NewtonianEuler::AnalyticData::ShuOsherTube Class Reference

Initial data for the Shu-Osher oscillatory shock tube [180]. More...

#include <ShuOsherTube.hpp>

Classes
struct	Epsilon
struct	JumpPosition
	Initial postition of the discontinuity. More...
struct	Lambda
struct	LeftMassDensity
struct	LeftPressure
struct	LeftVelocity
struct	RightPressure
struct	RightVelocity

Public Types
using	equation_of_state_type = EquationsOfState::IdealFluid< false >
using	options = implementation defined

Public Member Functions
	ShuOsherTube (double jump_position, double mass_density_l, double velocity_l, double pressure_l, double velocity_r, double pressure_r, double epsilon, double lambda)
	ShuOsherTube (const ShuOsherTube &)=default
ShuOsherTube &	operator= (const ShuOsherTube &)=default
	ShuOsherTube (ShuOsherTube &&)=default
ShuOsherTube &	operator= (ShuOsherTube &&)=default
auto	get_clone () const -> std::unique_ptr< evolution::initial_data::InitialData > override
template<typename DataType , typename... Tags>
tuples::TaggedTuple< Tags... >	variables (const tnsr::I< DataType, 1, Frame::Inertial > &x, tmpl::list< Tags... >) const
const EquationsOfState::IdealFluid< false > &	equation_of_state () const
void	pup (PUP::er &p) override
virtual auto	get_clone () const -> std::unique_ptr< InitialData >=0

Static Public Attributes
static constexpr Options::String	help

Friends
bool	operator== (const ShuOsherTube &lhs, const ShuOsherTube &rhs)

Detailed Description

Initial data for the Shu-Osher oscillatory shock tube [180].

The general initial data is given by:

$$\begin{array}{r}\{\rho ,v^x,p\}=\{\begin{array}{ll}\{{\rho }_L,v_L^x,p_L\}& \mathrm{i}\mathrm{f}x<\mathrm{\Delta }\\ \{1+ϵ\sin \left(\lambda x\right),v_R^x,p_R\}& \mathrm{i}\mathrm{f}x\ge \mathrm{\Delta }\end{array}\end{array}$$

with the adiabatic index being 1.4.

With the standard parameters given below, this is a Mach-3 shock moving into a sinusoidal density profile.

$$\begin{array}{r}\{\rho ,v^x,p\}=\{\begin{array}{ll}\{3.857143,2.629369,10.33333\}& \mathrm{i}\mathrm{f}x<-4\\ \{1+0.2\sin \left(5x\right),0,1\}& \mathrm{i}\mathrm{f}x\ge -4\end{array}\end{array}$$

With these values the usual final time is 1.8.

Member Function Documentation

get_clone()

auto NewtonianEuler::AnalyticData::ShuOsherTube::get_clone ( ) const -> std::unique_ptr< evolution::initial_data::InitialData >

overridevirtual

Implements evolution::initial_data::InitialData.

Member Data Documentation

help

constexpr Options::String NewtonianEuler::AnalyticData::ShuOsherTube::help

staticconstexpr

Initial value:

= {
      "1D Shu-Osher oscillatory shock tube."}

---

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

• src/PointwiseFunctions/AnalyticData/NewtonianEuler/ShuOsherTube.hpp