- 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
9 years ago
Last modified on 01/18/17 11:38:38
Note:
See TracWiki
for help on using the wiki.








