shamphys::EOS_Config_LocallyIsothermalDisc_ExtendedFarris2014< Tscal > Struct Template Reference

Configuration struct for the locally isothermal equation of state extended from Farris 2014 to include for the q index of the disc. More...

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

Public Attributes
Tscal	cs0
	Soundspeed at the reference radius.
Tscal	q
	Power exponent of the soundspeed profile.
Tscal	r0
	Reference radius.
u32	n_sinks
	Number of sinks to consider for the equation of state.

Detailed Description

template<class Tscal>
struct shamphys::EOS_Config_LocallyIsothermalDisc_ExtendedFarris2014< Tscal >

Configuration struct for the locally isothermal equation of state extended from Farris 2014 to include for the q index of the disc.

This EOS should match with ieos 13 and 14 of phantom.

The equation in phantom is a bit weird so re-derived it here.

Farris 2014 EOS which only corresponds to q=1/2:

\( c_s = \frac{H_0}{r_0}\left(\frac{G M_1}{r_1} + \frac{G M_2}{r_2}\right) \)

However the extension of that EOS to q != 1/2 was only introduced in Ragussa et al 2016, if I'm right with:

\( c_s = \frac{H_0}{r_0} \left(\frac{G M_1}{r_1} + \frac{G M_2}{r_2}\right)^{q} \)

But as is the units are broken if q is not 1/2 so you need to compensate with \(r_0 \Omega_0\)

\(c_s = \frac{H_0}{r_0}\frac{1}{(r_0 \Omega_0)^{2q - 1}}\left(\frac{G M_1}{r_1} + \frac{G M_2}{r_2}\right)^{q} \)

\(= c_{s0} \frac{1}{(r_0 \Omega_0)^{q}}\left(\frac{G M_1}{r_1} + \frac{G M_2}{r_2}\right)^{q}\)

\(= c_{s0}\frac{1}{r_0^{q}}\left[\frac{1}{\sum_i M_i}\sum_i \frac{M_i}{r_i}\right]^{q}\)

Template Parameters

Tscal

Scalar type

Definition at line 259 of file eos_config.hpp.

Member Data Documentation

cs0

template<class Tscal >

Tscal shamphys::EOS_Config_LocallyIsothermalDisc_ExtendedFarris2014< Tscal >::cs0

Soundspeed at the reference radius.

Definition at line 261 of file eos_config.hpp.

n_sinks

template<class Tscal >

u32 shamphys::EOS_Config_LocallyIsothermalDisc_ExtendedFarris2014< Tscal >::n_sinks

Number of sinks to consider for the equation of state.

Definition at line 270 of file eos_config.hpp.

q

template<class Tscal >

Tscal shamphys::EOS_Config_LocallyIsothermalDisc_ExtendedFarris2014< Tscal >::q

Power exponent of the soundspeed profile.

Definition at line 264 of file eos_config.hpp.

r0

template<class Tscal >

Tscal shamphys::EOS_Config_LocallyIsothermalDisc_ExtendedFarris2014< Tscal >::r0

Reference radius.

Definition at line 267 of file eos_config.hpp.

---

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

• shamphys/include/shamphys/eos_config.hpp