Aluminum Cooling

1. Original Cooling Table

0d0 2d0 1d-1 16d0 19d0 1d0 1 1
9.5013460857276d-21 9.452949522506d-22 8.2247033959d-23 7.816645303572d-24
1.880126789552d-20 2.383444755056d-21 2.1863738552096d-22 1.7507272234946d-23
2.33804509917d-20 3.60591383635d-21 4.704545192328d-22 4.099141290076d-23
4.040940257088d-20 5.685016012453d-21 6.956451960545d-22 7.798493275244d-23
1.5556940669088d-20 5.474266076639d-21 1.0478264664546d-21 1.1746353664656d-22
2.7098495276281d-21 1.4346129044304d-21 8.019099665396d-22 1.6389276429656d-22
6.252357428359d-22 3.557024783101d-22 2.8188636111456d-22 1.3623718647462d-22
1.9927303071d-22 1.1807499686052d-22 1.0289658697416d-22 7.129151345264d-23
8.531070138368d-23 5.386946917198d-23 4.871510796939d-23 3.873437631258d-23
8.53258939671d-23 7.949898535722d-23 7.62189248437d-23 5.223399320165d-23
5.801904750093d-22 1.0513297958832d-21 9.667645791854d-22 3.570565232044d-22
7.09672745925d-21 1.028642200642d-20 7.707712486531d-21 2.030515975311d-21
3.1222912043725d-20 3.438183443793d-20 2.5489688126791d-20 6.371784696976d-21
8.7152142903d-20 8.728794090588d-20 5.91507595616d-20 1.3956125144832d-20
1.812743479808d-19 1.6886600316919d-19 1.0236362150679d-19 2.2102807566645d-20
2.9830084425126d-19 2.6880719552326d-19 1.3382971945944d-19 2.4675590299233d-20
4.3075068682608d-19 3.4364315479333d-19 1.2058713101736d-19 2.0951549263728d-20
5.0458575060733d-19 3.0502142941424d-19 8.6011375813842d-20 1.9502579927838d-20
4.8549660579714d-19 2.4183313538395d-19 6.7458189069d-20 2.118855768717d-20
3.641803190773d-19 1.741573040681d-19 5.0134084958848d-20 2.19771341513d-20
2.038953583702d-19 1.046903674496d-19 3.661435334144d-20 1.982087435244d-20

Temperature range	1.16045e4 - 1.16045e6 K
electron number Density range	10^{16}-10^{19} in cgs
number density range	10^{16}-10^{19} in cgs
BC	extrapolate

2. Updated to the new form

Dimension low_logT high_logT d_logT l_bc r_bc low_logRho high_logRho d_logRho l_bc r_bc

lambda(rho,T) …

2
0d0 2d0 1d-1 1 1 
16d0 19d0 1d0 1 1
9.5013460857276d-21
9.452949522506d-22
8.2247033959d-23
7.816645303572d-24
1.880126789552d-20
2.383444755056d-21
2.1863738552096d-22
1.7507272234946d-23
2.33804509917d-20
3.60591383635d-21
4.704545192328d-22
4.099141290076d-23
4.040940257088d-20
...

3. 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	10^{19} cgs or 100 cu
Wind Density	1 cu	10^{17} cgs or 1 cu
Ambient Density	0.01 cu	10^{15} 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
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	;

2. 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

2. 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*cm^3/s/K.5)

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

A_beta=beta_cooling

3. 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

4. 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