- 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
- hydro data without cooling
frame 11 | ||
frame 12 | ||
frame 14 | ||
Wire Temperature at x=1.15 | ||
Gas Temperature at x=1.08 | ||
Gas Temperature at x=1.09 | ||
Wire Temperature at x=1.15 | ||
internal energy at x=1.1 | ||
Vy at x=1.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 |
- 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
- 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 |
- 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.