A calorically imperfect gas fluid property class This fluid property assumes that internal energy is a general monotonic function of temperature; behaves like an ideal gas otherwise. More...

#include <CaloricallyImperfectGas.h>

Inheritance diagram for CaloricallyImperfectGas:

[legend]

Public Types
typedef DataFileName	DataFileParameterType

Public Member Functions
	CaloricallyImperfectGas (const InputParameters &parameters)

virtual void	initialSetup () override

virtual Real	p_from_v_e (Real v, Real e) const override

virtual ADReal	p_from_v_e (const ADReal &v, const ADReal &e) const override

virtual void	p_from_v_e (Real v, Real e, Real &p, Real &dp_dv, Real &dp_de) const override

virtual Real	T_from_v_e (Real v, Real e) const override

virtual ADReal	T_from_v_e (const ADReal &v, const ADReal &e) const override

virtual void	T_from_v_e (Real v, Real e, Real &T, Real &dT_dv, Real &dT_de) const override

virtual Real	cp_from_v_e (Real v, Real e) const override

virtual void	cp_from_v_e (Real v, Real e, Real &cp, Real &dcp_dv, Real &dcp_de) const override

virtual Real	cv_from_v_e (Real v, Real e) const override

virtual ADReal	cv_from_v_e (ADReal v, ADReal e) const

virtual void	cv_from_v_e (Real v, Real e, Real &cv, Real &dcv_dv, Real &dcv_de) const override

virtual Real	mu_from_v_e (Real v, Real e) const override

virtual void	mu_from_v_e (Real v, Real e, Real &mu, Real &dmu_dv, Real &dmu_de) const override

virtual Real	k_from_v_e (Real v, Real e) const override

virtual void	k_from_v_e (Real v, Real e, Real &k, Real &dk_dv, Real &dk_de) const override

virtual Real	s_from_v_e (Real v, Real e) const override

virtual void	s_from_v_e (Real v, Real e, Real &s, Real &ds_dv, Real &ds_de) const override

virtual Real	s_from_p_T (Real p, Real T) const override

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

virtual Real	s_from_h_p (Real h, Real p) const override

virtual void	s_from_h_p (Real h, Real p, Real &s, Real &ds_dh, Real &ds_dp) const override

virtual Real	rho_from_p_s (Real p, Real s) const override

virtual Real	e_from_v_h (Real v, Real h) const override

virtual void	e_from_v_h (Real v, Real h, Real &e, Real &de_dv, Real &de_dh) const override

virtual Real	rho_from_p_T (Real p, Real T) const override

virtual ADReal	rho_from_p_T (const ADReal &p, const ADReal &T) const override

virtual void	rho_from_p_T (Real p, Real T, Real &rho, Real &drho_dp, Real &drho_dT) const override

virtual void	rho_from_p_T (const ADReal &p, const ADReal &T, ADReal &rho, ADReal &drho_dp, ADReal &drho_dT) const override

	propfuncWithDefinitionOverride (e, p, rho)

virtual Real	e_from_T_v (Real T, Real v) const override

virtual void	e_from_T_v (Real T, Real v, Real &e, Real &de_dT, Real &de_dv) const override

virtual ADReal	e_from_T_v (const ADReal &T, const ADReal &v) const override

virtual Real	p_from_T_v (Real T, Real v) const override

virtual void	p_from_T_v (Real T, Real v, Real &p, Real &dp_dT, Real &dp_dv) const override

virtual Real	h_from_T_v (Real T, Real v) const override

virtual void	h_from_T_v (Real T, Real v, Real &h, Real &dh_dT, Real &dh_dv) const override

virtual Real	s_from_T_v (Real T, Real v) const override

virtual void	s_from_T_v (Real T, Real v, Real &s, Real &ds_dT, Real &ds_dv) const override

virtual Real	cv_from_T_v (Real T, Real v) const override

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

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

virtual Real	h_from_p_T (Real p, Real T) const override

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

virtual Real	e_from_p_T (Real p, Real T) const override

virtual ADReal	e_from_p_T (ADReal p, ADReal T) const

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

virtual Real	T_from_p_h (Real p, Real h) const override

virtual void	T_from_p_h (Real p, Real h, Real &T, Real &dT_dp, Real &dT_dh) const override

virtual Real	cv_from_p_T (Real p, Real T) const override

virtual void	cv_from_p_T (Real p, Real T, Real &cv, Real &dcv_dp, Real &dcv_dT) const override

virtual Real	cp_from_p_T (Real p, Real T) const override

virtual void	cp_from_p_T (Real p, Real T, Real &cp, Real &dcp_dp, Real &dcp_dT) const override

virtual Real	mu_from_p_T (Real p, Real T) const override

virtual void	mu_from_p_T (Real p, Real T, Real &mu, Real &dmu_dp, Real &dmu_dT) 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::string	fluidName () const override
	Fluid name.

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

virtual Real	criticalTemperature () const override
	Critical temperature.

virtual Real	criticalDensity () const override
	Critical density.

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

virtual Real	gamma_from_v_e (Real v, Real e) const override

virtual ADReal	gamma_from_v_e (ADReal v, ADReal e) const

virtual void	gamma_from_v_e (Real v, Real e, Real &gamma, Real &dgamma_dv, Real &dgamma_de) const override

virtual Real	gamma_from_p_T (Real p, Real T) const override

virtual ADReal	gamma_from_p_T (ADReal p, ADReal T) const

virtual void	gamma_from_p_T (Real p, Real T, Real &gamma, Real &dgamma_dp, Real &dgamma_dT) const override

virtual Real	g_from_v_e (Real v, Real e) const override

virtual Real	p_from_h_s (Real h, Real s) const override

virtual void	p_from_h_s (Real h, Real s, Real &p, Real &dp_dh, Real &dp_ds) const override

virtual Real	c_from_v_e (Real v, Real e) const override

virtual ADReal	c_from_v_e (const ADReal &v, const ADReal &e) const override

virtual void	c_from_v_e (Real v, Real e, Real &c, Real &dc_dv, Real &dc_de) const override

virtual Real	c_from_p_T (Real p, Real T) const override

virtual ADReal	c_from_p_T (const ADReal &p, const ADReal &T) const override

virtual void	c_from_p_T (Real, Real T, Real &c, Real &dc_dp, Real &dc_dT) const override

virtual Real	criticalPressure () const
	Critical pressure.

virtual Real	triplePointPressure () const
	Triple point pressure.

virtual Real	triplePointTemperature () const
	Triple point temperature.

virtual Real	vaporPressure (Real T) const
	Vapor pressure.

virtual void	vaporPressure (Real T, Real &psat, Real &dpsat_dT) const

virtual ADReal	vaporPressure (const ADReal &T) const

virtual Real	vaporTemperature (Real p) const
	Vapor temperature.

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

virtual ADReal	vaporTemperature (const ADReal &p) const

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

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 (Real p, Real T, Real &rho, Real &mu) const
	Combined methods.

virtual void	rho_mu_from_p_T (Real p, Real T, Real &rho, Real &drho_dp, Real &drho_dT, Real &mu, Real &dmu_dp, Real &dmu_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.

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.

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.

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

virtual Real	spatialValue (const Point &) const

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

void	setPrimaryThreadCopy (UserObject *primary)

UserObject *	primaryThreadCopy ()

SubProblem &	getSubProblem () const

bool	shouldDuplicateInitialExecution () 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)

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

bool	isKokkosObject () const

MooseApp &	getMooseApp () const

const std::string &	type () const

const std::string &	name () const

std::string	typeAndName () const

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

MooseObjectName	uniqueName () const

const InputParameters &	parameters () const

const hit::Node *	getHitNode () const

bool	hasBase () const

const std::string &	getBase () 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 &name) const

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

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

void	paramWarning (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

std::string	messagePrefix (const bool hit_prefix=true) const

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

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

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

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

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

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

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

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

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

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

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

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

void	callMooseError (std::string msg, const bool with_prefix, const hit::Node *node=nullptr, const bool show_trace=true) 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	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 UserObjectBase &	getUserObjectBase (const std::string &param_name, bool is_dependency=true) const

const UserObjectBase &	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)

Moose::Kokkos::MaterialProperty< T, dimension >	getKokkosMaterialPropertyByName (const std::string &prop_name_in)

Moose::Kokkos::MaterialProperty< T, dimension >	getKokkosMaterialPropertyOldByName (const std::string &prop_name)

Moose::Kokkos::MaterialProperty< T, dimension >	getKokkosMaterialPropertyOlderByName (const std::string &prop_name)

Moose::Kokkos::MaterialProperty< T, dimension >	getKokkosMaterialProperty (const std::string &name)

Moose::Kokkos::MaterialProperty< T, dimension >	getKokkosMaterialPropertyOld (const std::string &name)

Moose::Kokkos::MaterialProperty< T, dimension >	getKokkosMaterialPropertyOlder (const std::string &name)

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

std::pair< Moose::Kokkos::MaterialProperty< T, dimension >, std::set< SubdomainID > >	getKokkosBlockMaterialProperty (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)

virtual 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	hasKokkosMaterialProperty (const std::string &name)

bool	hasKokkosMaterialProperty (const std::string &name)

bool	hasKokkosMaterialPropertyByName (const std::string &name)

bool	hasKokkosMaterialPropertyByName (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

Moose::Kokkos::Function	getKokkosFunction (const std::string &name) const

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

Moose::Kokkos::Function	getKokkosFunctionByName (const FunctionName &name) const

const T &	getKokkosFunctionByName (const FunctionName &name) const

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

bool	hasKokkosFunctionByName (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	callMooseError (MooseApp const app, const InputParameters &params, std::string msg, const bool with_prefix, const hit::Node node, const bool show_trace=true)

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, NameFunc &&name_func)

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

Public Attributes
	usingCombinedWarningSolutionWarnings

const ConsoleStream	_console

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

static const std::string	type_param

static const std::string	name_param

static const std::string	unique_name_param

static const std::string	app_param

static const std::string	moose_base_param

static const std::string	kokkos_object_param

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 Types
enum	NaNMessage { NAN_MESSAGE_NONE = 0 , NAN_MESSAGE_WARNING = 1 , NAN_MESSAGE_EXCEPTION = 2 , NAN_MESSAGE_ERROR = 3 }

Protected Member Functions
void	setupLookupTables ()
	sets up the T(e) reverse lookup table

void	outOfBounds (const std::string &function, Real value, Real min, Real max) const
	function that handles exceeding parameter limits

virtual void	addPostprocessorDependencyHelper (const PostprocessorName &name) const override

virtual void	addVectorPostprocessorDependencyHelper (const VectorPostprocessorName &name) const override

virtual void	addUserObjectDependencyHelper (const UserObjectBase &uo) const override

void	addReporterDependencyHelper (const ReporterName &reporter_name) override

void	flagInvalidSolutionInternal (const InvalidSolutionID invalid_solution_id) const

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

const ReporterContextBase &	getReporterContextBaseByName (const ReporterName &reporter_name) const

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)

virtual void	getKokkosMaterialPropertyHook (const std::string &, const unsigned int)

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)

Real	getNaN () const
	Throws an error or returns a NaN with or without a warning, with a default message.

template<typename... Args>
Real	getNaN (Args &&... args) const
	Throws an error or returns a NaN with or without a warning.

std::vector< Real >	getNaNVector (const unsigned int &n) const
	Throws an error or returns NaNs with or without a warning, with a default message.

template<typename... Args>
std::vector< Real >	getNaNVector (const unsigned int &n, Args &&... args) const
	Throws an error or returns NaNs with or without a warning.


Real	e_from_T (Real T) const
	helper functions for e(T) and h(T), T(e), T(h), cv(T), cp(T)

Real	h_from_T (Real T) const

Real	T_from_e (Real e) const

Real	T_from_h (Real h) const

template<typename CppType >
Real	cv_from_T (const CppType &T) const

Real	cp_from_T (Real T) const

void	cv_from_T (Real T, Real &cv, Real &dcv_dT) const

void	cp_from_T (Real T, Real &cp, Real &dcp_dT) const

Real	Z_from_T (Real T) 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.

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)

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

Protected Attributes
const Real &	_molar_mass
	molar mass

const Real	_R_specific
	Specific gas constant (R / molar mass)

Real	_T_c

Real	_rho_c

Real	_e_c

Real	_delta_T
	temperature interval in lookup tables

const bool	_out_of_bound_error
	Flag to error if out of bounds.

Real	_tol

const Function *	_e_T
	Internal energy as a function of temperature.

const Function *	_mu_T
	Dynamic viscosity.

const Function *	_k_T
	Thermal conductivity.

std::vector< Real >	_Z_T_lookup
	Z(T) lookup table on uniform grid between _min_temperature and _max_temperature.

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

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

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

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

const bool	_verbose_newton
	Whether to output information about newton solves to console.

const Real	_T_c2k
	Conversion of temperature from Celsius to Kelvin.

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

const Moose::CoordinateSystemType &	_coord_sys

const THREAD_ID	_tid

SubProblem &	_subproblem

FEProblemBase &	_fe_problem

SystemBase &	_sys

Assembly &	_assembly

const bool	_duplicate_initial_execution

std::set< std::string >	_depend_uo

const bool &	_enabled

MooseApp &	_app

Factory &	_factory

ActionFactory &	_action_factory

const std::string &	_type

const std::string &	_name

const InputParameters &	_pars

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 bool	_is_kokkos_object

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

const MooseObject *const	_moose_object

enum NaNMessage	_emit_on_nan
	Raise mooseWarning or mooseError?


Real	_min_temperature
	temperature limits when creating lookup tables

Real	_max_temperature


Real	_min_e
	internal energy and enthalpy limits when creating lookup tables

Real	_max_e

Real	_min_h

Real	_max_h


Real	_delta_e
	step size in internal energy and enthalpy

Real	_delta_h


std::vector< Real >	_T_h_lookup
	inverse lookup table data

std::vector< Real >	_T_e_lookup

Static Protected Attributes
static const std::string	_interpolated_old

static const std::string	_interpolated_older

Private Member Functions
void	unimplementedDerivativeMethod (const std::string &property_function_name) const

const Moose::FunctionBase *	getKokkosFunctionByNameHelper (const FunctionName &name) const

const UserObjectBase &	getUserObjectFromFEProblem (const UserObjectName &object_name, const THREAD_ID tid=0) const

const T &	castUserObject (const UserObjectBase &uo_base, const std::string &param_name="") const

void	mooseObjectError (const std::string &param_name, std::stringstream &oss) const

const std::string &	userObjectType (const UserObjectBase &uo) const

const std::string &	userObjectName (const UserObjectBase &uo) const

const PostprocessorName &	getPostprocessorNameInternal (const std::string &param_name, const unsigned int index, const bool allow_default_value=true) const

bool	isDefaultPostprocessorValueByName (const PostprocessorName &name) const

PostprocessorValue	getDefaultPostprocessorValueByName (const PostprocessorName &name) const

void	checkParam (const std::string &param_name, const unsigned int index=std::numeric_limits< unsigned int >::max()) const

bool	postprocessorsAdded () const

const VectorPostprocessorValue &	getVectorPostprocessorByNameHelper (const VectorPostprocessorName &name, const std::string &vector_name, bool broadcast, std::size_t t_index) const

const VectorPostprocessorContext< VectorPostprocessorValue > &	getVectorPostprocessorContextByNameHelper (const VectorPostprocessorName &name, const std::string &vector_name) const

bool	vectorPostprocessorsAdded () const

bool	reportersAdded () const

void	possiblyCheckHasReporter (const ReporterName &reporter_name, const std::string &param_name="") const

RestartableDataValue &	registerRestartableDataOnApp (std::unique_ptr< RestartableDataValue > data, THREAD_ID tid) const

void	registerRestartableNameWithFilterOnApp (const std::string &name, Moose::RESTARTABLE_FILTER filter)

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

virtual std::string	meshPropertyPrefix (const std::string &data_name) const

const RestartableDataValue &	getMeshPropertyInternal (const std::string &data_name, const std::string &prefix) const

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

const VariableValue *	getDefaultValue (const std::string &var_name) const

const ADVariableValue *	getADDefaultValue (const std::string &var_name) const

void	checkVar (const std::string &var_name) const

void	validateExecutionerType (const std::string &name, const std::string &fn_name) const

Moose::MaterialDataType	getMaterialDataType (const std::set< BoundaryID > &boundary_ids) const

unsigned int	getMaxQps () const

void	addConsumedPropertyName (const MooseObjectName &obj_name, const std::string &prop_name)

const PostprocessorValue &	getPostprocessorValueInternal (const std::string &param_name, unsigned int index, std::size_t t_index) const

const PostprocessorValue &	getPostprocessorValueInternal (const std::string &param_name, unsigned int index, std::size_t t_index) const

const PostprocessorValue &	getPostprocessorValueByNameInternal (const PostprocessorName &name, std::size_t t_index) const

const PostprocessorValue &	getPostprocessorValueByNameInternal (const PostprocessorName &name, std::size_t t_index) const

void	possiblyCheckHasVectorPostprocessor (const std::string &param_name, const std::string &vector_name) const

void	possiblyCheckHasVectorPostprocessor (const std::string &param_name, const std::string &vector_name) const

void	possiblyCheckHasVectorPostprocessorByName (const VectorPostprocessorName &name, const std::string &vector_name) const

void	possiblyCheckHasVectorPostprocessorByName (const VectorPostprocessorName &name, const std::string &vector_name) const

Static Private Member Functions
static const hit::Node *	getHitNode (const InputParameters &params)

static std::string	messagePrefix (const InputParameters &params, const bool hit_prefix)

Private Attributes
UserObject *	_primary_thread_copy

std::set< std::string >	_supplied_uo

const ParallelParamObject &	_parent

const MooseBase &	_si_moose_base

const FEProblemBase *	_si_problem

ExecFlagEnum	_empty_execute_enum

const MooseObject &	_fni_object

const InputParameters &	_fni_params

FEProblemBase &	_fni_feproblem

const THREAD_ID	_fni_tid

const MooseObject &	_uoi_moose_object

const FEProblemBase &	_uoi_feproblem

const THREAD_ID	_uoi_tid

const MooseObject &	_ppi_moose_object

const InputParameters &	_ppi_params

const FEProblemBase &	_ppi_feproblem

std::map< PostprocessorName, std::unique_ptr< PostprocessorValue > >	_default_values

const bool	_broadcast_by_default

const MooseObject &	_vpi_moose_object

const FEProblemBase &	_vpi_feproblem

const THREAD_ID	_vpi_tid

const InputParameters &	_ri_params

FEProblemBase &	_ri_fe_problem_base

const ReporterData &	_ri_reporter_data

const MooseObject &	_ri_moose_object

const InputParameters &	_dni_params

FEProblemBase &	_dni_feproblem

const MooseObject *const	_dni_moose_object_ptr

const InputParameters &	_si_params

FEProblemBase &	_si_feproblem

THREAD_ID	_si_tid

const RestartableDataMapName	_metaname

std::string	_restartable_name

MooseApp &	_meta_data_app

const MooseObject *const	_meta_data_object

const InputParameters &	_sc_parameters

const std::string &	_sc_name

const bool	_sc_is_implicit

std::unordered_map< std::string, std::vector< MooseVariableScalar * > >	_coupled_scalar_vars

std::unordered_map< std::string, std::unique_ptr< VariableValue > >	_default_value

std::unordered_map< std::string, std::unique_ptr< ADVariableValue > >	_dual_default_value

std::vector< MooseVariableScalar * >	_coupled_moose_scalar_vars

std::unordered_map< std::string, std::vector< MooseVariableFieldBase * > >	_sc_coupled_vars

std::set< TagID >	_sc_coupleable_vector_tags

std::set< TagID >	_sc_coupleable_matrix_tags

const MooseObject &	_mi_moose_object

const bool	_mi_boundary_restricted

const std::set< SubdomainID > &	_mi_block_ids

const std::set< BoundaryID > &	_mi_boundary_ids

std::vector< std::unique_ptr< OptionalMaterialPropertyProxyBase< MaterialPropertyInterface > > >	_optional_property_proxies

const std::string	_ti_name

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

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

A calorically imperfect gas fluid property class This fluid property assumes that internal energy is a general monotonic function of temperature; behaves like an ideal gas otherwise.

In particular, it assumes: p v = R_s T

Definition at line 26 of file CaloricallyImperfectGas.h.

Member Enumeration Documentation

◆ NaNMessage

enum NaNInterface::NaNMessage

protectedinherited

Enumerator
NAN_MESSAGE_NONE
NAN_MESSAGE_WARNING
NAN_MESSAGE_EXCEPTION
NAN_MESSAGE_ERROR

Definition at line 30 of file NaNInterface.h.

  {
    NAN_MESSAGE_NONE = 0,
    NAN_MESSAGE_WARNING = 1,
    NAN_MESSAGE_EXCEPTION = 2,
    NAN_MESSAGE_ERROR = 3
  };

Constructor & Destructor Documentation

◆ CaloricallyImperfectGas()

CaloricallyImperfectGas::CaloricallyImperfectGas ( const InputParameters & parameters )

Definition at line 42 of file CaloricallyImperfectGas.C.

  : SinglePhaseFluidProperties(parameters),
    NaNInterface(this),
    _molar_mass(getParam<Real>("molar_mass")),
    _R_specific(_R / _molar_mass),
    _T_c(getParam<Real>("T_c")),
    _rho_c(getParam<Real>("rho_c")),
    _e_c(getParam<Real>("e_c")),
    _min_temperature(getParam<Real>("min_temperature")),
    _max_temperature(getParam<Real>("max_temperature")),
    _delta_T(getParam<Real>("temperature_resolution")),
    _out_of_bound_error(getParam<bool>("out_of_bound_error")),
    _tol(1e-4)
{
}

Member Function Documentation

◆ c_from_p_T() [1/3]

ADReal CaloricallyImperfectGas::c_from_p_T	(	const ADReal &	p,
		const ADReal &	T
	)		const

overridevirtual

Definition at line 369 of file CaloricallyImperfectGas.C.

{
  using std::sqrt;
  return sqrt(gamma_from_p_T(p, T) * _R_specific * T);
}

◆ c_from_p_T() [2/3]

Real CaloricallyImperfectGas::c_from_p_T	(	Real	p,
		Real	T
	)		const

overridevirtual

Definition at line 363 of file CaloricallyImperfectGas.C.

{
  return std::sqrt(gamma_from_p_T(p, T) * _R_specific * T);
}

Referenced by TEST_F().

◆ c_from_p_T() [3/3]

void CaloricallyImperfectGas::c_from_p_T	(	Real	p,
		Real	T,
		Real &	c,
		Real &	dc_dp,
		Real &	dc_dT
	)		const

overridevirtual

Definition at line 376 of file CaloricallyImperfectGas.C.

{
  Real gamma, dgamma_dp, dgamma_dT;
  gamma_from_p_T(p, T, gamma, dgamma_dp, dgamma_dT);
  c = std::sqrt(gamma * _R_specific * T);
  dc_dp = 0;
  dc_dT = 0.5 / c * _R_specific * (gamma + T * dgamma_dT);
}

◆ c_from_v_e() [1/3]

ADReal CaloricallyImperfectGas::c_from_v_e	(	const ADReal &	v,
		const ADReal &	e
	)		const

overridevirtual

Definition at line 334 of file CaloricallyImperfectGas.C.

{
  using std::sqrt;
  const auto T = T_from_v_e(v, e);
 
  const auto c2 = gamma_from_v_e(v, e) * _R_specific * T;
  if (MetaPhysicL::raw_value(c2) < 0)
    return getNaN("Sound speed squared (gamma * R * T) is negative: c2 = " +
                  Moose::stringify(MetaPhysicL::raw_value(c2)) + ".");
 
  return sqrt(c2);
}

◆ c_from_v_e() [2/3]

Real CaloricallyImperfectGas::c_from_v_e	(	Real	v,
		Real	e
	)		const

overridevirtual

Definition at line 321 of file CaloricallyImperfectGas.C.

{
  Real T = T_from_v_e(v, e);
 
  const Real c2 = gamma_from_v_e(v, e) * _R_specific * T;
  if (c2 < 0)
    return getNaN("Sound speed squared (gamma * R * T) is negative: c2 = " + Moose::stringify(c2) +
                  ".");
 
  return std::sqrt(c2);
}

Referenced by TEST_F().

◆ c_from_v_e() [3/3]

void CaloricallyImperfectGas::c_from_v_e	(	Real	v,
		Real	e,
		Real &	c,
		Real &	dc_dv,
		Real &	dc_de
	)		const

overridevirtual

Definition at line 348 of file CaloricallyImperfectGas.C.

{
  Real T, dT_dv, dT_de;
  T_from_v_e(v, e, T, dT_dv, dT_de);
 
  Real gamma, dgamma_dv, dgamma_dT;
  gamma_from_v_e(v, e, gamma, dgamma_dv, dgamma_dT);
  c = std::sqrt(gamma * _R_specific * T);
 
  const Real dc_dT = 0.5 / c * _R_specific * (gamma + T * dgamma_dT);
  dc_dv = dc_dT * dT_dv;
  dc_de = dc_dT * dT_de;
}

◆ cp_from_p_T() [1/2]

Real CaloricallyImperfectGas::cp_from_p_T	(	Real	p,
		Real	T
	)		const

overridevirtual

Definition at line 868 of file CaloricallyImperfectGas.C.

{
  return cp_from_T(temperature);
}

Referenced by gamma_from_p_T(), gamma_from_p_T(), h_from_p_T(), s_from_p_T(), T_from_p_h(), and TEST_F().

◆ cp_from_p_T() [2/2]

void CaloricallyImperfectGas::cp_from_p_T	(	Real	p,
		Real	T,
		Real &	cp,
		Real &	dcp_dp,
		Real &	dcp_dT
	)		const

overridevirtual

Definition at line 874 of file CaloricallyImperfectGas.C.

{
  cp_from_T(T, cp, dcp_dT);
  dcp_dp = 0.0;
}

◆ cp_from_T() [1/2]

Real CaloricallyImperfectGas::cp_from_T ( Real T ) const

protected

Definition at line 177 of file CaloricallyImperfectGas.C.

{
  if (T < _min_temperature || T > _max_temperature)
    outOfBounds("cp_from_T", T, _min_temperature, _max_temperature);
 
  return _e_T->timeDerivative(T, Point()) + _R_specific;
}

Referenced by cp_from_p_T(), cp_from_p_T(), cp_from_T(), cp_from_v_e(), cp_from_v_e(), and h_from_T_v().

◆ cp_from_T() [2/2]

void CaloricallyImperfectGas::cp_from_T	(	Real	T,
		Real &	cp,
		Real &	dcp_dT
	)		const

protected

Definition at line 197 of file CaloricallyImperfectGas.C.

{
  if (T < _min_temperature || T > _max_temperature)
    outOfBounds("cp_from_T (3 args)", T, _min_temperature, _max_temperature);
  Real pert = 1.0e-7;
  cp = cp_from_T(T);
  Real cp_pert = cp_from_T(T * (1 + pert));
  dcp_dT = (cp_pert - cp) / (T * pert);
}

◆ cp_from_v_e() [1/2]

Real CaloricallyImperfectGas::cp_from_v_e	(	Real	v,
		Real	e
	)		const

overridevirtual

Definition at line 387 of file CaloricallyImperfectGas.C.

{
  Real T = T_from_v_e(v, e);
  return cp_from_T(T);
}

Referenced by e_from_v_h(), gamma_from_v_e(), gamma_from_v_e(), and TEST_F().

◆ cp_from_v_e() [2/2]

void CaloricallyImperfectGas::cp_from_v_e	(	Real	v,
		Real	e,
		Real &	cp,
		Real &	dcp_dv,
		Real &	dcp_de
	)		const

overridevirtual

Definition at line 394 of file CaloricallyImperfectGas.C.

{
  Real T = T_from_v_e(v, e);
  Real dcp_dT;
  cp_from_T(T, cp, dcp_dT);
  dcp_dv = 0.0;
  dcp_de = dcp_dT / cv_from_T(T);
}

◆ criticalDensity()

Real CaloricallyImperfectGas::criticalDensity ( ) const

overridevirtual

Critical density.

Returns: critical density (kg/m^3)

Reimplemented from SinglePhaseFluidProperties.

Definition at line 828 of file CaloricallyImperfectGas.C.

{
  return _rho_c;
}

◆ criticalInternalEnergy()

Real CaloricallyImperfectGas::criticalInternalEnergy ( ) const

overridevirtual

Critical specific internal energy.

Returns: specific internal energy (J/kg)

Reimplemented from SinglePhaseFluidProperties.

Definition at line 834 of file CaloricallyImperfectGas.C.

{
  return _e_c;
}

◆ criticalPressure()

Real SinglePhaseFluidProperties::criticalPressure ( ) const

virtualinherited

Critical pressure.

Returns: critical pressure (Pa)

Reimplemented in CO2FluidProperties, HydrogenFluidProperties, MethaneFluidProperties, NaClFluidProperties, NitrogenFluidProperties, TabulatedFluidProperties, Water97FluidProperties, and TestSinglePhaseFluidProperties.

Definition at line 302 of file SinglePhaseFluidProperties.C.

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

Referenced by SinglePhaseFluidProperties::criticalDensity(), SinglePhaseFluidProperties::criticalInternalEnergy(), and TabulatedFluidProperties::criticalPressure().

◆ criticalTemperature()

Real CaloricallyImperfectGas::criticalTemperature ( ) const

overridevirtual

Critical temperature.

Returns: critical temperature (K)

Reimplemented from SinglePhaseFluidProperties.

Definition at line 822 of file CaloricallyImperfectGas.C.

{
  return _T_c;
}

◆ cv_from_p_T() [1/2]

Real CaloricallyImperfectGas::cv_from_p_T	(	Real	p,
		Real	T
	)		const

overridevirtual

Definition at line 854 of file CaloricallyImperfectGas.C.

{
  return cv_from_T(temperature);
}

Referenced by e_from_p_T(), gamma_from_p_T(), gamma_from_p_T(), k_from_v_e(), mu_from_v_e(), and TEST_F().

◆ cv_from_p_T() [2/2]

void CaloricallyImperfectGas::cv_from_p_T	(	Real	p,
		Real	T,
		Real &	cv,
		Real &	dcv_dp,
		Real &	dcv_dT
	)		const

overridevirtual

Definition at line 860 of file CaloricallyImperfectGas.C.

{
  cv_from_T(T, cv, dcv_dT);
  dcv_dp = 0.0;
}

◆ cv_from_T() [1/2]

template<typename CppType >

Real CaloricallyImperfectGas::cv_from_T ( const CppType & T ) const

protected

Definition at line 218 of file CaloricallyImperfectGas.h.

{
  const auto raw_T = MetaPhysicL::raw_value(T);
  if (raw_T < _min_temperature || raw_T > _max_temperature)
    outOfBounds("cv_from_T", raw_T, _min_temperature, _max_temperature);
 
  return _e_T->timeDerivative(raw_T, Point());
}

Referenced by cp_from_v_e(), cv_from_p_T(), cv_from_p_T(), cv_from_T(), cv_from_T_v(), cv_from_v_e(), cv_from_v_e(), p_from_v_e(), and setupLookupTables().

◆ cv_from_T() [2/2]

void CaloricallyImperfectGas::cv_from_T	(	Real	T,
		Real &	cv,
		Real &	dcv_dT
	)		const

protected

Definition at line 186 of file CaloricallyImperfectGas.C.

{
  if (T < _min_temperature || T > _max_temperature)
    outOfBounds("cv_from_T (3 args)", T, _min_temperature, _max_temperature);
  Real pert = 1.0e-7;
  cv = cv_from_T(T);
  Real cv_pert = cv_from_T(T * (1 + pert));
  dcv_dT = (cv_pert - cv) / (T * pert);
}

◆ cv_from_T_v()

Real CaloricallyImperfectGas::cv_from_T_v	(	Real	T,
		Real	v
	)		const

overridevirtual

Definition at line 754 of file CaloricallyImperfectGas.C.

{
  return cv_from_T(T);
}

Referenced by e_from_T_v().

◆ cv_from_v_e() [1/3]

ADReal CaloricallyImperfectGas::cv_from_v_e	(	ADReal	v,
		ADReal	e
	)		const

virtual

Definition at line 411 of file CaloricallyImperfectGas.C.

{
  Real x = 0;
  Real raw1 = v.value();
  Real raw2 = e.value();
  Real dxd1 = 0;
  Real dxd2 = 0;
  cv_from_v_e(raw1, raw2, x, dxd1, dxd2);
 
  ADReal result = x;
  result.derivatives() = v.derivatives() * dxd1 + e.derivatives() * dxd2;
  return result;
}

◆ cv_from_v_e() [2/3]

Real CaloricallyImperfectGas::cv_from_v_e	(	Real	v,
		Real	e
	)		const

overridevirtual

Definition at line 404 of file CaloricallyImperfectGas.C.

{
  Real T = T_from_v_e(v, e);
  return cv_from_T(T);
}

Referenced by cv_from_v_e(), e_from_v_h(), gamma_from_v_e(), gamma_from_v_e(), T_from_v_e(), and TEST_F().

◆ cv_from_v_e() [3/3]

void CaloricallyImperfectGas::cv_from_v_e	(	Real	v,
		Real	e,
		Real &	cv,
		Real &	dcv_dv,
		Real &	dcv_de
	)		const

overridevirtual

Definition at line 426 of file CaloricallyImperfectGas.C.

{
  Real T = T_from_v_e(v, e);
  Real dcv_dT;
  cv_from_T(T, cv, dcv_dT);
  dcv_dv = 0.0;
  dcv_de = dcv_dT / cv;
}

◆ e_from_p_T() [1/3]

ADReal CaloricallyImperfectGas::e_from_p_T	(	ADReal	p,
		ADReal	T
	)		const

virtual

Definition at line 793 of file CaloricallyImperfectGas.C.

{
  Real x = 0;
  Real raw1 = p.value();
  Real raw2 = T.value();
  Real dxd1 = 0;
  Real dxd2 = 0;
  e_from_p_T(raw1, raw2, x, dxd1, dxd2);
 
  ADReal result = x;
  result.derivatives() = p.derivatives() * dxd1 + T.derivatives() * dxd2;
  return result;
}

◆ e_from_p_T() [2/3]

Real CaloricallyImperfectGas::e_from_p_T	(	Real	p,
		Real	T
	)		const

overridevirtual

Definition at line 787 of file CaloricallyImperfectGas.C.

{
  return e_from_T(T);
}

Referenced by e_from_p_T(), e_from_p_T(), setupLookupTables(), and TEST_F().

◆ e_from_p_T() [3/3]

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

overridevirtual

Definition at line 808 of file CaloricallyImperfectGas.C.

{
  e = e_from_p_T(p, T);
  de_dp = 0.0;
  de_dT = cv_from_p_T(p, T);
}

◆ e_from_T()

Real CaloricallyImperfectGas::e_from_T ( Real T ) const

protected

helper functions for e(T) and h(T), T(e), T(h), cv(T), cp(T)

Definition at line 160 of file CaloricallyImperfectGas.C.

{
  if (T < _min_temperature || T > _max_temperature)
    outOfBounds("e_from_T", T, _min_temperature, _max_temperature);
 
  return _e_T->value(T, Point());
}

Referenced by e_from_p_T(), e_from_T_v(), e_from_v_h(), and p_from_h_s().

◆ e_from_T_v() [1/3]

ADReal CaloricallyImperfectGas::e_from_T_v	(	const ADReal &	T,
		const ADReal &	v
	)		const

overridevirtual

Definition at line 693 of file CaloricallyImperfectGas.C.

{
  Real x = 0;
  Real raw1 = T.value();
  Real raw2 = v.value();
  Real dxd1 = 0;
  Real dxd2 = 0;
  e_from_T_v(raw1, raw2, x, dxd1, dxd2);
 
  ADReal result = x;
  result.derivatives() = T.derivatives() * dxd1 + v.derivatives() * dxd2;
  return result;
}

◆ e_from_T_v() [2/3]

Real CaloricallyImperfectGas::e_from_T_v	(	Real	T,
		Real	v
	)		const

overridevirtual

Definition at line 679 of file CaloricallyImperfectGas.C.

{
  return e_from_T(T);
}

Referenced by e_from_T_v(), e_from_T_v(), and TEST_F().

◆ e_from_T_v() [3/3]

void CaloricallyImperfectGas::e_from_T_v	(	Real	T,
		Real	v,
		Real &	e,
		Real &	de_dT,
		Real &	de_dv
	)		const

overridevirtual

Definition at line 685 of file CaloricallyImperfectGas.C.

{
  e = e_from_T_v(T, v);
  de_dT = cv_from_T_v(T, v);
  de_dv = 0.0;
}

◆ e_from_v_h() [1/2]

Real CaloricallyImperfectGas::e_from_v_h	(	Real	v,
		Real	h
	)		const

overridevirtual

Definition at line 606 of file CaloricallyImperfectGas.C.

{
  Real T = T_from_h(h);
  return e_from_T(T);
}

Referenced by e_from_v_h(), and TEST_F().

◆ e_from_v_h() [2/2]

void CaloricallyImperfectGas::e_from_v_h	(	Real	v,
		Real	h,
		Real &	e,
		Real &	de_dv,
		Real &	de_dh
	)		const

overridevirtual

Definition at line 613 of file CaloricallyImperfectGas.C.

{
  e = e_from_v_h(v, h);
  Real cv = cv_from_v_e(v, e);
  Real cp = cp_from_v_e(v, e);
  de_dv = 0.0;
  de_dh = cv / cp;
}

◆ e_spndl_from_v()

Real CaloricallyImperfectGas::e_spndl_from_v ( Real v ) const

overridevirtual

Specific internal energy from temperature and specific volume.

Parameters

[in]	T	temperature
[in]	v	specific volume

Reimplemented from SinglePhaseFluidProperties.

Definition at line 760 of file CaloricallyImperfectGas.C.

{
  return _e_c;
}

◆ 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 CaloricallyImperfectGas::fluidName ( ) const

overridevirtual

Fluid name.

Returns: string representing fluid name

Reimplemented from SinglePhaseFluidProperties.

Definition at line 154 of file CaloricallyImperfectGas.C.

{
  return "caloric_imperfect_gas";
}

◆ g_from_v_e()

Real CaloricallyImperfectGas::g_from_v_e	(	Real	v,
		Real	e
	)		const

overridevirtual

Definition at line 914 of file CaloricallyImperfectGas.C.

{
  // use the definition of the Gibbs free energy
  Real T = T_from_v_e(v, e);
  Real p = p_from_v_e(v, e);
  return h_from_p_T(p, T) - T * s_from_p_T(p, T);
}

◆ gamma_from_p_T() [1/3]

ADReal CaloricallyImperfectGas::gamma_from_p_T	(	ADReal	p,
		ADReal	T
	)		const

virtual

Definition at line 489 of file CaloricallyImperfectGas.C.

{
  Real x = 0;
  Real raw1 = p.value();
  Real raw2 = T.value();
  Real dxd1 = 0;
  Real dxd2 = 0;
  gamma_from_p_T(raw1, raw2, x, dxd1, dxd2);
  ADReal result = x;
  result.derivatives() = p.derivatives() * dxd1 + T.derivatives() * dxd2;
  return result;
}

◆ gamma_from_p_T() [2/3]

Real CaloricallyImperfectGas::gamma_from_p_T	(	Real	p,
		Real	T
	)		const

overridevirtual

Definition at line 483 of file CaloricallyImperfectGas.C.

{
  return cp_from_p_T(p, T) / cv_from_p_T(p, T);
}

Referenced by c_from_p_T(), c_from_p_T(), c_from_p_T(), gamma_from_p_T(), and TEST_F().

◆ gamma_from_p_T() [3/3]

void CaloricallyImperfectGas::gamma_from_p_T	(	Real	p,
		Real	T,
		Real &	gamma,
		Real &	dgamma_dp,
		Real &	dgamma_dT
	)		const

overridevirtual

Definition at line 470 of file CaloricallyImperfectGas.C.

{
  Real cp, dcp_dp, dcp_dT;
  cp_from_p_T(p, T, cp, dcp_dp, dcp_dT);
  Real cv, dcv_dp, dcv_dT;
  cv_from_p_T(p, T, cv, dcv_dp, dcv_dT);
  gamma = cp / cv;
  dgamma_dp = 0.0;
  dgamma_dT = 1.0 / cv * dcp_dT - gamma / cv * dcv_dT;
}

◆ gamma_from_v_e() [1/3]

ADReal CaloricallyImperfectGas::gamma_from_v_e	(	ADReal	v,
		ADReal	e
	)		const

virtual

Definition at line 442 of file CaloricallyImperfectGas.C.

{
  Real x = 0;
  Real raw1 = v.value();
  Real raw2 = e.value();
  Real dxd1 = 0;
  Real dxd2 = 0;
  gamma_from_v_e(raw1, raw2, x, dxd1, dxd2);
 
  ADReal result = x;
  result.derivatives() = v.derivatives() * dxd1 + e.derivatives() * dxd2;
  return result;
}

◆ gamma_from_v_e() [2/3]

Real CaloricallyImperfectGas::gamma_from_v_e	(	Real	v,
		Real	e
	)		const

overridevirtual

Definition at line 436 of file CaloricallyImperfectGas.C.

{
  return cp_from_v_e(v, e) / cv_from_v_e(v, e);
}

Referenced by c_from_v_e(), c_from_v_e(), c_from_v_e(), gamma_from_v_e(), and TEST_F().

◆ gamma_from_v_e() [3/3]

void CaloricallyImperfectGas::gamma_from_v_e	(	Real	v,
		Real	e,
		Real &	gamma,
		Real &	dgamma_dv,
		Real &	dgamma_de
	)		const

overridevirtual

Definition at line 457 of file CaloricallyImperfectGas.C.

{
  Real cp, dcp_dv, dcp_de;
  cp_from_v_e(v, e, cp, dcp_dv, dcp_de);
  Real cv, dcv_dv, dcv_de;
  cv_from_v_e(v, e, cv, dcv_dv, dcv_de);
  gamma = cp / cv;
  dgamma_dv = 0.0;
  dgamma_de = 1.0 / cv * dcp_de - gamma / cv * dcv_de;
}

◆ getNaN() [1/2]

Real NaNInterface::getNaN ( ) const

inlineprotectedinherited

Throws an error or returns a NaN with or without a warning, with a default message.

Definition at line 46 of file NaNInterface.h.

46{ return getNaN("A NaN was produced."); }

Referenced by c_from_v_e(), c_from_v_e(), StiffenedGasFluidProperties::c_from_v_e(), StiffenedGasFluidProperties::c_from_v_e(), ADNumericalFlux3EqnHLLC::calcFlux(), NumericalFluxGasMixHLLC::calcFlux(), NaNInterface::getNaN(), p_from_v_e(), IdealGasFluidProperties::p_from_v_e(), p_from_v_e(), IdealGasFluidProperties::p_from_v_e(), NaNInterfaceTestFluidProperties::p_from_v_e(), StiffenedGasFluidProperties::rho_from_p_s(), IdealGasFluidProperties::s_from_h_p(), StiffenedGasFluidProperties::s_from_h_p(), StiffenedGasFluidProperties::s_from_h_p(), IdealGasFluidProperties::s_from_p_T(), IdealGasFluidProperties::s_from_p_T(), IdealGasFluidProperties::s_from_v_e(), IdealGasFluidProperties::s_from_v_e(), T_from_v_e(), and IdealRealGasMixtureFluidProperties::xs_prim_from_p_T().

◆ getNaN() [2/2]

template<typename... Args>

Real NaNInterface::getNaN ( Args &&... args ) const

inlineprotectedinherited

Throws an error or returns a NaN with or without a warning.

Definition at line 62 of file NaNInterface.h.

  {
    switch (_emit_on_nan)
    {
      case (NAN_MESSAGE_WARNING):
        mooseWarning(_moose_object->name(), ": ", std::forward<Args>(args)...);
        break;
      case (NAN_MESSAGE_EXCEPTION):
        mooseException(_moose_object->name(), ": ", std::forward<Args>(args)...);
        break;
      case (NAN_MESSAGE_ERROR):
        mooseError(_moose_object->name(), ": ", std::forward<Args>(args)...);
        break;
      default:
        break;
    }
    // return a quiet NaN
    return std::nan("");
  }

◆ getNaNVector() [1/2]

std::vector< Real > NaNInterface::getNaNVector ( const unsigned int & n ) const

inlineprotectedinherited

Throws an error or returns NaNs with or without a warning, with a default message.

Parameters

[in] n Vector size

Definition at line 53 of file NaNInterface.h.

  {
    return getNaNVector(n, "A NaN was produced.");
  }

Referenced by NaNInterface::getNaNVector(), and NaNInterfaceTestFluidProperties::returnNaNVector().

◆ getNaNVector() [2/2]

template<typename... Args>

std::vector< Real > NaNInterface::getNaNVector	(	const unsigned int &	n,
		Args &&...	args
	)		const

inlineprotectedinherited

Throws an error or returns NaNs with or without a warning.

Parameters

[in] n Vector size

Definition at line 88 of file NaNInterface.h.

  {
    switch (_emit_on_nan)
    {
      case (NAN_MESSAGE_WARNING):
        mooseWarning(_moose_object->name(), ": ", std::forward<Args>(args)...);
        break;
      case (NAN_MESSAGE_EXCEPTION):
        mooseException(_moose_object->name(), ": ", std::forward<Args>(args)...);
        break;
      case (NAN_MESSAGE_ERROR):
        mooseError(_moose_object->name(), ": ", std::forward<Args>(args)...);
        break;
      default:
        break;
    }
    // return quiet NaNs
    return std::vector<Real>(n, std::nan(""));
  }

◆ h_from_p_T() [1/2]

Real CaloricallyImperfectGas::h_from_p_T	(	Real	p,
		Real	T
	)		const

overridevirtual

Definition at line 773 of file CaloricallyImperfectGas.C.

{
  return h_from_T(T);
}

Referenced by g_from_v_e(), h_from_p_T(), setupLookupTables(), and TEST_F().

◆ h_from_p_T() [2/2]

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

overridevirtual

Definition at line 779 of file CaloricallyImperfectGas.C.

{
  h = h_from_p_T(p, T);
  dh_dp = 0.0;
  dh_dT = cp_from_p_T(p, T);
}

◆ h_from_T()

Real CaloricallyImperfectGas::h_from_T ( Real T ) const

protected

Definition at line 169 of file CaloricallyImperfectGas.C.

{
  if (T < _min_temperature || T > _max_temperature)
    outOfBounds("h_from_T", T, _min_temperature, _max_temperature);
  return _e_T->value(T, Point()) + _R_specific * T;
}

Referenced by h_from_p_T(), h_from_T_v(), and h_from_T_v().

◆ h_from_T_v() [1/2]

Real CaloricallyImperfectGas::h_from_T_v	(	Real	T,
		Real	v
	)		const

overridevirtual

Definition at line 722 of file CaloricallyImperfectGas.C.

{
  return h_from_T(T);
}

◆ h_from_T_v() [2/2]

void CaloricallyImperfectGas::h_from_T_v	(	Real	T,
		Real	v,
		Real &	h,
		Real &	dh_dT,
		Real &	dh_dv
	)		const

overridevirtual

Definition at line 728 of file CaloricallyImperfectGas.C.

{
  h = h_from_T(T);
  dh_dT = cp_from_T(T);
  dh_dv = 0.0;
}

◆ henryCoefficients()

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

virtualinherited

Henry's law coefficients for dissolution in water.

Returns: Henry's constant coefficients

Reimplemented in CO2FluidProperties, HydrogenFluidProperties, MethaneFluidProperties, NitrogenFluidProperties, and TabulatedFluidProperties.

Definition at line 378 of file SinglePhaseFluidProperties.C.

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

Referenced by TabulatedFluidProperties::henryCoefficients().

◆ initialize()

virtual void FluidProperties::initialize ( )

inlinefinalvirtualinherited

Implements ThreadedGeneralUserObject.

Definition at line 34 of file FluidProperties.h.

34{}

◆ initialSetup()

void CaloricallyImperfectGas::initialSetup ( )

overridevirtual

Reimplemented from ThreadedGeneralUserObject.

Definition at line 59 of file CaloricallyImperfectGas.C.

{
  SinglePhaseFluidProperties::initialSetup();
  _e_T = &getFunction("e");
  _mu_T = &getFunction("mu");
  _k_T = &getFunction("k");
 
  setupLookupTables();
}

◆ k_from_p_T() [1/2]

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

overridevirtual

Definition at line 898 of file CaloricallyImperfectGas.C.

{
  return _k_T->value(temperature, Point());
}

Referenced by k_from_p_T(), k_from_v_e(), and k_from_v_e().

◆ k_from_p_T() [2/2]

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

overridevirtual

Definition at line 904 of file CaloricallyImperfectGas.C.

{
  k = k_from_p_T(p, T);
  dk_dp = 0.0;
  Real pert = 1.0e-7;
  Real k2 = this->k_from_p_T(p, T * (1 + pert));
  dk_dT = (k2 - k) / (T * pert);
}

◆ k_from_v_e() [1/2]

Real CaloricallyImperfectGas::k_from_v_e	(	Real	v,
		Real	e
	)		const

overridevirtual

Definition at line 524 of file CaloricallyImperfectGas.C.

{
  const Real T = T_from_v_e(v, e);
  const Real p = p_from_v_e(v, e);
  return k_from_p_T(p, T);
}

◆ k_from_v_e() [2/2]

void CaloricallyImperfectGas::k_from_v_e	(	Real	v,
		Real	e,
		Real &	k,
		Real &	dk_dv,
		Real &	dk_de
	)		const

overridevirtual

Definition at line 532 of file CaloricallyImperfectGas.C.

{
  const Real T = T_from_v_e(v, e);
  const Real p = p_from_v_e(v, e);
  Real dk_dp, dk_dT;
  k_from_p_T(p, T, k, dk_dp, dk_dT);
  // only dk_de = dk/dT * dT/de = 1/cv * dk/dT is nonzero
  // (dk/dv)_e = 0 because e only depends on T and e is constant
  dk_dv = 0.0;
  dk_de = dk_dT / cv_from_p_T(p, T);
}

◆ makeZeroAndOne() [1/2]

template<>

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

inlineinherited

Definition at line 485 of file SinglePhaseFluidProperties.h.

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

◆ makeZeroAndOne() [2/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 478 of file SinglePhaseFluidProperties.h.

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

Referenced by SinglePhaseFluidProperties::xyDerivatives().

◆ molarMass()

Real CaloricallyImperfectGas::molarMass ( ) const

overridevirtual

Molar mass [kg/mol].

Returns: molar mass

Reimplemented from SinglePhaseFluidProperties.

Definition at line 816 of file CaloricallyImperfectGas.C.

{
  return _molar_mass;
}

◆ mu_from_p_T() [1/2]

Real CaloricallyImperfectGas::mu_from_p_T	(	Real	p,
		Real	T
	)		const

overridevirtual

Definition at line 882 of file CaloricallyImperfectGas.C.

{
  return _mu_T->value(temperature, Point());
}

Referenced by mu_from_p_T(), mu_from_v_e(), and mu_from_v_e().

◆ mu_from_p_T() [2/2]

void CaloricallyImperfectGas::mu_from_p_T	(	Real	p,
		Real	T,
		Real &	mu,
		Real &	dmu_dp,
		Real &	dmu_dT
	)		const

overridevirtual

Definition at line 888 of file CaloricallyImperfectGas.C.

{
  mu = this->mu_from_p_T(p, T);
  dmu_dp = 0.0;
  Real pert = 1.0e-7;
  Real mu2 = this->mu_from_p_T(p, T * (1 + pert));
  dmu_dT = (mu2 - mu) / (T * pert);
}

◆ mu_from_v_e() [1/2]

Real CaloricallyImperfectGas::mu_from_v_e	(	Real	v,
		Real	e
	)		const

overridevirtual

Definition at line 503 of file CaloricallyImperfectGas.C.

{
  const Real T = T_from_v_e(v, e);
  const Real p = p_from_v_e(v, e);
  return mu_from_p_T(p, T);
}

◆ mu_from_v_e() [2/2]

void CaloricallyImperfectGas::mu_from_v_e	(	Real	v,
		Real	e,
		Real &	mu,
		Real &	dmu_dv,
		Real &	dmu_de
	)		const

overridevirtual

Definition at line 511 of file CaloricallyImperfectGas.C.

{
  const Real T = T_from_v_e(v, e);
  const Real p = p_from_v_e(v, e);
  Real dmu_dp, dmu_dT;
  mu_from_p_T(p, T, mu, dmu_dp, dmu_dT);
  // only dmu_de = dmu/dT * dT/de = 1/cv * dmu/dT is nonzero
  // (dmu/dv)_e = 0 because e only depends on T and e is constant
  dmu_dv = 0.0;
  dmu_de = dmu_dT / cv_from_p_T(p, T);
}

◆ outOfBounds()

void CaloricallyImperfectGas::outOfBounds	(	const std::string &	function,
		Real	value,
		Real	min,
		Real	max
	)		const

protected

function that handles exceeding parameter limits

Definition at line 944 of file CaloricallyImperfectGas.C.

{
  std::stringstream ss;
  ss << "Function " << function << " encountered argument value of " << value
     << " which is outside of the bounds of " << min << " .. " << max;
  if (_out_of_bound_error)
    mooseError(ss.str());
  else
    mooseWarning(ss.str());
}

Referenced by cp_from_T(), cp_from_T(), cv_from_T(), cv_from_T(), e_from_T(), h_from_T(), T_from_e(), T_from_h(), and Z_from_T().

◆ p_from_h_s() [1/2]

Real CaloricallyImperfectGas::p_from_h_s	(	Real	h,
		Real	s
	)		const

overridevirtual

Definition at line 923 of file CaloricallyImperfectGas.C.

{
  Real T = T_from_h(h);
  Real e = e_from_T(T);
  Real Z = Z_from_T(T);
  Real v = std::exp((s - Z) / _R_specific);
  return p_from_v_e(v, e);
}

Referenced by p_from_h_s().

◆ p_from_h_s() [2/2]

void CaloricallyImperfectGas::p_from_h_s	(	Real	h,
		Real	s,
		Real &	p,
		Real &	dp_dh,
		Real &	dp_ds
	)		const

overridevirtual

Definition at line 933 of file CaloricallyImperfectGas.C.

{
  p = p_from_h_s(h, s);
  Real pert = 1e-7;
  Real p_pert = p_from_h_s(h * (1 + pert), s);
  dp_dh = (p_pert - p) / (h * pert);
  p_pert = p_from_h_s(h, s * (1 + pert));
  dp_ds = (p_pert - p) / (s * pert);
}

◆ p_from_T_v() [1/2]

Real CaloricallyImperfectGas::p_from_T_v	(	Real	T,
		Real	v
	)		const

overridevirtual

Definition at line 708 of file CaloricallyImperfectGas.C.

{
  return _R_specific * T / v;
}

Referenced by s_from_T_v(), s_from_T_v(), and TEST_F().

◆ p_from_T_v() [2/2]

void CaloricallyImperfectGas::p_from_T_v	(	Real	T,
		Real	v,
		Real &	p,
		Real &	dp_dT,
		Real &	dp_dv
	)		const

overridevirtual

Definition at line 714 of file CaloricallyImperfectGas.C.

{
  p = _R_specific * T / v;
  dp_dT = _R_specific / v;
  dp_dv = -_R_specific * T / (v * v);
}

◆ p_from_v_e() [1/3]

ADReal CaloricallyImperfectGas::p_from_v_e	(	const ADReal &	v,
		const ADReal &	e
	)		const

overridevirtual

Definition at line 260 of file CaloricallyImperfectGas.C.

{
  if (v.value() == 0.0)
    return getNaN("Invalid value of specific volume detected (v = " + Moose::stringify(v.value()) +
                  ").");
 
  Real x = 0;
  Real raw1 = v.value();
  Real raw2 = e.value();
  Real dxd1 = 0;
  Real dxd2 = 0;
  p_from_v_e(raw1, raw2, x, dxd1, dxd2);
 
  ADReal result = x;
  result.derivatives() = v.derivatives() * dxd1 + e.derivatives() * dxd2;
  return result;
}

◆ p_from_v_e() [2/3]

Real CaloricallyImperfectGas::p_from_v_e	(	Real	v,
		Real	e
	)		const

overridevirtual

Definition at line 251 of file CaloricallyImperfectGas.C.

{
  if (v == 0.0)
    return getNaN("Invalid value of specific volume detected (v = " + Moose::stringify(v) + ").");
 
  return _R_specific * T_from_e(e) / v;
}

Referenced by g_from_v_e(), k_from_v_e(), k_from_v_e(), mu_from_v_e(), mu_from_v_e(), p_from_h_s(), p_from_v_e(), p_from_v_e(), and TEST_F().

◆ p_from_v_e() [3/3]

void CaloricallyImperfectGas::p_from_v_e	(	Real	v,
		Real	e,
		Real &	p,
		Real &	dp_dv,
		Real &	dp_de
	)		const

overridevirtual

Definition at line 279 of file CaloricallyImperfectGas.C.

{
  p = p_from_v_e(v, e);
  Real T = T_from_e(e);
  dp_dv = -_R_specific * T / v / v;
  dp_de = _R_specific / v / cv_from_T(T);
}

◆ 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 599 of file SinglePhaseFluidProperties.h.

{
  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,
                                        "p_T_from_h_s",
                                        _max_newton_its,
                                        _verbose_newton);
    conversion_succeeded = true;
  }
  catch (MooseException &)
  {
    conversion_succeeded = false;
  }
 
  if (!conversion_succeeded)
    mooseDoOnce(mooseWarning("Conversion from (h, s)=(", h, ", ", s, ") to (p, T) failed"));
}

Referenced by TabulatedFluidProperties::T_from_h_s().

◆ 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 513 of file SinglePhaseFluidProperties.h.

{
  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,
                                        "p_T_from_v_e",
                                        _max_newton_its,
                                        _verbose_newton);
    conversion_succeeded = true;
  }
  catch (MooseException &)
  {
    conversion_succeeded = false;
  }
 
  if (!conversion_succeeded)
    mooseDoOnce(mooseWarning("Conversion from (v, e)=(", v, ", ", e, ") to (p, T) failed"));
}

Referenced by TabulatedBicubicFluidProperties::constructInterpolation(), TemperaturePressureFunctionFluidProperties::cp_from_v_e(), TemperaturePressureFunctionFluidProperties::cv_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().

◆ 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 559 of file SinglePhaseFluidProperties.h.

{
  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,
                                        "p_T_from_v_h",
                                        _max_newton_its,
                                        _verbose_newton);
    conversion_succeeded = true;
  }
  catch (MooseException &)
  {
    conversion_succeeded = false;
  }
 
  if (!conversion_succeeded)
    mooseDoOnce(mooseWarning("Conversion from (v, h)=(", v, ", ", h, ") to (p, T) failed"));
}

Referenced by TabulatedBicubicFluidProperties::constructInterpolation().

◆ propfunc() [1/12]

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

inherited

◆ propfunc() [2/12]

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

inherited

◆ propfunc() [3/12]

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

inherited

◆ propfunc() [4/12]

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

inherited

◆ propfunc() [5/12]

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

inherited

◆ propfunc() [6/12]

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

inherited

◆ propfunc() [7/12]

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

inherited

◆ propfunc() [8/12]

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

inherited

◆ propfunc() [9/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:

@begincode 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:

@begincode // 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:

@begincode 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() [10/12]

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

inherited

◆ propfunc() [11/12]

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

inherited

◆ propfunc() [12/12]

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

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 T rho v v T e p T T SinglePhaseFluidProperties::propfuncWithDefault	(	gamma	,
		v	,
		e
	)

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 SinglePhaseFluidProperties::propfuncWithDefault	(	mu	,
		p	,
		T
	)

inherited

◆ propfuncWithDefault() [5/7]

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

inherited

◆ propfuncWithDefault() [6/7]

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

inherited

◆ propfuncWithDefault() [7/7]

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

inherited

◆ propfuncWithDefinitionOverride()

CaloricallyImperfectGas::propfuncWithDefinitionOverride	(	e	,
		p	,
		rho
	)

◆ 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 438 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_s()

Real CaloricallyImperfectGas::rho_from_p_s	(	Real	p,
		Real	s
	)		const

overridevirtual

Definition at line 595 of file CaloricallyImperfectGas.C.

{
  auto s_diff = [&p, &s, this](Real x) { return this->s_from_p_T(p, x) - s; };
  Real min = _min_temperature;
  Real max = _max_temperature;
  BrentsMethod::bracket(s_diff, min, max);
  Real T = BrentsMethod::root(s_diff, min, max);
  return rho_from_p_T(p, T);
}

◆ rho_from_p_T() [1/4]

ADReal CaloricallyImperfectGas::rho_from_p_T	(	const ADReal &	p,
		const ADReal &	T
	)		const

overridevirtual

Definition at line 629 of file CaloricallyImperfectGas.C.

{
  return p * _molar_mass / (_R * T);
}

◆ rho_from_p_T() [2/4]

void CaloricallyImperfectGas::rho_from_p_T	(	const ADReal &	p,
		const ADReal &	T,
		ADReal &	rho,
		ADReal &	drho_dp,
		ADReal &	drho_dT
	)		const

overridevirtual

Definition at line 635 of file CaloricallyImperfectGas.C.

{
  rho = rho_from_p_T(p, T);
  drho_dp = _molar_mass / (_R * T);
  drho_dT = -p * _molar_mass / (_R * T * T);
}

◆ rho_from_p_T() [3/4]

Real CaloricallyImperfectGas::rho_from_p_T	(	Real	p,
		Real	T
	)		const

overridevirtual

Definition at line 623 of file CaloricallyImperfectGas.C.

{
  return p * _molar_mass / (_R * T);
}

Referenced by rho_from_p_s(), rho_from_p_T(), rho_from_p_T(), s_from_h_p(), s_from_p_T(), and TEST_F().

◆ rho_from_p_T() [4/4]

void CaloricallyImperfectGas::rho_from_p_T	(	Real	p,
		Real	T,
		Real &	rho,
		Real &	drho_dp,
		Real &	drho_dT
	)		const

overridevirtual

Definition at line 644 of file CaloricallyImperfectGas.C.

{
  rho = rho_from_p_T(p, T);
  drho_dp = _molar_mass / (_R * T);
  drho_dT = -p * _molar_mass / (_R * T * T);
}

◆ rho_mu_from_p_T() [1/3]

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

virtualinherited

Definition at line 473 of file SinglePhaseFluidProperties.C.

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

◆ rho_mu_from_p_T() [2/3]

void SinglePhaseFluidProperties::rho_mu_from_p_T	(	Real	p,
		Real	T,
		Real &	rho,
		Real &	drho_dp,
		Real &	drho_dT,
		Real &	mu,
		Real &	dmu_dp,
		Real &	dmu_dT
	)		const

virtualinherited

Reimplemented in CO2FluidProperties, HydrogenFluidProperties, NitrogenFluidProperties, and Water97FluidProperties.

Definition at line 459 of file SinglePhaseFluidProperties.C.

{
  rho_from_p_T(p, T, rho, drho_dp, drho_dT);
  mu_from_p_T(p, T, mu, dmu_dp, dmu_dT);
}

◆ rho_mu_from_p_T() [3/3]

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

virtualinherited

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 in CO2FluidProperties, HydrogenFluidProperties, NitrogenFluidProperties, and Water97FluidProperties.

Definition at line 452 of file SinglePhaseFluidProperties.C.

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

Referenced by PorousFlowBrineCO2::gasProperties(), PorousFlowWaterNCG::gasProperties(), and PorousFlowWaterNCG::liquidProperties().

◆ s_from_h_p() [1/2]

Real CaloricallyImperfectGas::s_from_h_p	(	Real	h,
		Real	p
	)		const

overridevirtual

Definition at line 578 of file CaloricallyImperfectGas.C.

{
  Real T = T_from_p_h(p, h);
  Real v = 1.0 / rho_from_p_T(p, T);
  Real Z = Z_from_T(T);
  return Z + _R_specific * std::log(v);
}

Referenced by s_from_h_p().

◆ s_from_h_p() [2/2]

void CaloricallyImperfectGas::s_from_h_p	(	Real	h,
		Real	p,
		Real &	s,
		Real &	ds_dh,
		Real &	ds_dp
	)		const

overridevirtual

Definition at line 587 of file CaloricallyImperfectGas.C.

{
  s = s_from_h_p(h, p);
  ds_dh = 1.0 / T_from_p_h(p, h);
  ds_dp = -_R_specific / p;
}

◆ s_from_p_T() [1/2]

Real CaloricallyImperfectGas::s_from_p_T	(	Real	p,
		Real	T
	)		const

overridevirtual

Definition at line 562 of file CaloricallyImperfectGas.C.

{
  Real Z = Z_from_T(T);
  Real v = 1.0 / rho_from_p_T(p, T);
  return Z + _R_specific * std::log(v);
}

Referenced by g_from_v_e(), rho_from_p_s(), s_from_p_T(), s_from_T_v(), and s_from_T_v().

◆ s_from_p_T() [2/2]

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

overridevirtual

Definition at line 570 of file CaloricallyImperfectGas.C.

{
  s = s_from_p_T(p, T);
  ds_dp = -_R_specific / p;
  ds_dT = cp_from_p_T(p, T) / T;
}

◆ s_from_T_v() [1/2]

Real CaloricallyImperfectGas::s_from_T_v	(	Real	T,
		Real	v
	)		const

overridevirtual

Definition at line 736 of file CaloricallyImperfectGas.C.

{
  Real p = p_from_T_v(T, v);
  return s_from_p_T(p, T);
}

◆ s_from_T_v() [2/2]

void CaloricallyImperfectGas::s_from_T_v	(	Real	T,
		Real	v,
		Real &	s,
		Real &	ds_dT,
		Real &	ds_dv
	)		const

overridevirtual

Definition at line 743 of file CaloricallyImperfectGas.C.

{
  Real p, dp_dT_v, dp_dv_T;
  Real ds_dp_T, ds_dT_p;
  p_from_T_v(T, v, p, dp_dT_v, dp_dv_T);
  s_from_p_T(p, T, s, ds_dp_T, ds_dT_p);
  ds_dT = ds_dT_p + ds_dp_T * dp_dT_v;
  ds_dv = ds_dp_T * dp_dv_T;
}

◆ s_from_v_e() [1/2]

Real CaloricallyImperfectGas::s_from_v_e	(	Real	v,
		Real	e
	)		const

overridevirtual

Definition at line 545 of file CaloricallyImperfectGas.C.

{
  Real T = T_from_v_e(v, e);
  Real Z = Z_from_T(T);
  return Z + _R_specific * std::log(v);
}

Referenced by s_from_v_e().

◆ s_from_v_e() [2/2]

void CaloricallyImperfectGas::s_from_v_e	(	Real	v,
		Real	e,
		Real &	s,
		Real &	ds_dv,
		Real &	ds_de
	)		const

overridevirtual

Definition at line 553 of file CaloricallyImperfectGas.C.

{
  s = s_from_v_e(v, e);
  // see documentation for derivation
  ds_dv = _R_specific / v;
  ds_de = 1.0 / T_from_v_e(v, e);
}

◆ setupLookupTables()

void CaloricallyImperfectGas::setupLookupTables ( )

protected

sets up the T(e) reverse lookup table

Definition at line 70 of file CaloricallyImperfectGas.C.

{
  // estimate number of points in inverse lookup & then adjust
  // the _delta_T so (n - 1) * _delta_T = _max_temperature - _min_temperature
  unsigned int n = std::floor((_max_temperature - _min_temperature) / _delta_T) + 1;
  _delta_T = (_max_temperature - _min_temperature) / ((Real)n - 1.0);
 
  // ensure e(T) is monotonic
  for (unsigned int j = 0; j < n; ++j)
  {
    Real temperature = _min_temperature + j * _delta_T;
 
    // Note that the function behavior at the end points
    // can lead to de/dT = 0 & we want to allow that
    Real cv = _e_T->timeDerivative(temperature, Point());
    if (cv < 0 || (cv == 0 && j > 0 && j < n - 1))
      mooseError("e(T) is not monotonically increasing with T");
  }
 
  // backward lookup tables
  _min_e = _e_T->value(_min_temperature, Point());
  _max_e = _e_T->value(_max_temperature, Point());
  _delta_e = (_max_e - _min_e) / ((Real)n - 1.0);
  _min_h = _min_e + _R_specific * _min_temperature;
  _max_h = _max_e + _R_specific * _max_temperature;
  _delta_h = (_max_h - _min_h) / ((Real)n - 1.0);
 
  _T_e_lookup.resize(n);
  _T_h_lookup.resize(n);
  _Z_T_lookup.resize(n);
  for (unsigned int j = 0; j < n; ++j)
  {
    // internal energy
    {
      Real internal_energy = _min_e + j * _delta_e;
 
      // guarding against roundoff when summing
      // at least once I saw roundoff cause an error here
      internal_energy = std::min(internal_energy, _max_e);
 
      // The temperature is found by e - e(T) = 0
      Real min = _min_temperature;
      Real max = _max_temperature;
      auto e_diff = [&internal_energy, this](Real x)
      { return this->e_from_p_T(0.0, x) - internal_energy; };
      BrentsMethod::bracket(e_diff, min, max);
      Real temperature = BrentsMethod::root(e_diff, min, max);
      _T_e_lookup[j] = temperature;
    }
 
    // enthalpy
    {
      Real enthalpy = _min_h + j * _delta_h;
 
      // guarding against roundoff when summing
      // at least once I saw roundoff cause an error here
      enthalpy = std::min(enthalpy, _max_h);
 
      // The temperature is found by h - h(T) = 0
      Real min = _min_temperature;
      Real max = _max_temperature;
      auto h_diff = [&enthalpy, this](Real x) { return this->h_from_p_T(0.0, x) - enthalpy; };
      BrentsMethod::bracket(h_diff, min, max);
      Real temperature = BrentsMethod::root(h_diff, min, max);
      _T_h_lookup[j] = temperature;
    }
 
    // Z(T)
    {
      if (j == 0)
        _Z_T_lookup[j] = 0;
      else
      {
        Real temperature = _min_temperature + j * _delta_T;
        Real temperature_prev = _min_temperature + (j - 1) * _delta_T;
        Real f1 = cv_from_T(temperature) / temperature;
        Real f2 = cv_from_T(temperature_prev) / temperature_prev;
        _Z_T_lookup[j] = _Z_T_lookup[j - 1] + 0.5 * _delta_T * (f1 + f2);
      }
    }
  }
}

Referenced by initialSetup().

◆ subdomainSetup()

virtual void FluidProperties::subdomainSetup ( )

inlinefinalvirtualinherited

Reimplemented from ThreadedGeneralUserObject.

Definition at line 38 of file FluidProperties.h.

38{}

◆ T_from_e()

Real CaloricallyImperfectGas::T_from_e ( Real e ) const

protected

Definition at line 208 of file CaloricallyImperfectGas.C.

{
  if (e < _min_e || e > _max_e)
    outOfBounds("T_from_e", e, _min_e, _max_e);
 
  if (e < _min_e)
    return _T_e_lookup[0];
  else if (e > _max_e)
    return _T_e_lookup.back();
 
  unsigned int index = std::floor((e - _min_e) / _delta_e);
  Real x = (e - _min_e - index * _delta_e) / _delta_e;
  return x * _T_e_lookup[index + 1] + (1 - x) * _T_e_lookup[index];
}

Referenced by p_from_v_e(), p_from_v_e(), T_from_v_e(), and T_from_v_e().

◆ T_from_h()

Real CaloricallyImperfectGas::T_from_h ( Real h ) const

protected

Definition at line 224 of file CaloricallyImperfectGas.C.

{
  if (h < _min_h || h > _max_h)
    outOfBounds("h_from_e", h, _min_h, _max_h);
 
  if (h < _min_h)
    return _T_h_lookup[0];
  else if (h > _max_h)
    return _T_h_lookup.back();
 
  unsigned int index = std::floor((h - _min_h) / _delta_h);
  Real x = (h - _min_h - index * _delta_h) / _delta_h;
  return x * _T_h_lookup[index + 1] + (1 - x) * _T_h_lookup[index];
}

Referenced by e_from_v_h(), p_from_h_s(), and T_from_p_h().

◆ T_from_p_h() [1/2]

Real CaloricallyImperfectGas::T_from_p_h	(	Real	p,
		Real	h
	)		const

overridevirtual

Definition at line 840 of file CaloricallyImperfectGas.C.

{
  return T_from_h(h);
}

Referenced by s_from_h_p(), s_from_h_p(), and T_from_p_h().

◆ T_from_p_h() [2/2]

void CaloricallyImperfectGas::T_from_p_h	(	Real	p,
		Real	h,
		Real &	T,
		Real &	dT_dp,
		Real &	dT_dh
	)		const

overridevirtual

Definition at line 846 of file CaloricallyImperfectGas.C.

{
  T = T_from_p_h(p, h);
  dT_dp = 0;
  dT_dh = 1.0 / cp_from_p_T(p, T);
}

◆ T_from_v_e() [1/3]

ADReal CaloricallyImperfectGas::T_from_v_e	(	const ADReal &	v,
		const ADReal &	e
	)		const

overridevirtual

Definition at line 294 of file CaloricallyImperfectGas.C.

{
  if (v.value() == 0.0)
    return getNaN("Invalid value of specific volume detected (v = " + Moose::stringify(v.value()) +
                  ").");
 
  Real x = 0;
  Real raw1 = v.value();
  Real raw2 = e.value();
  Real dxd1 = 0;
  Real dxd2 = 0;
  T_from_v_e(raw1, raw2, x, dxd1, dxd2);
 
  ADReal result = x;
  result.derivatives() = v.derivatives() * dxd1 + e.derivatives() * dxd2;
  return result;
}

◆ T_from_v_e() [2/3]

Real CaloricallyImperfectGas::T_from_v_e	(	Real	v,
		Real	e
	)		const

overridevirtual

Definition at line 288 of file CaloricallyImperfectGas.C.

{
  return T_from_e(e);
}

Referenced by c_from_v_e(), c_from_v_e(), c_from_v_e(), cp_from_v_e(), cp_from_v_e(), cv_from_v_e(), cv_from_v_e(), g_from_v_e(), k_from_v_e(), k_from_v_e(), mu_from_v_e(), mu_from_v_e(), s_from_v_e(), s_from_v_e(), T_from_v_e(), and TEST_F().

◆ T_from_v_e() [3/3]

void CaloricallyImperfectGas::T_from_v_e	(	Real	v,
		Real	e,
		Real &	T,
		Real &	dT_dv,
		Real &	dT_de
	)		const

overridevirtual

Definition at line 313 of file CaloricallyImperfectGas.C.

{
  T = T_from_e(e);
  dT_dv = 0.0;
  dT_de = 1.0 / cv_from_v_e(v, e);
}

◆ threadJoin()

virtual void FluidProperties::threadJoin ( const UserObject & )

inlinefinalvirtualinherited

Reimplemented from ThreadedGeneralUserObject.

Definition at line 37 of file FluidProperties.h.

37{}

◆ triplePointPressure()

Real SinglePhaseFluidProperties::triplePointPressure ( ) const

virtualinherited

Triple point pressure.

Returns: triple point pressure (Pa)

Reimplemented in CO2FluidProperties, HydrogenFluidProperties, MethaneFluidProperties, NaClFluidProperties, NitrogenFluidProperties, TabulatedFluidProperties, and Water97FluidProperties.

Definition at line 326 of file SinglePhaseFluidProperties.C.

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

Referenced by TabulatedFluidProperties::triplePointPressure().

◆ triplePointTemperature()

Real SinglePhaseFluidProperties::triplePointTemperature ( ) const

virtualinherited

Triple point temperature.

Returns: triple point temperature (K)

Reimplemented in CO2FluidProperties, HydrogenFluidProperties, MethaneFluidProperties, NaClFluidProperties, NitrogenFluidProperties, TabulatedFluidProperties, and Water97FluidProperties.

Definition at line 332 of file SinglePhaseFluidProperties.C.

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

Referenced by TabulatedFluidProperties::triplePointTemperature().

◆ unimplementedDerivativeMethod()

void SinglePhaseFluidProperties::unimplementedDerivativeMethod ( const std::string & property_function_name ) const

inlineprivateinherited

Definition at line 454 of file SinglePhaseFluidProperties.h.

  {
    const std::string message =
        "The fluid properties class '" + type() +
        "' has not implemented the method below, which computes derivatives of fluid properties "
        "with regards to the flow variables. If your application requires this "
        "method, you must either implement it or use a different fluid properties "
        " class.\n\n" +
        property_function_name;
 
    if (_allow_imperfect_jacobians)
      mooseDoOnce(mooseWarning(message + "\nThe unimplemented derivatives for this fluid property "
                                         "are currently neglected, set to 0."));
    else
      mooseError(message + "\n\nYou can avoid this error by neglecting the "
                           "unimplemented derivatives of fluid properties by setting the "
                           "'allow_imperfect_jacobians' parameter");
  }

Referenced by SinglePhaseFluidProperties::vaporPressure(), and SinglePhaseFluidProperties::vaporTemperature().

◆ 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 &	dv_dp,
		CppType &	dv_dT,
		CppType &	e,
		CppType &	de_dp,
		CppType &	de_dT
	)		const

inherited

Definition at line 659 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_from_p_T() [2/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 639 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);
  }
}

Referenced by TabulatedFluidProperties::h_from_p_T(), TabulatedFluidProperties::h_from_p_T(), TabulatedFluidProperties::h_from_p_T(), and TabulatedFluidProperties::s_from_p_T().

◆ v_e_spndl_from_T()

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

overridevirtual

Specific internal energy from temperature and specific volume.

Parameters

[in]	T	temperature
[in]	v	specific volume

Reimplemented from SinglePhaseFluidProperties.

Definition at line 766 of file CaloricallyImperfectGas.C.

{
  v = 1. / _rho_c;
  e = _e_c;
}

◆ validParams()

InputParameters CaloricallyImperfectGas::validParams ( )

static

Definition at line 19 of file CaloricallyImperfectGas.C.

{
  InputParameters params = SinglePhaseFluidProperties::validParams();
  params += NaNInterface::validParams();
 
  params.addRequiredParam<Real>("molar_mass", "Constant molar mass of the fluid [kg/mol]");
  params.addRequiredParam<FunctionName>("mu", "Dynamic viscosity, [Pa-s]");
  params.addRequiredParam<FunctionName>("k", "Thermal conductivity, [W/(m-K)]");
  params.addParam<Real>("T_c", 0, "Critical temperature, [K]");
  params.addParam<Real>("rho_c", 0, "Critical density, [kg/m3]");
  params.addParam<Real>("e_c", 0, "Specific internal energy at the critical point, [J/kg]");
  params.addRequiredParam<FunctionName>(
      "e", "Specific internal energy as a function of temperature [J/kg]");
  params.addRequiredParam<Real>("min_temperature", "Smallest temperature for lookup tables [K]");
  params.addRequiredParam<Real>("max_temperature", "Largest temperature for lookup tables [K]");
  params.addParam<Real>(
      "temperature_resolution", 1, "Step size in temperature for creating the inverse lookup T(e)");
  params.addParam<bool>("out_of_bound_error", true, "Error if out of bounds");
  params.addClassDescription("Fluid properties for an ideal gas with imperfect caloric behavior.");
 
  return params;
}

◆ vaporPressure() [1/3]

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

virtualinherited

Definition at line 393 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;
}

◆ vaporPressure() [2/3]

Real SinglePhaseFluidProperties::vaporPressure ( Real T ) const

virtualinherited

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 in CO2FluidProperties, HydrogenFluidProperties, MethaneFluidProperties, NitrogenFluidProperties, TabulatedFluidProperties, and Water97FluidProperties.

Definition at line 372 of file SinglePhaseFluidProperties.C.

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

Referenced by PorousFlowWaterNCG::equilibriumMassFractions(), PorousFlowWaterNCG::gasDensity(), PorousFlowWaterNCG::gasProperties(), SinglePhaseFluidProperties::vaporPressure(), SinglePhaseFluidProperties::vaporPressure(), TabulatedFluidProperties::vaporPressure(), TabulatedFluidProperties::vaporPressure(), and BrineFluidProperties::vaporPressure().

◆ vaporPressure() [3/3]

void SinglePhaseFluidProperties::vaporPressure	(	Real	T,
		Real &	psat,
		Real &	dpsat_dT
	)		const

virtualinherited

Reimplemented in CO2FluidProperties, HydrogenFluidProperties, MethaneFluidProperties, NitrogenFluidProperties, TabulatedFluidProperties, and Water97FluidProperties.

Definition at line 384 of file SinglePhaseFluidProperties.C.

{
  unimplementedDerivativeMethod(__PRETTY_FUNCTION__);
 
  dp_dT = 0.0;
  p = vaporPressure(T);
}

◆ vaporTemperature() [1/3]

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

virtualinherited

Definition at line 423 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;
}

◆ vaporTemperature() [2/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 TabulatedFluidProperties, and Water97FluidProperties.

Definition at line 408 of file SinglePhaseFluidProperties.C.

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

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

◆ vaporTemperature() [3/3]

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

virtualinherited

Reimplemented in TabulatedFluidProperties, and Water97FluidProperties.

Definition at line 414 of file SinglePhaseFluidProperties.C.

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

◆ 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 492 of file SinglePhaseFluidProperties.h.

{
  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];
}

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().

◆ Z_from_T()

Real CaloricallyImperfectGas::Z_from_T ( Real T ) const

protected

Definition at line 240 of file CaloricallyImperfectGas.C.

{
  if (T < _min_temperature || T > _max_temperature)
    outOfBounds("Z_from_T", T, _min_temperature, _max_temperature);
 
  unsigned int index = std::floor((T - _min_temperature) / _delta_T);
  Real x = (T - _min_temperature - index * _delta_T) / _delta_T;
  return x * _Z_T_lookup[index + 1] + (1 - x) * _Z_T_lookup[index];
}

Referenced by p_from_h_s(), s_from_h_p(), s_from_p_T(), and s_from_v_e().

Member Data Documentation

◆ _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(), Water97FluidProperties::T_from_p_h_ad(), and SinglePhaseFluidProperties::unimplementedDerivativeMethod().

◆ _delta_e

Real CaloricallyImperfectGas::_delta_e

protected

step size in internal energy and enthalpy

Definition at line 181 of file CaloricallyImperfectGas.h.

Referenced by setupLookupTables(), and T_from_e().

◆ _delta_h

Real CaloricallyImperfectGas::_delta_h

protected

Definition at line 182 of file CaloricallyImperfectGas.h.

Referenced by setupLookupTables(), and T_from_h().

◆ _delta_T

Real CaloricallyImperfectGas::_delta_T

protected

temperature interval in lookup tables

Definition at line 158 of file CaloricallyImperfectGas.h.

Referenced by setupLookupTables(), and Z_from_T().

◆ _e_c

Real CaloricallyImperfectGas::_e_c

protected

Definition at line 150 of file CaloricallyImperfectGas.h.

Referenced by criticalInternalEnergy(), e_spndl_from_v(), and v_e_spndl_from_T().

◆ _e_T

const Function* CaloricallyImperfectGas::_e_T

protected

Internal energy as a function of temperature.

Definition at line 167 of file CaloricallyImperfectGas.h.

Referenced by cp_from_T(), cv_from_T(), e_from_T(), h_from_T(), initialSetup(), and setupLookupTables().

◆ _emit_on_nan

enum NaNMessage NaNInterface::_emit_on_nan

protectedinherited

Raise mooseWarning or mooseError?

Definition at line 41 of file NaNInterface.h.

Referenced by NaNInterface::getNaN(), and NaNInterface::getNaNVector().

◆ _k_T

const Function* CaloricallyImperfectGas::_k_T

protected

Thermal conductivity.

Definition at line 171 of file CaloricallyImperfectGas.h.

Referenced by initialSetup(), and k_from_p_T().

◆ _max_e

Real CaloricallyImperfectGas::_max_e

protected

Definition at line 175 of file CaloricallyImperfectGas.h.

Referenced by setupLookupTables(), and T_from_e().

◆ _max_h

Real CaloricallyImperfectGas::_max_h

protected

Definition at line 177 of file CaloricallyImperfectGas.h.

Referenced by setupLookupTables(), and T_from_h().

◆ _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.

◆ _max_temperature

Real CaloricallyImperfectGas::_max_temperature

protected

Definition at line 154 of file CaloricallyImperfectGas.h.

Referenced by cp_from_T(), cp_from_T(), cv_from_T(), cv_from_T(), e_from_T(), h_from_T(), rho_from_p_s(), setupLookupTables(), and Z_from_T().

◆ _min_e

Real CaloricallyImperfectGas::_min_e

protected

internal energy and enthalpy limits when creating lookup tables

Definition at line 174 of file CaloricallyImperfectGas.h.

Referenced by setupLookupTables(), and T_from_e().

◆ _min_h

Real CaloricallyImperfectGas::_min_h

protected

Definition at line 176 of file CaloricallyImperfectGas.h.

Referenced by setupLookupTables(), and T_from_h().

◆ _min_temperature

Real CaloricallyImperfectGas::_min_temperature

protected

temperature limits when creating lookup tables

Definition at line 153 of file CaloricallyImperfectGas.h.

Referenced by cp_from_T(), cp_from_T(), cv_from_T(), cv_from_T(), e_from_T(), h_from_T(), rho_from_p_s(), setupLookupTables(), and Z_from_T().

◆ _molar_mass

const Real& CaloricallyImperfectGas::_molar_mass

protected

molar mass

Definition at line 142 of file CaloricallyImperfectGas.h.

Referenced by molarMass(), rho_from_p_T(), rho_from_p_T(), rho_from_p_T(), and rho_from_p_T().

◆ _moose_object

const MooseObject* const NaNInterface::_moose_object

protectedinherited

Definition at line 38 of file NaNInterface.h.

Referenced by NaNInterface::getNaN(), and NaNInterface::getNaNVector().

◆ _mu_T

const Function* CaloricallyImperfectGas::_mu_T

protected

Dynamic viscosity.

Definition at line 169 of file CaloricallyImperfectGas.h.

Referenced by initialSetup(), and mu_from_p_T().

◆ _out_of_bound_error

const bool CaloricallyImperfectGas::_out_of_bound_error

protected

Flag to error if out of bounds.

Definition at line 161 of file CaloricallyImperfectGas.h.

Referenced by outOfBounds().

◆ _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(), TemperaturePressureFunctionFluidProperties::cv_from_v_e(), TabulatedFluidProperties::e_from_v_h(), 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().

◆ _R

const Real FluidProperties::_R = 8.3144598

staticinherited

Universal gas constant (J/mol/K)

Definition at line 41 of file FluidProperties.h.

◆ _R_specific

const Real CaloricallyImperfectGas::_R_specific

protected

Specific gas constant (R / molar mass)

Definition at line 145 of file CaloricallyImperfectGas.h.

Referenced by c_from_p_T(), c_from_p_T(), c_from_p_T(), c_from_v_e(), c_from_v_e(), c_from_v_e(), cp_from_T(), h_from_T(), p_from_h_s(), p_from_T_v(), p_from_T_v(), p_from_v_e(), p_from_v_e(), s_from_h_p(), s_from_h_p(), s_from_p_T(), s_from_p_T(), s_from_v_e(), s_from_v_e(), and setupLookupTables().

◆ _rho_c

Real CaloricallyImperfectGas::_rho_c

protected

Definition at line 149 of file CaloricallyImperfectGas.h.

Referenced by criticalDensity(), and v_e_spndl_from_T().

◆ _T_c

Real CaloricallyImperfectGas::_T_c

protected

Definition at line 148 of file CaloricallyImperfectGas.h.

Referenced by criticalTemperature().

◆ _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 NaClFluidProperties::cp_from_p_T(), NaKFluidProperties::cp_from_p_T(), NaClFluidProperties::cp_from_p_T(), NaKFluidProperties::cp_from_p_T(), BrineFluidProperties::cp_from_p_T_X(), NaClFluidProperties::h_from_p_T(), NaKFluidProperties::h_from_p_T(), NaClFluidProperties::h_from_p_T(), NaKFluidProperties::h_from_p_T(), BrineFluidProperties::h_from_p_T_X(), BrineFluidProperties::haliteSolubility(), NaClFluidProperties::k_from_p_T(), NaKFluidProperties::k_from_p_T(), NaClFluidProperties::k_from_p_T(), NaKFluidProperties::k_from_p_T(), BrineFluidProperties::k_from_p_T_X(), NaKFluidProperties::mu_from_p_T(), NaKFluidProperties::mu_from_p_T(), BrineFluidProperties::mu_from_p_T_X(), BrineFluidProperties::mu_from_p_T_X(), CO2FluidProperties::partialDensity(), NaClFluidProperties::rho_from_p_T(), NaKFluidProperties::rho_from_p_T(), NaClFluidProperties::rho_from_p_T(), NaKFluidProperties::rho_from_p_T(), BrineFluidProperties::rho_from_p_T_X(), and NaKFluidProperties::T_from_p_h().

◆ _T_e_lookup

std::vector<Real> CaloricallyImperfectGas::_T_e_lookup

protected

Definition at line 187 of file CaloricallyImperfectGas.h.

Referenced by setupLookupTables(), and T_from_e().

◆ _T_h_lookup

std::vector<Real> CaloricallyImperfectGas::_T_h_lookup

protected

inverse lookup table data

Definition at line 186 of file CaloricallyImperfectGas.h.

Referenced by setupLookupTables(), and T_from_h().

◆ _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.

◆ _tol

Real CaloricallyImperfectGas::_tol

protected

Definition at line 164 of file CaloricallyImperfectGas.h.

◆ _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.

◆ _verbose_newton

const bool SinglePhaseFluidProperties::_verbose_newton

protectedinherited

Whether to output information about newton solves to console.

Definition at line 451 of file SinglePhaseFluidProperties.h.

◆ _Z_T_lookup

std::vector<Real> CaloricallyImperfectGas::_Z_T_lookup

protected

Z(T) lookup table on uniform grid between _min_temperature and _max_temperature.

Definition at line 191 of file CaloricallyImperfectGas.h.

Referenced by setupLookupTables(), and Z_from_T().

◆ 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/CaloricallyImperfectGas.h
fluid_properties/src/fluidproperties/CaloricallyImperfectGas.C

Public Types

Public Member Functions

Static Public Member Functions

Public Attributes

Static Public Attributes

Protected Types

Protected Member Functions

Static Protected Member Functions

Protected Attributes

Static Protected Attributes

Private Member Functions

Static Private Member Functions

Private Attributes

Detailed Description

Member Enumeration Documentation

◆ NaNMessage

Constructor & Destructor Documentation

◆ CaloricallyImperfectGas()

Member Function Documentation

◆ c_from_p_T() [1/3]

◆ c_from_p_T() [2/3]

◆ c_from_p_T() [3/3]

◆ c_from_v_e() [1/3]

◆ c_from_v_e() [2/3]

◆ c_from_v_e() [3/3]

◆ cp_from_p_T() [1/2]

◆ cp_from_p_T() [2/2]

◆ cp_from_T() [1/2]

◆ cp_from_T() [2/2]

◆ cp_from_v_e() [1/2]

◆ cp_from_v_e() [2/2]

◆ criticalDensity()

◆ criticalInternalEnergy()

◆ criticalPressure()

◆ criticalTemperature()

◆ cv_from_p_T() [1/2]

◆ cv_from_p_T() [2/2]

◆ cv_from_T() [1/2]

◆ cv_from_T() [2/2]

◆ cv_from_T_v()

◆ cv_from_v_e() [1/3]

◆ cv_from_v_e() [2/3]

◆ cv_from_v_e() [3/3]

◆ e_from_p_T() [1/3]

◆ e_from_p_T() [2/3]

◆ e_from_p_T() [3/3]

◆ e_from_T()

◆ e_from_T_v() [1/3]

◆ e_from_T_v() [2/3]

◆ e_from_T_v() [3/3]

◆ e_from_v_h() [1/2]

◆ e_from_v_h() [2/2]

◆ e_spndl_from_v()

◆ execute()

◆ finalize()

◆ fluidName()

◆ g_from_v_e()

◆ gamma_from_p_T() [1/3]

◆ gamma_from_p_T() [2/3]

◆ gamma_from_p_T() [3/3]

◆ gamma_from_v_e() [1/3]

◆ gamma_from_v_e() [2/3]

◆ gamma_from_v_e() [3/3]

◆ getNaN() [1/2]

◆ getNaN() [2/2]

◆ getNaNVector() [1/2]

◆ getNaNVector() [2/2]

◆ h_from_p_T() [1/2]

◆ h_from_p_T() [2/2]

◆ h_from_T()

◆ h_from_T_v() [1/2]

◆ h_from_T_v() [2/2]

◆ henryCoefficients()

◆ initialize()

◆ initialSetup()

◆ k_from_p_T() [1/2]

◆ k_from_p_T() [2/2]

◆ k_from_v_e() [1/2]

◆ k_from_v_e() [2/2]

◆ makeZeroAndOne() [1/2]