https://mooseframework.inl.gov
Loading...
Searching...
No Matches
Public Types | Public Member Functions | Static Public Member Functions | Public Attributes | Static Public Attributes | Protected Member Functions | Static Protected Member Functions | Protected Attributes | Static Protected Attributes | Private Member Functions | Static Private Member Functions | Private Attributes | List of all members
PorousFlowWaterVapor Class Reference

Specialized class for water and vapor mixture using pressure and enthalpy. More...

#include <PorousFlowWaterVapor.h>

Inheritance diagram for PorousFlowWaterVapor:
[legend]

Public Types

typedef DataFileName DataFileParameterType
 

Public Member Functions

 PorousFlowWaterVapor (const InputParameters &parameters)
 
virtual std::string fluidStateName () const override
 Name of FluidState.
 
void thermophysicalProperties (Real pressure, Real enthalpy, unsigned int qp, std::vector< FluidStateProperties > &fsp) const override
 Determines the complete thermophysical state of the system for a given set of primary variables.
 
void thermophysicalProperties (const ADReal &pressure, const ADReal &enthalpy, unsigned int qp, std::vector< FluidStateProperties > &fsp) const override
 
void thermophysicalProperties (const ADReal &pressure, const ADReal &enthalpy, unsigned int qp, FluidStatePhaseEnum &phase_state, std::vector< FluidStateProperties > &fsp) const
 
unsigned int aqueousComponentIndex () const
 The index of the aqueous fluid component.
 
unsigned int gasComponentIndex () const
 The index of the gas fluid component.
 
unsigned int getPressureIndex () const
 
unsigned int getEnthalpyIndex () const
 
void initialize () final
 
void execute () final
 
void finalize () final
 
unsigned int numPhases () const
 The maximum number of phases in this model.
 
unsigned int numComponents () const
 The maximum number of components in this model.
 
unsigned int aqueousPhaseIndex () const
 The index of the aqueous phase.
 
unsigned int gasPhaseIndex () const
 The index of the gas phase.
 
void clearFluidStateProperties (std::vector< FluidStateProperties > &fsp) const
 Clears the contents of the FluidStateProperties data structure.
 
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 needThreadedCopy () 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 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 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
 
virtual void threadJoin (const UserObject &) override
 
virtual void threadJoin (const UserObject &) override
 
virtual void subdomainSetup () override
 
virtual void subdomainSetup () override
 
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 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 Member Functions

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)
 

Static Protected Member Functions

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

Protected Attributes

const SinglePhaseFluidProperties & _water_fp
 Fluid properties UserObject for water.
 
const Real _Mh2o
 Molar mass of water (kg/mol)
 
const Real _p_triple
 Triple point pressure of water (Pa)
 
const Real _p_critical
 Critical pressure of water (Pa)
 
const Real _T_triple
 Triple point temperature of water (K)
 
const Real _T_critical
 Critical temperature of water (K)
 
const unsigned int _fluid_component
 Fluid component number (only one fluid component in all phases)
 
const unsigned int _pidx
 Index of derivative wrt pressure.
 
const unsigned int _hidx
 Index of derivative wrt enthalpy.
 
const Real _dT
 Perturbation applied to saturation temperature to move to gas/liquid phase.
 
unsigned int _num_phases
 Number of phases.
 
unsigned int _num_components
 Number of components.
 
const unsigned int _aqueous_phase_number
 Phase number of the aqueous phase.
 
unsigned int _gas_phase_number
 Phase number of the gas phase.
 
const Real _R
 Universal gas constant (J/mol/K)
 
const Real _T_c2k
 Conversion from C to K.
 
const PorousFlowCapillaryPressure & _pc
 Capillary pressure UserObject.
 
FluidStateProperties _empty_fsp
 Empty FluidStateProperties object.
 
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
 

Static Protected Attributes

static const std::string _interpolated_old
 
static const std::string _interpolated_older
 

Private Member Functions

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
 

Detailed Description

Specialized class for water and vapor mixture using pressure and enthalpy.

Note: It is important to note that this class must be used in a non-isothernal setting (i.e., with both fluid and heat transport) in order for the problem to be well posed.

Definition at line 22 of file PorousFlowWaterVapor.h.

Constructor & Destructor Documentation

◆ PorousFlowWaterVapor()

PorousFlowWaterVapor::PorousFlowWaterVapor ( const InputParameters &  parameters)

Definition at line 26 of file PorousFlowWaterVapor.C.

27 : PorousFlowFluidStateSingleComponentBase(parameters),
28 _water_fp(getUserObject<SinglePhaseFluidProperties>("water_fp")),
29 _Mh2o(_water_fp.molarMass()),
30 _p_triple(_water_fp.triplePointPressure()),
31 _p_critical(_water_fp.criticalPressure()),
32 _T_triple(_water_fp.triplePointTemperature()),
33 _T_critical(_water_fp.criticalTemperature())
34{
35 // Check that the correct FluidProperties UserObjects have been provided
36 if (_water_fp.fluidName() != "water")
37 paramError("water_fp", "A valid water FluidProperties UserObject must be provided in water_fp");
38
39 // Set the number of phases and components, and their indexes
40 _num_phases = 2;
41 _num_components = 1;
42 _gas_phase_number = 1 - _aqueous_phase_number;
43
44 // Check that _aqueous_phase_number is <= total number of phases
45 if (_aqueous_phase_number >= _num_phases)
46 paramError("liquid_phase_number",
47 "This value is larger than the possible number of phases ",
48 _num_phases);
49
50 // Check that _fluid_component is <= total number of fluid components
51 if (_fluid_component >= _num_components)
52 paramError("fluid_component",
53 "This value is larger than the possible number of fluid components",
54 _num_components);
55
56 _empty_fsp = FluidStateProperties(_num_components);
57}
GeneralUserObject::parameters
const InputParameters & parameters() const
GeneralUserObject::paramError
void paramError(const std::string &param, Args... args) const
PorousFlowFluidStateBase::_num_components
unsigned int _num_components
Number of components.
Definition PorousFlowFluidStateBase.h:102
PorousFlowFluidStateBase::_gas_phase_number
unsigned int _gas_phase_number
Phase number of the gas phase.
Definition PorousFlowFluidStateBase.h:106
PorousFlowFluidStateBase::_aqueous_phase_number
const unsigned int _aqueous_phase_number
Phase number of the aqueous phase.
Definition PorousFlowFluidStateBase.h:104
PorousFlowFluidStateBase::_empty_fsp
FluidStateProperties _empty_fsp
Empty FluidStateProperties object.
Definition PorousFlowFluidStateBase.h:114
PorousFlowFluidStateBase::_num_phases
unsigned int _num_phases
Number of phases.
Definition PorousFlowFluidStateBase.h:100
PorousFlowFluidStateSingleComponentBase
Base class for miscible multiphase flow classes with a single fluid component using a pressure and en...
Definition PorousFlowFluidStateSingleComponentBase.h:19
PorousFlowFluidStateSingleComponentBase::_fluid_component
const unsigned int _fluid_component
Fluid component number (only one fluid component in all phases)
Definition PorousFlowFluidStateSingleComponentBase.h:61
PorousFlowWaterVapor::_T_triple
const Real _T_triple
Triple point temperature of water (K)
Definition PorousFlowWaterVapor.h:57
PorousFlowWaterVapor::_Mh2o
const Real _Mh2o
Molar mass of water (kg/mol)
Definition PorousFlowWaterVapor.h:51
PorousFlowWaterVapor::_p_triple
const Real _p_triple
Triple point pressure of water (Pa)
Definition PorousFlowWaterVapor.h:53
PorousFlowWaterVapor::_water_fp
const SinglePhaseFluidProperties & _water_fp
Fluid properties UserObject for water.
Definition PorousFlowWaterVapor.h:49
PorousFlowWaterVapor::_T_critical
const Real _T_critical
Critical temperature of water (K)
Definition PorousFlowWaterVapor.h:59
PorousFlowWaterVapor::_p_critical
const Real _p_critical
Critical pressure of water (Pa)
Definition PorousFlowWaterVapor.h:55
SinglePhaseFluidProperties::triplePointPressure
virtual Real triplePointPressure() const
Triple point pressure.
Definition SinglePhaseFluidProperties.C:326
SinglePhaseFluidProperties::criticalPressure
virtual Real criticalPressure() const
Critical pressure.
Definition SinglePhaseFluidProperties.C:302
SinglePhaseFluidProperties::molarMass
virtual Real molarMass() const
Molar mass [kg/mol].
Definition SinglePhaseFluidProperties.C:290
SinglePhaseFluidProperties::criticalTemperature
virtual Real criticalTemperature() const
Critical temperature.
Definition SinglePhaseFluidProperties.C:308
SinglePhaseFluidProperties::triplePointTemperature
virtual Real triplePointTemperature() const
Triple point temperature.
Definition SinglePhaseFluidProperties.C:332
SinglePhaseFluidProperties::fluidName
e e e e s T T T T T rho v v T e p T T virtual T std::string fluidName() const
Fluid name.
Definition SinglePhaseFluidProperties.C:296
FluidStateProperties
AD data structure to pass calculated thermophysical properties.
Definition PorousFlowFluidStateBase.h:25

Member Function Documentation

◆ aqueousComponentIndex()

unsigned int PorousFlowFluidStateSingleComponentBase::aqueousComponentIndex ( ) const
inlineinherited

The index of the aqueous fluid component.

Returns
aqueous fluid component number

Definition at line 29 of file PorousFlowFluidStateSingleComponentBase.h.

29{ return _fluid_component; };

◆ aqueousPhaseIndex()

unsigned int PorousFlowFluidStateBase::aqueousPhaseIndex ( ) const
inlineinherited

The index of the aqueous phase.

Returns
aqueous phase number

Definition at line 79 of file PorousFlowFluidStateBase.h.

79{ return _aqueous_phase_number; };

◆ clearFluidStateProperties()

void PorousFlowFluidStateBase::clearFluidStateProperties ( std::vector< FluidStateProperties > &  fsp) const
inherited

Clears the contents of the FluidStateProperties data structure.

Parameters
[out]fspFluidStateProperties data structure with all data initialized to 0

Definition at line 33 of file PorousFlowFluidStateBase.C.

34{
35 std::fill(fsp.begin(), fsp.end(), _empty_fsp);
36}

Referenced by PorousFlowBrineCO2::thermophysicalProperties(), and PorousFlowWaterNCG::thermophysicalProperties().

◆ execute()

void PorousFlowFluidStateBase::execute ( )
inlinefinalvirtualinherited

Implements GeneralUserObject.

Definition at line 60 of file PorousFlowFluidStateBase.h.

60{};

◆ finalize()

void PorousFlowFluidStateBase::finalize ( )
inlinefinalvirtualinherited

Implements GeneralUserObject.

Definition at line 61 of file PorousFlowFluidStateBase.h.

61{};

◆ fluidStateName()

std::string PorousFlowWaterVapor::fluidStateName ( ) const
overridevirtual

Name of FluidState.

Implements PorousFlowFluidStateBase.

Definition at line 60 of file PorousFlowWaterVapor.C.

61{
62 return "water-vapor";
63}

◆ gasComponentIndex()

unsigned int PorousFlowFluidStateSingleComponentBase::gasComponentIndex ( ) const
inlineinherited

The index of the gas fluid component.

Returns
gas fluid component number

Definition at line 35 of file PorousFlowFluidStateSingleComponentBase.h.

35{ return _fluid_component; };

◆ gasPhaseIndex()

unsigned int PorousFlowFluidStateBase::gasPhaseIndex ( ) const
inlineinherited

The index of the gas phase.

Returns
gas phase number

Definition at line 85 of file PorousFlowFluidStateBase.h.

85{ return _gas_phase_number; };

◆ getEnthalpyIndex()

unsigned int PorousFlowFluidStateSingleComponentBase::getEnthalpyIndex ( ) const
inlineinherited

Definition at line 57 of file PorousFlowFluidStateSingleComponentBase.h.

57{ return _hidx; };
PorousFlowFluidStateSingleComponentBase::_hidx
const unsigned int _hidx
Index of derivative wrt enthalpy.
Definition PorousFlowFluidStateSingleComponentBase.h:65

◆ getPressureIndex()

unsigned int PorousFlowFluidStateSingleComponentBase::getPressureIndex ( ) const
inlineinherited

Definition at line 56 of file PorousFlowFluidStateSingleComponentBase.h.

56{ return _pidx; };
PorousFlowFluidStateSingleComponentBase::_pidx
const unsigned int _pidx
Index of derivative wrt pressure.
Definition PorousFlowFluidStateSingleComponentBase.h:63

◆ initialize()

void PorousFlowFluidStateBase::initialize ( )
inlinefinalvirtualinherited

Implements GeneralUserObject.

Definition at line 59 of file PorousFlowFluidStateBase.h.

59{};

◆ numComponents()

unsigned int PorousFlowFluidStateBase::numComponents ( ) const
inlineinherited

The maximum number of components in this model.

Returns
number of components

Definition at line 73 of file PorousFlowFluidStateBase.h.

73{ return _num_components; };

◆ numPhases()

unsigned int PorousFlowFluidStateBase::numPhases ( ) const
inlineinherited

The maximum number of phases in this model.

Returns
number of phases

Definition at line 67 of file PorousFlowFluidStateBase.h.

67{ return _num_phases; };

Referenced by PorousFlowFluidStateSingleComponentTempl< is_ad >::PorousFlowFluidStateSingleComponentTempl(), and PorousFlowFluidStateTempl< is_ad >::PorousFlowFluidStateTempl().

◆ thermophysicalProperties() [1/3]

void PorousFlowWaterVapor::thermophysicalProperties ( const ADReal &  pressure,
const ADReal &  enthalpy,
unsigned int  qp,
FluidStatePhaseEnum &  phase_state,
std::vector< FluidStateProperties > &  fsp 
) const

Definition at line 92 of file PorousFlowWaterVapor.C.

97{
98 FluidStateProperties & liquid = fsp[_aqueous_phase_number];
99 FluidStateProperties & gas = fsp[_gas_phase_number];
100
101 ADReal Tsat = 0.0;
102 ADReal hl = 0.0;
103 ADReal hv = 0.0;
104
105 // Determine the phase state of the system
106 if (pressure.value() >= _p_triple && pressure.value() <= _p_critical)
107 {
108 // Saturation temperature at the given pressure
109 Tsat = _water_fp.vaporTemperature(pressure);
110
111 // Enthalpy of saturated liquid and saturated vapor
112 hl = _water_fp.h_from_p_T(pressure, Tsat - _dT);
113 hv = _water_fp.h_from_p_T(pressure, Tsat + _dT);
114
115 if (enthalpy.value() < hl.value())
116 phase_state = FluidStatePhaseEnum::LIQUID;
117
118 else if (enthalpy.value() >= hl.value() && enthalpy.value() <= hv.value())
119 phase_state = FluidStatePhaseEnum::TWOPHASE;
120
121 else // h > hv
122 phase_state = FluidStatePhaseEnum::GAS;
123 }
124 else // p.value() > _p_critical
125 {
126 // Check whether the phase point is in the liquid or vapor state
127 const ADReal T = _water_fp.T_from_p_h(pressure, enthalpy);
128
129 if (T.value() <= _T_critical)
130 phase_state = FluidStatePhaseEnum::LIQUID;
131 else
132 {
133 // The supercritical state is treated as a gas
134 phase_state = FluidStatePhaseEnum::GAS;
135 }
136 }
137
138 // Calculate the properties for each phase as required
139 switch (phase_state)
140 {
141 case FluidStatePhaseEnum::GAS:
142 {
143 gas.pressure = pressure + _pc.capillaryPressure(0.0, qp);
144 gas.saturation = 1.0;
145
146 const ADReal T = _water_fp.T_from_p_h(gas.pressure, enthalpy);
147
148 gas.temperature = T;
149 liquid.temperature = T;
150
151 gas.density = _water_fp.rho_from_p_T(gas.pressure, T);
152 gas.viscosity = _water_fp.mu_from_p_T(gas.pressure, T);
153 gas.enthalpy = enthalpy;
154 gas.internal_energy = _water_fp.e_from_p_T(gas.pressure, T);
155 gas.mass_fraction[_fluid_component] = 1.0;
156
157 break;
158 }
159
160 case FluidStatePhaseEnum::LIQUID:
161 {
162 const ADReal T = _water_fp.T_from_p_h(pressure, enthalpy);
163
164 liquid.pressure = pressure;
165 liquid.temperature = T;
166 liquid.density = _water_fp.rho_from_p_T(pressure, T);
167 liquid.viscosity = _water_fp.mu_from_p_T(pressure, T);
168 liquid.enthalpy = enthalpy;
169 liquid.internal_energy = _water_fp.e_from_p_T(pressure, T);
170 liquid.saturation = 1.0;
171 liquid.mass_fraction[_fluid_component] = 1.0;
172
173 break;
174 }
175
176 case FluidStatePhaseEnum::TWOPHASE:
177 {
178 // Latent heat of vaporization
179 const ADReal hvl = hv - hl;
180
181 // Vapor quality
182 const ADReal X = (enthalpy - hl) / hvl;
183
184 // Perturbed saturation temperature to ensure that the correct
185 // phase properties are calculated
186 const ADReal Tsatl = Tsat - _dT;
187 const ADReal Tsatv = Tsat + _dT;
188
189 // Density
190 const ADReal rhol = _water_fp.rho_from_p_T(pressure, Tsatl);
191 const ADReal rhov = _water_fp.rho_from_p_T(pressure, Tsatv);
192
193 // Vapor (gas) saturation
194 const ADReal satv = X * rhol / (rhov + X * (rhol - rhov));
195
196 gas.temperature = Tsat;
197 gas.density = rhov;
198 gas.viscosity = _water_fp.mu_from_p_T(pressure, Tsatv);
199 gas.enthalpy = hv;
200 gas.internal_energy = _water_fp.e_from_p_T(pressure, Tsatv);
201 gas.saturation = satv;
202
203 liquid.temperature = Tsat;
204 liquid.density = rhol;
205 liquid.viscosity = _water_fp.mu_from_p_T(pressure, Tsatl);
206 liquid.enthalpy = hl;
207 liquid.internal_energy = _water_fp.e_from_p_T(pressure, Tsatl);
208 liquid.saturation = 1.0 - satv;
209
210 liquid.pressure = pressure;
211 gas.pressure = pressure + _pc.capillaryPressure(liquid.saturation, qp);
212
213 gas.mass_fraction[_fluid_component] = 1.0;
214 liquid.mass_fraction[_fluid_component] = 1.0;
215
216 break;
217 }
218 }
219}
ADReal
DualNumber< Real, DNDerivativeType, true > ADReal
T
const double T
Definition GasLiquidMassTransferTest.C:19
PorousFlowCapillaryPressure::capillaryPressure
virtual Real capillaryPressure(Real saturation, unsigned qp=0) const
Capillary pressure is calculated as a function of true saturation.
Definition PorousFlowCapillaryPressure.C:63
PorousFlowFluidStateBase::_pc
const PorousFlowCapillaryPressure & _pc
Capillary pressure UserObject.
Definition PorousFlowFluidStateBase.h:112
PorousFlowFluidStateSingleComponentBase::_dT
const Real _dT
Perturbation applied to saturation temperature to move to gas/liquid phase.
Definition PorousFlowFluidStateSingleComponentBase.h:67
SinglePhaseFluidProperties::vaporTemperature
virtual Real vaporTemperature(Real p) const
Vapor temperature.
Definition SinglePhaseFluidProperties.C:408
FluidStateProperties::mass_fraction
std::vector< ADReal > mass_fraction
Definition PorousFlowFluidStateBase.h:46

◆ thermophysicalProperties() [2/3]

void PorousFlowWaterVapor::thermophysicalProperties ( const ADReal &  pressure,
const ADReal &  enthalpy,
unsigned int  qp,
std::vector< FluidStateProperties > &  fsp 
) const
overridevirtual

Implements PorousFlowFluidStateSingleComponentBase.

Definition at line 81 of file PorousFlowWaterVapor.C.

85{
86 FluidStatePhaseEnum phase_state;
87
88 thermophysicalProperties(pressure, enthalpy, qp, phase_state, fsp);
89}
FluidStatePhaseEnum
FluidStatePhaseEnum
Phase state enum.
Definition PorousFlowFluidStateBase.h:17
PorousFlowWaterVapor::thermophysicalProperties
void thermophysicalProperties(Real pressure, Real enthalpy, unsigned int qp, std::vector< FluidStateProperties > &fsp) const override
Determines the complete thermophysical state of the system for a given set of primary variables.
Definition PorousFlowWaterVapor.C:66

◆ thermophysicalProperties() [3/3]

void PorousFlowWaterVapor::thermophysicalProperties ( Real  pressure,
Real  enthalpy,
unsigned int  qp,
std::vector< FluidStateProperties > &  fsp 
) const
overridevirtual

Determines the complete thermophysical state of the system for a given set of primary variables.

Parameters
pressuregas phase pressure (Pa)
enthalpyfluid enthalpy (J/kg)
qpquadpoint index
[out]fspthe FluidStateProperties struct containing all properties

Implements PorousFlowFluidStateSingleComponentBase.

Definition at line 66 of file PorousFlowWaterVapor.C.

70{
71 // AD versions of primary variables
72 ADReal p = pressure;
73 Moose::derivInsert(p.derivatives(), _pidx, 1.0);
74 ADReal h = enthalpy;
75 Moose::derivInsert(h.derivatives(), _hidx, 1.0);
76
77 thermophysicalProperties(p, h, qp, fsp);
78}
p
const Real p
Definition GasLiquidMassTransferTest.C:18
Moose::derivInsert
void derivInsert(SemiDynamicSparseNumberArray< Real, libMesh::dof_id_type, NWrapper< N > > &derivs, libMesh::dof_id_type index, Real value)

Referenced by thermophysicalProperties(), and thermophysicalProperties().

◆ validParams()

InputParameters PorousFlowWaterVapor::validParams ( )
static

Definition at line 18 of file PorousFlowWaterVapor.C.

19{
20 InputParameters params = PorousFlowFluidStateSingleComponentBase::validParams();
21 params.addRequiredParam<UserObjectName>("water_fp", "The name of the user object for water");
22 params.addClassDescription("Fluid state class for water and vapor");
23 return params;
24}
InputParameters::addRequiredParam
void addRequiredParam(const std::string &name, const std::string &doc_string)
InputParameters::addClassDescription
void addClassDescription(const std::string &doc_string)

Member Data Documentation

◆ _aqueous_phase_number

const unsigned int PorousFlowFluidStateBase::_aqueous_phase_number
protectedinherited

Phase number of the aqueous phase.

Definition at line 104 of file PorousFlowFluidStateBase.h.

Referenced by PorousFlowFluidStateBase::aqueousPhaseIndex(), PorousFlowWaterNCG::gasDensity(), PorousFlowWaterNCG::gasProperties(), PorousFlowBrineCO2::liquidProperties(), PorousFlowWaterNCG::liquidProperties(), PorousFlowBrineCO2::massFractions(), PorousFlowWaterNCG::massFractions(), PorousFlowBrineCO2::PorousFlowBrineCO2(), PorousFlowWaterNCG::PorousFlowWaterNCG(), PorousFlowWaterVapor(), PorousFlowBrineCO2::saturation(), PorousFlowWaterNCG::saturation(), thermophysicalProperties(), PorousFlowBrineCO2::thermophysicalProperties(), PorousFlowWaterNCG::thermophysicalProperties(), PorousFlowBrineCO2::totalMassFraction(), and PorousFlowWaterNCG::totalMassFraction().

◆ _dT

const Real PorousFlowFluidStateSingleComponentBase::_dT
protectedinherited

Perturbation applied to saturation temperature to move to gas/liquid phase.

Definition at line 67 of file PorousFlowFluidStateSingleComponentBase.h.

Referenced by thermophysicalProperties().

◆ _empty_fsp

FluidStateProperties PorousFlowFluidStateBase::_empty_fsp
protectedinherited

Empty FluidStateProperties object.

Definition at line 114 of file PorousFlowFluidStateBase.h.

Referenced by PorousFlowFluidStateBase::clearFluidStateProperties(), PorousFlowBrineCO2::PorousFlowBrineCO2(), PorousFlowWaterNCG::PorousFlowWaterNCG(), and PorousFlowWaterVapor().

◆ _fluid_component

const unsigned int PorousFlowFluidStateSingleComponentBase::_fluid_component
protectedinherited

Fluid component number (only one fluid component in all phases)

Definition at line 61 of file PorousFlowFluidStateSingleComponentBase.h.

Referenced by PorousFlowFluidStateSingleComponentBase::aqueousComponentIndex(), PorousFlowFluidStateSingleComponentBase::gasComponentIndex(), PorousFlowWaterVapor(), and thermophysicalProperties().

◆ _gas_phase_number

unsigned int PorousFlowFluidStateBase::_gas_phase_number
protectedinherited

Phase number of the gas phase.

Definition at line 106 of file PorousFlowFluidStateBase.h.

Referenced by PorousFlowWaterNCG::gasDensity(), PorousFlowFluidStateBase::gasPhaseIndex(), PorousFlowBrineCO2::gasProperties(), PorousFlowWaterNCG::gasProperties(), PorousFlowBrineCO2::massFractions(), PorousFlowWaterNCG::massFractions(), PorousFlowBrineCO2::PorousFlowBrineCO2(), PorousFlowWaterNCG::PorousFlowWaterNCG(), PorousFlowWaterVapor(), PorousFlowBrineCO2::saturation(), PorousFlowWaterNCG::saturation(), thermophysicalProperties(), PorousFlowBrineCO2::thermophysicalProperties(), PorousFlowWaterNCG::thermophysicalProperties(), PorousFlowBrineCO2::totalMassFraction(), PorousFlowWaterNCG::totalMassFraction(), PorousFlowBrineCO2::twoPhaseProperties(), and PorousFlowWaterNCG::twoPhaseProperties().

◆ _hidx

const unsigned int PorousFlowFluidStateSingleComponentBase::_hidx
protectedinherited

Index of derivative wrt enthalpy.

Definition at line 65 of file PorousFlowFluidStateSingleComponentBase.h.

Referenced by PorousFlowFluidStateSingleComponentBase::getEnthalpyIndex(), and thermophysicalProperties().

◆ _Mh2o

const Real PorousFlowWaterVapor::_Mh2o
protected

Molar mass of water (kg/mol)

Definition at line 51 of file PorousFlowWaterVapor.h.

◆ _num_components

unsigned int PorousFlowFluidStateBase::_num_components
protectedinherited

Number of components.

Definition at line 102 of file PorousFlowFluidStateBase.h.

Referenced by PorousFlowFluidStateBase::numComponents(), PorousFlowBrineCO2::PorousFlowBrineCO2(), PorousFlowWaterNCG::PorousFlowWaterNCG(), PorousFlowWaterVapor(), PorousFlowBrineCO2::totalMassFraction(), and PorousFlowWaterNCG::totalMassFraction().

◆ _num_phases

unsigned int PorousFlowFluidStateBase::_num_phases
protectedinherited

Number of phases.

Definition at line 100 of file PorousFlowFluidStateBase.h.

Referenced by PorousFlowFluidStateBase::numPhases(), PorousFlowBrineCO2::PorousFlowBrineCO2(), PorousFlowWaterNCG::PorousFlowWaterNCG(), PorousFlowWaterVapor(), PorousFlowBrineCO2::totalMassFraction(), and PorousFlowWaterNCG::totalMassFraction().

◆ _p_critical

const Real PorousFlowWaterVapor::_p_critical
protected

Critical pressure of water (Pa)

Definition at line 55 of file PorousFlowWaterVapor.h.

Referenced by thermophysicalProperties().

◆ _p_triple

const Real PorousFlowWaterVapor::_p_triple
protected

Triple point pressure of water (Pa)

Definition at line 53 of file PorousFlowWaterVapor.h.

Referenced by thermophysicalProperties().

◆ _pc

const PorousFlowCapillaryPressure& PorousFlowFluidStateBase::_pc
protectedinherited

Capillary pressure UserObject.

Definition at line 112 of file PorousFlowFluidStateBase.h.

Referenced by thermophysicalProperties(), PorousFlowBrineCO2::thermophysicalProperties(), PorousFlowWaterNCG::thermophysicalProperties(), PorousFlowBrineCO2::totalMassFraction(), PorousFlowWaterNCG::totalMassFraction(), PorousFlowBrineCO2::twoPhaseProperties(), and PorousFlowWaterNCG::twoPhaseProperties().

◆ _pidx

const unsigned int PorousFlowFluidStateSingleComponentBase::_pidx
protectedinherited

Index of derivative wrt pressure.

Definition at line 63 of file PorousFlowFluidStateSingleComponentBase.h.

Referenced by PorousFlowFluidStateSingleComponentBase::getPressureIndex(), and thermophysicalProperties().

◆ _R

const Real PorousFlowFluidStateBase::_R
protectedinherited

Universal gas constant (J/mol/K)

Definition at line 108 of file PorousFlowFluidStateBase.h.

Referenced by PorousFlowWaterNCG::enthalpyOfDissolution(), and PorousFlowBrineCO2::enthalpyOfDissolutionGas().

◆ _T_c2k

const Real PorousFlowFluidStateBase::_T_c2k
protectedinherited

Conversion from C to K.

Definition at line 110 of file PorousFlowFluidStateBase.h.

Referenced by PorousFlowBrineCO2::equilibriumConstantCO2(), PorousFlowBrineCO2::equilibriumConstantH2O(), PorousFlowBrineCO2::funcABHighTemp(), PorousFlowBrineCO2::henryConstant(), and PorousFlowBrineCO2::partialDensityCO2().

◆ _T_critical

const Real PorousFlowWaterVapor::_T_critical
protected

Critical temperature of water (K)

Definition at line 59 of file PorousFlowWaterVapor.h.

Referenced by thermophysicalProperties().

◆ _T_triple

const Real PorousFlowWaterVapor::_T_triple
protected

Triple point temperature of water (K)

Definition at line 57 of file PorousFlowWaterVapor.h.

◆ _water_fp

const SinglePhaseFluidProperties& PorousFlowWaterVapor::_water_fp
protected

Fluid properties UserObject for water.

Definition at line 49 of file PorousFlowWaterVapor.h.

Referenced by PorousFlowWaterVapor(), and thermophysicalProperties().

The documentation for this class was generated from the following files:
Generated on Fri Aug 21 2026 13:43:40 for https://mooseframework.inl.gov by doxygen 1.9.8