26 #ifndef GRINS_REACTING_LOW_MACH_NAVIER_STOKES_ABSTRACT_H
27 #define GRINS_REACTING_LOW_MACH_NAVIER_STOKES_ABSTRACT_H
30 #include "grins_config.h"
64 std::vector<libMesh::Real>& mass_fracs )
const;
66 libMesh::Real
rho( libMesh::Real
T, libMesh::Real p0, libMesh::Real R_mix)
const;
87 libMesh::UniquePtr<ThermoPressureFEVariable>
_p0_var;
128 std::vector<libMesh::Real>& mass_fracs )
const
130 libmesh_assert_equal_to(mass_fracs.size(), this->
_n_species);
132 for(
unsigned int var = 0; var < this->
_n_species; var++ )
141 libMesh::Real R_mix)
const
143 libMesh::Real value = 0;
147 value = p0/(R_mix*
T);
154 unsigned int qp )
const
159 p0 = c.interior_value(
_p0_var->p0(), qp );
170 unsigned int qp )
const
175 p0 = c.side_value(
_p0_var->p0(), qp );
186 const libMesh::Point& p )
const
191 p0 = c.point_value(
_p0_var->p0(), p );
202 unsigned int qp )
const
207 p0 = c.fixed_interior_value(
_p0_var->p0(), qp );
218 #endif //GRINS_REACTING_LOW_MACH_NAVIER_STOKES_ABSTRACT_H
unsigned int n_species() const
libMesh::Real get_p0_steady_side(const AssemblyContext &c, unsigned int qp) const
bool _enable_thermo_press_calc
Flag to enable thermodynamic pressure calculation.
PrimitiveTempFEVariables _temp_vars
VelocityFEVariables _flow_vars
Physics abstract base class. Defines API for physics to be added to MultiphysicsSystem.
unsigned int _n_species
Number of species.
libMesh::Real T(const libMesh::Point &p, const AssemblyContext &c) const
void mass_fractions(const libMesh::Point &p, const AssemblyContext &c, std::vector< libMesh::Real > &mass_fracs) const
SpeciesMassFractionsFEVariables _species_vars
unsigned int _dim
Physical dimension of problem.
libMesh::UniquePtr< ThermoPressureFEVariable > _p0_var
libMesh::Real rho(libMesh::Real T, libMesh::Real p0, libMesh::Real R_mix) const
void register_variables()
virtual void init_context(AssemblyContext &context)
Initialize context for added physics variables.
void read_input_options(const GetPot &input)
Read options from GetPot input file.
virtual ~ReactingLowMachNavierStokesAbstract()
PressureFEVariable _press_var
virtual void set_time_evolving_vars(libMesh::FEMSystem *system)
Sets velocity variables to be time-evolving.
libMesh::Real get_p0_transient(const AssemblyContext &c, unsigned int qp) const
VariableIndex species(unsigned int species) const
virtual void init_variables(libMesh::FEMSystem *system)
Initialize variables for this physics.
libMesh::Real get_p0_steady(const AssemblyContext &c, unsigned int qp) const
libMesh::Real _fixed_rho_value
libMesh::Point _g
Gravity vector.
ReactingLowMachNavierStokesAbstract()