Changes between Version 82 and Version 83 of u/bliu/MHDJets

Timestamp:

06/22/20 12:16:32 (5 years ago)

Author:

Baowei Liu

Comment:

—

u/bliu/MHDJets

@@ +1 @@
+= New analytic cooling 6.20.2020 =
+
+
+{{{
+gamma=alpha*(T/T0)^beta
+alpha = 1e-21
+beta = [-2.5 -0.2 0 -.2 2.5]
+T0=10^5 K
+
+Cooling only applies to Ambient trailer < 1/3 which is the initial ambient density
+
+Current version of code -- need better ways of handling Temp0 and trailer!!!
+
+  FUNCTION AnalyticCoolingStrength2(q, Temp)
+     REAL(KIND=qPREC) :: AnalyticCoolingStrength2
+     REAL(KIND=qPrec) :: q(:)
+     ! Local declarations
+     REAL(KIND=qPrec) :: Temp, Temp0, T0
+
+     Temp0 = 100000  !10^5 Kelvin
+     T0 = Temp/Temp0
+     AnalyticCoolingStrength2=q(1)**2 * A_alpha*T0**A_beta*ScaleCool
+
+     !! only for ambient traer > 1/3 which is the scaless value for initial ambient density
+     if( q(9) <= 1d0/3d0 ) AnalyticCoolingStrength2=0d0
+
+
+  END FUNCTION AnalyticCoolingStrength2
+
+}}}
+
+|| rho_amb || 1e17 1/cm3 ||
+|| T_amb || 4320 K ||
+|| rho_jet || 6e17 1/cm3 ||
+|| T_jet || 720 K ||  
+|| T_ps || 183920 K ||
+|| V_jet || 44.982 km/s ||
+|| Mach || 14.2404 ||
+|| T_floor || 720 K ||
+
+results:
+
+|| beta=-2.5 || ||
+|| beta=-0.2 || ||
+|| beta=0.0 || ||
+|| beta=0.2 || ||
+|| beta=2.5 || ||
+
+
+----
+
+
 = Summaries 5.4.2020 =