def critL(pre0, nh0, r_inj, vinf): g = 5.0 / 3.0 #T = 147910.0 mu = 0.5 #T = 1.0e8 #c = math.sqrt(R_GAS * T / mu) c = math.sqrt(g * pre0 / (mh * nh0)) T = pre0 * mu / (R_GAS * mh * nh0) alpha = 0.28 #COOL = 3.185e-15 * math.pow(T, -0.63) * (1.0 - math.exp(-math.pow(T / 1.0e5, 1.63))) COOL = henney.coolTotal(1.0, nh0, T) coeff1 = 6.0 * (g - 1.0) * math.pi * alpha**2 * mu**2 * mh**2 * r_inj * vinf**4 coeff2 = (g + 1.0) * COOL return (coeff1 / coeff2) * ((vinf**2 / 2.0) - (c**2 / (g - 1.0)))
def coolTime(pre0, nh0): mu = 0.5 T = 1.0e8 T = pre0 * mu / (R_GAS * mh * nh0) COOL = henney.coolTotal(1.0, nh0, T) return 1.5 * nh0 * KB * T / (COOL)