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
BrineFluidProperties Class Reference

Brine (NaCl in H2O) fluid properties as a function of pressure (Pa), temperature (K) and NaCl mass fraction. More...

#include <BrineFluidProperties.h>

Inheritance diagram for BrineFluidProperties:
[legend]

Public Types

typedef DataFileName DataFileParameterType
 

Public Member Functions

 BrineFluidProperties (const InputParameters &parameters)
 
virtual ~BrineFluidProperties ()
 
virtual std::string fluidName () const override
 Fluid name.
 
Real molarMass (Real xnacl) const
 Average molar mass of brine.
 
FPADReal molarMass (const FPADReal &xnacl) const
 
Real molarMassNaCl () const
 NaCl molar mass.
 
Real molarMassH2O () const
 H2O molar mass.
 
virtual Real rho_from_p_T_X (Real pressure, Real temperature, Real xnacl) const override
 
FPADReal rho_from_p_T_X (const FPADReal &pressure, const FPADReal &temperature, const FPADReal &xnacl) const
 
virtual void rho_from_p_T_X (Real pressure, Real temperature, Real xnacl, Real &rho, Real &drho_dp, Real &drho_dT, Real &drho_dx) const override
 
virtual Real mu_from_p_T_X (Real pressure, Real temperature, Real xnacl) const override
 
virtual void mu_from_p_T_X (Real pressure, Real temperature, Real xnacl, Real &mu, Real &dmu_dp, Real &dmu_dT, Real &dmu_dx) const override
 
FPADReal h_from_p_T_X (const FPADReal &pressure, const FPADReal &temperature, const FPADReal &xnacl) const
 
virtual Real h_from_p_T_X (Real pressure, Real temperature, Real xnacl) const override
 
virtual void h_from_p_T_X (Real pressure, Real temperature, Real xnacl, Real &h, Real &dh_dp, Real &dh_dT, Real &dh_dx) const override
 
virtual Real cp_from_p_T_X (Real pressure, Real temperature, Real xnacl) const override
 
FPADReal e_from_p_T_X (const FPADReal &pressure, const FPADReal &temperature, const FPADReal &xnacl) const
 
virtual Real e_from_p_T_X (Real pressure, Real temperature, Real xnacl) const override
 
virtual void e_from_p_T_X (Real pressure, Real temperature, Real xnacl, Real &e, Real &de_dp, Real &de_dT, Real &de_dx) const override
 
virtual Real k_from_p_T_X (Real pressure, Real temperature, Real xnacl) const override
 
Real vaporPressure (Real temperature, Real xnacl) const
 Brine vapour pressure From Haas, Physical properties of the coexisting phases and thermochemical properties of the H2O component in boiling NaCl solutions, Geological Survey Bulletin, 1421-A (1976).
 
Real haliteSolubility (Real temperature) const
 Solubility of halite (solid NaCl) in water Originally from Potter et al., A new method for determining the solubility of salts in aqueous solutions at elevated temperatures, J.
 
Real henryConstant (Real temperature, const std::vector< Real > &coeffs) const
 IAPWS formulation of Henry's law constant for dissolution in water (implemented in water FluidProperties userobject)
 
void henryConstant (Real temperature, const std::vector< Real > &coeffs, Real &Kh, Real &dKh_dT) const
 
ADReal henryConstant (const ADReal &temperature, const std::vector< Real > &coeffs) const
 
virtual const SinglePhaseFluidProperties & getComponent (unsigned int component) const override
 Get UserObject for specified component.
 
virtual void rho_mu_from_p_T_X (Real pressure, Real temperature, Real xmass, Real &rho, Real &mu) const
 Density and viscosity.
 
virtual void rho_mu_from_p_T_X (ADReal pressure, ADReal temperature, ADReal xmass, ADReal &rho, ADReal &mu) const
 
virtual void rho_mu_from_p_T_X (Real pressure, Real temperature, Real xmass, Real &rho, Real &drho_dp, Real &drho_dT, Real &drho_dx, Real &mu, Real &dmu_dp, Real &dmu_dT, Real &dmu_dx) const
 Density and viscosity and their derivatives wrt pressure, temperature and mass fraction.
 
virtual void execute () final
 
virtual void initialize () final
 
virtual void finalize () final
 
virtual void threadJoin (const UserObject &) final
 
virtual void subdomainSetup () final
 
bool needThreadedCopy () const override final
 
virtual Real spatialValue (const Point &) const
 
virtual const std::vector< Point > spatialPoints () const
 
void setPrimaryThreadCopy (UserObject *primary)
 
UserObject * primaryThreadCopy ()
 
SubProblem & getSubProblem () const
 
bool shouldDuplicateInitialExecution () const
 
void gatherSum (T &value)
 
void gatherMax (T &value)
 
void gatherMin (T &value)
 
void gatherProxyValueMax (T1 &proxy, T2 &value)
 
void gatherProxyValueMin (T1 &proxy, T2 &value)
 
std::set< UserObjectName > getDependObjects () const
 
const std::set< std::string > & getRequestedItems () override
 
const std::set< std::string > & getSuppliedItems () override
 
unsigned int systemNumber () const
 
virtual bool enabled () const
 
std::shared_ptr< MooseObject > getSharedPtr ()
 
std::shared_ptr< const MooseObject > getSharedPtr () const
 
bool isKokkosObject () const
 
MooseApp & getMooseApp () const
 
const std::string & type () const
 
const std::string & name () const
 
std::string typeAndName () const
 
MooseObjectParameterName uniqueParameterName (const std::string &parameter_name) const
 
MooseObjectName uniqueName () const
 
const InputParameters & parameters () const
 
const hit::Node * getHitNode () const
 
bool hasBase () const
 
const std::string & getBase () const
 
const T & getParam (const std::string &name) const
 
std::vector< std::pair< T1, T2 > > getParam (const std::string &param1, const std::string &param2) const
 
const T * queryParam (const std::string &name) const
 
const T & getRenamedParam (const std::string &old_name, const std::string &new_name) const
 
T getCheckedPointerParam (const std::string &name, const std::string &error_string="") const
 
bool isParamValid (const std::string &name) const
 
bool isParamSetByUser (const std::string &name) const
 
void connectControllableParams (const std::string &parameter, const std::string &object_type, const std::string &object_name, const std::string &object_parameter) const
 
void paramError (const std::string &param, Args... args) const
 
void paramWarning (const std::string &param, Args... args) const
 
void paramWarning (const std::string &param, Args... args) const
 
void paramInfo (const std::string &param, Args... args) const
 
std::string messagePrefix (const bool hit_prefix=true) const
 
std::string errorPrefix (const std::string &) const
 
void mooseError (Args &&... args) const
 
void mooseDocumentedError (const std::string &repo_name, const unsigned int issue_num, Args &&... args) const
 
void mooseErrorNonPrefixed (Args &&... args) const
 
void mooseWarning (Args &&... args) const
 
void mooseWarning (Args &&... args) const
 
void mooseWarningNonPrefixed (Args &&... args) const
 
void mooseWarningNonPrefixed (Args &&... args) const
 
void mooseDeprecated (Args &&... args) const
 
void mooseDeprecated (Args &&... args) const
 
void mooseDeprecatedNoTrace (Args &&... args) const
 
void mooseInfo (Args &&... args) const
 
void callMooseError (std::string msg, const bool with_prefix, const hit::Node *node=nullptr, const bool show_trace=true) const
 
std::string getDataFileName (const std::string &param) const
 
std::string getDataFileNameByName (const std::string &relative_path) const
 
std::string getDataFilePath (const std::string &relative_path) const
 
virtual void 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
 
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 unsigned int WATER = 0
 Fluid component numbers for water and NaCl.
 
static const unsigned int NACL = 1
 
static const Real _R = 8.3144598
 Universal gas constant (J/mol/K)
 
static const std::string type_param
 
static const std::string name_param
 
static const std::string unique_name_param
 
static const std::string app_param
 
static const std::string moose_base_param
 
static const std::string kokkos_object_param
 
static constexpr PropertyValue::id_type default_property_id
 
static constexpr PropertyValue::id_type zero_property_id
 
static constexpr auto SYSTEM
 
static constexpr auto NAME
 

Protected Member Functions

Real massFractionToMolalConc (Real xnacl) const
 Conversion from mass fraction to molal concentration (molality)
 
Real massFractionToMoleFraction (Real xnacl) const
 Conversion from mass fraction to mole fraction.
 
FPADReal massFractionToMoleFraction (const FPADReal &xnacl) const
 
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 Water97FluidProperties * _water97_fp
 Water97FluidProperties UserObject (for Henry's law)
 
const SinglePhaseFluidProperties * _water_fp
 Water97FluidProperties UserObject.
 
const SinglePhaseFluidProperties * _nacl_fp
 NaClFluidProperties UserObject.
 
Real _Mnacl
 Molar mass of NaCl (kg/mol)
 
Real _Mh2o
 Molar mass of water (H2O) (kg/mol)
 
bool _water_fp_derivs
 Flag to indicate whether to calculate derivatives in water_fp.
 
const Real _T_c2k
 Conversion of temperature from Celsius to Kelvin.
 
const bool _allow_imperfect_jacobians
 Flag to set unimplemented Jacobian entries to zero.
 
const Moose::CoordinateSystemType & _coord_sys
 
const THREAD_ID _tid
 
SubProblem & _subproblem
 
FEProblemBase & _fe_problem
 
SystemBase & _sys
 
Assembly & _assembly
 
const bool _duplicate_initial_execution
 
std::set< std::string > _depend_uo
 
const bool & _enabled
 
MooseApp & _app
 
Factory & _factory
 
ActionFactory & _action_factory
 
const std::string & _type
 
const std::string & _name
 
const InputParameters & _pars
 
const ExecFlagEnum & _execute_enum
 
const ExecFlagType & _current_execute_flag
 
MooseApp & _restartable_app
 
const std::string _restartable_system_name
 
const THREAD_ID _restartable_tid
 
const bool _restartable_read_only
 
FEProblemBase & _mci_feproblem
 
FEProblemBase & _mdi_feproblem
 
MooseApp & _pg_moose_app
 
const std::string _prefix
 
FEProblemBase & _sc_fe_problem
 
const THREAD_ID _sc_tid
 
const Real & _real_zero
 
const VariableValue & _scalar_zero
 
const Point & _point_zero
 
const InputParameters & _mi_params
 
const std::string _mi_name
 
const MooseObjectName _mi_moose_object_name
 
FEProblemBase & _mi_feproblem
 
SubProblem & _mi_subproblem
 
const THREAD_ID _mi_tid
 
const bool _is_kokkos_object
 
const Moose::MaterialDataType _material_data_type
 
MaterialData & _material_data
 
bool _stateful_allowed
 
bool _get_material_property_called
 
std::vector< std::unique_ptr< PropertyValue > > _default_properties
 
std::unordered_set< unsigned int > _material_property_dependencies
 
const MaterialPropertyName _get_suffix
 
const bool _use_interpolated_state
 
const InputParameters & _ti_params
 
FEProblemBase & _ti_feproblem
 
bool _is_implicit
 
Real & _t
 
const Real & _t_old
 
int & _t_step
 
Real & _dt
 
Real & _dt_old
 
bool _is_transient
 
const Parallel::Communicator & _communicator
 

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
 
 propfunc (rho, p, T, X) propfunc(mu
 Compute a fluid property given for the state defined by three given properties.
 
X propfunc (h, p, T, X) propfunc(cp
 
X X propfunc (e, p, T, X) propfunc(k
 
 p
 
X p
 
X X p
 
 T
 
X T
 
X X T
 

Detailed Description

Brine (NaCl in H2O) fluid properties as a function of pressure (Pa), temperature (K) and NaCl mass fraction.

Most properties from: Driesner, The system H2O-NaCl. Part II: Correlations for molar volume, enthalpy, and isobaric heat capacity from 0 to 1000 C, 1 to 5000 bar, and 0 to 1 Xnacl, Geochimica et Cosmochimica Acta 71, 4902-4919 (2007)

Viscosity and thermal conductivity from: Phillips et al, A technical databook for geothermal energy utilization, LbL-12810 (1981) Note: uses water thermal conductivity from IAPWS rather than the correlation given by Phillips et al.

Definition at line 33 of file BrineFluidProperties.h.

Constructor & Destructor Documentation

◆ BrineFluidProperties()

BrineFluidProperties::BrineFluidProperties ( const InputParameters &  parameters)

Definition at line 24 of file BrineFluidProperties.C.

25 : MultiComponentFluidProperties(parameters), _water_fp_derivs(true)
26{
27 // There are two possibilities to consider:
28 // 1) No water_fp has been supplied (in which case one is constructed)
29 // 2) A water_fp hase been supplied (in which case it is used)
30 // In both cases, though, a Water97FluidProperties UserObject must be added
31 // Note: this UserObject is only used to gain access to the Henry's constant
32 // formulation. All property calculations are performed using _water_fp
33 const std::string water_name = name() + ":water";
34 {
35 const std::string class_name = "Water97FluidProperties";
36 InputParameters params = _app.getFactory().getValidParams(class_name);
37 if (_tid == 0)
38 _fe_problem.addUserObject(class_name, water_name, params);
39 }
40 _water97_fp = &_fe_problem.getUserObject<Water97FluidProperties>(water_name);
41
42 if (parameters.isParamSetByUser("water_fp"))
43 {
44 // SinglePhaseFluidPropertiesPT UserObject for water
45 _water_fp = &getUserObject<SinglePhaseFluidProperties>("water_fp");
46
47 // Check that a water userobject has actually been supplied
48 if (_water_fp->fluidName() != "water")
49 paramError("water_fp", "A water FluidProperties UserObject must be supplied");
50 }
51 else
52 {
53 // Construct a SinglePhaseFluidProperties UserObject for water
54 _water_fp = &_fe_problem.getUserObject<SinglePhaseFluidProperties>(water_name);
55 }
56
57 // SinglePhaseFluidProperties UserObject for NaCl
58 const std::string nacl_name = name() + ":nacl";
59 {
60 const std::string class_name = "NaClFluidProperties";
61 InputParameters params = _app.getFactory().getValidParams(class_name);
62 if (_tid == 0)
63 _fe_problem.addUserObject(class_name, nacl_name, params);
64 }
65 _nacl_fp = &_fe_problem.getUserObject<SinglePhaseFluidProperties>(nacl_name);
66
67 // Molar mass of NaCl and H20
68 _Mnacl = _nacl_fp->molarMass();
69 _Mh2o = _water_fp->molarMass();
70}
BrineFluidProperties::_nacl_fp
const SinglePhaseFluidProperties * _nacl_fp
NaClFluidProperties UserObject.
Definition BrineFluidProperties.h:194
BrineFluidProperties::_Mh2o
Real _Mh2o
Molar mass of water (H2O) (kg/mol)
Definition BrineFluidProperties.h:199
BrineFluidProperties::_water_fp
const SinglePhaseFluidProperties * _water_fp
Water97FluidProperties UserObject.
Definition BrineFluidProperties.h:192
BrineFluidProperties::_water_fp_derivs
bool _water_fp_derivs
Flag to indicate whether to calculate derivatives in water_fp.
Definition BrineFluidProperties.h:201
BrineFluidProperties::_water97_fp
const Water97FluidProperties * _water97_fp
Water97FluidProperties UserObject (for Henry's law)
Definition BrineFluidProperties.h:190
BrineFluidProperties::_Mnacl
Real _Mnacl
Molar mass of NaCl (kg/mol)
Definition BrineFluidProperties.h:197
FEProblemBase::getUserObject
T & getUserObject(const std::string &name, unsigned int tid=0) const
FEProblemBase::addUserObject
virtual std::vector< std::shared_ptr< UserObject > > addUserObject(const std::string &user_object_name, const std::string &name, InputParameters &parameters)
Factory::getValidParams
InputParameters getValidParams(const std::string &name) const
InputParameters::isParamSetByUser
bool isParamSetByUser(const std::string &name) const
MooseApp::getFactory
Factory & getFactory()
ThreadedGeneralUserObject::parameters
const InputParameters & parameters() const
ThreadedGeneralUserObject::name
const std::string & name() const
ThreadedGeneralUserObject::paramError
void paramError(const std::string &param, Args... args) const
MultiComponentFluidProperties
Common class for multiple component fluid properties using a pressure and temperature formulation.
Definition MultiComponentFluidProperties.h:65
SinglePhaseFluidProperties
Common class for single phase fluid properties.
Definition SinglePhaseFluidProperties.h:115
SinglePhaseFluidProperties::molarMass
virtual Real molarMass() const
Molar mass [kg/mol].
Definition SinglePhaseFluidProperties.C:290
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
ThreadedGeneralUserObject::_fe_problem
FEProblemBase & _fe_problem
Water97FluidProperties
Water (H2O) fluid properties as a function of pressure (Pa) and temperature (K) from IAPWS-IF97: Revi...
Definition Water97FluidProperties.h:39

◆ ~BrineFluidProperties()

BrineFluidProperties::~BrineFluidProperties ( )
virtual

Definition at line 72 of file BrineFluidProperties.C.

72{}

Member Function Documentation

◆ cp_from_p_T_X()

Real BrineFluidProperties::cp_from_p_T_X ( Real  pressure,
Real  temperature,
Real  xnacl 
) const
overridevirtual

Definition at line 354 of file BrineFluidProperties.C.

355{
356 Real q1, q2, q10, q11, q12, q20, q21, q22, q23, q1x1, q2x1, Th;
357
358 // The correlation requires the pressure in bar, not Pa.
359 Real pbar = pressure * 1.0e-5;
360 Real pbar2 = pbar * pbar;
361
362 // The correlation requires mole fraction
363 Real Xnacl = massFractionToMoleFraction(xnacl);
364
365 q11 = -32.1724 + 0.0621255 * pbar;
366 q21 = -1.69513 - 4.52781e-4 * pbar - 6.04279e-8 * pbar2;
367 q22 = 0.0612567 + 1.88082e-5 * pbar;
368
369 q1x1 = 47.9048 - 9.36994e-3 * pbar + 6.51059e-6 * pbar2;
370 q2x1 = 0.241022 + 3.45087e-5 * pbar - 4.28356e-9 * pbar2;
371
372 q12 = -q11 - q1x1;
373 q10 = q1x1;
374
375 q20 = 1.0 - q21 * std::sqrt(q22);
376 q23 = q2x1 - q20 - q21 * std::sqrt(1.0 + q22);
377
378 q1 = q10 + q11 * (1.0 - Xnacl) + q12 * (1.0 - Xnacl) * (1.0 - Xnacl);
379 q2 = q20 + q21 * std::sqrt(Xnacl + q22) + q23 * Xnacl;
380 // The temperature Th where the brine has the same isobaric heat capacity
381 // as pure water. Note: correlation uses temperature in Celcius
382 Th = q1 + q2 * (temperature - _T_c2k);
383
384 // The brine isobaric heat capacity is then given by the isobaric heat
385 // capacity of water at temperature Th multiplied by q2
386 // Note: water isobaric heat capacity requires temperature in Kelvin
387 return q2 * _water_fp->cp_from_p_T(pressure, Th + _T_c2k);
388}
BrineFluidProperties::massFractionToMoleFraction
Real massFractionToMoleFraction(Real xnacl) const
Conversion from mass fraction to mole fraction.
Definition BrineFluidProperties.C:488
FluidProperties::_T_c2k
const Real _T_c2k
Conversion of temperature from Celsius to Kelvin.
Definition FluidProperties.h:45
libMesh::Real
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real

◆ e_from_p_T_X() [1/3]

FPADReal BrineFluidProperties::e_from_p_T_X ( const FPADReal &  pressure,
const FPADReal &  temperature,
const FPADReal &  xnacl 
) const

Definition at line 391 of file BrineFluidProperties.C.

394{
395 FPADReal enthalpy = h_from_p_T_X(pressure, temperature, xnacl);
396 FPADReal density = rho_from_p_T_X(pressure, temperature, xnacl);
397
398 return enthalpy - pressure / density;
399}
FPADReal
DualNumber< Real, DNDerivativeSize< 5 > > FPADReal
Definition FluidProperties.h:23
BrineFluidProperties::rho_from_p_T_X
virtual Real rho_from_p_T_X(Real pressure, Real temperature, Real xnacl) const override
Definition BrineFluidProperties.C:171
BrineFluidProperties::h_from_p_T_X
FPADReal h_from_p_T_X(const FPADReal &pressure, const FPADReal &temperature, const FPADReal &xnacl) const
Definition BrineFluidProperties.C:264

Referenced by PorousFlowBrine::computeQpProperties(), e_from_p_T_X(), and PorousFlowBrine::initQpStatefulProperties().

◆ e_from_p_T_X() [2/3]

Real BrineFluidProperties::e_from_p_T_X ( Real  pressure,
Real  temperature,
Real  xnacl 
) const
overridevirtual

Definition at line 402 of file BrineFluidProperties.C.

403{
404 Real enthalpy = h_from_p_T_X(pressure, temperature, xnacl);
405 Real density = rho_from_p_T_X(pressure, temperature, xnacl);
406
407 return enthalpy - pressure / density;
408}

◆ e_from_p_T_X() [3/3]

void BrineFluidProperties::e_from_p_T_X ( Real  pressure,
Real  temperature,
Real  xnacl,
Real &  e,
Real &  de_dp,
Real &  de_dT,
Real &  de_dx 
) const
overridevirtual

Definition at line 411 of file BrineFluidProperties.C.

418{
419 // Initialise the AD value and derivatives
420 FPADReal p = pressure;
421 Moose::derivInsert(p.derivatives(), 0, 1.0);
422 FPADReal T = temperature;
423 Moose::derivInsert(T.derivatives(), 1, 1.0);
424 FPADReal x = xnacl;
425 Moose::derivInsert(x.derivatives(), 2, 1.0);
426
427 _water_fp_derivs = true;
428 FPADReal ad_e = e_from_p_T_X(p, T, x);
429
430 e = ad_e.value();
431 de_dp = ad_e.derivatives()[0];
432 de_dT = ad_e.derivatives()[1];
433 de_dx = ad_e.derivatives()[2];
434}
x
const std::vector< double > x
Definition EquilibriumConstantFitTest.C:17
BrineFluidProperties::e_from_p_T_X
FPADReal e_from_p_T_X(const FPADReal &pressure, const FPADReal &temperature, const FPADReal &xnacl) const
Definition BrineFluidProperties.C:391
Moose::derivInsert
void derivInsert(SemiDynamicSparseNumberArray< Real, libMesh::dof_id_type, NWrapper< N > > &derivs, libMesh::dof_id_type index, Real value)

◆ execute()

virtual void FluidProperties::execute ( )
inlinefinalvirtualinherited

Implements ThreadedGeneralUserObject.

Definition at line 33 of file FluidProperties.h.

33{}

◆ finalize()

virtual void FluidProperties::finalize ( )
inlinefinalvirtualinherited

Implements ThreadedGeneralUserObject.

Definition at line 35 of file FluidProperties.h.

35{}

◆ fluidName()

std::string BrineFluidProperties::fluidName ( ) const
overridevirtual

Fluid name.

Returns
"brine"

Reimplemented from MultiComponentFluidProperties.

Definition at line 91 of file BrineFluidProperties.C.

92{
93 return "brine";
94}

Referenced by PorousFlowBrineCO2::PorousFlowBrineCO2().

◆ getComponent()

const SinglePhaseFluidProperties & BrineFluidProperties::getComponent ( unsigned int  component) const
overridevirtual

Get UserObject for specified component.

Parameters
componentfluid component
Returns
reference to SinglePhaseFluidPropertiesPT UserObject for component

Reimplemented from MultiComponentFluidProperties.

Definition at line 75 of file BrineFluidProperties.C.

76{
77 switch (component)
78 {
79 case WATER:
80 return *_water_fp;
81
82 case NACL:
83 return *_nacl_fp;
84
85 default:
86 mooseError("BrineFluidProperties::getComponent has been provided an incorrect component");
87 }
88}
BrineFluidProperties::NACL
static const unsigned int NACL
Definition BrineFluidProperties.h:169
BrineFluidProperties::WATER
static const unsigned int WATER
Fluid component numbers for water and NaCl.
Definition BrineFluidProperties.h:168
ThreadedGeneralUserObject::mooseError
void mooseError(Args &&... args) const

Referenced by BrineFluidPropertiesTest::buildObjects().

◆ h_from_p_T_X() [1/3]

FPADReal BrineFluidProperties::h_from_p_T_X ( const FPADReal &  pressure,
const FPADReal &  temperature,
const FPADReal &  xnacl 
) const

Definition at line 264 of file BrineFluidProperties.C.

267{
268 using std::sqrt;
269
270 FPADReal q1, q2, q10, q11, q12, q20, q21, q22, q23, q1x1, q2x1, Th;
271
272 // The correlation requires the pressure in bar, not Pa.
273 const FPADReal pbar = pressure * 1.0e-5;
274 const FPADReal pbar2 = pbar * pbar;
275
276 // The correlation requires mole fraction
277 const FPADReal Xnacl = massFractionToMoleFraction(xnacl);
278
279 q11 = -32.1724 + 0.0621255 * pbar;
280 q21 = -1.69513 - 4.52781e-4 * pbar - 6.04279e-8 * pbar2;
281 q22 = 0.0612567 + 1.88082e-5 * pbar;
282
283 q1x1 = 47.9048 - 9.36994e-3 * pbar + 6.51059e-6 * pbar2;
284 q2x1 = 0.241022 + 3.45087e-5 * pbar - 4.28356e-9 * pbar2;
285
286 q12 = -q11 - q1x1;
287 q10 = q1x1;
288
289 q20 = 1.0 - q21 * sqrt(q22);
290 q23 = q2x1 - q20 - q21 * sqrt(1.0 + q22);
291
292 q1 = q10 + q11 * (1.0 - Xnacl) + q12 * (1.0 - Xnacl) * (1.0 - Xnacl);
293 q2 = q20 + q21 * sqrt(Xnacl + q22) + q23 * Xnacl;
294 // The temperature Th where the brine has the same enthalpy as pure water
295 // Note: correlation uses temperature in Celcius
296 Th = q1 + q2 * (temperature - _T_c2k);
297
298 // The brine enthalpy is then given by the enthalpy of water at temperature Th
299 // Note: water enthalpy requires temperature in Kelvin
300 FPADReal enthalpy;
301 if (_water_fp_derivs)
302 {
303 Real h, dh_dp, dh_dT;
304 _water_fp->h_from_p_T(pressure.value(), Th.value() + _T_c2k, h, dh_dp, dh_dT);
305 enthalpy = h;
306
307 enthalpy.derivatives() = pressure.derivatives() * dh_dp + Th.derivatives() * dh_dT;
308 }
309 else
310 enthalpy = _water_fp->h_from_p_T(pressure.value(), Th.value() + _T_c2k);
311
312 return enthalpy;
313}
CompileTimeDerivatives::sqrt
CTSub CT_OPERATOR_BINARY CTMul CTCompareLess CTCompareGreater CTCompareEqual _arg template * sqrt(_arg)) *_arg.template D< dtag >()) CT_SIMPLE_UNARY_FUNCTION(tanh

Referenced by PorousFlowBrine::computeQpProperties(), e_from_p_T_X(), e_from_p_T_X(), h_from_p_T_X(), h_from_p_T_X(), PorousFlowBrine::initQpStatefulProperties(), and PorousFlowBrineCO2::liquidProperties().

◆ h_from_p_T_X() [2/3]

Real BrineFluidProperties::h_from_p_T_X ( Real  pressure,
Real  temperature,
Real  xnacl 
) const
overridevirtual

Definition at line 316 of file BrineFluidProperties.C.

317{
318 // Initialise the AD value (no derivatives required)
319 FPADReal p = pressure;
320 FPADReal T = temperature;
321 FPADReal x = xnacl;
322
323 _water_fp_derivs = false;
324 return h_from_p_T_X(p, T, x).value();
325}

◆ h_from_p_T_X() [3/3]

void BrineFluidProperties::h_from_p_T_X ( Real  pressure,
Real  temperature,
Real  xnacl,
Real &  h,
Real &  dh_dp,
Real &  dh_dT,
Real &  dh_dx 
) const
overridevirtual

Definition at line 328 of file BrineFluidProperties.C.

335{
336 // Initialise the AD value and derivatives
337 FPADReal p = pressure;
338 Moose::derivInsert(p.derivatives(), 0, 1.0);
339 FPADReal T = temperature;
340 Moose::derivInsert(T.derivatives(), 1, 1.0);
341 FPADReal x = xnacl;
342 Moose::derivInsert(x.derivatives(), 2, 1.0);
343
344 _water_fp_derivs = true;
345 FPADReal ad_h = h_from_p_T_X(p, T, x);
346
347 h = ad_h.value();
348 dh_dp = ad_h.derivatives()[0];
349 dh_dT = ad_h.derivatives()[1];
350 dh_dx = ad_h.derivatives()[2];
351}

◆ haliteSolubility()

Real BrineFluidProperties::haliteSolubility ( Real  temperature) const

Solubility of halite (solid NaCl) in water Originally from Potter et al., A new method for determining the solubility of salts in aqueous solutions at elevated temperatures, J.

Res. U.S. Geol. Surv., 5, 389-395 (1977). Equation describing halite solubility is repeated in Chou, Phase relations in the system NaCI-KCI-H2O. III: Solubilities of halite in vapor-saturated liquids above 445 C and redetermination of phase equilibrium properties in the system NaCI-HzO to 1000 C and 1500 bars, Geochimica et Cosmochimica Acta 51, 1965-1975 (1987). Note: this correlation is valid for 0 <= T <= 424.5 C

Parameters
temperaturetemperature (K)
Returns
halite solubility (kg/kg)

Definition at line 473 of file BrineFluidProperties.C.

474{
475 // This correlation requires temperature in Celcius
476 Real Tc = temperature - _T_c2k;
477
478 return (26.18 + 7.2e-3 * Tc + 1.06e-4 * Tc * Tc) / 100.0;
479}

◆ henryConstant() [1/3]

ADReal BrineFluidProperties::henryConstant ( const ADReal &  temperature,
const std::vector< Real > &  coeffs 
) const

Definition at line 521 of file BrineFluidProperties.C.

523{
524 Real Kh, dKh_dT;
525 henryConstant(temperature.value(), coeffs, Kh, dKh_dT);
526
527 ADReal henry = Kh;
528 henry.derivatives() = temperature.derivatives() * dKh_dT;
529
530 return henry;
531}
ADReal
DualNumber< Real, DNDerivativeType, true > ADReal
BrineFluidProperties::henryConstant
Real henryConstant(Real temperature, const std::vector< Real > &coeffs) const
IAPWS formulation of Henry's law constant for dissolution in water (implemented in water FluidPropert...
Definition BrineFluidProperties.C:506

◆ henryConstant() [2/3]

Real BrineFluidProperties::henryConstant ( Real  temperature,
const std::vector< Real > &  coeffs 
) const

IAPWS formulation of Henry's law constant for dissolution in water (implemented in water FluidProperties userobject)

Parameters
Tfluid temperature (K)
coeffsHenry's constant coefficients of gas
[out]KhHenry's constant
[out]dKh_dTderivative of Kh wrt temperature

Definition at line 506 of file BrineFluidProperties.C.

507{
508 return _water97_fp->henryConstant(temperature, coeffs);
509}
Water97FluidProperties::henryConstant
Real henryConstant(Real temperature, const std::vector< Real > &coeffs) const
IAPWS formulation of Henry's law constant for dissolution in water From Guidelines on the Henry's con...
Definition Water97FluidProperties.C:950

Referenced by PorousFlowBrineCO2::henryConstant(), and henryConstant().

◆ henryConstant() [3/3]

void BrineFluidProperties::henryConstant ( Real  temperature,
const std::vector< Real > &  coeffs,
Real &  Kh,
Real &  dKh_dT 
) const

Definition at line 512 of file BrineFluidProperties.C.

516{
517 _water97_fp->henryConstant(temperature, coeffs, Kh, dKh_dT);
518}

◆ initialize()

virtual void FluidProperties::initialize ( )
inlinefinalvirtualinherited

Implements ThreadedGeneralUserObject.

Definition at line 34 of file FluidProperties.h.

34{}

◆ k_from_p_T_X()

Real BrineFluidProperties::k_from_p_T_X ( Real  pressure,
Real  temperature,
Real  xnacl 
) const
overridevirtual

Definition at line 437 of file BrineFluidProperties.C.

438{
439 // Correlation requires molal concentration (mol/kg)
440 Real mol = massFractionToMolalConc(xnacl);
441 // Correlation requires temperature in C
442 Real Tc = temperature - _T_c2k;
443
444 Real S = 100.0 * _Mnacl * mol / (1.0 + _Mnacl * mol);
445 Real lambdaw = _water_fp->k_from_p_T(pressure, temperature);
446 Real lambda = 1.0 - (2.3434e-3 - 7.924e-6 * Tc + 3.924e-8 * Tc * Tc) * S +
447 (1.06e-5 - 2.0e-8 * Tc - 1.2e-10 * Tc * Tc) * S * S;
448
449 return lambda * lambdaw;
450}
BrineFluidProperties::massFractionToMolalConc
Real massFractionToMolalConc(Real xnacl) const
Conversion from mass fraction to molal concentration (molality)
Definition BrineFluidProperties.C:482

◆ massFractionToMolalConc()

Real BrineFluidProperties::massFractionToMolalConc ( Real  xnacl) const
protected

Conversion from mass fraction to molal concentration (molality)

Parameters
xnaclNaCl mass fraction (kg/kg)
Returns
molal concentration (mol/kg)

Definition at line 482 of file BrineFluidProperties.C.

483{
484 return xnacl / ((1.0 - xnacl) * _Mnacl);
485}

Referenced by k_from_p_T_X(), mu_from_p_T_X(), mu_from_p_T_X(), and vaporPressure().

◆ massFractionToMoleFraction() [1/2]

FPADReal BrineFluidProperties::massFractionToMoleFraction ( const FPADReal &  xnacl) const
protected

Definition at line 497 of file BrineFluidProperties.C.

498{
499 // The average molar mass of brine from the mass fraction
500 FPADReal Mbrine = molarMass(xnacl);
501 // The mole fraction is then
502 return xnacl * Mbrine / _Mnacl;
503}
BrineFluidProperties::molarMass
Real molarMass(Real xnacl) const
Average molar mass of brine.
Definition BrineFluidProperties.C:103

◆ massFractionToMoleFraction() [2/2]

Real BrineFluidProperties::massFractionToMoleFraction ( Real  xnacl) const
protected

Conversion from mass fraction to mole fraction.

Parameters
xnaclNaCl mass fraction (kg/kg)
Returns
mole fraction (mol/mol)

Definition at line 488 of file BrineFluidProperties.C.

489{
490 // The average molar mass of brine from the mass fraction
491 Real Mbrine = molarMass(xnacl);
492 // The mole fraction is then
493 return xnacl * Mbrine / _Mnacl;
494}

Referenced by cp_from_p_T_X(), h_from_p_T_X(), and rho_from_p_T_X().

◆ molarMass() [1/2]

FPADReal BrineFluidProperties::molarMass ( const FPADReal &  xnacl) const

Definition at line 97 of file BrineFluidProperties.C.

98{
99 return 1.0 / (xnacl / _Mnacl + (1.0 - xnacl) / _Mh2o);
100}

◆ molarMass() [2/2]

Real BrineFluidProperties::molarMass ( Real  xnacl) const

Average molar mass of brine.

Parameters
xnaclNaCl mass fraction (-)
Returns
average molar mass (kg/mol)

Definition at line 103 of file BrineFluidProperties.C.

104{
105 return 1.0 / (xnacl / _Mnacl + (1.0 - xnacl) / _Mh2o);
106}

Referenced by massFractionToMoleFraction(), massFractionToMoleFraction(), and rho_from_p_T_X().

◆ molarMassH2O()

Real BrineFluidProperties::molarMassH2O ( ) const

H2O molar mass.

Returns
molar mass of H2O (kg/mol)

Definition at line 115 of file BrineFluidProperties.C.

116{
117 return _Mh2o;
118}

◆ molarMassNaCl()

Real BrineFluidProperties::molarMassNaCl ( ) const

NaCl molar mass.

Returns
molar mass of NaCl (kg/mol)

Definition at line 109 of file BrineFluidProperties.C.

110{
111 return _Mnacl;
112}

◆ mu_from_p_T_X() [1/2]

Real BrineFluidProperties::mu_from_p_T_X ( Real  pressure,
Real  temperature,
Real  xnacl 
) const
overridevirtual

Definition at line 211 of file BrineFluidProperties.C.

212{
213 // Correlation requires molal concentration (mol/kg)
214 const Real mol = massFractionToMolalConc(xnacl);
215 const Real mol2 = mol * mol;
216 const Real mol3 = mol2 * mol;
217
218 // Correlation requires temperature in C
219 const Real Tc = temperature - _T_c2k;
220
221 const Real a = 1.0 + 0.0816 * mol + 0.0122 * mol2 + 0.128e-3 * mol3 +
222 0.629e-3 * Tc * (1.0 - std::exp(-0.7 * mol));
223
224 const Real water_viscosity = _water_fp->mu_from_p_T(pressure, temperature);
225
226 return a * water_viscosity;
227}

Referenced by PorousFlowBrine::computeQpProperties(), and PorousFlowBrineCO2::liquidProperties().

◆ mu_from_p_T_X() [2/2]

void BrineFluidProperties::mu_from_p_T_X ( Real  pressure,
Real  temperature,
Real  xnacl,
Real &  mu,
Real &  dmu_dp,
Real &  dmu_dT,
Real &  dmu_dx 
) const
overridevirtual

Definition at line 230 of file BrineFluidProperties.C.

237{
238 // Viscosity of water and derivatives wrt pressure and temperature
239 Real muw, dmuw_dp, dmuw_dT;
240 _water_fp->mu_from_p_T(pressure, temperature, muw, dmuw_dp, dmuw_dT);
241
242 // Correlation requires molal concentration (mol/kg)
243 Real mol = massFractionToMolalConc(xnacl);
244 Real dmol_dx = 1.0 / ((1.0 - xnacl) * (1.0 - xnacl) * _Mnacl);
245 Real mol2 = mol * mol;
246 Real mol3 = mol2 * mol;
247
248 // Correlation requires temperature in C
249 Real Tc = temperature - _T_c2k;
250
251 Real a = 1.0 + 0.0816 * mol + 0.0122 * mol2 + 0.128e-3 * mol3 +
252 0.629e-3 * Tc * (1.0 - std::exp(-0.7 * mol));
253 Real da_dx =
254 (0.0816 + 0.0244 * mol + 3.84e-4 * mol2 + 4.403e-4 * Tc * std::exp(-0.7 * mol)) * dmol_dx;
255 Real da_dT = 0.629e-3 * (1.0 - std::exp(-0.7 * mol));
256
257 mu = a * muw;
258 dmu_dp = a * dmuw_dp;
259 dmu_dx = da_dx * muw;
260 dmu_dT = da_dT * muw + a * dmuw_dT;
261}
mu
const double mu
Definition GasLiquidMassTransferTest.C:22

◆ propfunc() [1/3]

X X MultiComponentFluidProperties::propfunc ( e  ,
p  ,
T  ,
X   
)
inherited

◆ propfunc() [2/3]

X MultiComponentFluidProperties::propfunc ( h  ,
p  ,
T  ,
X   
)
inherited

◆ propfunc() [3/3]

MultiComponentFluidProperties::propfunc ( rho  ,
p  ,
T  ,
X   
)
inherited

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

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

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

@begincode p pressure [Pa] T temperature [K] X solute mass fraction [-] e specific internal energy [J/kg] rho density [kg/m^3] h specific enthalpy [J/kg] mu viscosity [Pa*s] k thermal conductivity [W/(m*K)] c speed of sound [m/s] cp constant-pressure specific heat [J/K] cv constant-volume specific heat [J/K]

As an example:

@begincode // calculate desnity given pressure, temperature and solute mass fraction: auto density = your_fluid_properties_object.rho_from_p_T_X(p, T, X);

// or use the derivative variant: Real rho = 0; // density will be stored into here Real drho_dp = 0; // derivative will be stored into here Real drho_dT = 0; // derivative will be stored into here Real drho_dX = 0; // derivative will be stored into here your_fluid_properties_object.rho_from_p_T_X(p, T, X, rho, drho_dp, drho_dT, drho_dX);

Automatic differentiation (AD) support is provided through prop_from_p_T_X(ADReal p, ADReal T, ADReal X) versions of the functions where p, T and X must be ADReal/DualNumber's calculated using all AD-supporting values.

◆ rho_from_p_T_X() [1/3]

FPADReal BrineFluidProperties::rho_from_p_T_X ( const FPADReal &  pressure,
const FPADReal &  temperature,
const FPADReal &  xnacl 
) const

Definition at line 121 of file BrineFluidProperties.C.

124{
125 using std::exp, std::sqrt;
126 // The correlation requires the pressure in bar, not Pa.
127 FPADReal pbar = pressure * 1.0e-5;
128 FPADReal pbar2 = pbar * pbar;
129 FPADReal pbar3 = pbar2 * pbar;
130
131 // The correlation requires mole fraction
132 const FPADReal Xnacl = massFractionToMoleFraction(xnacl);
133
134 const FPADReal n11 = -54.2958 - 45.7623 * exp(-9.44785e-4 * pbar);
135 const FPADReal n21 = -2.6142 - 2.39092e-4 * pbar;
136 const FPADReal n22 = 0.0356828 + 4.37235e-6 * pbar + 2.0566e-9 * pbar2;
137 const FPADReal n1x1 =
138 330.47 + 0.942876 * sqrt(pbar) + 0.0817193 * pbar - 2.47556e-8 * pbar2 + 3.45052e-10 * pbar3;
139 const FPADReal n2x1 = -0.0370751 + 0.00237723 * sqrt(pbar) + 5.42049e-5 * pbar +
140 5.84709e-9 * pbar2 - 5.99373e-13 * pbar3;
141 const FPADReal n12 = -n1x1 - n11;
142 const FPADReal n20 = 1.0 - n21 * sqrt(n22);
143 const FPADReal n23 = n2x1 - n20 - n21 * sqrt(1.0 + n22);
144
145 // The temperature Tv where the brine has the same molar volume as pure water
146 // Note: correlation uses temperature in Celcius
147 const FPADReal n1 = n1x1 + n11 * (1.0 - Xnacl) + n12 * (1.0 - Xnacl) * (1.0 - Xnacl);
148 const FPADReal n2 = n20 + n21 * sqrt(Xnacl + n22) + n23 * Xnacl;
149 const FPADReal Tv = n1 + n2 * (temperature - _T_c2k);
150
151 // The density of water at temperature Tv
152 // Note: convert Tv to Kelvin to calculate water density
153 FPADReal water_density;
154 if (_water_fp_derivs)
155 {
156 Real rho, drho_dp, drho_dT;
157 _water_fp->rho_from_p_T(pressure.value(), Tv.value() + _T_c2k, rho, drho_dp, drho_dT);
158 water_density = rho;
159
160 water_density.derivatives() = pressure.derivatives() * drho_dp + Tv.derivatives() * drho_dT;
161 }
162 else
163 water_density = _water_fp->rho_from_p_T(pressure.value(), Tv.value() + _T_c2k);
164
165 // The brine density is given by the water density scaled by the ratio of
166 // brine molar mass to pure water molar mass
167 return water_density * molarMass(xnacl) / _Mh2o;
168}
rho
const double rho
Definition GasLiquidMassTransferTest.C:27
CompileTimeDerivatives::exp
auto exp(const T &)

◆ rho_from_p_T_X() [2/3]

Real BrineFluidProperties::rho_from_p_T_X ( Real  pressure,
Real  temperature,
Real  xnacl 
) const
overridevirtual

Definition at line 171 of file BrineFluidProperties.C.

172{
173 // Initialise the AD value (no derivatives required)
174 FPADReal p = pressure;
175 FPADReal T = temperature;
176 FPADReal x = xnacl;
177
178 _water_fp_derivs = false;
179 FPADReal ad_rho = this->rho_from_p_T_X(p, T, x);
180
181 return ad_rho.value();
182}

Referenced by PorousFlowBrine::computeQpProperties(), e_from_p_T_X(), e_from_p_T_X(), PorousFlowBrine::initQpStatefulProperties(), PorousFlowBrineCO2::liquidProperties(), rho_from_p_T_X(), rho_from_p_T_X(), and PorousFlowBrineCO2::saturation().

◆ rho_from_p_T_X() [3/3]

void BrineFluidProperties::rho_from_p_T_X ( Real  pressure,
Real  temperature,
Real  xnacl,
Real &  rho,
Real &  drho_dp,
Real &  drho_dT,
Real &  drho_dx 
) const
overridevirtual

Definition at line 185 of file BrineFluidProperties.C.

192{
193 // Initialise the AD value and derivatives
194 FPADReal p = pressure;
195 Moose::derivInsert(p.derivatives(), 0, 1.0);
196 FPADReal T = temperature;
197 Moose::derivInsert(T.derivatives(), 1, 1.0);
198 FPADReal x = xnacl;
199 Moose::derivInsert(x.derivatives(), 2, 1.0);
200
201 _water_fp_derivs = true;
202 FPADReal ad_rho = this->rho_from_p_T_X(p, T, x);
203
204 rho = ad_rho.value();
205 drho_dp = ad_rho.derivatives()[0];
206 drho_dT = ad_rho.derivatives()[1];
207 drho_dx = ad_rho.derivatives()[2];
208}

◆ rho_mu_from_p_T_X() [1/3]

void MultiComponentFluidProperties::rho_mu_from_p_T_X ( ADReal  pressure,
ADReal  temperature,
ADReal  xmass,
ADReal &  rho,
ADReal &  mu 
) const
virtualinherited

Definition at line 41 of file MultiComponentFluidProperties.C.

43{
44 rho = rho_from_p_T_X(pressure, temperature, xmass);
45 mu = mu_from_p_T_X(pressure, temperature, xmass);
46}

◆ rho_mu_from_p_T_X() [2/3]

void MultiComponentFluidProperties::rho_mu_from_p_T_X ( Real  pressure,
Real  temperature,
Real  xmass,
Real &  rho,
Real &  drho_dp,
Real &  drho_dT,
Real &  drho_dx,
Real &  mu,
Real &  dmu_dp,
Real &  dmu_dT,
Real &  dmu_dx 
) const
virtualinherited

Density and viscosity and their derivatives wrt pressure, temperature and mass fraction.

Parameters
pressurefluid pressure (Pa)
temperaturefluid temperature (K)
xmassmass fraction (-)
[out]rhodensity (kg/m^3)
[out]drho_dpderivative of density wrt pressure
[out]drho_dTderivative of density wrt temperature
[out]drho_dxderivative of density wrt mass fraction
[out]muviscosity (Pa.s)
[out]dmu_dpderivative of viscosity wrt pressure
[out]dmu_dTderivative of viscosity wrt temperature
[out]dmu_dxderivative of viscosity wrt mass fraction

Definition at line 49 of file MultiComponentFluidProperties.C.

60{
61 rho_from_p_T_X(pressure, temperature, xmass, rho, drho_dp, drho_dT, drho_dx);
62 mu_from_p_T_X(pressure, temperature, xmass, mu, dmu_dp, dmu_dT, dmu_dx);
63}

◆ rho_mu_from_p_T_X() [3/3]

void MultiComponentFluidProperties::rho_mu_from_p_T_X ( Real  pressure,
Real  temperature,
Real  xmass,
Real &  rho,
Real &  mu 
) const
virtualinherited

Density and viscosity.

Parameters
pressurefluid pressure (Pa)
temperaturefluid temperature (K)
xmassmass fraction (-)
[out]rhodensity (kg/m^3)

Definition at line 33 of file MultiComponentFluidProperties.C.

35{
36 rho = rho_from_p_T_X(pressure, temperature, xmass);
37 mu = mu_from_p_T_X(pressure, temperature, xmass);
38}

◆ subdomainSetup()

virtual void FluidProperties::subdomainSetup ( )
inlinefinalvirtualinherited

Reimplemented from ThreadedGeneralUserObject.

Definition at line 38 of file FluidProperties.h.

38{}

◆ threadJoin()

virtual void FluidProperties::threadJoin ( const UserObject &  )
inlinefinalvirtualinherited

Reimplemented from ThreadedGeneralUserObject.

Definition at line 37 of file FluidProperties.h.

37{}

◆ validParams()

InputParameters BrineFluidProperties::validParams ( )
static

Definition at line 15 of file BrineFluidProperties.C.

16{
17 InputParameters params = MultiComponentFluidProperties::validParams();
18 params.addParam<UserObjectName>("water_fp",
19 "The name of the FluidProperties UserObject for water");
20 params.addClassDescription("Fluid properties for brine");
21 return params;
22}
InputParameters::addParam
void addParam(const std::string &name, const std::initializer_list< typename T::value_type > &value, const std::string &doc_string)
InputParameters::addClassDescription
void addClassDescription(const std::string &doc_string)

◆ vaporPressure()

Real BrineFluidProperties::vaporPressure ( Real  temperature,
Real  xnacl 
) const

Brine vapour pressure From Haas, Physical properties of the coexisting phases and thermochemical properties of the H2O component in boiling NaCl solutions, Geological Survey Bulletin, 1421-A (1976).

Parameters
temperaturebrine temperature (K)
xnaclsalt mass fraction (-)
Returns
brine vapour pressure (Pa)

Definition at line 453 of file BrineFluidProperties.C.

454{
455 // Correlation requires molal concentration (mol/kg)
456 Real mol = massFractionToMolalConc(xnacl);
457 Real mol2 = mol * mol;
458 Real mol3 = mol2 * mol;
459
460 Real a = 1.0 + 5.93582e-6 * mol - 5.19386e-5 * mol2 + 1.23156e-5 * mol3;
461 Real b = 1.1542e-6 * mol + 1.41254e-7 * mol2 - 1.92476e-8 * mol3 - 1.70717e-9 * mol * mol3 +
462 1.0539e-10 * mol2 * mol3;
463
464 // The temperature of pure water at the same pressure as the brine is given by
465 Real th20 = std::exp(std::log(temperature) / (a + b * temperature));
466
467 // The brine vapour pressure is then found by evaluating the saturation pressure for pure water
468 // using this effective temperature
469 return _water_fp->vaporPressure(th20);
470}
SinglePhaseFluidProperties::vaporPressure
virtual Real vaporPressure(Real T) const
Vapor pressure.
Definition SinglePhaseFluidProperties.C:372

Referenced by PorousFlowBrineCO2::solveEquilibriumMoleFractionHighTemp().

Member Data Documentation

◆ _allow_imperfect_jacobians

const bool FluidProperties::_allow_imperfect_jacobians
protectedinherited

Flag to set unimplemented Jacobian entries to zero.

Definition at line 47 of file FluidProperties.h.

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

◆ _Mh2o

Real BrineFluidProperties::_Mh2o
protected

Molar mass of water (H2O) (kg/mol)

Definition at line 199 of file BrineFluidProperties.h.

Referenced by BrineFluidProperties(), molarMass(), molarMass(), molarMassH2O(), and rho_from_p_T_X().

◆ _Mnacl

Real BrineFluidProperties::_Mnacl
protected

Molar mass of NaCl (kg/mol)

Definition at line 197 of file BrineFluidProperties.h.

Referenced by BrineFluidProperties(), k_from_p_T_X(), massFractionToMolalConc(), massFractionToMoleFraction(), massFractionToMoleFraction(), molarMass(), molarMass(), molarMassNaCl(), and mu_from_p_T_X().

◆ _nacl_fp

const SinglePhaseFluidProperties* BrineFluidProperties::_nacl_fp
protected

NaClFluidProperties UserObject.

Definition at line 194 of file BrineFluidProperties.h.

Referenced by BrineFluidProperties(), and getComponent().

◆ _R

const Real FluidProperties::_R = 8.3144598
staticinherited

Universal gas constant (J/mol/K)

Definition at line 41 of file FluidProperties.h.

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

◆ _T_c2k

const Real FluidProperties::_T_c2k
protectedinherited

Conversion of temperature from Celsius to Kelvin.

Definition at line 45 of file FluidProperties.h.

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

◆ _water97_fp

const Water97FluidProperties* BrineFluidProperties::_water97_fp
protected

Water97FluidProperties UserObject (for Henry's law)

Definition at line 190 of file BrineFluidProperties.h.

Referenced by BrineFluidProperties(), henryConstant(), and henryConstant().

◆ _water_fp

const SinglePhaseFluidProperties* BrineFluidProperties::_water_fp
protected

Water97FluidProperties UserObject.

Definition at line 192 of file BrineFluidProperties.h.

Referenced by BrineFluidProperties(), cp_from_p_T_X(), getComponent(), h_from_p_T_X(), k_from_p_T_X(), mu_from_p_T_X(), mu_from_p_T_X(), rho_from_p_T_X(), and vaporPressure().

◆ _water_fp_derivs

bool BrineFluidProperties::_water_fp_derivs
mutableprotected

Flag to indicate whether to calculate derivatives in water_fp.

Definition at line 201 of file BrineFluidProperties.h.

Referenced by e_from_p_T_X(), h_from_p_T_X(), h_from_p_T_X(), h_from_p_T_X(), rho_from_p_T_X(), rho_from_p_T_X(), and rho_from_p_T_X().

◆ NACL

const unsigned int BrineFluidProperties::NACL = 1
static

Definition at line 169 of file BrineFluidProperties.h.

Referenced by getComponent(), and TEST_F().

◆ p [1/3]

MultiComponentFluidProperties::p
inherited

Definition at line 120 of file MultiComponentFluidProperties.h.

Referenced by e_from_p_T_X(), h_from_p_T_X(), h_from_p_T_X(), rho_from_p_T_X(), and rho_from_p_T_X().

◆ p [2/3]

X MultiComponentFluidProperties::p
inherited

Definition at line 122 of file MultiComponentFluidProperties.h.

◆ p [3/3]

X X MultiComponentFluidProperties::p
inherited

Definition at line 124 of file MultiComponentFluidProperties.h.

◆ T [1/3]

MultiComponentFluidProperties::T
inherited

Definition at line 120 of file MultiComponentFluidProperties.h.

Referenced by e_from_p_T_X(), h_from_p_T_X(), h_from_p_T_X(), rho_from_p_T_X(), and rho_from_p_T_X().

◆ T [2/3]

X MultiComponentFluidProperties::T
inherited

Definition at line 122 of file MultiComponentFluidProperties.h.

◆ T [3/3]

X X MultiComponentFluidProperties::T
inherited

Definition at line 124 of file MultiComponentFluidProperties.h.

◆ WATER

const unsigned int BrineFluidProperties::WATER = 0
static

Fluid component numbers for water and NaCl.

Definition at line 168 of file BrineFluidProperties.h.

Referenced by BrineFluidPropertiesTest::buildObjects(), getComponent(), and TEST_F().

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