![]() |
ryujin 2.1.1 revision b04fb4aba7bd97ddfa5994a6b4b6e8586eb15db2
|
#include <source/euler_barotropic/hyperbolic_system.h>
Public Member Functions | |
Constructor and setup | |
| HyperbolicSystemView (const HyperbolicSystem &hyperbolic_system) | |
Access to runtime parameters | |
| DEAL_II_ALWAYS_INLINE const std::string & | barotropic_equation_of_state () const |
| DEAL_II_HOST_DEVICE_ALWAYS_INLINE ScalarNumber | reference_density () const |
| DEAL_II_HOST_DEVICE_ALWAYS_INLINE ScalarNumber | vacuum_state_relaxation_small () const |
| DEAL_II_HOST_DEVICE_ALWAYS_INLINE ScalarNumber | vacuum_state_relaxation_large () const |
Low-level access to the selected equation of state | |
| |
| DEAL_II_ALWAYS_INLINE Number | beos_specific_internal_energy (const Number &rho) const |
| DEAL_II_ALWAYS_INLINE Number | beos_pressure (const Number &rho) const |
| DEAL_II_ALWAYS_INLINE Number | beos_speed_of_sound (const Number &rho) const |
Special functions for boundary states | |
| template<int component> | |
| state_type | prescribe_riemann_characteristic (const state_type &U, const Number &p, const state_type &U_bar, const Number &p_bar, const dealii::Tensor< 1, dim, Number > &normal) const |
| template<typename Lambda > | |
| state_type | apply_boundary_conditions (const dealii::types::boundary_id id, const state_type &U, const dealii::Tensor< 1, dim, Number > &normal, const Lambda &get_dirichlet_data) const |
State transformations | |
| template<typename ST > | |
| DEAL_II_HOST_DEVICE state_type | expand_state (const ST &state) const |
| template<typename ST > | |
| DEAL_II_HOST_DEVICE state_type | from_initial_state (const ST &initial_state) const |
| DEAL_II_HOST_DEVICE state_type | from_primitive_state (const state_type &primitive_state) const |
| DEAL_II_HOST_DEVICE state_type | to_primitive_state (const state_type &state) const |
| template<typename Lambda > | |
| DEAL_II_HOST_DEVICE state_type | apply_galilei_transform (const state_type &state, const Lambda &lambda) const |
Static Public Attributes | |
Constexpr booleans used in the EulerInitialStates namespace | |
| static constexpr bool | have_gamma = false |
| static constexpr bool | have_covolume_constant = false |
| static constexpr bool | have_energy_equation = false |
Typedefs and constexpr constants | |
| using | ScalarNumber = typename get_value_type< Number >::type |
| using | state_type = dealii::Tensor< 1, problem_dimension, Number > |
| using | flux_type = dealii::Tensor< 1, problem_dimension, dealii::Tensor< 1, dim, Number > > |
| using | flux_contribution_type = flux_type |
| using | precomputed_type = std::array< Number, n_precomputed_values > |
| using | initial_precomputed_type = std::array< Number, n_initial_precomputed_values > |
| template<int n_comp> | |
| using | VectorView = MultiComponentVectorView< ScalarNumber, n_comp, MemorySpace, false > |
| using | PrecomputedVectorView = VectorView< n_precomputed_values > |
| using | InitialPrecomputedVectorView = VectorView< n_initial_precomputed_values > |
| static constexpr unsigned int | problem_dimension = 1 + dim |
| static constexpr unsigned int | n_precomputed_values = 3 |
| static constexpr unsigned int | n_initial_precomputed_values = 0 |
Flux computations | |
| static constexpr bool | have_high_order_flux = false |
| DEAL_II_HOST_DEVICE flux_type | f (const state_type &U, const Number &p) const |
| DEAL_II_HOST_DEVICE flux_contribution_type | flux_contribution (const PrecomputedVectorView &pv, const InitialPrecomputedVectorView &piv, const unsigned int i, const state_type &U_i) const |
| DEAL_II_HOST_DEVICE flux_contribution_type | flux_contribution (const PrecomputedVectorView &pv, const InitialPrecomputedVectorView &piv, const unsigned int *js, const state_type &U_j) const |
| DEAL_II_HOST_DEVICE state_type | flux_divergence (const flux_contribution_type &flux_i, const flux_contribution_type &flux_j, const dealii::Tensor< 1, dim, Number > &c_ij) const |
| DEAL_II_HOST_DEVICE state_type | high_order_flux_divergence (const flux_contribution_type &flux_i, const flux_contribution_type &flux_j, const dealii::Tensor< 1, dim, Number > &c_ij) const =delete |
Computing stencil source terms | |
| static constexpr bool | have_source_terms = false |
| DEAL_II_HOST_DEVICE state_type | nodal_source (const PrecomputedVectorView &pv, const unsigned int i, const state_type &U_i, const ScalarNumber tau) const =delete |
| DEAL_II_HOST_DEVICE state_type | nodal_source (const PrecomputedVectorView &pv, const unsigned int *js, const state_type &U_j, const ScalarNumber tau) const =delete |
Computing derived physical quantities | |
| DEAL_II_HOST_DEVICE Number | filter_vacuum_density (const Number &rho) const |
| DEAL_II_HOST_DEVICE Number | total_energy (const state_type &U, const Number &specific_internal_energy) const |
| DEAL_II_HOST_DEVICE state_type | total_energy_derivative (const state_type &U, const Number &specific_internal_energy, const Number &pressure) const |
| DEAL_II_HOST_DEVICE bool | is_admissible (const state_type &U) const |
| static DEAL_II_HOST_DEVICE Number | density (const state_type &U) |
| static DEAL_II_HOST_DEVICE dealii::Tensor< 1, dim, Number > | momentum (const state_type &U) |
A view of the HyperbolicSystem that makes methods available for a given dimension dim and choice of number type Number (which can be a scalar float, or double, as well as a VectorizedArray holding packed scalars.
Intended usage:
Definition at line 168 of file hyperbolic_system.h.
| using ryujin::EulerBarotropic::HyperbolicSystemView< dim, Number, MemorySpace >::ScalarNumber = typename get_value_type<Number>::type |
The underlying scalar number type.
Definition at line 184 of file hyperbolic_system.h.
| using ryujin::EulerBarotropic::HyperbolicSystemView< dim, Number, MemorySpace >::state_type = dealii::Tensor<1, problem_dimension, Number> |
Storage type for a (conserved) state vector \(\boldsymbol U\).
Definition at line 194 of file hyperbolic_system.h.
| using ryujin::EulerBarotropic::HyperbolicSystemView< dim, Number, MemorySpace >::flux_type = dealii::Tensor<1, problem_dimension, dealii::Tensor<1, dim, Number> > |
Storage type for the flux \(\mathbf{f}\).
Definition at line 199 of file hyperbolic_system.h.
| using ryujin::EulerBarotropic::HyperbolicSystemView< dim, Number, MemorySpace >::flux_contribution_type = flux_type |
The storage type used for flux contributions.
Definition at line 205 of file hyperbolic_system.h.
| using ryujin::EulerBarotropic::HyperbolicSystemView< dim, Number, MemorySpace >::precomputed_type = std::array<Number, n_precomputed_values> |
Array type used for precomputed values.
Definition at line 215 of file hyperbolic_system.h.
| using ryujin::EulerBarotropic::HyperbolicSystemView< dim, Number, MemorySpace >::initial_precomputed_type = std::array<Number, n_initial_precomputed_values> |
Array type used for precomputed initial values.
Definition at line 225 of file hyperbolic_system.h.
| using ryujin::EulerBarotropic::HyperbolicSystemView< dim, Number, MemorySpace >::VectorView = MultiComponentVectorView<ScalarNumber, n_comp, MemorySpace, false> |
Definition at line 229 of file hyperbolic_system.h.
| using ryujin::EulerBarotropic::HyperbolicSystemView< dim, Number, MemorySpace >::PrecomputedVectorView = VectorView<n_precomputed_values> |
MulticomponentVectorView for accessing a vector of precomputed states:
Definition at line 236 of file hyperbolic_system.h.
| using ryujin::EulerBarotropic::HyperbolicSystemView< dim, Number, MemorySpace >::InitialPrecomputedVectorView = VectorView<n_initial_precomputed_values> |
MulticomponentVectorView for accessing a vector of precomputed initial states:
Definition at line 242 of file hyperbolic_system.h.
|
inline |
Constructor taking a reference to the underlying HyperbolicSystem
Definition at line 255 of file hyperbolic_system.h.
|
inline |
Return the name of the selected barotropic equation of state.
Definition at line 275 of file hyperbolic_system.h.
|
inline |
Definition at line 280 of file hyperbolic_system.h.
References ryujin::EulerBarotropic::HyperbolicSystem::Parameters::reference_density.
|
inline |
Definition at line 286 of file hyperbolic_system.h.
References ryujin::EulerBarotropic::HyperbolicSystem::Parameters::vacuum_state_relaxation_small.
|
inline |
Definition at line 292 of file hyperbolic_system.h.
References ryujin::EulerBarotropic::HyperbolicSystem::Parameters::vacuum_state_relaxation_large.
|
inline |
For a given density \(\rho\) return the specific internal energy \(e\).
Definition at line 313 of file hyperbolic_system.h.
|
inline |
For a given density \(\rho\) return the pressure \(p\).
Definition at line 331 of file hyperbolic_system.h.
|
inline |
For a given density \(\rho\) and specific internal energy \(e\) return the sound speed \(a\).
Definition at line 351 of file hyperbolic_system.h.
|
static |
For a given (1+dim dimensional) state vector U, return the density U[0]
Definition at line 823 of file hyperbolic_system.h.
| DEAL_II_HOST_DEVICE_ALWAYS_INLINE Number ryujin::EulerBarotropic::HyperbolicSystemView< dim, Number, MemorySpace >::filter_vacuum_density | ( | const Number & | rho | ) | const |
Given a density rho this function returns 0 if the magnitude of rho is smaller or equal than relaxation_large * rho_cutoff. Otherwise rho is returned unmodified. Here, rho_cutoff is the reference density multiplied by eps.
Definition at line 831 of file hyperbolic_system.h.
|
static |
For a given (1+dim dimensional) state vector U, return the momentum vector [U[1], ..., U[1+dim]].
Definition at line 845 of file hyperbolic_system.h.
| DEAL_II_HOST_DEVICE_ALWAYS_INLINE Number ryujin::EulerBarotropic::HyperbolicSystemView< dim, Number, MemorySpace >::total_energy | ( | const state_type & | U, |
| const Number & | specific_internal_energy | ||
| ) | const |
For a given (1+dim dimensional) barotropic state vector U, compute and return the total energy of the system.
\[ \eta = \rho e(\rho) + \frac12\rho&{-1}|\vec m|^2 \]
Definition at line 857 of file hyperbolic_system.h.
| DEAL_II_HOST_DEVICE_ALWAYS_INLINE auto ryujin::EulerBarotropic::HyperbolicSystemView< dim, Number, MemorySpace >::total_energy_derivative | ( | const state_type & | U, |
| const Number & | specific_internal_energy, | ||
| const Number & | pressure | ||
| ) | const |
For a given (1+dim dimensional) barotropic state vector U, compute and return the derivative \(\eta'\) of the total energy.
Definition at line 871 of file hyperbolic_system.h.
| DEAL_II_HOST_DEVICE_ALWAYS_INLINE bool ryujin::EulerBarotropic::HyperbolicSystemView< dim, Number, MemorySpace >::is_admissible | ( | const state_type & | U | ) | const |
Returns whether the state U is admissible. If U is a vectorized state then U is admissible if all vectorized values are admissible.
Definition at line 894 of file hyperbolic_system.h.
| state_type ryujin::EulerBarotropic::HyperbolicSystemView< dim, Number, MemorySpace >::prescribe_riemann_characteristic | ( | const state_type & | U, |
| const Number & | p, | ||
| const state_type & | U_bar, | ||
| const Number & | p_bar, | ||
| const dealii::Tensor< 1, dim, Number > & | normal | ||
| ) | const |
Decomposes a given state U into Riemann invariants and then replaces the first or second Riemann characteristic from the one taken from U_bar state. Note that the U_bar state is just the prescribed dirichlet values.
| state_type ryujin::EulerBarotropic::HyperbolicSystemView< dim, Number, MemorySpace >::apply_boundary_conditions | ( | const dealii::types::boundary_id | id, |
| const state_type & | U, | ||
| const dealii::Tensor< 1, dim, Number > & | normal, | ||
| const Lambda & | get_dirichlet_data | ||
| ) | const |
Apply boundary conditions.
For the compressible Euler equations we have:
| DEAL_II_HOST_DEVICE_ALWAYS_INLINE auto ryujin::EulerBarotropic::HyperbolicSystemView< dim, Number, MemorySpace >::f | ( | const state_type & | U, |
| const Number & | p | ||
| ) | const |
Given a state U and a pressure p compute the flux
\[ \begin{pmatrix} \textbf m \\ \textbf v\otimes \textbf m + p\mathbb{I}_d \end{pmatrix}, \]
Definition at line 1033 of file hyperbolic_system.h.
| DEAL_II_HOST_DEVICE_ALWAYS_INLINE auto ryujin::EulerBarotropic::HyperbolicSystemView< dim, Number, MemorySpace >::flux_contribution | ( | const PrecomputedVectorView & | pv, |
| const InitialPrecomputedVectorView & | piv, | ||
| const unsigned int | i, | ||
| const state_type & | U_i | ||
| ) | const |
Given a state U_i and an index i compute flux contributions.
Intended usage:
For the Euler equations we simply compute f(U_i).
Definition at line 1054 of file hyperbolic_system.h.
| DEAL_II_HOST_DEVICE_ALWAYS_INLINE auto ryujin::EulerBarotropic::HyperbolicSystemView< dim, Number, MemorySpace >::flux_contribution | ( | const PrecomputedVectorView & | pv, |
| const InitialPrecomputedVectorView & | piv, | ||
| const unsigned int * | js, | ||
| const state_type & | U_j | ||
| ) | const |
Definition at line 1068 of file hyperbolic_system.h.
| DEAL_II_HOST_DEVICE_ALWAYS_INLINE auto ryujin::EulerBarotropic::HyperbolicSystemView< dim, Number, MemorySpace >::flux_divergence | ( | const flux_contribution_type & | flux_i, |
| const flux_contribution_type & | flux_j, | ||
| const dealii::Tensor< 1, dim, Number > & | c_ij | ||
| ) | const |
Given flux contributions flux_i and flux_j compute the flux (-f(U_i) - f(U_j)
Definition at line 1082 of file hyperbolic_system.h.
References ryujin::add(), and ryujin::contract().
|
delete |
|
delete |
|
delete |
| DEAL_II_HOST_DEVICE state_type ryujin::EulerBarotropic::HyperbolicSystemView< dim, Number, MemorySpace >::expand_state | ( | const ST & | state | ) | const |
Given a state vector associated with a different spatial dimensions than the current one, return an "expanded" version of the state vector associated with dim spatial dimensions where the momentum vector of the conserved state state is expaned with zeros to a total length of dim entries.
| DEAL_II_HOST_DEVICE state_type ryujin::EulerBarotropic::HyperbolicSystemView< dim, Number, MemorySpace >::from_initial_state | ( | const ST & | initial_state | ) | const |
Given an initial state [rho, u_1, ..., u_d, p] return a conserved state [rho, m_1, ..., m_d, E]. Most notably, the specific equation of state oracle is queried to convert the pressure value into a specific internal energy.
| DEAL_II_HOST_DEVICE_ALWAYS_INLINE auto ryujin::EulerBarotropic::HyperbolicSystemView< dim, Number, MemorySpace >::from_primitive_state | ( | const state_type & | primitive_state | ) | const |
Given a primitive state [rho, v_1, ..., v_d] return a conserved state.
Definition at line 1127 of file hyperbolic_system.h.
| DEAL_II_HOST_DEVICE_ALWAYS_INLINE auto ryujin::EulerBarotropic::HyperbolicSystemView< dim, Number, MemorySpace >::to_primitive_state | ( | const state_type & | state | ) | const |
Given a conserved state return a primitive state [rho, v_1, ..., v_d]
Definition at line 1143 of file hyperbolic_system.h.
| DEAL_II_HOST_DEVICE state_type ryujin::EulerBarotropic::HyperbolicSystemView< dim, Number, MemorySpace >::apply_galilei_transform | ( | const state_type & | state, |
| const Lambda & | lambda | ||
| ) | const |
Transform the current state according to a given operator lambda acting on a dim dimensional momentum (or velocity) vector.
|
staticconstexpr |
The dimension of the state space.
Definition at line 189 of file hyperbolic_system.h.
|
staticconstexpr |
The number of precomputed values.
Definition at line 210 of file hyperbolic_system.h.
|
staticconstexpr |
The number of precomputed initial values.
Definition at line 220 of file hyperbolic_system.h.
|
staticconstexpr |
Definition at line 372 of file hyperbolic_system.h.
|
staticconstexpr |
Definition at line 373 of file hyperbolic_system.h.
|
staticconstexpr |
Definition at line 374 of file hyperbolic_system.h.
|
staticconstexpr |
The low-order and high-order fluxes are the same:
Definition at line 546 of file hyperbolic_system.h.
|
staticconstexpr |
We do not have source terms
Definition at line 561 of file hyperbolic_system.h.