Hydrogen (H2) fluid properties as a function of pressure (Pa) and temperature (K). More...

#include <HydrogenFluidProperties.h>

Inheritance diagram for HydrogenFluidProperties:

Public Types
typedef DataFileName	DataFileParameterType

Public Member Functions
	HydrogenFluidProperties (const InputParameters &parameters)

virtual std::string	fluidName () const override
	Fluid name. More...

virtual Real	molarMass () const override
	Molar mass [kg/mol]. More...

virtual Real	mu_from_rho_T (Real density, Real temperature) const override

void	mu_from_rho_T (Real density, Real temperature, Real ddensity_dT, Real &mu, Real &dmu_drho, Real &dmu_dT) const

virtual Real	mu_from_p_T (Real pressure, Real temperature) const override

virtual void	mu_from_p_T (Real pressure, Real temperature, Real &mu, Real &dmu_dp, Real &dmu_dT) const override

virtual void	rho_mu_from_p_T (Real pressure, Real temperature, Real &rho, Real &mu) const override
	Combined methods. More...

virtual void	rho_mu_from_p_T (Real pressure, Real temperature, Real &rho, Real &drho_dp, Real &drho_dT, Real &mu, Real &dmu_dp, Real &dmu_dT) const override

virtual Real	k_from_rho_T (Real density, Real temperature) const override

virtual Real	k_from_p_T (Real pressure, Real temperature) const override

virtual void	k_from_p_T (Real pressure, Real temperature, Real &k, Real &dk_dp, Real &dk_dT) const override

virtual std::vector< Real >	henryCoefficients () const override
	Henry's law coefficients for dissolution in water. More...

virtual Real	criticalPressure () const override
	Critical pressure. More...

virtual Real	criticalTemperature () const override
	Critical temperature. More...

virtual Real	criticalDensity () const override
	Critical density. More...

virtual Real	triplePointPressure () const override
	Triple point pressure. More...

virtual Real	triplePointTemperature () const override
	Triple point temperature. More...

virtual Real	vaporPressure (Real temperature) const override
	Vapor pressure. More...

virtual void	vaporPressure (Real temperature, Real &psat, Real &dpsat_dT) const override

virtual Real	rho_from_p_T (Real pressure, Real temperature) const override

virtual void	rho_from_p_T (Real pressure, Real temperature, Real &rho, Real &drho_dp, Real &drho_dT) const override

virtual Real	e_from_p_T (Real pressure, Real temperature) const override

virtual void	e_from_p_T (Real p, Real T, Real &e, Real &de_dp, Real &de_dT) const override

virtual Real	c_from_p_T (Real pressure, Real temperature) const override

virtual Real	cp_from_p_T (Real pressure, Real temperature) const override

virtual Real	cv_from_p_T (Real pressure, Real temperature) const override

virtual Real	s_from_p_T (Real pressure, Real temperature) const override

virtual void	s_from_p_T (Real p, Real T, Real &s, Real &ds_dp, Real &ds_dT) const override

virtual Real	h_from_p_T (Real pressure, Real temperature) const override

virtual void	h_from_p_T (Real p, Real T, Real &h, Real &dh_dp, Real &dh_dT) const override

virtual Real	T_from_p_h (Real pressure, Real enthalpy) const override

virtual Real	p_from_rho_T (Real rho, Real T) const
	Pressure as a function of density and temperature. More...

virtual Real	criticalInternalEnergy () const
	Critical specific internal energy. More...

virtual Real	e_spndl_from_v (Real v) const
	Specific internal energy from temperature and specific volume. More...

virtual void	v_e_spndl_from_T (Real T, Real &v, Real &e) const
	Specific internal energy from temperature and specific volume. More...

virtual ADReal	vaporPressure (const ADReal &T) const

virtual Real	vaporTemperature (Real p) const
	Vapor temperature. More...

virtual void	vaporTemperature (Real p, Real &Tsat, Real &dTsat_dp) const

virtual ADReal	vaporTemperature (const ADReal &p) const

template<typename CppType >
void	v_e_from_p_T (const CppType &p, const CppType &T, CppType &v, CppType &e) const

template<typename CppType >
void	v_e_from_p_T (const CppType &p, const CppType &T, CppType &v, CppType &dv_dp, CppType &dv_dT, CppType &e, CppType &de_dp, CppType &de_dT) const

virtual void	rho_mu_from_p_T (const ADReal &p, const ADReal &T, ADReal &rho, ADReal &mu) const

virtual void	rho_e_from_p_T (Real p, Real T, Real &rho, Real &drho_dp, Real &drho_dT, Real &e, Real &de_dp, Real &de_dT) const

template<typename CppType >
void	p_T_from_v_e (const CppType &v, const CppType &e, Real p0, Real T0, CppType &p, CppType &T, bool &conversion_succeeded) const
	Determines (p,T) from (v,e) using Newton Solve in 2D Useful for conversion between different sets of state variables. More...

template<typename T >
void	p_T_from_v_h (const T &v, const T &h, Real p0, Real T0, T &pressure, T &temperature, bool &conversion_succeeded) const
	Determines (p,T) from (v,h) using Newton Solve in 2D Useful for conversion between different sets of state variables. More...

template<typename T >
void	p_T_from_h_s (const T &h, const T &s, Real p0, Real T0, T &pressure, T &temperature, bool &conversion_succeeded) const
	Determines (p,T) from (h,s) using Newton Solve in 2D Useful for conversion between different sets of state variables. More...

template<>
std::pair< Real, Real >	makeZeroAndOne (const Real &)

virtual void	execute () final

virtual void	initialize () final

virtual void	finalize () final

virtual void	threadJoin (const UserObject &) final

virtual void	subdomainSetup () final

bool	needThreadedCopy () const override final

SubProblem &	getSubProblem () const

bool	shouldDuplicateInitialExecution () const

virtual Real	spatialValue (const Point &) const

virtual const std::vector< Point >	spatialPoints () const

void	gatherSum (T &value)

void	gatherMax (T &value)

void	gatherMin (T &value)

void	gatherProxyValueMax (T1 &proxy, T2 &value)

void	gatherProxyValueMin (T1 &proxy, T2 &value)

void	setPrimaryThreadCopy (UserObject *primary)

UserObject *	primaryThreadCopy ()

std::set< UserObjectName >	getDependObjects () const

const std::set< std::string > &	getRequestedItems () override

const std::set< std::string > &	getSuppliedItems () override

unsigned int	systemNumber () const

virtual bool	enabled () const

std::shared_ptr< MooseObject >	getSharedPtr ()

std::shared_ptr< const MooseObject >	getSharedPtr () const

MooseApp &	getMooseApp () const

const std::string &	type () const

virtual const std::string &	name () const

std::string	typeAndName () const

std::string	errorPrefix (const std::string &error_type) const

void	callMooseError (std::string msg, const bool with_prefix) const

MooseObjectParameterName	uniqueParameterName (const std::string &parameter_name) const

const InputParameters &	parameters () const

MooseObjectName	uniqueName () const

const T &	getParam (const std::string &name) const

std::vector< std::pair< T1, T2 > >	getParam (const std::string &param1, const std::string &param2) const

const T *	queryParam (const std::string &name) const

const T &	getRenamedParam (const std::string &old_name, const std::string &new_name) const

T	getCheckedPointerParam (const std::string &name, const std::string &error_string="") const

bool	isParamValid (const std::string &name) const

bool	isParamSetByUser (const std::string &nm) const

void	paramError (const std::string &param, Args... args) const

void	paramWarning (const std::string &param, Args... args) const

void	paramInfo (const std::string &param, Args... args) const

void	connectControllableParams (const std::string &parameter, const std::string &object_type, const std::string &object_name, const std::string &object_parameter) const

void	mooseError (Args &&... args) const

void	mooseErrorNonPrefixed (Args &&... args) const

void	mooseDocumentedError (const std::string &repo_name, const unsigned int issue_num, Args &&... args) const

void	mooseWarning (Args &&... args) const

void	mooseWarningNonPrefixed (Args &&... args) const

void	mooseDeprecated (Args &&... args) const

void	mooseInfo (Args &&... args) const

std::string	getDataFileName (const std::string &param) const

std::string	getDataFileNameByName (const std::string &relative_path) const

std::string	getDataFilePath (const std::string &relative_path) const

virtual void	initialSetup ()

virtual void	timestepSetup ()

virtual void	jacobianSetup ()

virtual void	residualSetup ()

virtual void	customSetup (const ExecFlagType &)

const ExecFlagEnum &	getExecuteOnEnum () const

UserObjectName	getUserObjectName (const std::string &param_name) const

const T &	getUserObject (const std::string &param_name, bool is_dependency=true) const

const T &	getUserObjectByName (const UserObjectName &object_name, bool is_dependency=true) const

const UserObject &	getUserObjectBase (const std::string &param_name, bool is_dependency=true) const

const UserObject &	getUserObjectBaseByName (const UserObjectName &object_name, bool is_dependency=true) const

const std::vector< MooseVariableScalar *> &	getCoupledMooseScalarVars ()

const std::set< TagID > &	getScalarVariableCoupleableVectorTags () const

const std::set< TagID > &	getScalarVariableCoupleableMatrixTags () const

const GenericMaterialProperty< T, is_ad > &	getGenericMaterialProperty (const std::string &name, MaterialData &material_data, const unsigned int state=0)

const GenericMaterialProperty< T, is_ad > &	getGenericMaterialProperty (const std::string &name, const unsigned int state=0)

const GenericMaterialProperty< T, is_ad > &	getGenericMaterialProperty (const std::string &name, const unsigned int state=0)

const MaterialProperty< T > &	getMaterialProperty (const std::string &name, MaterialData &material_data, const unsigned int state=0)

const MaterialProperty< T > &	getMaterialProperty (const std::string &name, const unsigned int state=0)

const MaterialProperty< T > &	getMaterialProperty (const std::string &name, const unsigned int state=0)

const ADMaterialProperty< T > &	getADMaterialProperty (const std::string &name, MaterialData &material_data)

const ADMaterialProperty< T > &	getADMaterialProperty (const std::string &name)

const ADMaterialProperty< T > &	getADMaterialProperty (const std::string &name)

const MaterialProperty< T > &	getMaterialPropertyOld (const std::string &name, MaterialData &material_data)

const MaterialProperty< T > &	getMaterialPropertyOld (const std::string &name)

const MaterialProperty< T > &	getMaterialPropertyOld (const std::string &name)

const MaterialProperty< T > &	getMaterialPropertyOlder (const std::string &name, MaterialData &material_data)

const MaterialProperty< T > &	getMaterialPropertyOlder (const std::string &name)

const MaterialProperty< T > &	getMaterialPropertyOlder (const std::string &name)

const GenericMaterialProperty< T, is_ad > &	getGenericMaterialPropertyByName (const MaterialPropertyName &name, MaterialData &material_data, const unsigned int state)

const GenericMaterialProperty< T, is_ad > &	getGenericMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)

const GenericMaterialProperty< T, is_ad > &	getGenericMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)

const MaterialProperty< T > &	getMaterialPropertyByName (const MaterialPropertyName &name, MaterialData &material_data, const unsigned int state=0)

const MaterialProperty< T > &	getMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)

const MaterialProperty< T > &	getMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)

const ADMaterialProperty< T > &	getADMaterialPropertyByName (const MaterialPropertyName &name, MaterialData &material_data)

const ADMaterialProperty< T > &	getADMaterialPropertyByName (const MaterialPropertyName &name)

const ADMaterialProperty< T > &	getADMaterialPropertyByName (const MaterialPropertyName &name)

const MaterialProperty< T > &	getMaterialPropertyOldByName (const MaterialPropertyName &name, MaterialData &material_data)

const MaterialProperty< T > &	getMaterialPropertyOldByName (const MaterialPropertyName &name)

const MaterialProperty< T > &	getMaterialPropertyOldByName (const MaterialPropertyName &name)

const MaterialProperty< T > &	getMaterialPropertyOlderByName (const MaterialPropertyName &name, MaterialData &material_data)

const MaterialProperty< T > &	getMaterialPropertyOlderByName (const MaterialPropertyName &name)

const MaterialProperty< T > &	getMaterialPropertyOlderByName (const MaterialPropertyName &name)

std::pair< const MaterialProperty< T > *, std::set< SubdomainID > >	getBlockMaterialProperty (const MaterialPropertyName &name)

const GenericMaterialProperty< T, is_ad > &	getGenericZeroMaterialProperty (const std::string &name)

const GenericMaterialProperty< T, is_ad > &	getGenericZeroMaterialProperty ()

const GenericMaterialProperty< T, is_ad > &	getGenericZeroMaterialPropertyByName (const std::string &prop_name)

const MaterialProperty< T > &	getZeroMaterialProperty (Ts... args)

std::set< SubdomainID >	getMaterialPropertyBlocks (const std::string &name)

std::vector< SubdomainName >	getMaterialPropertyBlockNames (const std::string &name)

std::set< BoundaryID >	getMaterialPropertyBoundaryIDs (const std::string &name)

std::vector< BoundaryName >	getMaterialPropertyBoundaryNames (const std::string &name)

void	checkBlockAndBoundaryCompatibility (std::shared_ptr< MaterialBase > discrete)

std::unordered_map< SubdomainID, std::vector< MaterialBase *> >	buildRequiredMaterials (bool allow_stateful=true)

void	statefulPropertiesAllowed (bool)

bool	getMaterialPropertyCalled () const

virtual const std::unordered_set< unsigned int > &	getMatPropDependencies () const

virtual void	resolveOptionalProperties ()

const GenericMaterialProperty< T, is_ad > &	getPossiblyConstantGenericMaterialPropertyByName (const MaterialPropertyName &prop_name, MaterialData &material_data, const unsigned int state)

bool	isImplicit ()

Moose::StateArg	determineState () const

bool	hasUserObject (const std::string &param_name) const

bool	hasUserObject (const std::string &param_name) const

bool	hasUserObject (const std::string &param_name) const

bool	hasUserObject (const std::string &param_name) const

bool	hasUserObjectByName (const UserObjectName &object_name) const

bool	hasUserObjectByName (const UserObjectName &object_name) const

bool	hasUserObjectByName (const UserObjectName &object_name) const

bool	hasUserObjectByName (const UserObjectName &object_name) const

const GenericOptionalMaterialProperty< T, is_ad > &	getGenericOptionalMaterialProperty (const std::string &name, const unsigned int state=0)

const GenericOptionalMaterialProperty< T, is_ad > &	getGenericOptionalMaterialProperty (const std::string &name, const unsigned int state=0)

const OptionalMaterialProperty< T > &	getOptionalMaterialProperty (const std::string &name, const unsigned int state=0)

const OptionalMaterialProperty< T > &	getOptionalMaterialProperty (const std::string &name, const unsigned int state=0)

const OptionalADMaterialProperty< T > &	getOptionalADMaterialProperty (const std::string &name)

const OptionalADMaterialProperty< T > &	getOptionalADMaterialProperty (const std::string &name)

const OptionalMaterialProperty< T > &	getOptionalMaterialPropertyOld (const std::string &name)

const OptionalMaterialProperty< T > &	getOptionalMaterialPropertyOld (const std::string &name)

const OptionalMaterialProperty< T > &	getOptionalMaterialPropertyOlder (const std::string &name)

const OptionalMaterialProperty< T > &	getOptionalMaterialPropertyOlder (const std::string &name)

MaterialBase &	getMaterial (const std::string &name)

MaterialBase &	getMaterial (const std::string &name)

MaterialBase &	getMaterialByName (const std::string &name, bool no_warn=false)

MaterialBase &	getMaterialByName (const std::string &name, bool no_warn=false)

bool	hasMaterialProperty (const std::string &name)

bool	hasMaterialProperty (const std::string &name)

bool	hasMaterialPropertyByName (const std::string &name)

bool	hasMaterialPropertyByName (const std::string &name)

bool	hasADMaterialProperty (const std::string &name)

bool	hasADMaterialProperty (const std::string &name)

bool	hasADMaterialPropertyByName (const std::string &name)

bool	hasADMaterialPropertyByName (const std::string &name)

bool	hasGenericMaterialProperty (const std::string &name)

bool	hasGenericMaterialProperty (const std::string &name)

bool	hasGenericMaterialPropertyByName (const std::string &name)

bool	hasGenericMaterialPropertyByName (const std::string &name)

const Function &	getFunction (const std::string &name) const

const Function &	getFunctionByName (const FunctionName &name) const

bool	hasFunction (const std::string &param_name) const

bool	hasFunctionByName (const FunctionName &name) const

bool	isDefaultPostprocessorValue (const std::string &param_name, const unsigned int index=0) const

bool	hasPostprocessor (const std::string &param_name, const unsigned int index=0) const

bool	hasPostprocessorByName (const PostprocessorName &name) const

std::size_t	coupledPostprocessors (const std::string &param_name) const

const PostprocessorName &	getPostprocessorName (const std::string &param_name, const unsigned int index=0) const

const VectorPostprocessorValue &	getVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name) const

const VectorPostprocessorValue &	getVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name, bool needs_broadcast) const

const VectorPostprocessorValue &	getVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name) const

const VectorPostprocessorValue &	getVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name, bool needs_broadcast) const

const VectorPostprocessorValue &	getVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name) const

const VectorPostprocessorValue &	getVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name, bool needs_broadcast) const

const VectorPostprocessorValue &	getVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name) const

const VectorPostprocessorValue &	getVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name, bool needs_broadcast) const

const ScatterVectorPostprocessorValue &	getScatterVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name) const

const ScatterVectorPostprocessorValue &	getScatterVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name) const

const ScatterVectorPostprocessorValue &	getScatterVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name) const

const ScatterVectorPostprocessorValue &	getScatterVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name) const

bool	hasVectorPostprocessor (const std::string &param_name, const std::string &vector_name) const

bool	hasVectorPostprocessor (const std::string &param_name) const

bool	hasVectorPostprocessorByName (const VectorPostprocessorName &name, const std::string &vector_name) const

bool	hasVectorPostprocessorByName (const VectorPostprocessorName &name) const

const VectorPostprocessorName &	getVectorPostprocessorName (const std::string &param_name) const

T &	getSampler (const std::string &name)

Sampler &	getSampler (const std::string &name)

T &	getSamplerByName (const SamplerName &name)

Sampler &	getSamplerByName (const SamplerName &name)

virtual void	meshChanged ()

virtual void	meshDisplaced ()

PerfGraph &	perfGraph ()

const PostprocessorValue &	getPostprocessorValue (const std::string &param_name, const unsigned int index=0) const

const PostprocessorValue &	getPostprocessorValue (const std::string &param_name, const unsigned int index=0) const

const PostprocessorValue &	getPostprocessorValueOld (const std::string &param_name, const unsigned int index=0) const

const PostprocessorValue &	getPostprocessorValueOld (const std::string &param_name, const unsigned int index=0) const

const PostprocessorValue &	getPostprocessorValueOlder (const std::string &param_name, const unsigned int index=0) const

const PostprocessorValue &	getPostprocessorValueOlder (const std::string &param_name, const unsigned int index=0) const

virtual const PostprocessorValue &	getPostprocessorValueByName (const PostprocessorName &name) const

virtual const PostprocessorValue &	getPostprocessorValueByName (const PostprocessorName &name) const

const PostprocessorValue &	getPostprocessorValueOldByName (const PostprocessorName &name) const

const PostprocessorValue &	getPostprocessorValueOldByName (const PostprocessorName &name) const

const PostprocessorValue &	getPostprocessorValueOlderByName (const PostprocessorName &name) const

const PostprocessorValue &	getPostprocessorValueOlderByName (const PostprocessorName &name) const

bool	isVectorPostprocessorDistributed (const std::string &param_name) const

bool	isVectorPostprocessorDistributed (const std::string &param_name) const

bool	isVectorPostprocessorDistributedByName (const VectorPostprocessorName &name) const

bool	isVectorPostprocessorDistributedByName (const VectorPostprocessorName &name) const

const Distribution &	getDistribution (const std::string &name) const

const T &	getDistribution (const std::string &name) const

const Distribution &	getDistribution (const std::string &name) const

const T &	getDistribution (const std::string &name) const

const Distribution &	getDistributionByName (const DistributionName &name) const

const T &	getDistributionByName (const std::string &name) const

const Distribution &	getDistributionByName (const DistributionName &name) const

const T &	getDistributionByName (const std::string &name) const

const Parallel::Communicator &	comm () const

processor_id_type	n_processors () const

processor_id_type	processor_id () const

Static Public Member Functions
static InputParameters	validParams ()

static void	sort (typename std::vector< T > &vector)

static void	sortDFS (typename std::vector< T > &vector)

static void	cyclicDependencyError (CyclicDependencyException< T2 > &e, const std::string &header)

Public Attributes
const ConsoleStream	_console

Static Public Attributes
static const Real	_R = 8.3144598
	Universal gas constant (J/mol/K) More...

static constexpr PropertyValue::id_type	default_property_id

static constexpr PropertyValue::id_type	zero_property_id

static constexpr auto	SYSTEM

static constexpr auto	NAME

Protected Member Functions
virtual Real	alpha (Real delta, Real tau) const override
	Helmholtz free energy for H2 From Leachman et al (reference above) More...

virtual Real	dalpha_ddelta (Real delta, Real tau) const override
	Derivative of Helmholtz free energy wrt delta. More...

virtual Real	dalpha_dtau (Real delta, Real tau) const override
	Derivative of Helmholtz free energy wrt tau. More...

virtual Real	d2alpha_ddelta2 (Real delta, Real tau) const override
	Second derivative of Helmholtz free energy wrt delta. More...

virtual Real	d2alpha_dtau2 (Real delta, Real tau) const override
	Second derivative of Helmholtz free energy wrt tau. More...

virtual Real	d2alpha_ddeltatau (Real delta, Real tau) const override
	Second derivative of Helmholtz free energy wrt delta and tau. More...

virtual void	addPostprocessorDependencyHelper (const PostprocessorName &name) const override

virtual void	addVectorPostprocessorDependencyHelper (const VectorPostprocessorName &name) const override

virtual void	addUserObjectDependencyHelper (const UserObject &uo) const override

void	addReporterDependencyHelper (const ReporterName &reporter_name) override

const ReporterName &	getReporterName (const std::string &param_name) const

T &	declareRestartableData (const std::string &data_name, Args &&... args)

ManagedValue< T >	declareManagedRestartableDataWithContext (const std::string &data_name, void *context, Args &&... args)

const T &	getRestartableData (const std::string &data_name) const

T &	declareRestartableDataWithContext (const std::string &data_name, void *context, Args &&... args)

T &	declareRecoverableData (const std::string &data_name, Args &&... args)

T &	declareRestartableDataWithObjectName (const std::string &data_name, const std::string &object_name, Args &&... args)

T &	declareRestartableDataWithObjectNameWithContext (const std::string &data_name, const std::string &object_name, void *context, Args &&... args)

std::string	restartableName (const std::string &data_name) const

const T &	getMeshProperty (const std::string &data_name, const std::string &prefix)

const T &	getMeshProperty (const std::string &data_name)

bool	hasMeshProperty (const std::string &data_name, const std::string &prefix) const

bool	hasMeshProperty (const std::string &data_name, const std::string &prefix) const

bool	hasMeshProperty (const std::string &data_name) const

bool	hasMeshProperty (const std::string &data_name) const

std::string	meshPropertyName (const std::string &data_name) const

PerfID	registerTimedSection (const std::string &section_name, const unsigned int level) const

PerfID	registerTimedSection (const std::string &section_name, const unsigned int level, const std::string &live_message, const bool print_dots=true) const

std::string	timedSectionName (const std::string &section_name) const

bool	isCoupledScalar (const std::string &var_name, unsigned int i=0) const

unsigned int	coupledScalarComponents (const std::string &var_name) const

unsigned int	coupledScalar (const std::string &var_name, unsigned int comp=0) const

libMesh::Order	coupledScalarOrder (const std::string &var_name, unsigned int comp=0) const

const VariableValue &	coupledScalarValue (const std::string &var_name, unsigned int comp=0) const

const ADVariableValue &	adCoupledScalarValue (const std::string &var_name, unsigned int comp=0) const

const GenericVariableValue< is_ad > &	coupledGenericScalarValue (const std::string &var_name, unsigned int comp=0) const

const GenericVariableValue< false > &	coupledGenericScalarValue (const std::string &var_name, const unsigned int comp) const

const GenericVariableValue< true > &	coupledGenericScalarValue (const std::string &var_name, const unsigned int comp) const

const VariableValue &	coupledVectorTagScalarValue (const std::string &var_name, TagID tag, unsigned int comp=0) const

const VariableValue &	coupledMatrixTagScalarValue (const std::string &var_name, TagID tag, unsigned int comp=0) const

const VariableValue &	coupledScalarValueOld (const std::string &var_name, unsigned int comp=0) const

const VariableValue &	coupledScalarValueOlder (const std::string &var_name, unsigned int comp=0) const

const VariableValue &	coupledScalarDot (const std::string &var_name, unsigned int comp=0) const

const ADVariableValue &	adCoupledScalarDot (const std::string &var_name, unsigned int comp=0) const

const VariableValue &	coupledScalarDotDot (const std::string &var_name, unsigned int comp=0) const

const VariableValue &	coupledScalarDotOld (const std::string &var_name, unsigned int comp=0) const

const VariableValue &	coupledScalarDotDotOld (const std::string &var_name, unsigned int comp=0) const

const VariableValue &	coupledScalarDotDu (const std::string &var_name, unsigned int comp=0) const

const VariableValue &	coupledScalarDotDotDu (const std::string &var_name, unsigned int comp=0) const

const MooseVariableScalar *	getScalarVar (const std::string &var_name, unsigned int comp) const

virtual void	checkMaterialProperty (const std::string &name, const unsigned int state)

void	markMatPropRequested (const std::string &)

MaterialPropertyName	getMaterialPropertyName (const std::string &name) const

void	checkExecutionStage ()

const T &	getReporterValue (const std::string &param_name, const std::size_t time_index=0)

const T &	getReporterValue (const std::string &param_name, ReporterMode mode, const std::size_t time_index=0)

const T &	getReporterValue (const std::string &param_name, const std::size_t time_index=0)

const T &	getReporterValue (const std::string &param_name, ReporterMode mode, const std::size_t time_index=0)

const T &	getReporterValueByName (const ReporterName &reporter_name, const std::size_t time_index=0)

const T &	getReporterValueByName (const ReporterName &reporter_name, ReporterMode mode, const std::size_t time_index=0)

const T &	getReporterValueByName (const ReporterName &reporter_name, const std::size_t time_index=0)

const T &	getReporterValueByName (const ReporterName &reporter_name, ReporterMode mode, const std::size_t time_index=0)

bool	hasReporterValue (const std::string &param_name) const

bool	hasReporterValue (const std::string &param_name) const

bool	hasReporterValue (const std::string &param_name) const

bool	hasReporterValue (const std::string &param_name) const

bool	hasReporterValueByName (const ReporterName &reporter_name) const

bool	hasReporterValueByName (const ReporterName &reporter_name) const

bool	hasReporterValueByName (const ReporterName &reporter_name) const

bool	hasReporterValueByName (const ReporterName &reporter_name) const

const GenericMaterialProperty< T, is_ad > *	defaultGenericMaterialProperty (const std::string &name)

const GenericMaterialProperty< T, is_ad > *	defaultGenericMaterialProperty (const std::string &name)

const MaterialProperty< T > *	defaultMaterialProperty (const std::string &name)

const MaterialProperty< T > *	defaultMaterialProperty (const std::string &name)

const ADMaterialProperty< T > *	defaultADMaterialProperty (const std::string &name)

const ADMaterialProperty< T > *	defaultADMaterialProperty (const std::string &name)

void	flagInvalidSolutionInternal (const InvalidSolutionID invalid_solution_id) const

InvalidSolutionID	registerInvalidSolutionInternal (const std::string &message, const bool warning) const

Static Protected Member Functions
template<typename T , typename Functor >
static void	xyDerivatives (const T x, const T &y, T &z, T &dz_dx, T &dz_dy, const Functor &z_from_x_y)
	Computes the dependent variable z and its derivatives with respect to the independent variables x and y using the simple two parameter `z_from_x_y` functor. More...

template<typename T >
static std::pair< T, T >	makeZeroAndOne (const T &)
	Given a type example, this method returns zero and unity representations of that type (first and second members of returned pair respectively) More...

static std::string	meshPropertyName (const std::string &data_name, const std::string &prefix)

Protected Attributes
const Real	_Mh2
	Hydrogen molar mass (kg/mol) More...

const Real	_p_critical
	Critical pressure (Pa) More...

const Real	_T_critical
	Critical temperature (K) More...

const Real	_rho_molar_critical
	Critical molar density (mol/l) More...

const Real	_rho_critical
	Critical density (kg/m^3) More...

const Real	_p_triple
	Triple point pressure (Pa) More...

const Real	_T_triple
	Triple point temperature (K) More...

const std::array< Real, 5 >	_a {{1.616, -0.4117, -0.792, 0.758, 1.217}}
	Coefficients for ideal gas component of the Helmholtz free energy. More...

const std::array< Real, 5 >	_b

const std::array< Real, 7 >	_N1 {{-6.93643, 0.01, 2.1101, 4.52059, 0.732564, -1.34086, 0.130985}}
	Coefficients for residual component of the Helmholtz free energy. More...

const std::array< Real, 7 >	_t1 {{0.6844, 1.0, 0.989, 0.489, 0.803, 1.1444, 1.409}}

const std::array< unsigned int, 7 >	_d1 {{1, 4, 1, 1, 2, 2, 3}}

const std::array< Real, 2 >	_N2 {{-0.777414, 0.351944}}

const std::array< Real, 2 >	_t2 {{1.754, 1.311}}

const std::array< unsigned int, 2 >	_d2 {{1, 3}}

const std::array< Real, 5 >	_N3 {{-0.0211716, 0.0226312, 0.032187, -0.0231752, 0.0557346}}

const std::array< Real, 5 >	_t3 {{4.187, 5.646, 0.791, 7.249, 2.986}}

const std::array< unsigned int, 5 >	_d3 {{2, 1, 3, 1, 1}}

const std::array< Real, 5 >	_phi3 {{-1.685, -0.489, -0.103, -2.506, -1.607}}

const std::array< Real, 5 >	_beta3 {{-0.171, -0.2245, -0.1304, -0.2785, -0.3967}}

const std::array< Real, 5 >	_gamma3 {{0.7164, 1.3444, 1.4517, 0.7204, 1.5445}}

const std::array< Real, 5 >	_D3 {{1.506, 0.156, 1.736, 0.67, 1.662}}

const std::array< Real, 5 >	_amu {{2.09630e-1, -4.55274e-1, 1.423602e-1, -3.35325e-2, 2.76981e-3}}
	Coefficients for viscosity. More...

const std::array< Real, 7 >	_bmu {{-0.187, 2.4871, 3.7151, -11.0972, 9.0965, -3.8292, 0.5166}}

const std::array< Real, 6 >	_cmu

const std::array< Real, 7 >	_a1k
	Coefficients for thermal conductivity. More...

const std::array< Real, 4 >	_a2k {{1.38497e2, -2.21878e1, 4.57151, 1.0}}

const std::array< Real, 5 >	_b1k {{3.63081e-2, -2.07629e-2, 3.1481e-2, -1.43097e-2, 1.7498e-3}}

const std::array< Real, 5 >	_b2k {{1.8337e-3, -8.86716e-3, 1.5826e-2, -1.06283e-2, 2.80673e-3}}

const Real	_tolerance
	Newton's method may be used to convert between variable sets. More...

const Real	_T_initial_guess
	Initial guess for temperature (or temperature used to compute the initial guess) More...

const Real	_p_initial_guess
	Initial guess for pressure (or pressure used to compute the initial guess) More...

const unsigned int	_max_newton_its
	Maximum number of iterations for the variable conversion newton solves. More...

const Real	_T_c2k
	Conversion of temperature from Celsius to Kelvin. More...

const bool	_allow_imperfect_jacobians
	Flag to set unimplemented Jacobian entries to zero. More...

SubProblem &	_subproblem

FEProblemBase &	_fe_problem

SystemBase &	_sys

const THREAD_ID	_tid

Assembly &	_assembly

const Moose::CoordinateSystemType &	_coord_sys

const bool	_duplicate_initial_execution

std::set< std::string >	_depend_uo

const bool &	_enabled

MooseApp &	_app

const std::string	_type

const std::string	_name

const InputParameters &	_pars

Factory &	_factory

ActionFactory &	_action_factory

const ExecFlagEnum &	_execute_enum

const ExecFlagType &	_current_execute_flag

MooseApp &	_restartable_app

const std::string	_restartable_system_name

const THREAD_ID	_restartable_tid

const bool	_restartable_read_only

FEProblemBase &	_mci_feproblem

FEProblemBase &	_mdi_feproblem

MooseApp &	_pg_moose_app

const std::string	_prefix

FEProblemBase &	_sc_fe_problem

const THREAD_ID	_sc_tid

const Real &	_real_zero

const VariableValue &	_scalar_zero

const Point &	_point_zero

const InputParameters &	_mi_params

const std::string	_mi_name

const MooseObjectName	_mi_moose_object_name

FEProblemBase &	_mi_feproblem

SubProblem &	_mi_subproblem

const THREAD_ID	_mi_tid

const Moose::MaterialDataType	_material_data_type

MaterialData &	_material_data

bool	_stateful_allowed

bool	_get_material_property_called

std::vector< std::unique_ptr< PropertyValue > >	_default_properties

std::unordered_set< unsigned int >	_material_property_dependencies

const MaterialPropertyName	_get_suffix

const bool	_use_interpolated_state

const InputParameters &	_ti_params

FEProblemBase &	_ti_feproblem

bool	_is_implicit

Real &	_t

const Real &	_t_old

int &	_t_step

Real &	_dt

Real &	_dt_old

bool	_is_transient

const Parallel::Communicator &	_communicator

Static Protected Attributes
static const std::string	_interpolated_old

static const std::string	_interpolated_older

	propfunc (p, v, e) propfunc(T
	Compute a fluid property given for the state defined by two given properties. More...

e	propfunc (c, v, e) propfunc(cp

e e	propfunc (cv, v, e) propfunc(mu

e e e	propfunc (k, v, e) propfuncWithDefault(s

e e e e	propfunc (s, h, p) propfunc(rho

e e e e s	propfunc (e, v, h) propfuncWithDefault(s

e e e e s T	propfunc (pp_sat, p, T) propfunc(mu

e e e e s T T	propfunc (k, rho, T) propfuncWithDefault(c

e e e e s T T T T T	propfunc (rho, p, T) propfunc(e

e e e e s T T T T T rho	propfunc (e, T, v) propfunc(p

e e e e s T T T T T rho v	propfunc (h, T, v) propfunc(s

e e e e s T T T T T rho v v	propfunc (cv, T, v) propfunc(h

e e e e s T T T	propfuncWithDefault (cp, p, T) propfuncWithDefault(cv

e e e e s T T T T	propfuncWithDefault (mu, p, T) propfuncWithDefault(k

e e e e s T T T T T rho v v T	propfuncWithDefault (h, v, e) propfunc(g

e e e e s T T T T T rho v v T e	propfuncWithDefault (p, h, s) propfunc(T

e e e e s T T T T T rho v v T e p	propfuncWithDefault (T, p, h) propfuncWithDefault(beta

e e e e s T T T T T rho v v T e p T	propfuncWithDefault (v, p, T) propfuncWithDefault(e

e e e e s T T T T T rho v v T e p T T	propfuncWithDefault (gamma, v, e) propfuncWithDefault(gamma

	v

e	v

e e	v

e e e	v

e e e e s T T T T T rho v v T	v

e e e e	p

e e e e s	p

e e e e s T T	p

e e e e s T T T	p

e e e e s T T T T	p

e e e e s T T T T T	p

e e e e s T T T T T rho v v	p

e e e e s T T T T T rho v v T e p	p

e e e e s T T T T T rho v v T e p T	p

e e e e s T T T T T rho v v T e p T T	p

e e e e s T	rho

e e e e s T T T T T rho	T

e e e e s T T T T T rho v	T

e e e e s T T T T T rho v v T e	h

Detailed Description

Hydrogen (H2) fluid properties as a function of pressure (Pa) and temperature (K).

Thermodynamic properties calculated from: Leachman, Jacobsen, Penoncello,and Lemmon, Fundamental equations of state for parahydrogen, normal hydrogen, and orthohydrogen, Journal of Physical and Chemical Reference Data, 38, 721–748 (2009)

Viscosity from: Muzny, Huber and Kazakov, Correlation for the viscosity of normal hydrogen obtained from symbolic regression, Journal of Chemical and Engineering Data, 58, 969-979 (2013)

Thermal conductivity from: Assael, Assael, Huber, Perkins and Takata, Correlation of the thermal conductivity of normal and parahydrogen from the triple point to 1000 K and up to 100 Mpa, Journal of Physical and Chemical Reference Data, 40 (2011)

Definition at line 37 of file HydrogenFluidProperties.h.

Constructor & Destructor Documentation

◆ HydrogenFluidProperties()

HydrogenFluidProperties::HydrogenFluidProperties ( const InputParameters & parameters )

Definition at line 25 of file HydrogenFluidProperties.C.

   : HelmholtzFluidProperties(parameters),
     _Mh2(2.01588e-3),
     _p_critical(1.315e6),
     _T_critical(33.19),
     _rho_molar_critical(15.508),
     _rho_critical(1000.0 * _rho_molar_critical * _Mh2),
     _p_triple(7.7e3),
     _T_triple(13.952)
 {
 }

Member Function Documentation

◆ alpha()

Real HydrogenFluidProperties::alpha	(	Real	delta,
		Real	tau
	)		const

overrideprotectedvirtual

Helmholtz free energy for H2 From Leachman et al (reference above)

Parameters

delta	scaled density (-)
tau	scaled temperature (-)

Returns: alpha Helmholtz free energy

Implements HelmholtzFluidProperties.

Definition at line 299 of file HydrogenFluidProperties.C.

 {
   // Ideal gas component of the Helmholtz free energy
   Real alpha0 = std::log(delta) + 1.5 * std::log(tau) - 1.4579856475 + 1.888076782 * tau;
 
   for (std::size_t i = 0; i < _a.size(); ++i)
     alpha0 += _a[i] * std::log(1.0 - std::exp(_b[i] * tau));
 
   // Residual component of the Helmholtz free energy
   Real alphar = 0.0;
 
   for (std::size_t i = 0; i < _t1.size(); ++i)
     alphar += _N1[i] * MathUtils::pow(delta, _d1[i]) * std::pow(tau, _t1[i]);
 
   for (std::size_t i = 0; i < _t2.size(); ++i)
     alphar += _N2[i] * MathUtils::pow(delta, _d2[i]) * std::pow(tau, _t2[i]) * std::exp(-delta);
 
   for (std::size_t i = 0; i < _t3.size(); ++i)
     alphar += _N3[i] * MathUtils::pow(delta, _d3[i]) * std::pow(tau, _t3[i]) *
               std::exp(_phi3[i] * Utility::pow<2>(delta - _D3[i]) +
                        _beta3[i] * Utility::pow<2>(tau - _gamma3[i]));
 
   // The Helmholtz free energy is the sum of these two
   return alpha0 + alphar;
 }

◆ c_from_p_T()

Real HelmholtzFluidProperties::c_from_p_T	(	Real	pressure,
		Real	temperature
	)		const

overridevirtualinherited

Definition at line 98 of file HelmholtzFluidProperties.C.

 {
   // Require density first
   const Real density = rho_from_p_T(pressure, temperature);
   // Scale the input density and temperature
   const Real delta = density / criticalDensity();
   const Real tau = criticalTemperature() / temperature;
 
   const Real da_dd = dalpha_ddelta(delta, tau);
 
   Real w = 2.0 * delta * da_dd + delta * delta * d2alpha_ddelta2(delta, tau);
   w -= Utility::pow<2>(delta * da_dd - delta * tau * d2alpha_ddeltatau(delta, tau)) /
        (tau * tau * d2alpha_dtau2(delta, tau));
 
   return std::sqrt(_R * temperature * w / molarMass());
 }

◆ cp_from_p_T()

Real HelmholtzFluidProperties::cp_from_p_T	(	Real	pressure,
		Real	temperature
	)		const

overridevirtualinherited

Definition at line 116 of file HelmholtzFluidProperties.C.

 {
   // Require density first
   const Real density = rho_from_p_T(pressure, temperature);
   // Scale the input density and temperature
   const Real delta = density / criticalDensity();
   const Real tau = criticalTemperature() / temperature;
 
   const Real da_dd = dalpha_ddelta(delta, tau);
 
   const Real cp = _R *
                   (-tau * tau * d2alpha_dtau2(delta, tau) +
                    Utility::pow<2>(delta * da_dd - delta * tau * d2alpha_ddeltatau(delta, tau)) /
                        (2.0 * delta * da_dd + delta * delta * d2alpha_ddelta2(delta, tau))) /
                   molarMass();
 
   return cp;
 }

◆ criticalDensity()

Real HydrogenFluidProperties::criticalDensity ( ) const

overridevirtual

Critical density.

Returns: critical density (kg/m^3)

Reimplemented from SinglePhaseFluidProperties.

Definition at line 62 of file HydrogenFluidProperties.C.

 {
   return _rho_critical;
 }

◆ criticalInternalEnergy()

Real SinglePhaseFluidProperties::criticalInternalEnergy ( ) const

virtualinherited

Critical specific internal energy.

Returns: specific internal energy (J/kg)

Reimplemented in IdealGasFluidProperties, CaloricallyImperfectGas, and StiffenedGasFluidProperties.

Definition at line 316 of file SinglePhaseFluidProperties.C.

 {
   return e_from_p_rho(criticalPressure(), criticalDensity());
 }

◆ criticalPressure()

Real HydrogenFluidProperties::criticalPressure ( ) const

overridevirtual

Critical pressure.

Returns: critical pressure (Pa)

Reimplemented from SinglePhaseFluidProperties.

Definition at line 50 of file HydrogenFluidProperties.C.

 {
   return _p_critical;
 }

◆ criticalTemperature()

Real HydrogenFluidProperties::criticalTemperature ( ) const

overridevirtual

Critical temperature.

Returns: critical temperature (K)

Reimplemented from SinglePhaseFluidProperties.

Definition at line 56 of file HydrogenFluidProperties.C.

 {
   return _T_critical;
 }

◆ cv_from_p_T()

Real HelmholtzFluidProperties::cv_from_p_T	(	Real	pressure,
		Real	temperature
	)		const

overridevirtualinherited

Definition at line 136 of file HelmholtzFluidProperties.C.

 {
   // Require density first
   const Real density = rho_from_p_T(pressure, temperature);
   // Scale the input density and temperature
   const Real delta = density / criticalDensity();
   const Real tau = criticalTemperature() / temperature;
 
   return -_R * tau * tau * d2alpha_dtau2(delta, tau) / molarMass();
 }

◆ d2alpha_ddelta2()

Real HydrogenFluidProperties::d2alpha_ddelta2	(	Real	delta,
		Real	tau
	)		const

overrideprotectedvirtual

Second derivative of Helmholtz free energy wrt delta.

Parameters

delta	scaled density (-)
tau	scaled temperature (-)

Returns: second derivative of Helmholtz free energy wrt delta

Implements HelmholtzFluidProperties.

Definition at line 381 of file HydrogenFluidProperties.C.

 {
   // Ideal gas component of the Helmholtz free energy
   Real dalpha0 = -1.0 / delta / delta;
 
   // Residual component of the Helmholtz free energy
   Real dalphar = 0.0;
 
   for (std::size_t i = 0; i < _t1.size(); ++i)
     dalphar +=
         _N1[i] * _d1[i] * (_d1[i] - 1.0) * MathUtils::pow(delta, _d1[i]) * std::pow(tau, _t1[i]);
 
   for (std::size_t i = 0; i < _t2.size(); ++i)
     dalphar += _N2[i] * MathUtils::pow(delta, _d2[i]) * std::pow(tau, _t2[i]) * std::exp(-delta) *
                (delta * delta - 2.0 * _d2[i] * delta + _d2[i] * (_d2[i] - 1.0));
 
   for (std::size_t i = 0; i < _t3.size(); ++i)
     dalphar += _N3[i] * MathUtils::pow(delta, _d3[i]) * std::pow(tau, _t3[i]) *
                std::exp(_phi3[i] * Utility::pow<2>(delta - _D3[i]) +
                         _beta3[i] * Utility::pow<2>(tau - _gamma3[i])) *
                (_d3[i] * _d3[i] +
                 2.0 * delta * delta * _phi3[i] *
                     (1.0 + 2.0 * _phi3[i] * (_D3[i] - delta) * (_D3[i] - delta)) +
                 _d3[i] * (4.0 * delta * _phi3[i] * (delta - _D3[i]) - 1.0));
 
   // The Helmholtz free energy is the sum of these two
   return dalpha0 + dalphar / delta / delta;
 }

◆ d2alpha_ddeltatau()

Real HydrogenFluidProperties::d2alpha_ddeltatau	(	Real	delta,
		Real	tau
	)		const

overrideprotectedvirtual

Second derivative of Helmholtz free energy wrt delta and tau.

Parameters

delta	scaled density (-)
tau	scaled temperature (-)

Returns: second derivative of Helmholtz free energy wrt delta and tau

Implements HelmholtzFluidProperties.

Definition at line 447 of file HydrogenFluidProperties.C.

 {
   // Residual component of the Helmholtz free energy
   Real dalphar = 0.0;
 
   for (std::size_t i = 0; i < _t1.size(); ++i)
     dalphar += _N1[i] * _d1[i] * _t1[i] * std::pow(delta, _d1[i]) * std::pow(tau, _t1[i]);
 
   for (std::size_t i = 0; i < _t2.size(); ++i)
     dalphar += _N2[i] * _t2[i] * std::pow(delta, _d2[i]) * std::pow(tau, _t2[i]) *
                std::exp(-delta) * (_d2[i] - delta);
 
   for (std::size_t i = 0; i < _t3.size(); ++i)
     dalphar += _N3[i] * std::pow(delta, _d3[i]) * std::pow(tau, _t3[i]) *
                std::exp(_phi3[i] * Utility::pow<2>(delta - _D3[i]) +
                         _beta3[i] * Utility::pow<2>(tau - _gamma3[i])) *
                (_d3[i] + delta * (2.0 * _phi3[i] * (delta - _D3[i]))) *
                (_t3[i] + 2.0 * _beta3[i] * tau * (tau - _gamma3[i]));
 
   // The Helmholtz free energy is the sum of these two
   return dalphar / delta / tau;
 }

◆ d2alpha_dtau2()

Real HydrogenFluidProperties::d2alpha_dtau2	(	Real	delta,
		Real	tau
	)		const

overrideprotectedvirtual

Second derivative of Helmholtz free energy wrt tau.

Parameters

delta	scaled density (-)
tau	scaled temperature (-)

Returns: second derivative of Helmholtz free energy wrt tau

Implements HelmholtzFluidProperties.

Definition at line 411 of file HydrogenFluidProperties.C.

 {
   // Ideal gas component of the Helmholtz free energy
   Real dalpha0 = -1.5 / tau / tau;
 
   for (std::size_t i = 0; i < _a.size(); ++i)
   {
     Real exptau = std::exp(_b[i] * tau);
     dalpha0 -= _a[i] * (_b[i] * _b[i] * exptau / (1.0 - exptau) * (exptau / (1.0 - exptau) + 1.0));
   }
 
   // Residual component of the Helmholtz free energy
   Real dalphar = 0.0;
 
   for (std::size_t i = 0; i < _t1.size(); ++i)
     dalphar +=
         _N1[i] * _t1[i] * (_t1[i] - 1.0) * MathUtils::pow(delta, _d1[i]) * std::pow(tau, _t1[i]);
 
   for (std::size_t i = 0; i < _t2.size(); ++i)
     dalphar += _N2[i] * _t2[i] * (_t2[i] - 1.0) * MathUtils::pow(delta, _d2[i]) *
                std::pow(tau, _t2[i]) * std::exp(-delta);
 
   for (std::size_t i = 0; i < _t3.size(); ++i)
     dalphar += _N3[i] * MathUtils::pow(delta, _d3[i]) * std::pow(tau, _t3[i]) *
                std::exp(_phi3[i] * Utility::pow<2>(delta - _D3[i]) +
                         _beta3[i] * Utility::pow<2>(tau - _gamma3[i])) *
                (_t3[i] * _t3[i] +
                 2.0 * _beta3[i] * tau * tau *
                     (1.0 + 2.0 * _beta3[i] * MathUtils::pow(tau - _gamma3[i], 2)) -
                 _t3[i] * (1.0 + 4.0 * _beta3[i] * tau * (tau - _gamma3[i])));
 
   // The Helmholtz free energy is the sum of these two
   return dalpha0 + dalphar / tau / tau;
 }

◆ dalpha_ddelta()

Real HydrogenFluidProperties::dalpha_ddelta	(	Real	delta,
		Real	tau
	)		const

overrideprotectedvirtual

Derivative of Helmholtz free energy wrt delta.

Parameters

delta	scaled density (-)
tau	scaled temperature (-)

Returns: derivative of Helmholtz free energy wrt delta

Implements HelmholtzFluidProperties.

Definition at line 326 of file HydrogenFluidProperties.C.

 {
   // Ideal gas component of the Helmholtz free energy
   Real dalpha0 = 1.0 / delta;
 
   // Residual component of the Helmholtz free energy
   Real dalphar = 0.0;
 
   for (std::size_t i = 0; i < _t1.size(); ++i)
     dalphar += _N1[i] * _d1[i] * MathUtils::pow(delta, _d1[i]) * std::pow(tau, _t1[i]);
 
   for (std::size_t i = 0; i < _t2.size(); ++i)
     dalphar += _N2[i] * MathUtils::pow(delta, _d2[i]) * std::pow(tau, _t2[i]) * std::exp(-delta) *
                (_d2[i] - delta);
 
   for (std::size_t i = 0; i < _t3.size(); ++i)
     dalphar += _N3[i] * MathUtils::pow(delta, _d3[i]) * std::pow(tau, _t3[i]) *
                std::exp(_phi3[i] * Utility::pow<2>(delta - _D3[i]) +
                         _beta3[i] * Utility::pow<2>(tau - _gamma3[i])) *
                (_d3[i] + delta * (2.0 * _phi3[i] * (delta - _D3[i])));
 
   // The Helmholtz free energy is the sum of these two
   return dalpha0 + dalphar / delta;
 }

◆ dalpha_dtau()

Real HydrogenFluidProperties::dalpha_dtau	(	Real	delta,
		Real	tau
	)		const

overrideprotectedvirtual

Derivative of Helmholtz free energy wrt tau.

Parameters

delta	scaled density (-)
tau	scaled temperature (-)

Returns: derivative of Helmholtz free energy wrt tau

Implements HelmholtzFluidProperties.

Definition at line 352 of file HydrogenFluidProperties.C.

 {
   // Ideal gas component of the Helmholtz free energy
   Real dalpha0 = 1.5 / tau + 1.888076782;
 
   for (std::size_t i = 0; i < _a.size(); ++i)
     dalpha0 += _a[i] * _b[i] * (1.0 - 1.0 / (1.0 - std::exp(_b[i] * tau)));
 
   // Residual component of the Helmholtz free energy
   Real dalphar = 0.0;
 
   for (std::size_t i = 0; i < _t1.size(); ++i)
     dalphar += _N1[i] * _t1[i] * MathUtils::pow(delta, _d1[i]) * std::pow(tau, _t1[i]);
 
   for (std::size_t i = 0; i < _t2.size(); ++i)
     dalphar +=
         _N2[i] * _t2[i] * MathUtils::pow(delta, _d2[i]) * std::pow(tau, _t2[i]) * std::exp(-delta);
 
   for (std::size_t i = 0; i < _t3.size(); ++i)
     dalphar += _N3[i] * MathUtils::pow(delta, _d3[i]) * std::pow(tau, _t3[i]) *
                std::exp(_phi3[i] * Utility::pow<2>(delta - _D3[i]) +
                         _beta3[i] * Utility::pow<2>(tau - _gamma3[i])) *
                (_t3[i] + tau * (2.0 * _beta3[i] * (tau - _gamma3[i])));
 
   // The Helmholtz free energy is the sum of these two
   return dalpha0 + dalphar / tau;
 }

◆ e_from_p_T() [1/2]

Real HelmholtzFluidProperties::e_from_p_T	(	Real	pressure,
		Real	temperature
	)		const

overridevirtualinherited

Definition at line 63 of file HelmholtzFluidProperties.C.

Referenced by HelmholtzFluidProperties::e_from_p_T().

 {
   // Require density first
   const Real density = rho_from_p_T(pressure, temperature);
   // Scale the input density and temperature
   const Real delta = density / criticalDensity();
   const Real tau = criticalTemperature() / temperature;
 
   return _R * temperature * tau * dalpha_dtau(delta, tau) / molarMass();
 }

◆ e_from_p_T() [2/2]

void HelmholtzFluidProperties::e_from_p_T	(	Real	p,
		Real	T,
		Real &	e,
		Real &	de_dp,
		Real &	de_dT
	)		const

overridevirtualinherited

Definition at line 75 of file HelmholtzFluidProperties.C.

 {
   e = this->e_from_p_T(pressure, temperature);
 
   // Require density first
   const Real density = rho_from_p_T(pressure, temperature);
   // Scale the input density and temperature
   const Real delta = density / criticalDensity();
   const Real tau = criticalTemperature() / temperature;
 
   const Real da_dd = dalpha_ddelta(delta, tau);
   const Real d2a_dd2 = d2alpha_ddelta2(delta, tau);
   const Real d2a_ddt = d2alpha_ddeltatau(delta, tau);
 
   de_dp = tau * d2a_ddt / (density * (2.0 * da_dd + delta * d2a_dd2));
   de_dT = -_R *
           (delta * tau * d2a_ddt * (da_dd - tau * d2a_ddt) / (2.0 * da_dd + delta * d2a_dd2) +
            tau * tau * d2alpha_dtau2(delta, tau)) /
           molarMass();
 }

◆ e_spndl_from_v()

Real SinglePhaseFluidProperties::e_spndl_from_v ( Real v ) const

virtualinherited

Specific internal energy from temperature and specific volume.

Parameters

[in]	T	temperature
[in]	v	specific volume

Reimplemented in IdealGasFluidProperties, CaloricallyImperfectGas, and StiffenedGasFluidProperties.

Definition at line 479 of file SinglePhaseFluidProperties.C.

 {
   mooseError(__PRETTY_FUNCTION__, " not implemented.");
 }

◆ execute()

virtual void FluidProperties::execute ( )

inlinefinalvirtualinherited

Implements ThreadedGeneralUserObject.

Definition at line 33 of file FluidProperties.h.

33 {}

◆ finalize()

virtual void FluidProperties::finalize ( )

inlinefinalvirtualinherited

Implements ThreadedGeneralUserObject.

Definition at line 35 of file FluidProperties.h.

35 {}

◆ fluidName()

std::string HydrogenFluidProperties::fluidName ( ) const

overridevirtual

Fluid name.

Returns: string representing fluid name

Reimplemented from SinglePhaseFluidProperties.

Definition at line 38 of file HydrogenFluidProperties.C.

 {
   return "hydrogen";
 }

◆ h_from_p_T() [1/2]

Real HelmholtzFluidProperties::h_from_p_T	(	Real	pressure,
		Real	temperature
	)		const

overridevirtualinherited

Definition at line 183 of file HelmholtzFluidProperties.C.

Referenced by HelmholtzFluidProperties::h_from_p_T(), and HelmholtzFluidProperties::T_from_p_h().

 {
   // Require density first
   const Real density = rho_from_p_T(pressure, temperature);
   // Scale the input density and temperature
   const Real delta = density / criticalDensity();
   const Real tau = criticalTemperature() / temperature;
 
   return _R * temperature * (tau * dalpha_dtau(delta, tau) + delta * dalpha_ddelta(delta, tau)) /
          molarMass();
 }

◆ h_from_p_T() [2/2]

void HelmholtzFluidProperties::h_from_p_T	(	Real	p,
		Real	T,
		Real &	h,
		Real &	dh_dp,
		Real &	dh_dT
	)		const

overridevirtualinherited

Definition at line 196 of file HelmholtzFluidProperties.C.

 {
   h = this->h_from_p_T(pressure, temperature);
 
   // Require density first
   const Real density = rho_from_p_T(pressure, temperature);
   // Scale the input density and temperature
   const Real delta = density / criticalDensity();
   const Real tau = criticalTemperature() / temperature;
 
   const Real da_dd = dalpha_ddelta(delta, tau);
   const Real d2a_dd2 = d2alpha_ddelta2(delta, tau);
   const Real d2a_ddt = d2alpha_ddeltatau(delta, tau);
 
   dh_dp = (da_dd + delta * d2a_dd2 + tau * d2a_ddt) / (density * (2.0 * da_dd + delta * d2a_dd2));
   dh_dT = _R *
           (delta * da_dd * (1.0 - tau * d2a_ddt / da_dd) * (1.0 - tau * d2a_ddt / da_dd) /
                (2.0 + delta * d2a_dd2 / da_dd) -
            tau * tau * d2alpha_dtau2(delta, tau)) /
           molarMass();
 }

◆ henryCoefficients()

std::vector< Real > HydrogenFluidProperties::henryCoefficients ( ) const

overridevirtual

Henry's law coefficients for dissolution in water.

Returns: Henry's constant coefficients

Reimplemented from SinglePhaseFluidProperties.

Definition at line 271 of file HydrogenFluidProperties.C.

 {
   return {-4.73284, 6.08954, 6.06066};
 }

◆ initialize()

virtual void FluidProperties::initialize ( )

inlinefinalvirtualinherited

Implements ThreadedGeneralUserObject.

Definition at line 34 of file FluidProperties.h.

34 {}

◆ k_from_p_T() [1/2]

Real HydrogenFluidProperties::k_from_p_T	(	Real	pressure,
		Real	temperature
	)		const

overridevirtual

Definition at line 249 of file HydrogenFluidProperties.C.

Referenced by k_from_p_T().

 {
   // Require density first
   const Real density = rho_from_p_T(pressure, temperature);
   return k_from_rho_T(density, temperature);
 }

◆ k_from_p_T() [2/2]

void HydrogenFluidProperties::k_from_p_T	(	Real	pressure,
		Real	temperature,
		Real &	k,
		Real &	dk_dp,
		Real &	dk_dT
	)		const

overridevirtual

Definition at line 257 of file HydrogenFluidProperties.C.

 {
   k = this->k_from_p_T(pressure, temperature);
   // Calculate derivatives using finite differences
   const Real eps = 1.0e-6;
   const Real peps = pressure * eps;
   const Real Teps = temperature * eps;
 
   dk_dp = (this->k_from_p_T(pressure + peps, temperature) - k) / peps;
   dk_dT = (this->k_from_p_T(pressure, temperature + Teps) - k) / Teps;
 }

◆ k_from_rho_T()

Real HydrogenFluidProperties::k_from_rho_T	(	Real	density,
		Real	temperature
	)		const

overridevirtual

Definition at line 218 of file HydrogenFluidProperties.C.

Referenced by k_from_p_T().

 {
   // // Scaled variables
   const Real Tr = temperature / 33.145;
   const Real rhor = density / 31.262;
 
   // The ideal gas component
   Real sum1 = 0.0;
   for (std::size_t i = 0; i < _a1k.size(); ++i)
     sum1 += _a1k[i] * MathUtils::pow(Tr, i);
 
   Real sum2 = 0.0;
   for (std::size_t i = 0; i < _a2k.size(); ++i)
     sum2 += _a2k[i] * MathUtils::pow(Tr, i);
 
   const Real lambda0 = sum1 / sum2;
 
   // The excess contribution due to density
   Real lambdah = 0.0;
   for (std::size_t i = 0; i < _b1k.size(); ++i)
     lambdah += (_b1k[i] + _b2k[i] * Tr) * MathUtils::pow(rhor, i + 1);
 
   // The critical enhancement
   const Real lambdac = 6.24e-4 / (-2.58e-7 + std::abs(Tr - 1.0)) *
                        std::exp(-MathUtils::pow(0.837 * (rhor - 1.0), 2));
 
   // The thermal conductivity
   return lambda0 + lambdah + lambdac;
 }

◆ makeZeroAndOne() [1/2]

template<typename T >

std::pair< T, T > SinglePhaseFluidProperties::makeZeroAndOne ( const T & )

staticprotectedinherited

Given a type example, this method returns zero and unity representations of that type (first and second members of returned pair respectively)

Definition at line 476 of file SinglePhaseFluidProperties.h.

Referenced by SinglePhaseFluidProperties::xyDerivatives().

 {
   return {T{0, 0}, T{1, 0}};
 }

◆ makeZeroAndOne() [2/2]

template<>

std::pair<Real, Real> SinglePhaseFluidProperties::makeZeroAndOne ( const Real & )

inlineinherited

Definition at line 483 of file SinglePhaseFluidProperties.h.

 {
   return {Real{0}, Real{1}};
 }

◆ molarMass()

Real HydrogenFluidProperties::molarMass ( ) const

overridevirtual

Molar mass [kg/mol].

Returns: molar mass

Reimplemented from SinglePhaseFluidProperties.

Definition at line 44 of file HydrogenFluidProperties.C.

 {
   return _Mh2;
 }

◆ mu_from_p_T() [1/2]

Real HydrogenFluidProperties::mu_from_p_T	(	Real	pressure,
		Real	temperature
	)		const

overridevirtual

Definition at line 172 of file HydrogenFluidProperties.C.

 {
   // Require density first
   const Real density = rho_from_p_T(pressure, temperature);
   return mu_from_rho_T(density, temperature);
 }

◆ mu_from_p_T() [2/2]

void HydrogenFluidProperties::mu_from_p_T	(	Real	pressure,
		Real	temperature,
		Real &	mu,
		Real &	dmu_dp,
		Real &	dmu_dT
	)		const

overridevirtual

Definition at line 180 of file HydrogenFluidProperties.C.

 {
   Real rho, drho_dp, drho_dT;
   rho_from_p_T(pressure, temperature, rho, drho_dp, drho_dT);
 
   Real dmu_drho;
   mu_from_rho_T(rho, temperature, drho_dT, mu, dmu_drho, dmu_dT);
   dmu_dp = dmu_drho * drho_dp;
 }

◆ mu_from_rho_T() [1/2]

Real HydrogenFluidProperties::mu_from_rho_T	(	Real	density,
		Real	temperature
	)		const

overridevirtual

Definition at line 80 of file HydrogenFluidProperties.C.

Referenced by mu_from_p_T(), and rho_mu_from_p_T().

 {
   // Scaled variables
   const Real Tstar = temperature / 30.41;
   const Real logTstar = std::log(Tstar);
   const Real Tr = temperature / _T_critical;
   const Real rhor = density / 90.5;
 
   // Ideal gas component
   Real sum = 0.0;
   for (std::size_t i = 0; i < _amu.size(); ++i)
     sum += _amu[i] * MathUtils::pow(logTstar, i);
 
   const Real mu0 =
       0.021357 * std::sqrt(1000.0 * _Mh2 * temperature) / (0.297 * 0.297 * std::exp(sum));
 
   // The excess contribution due to density
   Real sumr = 0.0;
   for (std::size_t i = 0; i < _bmu.size(); ++i)
     sumr += _bmu[i];
 
   const Real mu1 = MathUtils::pow(0.297, 3) * sumr * mu0 / Tstar;
 
   // The viscosity is then
   const Real mu =
       mu0 + mu1 * density +
       _cmu[0] * rhor * rhor *
           std::exp(_cmu[1] * Tr + _cmu[2] / Tr + _cmu[3] * rhor * rhor / (_cmu[4] + Tr) +
                    _cmu[5] * MathUtils::pow(rhor, 6));
 
   return mu * 1.0e-6;
 }

◆ mu_from_rho_T() [2/2]

void HydrogenFluidProperties::mu_from_rho_T	(	Real	density,
		Real	temperature,
		Real	ddensity_dT,
		Real &	mu,
		Real &	dmu_drho,
		Real &	dmu_dT
	)		const

Definition at line 114 of file HydrogenFluidProperties.C.

 {
   // Scaled variables
   const Real Tstar = temperature / 30.41;
   const Real logTstar = std::log(Tstar);
   const Real Tr = temperature / _T_critical;
   const Real rhor = density / 90.5;
   const Real drhor_drho = 1.0 / 90.5;
   const Real dTr_dT = 1.0 / _T_critical;
 
   // The dilute gas component
   Real sum = 0.0, dsum_dT = 0.0;
   for (std::size_t i = 0; i < _amu.size(); ++i)
   {
     sum += _amu[i] * MathUtils::pow(logTstar, i);
     dsum_dT += i * _amu[i] * MathUtils::pow(logTstar, i) / (temperature * logTstar);
   }
 
   const Real mu0 =
       0.021357 * std::sqrt(1000.0 * _Mh2 * temperature) / (0.297 * 0.297 * std::exp(sum));
   const Real dmu0_dT = 21.357 * _Mh2 * (1.0 - 2.0 * temperature * dsum_dT) * std::exp(-sum) /
                        (2.0 * std::sqrt(1000.0 * _Mh2 * temperature) * 0.297 * 0.297);
 
   // The excess contribution due to density
   Real sumr = 0.0;
   for (std::size_t i = 0; i < _bmu.size(); ++i)
     sumr += _bmu[i];
 
   const Real mu1 = MathUtils::pow(0.297, 3) * sumr * mu0 / Tstar;
   const Real dmu1_dT =
       MathUtils::pow(0.297, 3) * sumr * (dmu0_dT / Tstar - mu0 / (30.41 * Tstar * Tstar));
 
   // The viscosity and derivatives are then
   const Real exponent = _cmu[1] * Tr + _cmu[2] / Tr + _cmu[3] * rhor * rhor / (_cmu[4] + Tr) +
                         _cmu[5] * MathUtils::pow(rhor, 6);
   const Real dexponent_drho =
       (2.0 * _cmu[3] * rhor / (_cmu[4] + Tr) + 6.0 * _cmu[5] * MathUtils::pow(rhor, 5)) *
       drhor_drho;
   const Real dexponent_dT =
       (_cmu[1] - _cmu[2] / Tr / Tr - _cmu[3] * rhor * rhor / (_cmu[4] + Tr) / (_cmu[4] + Tr)) *
       dTr_dT;
 
   mu = (mu0 + mu1 * density + _cmu[0] * rhor * rhor * std::exp(exponent)) * 1.0e-6;
   dmu_drho =
       (mu1 + _cmu[0] * rhor * std::exp(exponent) * (2.0 * drhor_drho + rhor * dexponent_drho)) *
       1.0e-6;
   dmu_dT =
       (dmu0_dT + density * dmu1_dT + _cmu[0] * rhor * rhor * dexponent_dT * std::exp(exponent)) *
           1.0e-6 +
       dmu_drho * ddensity_dT;
 }

◆ p_from_rho_T()

Real HelmholtzFluidProperties::p_from_rho_T	(	Real	rho,
		Real	T
	)		const

virtualinherited

Pressure as a function of density and temperature.

Parameters

rho	density (kg/m^3)
T	temperature (K)

Returns: pressure (Pa)

Reimplemented in CO2FluidProperties.

Definition at line 238 of file HelmholtzFluidProperties.C.

Referenced by CO2FluidProperties::p_from_rho_T(), and HelmholtzFluidProperties::rho_from_p_T().

 {
   // Scale the input density and temperature
   const Real delta = density / criticalDensity();
   const Real tau = criticalTemperature() / temperature;
 
   return _R * density * temperature * delta * dalpha_ddelta(delta, tau) / molarMass();
 }

◆ p_T_from_h_s()

template<typename T >

void SinglePhaseFluidProperties::p_T_from_h_s	(	const T &	h,
		const T &	s,
		Real	p0,
		Real	T0,
		T &	pressure,
		T &	temperature,
		bool &	conversion_succeeded
	)		const

inherited

Determines (p,T) from (h,s) using Newton Solve in 2D Useful for conversion between different sets of state variables.

Parameters

[in]	h	specific enthalpy (J / kg)
[in]	s	specific entropy (J/K*kg)
[in]	p0	initial guess for pressure (Pa / kg)
[in]	T0	initial guess for temperature (K)
[out]	fluid	pressure (Pa / kg)
[out]	Temperature	(K)

Definition at line 575 of file SinglePhaseFluidProperties.h.

Referenced by TabulatedFluidProperties::T_from_h_s().

 {
   auto h_lambda = [&](const T & pressure, const T & temperature, T & new_h, T & dh_dp, T & dh_dT)
   { h_from_p_T(pressure, temperature, new_h, dh_dp, dh_dT); };
   auto s_lambda = [&](const T & pressure, const T & temperature, T & new_s, T & ds_dp, T & ds_dT)
   { s_from_p_T(pressure, temperature, new_s, ds_dp, ds_dT); };
   try
   {
     FluidPropertiesUtils::NewtonSolve2D(
         h, s, p0, T0, pressure, temperature, _tolerance, _tolerance, h_lambda, s_lambda);
     conversion_succeeded = true;
   }
   catch (MooseException &)
   {
     conversion_succeeded = false;
   }
 
   if (!conversion_succeeded)
     mooseDoOnce(mooseWarning("Conversion from (h, s)=(", h, ", ", s, ") to (p, T) failed"));
 }

◆ p_T_from_v_e()

template<typename CppType >

void SinglePhaseFluidProperties::p_T_from_v_e	(	const CppType &	v,
		const CppType &	e,
		Real	p0,
		Real	T0,
		CppType &	p,
		CppType &	T,
		bool &	conversion_succeeded
	)		const

inherited

Determines (p,T) from (v,e) using Newton Solve in 2D Useful for conversion between different sets of state variables.

Parameters

[in]	v	specific volume (m^3 / kg)
[in]	e	specific internal energy (J / kg)
[in]	p0	initial guess for pressure (Pa / kg)
[in]	T0	initial guess for temperature (K)
[out]	fluid	pressure (Pa / kg)
[out]	Temperature	(K)

Definition at line 511 of file SinglePhaseFluidProperties.h.

Referenced by TabulatedBicubicFluidProperties::constructInterpolation(), TemperaturePressureFunctionFluidProperties::cp_from_v_e(), TemperaturePressureFunctionFluidProperties::cv_from_v_e(), TabulatedFluidProperties::g_from_v_e(), TemperaturePressureFunctionFluidProperties::k_from_v_e(), TemperaturePressureFunctionFluidProperties::mu_from_v_e(), and TemperaturePressureFunctionFluidProperties::T_from_v_e().

 {
   auto v_lambda = [&](const CppType & pressure,
                       const CppType & temperature,
                       CppType & new_v,
                       CppType & dv_dp,
                       CppType & dv_dT) { v_from_p_T(pressure, temperature, new_v, dv_dp, dv_dT); };
   auto e_lambda = [&](const CppType & pressure,
                       const CppType & temperature,
                       CppType & new_e,
                       CppType & de_dp,
                       CppType & de_dT) { e_from_p_T(pressure, temperature, new_e, de_dp, de_dT); };
   try
   {
     FluidPropertiesUtils::NewtonSolve2D(
         v, e, p0, T0, p, T, _tolerance, _tolerance, v_lambda, e_lambda);
     conversion_succeeded = true;
   }
   catch (MooseException &)
   {
     conversion_succeeded = false;
   }
 
   if (!conversion_succeeded)
     mooseDoOnce(mooseWarning("Conversion from (v, e)=(", v, ", ", e, ") to (p, T) failed"));
 }

◆ p_T_from_v_h()

template<typename T >

void SinglePhaseFluidProperties::p_T_from_v_h	(	const T &	v,
		const T &	h,
		Real	p0,
		Real	T0,
		T &	pressure,
		T &	temperature,
		bool &	conversion_succeeded
	)		const

inherited

Determines (p,T) from (v,h) using Newton Solve in 2D Useful for conversion between different sets of state variables.

Parameters

[in]	v	specific volume (m^3 / kg)
[in]	h	specific enthalpy (J / kg)
[in]	p0	initial guess for pressure (Pa / kg)
[in]	T0	initial guess for temperature (K)
[out]	fluid	pressure (Pa / kg)
[out]	Temperature	(K)

Definition at line 546 of file SinglePhaseFluidProperties.h.

Referenced by TabulatedBicubicFluidProperties::constructInterpolation().

 {
   auto v_lambda = [&](const T & pressure, const T & temperature, T & new_v, T & dv_dp, T & dv_dT)
   { v_from_p_T(pressure, temperature, new_v, dv_dp, dv_dT); };
   auto h_lambda = [&](const T & pressure, const T & temperature, T & new_h, T & dh_dp, T & dh_dT)
   { h_from_p_T(pressure, temperature, new_h, dh_dp, dh_dT); };
   try
   {
     FluidPropertiesUtils::NewtonSolve2D(
         v, h, p0, T0, pressure, temperature, _tolerance, _tolerance, v_lambda, h_lambda);
     conversion_succeeded = true;
   }
   catch (MooseException &)
   {
     conversion_succeeded = false;
   }
 
   if (!conversion_succeeded)
     mooseDoOnce(mooseWarning("Conversion from (v, h)=(", v, ", ", h, ") to (p, T) failed"));
 }

◆ propfunc() [1/12]

SinglePhaseFluidProperties::propfunc	(	p	,
		v	,
		e
	)

inherited

Compute a fluid property given for the state defined by two given properties.

For all functions, the first two arguments are the given properties that define the fluid state. For the two-argument variants, the desired property is the return value. The five-argument variants also provide partial derivatives dx/da and dx/db where x is the desired property being computed, a is the first given property, and b is the second given property. The desired property, dx/da, and dx/db are stored into the 3rd, 4th, and 5th arguments respectively.

Properties/parameters used in these function are listed below with their units:

p pressure [Pa] T temperature [K] e specific internal energy [J/kg] v specific volume [m^3/kg] rho density [kg/m^3] h specific enthalpy [J/kg] s specific entropy [J/(kg*K)] mu viscosity [Pa*s] k thermal conductivity [W/(m*K)] c speed of sound [m/s] cp constant-pressure specific heat [J/K] cv constant-volume specific heat [J/K] beta volumetric thermal expansion coefficient [1/K] g Gibbs free energy [J] pp_sat partial pressure at saturation [Pa] gamma Adiabatic ratio (cp/cv) [-]

As an example:

// calculate pressure given specific vol and energy: auto pressure = your_fluid_properties_object.p_from_v_e(specific_vol, specific_energy);

// or use the derivative variant: Real dp_dv = 0; // derivative will be stored into here Real dp_de = 0; // derivative will be stored into here your_fluid_properties_object.p_from_v_e(specific_vol, specific_energy, pressure, dp_dv, dp_de);

Automatic differentiation (AD) support is provided through x_from_a_b(ADReal a, ADReal b) and x_from_a_b(ADReal a, ADReal b, ADReal x, ADReal dx_da, ADReal dx_db) versions of the functions where a and b must be ADReal/DualNumber's calculated using all AD-supporting values:

auto v = 1/rho; // rho must be an AD non-linear variable. auto e = rhoE/rho - vel_energy; // rhoE and vel_energy must be AD variables/numbers also. auto pressure = your_fluid_properties_object.p_from_v_e(v, e); // pressure now contains partial derivatives w.r.t. all degrees of freedom

◆ propfunc() [2/12]

e SinglePhaseFluidProperties::propfunc	(	c	,
		v	,
		e
	)

inherited

◆ propfunc() [3/12]

e e SinglePhaseFluidProperties::propfunc	(	cv	,
		v	,
		e
	)

inherited

◆ propfunc() [4/12]

e e e SinglePhaseFluidProperties::propfunc	(	k	,
		v	,
		e
	)

inherited

◆ propfunc() [5/12]

e e e e SinglePhaseFluidProperties::propfunc	(	s	,
		h	,
		p
	)

inherited

◆ propfunc() [6/12]

e e e e s SinglePhaseFluidProperties::propfunc	(	e	,
		v	,
		h
	)

inherited

◆ propfunc() [7/12]

e e e e s T SinglePhaseFluidProperties::propfunc	(	pp_sat	,
		p	,
		T
	)

inherited

◆ propfunc() [8/12]

e e e e s T T SinglePhaseFluidProperties::propfunc	(	k	,
		rho	,
		T
	)

inherited

◆ propfunc() [9/12]

e e e e s T T T T T SinglePhaseFluidProperties::propfunc	(	rho	,
		p	,
		T
	)

inherited

◆ propfunc() [10/12]

e e e e s T T T T T rho SinglePhaseFluidProperties::propfunc	(	e	,
		T	,
		v
	)

inherited

◆ propfunc() [11/12]

e e e e s T T T T T rho v SinglePhaseFluidProperties::propfunc	(	h	,
		T	,
		v
	)

inherited

◆ propfunc() [12/12]

e e e e s T T T T T rho v v SinglePhaseFluidProperties::propfunc	(	cv	,
		T	,
		v
	)

inherited

◆ propfuncWithDefault() [1/7]

e e e e s T T T SinglePhaseFluidProperties::propfuncWithDefault	(	cp	,
		p	,
		T
	)

inherited

◆ propfuncWithDefault() [2/7]

e e e e s T T T T SinglePhaseFluidProperties::propfuncWithDefault	(	mu	,
		p	,
		T
	)

inherited

◆ propfuncWithDefault() [3/7]

e e e e s T T T T T rho v v T SinglePhaseFluidProperties::propfuncWithDefault	(	h	,
		v	,
		e
	)

inherited

◆ propfuncWithDefault() [4/7]

e e e e s T T T T T rho v v T e SinglePhaseFluidProperties::propfuncWithDefault	(	p	,
		h	,
		s
	)

inherited

◆ propfuncWithDefault() [5/7]

e e e e s T T T T T rho v v T e p SinglePhaseFluidProperties::propfuncWithDefault	(	T	,
		p	,
		h
	)

inherited

◆ propfuncWithDefault() [6/7]

e e e e s T T T T T rho v v T e p T SinglePhaseFluidProperties::propfuncWithDefault	(	v	,
		p	,
		T
	)

inherited

◆ propfuncWithDefault() [7/7]

e e e e s T T T T T rho v v T e p T T SinglePhaseFluidProperties::propfuncWithDefault	(	gamma	,
		v	,
		e
	)

inherited

◆ rho_e_from_p_T()

void SinglePhaseFluidProperties::rho_e_from_p_T	(	Real	p,
		Real	T,
		Real &	rho,
		Real &	drho_dp,
		Real &	drho_dT,
		Real &	e,
		Real &	de_dp,
		Real &	de_dT
	)		const

virtualinherited

Definition at line 434 of file SinglePhaseFluidProperties.C.

 {
   rho_from_p_T(p, T, rho, drho_dp, drho_dT);
   e_from_p_T(p, T, e, de_dp, de_dT);
 }

◆ rho_from_p_T() [1/2]

Real HelmholtzFluidProperties::rho_from_p_T	(	Real	pressure,
		Real	temperature
	)		const

overridevirtualinherited

Reimplemented in CO2FluidProperties.

Definition at line 28 of file HelmholtzFluidProperties.C.

Referenced by HelmholtzFluidProperties::c_from_p_T(), HelmholtzFluidProperties::cp_from_p_T(), HelmholtzFluidProperties::cv_from_p_T(), HelmholtzFluidProperties::e_from_p_T(), HelmholtzFluidProperties::h_from_p_T(), NitrogenFluidProperties::k_from_p_T(), k_from_p_T(), NitrogenFluidProperties::mu_from_p_T(), CO2FluidProperties::mu_from_p_T(), mu_from_p_T(), HelmholtzFluidProperties::rho_from_p_T(), CO2FluidProperties::rho_from_p_T(), NitrogenFluidProperties::rho_mu_from_p_T(), rho_mu_from_p_T(), and HelmholtzFluidProperties::s_from_p_T().

 {
   Real density;
   // Initial estimate of a bracketing interval for the density
   Real lower_density = 1.0e-2;
   Real upper_density = 100.0;
 
   // The density is found by finding the zero of the pressure
   auto pressure_diff = [&pressure, &temperature, this](Real x)
   { return this->p_from_rho_T(x, temperature) - pressure; };
 
   BrentsMethod::bracket(pressure_diff, lower_density, upper_density);
   density = BrentsMethod::root(pressure_diff, lower_density, upper_density);
 
   return density;
 }

◆ rho_from_p_T() [2/2]

void HelmholtzFluidProperties::rho_from_p_T	(	Real	pressure,
		Real	temperature,
		Real &	rho,
		Real &	drho_dp,
		Real &	drho_dT
	)		const

overridevirtualinherited

Reimplemented in CO2FluidProperties.

Definition at line 46 of file HelmholtzFluidProperties.C.

 {
   rho = this->rho_from_p_T(pressure, temperature);
 
   // Scale the density and temperature
   const Real delta = rho / criticalDensity();
   const Real tau = criticalTemperature() / temperature;
   const Real da_dd = dalpha_ddelta(delta, tau);
   const Real d2a_dd2 = d2alpha_ddelta2(delta, tau);
 
   drho_dp = molarMass() / (_R * temperature * delta * (2.0 * da_dd + delta * d2a_dd2));
   drho_dT = rho * (tau * d2alpha_ddeltatau(delta, tau) - da_dd) / temperature /
             (2.0 * da_dd + delta * d2a_dd2);
 }

◆ rho_mu_from_p_T() [1/3]

void HydrogenFluidProperties::rho_mu_from_p_T	(	Real	p,
		Real	T,
		Real &	rho,
		Real &	mu
	)		const

overridevirtual

Combined methods.

These methods are particularly useful for the PorousFlow module, where density and viscosity are typically both computed everywhere. The combined methods allow the most efficient means of calculating both properties, especially where rho(p, T) and mu(rho, T). In this case, an extra density calculation would be required to calculate mu(p, T). All property names are described above.

Reimplemented from SinglePhaseFluidProperties.

Definition at line 192 of file HydrogenFluidProperties.C.

 {
   rho = rho_from_p_T(pressure, temperature);
   mu = mu_from_rho_T(rho, temperature);
 }

◆ rho_mu_from_p_T() [2/3]

void HydrogenFluidProperties::rho_mu_from_p_T	(	Real	pressure,
		Real	temperature,
		Real &	rho,
		Real &	drho_dp,
		Real &	drho_dT,
		Real &	mu,
		Real &	dmu_dp,
		Real &	dmu_dT
	)		const

overridevirtual

Reimplemented from SinglePhaseFluidProperties.

Definition at line 202 of file HydrogenFluidProperties.C.

 {
   rho_from_p_T(pressure, temperature, rho, drho_dp, drho_dT);
   Real dmu_drho;
   mu_from_rho_T(rho, temperature, drho_dT, mu, dmu_drho, dmu_dT);
   dmu_dp = dmu_drho * drho_dp;
 }

◆ rho_mu_from_p_T() [3/3]

void SinglePhaseFluidProperties::rho_mu_from_p_T	(	const ADReal &	p,
		const ADReal &	T,
		ADReal &	rho,
		ADReal &	mu
	)		const

virtualinherited

Definition at line 469 of file SinglePhaseFluidProperties.C.

 {
   rho = rho_from_p_T(p, T);
   mu = mu_from_p_T(p, T);
 }

◆ s_from_p_T() [1/2]

Real HelmholtzFluidProperties::s_from_p_T	(	Real	pressure,
		Real	temperature
	)		const

overridevirtualinherited

Definition at line 148 of file HelmholtzFluidProperties.C.

Referenced by HelmholtzFluidProperties::s_from_p_T().

 {
   // Require density first
   const Real density = rho_from_p_T(pressure, temperature);
   // Scale the input density and temperature
   const Real delta = density / criticalDensity();
   const Real tau = criticalTemperature() / temperature;
 
   return _R * (tau * dalpha_dtau(delta, tau) - alpha(delta, tau)) / molarMass();
 }

◆ s_from_p_T() [2/2]

void HelmholtzFluidProperties::s_from_p_T	(	Real	p,
		Real	T,
		Real &	s,
		Real &	ds_dp,
		Real &	ds_dT
	)		const

overridevirtualinherited

Definition at line 160 of file HelmholtzFluidProperties.C.

 {
   s = this->s_from_p_T(pressure, temperature);
 
   // Require density first
   const Real density = rho_from_p_T(pressure, temperature);
   // Scale the input density and temperature
   const Real delta = density / criticalDensity();
   const Real tau = criticalTemperature() / temperature;
 
   const Real da_dd = dalpha_ddelta(delta, tau);
   const Real da_dt = dalpha_dtau(delta, tau);
   const Real d2a_dd2 = d2alpha_ddelta2(delta, tau);
   const Real d2a_dt2 = d2alpha_dtau2(delta, tau);
   const Real d2a_ddt = d2alpha_ddeltatau(delta, tau);
 
   ds_dp = tau * (d2a_ddt - da_dd) / (density * temperature * (2.0 * da_dd + delta * d2a_dd2));
   ds_dT = -_R * tau * (da_dt - alpha(delta, tau) + tau * (d2a_dt2 - da_dt)) /
           (molarMass() * temperature);
 }

◆ subdomainSetup()

virtual void FluidProperties::subdomainSetup ( )

inlinefinalvirtualinherited

Reimplemented from ThreadedGeneralUserObject.

Definition at line 38 of file FluidProperties.h.

38 {}

◆ T_from_p_h()

Real HelmholtzFluidProperties::T_from_p_h	(	Real	pressure,
		Real	enthalpy
	)		const

overridevirtualinherited

Definition at line 220 of file HelmholtzFluidProperties.C.

 {
   auto lambda = [&](Real pressure, Real current_T, Real & new_h, Real & dh_dp, Real & dh_dT)
   { h_from_p_T(pressure, current_T, new_h, dh_dp, dh_dT); };
   Real T = FluidPropertiesUtils::NewtonSolve(
                pressure, enthalpy, _T_initial_guess, _tolerance, lambda, name() + "::T_from_p_h")
                .first;
   // check for nans
   if (std::isnan(T))
     mooseError("Conversion from enthalpy (h = ",
                enthalpy,
                ") and pressure (p = ",
                pressure,
                ") to temperature failed to converge.");
   return T;
 }

◆ threadJoin()

virtual void FluidProperties::threadJoin ( const UserObject & )

inlinefinalvirtualinherited

Reimplemented from ThreadedGeneralUserObject.

Definition at line 37 of file FluidProperties.h.

37 {}

◆ triplePointPressure()

Real HydrogenFluidProperties::triplePointPressure ( ) const

overridevirtual

Triple point pressure.

Returns: triple point pressure (Pa)

Reimplemented from SinglePhaseFluidProperties.

Definition at line 68 of file HydrogenFluidProperties.C.

 {
   return _p_triple;
 }

◆ triplePointTemperature()

Real HydrogenFluidProperties::triplePointTemperature ( ) const

overridevirtual

Triple point temperature.

Returns: triple point temperature (K)

Reimplemented from SinglePhaseFluidProperties.

Definition at line 74 of file HydrogenFluidProperties.C.

 {
   return _T_triple;
 }

◆ v_e_from_p_T() [1/2]

template<typename CppType >

void SinglePhaseFluidProperties::v_e_from_p_T	(	const CppType &	p,
		const CppType &	T,
		CppType &	v,
		CppType &	e
	)		const

inherited

Definition at line 604 of file SinglePhaseFluidProperties.h.

 {
   const CppType rho = rho_from_p_T(p, T);
   v = 1.0 / rho;
   try
   {
     // more likely to not involve a Newton search
     e = e_from_p_T(p, T);
   }
   catch (...)
   {
     e = e_from_p_rho(p, rho);
   }
 }

◆ v_e_from_p_T() [2/2]

template<typename CppType >

void SinglePhaseFluidProperties::v_e_from_p_T	(	const CppType &	p,
		const CppType &	T,
		CppType &	v,
		CppType &	dv_dp,
		CppType &	dv_dT,
		CppType &	e,
		CppType &	de_dp,
		CppType &	de_dT
	)		const

inherited

Definition at line 624 of file SinglePhaseFluidProperties.h.

 {
   CppType rho, drho_dp, drho_dT;
   rho_from_p_T(p, T, rho, drho_dp, drho_dT);
 
   v = 1.0 / rho;
   const CppType dv_drho = -1.0 / (rho * rho);
   dv_dp = dv_drho * drho_dp;
   dv_dT = dv_drho * drho_dT;
 
   CppType de_dp_partial, de_drho;
   e_from_p_rho(p, rho, e, de_dp_partial, de_drho);
   de_dp = de_dp_partial + de_drho * drho_dp;
   de_dT = de_drho * drho_dT;
 }

◆ v_e_spndl_from_T()

void SinglePhaseFluidProperties::v_e_spndl_from_T	(	Real	T,
		Real &	v,
		Real &	e
	)		const

virtualinherited

Specific internal energy from temperature and specific volume.

Parameters

[in]	T	temperature
[in]	v	specific volume

Reimplemented in IdealGasFluidProperties, CaloricallyImperfectGas, and StiffenedGasFluidProperties.

Definition at line 485 of file SinglePhaseFluidProperties.C.

 {
   mooseError(__PRETTY_FUNCTION__, " not implemented.");
 }

◆ validParams()

InputParameters HydrogenFluidProperties::validParams ( )

static

Definition at line 18 of file HydrogenFluidProperties.C.

 {
   InputParameters params = HelmholtzFluidProperties::validParams();
   params.addClassDescription("Fluid properties for Hydrogen (H2)");
   return params;
 }

◆ vaporPressure() [1/3]

Real HydrogenFluidProperties::vaporPressure ( Real T ) const

overridevirtual

Vapor pressure.

Used to delineate liquid and gas phases. Valid for temperatures between the triple point temperature and the critical temperature

Parameters

	T	fluid temperature (K)
[out]	saturation	pressure (Pa)
[out]	derivative	of saturation pressure wrt temperature (Pa/K)

Reimplemented from SinglePhaseFluidProperties.

Definition at line 277 of file HydrogenFluidProperties.C.

 {
   if (temperature < _T_triple || temperature > _T_critical)
     throw MooseException("Temperature is out of range in " + name() + ": vaporPressure()");
 
   const Real Tr = temperature / _T_critical;
   const Real theta = 1.0 - Tr;
 
   const Real logpressure = (-4.89789 * theta + 0.988588 * std::pow(theta, 1.5) +
                             0.349689 * Utility::pow<2>(theta) + 0.499356 * std::pow(theta, 2.85)) /
                            Tr;
 
   return _p_critical * std::exp(logpressure);
 }

◆ vaporPressure() [2/3]

void HydrogenFluidProperties::vaporPressure	(	Real	temperature,
		Real &	psat,
		Real &	dpsat_dT
	)		const

overridevirtual

Reimplemented from SinglePhaseFluidProperties.

Definition at line 293 of file HydrogenFluidProperties.C.

 {
   mooseError("vaporPressure() is not implemented");
 }

◆ vaporPressure() [3/3]

ADReal SinglePhaseFluidProperties::vaporPressure ( const ADReal & T ) const

virtualinherited

Definition at line 389 of file SinglePhaseFluidProperties.C.

 {
   Real p = 0.0;
   Real temperature = T.value();
   Real dpdT = 0.0;
 
   vaporPressure(temperature, p, dpdT);
 
   ADReal result = p;
   result.derivatives() = T.derivatives() * dpdT;
 
   return result;
 }

◆ vaporTemperature() [1/3]

Real SinglePhaseFluidProperties::vaporTemperature ( Real p ) const

virtualinherited

Vapor temperature.

Used to delineate liquid and gas phases. Valid for pressures between the triple point pressure and the critical pressure

Parameters

	p	fluid pressure (Pa)
[out]	saturation	temperature (K)
[out]	derivative	of saturation temperature wrt pressure

Reimplemented in Water97FluidProperties, and TabulatedFluidProperties.

Definition at line 404 of file SinglePhaseFluidProperties.C.

Referenced by PorousFlowWaterVapor::thermophysicalProperties(), TabulatedFluidProperties::vaporTemperature(), and SinglePhaseFluidProperties::vaporTemperature().

 {
   mooseError(__PRETTY_FUNCTION__, " not implemented.");
 }

◆ vaporTemperature() [2/3]

void SinglePhaseFluidProperties::vaporTemperature	(	Real	p,
		Real &	Tsat,
		Real &	dTsat_dp
	)		const

virtualinherited

Reimplemented in Water97FluidProperties, and TabulatedFluidProperties.

Definition at line 410 of file SinglePhaseFluidProperties.C.

 {
   unimplementedDerivativeMethod(__PRETTY_FUNCTION__);
 
   dT_dp = 0.0;
   T = vaporTemperature(p);
 }

◆ vaporTemperature() [3/3]

ADReal SinglePhaseFluidProperties::vaporTemperature ( const ADReal & p ) const

virtualinherited

Definition at line 419 of file SinglePhaseFluidProperties.C.

 {
   Real T = 0.0;
   Real pressure = p.value();
   Real dTdp = 0.0;
 
   vaporTemperature(pressure, T, dTdp);
 
   ADReal result = T;
   result.derivatives() = p.derivatives() * dTdp;
 
   return result;
 }

◆ xyDerivatives()

template<typename T , typename Functor >

void SinglePhaseFluidProperties::xyDerivatives	(	const T	x,
		const T &	y,
		T &	z,
		T &	dz_dx,
		T &	dz_dy,
		const Functor &	z_from_x_y
	)

staticprotectedinherited

Computes the dependent variable z and its derivatives with respect to the independent variables x and y using the simple two parameter z_from_x_y functor.

The derivatives are computed using a compound automatic differentiation type

Definition at line 490 of file SinglePhaseFluidProperties.h.

Referenced by Water97FluidProperties::e_from_p_rho_template(), Water97FluidProperties::h_from_p_T_template(), Water97FluidProperties::rho_from_p_T_template(), and Water97FluidProperties::v_from_p_T_template().

 {
   typedef MetaPhysicL::DualNumber<T, MetaPhysicL::NumberArray<2, T>> CompoundType;
   const auto [zero, one] = makeZeroAndOne(x);
 
   CompoundType x_c(x, zero);
   auto & x_cd = x_c.derivatives();
   x_cd[0] = one;
   CompoundType y_c(y, zero);
   auto & y_cd = y_c.derivatives();
   y_cd[1] = one;
 
   const auto z_c = z_from_x_y(x_c, y_c);
   z = z_c.value();
   dz_dx = z_c.derivatives()[0];
   dz_dy = z_c.derivatives()[1];
 }

Member Data Documentation

◆ _a

const std::array<Real, 5> HydrogenFluidProperties::_a {{1.616, -0.4117, -0.792, 0.758, 1.217}}

protected

Coefficients for ideal gas component of the Helmholtz free energy.

Definition at line 169 of file HydrogenFluidProperties.h.

Referenced by alpha(), d2alpha_dtau2(), and dalpha_dtau().

◆ _a1k

const std::array<Real, 7> HydrogenFluidProperties::_a1k

protected

Initial value:

{

{-3.40976e-1, 4.5882, -1.4508, 3.26394e-1, 3.16939e-3, 1.90592e-4, -1.139e-6}}

Coefficients for thermal conductivity.

Definition at line 196 of file HydrogenFluidProperties.h.

Referenced by k_from_rho_T().

◆ _a2k

const std::array<Real, 4> HydrogenFluidProperties::_a2k {{1.38497e2, -2.21878e1, 4.57151, 1.0}}

protected

Definition at line 198 of file HydrogenFluidProperties.h.

Referenced by k_from_rho_T().

◆ _allow_imperfect_jacobians

const bool FluidProperties::_allow_imperfect_jacobians

protectedinherited

Flag to set unimplemented Jacobian entries to zero.

Definition at line 47 of file FluidProperties.h.

Referenced by VaporMixtureFluidProperties::imperfectJacobianMessage(), and SinglePhaseFluidProperties::unimplementedDerivativeMethod().

◆ _amu

const std::array<Real, 5> HydrogenFluidProperties::_amu {{2.09630e-1, -4.55274e-1, 1.423602e-1, -3.35325e-2, 2.76981e-3}}

protected

Coefficients for viscosity.

Definition at line 190 of file HydrogenFluidProperties.h.

Referenced by mu_from_rho_T().

◆ _b

const std::array<Real, 5> HydrogenFluidProperties::_b

protected

Initial value:

{

{-16.0205159149, -22.6580178006, -60.0090511389, -74.9434303817, -206.9392065168}}

Definition at line 170 of file HydrogenFluidProperties.h.

Referenced by alpha(), d2alpha_dtau2(), and dalpha_dtau().

◆ _b1k

const std::array<Real, 5> HydrogenFluidProperties::_b1k {{3.63081e-2, -2.07629e-2, 3.1481e-2, -1.43097e-2, 1.7498e-3}}

protected

Definition at line 199 of file HydrogenFluidProperties.h.

Referenced by k_from_rho_T().

◆ _b2k

const std::array<Real, 5> HydrogenFluidProperties::_b2k {{1.8337e-3, -8.86716e-3, 1.5826e-2, -1.06283e-2, 2.80673e-3}}

protected

Definition at line 200 of file HydrogenFluidProperties.h.

Referenced by k_from_rho_T().

◆ _beta3

const std::array<Real, 5> HydrogenFluidProperties::_beta3 {{-0.171, -0.2245, -0.1304, -0.2785, -0.3967}}

protected

Definition at line 185 of file HydrogenFluidProperties.h.

Referenced by alpha(), d2alpha_ddelta2(), d2alpha_ddeltatau(), d2alpha_dtau2(), dalpha_ddelta(), and dalpha_dtau().

◆ _bmu

const std::array<Real, 7> HydrogenFluidProperties::_bmu {{-0.187, 2.4871, 3.7151, -11.0972, 9.0965, -3.8292, 0.5166}}

protected

Definition at line 191 of file HydrogenFluidProperties.h.

Referenced by mu_from_rho_T().

◆ _cmu

const std::array<Real, 6> HydrogenFluidProperties::_cmu

protected

Initial value:

{

{6.43449673, 4.56334068e-2, 2.32797868e-1, 9.5832612e-1, 1.27941189e-1, 3.63576595e-1}}

Definition at line 192 of file HydrogenFluidProperties.h.

Referenced by mu_from_rho_T().

◆ _d1

const std::array<unsigned int, 7> HydrogenFluidProperties::_d1 {{1, 4, 1, 1, 2, 2, 3}}

protected

Definition at line 175 of file HydrogenFluidProperties.h.

Referenced by alpha(), d2alpha_ddelta2(), d2alpha_ddeltatau(), d2alpha_dtau2(), dalpha_ddelta(), and dalpha_dtau().

◆ _d2

const std::array<unsigned int, 2> HydrogenFluidProperties::_d2 {{1, 3}}

protected

Definition at line 179 of file HydrogenFluidProperties.h.

Referenced by alpha(), d2alpha_ddelta2(), d2alpha_ddeltatau(), d2alpha_dtau2(), dalpha_ddelta(), and dalpha_dtau().

◆ _d3

const std::array<unsigned int, 5> HydrogenFluidProperties::_d3 {{2, 1, 3, 1, 1}}

protected

Definition at line 183 of file HydrogenFluidProperties.h.

Referenced by alpha(), d2alpha_ddelta2(), d2alpha_ddeltatau(), d2alpha_dtau2(), dalpha_ddelta(), and dalpha_dtau().

◆ _D3

const std::array<Real, 5> HydrogenFluidProperties::_D3 {{1.506, 0.156, 1.736, 0.67, 1.662}}

protected

Definition at line 187 of file HydrogenFluidProperties.h.

Referenced by alpha(), d2alpha_ddelta2(), d2alpha_ddeltatau(), d2alpha_dtau2(), dalpha_ddelta(), and dalpha_dtau().

◆ _gamma3

const std::array<Real, 5> HydrogenFluidProperties::_gamma3 {{0.7164, 1.3444, 1.4517, 0.7204, 1.5445}}

protected

Definition at line 186 of file HydrogenFluidProperties.h.

Referenced by alpha(), d2alpha_ddelta2(), d2alpha_ddeltatau(), d2alpha_dtau2(), dalpha_ddelta(), and dalpha_dtau().

◆ _max_newton_its

const unsigned int SinglePhaseFluidProperties::_max_newton_its

protectedinherited

Maximum number of iterations for the variable conversion newton solves.

Definition at line 449 of file SinglePhaseFluidProperties.h.

Referenced by TabulatedFluidProperties::e_from_v_h(), TemperaturePressureFunctionFluidProperties::T_from_p_h(), SimpleFluidProperties::T_from_p_h(), NaKFluidProperties::T_from_p_rho(), TabulatedFluidProperties::T_from_p_rho(), and TabulatedFluidProperties::T_from_p_s().

◆ _Mh2

const Real HydrogenFluidProperties::_Mh2

protected

Hydrogen molar mass (kg/mol)

Definition at line 154 of file HydrogenFluidProperties.h.

Referenced by molarMass(), and mu_from_rho_T().

◆ _N1

const std::array<Real, 7> HydrogenFluidProperties::_N1 {{-6.93643, 0.01, 2.1101, 4.52059, 0.732564, -1.34086, 0.130985}}

protected

Coefficients for residual component of the Helmholtz free energy.

Definition at line 173 of file HydrogenFluidProperties.h.

Referenced by alpha(), d2alpha_ddelta2(), d2alpha_ddeltatau(), d2alpha_dtau2(), dalpha_ddelta(), and dalpha_dtau().

◆ _N2

const std::array<Real, 2> HydrogenFluidProperties::_N2 {{-0.777414, 0.351944}}

protected

Definition at line 177 of file HydrogenFluidProperties.h.

Referenced by alpha(), d2alpha_ddelta2(), d2alpha_ddeltatau(), d2alpha_dtau2(), dalpha_ddelta(), and dalpha_dtau().

◆ _N3

const std::array<Real, 5> HydrogenFluidProperties::_N3 {{-0.0211716, 0.0226312, 0.032187, -0.0231752, 0.0557346}}

protected

Definition at line 181 of file HydrogenFluidProperties.h.

Referenced by alpha(), d2alpha_ddelta2(), d2alpha_ddeltatau(), d2alpha_dtau2(), dalpha_ddelta(), and dalpha_dtau().

◆ _p_critical

const Real HydrogenFluidProperties::_p_critical

protected

Critical pressure (Pa)

Definition at line 156 of file HydrogenFluidProperties.h.

Referenced by criticalPressure(), and vaporPressure().

◆ _p_initial_guess

const Real SinglePhaseFluidProperties::_p_initial_guess

protectedinherited

Initial guess for pressure (or pressure used to compute the initial guess)

Definition at line 447 of file SinglePhaseFluidProperties.h.

Referenced by TabulatedFluidProperties::checkInitialGuess(), TabulatedBicubicFluidProperties::constructInterpolation(), TemperaturePressureFunctionFluidProperties::cp_from_v_e(), TemperaturePressureFunctionFluidProperties::cv_from_v_e(), TabulatedFluidProperties::e_from_v_h(), TabulatedFluidProperties::g_from_v_e(), TemperaturePressureFunctionFluidProperties::k_from_v_e(), TemperaturePressureFunctionFluidProperties::mu_from_v_e(), TemperaturePressureFunctionFluidProperties::p_from_v_e(), Water97FluidProperties::p_from_v_e_template(), Water97FluidProperties::p_T_from_v_h(), TabulatedFluidProperties::T_from_h_s(), and TemperaturePressureFunctionFluidProperties::T_from_v_e().

◆ _p_triple

const Real HydrogenFluidProperties::_p_triple

protected

Triple point pressure (Pa)

Definition at line 164 of file HydrogenFluidProperties.h.

Referenced by triplePointPressure().

◆ _phi3

const std::array<Real, 5> HydrogenFluidProperties::_phi3 {{-1.685, -0.489, -0.103, -2.506, -1.607}}

protected

Definition at line 184 of file HydrogenFluidProperties.h.

Referenced by alpha(), d2alpha_ddelta2(), d2alpha_ddeltatau(), d2alpha_dtau2(), dalpha_ddelta(), and dalpha_dtau().

◆ _R

const Real FluidProperties::_R = 8.3144598

staticinherited

Universal gas constant (J/mol/K)

Definition at line 41 of file FluidProperties.h.

Referenced by HelmholtzFluidProperties::c_from_p_T(), IdealGasFluidProperties::c_from_p_T(), MassDiffusionEnergyGasMixDGKernel::computeComponentTotalEnthalpy(), BinaryDiffusionCoefMaterial::computeQpProperties(), HelmholtzFluidProperties::cp_from_p_T(), HelmholtzFluidProperties::cv_from_p_T(), HelmholtzFluidProperties::e_from_p_T(), HelmholtzFluidProperties::h_from_p_T(), LeadLithiumFluidProperties::mu_from_p_T(), LeadLithiumFluidProperties::mu_from_v_e(), HelmholtzFluidProperties::p_from_rho_T(), HelmholtzFluidProperties::rho_from_p_T(), CaloricallyImperfectGas::rho_from_p_T(), IdealGasFluidProperties::rho_from_p_T(), HelmholtzFluidProperties::s_from_p_T(), and TEST_F().

◆ _rho_critical

const Real HydrogenFluidProperties::_rho_critical

protected

Critical density (kg/m^3)

Definition at line 162 of file HydrogenFluidProperties.h.

Referenced by criticalDensity().

◆ _rho_molar_critical

const Real HydrogenFluidProperties::_rho_molar_critical

protected

Critical molar density (mol/l)

Definition at line 160 of file HydrogenFluidProperties.h.

◆ _t1

const std::array<Real, 7> HydrogenFluidProperties::_t1 {{0.6844, 1.0, 0.989, 0.489, 0.803, 1.1444, 1.409}}

protected

Definition at line 174 of file HydrogenFluidProperties.h.

Referenced by alpha(), d2alpha_ddelta2(), d2alpha_ddeltatau(), d2alpha_dtau2(), dalpha_ddelta(), and dalpha_dtau().

◆ _t2

const std::array<Real, 2> HydrogenFluidProperties::_t2 {{1.754, 1.311}}

protected

Definition at line 178 of file HydrogenFluidProperties.h.

Referenced by alpha(), d2alpha_ddelta2(), d2alpha_ddeltatau(), d2alpha_dtau2(), dalpha_ddelta(), and dalpha_dtau().

◆ _t3

const std::array<Real, 5> HydrogenFluidProperties::_t3 {{4.187, 5.646, 0.791, 7.249, 2.986}}

protected

Definition at line 182 of file HydrogenFluidProperties.h.

Referenced by alpha(), d2alpha_ddelta2(), d2alpha_ddeltatau(), d2alpha_dtau2(), dalpha_ddelta(), and dalpha_dtau().

◆ _T_c2k

const Real FluidProperties::_T_c2k

protectedinherited

Conversion of temperature from Celsius to Kelvin.

Definition at line 45 of file FluidProperties.h.

Referenced by NaKFluidProperties::cp_from_p_T(), NaClFluidProperties::cp_from_p_T(), BrineFluidProperties::cp_from_p_T_X(), NaKFluidProperties::h_from_p_T(), NaClFluidProperties::h_from_p_T(), BrineFluidProperties::h_from_p_T_X(), BrineFluidProperties::haliteSolubility(), NaKFluidProperties::k_from_p_T(), NaClFluidProperties::k_from_p_T(), BrineFluidProperties::k_from_p_T_X(), NaKFluidProperties::mu_from_p_T(), BrineFluidProperties::mu_from_p_T_X(), CO2FluidProperties::partialDensity(), NaKFluidProperties::rho_from_p_T(), NaClFluidProperties::rho_from_p_T(), BrineFluidProperties::rho_from_p_T_X(), and NaKFluidProperties::T_from_p_h().

◆ _T_critical

const Real HydrogenFluidProperties::_T_critical

protected

Critical temperature (K)

Definition at line 158 of file HydrogenFluidProperties.h.

Referenced by criticalTemperature(), mu_from_rho_T(), and vaporPressure().

◆ _T_initial_guess

const Real SinglePhaseFluidProperties::_T_initial_guess

protectedinherited

Initial guess for temperature (or temperature used to compute the initial guess)

Definition at line 445 of file SinglePhaseFluidProperties.h.

◆ _T_triple

const Real HydrogenFluidProperties::_T_triple

protected

Triple point temperature (K)

Definition at line 166 of file HydrogenFluidProperties.h.

Referenced by triplePointTemperature().

◆ _tolerance

const Real SinglePhaseFluidProperties::_tolerance

protectedinherited

Newton's method may be used to convert between variable sets.

Relative tolerance of the solves

Definition at line 443 of file SinglePhaseFluidProperties.h.

◆ h

e e e e s T T T T T rho v v T e SinglePhaseFluidProperties::h

inherited

Definition at line 212 of file SinglePhaseFluidProperties.h.

◆ p [1/10]

e e e e SinglePhaseFluidProperties::p

inherited

Definition at line 190 of file SinglePhaseFluidProperties.h.

◆ p [2/10]

e e e e s SinglePhaseFluidProperties::p

inherited

Definition at line 192 of file SinglePhaseFluidProperties.h.

◆ p [3/10]

e e e e s T T SinglePhaseFluidProperties::p

inherited

Definition at line 196 of file SinglePhaseFluidProperties.h.

◆ p [4/10]

e e e e s T T T SinglePhaseFluidProperties::p

inherited

Definition at line 198 of file SinglePhaseFluidProperties.h.

◆ p [5/10]

e e e e s T T T T SinglePhaseFluidProperties::p

inherited

Definition at line 200 of file SinglePhaseFluidProperties.h.

◆ p [6/10]

e e e e s T T T T T SinglePhaseFluidProperties::p

inherited

Definition at line 202 of file SinglePhaseFluidProperties.h.

◆ p [7/10]

e e e e s T T T T T rho v v SinglePhaseFluidProperties::p

inherited

Definition at line 208 of file SinglePhaseFluidProperties.h.

◆ p [8/10]

e e e e s T T T T T rho v v T e p SinglePhaseFluidProperties::p

inherited

Definition at line 214 of file SinglePhaseFluidProperties.h.

◆ p [9/10]

e e e e s T T T T T rho v v T e p T SinglePhaseFluidProperties::p

inherited

Definition at line 216 of file SinglePhaseFluidProperties.h.

◆ p [10/10]

e e e e s T T T T T rho v v T e p T T SinglePhaseFluidProperties::p

inherited

Definition at line 218 of file SinglePhaseFluidProperties.h.

◆ rho

e e e e s T SinglePhaseFluidProperties::rho

inherited

Definition at line 194 of file SinglePhaseFluidProperties.h.

◆ T [1/2]

e e e e s T T T T T rho SinglePhaseFluidProperties::T

inherited

Definition at line 204 of file SinglePhaseFluidProperties.h.

◆ T [2/2]

e e e e s T T T T T rho v SinglePhaseFluidProperties::T

inherited

Definition at line 206 of file SinglePhaseFluidProperties.h.

◆ v [1/5]

SinglePhaseFluidProperties::v

inherited

Definition at line 182 of file SinglePhaseFluidProperties.h.

◆ v [2/5]

e SinglePhaseFluidProperties::v

inherited

Definition at line 184 of file SinglePhaseFluidProperties.h.

◆ v [3/5]

e e SinglePhaseFluidProperties::v

inherited

Definition at line 186 of file SinglePhaseFluidProperties.h.

◆ v [4/5]

e e e SinglePhaseFluidProperties::v

inherited

Definition at line 188 of file SinglePhaseFluidProperties.h.

◆ v [5/5]

e e e e s T T T T T rho v v T SinglePhaseFluidProperties::v

inherited

Definition at line 210 of file SinglePhaseFluidProperties.h.

The documentation for this class was generated from the following files:

fluid_properties/include/fluidproperties/HydrogenFluidProperties.h
fluid_properties/src/fluidproperties/HydrogenFluidProperties.C

Public Types

Public Member Functions

Static Public Member Functions

Public Attributes

Static Public Attributes

Protected Member Functions

Static Protected Member Functions

Protected Attributes

Static Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ HydrogenFluidProperties()

Member Function Documentation

◆ alpha()

◆ c_from_p_T()

◆ cp_from_p_T()

◆ criticalDensity()

◆ criticalInternalEnergy()

◆ criticalPressure()

◆ criticalTemperature()

◆ cv_from_p_T()

◆ d2alpha_ddelta2()

◆ d2alpha_ddeltatau()

◆ d2alpha_dtau2()

◆ dalpha_ddelta()

◆ dalpha_dtau()

◆ e_from_p_T() [1/2]

◆ e_from_p_T() [2/2]

◆ e_spndl_from_v()

◆ execute()

◆ finalize()

◆ fluidName()

◆ h_from_p_T() [1/2]

◆ h_from_p_T() [2/2]

◆ henryCoefficients()

◆ initialize()

◆ k_from_p_T() [1/2]

◆ k_from_p_T() [2/2]

◆ k_from_rho_T()

◆ makeZeroAndOne() [1/2]

◆ makeZeroAndOne() [2/2]

◆ molarMass()

◆ mu_from_p_T() [1/2]

◆ mu_from_p_T() [2/2]

◆ mu_from_rho_T() [1/2]

◆ mu_from_rho_T() [2/2]

◆ p_from_rho_T()

◆ p_T_from_h_s()

◆ p_T_from_v_e()

◆ p_T_from_v_h()

◆ propfunc() [1/12]

◆ propfunc() [2/12]

◆ propfunc() [3/12]

◆ propfunc() [4/12]

◆ propfunc() [5/12]

◆ propfunc() [6/12]

◆ propfunc() [7/12]

◆ propfunc() [8/12]

◆ propfunc() [9/12]

◆ propfunc() [10/12]

◆ propfunc() [11/12]

◆ propfunc() [12/12]

◆ propfuncWithDefault() [1/7]

◆ propfuncWithDefault() [2/7]

◆ propfuncWithDefault() [3/7]

◆ propfuncWithDefault() [4/7]

◆ propfuncWithDefault() [5/7]

◆ propfuncWithDefault() [6/7]

◆ propfuncWithDefault() [7/7]

◆ rho_e_from_p_T()

◆ rho_from_p_T() [1/2]

◆ rho_from_p_T() [2/2]

◆ rho_mu_from_p_T() [1/3]

◆ rho_mu_from_p_T() [2/3]

◆ rho_mu_from_p_T() [3/3]

◆ s_from_p_T() [1/2]

◆ s_from_p_T() [2/2]

◆ subdomainSetup()

◆ T_from_p_h()

◆ threadJoin()