wiki:u/bliu/coolingLength
  1. Dimensions and Parameters
Original runs New runs with Aluminum Cooling Table
Dimensions 10x1x1 10x1x1
Pressure 1 10
gamma 1.001 1.001 for nonCooling
Grid Density 100 cu cgs or 100 cu
Wind Density 1 cu cgs or 1 cu
Ambient Density 0.01 cu cgs or 0.01 cu
Grid Temp 0.01 cu 0.01 cu or 1.5E2 K
Wind Temp 1 cu 1 cu or 1.5E5 K
Ambient Temp 100 cu 1.5E7 K
Sound speed of wind 1 cu 1 cu or 1E6
wind Velocity 20 cu 20 cu or 2E7 cgs
TimeScale 1 1E-6
LengthScale 1 1
TempScale 1.5E-8 1.5E4
pScale 1 207097.545
nScale 0.483451886595769408E+24 1E17
VelScale 1 1E6
FinalTime 1 1

Cooling Length and Analytical Cooling Parameters for Wire Turbulence

  1. hydro data without cooling
frame 11 http://www.pas.rochester.edu/~bliu/wireTurbulence/Cooling/frame11.png
frame 12 http://www.pas.rochester.edu/~bliu/wireTurbulence/Cooling/frame12.png
frame 14 http://www.pas.rochester.edu/~bliu/wireTurbulence/Cooling/frame14.png
Wire Temperature at x=1.15 http://www.pas.rochester.edu/~bliu/wireTurbulence/Cooling/wireT.png
Gas Temperature at x=1.08 http://www.pas.rochester.edu/~bliu/wireTurbulence/Cooling/x=1.08_f12.png
Gas Temperature at x=1.09 http://www.pas.rochester.edu/~bliu/wireTurbulence/Cooling/x=1.09_f12.png
Wire Temperature at x=1.15 http://www.pas.rochester.edu/~bliu/wireTurbulence/Cooling/wireT.png
internal energy at x=1.1 http://www.pas.rochester.edu/~bliu/wireTurbulence/Cooling/wire_e.png
Vy at x=1.1 http://www.pas.rochester.edu/~bliu/wireTurbulence/Cooling/wire_vy.png;http://www.pas.rochester.edu/~bliu/wireTurbulence/Cooling/wire_vyplot.png
  1. Low resolution data
frame Ave rho (number density) Wire Temp Post Shock Temp Ave e Max Vy
11 0.051 2.5E-4 115 5081 0.05
  1. Equations in the code

Cooling time t_cooling=e/(N*lambda(T)) for DM cooling

AnalyticCoolingStrength=rho2 * A_alpha*T0A_beta*ScaleCool for Analytic

A_alpha=alpha_cooling*4.76e-20*&! (ergs*cm3/s/K.5)

(3d0/16d0*Xmu*muH*amu/Boltzmann*(vy_cooling)2)(power)

A_beta=beta_cooling

  1. Cooling Length Estimate '
alpha_cooling 2.09E-46
beta_cooling 2.5
vy_cooling 0.05
A_alpha 0.2E-42
A_beta 2.5
nScale 4.8345E+23
pScale 1
TimeScale 1
scaleCool 2.3373E+47
Cooling Strength 2.5378E+02
Cooling Time 2.0021E+01
Cooling Length 1.0043E+00
  1. Aluminum Cooling
FUNCTION AlCoolingStrength(q,Temp)
      REAL(KIND=qPrec) :: q(:), AlCoolingStrength, Temp, ne, nAl, logne, logtemp
      ne = (6d0/7d0)*q(1)*nScale  ! number densities based on Z = 6
      nAl = (1d0/7d0)*q(1)*nScale
      logne = LOG10(ne)
      Temp = Temp/1.16045d4       ! K to eV
      logtemp = LOG10(Temp)
      AlCoolingStrength=ne*nAl*GetTableValue(AlTab, (/logtemp, logne/))/(nScale**2d0)
      AlCoolingStrength=AlCoolingStrength*ScaleCool
   END FUNCTION AlCoolingStrength

Last modified 8 years ago Last modified on 01/18/17 11:38:38
Note: See TracWiki for help on using the wiki.