M-9

Test 1 9.30.2015

• 2D

      if(lPulsed) then
        qjet(1)=qjet(1)*(abs(cos(Omega_Pulsed*time))**PulseModule)
      end if
  

The shape depends on the resolution. See the non-pulsed result with different resolution:

1AMR		​global.data
2AMR		​global.data

​Non Pulsed movie long

​Pulsed movie Period=150y Long

*3D

Non-rotate		​4AMR non Rotate movie;	​4AMR non Rotate 2D slice movie;
Rotate with 150y period		​4AMR Rotate movie;	​4AMR Rotate 2D slice movie;

Test 2

    if(lPulsed) then
      qjet(1)=qjet(1)*(abs(cos(Omega_Pulsed*time))**PulseModule)
    end if

&ProblemData
! MODEL: jet 0deg  (n=4e4, r=1000AU, T=100K) @ 200 km/s into 4e4 torus + AGB wind
! folder taper15n4e2v200namb4e4
!
!      BACKGROUND or “AMBIENT” SECTION. Values apply to origin
tamb = 1d3           ! ambient temp, 1cu = 0.1K (100K=1000cu)
namb = 4e4           ! ambient central density cm^-3. Usually 400 for 1/r^2 or torus.
stratified = t       ! true = add a 1/r^2 background 'AGB stellar wind'
torus      = f       ! true - add torus to the background
torusalpha = 0.7     ! alpha and beta specify the geometry
torusbeta  = 10d0    ! see Frank & Mellema, 1994ApJ...430..800F
rings      = f       ! true - add radial density modulations to AGB wind
!
!     FLOW  DESCRIPTION SECTION, values apply at origin at t=0
outflowType  = 2    ! TYPE OF FLOW    1 cyl jet, 2 conical wind, 3 is clump
njet  = 4d4         ! flow density at launch zone, 1cu = 1cm^-3
Rjet  = 2d0         ! flow radius at launch zone, 1cu = 500AU (outflowType=1 only)
vjet  = 2e7         ! flow velocity , 1cu = cm/s (100km/s=1e7cu)
tjet  = 1d3         ! flow temp, 1cu = 0.1K (100K=1000cu)
tt    = 0.0d0       ! flow accel time, 1cu = 8250y (0.02 = 165y)
open_angle = 15d0    ! conical flow open angle (deg)
tf    = 15d0        ! conical flow Gaussian taper (deg) for njet and vjet; 0= disable
sigma = 0d0         ! !toroidal.magnetic.energy / kinetic.energy, example 0.6

!
!     OTHER PARAMETERS
lcooling = t         ! radiative cooling?
buff  = 8            ! central refinement of a grid with a resolution 1/2

! for Taper
lObliqueCone = F    !T, 45 deg tamper
                    !F, non 45 deg tamper

! for M2-9 Pulsed jet (density)

lPulsed = T         !T, turn on pulsed jet
                    !F, turn off
PulsedPeriod = 150  ! year

PulseModule = 4     ! 1  abs(cos)
                    ! 2  abs(cos)^2
                    ! 4  abs(cos)^4
/

cos^4

​movie