Esempio n. 1
0
def calcStagnationTime(logQ, nh):
	ci = koo.soundSpeed(8000.0, 0.5)
	co = koo.soundSpeed(10.0, 1.0)
	cico = ci / co

	RS = koo.calcStromgrenRadius(logQ, nh)
	ts = RS / ci

	return (4.0 / 7.0) * ((3.0 / 4.0)**0.5) * (((4.0 / 3.0)**(7.0 / 6.0)) * (cico**(7.0 / 3.0)) - 1.0) * ts
Esempio n. 2
0
def windHiiVelRatio(isMVD):
	for row in range(9):
		mdot = mdotList[row]
		vinf = vinfList[row]
		times = mvd_times if isMVD else mvt_times

		tableline =  str(massList[row])
		for col in range(5):
			nh = mvd_iden[col][row] if isMVD else  mvt_iden[col][row]
			#nh = mvd_nhs[col] if isMVD else  mvt_nhs[col]

			R = koo.calcBubbleRadiusRB(mdot, vinf, mh * nh, times[col] * YR2S)
			v = R / (times[col] * YR2S)
			ci = koo.soundSpeed(8000.0, 0.5)
			tableline = tableline + " & " + str(v / ci)

		print tableline + "\\\\"