ryujin 2.1.1 revision feb53359f0c9a08baf43c3dfe847d8a9f7d6893a
|
#include <source/shallow_water/riemann_solver.h>
Public Types | |
using | View = HyperbolicSystemView< dim, Number > |
using | primitive_type = std::array< Number, riemann_data_size > |
using | state_type = typename View::state_type |
using | ScalarNumber = typename View::ScalarNumber |
using | Parameters = RiemannSolverParameters< ScalarNumber > |
Public Member Functions | |
Compute wavespeed estimates | |
RiemannSolver (const HyperbolicSystem &hyperbolic_system, const Parameters ¶meters, const MultiComponentVector< ScalarNumber, n_precomputed_values > &precomputed_values) | |
Number | compute (const primitive_type &riemann_data_i, const primitive_type &riemann_data_j) const |
Number | compute (const state_type &U_i, const state_type &U_j, const unsigned int i, const unsigned int *js, const dealii::Tensor< 1, dim, Number > &n_ij) const |
Static Public Attributes | |
static constexpr unsigned int | problem_dimension = View::problem_dimension |
static constexpr unsigned int | riemann_data_size = 3 |
static constexpr unsigned int | n_precomputed_values |
Internal functions used in the Riemann solver | |
Number | compute_h_star (const primitive_type &riemann_data_i, const primitive_type &riemann_data_j) const |
Number | f (const primitive_type &primitive_state, const Number &h_star) const |
Number | phi (const primitive_type &riemann_data_i, const primitive_type &riemann_data_j, const Number &h) const |
Number | lambda1_minus (const primitive_type &riemann_data, const Number h_star) const |
Number | lambda3_plus (const primitive_type &riemann_data, const Number h_star) const |
Number | compute_lambda (const primitive_type &riemann_data_i, const primitive_type &riemann_data_j, const Number h_star) const |
primitive_type | riemann_data_from_state (const state_type &U, const dealii::Tensor< 1, dim, Number > &n_ij) const |
A fast approximative solver for the associated 1D Riemann problem. The solver has to ensure that the estimate \(\lambda_{\text{max}}\) that is returned for the maximal wavespeed is a strict upper bound.
Definition at line 41 of file riemann_solver.h.
using ryujin::ShallowWater::RiemannSolver< dim, Number >::View = HyperbolicSystemView<dim, Number> |
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 47 of file riemann_solver.h.
using ryujin::ShallowWater::RiemannSolver< dim, Number >::primitive_type = std::array<Number, riemann_data_size> |
The array type to store the expanded primitive state for the Riemann solver \([\rho, v, p, a]\)
Definition at line 64 of file riemann_solver.h.
using ryujin::ShallowWater::RiemannSolver< dim, Number >::state_type = typename View::state_type |
The storage type used for a (conserved) state vector \(\boldsymbol U\).
Definition at line 69 of file riemann_solver.h.
using ryujin::ShallowWater::RiemannSolver< dim, Number >::ScalarNumber = typename View::ScalarNumber |
The underlying scalar number type.
Definition at line 80 of file riemann_solver.h.
using ryujin::ShallowWater::RiemannSolver< dim, Number >::Parameters = RiemannSolverParameters<ScalarNumber> |
Definition at line 85 of file riemann_solver.h.
|
inline |
Constructor taking a HyperbolicSystem instance as argument
Definition at line 95 of file riemann_solver.h.
|
inline |
For two given 1D primitive states riemann_data_i and riemann_data_j, compute an estimation of an upper bound for the maximum wavespeed lambda.
Definition at line 227 of file riemann_solver.template.h.
Number ryujin::ShallowWater::RiemannSolver< dim, Number >::compute | ( | const state_type & | U_i, |
const state_type & | U_j, | ||
const unsigned int | i, | ||
const unsigned int * | js, | ||
const dealii::Tensor< 1, dim, Number > & | n_ij | ||
) | const |
For two given states U_i a U_j and a (normalized) "direction" n_ij compute an estimation of an upper bound for lambda.
Definition at line 241 of file riemann_solver.template.h.
|
inlineprotected |
Definition at line 27 of file riemann_solver.template.h.
|
inlineprotected |
Definition at line 49 of file riemann_solver.template.h.
|
inlineprotected |
Definition at line 66 of file riemann_solver.template.h.
References ryujin::positive_part().
|
inlineprotected |
Definition at line 81 of file riemann_solver.template.h.
References ryujin::positive_part().
|
inlineprotected |
Definition at line 96 of file riemann_solver.template.h.
References ryujin::negative_part(), and ryujin::positive_part().
|
inline |
Definition at line 110 of file riemann_solver.template.h.
References ryujin::positive_part().
|
inlineprotected |
Definition at line 209 of file riemann_solver.template.h.
|
staticconstexpr |
The dimension of the state space.
Definition at line 52 of file riemann_solver.h.
|
staticconstexpr |
Number of components in a primitive state, we store \([\rho, v, p, a]\), thus, 4.
Definition at line 58 of file riemann_solver.h.
|
staticconstexpr |
The number of precomputed values.
Definition at line 74 of file riemann_solver.h.