Computes the Full Dyadic Second Moment which becomes the Fluid Momentum Flux Gamma = Full velocity space integral of vec{v}(A_{alpha} vec{v} f_{alpha}) d^{3}vec{v} or \Gamma = \iiint \boldsymbol{v} \left(A_{\alpha} \boldsymbol{v} f_{\alpha}\right) d^3\boldsymbol{v} where alpha is the species We find from moments of the Boltzmann equation: Gamma = A_{alpha}n_{alpha} vec{v}_{alpha}vec{v}_{alpha} + P_{alpha} or \Gamma = A_{\alpha}n_{\alpha} \boldsymbol{v}_{\alpha}\boldsymbol{v}_{\alpha} + P_{\alpha}.
More...
Computes the Full Dyadic Second Moment which becomes the Fluid Momentum Flux Gamma = Full velocity space integral of vec{v}(A_{alpha} vec{v} f_{alpha}) d^{3}vec{v} or \Gamma = \iiint \boldsymbol{v} \left(A_{\alpha} \boldsymbol{v} f_{\alpha}\right) d^3\boldsymbol{v} where alpha is the species We find from moments of the Boltzmann equation: Gamma = A_{alpha}n_{alpha} vec{v}_{alpha}vec{v}_{alpha} + P_{alpha} or \Gamma = A_{\alpha}n_{\alpha} \boldsymbol{v}_{\alpha}\boldsymbol{v}_{\alpha} + P_{\alpha}.
This moment is useful for using as a momentum flux boundary condition for the fluid (5N-Moment) from the Kinetic Model
Note: This becomes a symmetric tensor so instead of all 9 components, the unique 6 components are calculated: \iiint v_{x} \left(A_{\alpha} v_{x} f_{\alpha}\right) d^3\boldsymbol{v} \iiint v_{x} \left(A_{\alpha} v_{y} f_{\alpha}\right) d^3\boldsymbol{v} \iiint v_{x} \left(A_{\alpha} v_{z} f_{\alpha}\right) d^3\boldsymbol{v} \iiint v_{y} \left(A_{\alpha} v_{y} f_{\alpha}\right) d^3\boldsymbol{v} \iiint v_{y} \left(A_{\alpha} v_{z} f_{\alpha}\right) d^3\boldsymbol{v} \iiint v_{z} \left(A_{\alpha} v_{z} f_{\alpha}\right) d^3\boldsymbol{v}
|
| FullDyadicSecondMoment () |
| Constructor.
|
|
| ~FullDyadicSecondMoment () override |
| Destructor.
|
|
void | setup (const WxCryptSet &wxc) override |
| Setup.
|
|
const std::vector< int > & | getInputVariableIndexes (int flag) const override |
| Grab input variable indices.
|
|
const std::vector< int > & | getOutputVariableIndexes (int flag) const override |
| Grab output variable indices.
|
|
void | evaluate_function (const real *q, const real *aux, const solverVariables_t *pSV, real *result) const override |
| Redefinition of the function evaluator.
|
|
| WmApplication () |
|
virtual | ~WmApplication ()=default |
|
virtual void | setup (const WxCryptSet &wxc) |
|
virtual const std::vector< int > & | getAuxiliaryVariableIndexes (int flag=WMAPPLICATIONFLAG_NONE) const |
|
virtual const std::vector< int > & | getCrossVariableIndexes (int flag=WMAPPLICATIONFLAG_NONE) const |
|
const std::vector< std::string > & | getBoundaryNames () const |
|
bool | isOnBoundary (const std::string &boundaryName) const |
|
bool | has (int flag) const |
|
virtual real | numerical_flux (const real *q_l, const real *q_r, const real *aux_l, const real *aux_r, const solverVariables_t *pFV, real *numericalFlux) const |
|
virtual real | internal_flux (const real *q, const real *aux, const solverVariables_t *pSV, std::vector< std::vector< real > > &internalFlux) const |
|
virtual real | source (const real *q, const real *aux, const elementGeometry_t *pEG, real *source) const |
|
virtual void | bc_q (const real *q_in, const real *aux_in, const real *aux_out, const solverVariables_t *pFV, real *q_out) const |
| Boundary Condition Application which sets the boundary condition on ghost nodes.
|
|
virtual real | bcNumericalFlux (const real *q_l, const real *q_r, const real *aux_l, const real *aux_r, const solverVariables_t *pFV, real *numericalFlux) const |
|
virtual void | conserved_to_primitive (const real *q, const real *aux, real *w) const |
|
virtual void | primitive_to_conserved (const real *w, const real *aux, real *q) const |
|
virtual void | evaluate_function (const real *q, const real *aux, const solverVariables_t *pSV, real *result) const |
|
virtual void | bc_q_kinetic (const real *q_in, const real *aux_in, const solverVariables_t *pFV, real *q_out) const |
|
std::shared_ptr< std::string > | app_name () |
|
virtual const std::vector< int > & | getInputVariableIndexes (int flag=0) const |
|
virtual const std::vector< int > & | getOutputVariableIndexes (int flag=0) const |
|