https://mooseframework.inl.gov
Public Types | Public Member Functions | Static Public Member Functions | Protected Member Functions | Protected Attributes | List of all members
PorousFlowRelativePermeabilityCoreyTempl< is_ad > Class Template Reference

Material to calculate Corey-type relative permeability of an arbitrary phase given the effective saturation and Corey exponent of that phase. More...

#include <PorousFlowRelativePermeabilityCorey.h>

Inheritance diagram for PorousFlowRelativePermeabilityCoreyTempl< is_ad >:
[legend]

Public Types

typedef DerivativeMaterialPropertyNameInterface::SymbolName SymbolName
 

Public Member Functions

 PorousFlowRelativePermeabilityCoreyTempl (const InputParameters &parameters)
 
const GenericMaterialProperty< U, is_ad > & getDefaultMaterialProperty (const std::string &name)
 
const GenericMaterialProperty< U, is_ad > & getDefaultMaterialPropertyByName (const std::string &name)
 
void validateDerivativeMaterialPropertyBase (const std::string &base)
 
const MaterialPropertyName derivativePropertyName (const MaterialPropertyName &base, const std::vector< SymbolName > &c) const
 
const MaterialPropertyName derivativePropertyNameFirst (const MaterialPropertyName &base, const SymbolName &c1) const
 
const MaterialPropertyName derivativePropertyNameSecond (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2) const
 
const MaterialPropertyName derivativePropertyNameThird (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2, const SymbolName &c3) const
 
GenericMaterialProperty< U, is_ad > & declarePropertyDerivative (const std::string &base, const std::vector< VariableName > &c)
 
GenericMaterialProperty< U, is_ad > & declarePropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)
 
GenericMaterialProperty< U, is_ad > & declarePropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")
 
GenericMaterialProperty< U, is_ad > & declarePropertyDerivative (const std::string &base, const std::vector< VariableName > &c)
 
GenericMaterialProperty< U, is_ad > & declarePropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)
 
GenericMaterialProperty< U, is_ad > & declarePropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivative (const std::string &base, const std::vector< VariableName > &c)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, unsigned int v2, unsigned int v3=libMesh::invalid_uint)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivative (const std::string &base, unsigned int v1, unsigned int v2=libMesh::invalid_uint, unsigned int v3=libMesh::invalid_uint)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivative (const std::string &base, const std::vector< VariableName > &c)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, unsigned int v2, unsigned int v3=libMesh::invalid_uint)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivative (const std::string &base, unsigned int v1, unsigned int v2=libMesh::invalid_uint, unsigned int v3=libMesh::invalid_uint)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< VariableName > &c)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< SymbolName > &c)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< VariableName > &c)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< SymbolName > &c)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")
 
void validateCoupling (const MaterialPropertyName &base, const std::vector< VariableName > &c, bool validate_aux=true)
 
void validateCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")
 
void validateCoupling (const MaterialPropertyName &base, const std::vector< VariableName > &c, bool validate_aux=true)
 
void validateCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")
 
void validateNonlinearCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")
 
void validateNonlinearCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")
 
const MaterialPropertyName propertyName (const MaterialPropertyName &base, const std::vector< SymbolName > &c) const
 
const MaterialPropertyName propertyName (const MaterialPropertyName &base, const std::vector< SymbolName > &c) const
 
const MaterialPropertyName propertyNameFirst (const MaterialPropertyName &base, const SymbolName &c1) const
 
const MaterialPropertyName propertyNameFirst (const MaterialPropertyName &base, const SymbolName &c1) const
 
const MaterialPropertyName propertyNameSecond (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2) const
 
const MaterialPropertyName propertyNameSecond (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2) const
 
const MaterialPropertyName propertyNameThird (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2, const SymbolName &c3) const
 
const MaterialPropertyName propertyNameThird (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2, const SymbolName &c3) const
 

Static Public Member Functions

static InputParameters validParams ()
 

Protected Member Functions

virtual GenericReal< is_ad > relativePermeability (GenericReal< is_ad > seff) const override
 Relative permeability equation (must be overriden in derived class) More...
 
virtual Real dRelativePermeability (Real seff) const override
 Derivative of relative permeability with respect to effective saturation. More...
 
virtual void computeQpProperties () override
 
virtual GenericReal< is_ad > effectiveSaturation (GenericReal< is_ad > saturation) const
 Effective saturation of fluid phase. More...
 

Protected Attributes

const Real _n
 Corey exponent n for the specified phase. More...
 
const Real _scaling
 Relative permeability is multiplied by this quantity. More...
 
const GenericMaterialProperty< std::vector< Real >, is_ad > & _saturation
 Saturation material property. More...
 
GenericMaterialProperty< Real, is_ad > & _relative_permeability
 Relative permeability material property. More...
 
MaterialProperty< Real > *const _drelative_permeability_ds
 Derivative of relative permeability wrt phase saturation. More...
 
const Real _s_res
 Residual saturation of specified phase. More...
 
const Real _sum_s_res
 Sum of residual saturations over all phases. More...
 
const Real _dseff_ds
 Derivative of effective saturation with respect to saturation. More...
 
const unsigned int _phase_num
 Phase number of fluid. More...
 
const std::string _phase
 Stringified fluid phase number. More...
 

Detailed Description

template<bool is_ad>
class PorousFlowRelativePermeabilityCoreyTempl< is_ad >

Material to calculate Corey-type relative permeability of an arbitrary phase given the effective saturation and Corey exponent of that phase.

From Corey, A. T., The interrelation between gas and oil relative permeabilities, Prod. Month. 38-41 (1954)

Definition at line 22 of file PorousFlowRelativePermeabilityCorey.h.

Constructor & Destructor Documentation

◆ PorousFlowRelativePermeabilityCoreyTempl()

template<bool is_ad>
PorousFlowRelativePermeabilityCoreyTempl< is_ad >::PorousFlowRelativePermeabilityCoreyTempl ( const InputParameters &  parameters)

Definition at line 27 of file PorousFlowRelativePermeabilityCorey.C.

29  : PorousFlowRelativePermeabilityBaseTempl<is_ad>(parameters),
30  _n(this->template getParam<Real>("n"))
31 {
32 }
PorousFlowRelativePermeabilityCoreyTempl::_n
const Real _n
Corey exponent n for the specified phase.
Definition: PorousFlowRelativePermeabilityCorey.h:35
PorousFlowRelativePermeabilityBaseTempl
Base class for PorousFlow relative permeability materials.
Definition: PorousFlowRelativePermeabilityBase.h:20

Member Function Documentation

◆ computeQpProperties()

template<bool is_ad>
void PorousFlowRelativePermeabilityBaseTempl< is_ad >::computeQpProperties ( )
overrideprotectedvirtualinherited

Definition at line 65 of file PorousFlowRelativePermeabilityBase.C.

66 {
67  // Effective saturation
68  GenericReal<is_ad> seff = effectiveSaturation(_saturation[_qp][_phase_num]);
69  GenericReal<is_ad> relperm;
70  Real drelperm;
71 
72  if (seff < 0.0)
73  {
74  // Relative permeability is 0 for saturation less than residual
75  relperm = 0.0;
76  drelperm = 0.0;
77  }
78  else if (seff >= 0.0 && seff <= 1)
79  {
80  relperm = relativePermeability(seff);
81  drelperm = dRelativePermeability(MetaPhysicL::raw_value(seff));
82  }
83  else // seff > 1
84  {
85  // Relative permeability is 1 when fully saturated
86  relperm = 1.0;
87  drelperm = 0.0;
88  }
89 
90  _relative_permeability[_qp] = relperm * _scaling;
91 
92  if (!is_ad)
93  (*_drelative_permeability_ds)[_qp] = drelperm * _dseff_ds * _scaling;
94 }
PorousFlowRelativePermeabilityBaseTempl::dRelativePermeability
virtual Real dRelativePermeability(Real seff) const =0
Derivative of relative permeability with respect to effective saturation.
GenericReal
Moose::GenericType< Real, is_ad > GenericReal
PorousFlowRelativePermeabilityBaseTempl::_saturation
const GenericMaterialProperty< std::vector< Real >, is_ad > & _saturation
Saturation material property.
Definition: PorousFlowRelativePermeabilityBase.h:55
PorousFlowRelativePermeabilityBaseTempl::relativePermeability
virtual GenericReal< is_ad > relativePermeability(GenericReal< is_ad > seff) const =0
Relative permeability equation (must be overriden in derived class)
MetaPhysicL::raw_value
auto raw_value(const Eigen::Map< T > &in)
PorousFlowRelativePermeabilityBaseTempl::_scaling
const Real _scaling
Relative permeability is multiplied by this quantity.
Definition: PorousFlowRelativePermeabilityBase.h:52
PorousFlowMaterialBase::_phase_num
const unsigned int _phase_num
Phase number of fluid.
Definition: PorousFlowMaterialBase.h:31
PorousFlowRelativePermeabilityBaseTempl::effectiveSaturation
virtual GenericReal< is_ad > effectiveSaturation(GenericReal< is_ad > saturation) const
Effective saturation of fluid phase.
Definition: PorousFlowRelativePermeabilityBase.C:98
Real
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
PorousFlowRelativePermeabilityBaseTempl::_dseff_ds
const Real _dseff_ds
Derivative of effective saturation with respect to saturation.
Definition: PorousFlowRelativePermeabilityBase.h:70
PorousFlowRelativePermeabilityBaseTempl::_relative_permeability
GenericMaterialProperty< Real, is_ad > & _relative_permeability
Relative permeability material property.
Definition: PorousFlowRelativePermeabilityBase.h:58

◆ dRelativePermeability()

template<bool is_ad>
Real PorousFlowRelativePermeabilityCoreyTempl< is_ad >::dRelativePermeability ( Real  seff) const
overrideprotectedvirtual

Derivative of relative permeability with respect to effective saturation.

Parameters
seffeffective saturation
Returns
derivative of relative permeability wrt effective saturation

Implements PorousFlowRelativePermeabilityBaseTempl< is_ad >.

Definition at line 43 of file PorousFlowRelativePermeabilityCorey.C.

44 {
45  return _n * std::pow(seff, _n - 1.0);
46 }
PorousFlowRelativePermeabilityCoreyTempl::_n
const Real _n
Corey exponent n for the specified phase.
Definition: PorousFlowRelativePermeabilityCorey.h:35
std::pow
MooseUnits pow(const MooseUnits &, int)

◆ effectiveSaturation()

template<bool is_ad>
GenericReal< is_ad > PorousFlowRelativePermeabilityBaseTempl< is_ad >::effectiveSaturation ( GenericReal< is_ad >  saturation) const
protectedvirtualinherited

Effective saturation of fluid phase.

Parameters
saturationtrue saturation
Returns
effective saturation

Definition at line 98 of file PorousFlowRelativePermeabilityBase.C.

100 {
101  return (saturation - _s_res) / (1.0 - _sum_s_res);
102 }
PorousFlowRelativePermeabilityBaseTempl::_sum_s_res
const Real _sum_s_res
Sum of residual saturations over all phases.
Definition: PorousFlowRelativePermeabilityBase.h:67
PorousFlowRelativePermeabilityBaseTempl::_s_res
const Real _s_res
Residual saturation of specified phase.
Definition: PorousFlowRelativePermeabilityBase.h:64

◆ relativePermeability()

template<bool is_ad>
GenericReal< is_ad > PorousFlowRelativePermeabilityCoreyTempl< is_ad >::relativePermeability ( GenericReal< is_ad >  seff) const
overrideprotectedvirtual

Relative permeability equation (must be overriden in derived class)

Parameters
seffeffective saturation
Returns
relative permeability

Implements PorousFlowRelativePermeabilityBaseTempl< is_ad >.

Definition at line 36 of file PorousFlowRelativePermeabilityCorey.C.

37 {
38  return std::pow(seff, _n);
39 }
PorousFlowRelativePermeabilityCoreyTempl::_n
const Real _n
Corey exponent n for the specified phase.
Definition: PorousFlowRelativePermeabilityCorey.h:35
std::pow
MooseUnits pow(const MooseUnits &, int)

◆ validParams()

template<bool is_ad>
InputParameters PorousFlowRelativePermeabilityCoreyTempl< is_ad >::validParams ( )
static

Definition at line 17 of file PorousFlowRelativePermeabilityCorey.C.

18 {
19  InputParameters params = PorousFlowRelativePermeabilityBaseTempl<is_ad>::validParams();
20  params.addRequiredParam<Real>("n", "The Corey exponent of the phase.");
21  params.addClassDescription("This Material calculates relative permeability of the fluid phase, "
22  "using the simple Corey model ((S-S_res)/(1-sum(S_res)))^n");
23  return params;
24 }
InputParameters::addRequiredParam
void addRequiredParam(const std::string &name, const std::string &doc_string)
Real
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
InputParameters::addClassDescription
void addClassDescription(const std::string &doc_string)

Member Data Documentation

◆ _drelative_permeability_ds

template<bool is_ad>
MaterialProperty<Real>* const PorousFlowRelativePermeabilityBaseTempl< is_ad >::_drelative_permeability_ds
protectedinherited

Derivative of relative permeability wrt phase saturation.

Definition at line 61 of file PorousFlowRelativePermeabilityBase.h.

◆ _dseff_ds

template<bool is_ad>
const Real PorousFlowRelativePermeabilityBaseTempl< is_ad >::_dseff_ds
protectedinherited

Derivative of effective saturation with respect to saturation.

Definition at line 70 of file PorousFlowRelativePermeabilityBase.h.

◆ _n

template<bool is_ad>
const Real PorousFlowRelativePermeabilityCoreyTempl< is_ad >::_n
protected

Corey exponent n for the specified phase.

Definition at line 35 of file PorousFlowRelativePermeabilityCorey.h.

◆ _phase

const std::string PorousFlowMaterialBase::_phase
protectedinherited

Stringified fluid phase number.

Definition at line 34 of file PorousFlowMaterialBase.h.

◆ _phase_num

const unsigned int PorousFlowMaterialBase::_phase_num
protectedinherited

Phase number of fluid.

Definition at line 31 of file PorousFlowMaterialBase.h.

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

◆ _relative_permeability

template<bool is_ad>
GenericMaterialProperty<Real, is_ad>& PorousFlowRelativePermeabilityBaseTempl< is_ad >::_relative_permeability
protectedinherited

Relative permeability material property.

Definition at line 58 of file PorousFlowRelativePermeabilityBase.h.

◆ _s_res

template<bool is_ad>
const Real PorousFlowRelativePermeabilityBaseTempl< is_ad >::_s_res
protectedinherited

Residual saturation of specified phase.

Definition at line 64 of file PorousFlowRelativePermeabilityBase.h.

Referenced by PorousFlowRelativePermeabilityBaseTempl< is_ad >::PorousFlowRelativePermeabilityBaseTempl().

◆ _saturation

template<bool is_ad>
const GenericMaterialProperty<std::vector<Real>, is_ad>& PorousFlowRelativePermeabilityBaseTempl< is_ad >::_saturation
protectedinherited

Saturation material property.

Definition at line 55 of file PorousFlowRelativePermeabilityBase.h.

◆ _scaling

template<bool is_ad>
const Real PorousFlowRelativePermeabilityBaseTempl< is_ad >::_scaling
protectedinherited

Relative permeability is multiplied by this quantity.

Definition at line 52 of file PorousFlowRelativePermeabilityBase.h.

◆ _sum_s_res

template<bool is_ad>
const Real PorousFlowRelativePermeabilityBaseTempl< is_ad >::_sum_s_res
protectedinherited

Sum of residual saturations over all phases.

Definition at line 67 of file PorousFlowRelativePermeabilityBase.h.

Referenced by PorousFlowRelativePermeabilityBaseTempl< is_ad >::PorousFlowRelativePermeabilityBaseTempl().

The documentation for this class was generated from the following files:
Generated on Fri Jul 18 2025 13:50:17 for https://mooseframework.inl.gov by   doxygen 1.8.14