Computes fluid properties in (P, T) formulation. More...

#include <GeneralFluidProps.h>

Inheritance diagram for GeneralFluidProps:

Public Types
typedef DerivativeMaterialPropertyNameInterface::SymbolName	SymbolName

Public Member Functions
	GeneralFluidProps (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 void	computeQpProperties () override

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

Protected Attributes
const SinglePhaseFluidProperties &	_fluid

const VariableValue &	_eps
	Porosity.

const Real	_d
	Characteristic length $d$ used in computing the Reynolds number $Re=\rho_fVd/\mu_f$.

const ADMaterialProperty< Real > &	_pressure
	variables

const ADMaterialProperty< Real > &	_T_fluid

const ADMaterialProperty< Real > &	_rho

const ADMaterialProperty< Real > &	_speed

MaterialProperty< Real > &	_drho_dp
	Derivative of density with respect to pressure.

MaterialProperty< Real > &	_drho_dT
	Derivative of density with respect to temperature.

ADMaterialProperty< Real > &	_cp
	Isobaric specific heat capacity.

MaterialProperty< Real > &	_dcp_dp
	Derivative of isobaric specific heat with respect to pressure.

MaterialProperty< Real > &	_dcp_dT
	Derivative of isobaric specific heat with respect to temperature.

ADMaterialProperty< Real > &	_cv
	Isochoric specific heat capacity.

ADMaterialProperty< Real > &	_mu
	Dynamic viscosity.

MaterialProperty< Real > &	_dmu_dp
	Derivative of dynamic viscosity with respect to pressure.

MaterialProperty< Real > &	_dmu_dT
	Derivative of dynamic viscosity with respect to temperature.

ADMaterialProperty< Real > &	_k
	Thermal conductivity.

MaterialProperty< Real > &	_dk_dp
	Derivative of thermal conductivity with respect to pressure.

MaterialProperty< Real > &	_dk_dT
	Derivative of thermal conductivity with respect to temperature.

ADMaterialProperty< Real > &	_Pr
	Prandtl number.

MaterialProperty< Real > &	_dPr_dp
	Derivative of Prandtl number with respect to pressure.

MaterialProperty< Real > &	_dPr_dT
	Derivative of Prandtl number with respect to temperature.

ADMaterialProperty< Real > &	_Re
	Pore (particle) Reynolds number.

MaterialProperty< Real > &	_dRe_dp
	Derivative of pore Reynolds number with respect to pressure.

MaterialProperty< Real > &	_dRe_dT
	Derivative of pore Reynolds number with respect to temperature.

ADMaterialProperty< Real > &	_Re_h
	Hydraulic Reynolds number.

ADMaterialProperty< Real > &	_Re_i
	Interstitial Reynolds number.

Private Member Functions
bool	haveMaterialProperty (const std::string &prop_name)

std::vector< VariableName >	buildVariableVector (const VariableName &c1, const VariableName &c2, const VariableName &c3)

void	validateCouplingHelper (const MaterialPropertyName &base, const std::vector< VariableName > &c, const System &system, std::vector< VariableName > &missing)

bool	isNotObjectVariable (const VariableName &name)

Private Attributes
FEProblemBase &	_dmi_fe_problem

Detailed Description

Computes fluid properties in (P, T) formulation.

Definition at line 20 of file GeneralFluidProps.h.

Constructor & Destructor Documentation

◆ GeneralFluidProps()

GeneralFluidProps::GeneralFluidProps ( const InputParameters & parameters )

Definition at line 33 of file GeneralFluidProps.C.

  : DerivativeMaterialInterface<Material>(parameters),
    _fluid(UserObjectInterface::getUserObject<SinglePhaseFluidProperties>(NS::fluid)),
    _eps(coupledValue(NS::porosity)),
    _d(getParam<Real>("characteristic_length")),
 
    _pressure(getADMaterialProperty<Real>(NS::pressure)),
    _T_fluid(getADMaterialProperty<Real>(NS::T_fluid)),
    _rho(getADMaterialProperty<Real>(NS::density)),
    _speed(getADMaterialProperty<Real>(NS::speed)),
 
    _drho_dp(declarePropertyDerivative<Real>(NS::density, NS::pressure)),
    _drho_dT(declarePropertyDerivative<Real>(NS::density, NS::T_fluid)),
 
    _cp(declareADProperty<Real>(NS::cp)),
    _dcp_dp(declarePropertyDerivative<Real>(NS::cp, NS::pressure)),
    _dcp_dT(declarePropertyDerivative<Real>(NS::cp, NS::T_fluid)),
 
    _cv(declareADProperty<Real>(NS::cv)),
 
    _mu(declareADProperty<Real>(NS::mu)),
    _dmu_dp(declarePropertyDerivative<Real>(NS::mu, NS::pressure)),
    _dmu_dT(declarePropertyDerivative<Real>(NS::mu, NS::T_fluid)),
 
    _k(declareADProperty<Real>(NS::k)),
    _dk_dp(declarePropertyDerivative<Real>(NS::k, NS::pressure)),
    _dk_dT(declarePropertyDerivative<Real>(NS::k, NS::T_fluid)),
 
    _Pr(declareADProperty<Real>(NS::Prandtl)),
    _dPr_dp(declarePropertyDerivative<Real>(NS::Prandtl, NS::pressure)),
    _dPr_dT(declarePropertyDerivative<Real>(NS::Prandtl, NS::T_fluid)),
 
    _Re(declareADProperty<Real>(NS::Reynolds)),
    _dRe_dp(declarePropertyDerivative<Real>(NS::Reynolds, NS::pressure)),
    _dRe_dT(declarePropertyDerivative<Real>(NS::Reynolds, NS::T_fluid)),
 
    _Re_h(declareADProperty<Real>(NS::Reynolds_hydraulic)),
    _Re_i(declareADProperty<Real>(NS::Reynolds_interstitial))
{
}

Member Function Documentation

◆ computeQpProperties()

void GeneralFluidProps::computeQpProperties ( )

overrideprotectedvirtual

Definition at line 75 of file GeneralFluidProps.C.

{
  using std::max;
 
  auto raw_pressure = MetaPhysicL::raw_value(_pressure[_qp]);
  auto raw_T_fluid = MetaPhysicL::raw_value(_T_fluid[_qp]);
 
  // Density is not a material property because we will calculate it using
  // FluidDensityAux as needed.
  Real dummy = 0;
  _fluid.rho_from_p_T(raw_pressure, raw_T_fluid, dummy, _drho_dp[_qp], _drho_dT[_qp]);
 
  _cv[_qp] = _fluid.cv_from_p_T(_pressure[_qp], _T_fluid[_qp]);
  _cp[_qp] = _fluid.cp_from_p_T(_pressure[_qp], _T_fluid[_qp]);
  _mu[_qp] = _fluid.mu_from_p_T(_pressure[_qp], _T_fluid[_qp]);
  _k[_qp] = _fluid.k_from_p_T(_pressure[_qp], _T_fluid[_qp]);
  _fluid.cp_from_p_T(raw_pressure, raw_T_fluid, dummy, _dcp_dp[_qp], _dcp_dT[_qp]);
  _fluid.mu_from_p_T(raw_pressure, raw_T_fluid, dummy, _dmu_dp[_qp], _dmu_dT[_qp]);
  _fluid.k_from_p_T(raw_pressure, raw_T_fluid, dummy, _dk_dp[_qp], _dk_dT[_qp]);
 
  static constexpr Real small_number = 1e-8;
 
  _Pr[_qp] = HeatTransferUtils::prandtl(_cp[_qp], _mu[_qp], max(_k[_qp], small_number));
  _dPr_dp[_qp] = NS::prandtlPropertyDerivative(MetaPhysicL::raw_value(_mu[_qp]),
                                               MetaPhysicL::raw_value(_cp[_qp]),
                                               MetaPhysicL::raw_value(_k[_qp]),
                                               _dmu_dp[_qp],
                                               _dcp_dp[_qp],
                                               _dk_dp[_qp]);
  _dPr_dT[_qp] = NS::prandtlPropertyDerivative(MetaPhysicL::raw_value(_mu[_qp]),
                                               MetaPhysicL::raw_value(_cp[_qp]),
                                               MetaPhysicL::raw_value(_k[_qp]),
                                               _dmu_dT[_qp],
                                               _dcp_dT[_qp],
                                               _dk_dT[_qp]);
 
  // (pore / particle) Reynolds number based on superficial velocity and
  // characteristic length. Only call Reynolds() one time to compute all three so that
  // we don't redundantly check that viscosity is not too close to zero.
  _Re[_qp] = max(HeatTransferUtils::reynolds(
                     _rho[_qp], _eps[_qp] * _speed[_qp], _d, max(_mu[_qp], small_number)),
                 1.0);
  _dRe_dp[_qp] = NS::reynoldsPropertyDerivative(MetaPhysicL::raw_value(_Re[_qp]),
                                                MetaPhysicL::raw_value(_rho[_qp]),
                                                MetaPhysicL::raw_value(_mu[_qp]),
                                                _drho_dp[_qp],
                                                _dmu_dp[_qp]);
  _dRe_dT[_qp] = NS::reynoldsPropertyDerivative(MetaPhysicL::raw_value(_Re[_qp]),
                                                MetaPhysicL::raw_value(_rho[_qp]),
                                                MetaPhysicL::raw_value(_mu[_qp]),
                                                _drho_dT[_qp],
                                                _dmu_dT[_qp]);
 
  // (hydraulic) Reynolds number
  _Re_h[_qp] = _Re[_qp] / max(1 - _eps[_qp], small_number);
 
  // (interstitial) Reynolds number
  _Re_i[_qp] = _Re[_qp] / _eps[_qp];
}

◆ getUserObject()

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

protected

◆ validParams()

InputParameters GeneralFluidProps::validParams ( )

static

Definition at line 19 of file GeneralFluidProps.C.

{
  auto params = Material::validParams();
  params.addRequiredParam<UserObjectName>(NS::fluid, "Fluid properties userobject");
  params.addClassDescription("Computes fluid properties using a (P, T) formulation");
 
  params.addRequiredCoupledVar(NS::porosity, "porosity");
  params.addRequiredRangeCheckedParam<Real>(
      "characteristic_length",
      "characteristic_length > 0.0 ",
      "characteristic length for Reynolds number calculation");
  return params;
}

Member Data Documentation

◆ _cp

ADMaterialProperty<Real>& GeneralFluidProps::_cp

protected

Isobaric specific heat capacity.

Definition at line 53 of file GeneralFluidProps.h.

Referenced by computeQpProperties().

◆ _cv

ADMaterialProperty<Real>& GeneralFluidProps::_cv

protected

Isochoric specific heat capacity.

Definition at line 62 of file GeneralFluidProps.h.

Referenced by computeQpProperties().

◆ _d

const Real GeneralFluidProps::_d

protected

Characteristic length $d$ used in computing the Reynolds number $Re=\rho_fVd/\mu_f$.

Definition at line 38 of file GeneralFluidProps.h.

Referenced by computeQpProperties().

◆ _dcp_dp

MaterialProperty<Real>& GeneralFluidProps::_dcp_dp

protected

Derivative of isobaric specific heat with respect to pressure.

Definition at line 56 of file GeneralFluidProps.h.

Referenced by computeQpProperties().

◆ _dcp_dT

MaterialProperty<Real>& GeneralFluidProps::_dcp_dT

protected

Derivative of isobaric specific heat with respect to temperature.

Definition at line 59 of file GeneralFluidProps.h.

Referenced by computeQpProperties().

◆ _dk_dp

MaterialProperty<Real>& GeneralFluidProps::_dk_dp

protected

Derivative of thermal conductivity with respect to pressure.

Definition at line 77 of file GeneralFluidProps.h.

Referenced by computeQpProperties().

◆ _dk_dT

MaterialProperty<Real>& GeneralFluidProps::_dk_dT

protected

Derivative of thermal conductivity with respect to temperature.

Definition at line 80 of file GeneralFluidProps.h.

Referenced by computeQpProperties().

◆ _dmu_dp

MaterialProperty<Real>& GeneralFluidProps::_dmu_dp

protected

Derivative of dynamic viscosity with respect to pressure.

Definition at line 68 of file GeneralFluidProps.h.

Referenced by computeQpProperties().

◆ _dmu_dT

MaterialProperty<Real>& GeneralFluidProps::_dmu_dT

protected

Derivative of dynamic viscosity with respect to temperature.

Definition at line 71 of file GeneralFluidProps.h.

Referenced by computeQpProperties().

◆ _dPr_dp

MaterialProperty<Real>& GeneralFluidProps::_dPr_dp

protected

Derivative of Prandtl number with respect to pressure.

Definition at line 88 of file GeneralFluidProps.h.

Referenced by computeQpProperties().

◆ _dPr_dT

MaterialProperty<Real>& GeneralFluidProps::_dPr_dT

protected

Derivative of Prandtl number with respect to temperature.

Definition at line 93 of file GeneralFluidProps.h.

Referenced by computeQpProperties().

◆ _dRe_dp

MaterialProperty<Real>& GeneralFluidProps::_dRe_dp

protected

Derivative of pore Reynolds number with respect to pressure.

Definition at line 99 of file GeneralFluidProps.h.

Referenced by computeQpProperties().

◆ _dRe_dT

MaterialProperty<Real>& GeneralFluidProps::_dRe_dT

protected

Derivative of pore Reynolds number with respect to temperature.

Definition at line 102 of file GeneralFluidProps.h.

Referenced by computeQpProperties().

◆ _drho_dp

MaterialProperty<Real>& GeneralFluidProps::_drho_dp

protected

Derivative of density with respect to pressure.

Definition at line 47 of file GeneralFluidProps.h.

Referenced by computeQpProperties().

◆ _drho_dT

MaterialProperty<Real>& GeneralFluidProps::_drho_dT

protected

Derivative of density with respect to temperature.

Definition at line 50 of file GeneralFluidProps.h.

Referenced by computeQpProperties().

◆ _eps

const VariableValue& GeneralFluidProps::_eps

protected

Porosity.

Definition at line 32 of file GeneralFluidProps.h.

Referenced by computeQpProperties().

◆ _fluid

const SinglePhaseFluidProperties& GeneralFluidProps::_fluid

protected

Definition at line 29 of file GeneralFluidProps.h.

Referenced by computeQpProperties().

◆ _k

ADMaterialProperty<Real>& GeneralFluidProps::_k

protected

Thermal conductivity.

Definition at line 74 of file GeneralFluidProps.h.

Referenced by computeQpProperties().

◆ _mu

ADMaterialProperty<Real>& GeneralFluidProps::_mu

protected

Dynamic viscosity.

Definition at line 65 of file GeneralFluidProps.h.

Referenced by computeQpProperties().

◆ _Pr

ADMaterialProperty<Real>& GeneralFluidProps::_Pr

protected

Prandtl number.

Definition at line 83 of file GeneralFluidProps.h.

Referenced by computeQpProperties().

◆ _pressure

const ADMaterialProperty<Real>& GeneralFluidProps::_pressure

protected

variables

Definition at line 41 of file GeneralFluidProps.h.

Referenced by computeQpProperties().

◆ _Re

ADMaterialProperty<Real>& GeneralFluidProps::_Re

protected

Pore (particle) Reynolds number.

Definition at line 96 of file GeneralFluidProps.h.

Referenced by computeQpProperties().

◆ _Re_h

ADMaterialProperty<Real>& GeneralFluidProps::_Re_h

protected

Hydraulic Reynolds number.

Definition at line 105 of file GeneralFluidProps.h.

Referenced by computeQpProperties().

◆ _Re_i

ADMaterialProperty<Real>& GeneralFluidProps::_Re_i

protected

Interstitial Reynolds number.

Definition at line 108 of file GeneralFluidProps.h.

Referenced by computeQpProperties().

◆ _rho

const ADMaterialProperty<Real>& GeneralFluidProps::_rho

protected

Definition at line 43 of file GeneralFluidProps.h.

Referenced by computeQpProperties().

◆ _speed

const ADMaterialProperty<Real>& GeneralFluidProps::_speed

protected

Definition at line 44 of file GeneralFluidProps.h.

Referenced by computeQpProperties().

◆ _T_fluid

const ADMaterialProperty<Real>& GeneralFluidProps::_T_fluid

protected

Definition at line 42 of file GeneralFluidProps.h.

Referenced by computeQpProperties().

The documentation for this class was generated from the following files:

navier_stokes/include/materials/GeneralFluidProps.h
navier_stokes/src/materials/GeneralFluidProps.C

Public Types

Public Member Functions

Static Public Member Functions

Protected Member Functions

Protected Attributes

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ GeneralFluidProps()

Member Function Documentation

◆ computeQpProperties()

◆ getUserObject()

◆ validParams()

Member Data Documentation

◆ _cp

◆ _cv

◆ _d

◆ _dcp_dp

◆ _dcp_dT

◆ _dk_dp

◆ _dk_dT

◆ _dmu_dp

◆ _dmu_dT

◆ _dPr_dp

◆ _dPr_dT

◆ _dRe_dp

◆ _dRe_dT

◆ _drho_dp

◆ _drho_dT

◆ _eps

◆ _fluid

◆ _k

◆ _mu

◆ _Pr

◆ _pressure

◆ _Re

◆ _Re_h

◆ _Re_i

◆ _rho

◆ _speed

◆ _T_fluid