rotstars: module of rotating stars functions

PyHdust rotstars module: Rotating stars tools.

license:GNU GPL v3.0
pyhdust.rotstars.beta(par, is_ob=False)[source]

Calculate the \(\beta\) value from Espinosa-Lara for a given rotation rate \(w_{\rm frac} = \Omega/\Omega_c\)

If is_ob == True, it consider the param as ob (instead of \(w_{\rm frac}\)).

pyhdust.rotstars.ellips_th(th, rf)[source]

Ellipsoid radius

  • th – theta, in radians (0 = pole; pi/2 = equator).
  • rt – radius fraction (Req/Rp >= 1)

Read source generated with ref_estrela.txt.

OUTPUT: M, Req and TP (2*solar units and K).

pyhdust.rotstars.rochearea(wfrac, isW=False)[source]

Calculate the Roche area of a rigid rotator.

Equation 4.23 from Cranmer 1996 (thesis).

Area in (squared) radial unit (it must be multiplied to Rpole**2 to a physical size).

pyhdust.rotstars.rotStar(Tp=20000.0, M=10.3065, rp=5.38462, star='B', beta=0.25, wfrac=0.8, th_res=5001, quiet=False, LnotTp=False)[source]

Return the photospheric parameters of a rotating star.

LnotTp: the value of “Tp” is the Luminosity (in solar units).

Calculation of Von Zeipel’s Beta parameter as function of W: see math...

INPUT: th_res (theta resolution, integer)...

OUTPUT: printed status + (ob, Tp values, Area[cm2])

pyhdust.rotstars.rt(th, wfrac)[source]

Roche Rpole normalized radius as function of wfrac.

Parameters:th – theta, in radians (0 = pole; pi/2 = equator).

Returns the vrot value of a given source star.

OUTPUT: vrot in km/s.


Returns wfrac (Omega/Omega_crit) value from a W value.

Equation 1.23 de Faes (2015).

pyhdust.rotstars.wrot(par, is_ob=False)[source]

Converts \(w_{\rm frac} = \Omega/\Omega_c\) into \(W = vrot/vorb\).

If is_ob == True, it considers the param as the oblateness (instead of \(w_{\rm frac}\)).