agpy 0.1 documentation

Source code for agpy.hotwind

# formulae for hot wind computations
#
# Panagia & Felli 1975 / Wright & Barlow 1975:
from agpy.constants import *
from numpy import log10,exp,log

[docs]def mdot(snu, nu=10, Te=1e4, vinf=2000, mue=1.3, d=1): """ [snu] = mJy nu = 10 GHz (frequency) Te = 10^4 K (ionized gas temp) vinf = 2000 km/s (maximum wind speed) mue = 1.3 (mean atomic weight per free electron) d = distance (kpc) return: mdot in Msun/yr """ mdot = ( snu / 7.26 * (nu/10.)**-0.6 * (Te/1e4)**-0.1 * d**2 * (mue*vinf/100)**(4/3.) ) **(3/4.) * 1e-6 return mdot
[docs]def mdotvinfr(snu, nu=10, Te=1e4, vinf=2000, mue=1.3, d=1, R=25): """ [snu] = mJy nu = 10 GHz (frequency) Te = 10^4 K (ionized gas temp) vinf = 2000 km/s (maximum wind speed) mue = 1.3 (mean atomic weight per free electron) d = distance (kpc) return: Dmom in g*cm/s^2 Sp. type Log D0 x alphaP A I 14.22 2.41 2.64 0.47 0.38 0.07 Mid B I 17.07 1.05 1.95 0.20 0.51 0.05 Early B I 21.24 1.38 1.34 0.25 0.75 0.15 O I 20.69 1.04 1.51 0.18 0.66 0.06 O III, V 19.87 1.21 1.57 0.21 0.64 0.06 """ mdot = ( snu / 7.26 * (nu/10.)**-0.6 * (Te/1e4)**-0.1 * d**2 * (mue*vinf/100)**(4/3.) ) **(3/4.) * 1e-6 return mdot*msun/yr * (vinf*1e5) * (R)**0.5
[docs]def LofMdot(logdmom, x=1.51, D0=20.69, alphaP=0.66): """ Kudritzki 2000 log Dmom - log D0 = x log(L/Lsun); Dmom = Mdot * vinf *R/rsun returns L in log10(Lsun) """ L = ( (logdmom-D0) / x ) return L