Two-phase stiffened gas fluid properties. More...

#include <StiffenedGasTwoPhaseFluidProperties.h>

Inheritance diagram for StiffenedGasTwoPhaseFluidProperties:

Public Types
typedef DataFileName	DataFileParameterType

Public Member Functions
	StiffenedGasTwoPhaseFluidProperties (const InputParameters &parameters)

virtual Real	p_critical () const override
	Returns the critical pressure. More...

virtual Real	T_triple () const override
	Returns the triple-point temperature. More...

virtual Real	L_fusion () const override
	Returns the latent heat of fusion. More...

virtual Real	T_sat (Real pressure) const override
	Computes the saturation temperature at a pressure. More...

virtual Real	p_sat (Real temperature) const override
	Computes the saturation pressure at a temperature. More...

virtual Real	dT_sat_dp (Real pressure) const override
	Computes dT/dp along the saturation line. More...

virtual Real	sigma_from_T (Real T) const override
	Computes surface tension sigma of saturated liquid in contact with saturated vapor. More...

virtual Real	dsigma_dT_from_T (Real T) const override
	Computes dsigma/dT along the saturation line. More...

virtual bool	supportsPhaseChange () const override
	Returns true if phase change is supported, otherwise false. More...

virtual const UserObjectName &	getLiquidName () const
	Returns the name of the liquid single-phase fluid properties object. More...

virtual const UserObjectName &	getVaporName () const
	Returns the name of the vapor single-phase fluid properties object. More...

virtual ADReal	T_sat (const ADReal &p) const

virtual ADReal	p_sat (const ADReal &T) const

virtual Real	h_lat (Real p, Real T) const
	Computes latent heat of vaporization. More...

virtual ADReal	h_lat (const ADReal &p, const ADReal &T) const

virtual ADReal	sigma_from_T (const ADReal &T) const

virtual void	execute () final

virtual void	initialize () final

virtual void	finalize () final

virtual void	threadJoin (const UserObject &) final

virtual void	subdomainSetup () final

bool	needThreadedCopy () const override final

SubProblem &	getSubProblem () const

bool	shouldDuplicateInitialExecution () const

virtual Real	spatialValue (const Point &) const

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

void	gatherSum (T &value)

void	gatherMax (T &value)

void	gatherMin (T &value)

void	gatherProxyValueMax (T1 &proxy, T2 &value)

void	gatherProxyValueMin (T1 &proxy, T2 &value)

void	setPrimaryThreadCopy (UserObject *primary)

UserObject *	primaryThreadCopy ()

std::set< UserObjectName >	getDependObjects () const

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

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

unsigned int	systemNumber () const

virtual bool	enabled () const

std::shared_ptr< MooseObject >	getSharedPtr ()

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

bool	isKokkosObject (IsKokkosObjectKey &&) 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	mooseInfo (Args &&... args) const

void	callMooseError (std::string msg, const bool with_prefix, const hit::Node *node=nullptr) const

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

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

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

virtual void	initialSetup ()

virtual void	timestepSetup ()

virtual void	jacobianSetup ()

virtual void	residualSetup ()

virtual void	customSetup (const ExecFlagType &)

const ExecFlagEnum &	getExecuteOnEnum () const

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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)

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

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)

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

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

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

Public Attributes
	usingCombinedWarningSolutionWarnings

const ConsoleStream	_console

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

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
Real	compute_p_sat (const Real &T) const
	Computes saturation pressure value using Newton solve. More...

virtual void	addPostprocessorDependencyHelper (const PostprocessorName &name) const override

virtual void	addVectorPostprocessorDependencyHelper (const VectorPostprocessorName &name) const override

virtual void	addUserObjectDependencyHelper (const UserObject &uo) const override

void	addReporterDependencyHelper (const ReporterName &reporter_name) override

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. More...

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

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

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. More...

Static Protected Member Functions
static std::string	meshPropertyName (const std::string &data_name, const std::string &prefix)

Protected Attributes
const Real	_gamma_liquid

const Real	_cv_liquid

const Real	_cp_liquid

const Real	_q_liquid

const Real	_p_inf_liquid

const Real	_q_prime_liquid

const Real	_gamma_vapor

const Real	_cv_vapor

const Real	_cp_vapor

const Real	_q_vapor

const Real	_p_inf_vapor

const Real	_q_prime_vapor

const Real &	_T_c
	critical temperature More...

const Real &	_p_c
	critical pressure More...

const Real &	_T_triple
	Triple-point temperature. More...

const Real &	_L_fusion
	Latent heat of fusion. More...

const Real &	_sigma_A
	'A' constant used in surface tension correlation More...

const Real &	_sigma_B
	'B' constant used in surface tension correlation More...

const Real &	_sigma_C
	'C' constant used in surface tension correlation More...

const Real &	_T_sat_min
	Minimum temperature value in saturation curve. More...

const Real &	_T_sat_max
	Maximum temperature value in saturation curve. More...

const Real &	_p_sat_guess
	Initial guess for saturation pressure Newton solve. More...

const unsigned int &	_n_sat_samples
	Number of samples to take in saturation curve. More...

const Real	_dT_sat
	Temperature increments on saturation curve. More...

const Real	_A

const Real	_B

const Real	_C

const Real	_D

const Real	_newton_tol
	Newton solve tolerance. More...

const unsigned int	_newton_max_iter
	Newton max number of iterations. More...

std::vector< Real >	_T_vec

std::vector< Real >	_p_sat_vec

LinearInterpolation	_ipol_temp

LinearInterpolation	_ipol_pressure

const UserObjectName	_liquid_name
	The name of the user object that provides liquid phase fluid properties. More...

const UserObjectName	_vapor_name
	The name of the user object that provides vapor phase fluid properties. More...

const SinglePhaseFluidProperties *	_fp_liquid
	The user object that provides liquid phase fluid properties. More...

const SinglePhaseFluidProperties *	_fp_vapor
	The user object that provides vapor phase fluid properties. More...

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

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

SubProblem &	_subproblem

FEProblemBase &	_fe_problem

SystemBase &	_sys

const THREAD_ID	_tid

Assembly &	_assembly

const Moose::CoordinateSystemType &	_coord_sys

const bool	_duplicate_initial_execution

std::set< std::string >	_depend_uo

const bool &	_enabled

MooseApp &	_app

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? More...

Static Protected Attributes
static const std::string	_interpolated_old

static const std::string	_interpolated_older

Detailed Description

Two-phase stiffened gas fluid properties.

Definition at line 24 of file StiffenedGasTwoPhaseFluidProperties.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

◆ StiffenedGasTwoPhaseFluidProperties()

StiffenedGasTwoPhaseFluidProperties::StiffenedGasTwoPhaseFluidProperties ( const InputParameters & parameters )

Definition at line 65 of file StiffenedGasTwoPhaseFluidProperties.C.

   : TwoPhaseFluidProperties(parameters),
     NaNInterface(this),
 
     _gamma_liquid(getParam<Real>("gamma_liquid")),
     _cv_liquid(getParam<Real>("cv_liquid")),
     _cp_liquid(_gamma_liquid * _cv_liquid),
     _q_liquid(getParam<Real>("q_liquid")),
     _p_inf_liquid(getParam<Real>("p_inf_liquid")),
     _q_prime_liquid(getParam<Real>("q_prime_liquid")),
 
     _gamma_vapor(getParam<Real>("gamma_vapor")),
     _cv_vapor(getParam<Real>("cv_vapor")),
     _cp_vapor(_gamma_vapor * _cv_vapor),
     _q_vapor(getParam<Real>("q_vapor")),
     _p_inf_vapor(getParam<Real>("p_inf_vapor")),
     _q_prime_vapor(getParam<Real>("q_prime_vapor")),
 
     _T_c(getParam<Real>("T_c")),
     _p_c(getParam<Real>("p_c")),
     _T_triple(getParam<Real>("T_triple")),
     _L_fusion(getParam<Real>("L_fusion")),
 
     _sigma_A(getParam<Real>("sigma_A")),
     _sigma_B(getParam<Real>("sigma_B")),
     _sigma_C(getParam<Real>("sigma_C")),
 
     _T_sat_min(getParam<Real>("T_sat_min")),
     _T_sat_max(getParam<Real>("T_sat_max")),
     _p_sat_guess(getParam<Real>("p_sat_guess")),
     _n_sat_samples(getParam<unsigned int>("n_sat_samples")),
     _dT_sat((_T_sat_max - _T_sat_min) / (_n_sat_samples - 1)),
 
     _A((_cp_liquid - _cp_vapor + _q_prime_vapor - _q_prime_liquid) / (_cp_vapor - _cv_vapor)),
     _B((_q_liquid - _q_vapor) / (_cp_vapor - _cv_vapor)),
     _C((_cp_vapor - _cp_liquid) / (_cp_vapor - _cv_vapor)),
     _D((_cp_liquid - _cv_liquid) / (_cp_vapor - _cv_vapor)),
 
     _newton_tol(1.0e-8),
     _newton_max_iter(200)
 {
   if (_tid == 0)
   {
     std::string class_name = "StiffenedGasFluidProperties";
     InputParameters params = _app.getFactory().getValidParams(class_name);
     params.set<MooseEnum>("emit_on_nan") = getParam<MooseEnum>("emit_on_nan");
     params.set<bool>("allow_nonphysical_states") = false;
     params.set<Real>("gamma") = _gamma_liquid;
     params.set<Real>("cv") = _cv_liquid;
     params.set<Real>("q") = _q_liquid;
     params.set<Real>("p_inf") = _p_inf_liquid;
     params.set<Real>("q_prime") = _q_prime_liquid;
     params.set<Real>("k") = getParam<Real>("k_liquid");
     params.set<Real>("mu") = getParam<Real>("mu_liquid");
     params.set<Real>("M") = getParam<Real>("M_liquid");
     _fe_problem.addUserObject(class_name, _liquid_name, params);
   }
   _fp_liquid = &_fe_problem.getUserObject<SinglePhaseFluidProperties>(_liquid_name, _tid);
 
   if (_tid == 0)
   {
     std::string class_name = "StiffenedGasFluidProperties";
     InputParameters params = _app.getFactory().getValidParams(class_name);
     params.set<MooseEnum>("emit_on_nan") = getParam<MooseEnum>("emit_on_nan");
     params.set<bool>("allow_nonphysical_states") = false;
     params.set<Real>("gamma") = _gamma_vapor;
     params.set<Real>("cv") = _cv_vapor;
     params.set<Real>("q") = _q_vapor;
     params.set<Real>("p_inf") = _p_inf_vapor;
     params.set<Real>("q_prime") = _q_prime_vapor;
     params.set<Real>("k") = getParam<Real>("k_vapor");
     params.set<Real>("mu") = getParam<Real>("mu_vapor");
     params.set<Real>("M") = getParam<Real>("M_vapor");
     params.set<Real>("T_c") = getParam<Real>("T_c");
     params.set<Real>("rho_c") = getParam<Real>("rho_c");
     params.set<Real>("e_c") = getParam<Real>("e_c");
     _fe_problem.addUserObject(class_name, _vapor_name, params);
   }
   _fp_vapor = &_fe_problem.getUserObject<SinglePhaseFluidProperties>(_vapor_name, _tid);
 
   // calculate the saturation curve p(T) and store the data in two vectors for re-use
   {
     _T_vec.resize(_n_sat_samples);
     _p_sat_vec.resize(_n_sat_samples);
 
     // loop over sample temperatures, starting with the minimum
     Real T = _T_sat_min;
     for (unsigned int i = 0; i < _n_sat_samples; i++)
     {
       _T_vec[i] = T;
       _p_sat_vec[i] = compute_p_sat(T);
 
       // increment sample temperature
       T += _dT_sat;
     }
   }
 
   _ipol_pressure.setData(_T_vec, _p_sat_vec);
   _ipol_temp.setData(_p_sat_vec, _T_vec);
 }

Member Function Documentation

◆ compute_p_sat()

Real StiffenedGasTwoPhaseFluidProperties::compute_p_sat ( const Real & T ) const

protected

Computes saturation pressure value using Newton solve.

The process for determining the saturation pressure is given in the following reference:

Ray A. Berry, Richard Saurel, Olivier LeMetayer The discrete equation method (DEM) for fully compressible, two-phase flow in ducts of spatially varying cross-section Nuclear Engineering and Design 240 (2010) p. 3797-3818

The nonlinear equation to be solved is given by Equation (38) of this reference; it is obtained by reasoning that at thermodynamic equilibrium, the Gibbs free enthalpy of the phases must be equal at the interface.

Parameters

[in] T temperature at which saturation pressure is to be computed

Definition at line 198 of file StiffenedGasTwoPhaseFluidProperties.C.

Referenced by StiffenedGasTwoPhaseFluidProperties().

 {
   Real p_sat = _p_sat_guess;
   Real f_norm = 1.0e5;
   unsigned int iter = 1;
   while (std::fabs(f_norm / p_sat) > _newton_tol)
   {
     // check for maximum iteration
     if (iter > _newton_max_iter)
       mooseError("StiffenedGasTwoPhaseFluidProperties::compute_p_sat: ",
                  "Newton solve did not converge after ",
                  _newton_max_iter,
                  " iterations.");
 
     // evaluate nonlinear residual and its derivative w.r.t. p_sat
     Real f = std::log(p_sat + _p_inf_vapor) - _A - _B / T - _C * std::log(T) -
              _D * std::log(p_sat + _p_inf_liquid);
     Real f_prime = 1.0 / (p_sat + _p_inf_vapor) - _D / (p_sat + _p_inf_liquid);
 
     // take Newton step
     f_norm = f / f_prime;
     p_sat -= f_norm;
 
     iter++;
   }
   return p_sat;
 }

◆ dsigma_dT_from_T()

Real StiffenedGasTwoPhaseFluidProperties::dsigma_dT_from_T ( Real T ) const

overridevirtual

Computes dsigma/dT along the saturation line.

Parameters

[in] T temperature (K)

Reimplemented from TwoPhaseFluidProperties.

Definition at line 240 of file StiffenedGasTwoPhaseFluidProperties.C.

 {
   const Real aux = 1.0 - T / _T_c;
   const Real daux_dT = -1.0 / _T_c;
   const Real dsigma_daux =
       _sigma_A * (_sigma_B * std::pow(aux, _sigma_B - 1.0) * (1.0 - _sigma_C * aux) -
                   _sigma_C * std::pow(aux, _sigma_B));
   return dsigma_daux * daux_dT;
 }

◆ dT_sat_dp()

Real StiffenedGasTwoPhaseFluidProperties::dT_sat_dp ( Real p ) const

overridevirtual

Computes dT/dp along the saturation line.

Parameters

[in] p pressure

Implements TwoPhaseFluidProperties.

Definition at line 227 of file StiffenedGasTwoPhaseFluidProperties.C.

 {
   return _ipol_temp.sampleDerivative(pressure);
 }

◆ 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 {}

◆ getLiquidName()

virtual const UserObjectName& TwoPhaseFluidProperties::getLiquidName ( ) const

inlinevirtualinherited

Returns the name of the liquid single-phase fluid properties object.

Reimplemented in TwoPhaseNCGFluidProperties.

Definition at line 32 of file TwoPhaseFluidProperties.h.

Referenced by StiffenedGasTwoPhaseFluidPropertiesTest::buildObjects(), TwoPhaseNCGFluidProperties::getLiquidName(), and SaturationDensityFunction::initialSetup().

32 { return _liquid_name; }

TwoPhaseFluidProperties::_liquid_name

const UserObjectName _liquid_name

The name of the user object that provides liquid phase fluid properties.

Definition: TwoPhaseFluidProperties.h:109

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

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

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

NaNInterface::getNaN

Real getNaN() const

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

Definition: NaNInterface.h:46

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

Referenced by NaNInterfaceTestFluidProperties::returnNaNVector().

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

◆ 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(""));
   }

◆ getVaporName()

virtual const UserObjectName& TwoPhaseFluidProperties::getVaporName ( ) const

inlinevirtualinherited

Returns the name of the vapor single-phase fluid properties object.

Reimplemented in TwoPhaseNCGFluidProperties.

Definition at line 37 of file TwoPhaseFluidProperties.h.

Referenced by StiffenedGasTwoPhaseFluidPropertiesTest::buildObjects(), TwoPhaseNCGFluidProperties::getVaporName(), and SaturationDensityFunction::initialSetup().

37 { return _vapor_name; }

TwoPhaseFluidProperties::_vapor_name

const UserObjectName _vapor_name

The name of the user object that provides vapor phase fluid properties.

Definition: TwoPhaseFluidProperties.h:111

◆ h_lat() [1/2]

Real TwoPhaseFluidProperties::h_lat	(	Real	p,
		Real	T
	)		const

virtualinherited

Computes latent heat of vaporization.

Parameters

p	pressure
T	temperature

Reimplemented in TwoPhaseNCGFluidProperties.

Definition at line 92 of file TwoPhaseFluidProperties.C.

Referenced by FluidPropertiesInterrogator::compute2Phase(), and TwoPhaseNCGFluidProperties::h_lat().

 {
   return _fp_vapor->h_from_p_T(p, T) - _fp_liquid->h_from_p_T(p, T);
 }

◆ h_lat() [2/2]

ADReal TwoPhaseFluidProperties::h_lat	(	const ADReal &	p,
		const ADReal &	T
	)		const

virtualinherited

Reimplemented in TwoPhaseNCGFluidProperties.

Definition at line 98 of file TwoPhaseFluidProperties.C.

 {
   return _fp_vapor->h_from_p_T(p, T) - _fp_liquid->h_from_p_T(p, T);
 }

◆ initialize()

virtual void FluidProperties::initialize ( )

inlinefinalvirtualinherited

Implements ThreadedGeneralUserObject.

Definition at line 34 of file FluidProperties.h.

34 {}

◆ L_fusion()

Real StiffenedGasTwoPhaseFluidProperties::L_fusion ( ) const

overridevirtual

Returns the latent heat of fusion.

Reimplemented from TwoPhaseFluidProperties.

Definition at line 180 of file StiffenedGasTwoPhaseFluidProperties.C.

 {
   return _L_fusion;
 }

◆ p_critical()

Real StiffenedGasTwoPhaseFluidProperties::p_critical ( ) const

overridevirtual

Returns the critical pressure.

Implements TwoPhaseFluidProperties.

Definition at line 168 of file StiffenedGasTwoPhaseFluidProperties.C.

 {
   return _p_c;
 }

◆ p_sat() [1/2]

Real StiffenedGasTwoPhaseFluidProperties::p_sat ( Real T ) const

overridevirtual

Computes the saturation pressure at a temperature.

Parameters

[in] T temperature

Implements TwoPhaseFluidProperties.

Definition at line 192 of file StiffenedGasTwoPhaseFluidProperties.C.

Referenced by compute_p_sat().

 {
   return _ipol_pressure.sample(temperature);
 }

◆ p_sat() [2/2]

ADReal TwoPhaseFluidProperties::p_sat ( const ADReal & T ) const

virtualinherited

Reimplemented in TwoPhaseNCGFluidProperties, and TwoPhaseFluidPropertiesIndependent.

Definition at line 80 of file TwoPhaseFluidProperties.C.

 {
   const Real T_real = T.value();
   const Real p_sat_real = p_sat(T_real);
   const Real dp_sat_dT_real = 1.0 / dT_sat_dp(p_sat_real);
 
   ADReal p_sat = p_sat_real;
   p_sat.derivatives() = T.derivatives() * dp_sat_dT_real;
   return p_sat;
 }

◆ sigma_from_T() [1/2]

Real StiffenedGasTwoPhaseFluidProperties::sigma_from_T ( Real T ) const

overridevirtual

Computes surface tension sigma of saturated liquid in contact with saturated vapor.

Parameters

T	temperature

Reimplemented from TwoPhaseFluidProperties.

Definition at line 233 of file StiffenedGasTwoPhaseFluidProperties.C.

 {
   const Real aux = 1.0 - T / _T_c;
   return _sigma_A * std::pow(aux, _sigma_B) * (1.0 - _sigma_C * aux);
 }

◆ sigma_from_T() [2/2]

ADReal TwoPhaseFluidProperties::sigma_from_T ( const ADReal & T ) const

virtualinherited

Reimplemented in TwoPhaseFluidPropertiesIndependent.

Definition at line 115 of file TwoPhaseFluidProperties.C.

 {
   const Real T_real = T.value();
   const Real sigma_real = sigma_from_T(T_real);
   const Real dsigma_dT = dsigma_dT_from_T(T_real);
 
   ADReal sigma = sigma_real;
   sigma.derivatives() = T.derivatives() * dsigma_dT;
   return sigma;
 }

◆ subdomainSetup()

virtual void FluidProperties::subdomainSetup ( )

inlinefinalvirtualinherited

Reimplemented from ThreadedGeneralUserObject.

Definition at line 38 of file FluidProperties.h.

38 {}

◆ supportsPhaseChange()

virtual bool StiffenedGasTwoPhaseFluidProperties::supportsPhaseChange ( ) const

inlineoverridevirtual

Returns true if phase change is supported, otherwise false.

Implements TwoPhaseFluidProperties.

Definition at line 38 of file StiffenedGasTwoPhaseFluidProperties.h.

38 { return true; }

◆ T_sat() [1/2]

Real StiffenedGasTwoPhaseFluidProperties::T_sat ( Real p ) const

overridevirtual

Computes the saturation temperature at a pressure.

Parameters

[in] p pressure

Implements TwoPhaseFluidProperties.

Definition at line 186 of file StiffenedGasTwoPhaseFluidProperties.C.

 {
   return _ipol_temp.sample(pressure);
 }

◆ T_sat() [2/2]

ADReal TwoPhaseFluidProperties::T_sat ( const ADReal & p ) const

virtualinherited

Reimplemented in TwoPhaseNCGFluidProperties, and TwoPhaseFluidPropertiesIndependent.

Definition at line 68 of file TwoPhaseFluidProperties.C.

 {
   const Real p_real = p.value();
   const Real T_sat_real = T_sat(p_real);
   const Real dT_sat_dp_real = dT_sat_dp(p_real);
 
   ADReal T_sat = T_sat_real;
   T_sat.derivatives() = p.derivatives() * dT_sat_dp_real;
   return T_sat;
 }

◆ T_triple()

Real StiffenedGasTwoPhaseFluidProperties::T_triple ( ) const

overridevirtual

Returns the triple-point temperature.

Reimplemented from TwoPhaseFluidProperties.

Definition at line 174 of file StiffenedGasTwoPhaseFluidProperties.C.

 {
   return _T_triple;
 }

◆ threadJoin()

virtual void FluidProperties::threadJoin ( const UserObject & )

inlinefinalvirtualinherited

Reimplemented from ThreadedGeneralUserObject.

Definition at line 37 of file FluidProperties.h.

37 {}

◆ validParams()

InputParameters StiffenedGasTwoPhaseFluidProperties::validParams ( )

static

Definition at line 16 of file StiffenedGasTwoPhaseFluidProperties.C.

 {
   InputParameters params = TwoPhaseFluidProperties::validParams();
   params += NaNInterface::validParams();
 
   params.addParam<Real>("gamma_liquid", 2.35, "Liquid heat capacity ratio");
   params.addParam<Real>("cv_liquid", 1816.0, "Liquid isochoric specific heat capacity");
   params.addParam<Real>("q_liquid", -1.167e6, "Liquid reference specific internal energy");
   params.addParam<Real>("p_inf_liquid", 1.0e9, "Liquid stiffness pressure");
   params.addParam<Real>("q_prime_liquid", 0, "Liquid reference specific entropy");
   params.addParam<Real>("k_liquid", 0.5, "Liquid thermal conductivity");
   params.addParam<Real>("mu_liquid", 281.8e-6, "Liquid dynamic viscosity");
   params.addParam<Real>("M_liquid", 0.01801488, "Liquid molar mass");
 
   params.addParam<Real>("gamma_vapor", 1.43, "Vapor heat capacity ratio");
   params.addParam<Real>("cv_vapor", 1040.0, "Vapor isochoric specific heat capacity");
   params.addParam<Real>("q_vapor", 2.03e6, "Vapor reference specific internal energy");
   params.addParam<Real>("p_inf_vapor", 0.0, "Vapor stiffness pressure");
   params.addParam<Real>("q_prime_vapor", -2.3e4, "Vapor reference specific entropy");
   params.addParam<Real>("k_vapor", 0.026, "Vapor thermal conductivity");
   params.addParam<Real>("mu_vapor", 134.4e-7, "Vapor dynamic viscosity");
   params.addParam<Real>("M_vapor", 0.01801488, "Vapor molar mass");
 
   params.addParam<Real>("T_c", 647.096, "Critical temperature [K]");
   params.addParam<Real>("p_c", 22.09e6, "Critical pressure [Pa]");
   params.addParam<Real>("rho_c", 322.0, "Critical density [kg/m^3]");
   params.addParam<Real>("e_c", 2702979.84310559, "Critical specific internal energy [J/kg]");
   params.addParam<Real>("T_triple", 273.16, "Triple-point temperature [K]");
   params.addParam<Real>("L_fusion", 0.334, "Latent heat of fusion [J/kg]");
 
   params.addParam<Real>(
       "sigma_A", 0.2358, "'A' constant used in surface tension correlation [N/m]");
   params.addParam<Real>("sigma_B", 1.256, "'B' constant used in surface tension correlation");
   params.addParam<Real>("sigma_C", 0.625, "'C' constant used in surface tension correlation");
 
   // Default values correspond to water, from the freezing point to critical point, with 1 K
   // increments
   params.addParam<Real>("T_sat_min", 274.0, "Minimum temperature value in saturation curve [K]");
   params.addParam<Real>("T_sat_max", 647.0, "Maximum temperature value in saturation curve [K]");
   params.addParam<Real>(
       "p_sat_guess", 611.0, "Initial guess for saturation pressure Newton solve [Pa]");
   params.addParam<unsigned int>(
       "n_sat_samples", 374, "Number of samples to take in saturation curve");
 
   params.addClassDescription("Two-phase stiffened gas fluid properties");
 
   return params;
 }