26 #include "grins_config.h"
28 #ifdef GRINS_HAVE_ANTIOCH
37 #include "antioch/temp_cache.h"
38 #include "antioch/vector_utils.h"
43 : _antioch_mixture( mixture ),
44 _antioch_kinetics( mixture.reaction_set(), 0 ),
45 _antioch_cea_thermo( mixture.cea_mixture() )
49 const libMesh::Real rho,
50 const std::vector<libMesh::Real>& mass_fractions,
51 std::vector<libMesh::Real>& omega_dot )
53 Antioch::TempCache<libMesh::Real> temp_cache(T);
54 this->
omega_dot(temp_cache,rho,mass_fractions,omega_dot);
58 const libMesh::Real rho,
59 const std::vector<libMesh::Real>& mass_fractions,
60 std::vector<libMesh::Real>& omega_dot )
64 libmesh_assert_equal_to( mass_fractions.size(), n_species );
65 libmesh_assert_equal_to( omega_dot.size(), n_species );
67 std::vector<libMesh::Real> h_RT_minus_s_R(n_species, 0.0);
68 std::vector<libMesh::Real> molar_densities(n_species, 0.0);
82 #endif // GRINS_HAVE_ANTIOCH
Antioch::KineticsEvaluator< libMesh::Real > _antioch_kinetics
unsigned int n_species() const
void molar_densities(const libMesh::Real rho, const std::vector< libMesh::Real > &mass_fractions, std::vector< libMesh::Real > &molar_densities) const
Molar density for all species, [mol/m^3].
Antioch::CEAEvaluator< libMesh::Real > _antioch_cea_thermo
void omega_dot(const libMesh::Real &T, const libMesh::Real rho, const std::vector< libMesh::Real > &mass_fractions, std::vector< libMesh::Real > &omega_dot)
Wrapper class for storing state for Antioch thermo and kinetics.
const AntiochMixture & _antioch_mixture