Changes between Version 3 and Version 4 of u/bliu/pnStudy/M2-9

Timestamp:

10/01/15 13:55:13 (10 years ago)

Author:

Baowei Liu

Comment:

—

u/bliu/pnStudy/M2-9

@@ -7 +7 @@
 }}}
 
-|| [http://www.pas.rochester.edu/~bliu/pnStudy/M2-9/rho_nonPulsed.gif Non Pulsed movie] || [http://www.pas.rochester.edu/~bliu/pnStudy/M2-9/rho_pulsed_150y.gif Pulsed movie Period=150y] ||
 || [http://www.pas.rochester.edu/~bliu/pnStudy/M2-9/rho_nonPulsed_long.gif Non Pulsed movie long] || [http://www.pas.rochester.edu/~bliu/pnStudy/M2-9/rho_pulsed_150y_long.gif Pulsed movie Period=150y Long] ||
 
@@ -14 +13 @@
 {{{
     if(lPulsed) then
-      qjet(1)=qjet(1)*cos(Omega_Pulsed*time)
-      v=v*cos(Omega_Pulsed*time)
+      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$ || ||