Meeting update 06/22/2020

MHDCollidingFlows runs with the new analytic cooling and the ambient tracer. Need better way to handle parameters Temp0 and tracer in the cooling code though. Also need to double-check the scales..

​Mach15 result 

      A_alpha = alpha
      !cs=sqrt(gamma*Boltzmann*TempScale/Xmu/muH/amu)
      !power=.5d0*(1d0-2d0*beta)
      !A_alpha=alpha*4.76e-20*&! (ergs*cm^3/s/K^.5)
      !   (3d0/16d0*Xmu*muH*amu/Boltzmann*(cs*vx)**2)**(power)
      A_beta = beta


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

Comments

No comments.