Applies a pressure and density floor dependent on the distance to the origin.
More...
|
|
| RadiallyFallingFloor (double minimum_radius_at_which_to_apply_floor, double rest_mass_density_scale, double rest_mass_density_power, double pressure_scale, double pressure_power) |
|
| RadiallyFallingFloor (const RadiallyFallingFloor &)=default |
|
RadiallyFallingFloor & | operator= (const RadiallyFallingFloor &)=default |
|
| RadiallyFallingFloor (RadiallyFallingFloor &&)=default |
|
RadiallyFallingFloor & | operator= (RadiallyFallingFloor &&)=default |
|
void | pup (PUP::er &p) |
|
template<size_t ThermodynamicDim> |
| void | operator() (gsl::not_null< Scalar< DataVector > * > density, gsl::not_null< Scalar< DataVector > * > pressure, gsl::not_null< Scalar< DataVector > * > specific_internal_energy, gsl::not_null< Scalar< DataVector > * > temperature, gsl::not_null< Scalar< DataVector > * > electron_fraction, const tnsr::I< DataVector, Dim, Frame::Inertial > &coords, const EquationsOfState::EquationOfState< true, ThermodynamicDim > &equation_of_state) const |
template<size_t Dim>
class VariableFixing::RadiallyFallingFloor< Dim >
Applies a pressure and density floor dependent on the distance to the origin.
Applies the floors: \(\rho(r) \geq \rho_{\mathrm{fl}}(r) = C_\rho r^{k_\rho}\) and \(P(r) \geq P_{\mathrm{fl}}(r) = C_p r^{k_p}\) when \( r > r_{min}\), where \(C_\rho\) is given by the option ScaleDensityFloor, \(k_\rho\) is given by the option PowerDensityFloor, \(C_p\) is given by the option ScalePressureFloor, \(k_p\) is given by the option PowerPressureFloor, and \(r_{min}\) is given by the option MinimumRadius.
- Note
- In [173], the following floors are applied: \(\rho(r) \geq \rho_{\mathrm{fl}}(r) = 10^{-5}r^{-3/2}\) and \(P(r) \geq P_{\mathrm{fl}}(r) = \frac{1}{3} \times 10^{-7}r^{-5/2}\)
-
For 3d equations of state we treat the product \(\rho T\) as the "pressure" that we floor and then compute the pressure from the temperature. The floor is applied after flooring the density.