26 #include "grins_config.h"
28 #ifdef GRINS_HAVE_ANTIOCH
38 #include "libmesh/getpot.h"
39 #include "libmesh/parameter_multiaccessor.h"
42 #include "antioch/read_reaction_set_data.h"
43 #include "antioch/cea_mixture_ascii_parsing.h"
44 #include "antioch/stat_mech_thermo.h"
49 const std::string& material )
51 _reaction_set( new Antioch::ReactionSet<
libMesh::Real>( (*_antioch_gas.get()) ) ),
52 _cea_mixture( new Antioch::CEAThermoMixture<
libMesh::Real>( (*_antioch_gas.get()) ) )
56 bool verbose_read = input(
"screen-options/verbose_kinetics_read",
false );
58 Antioch::read_reaction_set_data_xml<libMesh::Real>( kinetics_data_filename, verbose_read, *
_reaction_set.get() );
60 std::string cea_data_filename = input(
"Materials/"+material+
"/GasMixture/Antioch/cea_data",
"default" );
61 if( cea_data_filename == std::string(
"default") )
62 cea_data_filename = Antioch::DefaultInstallFilename::thermo_data();
64 Antioch::read_cea_mixture_data_ascii( *
_cea_mixture.get(), cea_data_filename );
70 (
const std::string & param_name,
78 if (param_name.find(
"Antioch") == 0)
79 param_pointer.push_back
81 (*this->_reaction_set.get(), param_name));
87 Antioch::StatMechThermodynamics<libMesh::Real> thermo( *(this->
_antioch_gas.get()) );
91 for(
unsigned int s = 0; s < this->
n_species(); s++ )
99 #endif // GRINS_HAVE_ANTIOCH
libMesh::UniquePtr< Antioch::ChemicalMixture< libMesh::Real > > _antioch_gas
unsigned int n_species() const
libMesh::UniquePtr< Antioch::ReactionSet< libMesh::Real > > _reaction_set
Wrapper class for Antioch::ChemicalMixture.
void build_stat_mech_ref_correction()
virtual void register_parameter(const std::string ¶m_name, libMesh::ParameterMultiAccessor< libMesh::Number > ¶m_pointer) const
Each subclass will register its copy of an independent.
libMesh::UniquePtr< Antioch::CEAThermoMixture< libMesh::Real > > _cea_mixture
std::vector< libMesh::Real > _h_stat_mech_ref_correction
static std::string parse_chemical_kinetics_datafile_name(const GetPot &input, const std::string &material)
virtual void register_parameter(const std::string ¶m_name, libMesh::ParameterMultiAccessor< libMesh::Number > ¶m_pointer) const
Each subclass will register its copy of an independent.