|
|
| SodExplosion (const SodExplosion &)=default |
|
SodExplosion & | operator= (const SodExplosion &)=default |
|
| SodExplosion (SodExplosion &&)=default |
|
SodExplosion & | operator= (SodExplosion &&)=default |
| auto | get_clone () const -> std::unique_ptr< evolution::initial_data::InitialData > override |
|
| SodExplosion (double initial_radius, double inner_mass_density, double inner_pressure, double outer_mass_density, double outer_pressure, const Options::Context &context={}) |
|
template<typename... Tags> |
| tuples::TaggedTuple< Tags... > | variables (const tnsr::I< DataVector, Dim, Frame::Inertial > &x, tmpl::list< Tags... >) const |
| | Retrieve a collection of hydrodynamic variables at position x.
|
|
const equation_of_state_type & | equation_of_state () const |
|
void | pup (PUP::er &) override |
template<size_t Dim>
class NewtonianEuler::AnalyticData::SodExplosion< Dim >
A cylindrical or spherical Sod explosion [207] [190].
Common initial conditions are:
\begin{align*}(\rho, v^i, p) =
&\left\{
\begin{array}{ll}
(1 ,0, 1) & \mathrm{if} \; r \le 0.5 \\
(0.125 ,0, 0.1) & \mathrm{if} \; r > 0.5
\end{array}\right.
\end{align*}
where \(r\) is the cylindrical (2d) or spherical (3d) radius. This test problem uses an adiabatic index of 1.4. A reference solution can be computed in 1d by solving the Newtonian Euler equations in cylindrical or spherical symmetry. Note that the inner and outer density and pressure, as well as where the initial discontinuity is can be chosen arbitrarily.