Shamrock 2025.10.0
Astrophysical Code
Loading...
Searching...
No Matches
shamphys::EOS_Machida06< T > Struct Template Reference

Piecewise polytropic EOS from Machida et al. (2006). More...

#include <shamphys/include/shamphys/eos.hpp>

Static Public Member Functions

static constexpr T soundspeed (T P, T rho, T rho_c1, T rho_c2, T rho_c3)
static constexpr T temperature (T P, T rho, T mu, T mh, T kb)
static constexpr T pressure (T cs, T rho, T rho_c1, T rho_c2, T rho_c3)

Detailed Description

template<class T>
struct shamphys::EOS_Machida06< T >

Piecewise polytropic EOS from Machida et al. (2006).

Uses different gamma values across density thresholds for gravitational collapse modeling.

Sound speed: \( c_s = \sqrt{\frac{\gamma P}{\rho}} \) where \( \gamma \) depends on density:

\[\gamma = \begin{cases} 1.0 & \rho < \rho_{c1} \\ 7/5 & \rho_{c1} \leq \rho < \rho_{c2} \\ 1.1 & \rho_{c2} \leq \rho < \rho_{c3} \\ 5/3 & \rho \geq \rho_{c3} \end{cases} \]

Pressure (piecewise):

\[P = \begin{cases} c_s^2 \rho & \rho < \rho_{c1} \\ c_s^2 \rho_{c1} \left(\frac{\rho}{\rho_{c1}}\right)^{7/5} & \rho_{c1} \leq \rho < \rho_{c2} \\ c_s^2 \rho_{c1} \left(\frac{\rho_{c2}}{\rho_{c1}}\right)^{7/5} \left(\frac{\rho}{\rho_{c2}}\right)^{1.1} & \rho_{c2} \leq \rho < \rho_{c3} \\ c_s^2 \rho_{c1} \left(\frac{\rho_{c2}}{\rho_{c1}}\right)^{7/5} \left(\frac{\rho_{c3}}{\rho_{c2}}\right)^{1.1} \left(\frac{\rho}{\rho_{c3}}\right)^{5/3} & \rho \geq \rho_{c3} \end{cases} \]

Temperature: \( T = \frac{\mu m_H P}{\rho k_B} \)

Definition at line 133 of file eos.hpp.

Member Function Documentation

◆ pressure()

template<class T>
constexpr T shamphys::EOS_Machida06< T >::pressure ( T cs,
T rho,
T rho_c1,
T rho_c2,
T rho_c3 )
inlinestaticconstexpr

Definition at line 147 of file eos.hpp.

◆ soundspeed()

template<class T>
constexpr T shamphys::EOS_Machida06< T >::soundspeed ( T P,
T rho,
T rho_c1,
T rho_c2,
T rho_c3 )
inlinestaticconstexpr

Definition at line 135 of file eos.hpp.

◆ temperature()

template<class T>
constexpr T shamphys::EOS_Machida06< T >::temperature ( T P,
T rho,
T mu,
T mh,
T kb )
inlinestaticconstexpr

Definition at line 143 of file eos.hpp.


The documentation for this struct was generated from the following file: