PorousFlowTotalGravitationalDensityFullySaturatedFromPorosity Class Reference

Material designed to provide the density of the porous medium for the fully-saturated case. More...

#include <PorousFlowTotalGravitationalDensityFullySaturatedFromPorosity.h>

Inheritance diagram for PorousFlowTotalGravitationalDensityFullySaturatedFromPorosity:

Public Member Functions
	PorousFlowTotalGravitationalDensityFullySaturatedFromPorosity (const InputParameters &parameters)

Protected Member Functions
void	initQpStatefulProperties () override
void	computeQpProperties () override

Protected Attributes
const Real	_rho_s
	Solid density. More...
const MaterialProperty< std::vector< Real > > &	_rho_f_qp
	Fluid density at qps. More...
const MaterialProperty< Real > &	_porosity_qp
	Porosity at qps. More...
const MaterialProperty< std::vector< std::vector< Real > > > &	_drho_f_qp_dvar
	d(rho_f)/d(PorousFlow variable) More...
const MaterialProperty< std::vector< Real > > &	_dporosity_qp_dvar
	d(porosity)/d(PorousFlow variable) More...
MaterialProperty< Real > &	_gravdensity
	Computed density at quadpoints. More...
MaterialProperty< std::vector< Real > > &	_dgravdensity_dvar
	d(density)/d(PorousFlow variable) More...
const unsigned int	_num_phases
	Number of phases. More...
const unsigned int	_num_components
	Number of fluid components. More...
const unsigned int	_num_var
	Number of PorousFlow variables. More...

Detailed Description

Material designed to provide the density of the porous medium for the fully-saturated case.

Density is calculated as a weighted average of the fluid and solid densities: density = phi * rho_f + (1 - phi) * rho_s where phi is porosity, rho_f is fluid density and rho_s is solid density (assumed constant).

Definition at line 27 of file PorousFlowTotalGravitationalDensityFullySaturatedFromPorosity.h.

Constructor & Destructor Documentation

PorousFlowTotalGravitationalDensityFullySaturatedFromPorosity()

PorousFlowTotalGravitationalDensityFullySaturatedFromPorosity::PorousFlowTotalGravitationalDensityFullySaturatedFromPorosity

(

const InputParameters &

parameters

)

Definition at line 31 of file PorousFlowTotalGravitationalDensityFullySaturatedFromPorosity.C.

  : PorousFlowTotalGravitationalDensityBase(parameters),
    _rho_s(getParam<Real>("rho_s")),
    _rho_f_qp(getMaterialProperty<std::vector<Real>>("PorousFlow_fluid_phase_density_qp")),
    _porosity_qp(getMaterialProperty<Real>("PorousFlow_porosity_qp")),
    _drho_f_qp_dvar(getMaterialProperty<std::vector<std::vector<Real>>>(
        "dPorousFlow_fluid_phase_density_qp_dvar")),
    _dporosity_qp_dvar(getMaterialProperty<std::vector<Real>>("dPorousFlow_porosity_qp_dvar"))
{
  if (_num_phases != 1)
    mooseError("The Dictator proclaims that the number of phases is ",
               _dictator.numPhases(),
               " whereas PorousFlowThermalConductivityFromPorosity can only be used for 1-phase "
               "simulations. Be aware that the Dictator has noted your mistake.");
 
  if (_nodal_material == true)
    mooseError("PorousFlowTotalGravitationalDensityFullySaturatedFromPorosity is only defined for "
               "at_nodes = false");
}

Member Function Documentation

computeQpProperties()

void PorousFlowTotalGravitationalDensityFullySaturatedFromPorosity::computeQpProperties ( )

overrideprotected

Definition at line 60 of file PorousFlowTotalGravitationalDensityFullySaturatedFromPorosity.C.

{
  const unsigned ph = 0;
  _gravdensity[_qp] = _rho_s * (1.0 - _porosity_qp[_qp]) + _rho_f_qp[_qp][ph] * _porosity_qp[_qp];
 
  _dgravdensity_dvar[_qp].resize(_num_var);
  for (unsigned int v = 0; v < _num_var; ++v)
  {
    _dgravdensity_dvar[_qp][v] = _dporosity_qp_dvar[_qp][v] * (_rho_f_qp[_qp][ph] - _rho_s) +
                                 _drho_f_qp_dvar[_qp][ph][v] * _porosity_qp[_qp];
  }
}

initQpStatefulProperties()

void PorousFlowTotalGravitationalDensityFullySaturatedFromPorosity::initQpStatefulProperties ( )

overrideprotected

Definition at line 53 of file PorousFlowTotalGravitationalDensityFullySaturatedFromPorosity.C.

{
  const unsigned ph = 0;
  _gravdensity[_qp] = _rho_s * (1.0 - _porosity_qp[_qp]) + _rho_f_qp[_qp][ph] * _porosity_qp[_qp];
}

Member Data Documentation

_dgravdensity_dvar

MaterialProperty<std::vector<Real> >& PorousFlowTotalGravitationalDensityBase::_dgravdensity_dvar

protectedinherited

d(density)/d(PorousFlow variable)

Definition at line 32 of file PorousFlowTotalGravitationalDensityBase.h.

Referenced by computeQpProperties().

_dporosity_qp_dvar

const MaterialProperty<std::vector<Real> >& PorousFlowTotalGravitationalDensityFullySaturatedFromPorosity::_dporosity_qp_dvar

protected

d(porosity)/d(PorousFlow variable)

Definition at line 50 of file PorousFlowTotalGravitationalDensityFullySaturatedFromPorosity.h.

Referenced by computeQpProperties().

_drho_f_qp_dvar

const MaterialProperty<std::vector<std::vector<Real> > >& PorousFlowTotalGravitationalDensityFullySaturatedFromPorosity::_drho_f_qp_dvar

protected

d(rho_f)/d(PorousFlow variable)

Definition at line 47 of file PorousFlowTotalGravitationalDensityFullySaturatedFromPorosity.h.

Referenced by computeQpProperties().

_gravdensity

MaterialProperty<Real>& PorousFlowTotalGravitationalDensityBase::_gravdensity

protectedinherited

Computed density at quadpoints.

Definition at line 29 of file PorousFlowTotalGravitationalDensityBase.h.

Referenced by computeQpProperties(), and initQpStatefulProperties().

_num_components

const unsigned int PorousFlowMaterialVectorBase::_num_components

protectedinherited

Number of fluid components.

Definition at line 33 of file PorousFlowMaterialVectorBase.h.

Referenced by PorousFlowDiffusivityBase::computeQpProperties(), PorousFlowMassFractionAqueousEquilibriumChemistry::computeQpProperties(), PorousFlowMassFraction::computeQpProperties(), PorousFlowAqueousPreDisChemistry::PorousFlowAqueousPreDisChemistry(), PorousFlowDiffusivityBase::PorousFlowDiffusivityBase(), and PorousFlowMassFraction::PorousFlowMassFraction().

_num_phases

const unsigned int PorousFlowMaterialVectorBase::_num_phases

protectedinherited

Number of phases.

Definition at line 30 of file PorousFlowMaterialVectorBase.h.

Referenced by PorousFlowDiffusivityBase::computeQpProperties(), PorousFlowDiffusivityMillingtonQuirk::computeQpProperties(), PorousFlowEffectiveFluidPressure::computeQpProperties(), PorousFlowMassFraction::computeQpProperties(), PorousFlowJoiner::computeQpProperties(), PorousFlowEffectiveFluidPressure::initQpStatefulProperties(), PorousFlowJoiner::initQpStatefulProperties(), PorousFlowDiffusivityBase::PorousFlowDiffusivityBase(), PorousFlowDiffusivityConst::PorousFlowDiffusivityConst(), PorousFlowJoiner::PorousFlowJoiner(), PorousFlowMassFraction::PorousFlowMassFraction(), PorousFlowThermalConductivityFromPorosity::PorousFlowThermalConductivityFromPorosity(), PorousFlowThermalConductivityIdeal::PorousFlowThermalConductivityIdeal(), and PorousFlowTotalGravitationalDensityFullySaturatedFromPorosity().

_num_var

const unsigned int PorousFlowMaterialVectorBase::_num_var

protectedinherited

Number of PorousFlow variables.

Definition at line 36 of file PorousFlowMaterialVectorBase.h.

Referenced by PorousFlowDiffusivityBase::computeQpProperties(), PorousFlowPermeabilityConst::computeQpProperties(), PorousFlowVolumetricStrain::computeQpProperties(), PorousFlowDiffusivityMillingtonQuirk::computeQpProperties(), PorousFlowAqueousPreDisChemistry::computeQpProperties(), PorousFlowThermalConductivityIdeal::computeQpProperties(), PorousFlowEffectiveFluidPressure::computeQpProperties(), PorousFlowPermeabilityTensorFromVar::computeQpProperties(), PorousFlowPorosityConst::computeQpProperties(), PorousFlowMassFractionAqueousEquilibriumChemistry::computeQpProperties(), PorousFlowThermalConductivityFromPorosity::computeQpProperties(), PorousFlowMatrixInternalEnergy::computeQpProperties(), PorousFlowPorosityHMBiotModulus::computeQpProperties(), PorousFlowPermeabilityConstFromVar::computeQpProperties(), computeQpProperties(), PorousFlowPermeabilityExponential::computeQpProperties(), PorousFlowPermeabilityKozenyCarman::computeQpProperties(), PorousFlowMassFraction::computeQpProperties(), PorousFlowJoiner::computeQpProperties(), PorousFlowPorosityExponentialBase::computeQpProperties(), and PorousFlowAqueousPreDisChemistry::initQpStatefulProperties().

_porosity_qp

const MaterialProperty<Real>& PorousFlowTotalGravitationalDensityFullySaturatedFromPorosity::_porosity_qp

protected

Porosity at qps.

Definition at line 44 of file PorousFlowTotalGravitationalDensityFullySaturatedFromPorosity.h.

Referenced by computeQpProperties(), and initQpStatefulProperties().

_rho_f_qp

const MaterialProperty<std::vector<Real> >& PorousFlowTotalGravitationalDensityFullySaturatedFromPorosity::_rho_f_qp

protected

Fluid density at qps.

Definition at line 41 of file PorousFlowTotalGravitationalDensityFullySaturatedFromPorosity.h.

Referenced by computeQpProperties(), and initQpStatefulProperties().

_rho_s

const Real PorousFlowTotalGravitationalDensityFullySaturatedFromPorosity::_rho_s

protected

Solid density.

Definition at line 38 of file PorousFlowTotalGravitationalDensityFullySaturatedFromPorosity.h.

Referenced by computeQpProperties(), and initQpStatefulProperties().

---

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

• porous_flow/include/materials/PorousFlowTotalGravitationalDensityFullySaturatedFromPorosity.h
• porous_flow/src/materials/PorousFlowTotalGravitationalDensityFullySaturatedFromPorosity.C