Dispersive flux of component k in fluid phase alpha. More...

#include <PorousFlowDispersiveFlux.h>

Inheritance diagram for PorousFlowDispersiveFlux:

Public Member Functions
	PorousFlowDispersiveFlux (const InputParameters &parameters)

Protected Member Functions
virtual Real	computeQpResidual () override

virtual Real	computeQpJacobian () override

virtual Real	computeQpOffDiagJacobian (unsigned int jvar) override

Real	computeQpJac (unsigned int jvar) const
	Derivative of the residual with respect to the PorousFLow Variable with variable number jvar. More...

Protected Attributes
const MaterialProperty< std::vector< Real > > &	_fluid_density_qp
	Fluid density for each phase (at the qp) More...

const MaterialProperty< std::vector< std::vector< Real > > > &	_dfluid_density_qp_dvar
	Derivative of the fluid density for each phase wrt PorousFlow variables (at the qp) More...

const MaterialProperty< std::vector< std::vector< RealGradient > > > &	_grad_mass_frac
	Gradient of mass fraction of each component in each phase. More...

const MaterialProperty< std::vector< std::vector< std::vector< Real > > > > &	_dmass_frac_dvar
	Derivative of mass fraction wrt PorousFlow variables. More...

const MaterialProperty< Real > &	_porosity_qp
	Porosity at the qps. More...

const MaterialProperty< std::vector< Real > > &	_dporosity_qp_dvar
	Derivative of porosity wrt PorousFlow variables (at the qps) More...

const MaterialProperty< std::vector< Real > > &	_tortuosity
	Tortuosity tau_0 * tau_{alpha} for fluid phase alpha. More...

const MaterialProperty< std::vector< std::vector< Real > > > &	_dtortuosity_dvar
	Derivative of tortuosity wrt PorousFlow variables. More...

const MaterialProperty< std::vector< std::vector< Real > > > &	_diffusion_coeff
	Diffusion coefficients of component k in fluid phase alpha. More...

const MaterialProperty< std::vector< std::vector< std::vector< Real > > > > &	_ddiffusion_coeff_dvar
	Derivative of the diffusion coefficients wrt PorousFlow variables. More...

const PorousFlowDictator &	_dictator
	PorousFlowDictator UserObject. More...

const unsigned int	_fluid_component
	Index of the fluid component that this kernel acts on. More...

const unsigned int	_num_phases
	The number of fluid phases. More...

const RankTwoTensor	_identity_tensor
	Identity tensor. More...

const MaterialProperty< std::vector< Real > > &	_relative_permeability
	Relative permeability of each phase. More...

const MaterialProperty< std::vector< std::vector< Real > > > &	_drelative_permeability_dvar
	Derivative of relative permeability wrt PorousFlow variables. More...

const MaterialProperty< std::vector< Real > > &	_fluid_viscosity
	Viscosity of each component in each phase. More...

const MaterialProperty< std::vector< std::vector< Real > > > &	_dfluid_viscosity_dvar
	Derivative of viscosity wrt PorousFlow variables. More...

const MaterialProperty< RealTensorValue > &	_permeability
	Permeability of porous material. More...

const MaterialProperty< std::vector< RealTensorValue > > &	_dpermeability_dvar
	Derivative of permeability wrt PorousFlow variables. More...

const MaterialProperty< std::vector< std::vector< RealTensorValue > > > &	_dpermeability_dgradvar
	d(permeabiity)/d(grad(PorousFlow variable)) More...

const MaterialProperty< std::vector< RealGradient > > &	_grad_p
	Gradient of the pore pressure in each phase. More...

const MaterialProperty< std::vector< std::vector< Real > > > &	_dgrad_p_dgrad_var
	Derivative of Grad porepressure in each phase wrt grad(PorousFlow variables) More...

const MaterialProperty< std::vector< std::vector< RealGradient > > > &	_dgrad_p_dvar
	Derivative of Grad porepressure in each phase wrt PorousFlow variables. More...

const RealVectorValue	_gravity
	Gravitational acceleration. More...

const std::vector< Real >	_disp_long
	Longitudinal dispersivity for each phase. More...

const std::vector< Real >	_disp_trans
	Transverse dispersivity for each phase. More...

const bool	_perm_derivs
	Flag to check whether permeabiity derivatives are non-zero. More...

Detailed Description

Dispersive flux of component k in fluid phase alpha.

Includes the effects of both molecular diffusion and hydrodynamic dispersion.

Definition at line 25 of file PorousFlowDispersiveFlux.h.

Constructor & Destructor Documentation

◆ PorousFlowDispersiveFlux()

PorousFlowDispersiveFlux::PorousFlowDispersiveFlux ( const InputParameters & parameters )

Definition at line 36 of file PorousFlowDispersiveFlux.C.

   : Kernel(parameters),
  
     _fluid_density_qp(getMaterialProperty<std::vector<Real>>("PorousFlow_fluid_phase_density_qp")),
     _dfluid_density_qp_dvar(getMaterialProperty<std::vector<std::vector<Real>>>(
         "dPorousFlow_fluid_phase_density_qp_dvar")),
     _grad_mass_frac(getMaterialProperty<std::vector<std::vector<RealGradient>>>(
         "PorousFlow_grad_mass_frac_qp")),
     _dmass_frac_dvar(getMaterialProperty<std::vector<std::vector<std::vector<Real>>>>(
         "dPorousFlow_mass_frac_qp_dvar")),
     _porosity_qp(getMaterialProperty<Real>("PorousFlow_porosity_qp")),
     _dporosity_qp_dvar(getMaterialProperty<std::vector<Real>>("dPorousFlow_porosity_qp_dvar")),
     _tortuosity(getMaterialProperty<std::vector<Real>>("PorousFlow_tortuosity_qp")),
     _dtortuosity_dvar(
         getMaterialProperty<std::vector<std::vector<Real>>>("dPorousFlow_tortuosity_qp_dvar")),
     _diffusion_coeff(
         getMaterialProperty<std::vector<std::vector<Real>>>("PorousFlow_diffusion_coeff_qp")),
     _ddiffusion_coeff_dvar(getMaterialProperty<std::vector<std::vector<std::vector<Real>>>>(
         "dPorousFlow_diffusion_coeff_qp_dvar")),
     _dictator(getUserObject<PorousFlowDictator>("PorousFlowDictator")),
     _fluid_component(getParam<unsigned int>("fluid_component")),
     _num_phases(_dictator.numPhases()),
     _identity_tensor(RankTwoTensor::initIdentity),
     _relative_permeability(
         getMaterialProperty<std::vector<Real>>("PorousFlow_relative_permeability_qp")),
     _drelative_permeability_dvar(getMaterialProperty<std::vector<std::vector<Real>>>(
         "dPorousFlow_relative_permeability_qp_dvar")),
     _fluid_viscosity(getMaterialProperty<std::vector<Real>>("PorousFlow_viscosity_qp")),
     _dfluid_viscosity_dvar(
         getMaterialProperty<std::vector<std::vector<Real>>>("dPorousFlow_viscosity_qp_dvar")),
     _permeability(getMaterialProperty<RealTensorValue>("PorousFlow_permeability_qp")),
     _dpermeability_dvar(
         getMaterialProperty<std::vector<RealTensorValue>>("dPorousFlow_permeability_qp_dvar")),
     _dpermeability_dgradvar(getMaterialProperty<std::vector<std::vector<RealTensorValue>>>(
         "dPorousFlow_permeability_qp_dgradvar")),
     _grad_p(getMaterialProperty<std::vector<RealGradient>>("PorousFlow_grad_porepressure_qp")),
     _dgrad_p_dgrad_var(getMaterialProperty<std::vector<std::vector<Real>>>(
         "dPorousFlow_grad_porepressure_qp_dgradvar")),
     _dgrad_p_dvar(getMaterialProperty<std::vector<std::vector<RealGradient>>>(
         "dPorousFlow_grad_porepressure_qp_dvar")),
     _gravity(getParam<RealVectorValue>("gravity")),
     _disp_long(getParam<std::vector<Real>>("disp_long")),
     _disp_trans(getParam<std::vector<Real>>("disp_trans")),
     _perm_derivs(_dictator.usePermDerivs())
 {
   // Check that sufficient values of the dispersion coefficients have been entered
   if (_disp_long.size() != _num_phases)
     paramError(
         "disp_long",
         "The number of longitudinal dispersion coefficients is not equal to the number of phases");
  
   if (_disp_trans.size() != _num_phases)
     paramError("disp_trans",
                "The number of transverse dispersion coefficients disp_trans in is not equal to the "
                "number of phases");
 }

Member Function Documentation

◆ computeQpJac()

Real PorousFlowDispersiveFlux::computeQpJac ( unsigned int jvar ) const

protected

Derivative of the residual with respect to the PorousFLow Variable with variable number jvar.

This is used by both computeQpJacobian and computeQpOffDiagJacobian

Parameters

jvar	take the derivative wrt this variable number

Returns: dResidual_dVar

Definition at line 143 of file PorousFlowDispersiveFlux.C.

 {
   // If the variable is not a valid PorousFlow variable, set the Jacobian to 0
   if (_dictator.notPorousFlowVariable(jvar))
     return 0.0;
  
   const unsigned int pvar = _dictator.porousFlowVariableNum(jvar);
  
   RealVectorValue velocity;
   Real velocity_abs;
   RankTwoTensor v2;
   RankTwoTensor dispersion;
   dispersion.zero();
   Real diffusion;
   RealVectorValue flux = 0.0;
   RealVectorValue dflux = 0.0;
  
   for (unsigned int ph = 0; ph < _num_phases; ++ph)
   {
     // Diffusive component
     diffusion =
         _porosity_qp[_qp] * _tortuosity[_qp][ph] * _diffusion_coeff[_qp][ph][_fluid_component];
  
     // Calculate Darcy velocity
     velocity = (_permeability[_qp] * (_grad_p[_qp][ph] - _fluid_density_qp[_qp][ph] * _gravity) *
                 _relative_permeability[_qp][ph] / _fluid_viscosity[_qp][ph]);
     velocity_abs = std::sqrt(velocity * velocity);
  
     if (velocity_abs > 0.0)
     {
       v2.vectorOuterProduct(velocity, velocity);
  
       // Add longitudinal dispersion to diffusive component
       diffusion += _disp_trans[ph] * velocity_abs;
       dispersion = (_disp_long[ph] - _disp_trans[ph]) * v2 / velocity_abs;
     }
  
     // Derivative of Darcy velocity
     RealVectorValue dvelocity =
         _permeability[_qp] * (_grad_phi[_j][_qp] * _dgrad_p_dgrad_var[_qp][ph][pvar] -
                               _phi[_j][_qp] * _dfluid_density_qp_dvar[_qp][ph][pvar] * _gravity);
     dvelocity += _permeability[_qp] * (_dgrad_p_dvar[_qp][ph][pvar] * _phi[_j][_qp]);
  
     if (_perm_derivs)
     {
       dvelocity += _dpermeability_dvar[_qp][pvar] * _phi[_j][_qp] *
                    (_grad_p[_qp][ph] - _fluid_density_qp[_qp][ph] * _gravity);
  
       for (unsigned int i = 0; i < LIBMESH_DIM; ++i)
         dvelocity += _dpermeability_dgradvar[_qp][i][pvar] * _grad_phi[_j][_qp](i) *
                      (_grad_p[_qp][ph] - _fluid_density_qp[_qp][ph] * _gravity);
     }
  
     dvelocity = dvelocity * _relative_permeability[_qp][ph] / _fluid_viscosity[_qp][ph] +
                 (_permeability[_qp] * (_grad_p[_qp][ph] - _fluid_density_qp[_qp][ph] * _gravity)) *
                     (_drelative_permeability_dvar[_qp][ph][pvar] / _fluid_viscosity[_qp][ph] -
                      _relative_permeability[_qp][ph] * _dfluid_viscosity_dvar[_qp][ph][pvar] /
                          std::pow(_fluid_viscosity[_qp][ph], 2)) *
                     _phi[_j][_qp];
  
     Real dvelocity_abs = 0.0;
     if (velocity_abs > 0.0)
       dvelocity_abs = velocity * dvelocity / velocity_abs;
  
     // Derivative of diffusion term (note: dispersivity is assumed constant)
     Real ddiffusion = _phi[_j][_qp] * _dporosity_qp_dvar[_qp][pvar] * _tortuosity[_qp][ph] *
                       _diffusion_coeff[_qp][ph][_fluid_component];
     ddiffusion += _phi[_j][_qp] * _porosity_qp[_qp] * _dtortuosity_dvar[_qp][ph][pvar] *
                   _diffusion_coeff[_qp][ph][_fluid_component];
     ddiffusion += _phi[_j][_qp] * _porosity_qp[_qp] * _tortuosity[_qp][ph] *
                   _ddiffusion_coeff_dvar[_qp][ph][_fluid_component][pvar];
     ddiffusion += _disp_trans[ph] * dvelocity_abs;
  
     // Derivative of dispersion term (note: dispersivity is assumed constant)
     RankTwoTensor ddispersion;
     ddispersion.zero();
     if (velocity_abs > 0.0)
     {
       RankTwoTensor dv2a, dv2b;
       dv2a.vectorOuterProduct(velocity, dvelocity);
       dv2b.vectorOuterProduct(dvelocity, velocity);
       ddispersion = (_disp_long[ph] - _disp_trans[ph]) * (dv2a + dv2b) / velocity_abs;
       ddispersion -=
           (_disp_long[ph] - _disp_trans[ph]) * v2 * dvelocity_abs / velocity_abs / velocity_abs;
     }
  
     dflux += _phi[_j][_qp] * _dfluid_density_qp_dvar[_qp][ph][pvar] *
              (diffusion * _identity_tensor + dispersion) *
              _grad_mass_frac[_qp][ph][_fluid_component];
     dflux += _fluid_density_qp[_qp][ph] * (ddiffusion * _identity_tensor + ddispersion) *
              _grad_mass_frac[_qp][ph][_fluid_component];
  
     // NOTE: Here we assume that d(grad_mass_frac)/d(var) = d(mass_frac)/d(var) * grad_phi
     //       This is true for most PorousFlow scenarios, but not for chemical reactions
     //       where mass_frac is a nonlinear function of the primary MOOSE Variables
     dflux += _fluid_density_qp[_qp][ph] * (diffusion * _identity_tensor + dispersion) *
              _dmass_frac_dvar[_qp][ph][_fluid_component][pvar] * _grad_phi[_j][_qp];
   }
  
   return _grad_test[_i][_qp] * dflux;
 }

Referenced by computeQpJacobian(), and computeQpOffDiagJacobian().

◆ computeQpJacobian()

Real PorousFlowDispersiveFlux::computeQpJacobian ( )

overrideprotectedvirtual

Definition at line 131 of file PorousFlowDispersiveFlux.C.

 {
   return computeQpJac(_var.number());
 }

◆ computeQpOffDiagJacobian()

Real PorousFlowDispersiveFlux::computeQpOffDiagJacobian ( unsigned int jvar )

overrideprotectedvirtual

Definition at line 137 of file PorousFlowDispersiveFlux.C.

 {
   return computeQpJac(jvar);
 }

◆ computeQpResidual()

Real PorousFlowDispersiveFlux::computeQpResidual ( )

overrideprotectedvirtual

Definition at line 94 of file PorousFlowDispersiveFlux.C.

 {
   RealVectorValue flux = 0.0;
   RealVectorValue velocity;
   Real velocity_abs;
   RankTwoTensor v2;
   RankTwoTensor dispersion;
   dispersion.zero();
   Real diffusion;
  
   for (unsigned int ph = 0; ph < _num_phases; ++ph)
   {
     // Diffusive component
     diffusion =
         _porosity_qp[_qp] * _tortuosity[_qp][ph] * _diffusion_coeff[_qp][ph][_fluid_component];
  
     // Calculate Darcy velocity
     velocity = (_permeability[_qp] * (_grad_p[_qp][ph] - _fluid_density_qp[_qp][ph] * _gravity) *
                 _relative_permeability[_qp][ph] / _fluid_viscosity[_qp][ph]);
     velocity_abs = std::sqrt(velocity * velocity);
  
     if (velocity_abs > 0.0)
     {
       v2.vectorOuterProduct(velocity, velocity);
  
       // Add longitudinal dispersion to diffusive component
       diffusion += _disp_trans[ph] * velocity_abs;
       dispersion = (_disp_long[ph] - _disp_trans[ph]) * v2 / velocity_abs;
     }
  
     flux += _fluid_density_qp[_qp][ph] * (diffusion * _identity_tensor + dispersion) *
             _grad_mass_frac[_qp][ph][_fluid_component];
   }
   return _grad_test[_i][_qp] * flux;
 }

Member Data Documentation

◆ _ddiffusion_coeff_dvar

const MaterialProperty<std::vector<std::vector<std::vector<Real> > > >& PorousFlowDispersiveFlux::_ddiffusion_coeff_dvar

protected

Derivative of the diffusion coefficients wrt PorousFlow variables.

Definition at line 74 of file PorousFlowDispersiveFlux.h.

Referenced by computeQpJac().

◆ _dfluid_density_qp_dvar

const MaterialProperty<std::vector<std::vector<Real> > >& PorousFlowDispersiveFlux::_dfluid_density_qp_dvar

protected

Derivative of the fluid density for each phase wrt PorousFlow variables (at the qp)

Definition at line 50 of file PorousFlowDispersiveFlux.h.

Referenced by computeQpJac().

◆ _dfluid_viscosity_dvar

const MaterialProperty<std::vector<std::vector<Real> > >& PorousFlowDispersiveFlux::_dfluid_viscosity_dvar

protected

Derivative of viscosity wrt PorousFlow variables.

Definition at line 98 of file PorousFlowDispersiveFlux.h.

Referenced by computeQpJac().

◆ _dgrad_p_dgrad_var

const MaterialProperty<std::vector<std::vector<Real> > >& PorousFlowDispersiveFlux::_dgrad_p_dgrad_var

protected

Derivative of Grad porepressure in each phase wrt grad(PorousFlow variables)

Definition at line 113 of file PorousFlowDispersiveFlux.h.

Referenced by computeQpJac().

◆ _dgrad_p_dvar

const MaterialProperty<std::vector<std::vector<RealGradient> > >& PorousFlowDispersiveFlux::_dgrad_p_dvar

protected

Derivative of Grad porepressure in each phase wrt PorousFlow variables.

Definition at line 116 of file PorousFlowDispersiveFlux.h.

Referenced by computeQpJac().

◆ _dictator

const PorousFlowDictator& PorousFlowDispersiveFlux::_dictator

protected

PorousFlowDictator UserObject.

Definition at line 77 of file PorousFlowDispersiveFlux.h.

Referenced by computeQpJac().

◆ _diffusion_coeff

const MaterialProperty<std::vector<std::vector<Real> > >& PorousFlowDispersiveFlux::_diffusion_coeff

protected

Diffusion coefficients of component k in fluid phase alpha.

Definition at line 71 of file PorousFlowDispersiveFlux.h.

Referenced by computeQpJac(), and computeQpResidual().

◆ _disp_long

const std::vector<Real> PorousFlowDispersiveFlux::_disp_long

protected

Longitudinal dispersivity for each phase.

Definition at line 122 of file PorousFlowDispersiveFlux.h.

Referenced by computeQpJac(), computeQpResidual(), and PorousFlowDispersiveFlux().

◆ _disp_trans

const std::vector<Real> PorousFlowDispersiveFlux::_disp_trans

protected

Transverse dispersivity for each phase.

Definition at line 125 of file PorousFlowDispersiveFlux.h.

Referenced by computeQpJac(), computeQpResidual(), and PorousFlowDispersiveFlux().

◆ _dmass_frac_dvar

const MaterialProperty<std::vector<std::vector<std::vector<Real> > > >& PorousFlowDispersiveFlux::_dmass_frac_dvar

protected

Derivative of mass fraction wrt PorousFlow variables.

Definition at line 56 of file PorousFlowDispersiveFlux.h.

Referenced by computeQpJac().

◆ _dpermeability_dgradvar

const MaterialProperty<std::vector<std::vector<RealTensorValue> > >& PorousFlowDispersiveFlux::_dpermeability_dgradvar

protected

d(permeabiity)/d(grad(PorousFlow variable))

Definition at line 107 of file PorousFlowDispersiveFlux.h.

Referenced by computeQpJac().

◆ _dpermeability_dvar

const MaterialProperty<std::vector<RealTensorValue> >& PorousFlowDispersiveFlux::_dpermeability_dvar

protected

Derivative of permeability wrt PorousFlow variables.

Definition at line 104 of file PorousFlowDispersiveFlux.h.

Referenced by computeQpJac().

◆ _dporosity_qp_dvar

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

protected

Derivative of porosity wrt PorousFlow variables (at the qps)

Definition at line 62 of file PorousFlowDispersiveFlux.h.

Referenced by computeQpJac().

◆ _drelative_permeability_dvar

const MaterialProperty<std::vector<std::vector<Real> > >& PorousFlowDispersiveFlux::_drelative_permeability_dvar

protected

Derivative of relative permeability wrt PorousFlow variables.

Definition at line 92 of file PorousFlowDispersiveFlux.h.

Referenced by computeQpJac().

◆ _dtortuosity_dvar

const MaterialProperty<std::vector<std::vector<Real> > >& PorousFlowDispersiveFlux::_dtortuosity_dvar

protected

Derivative of tortuosity wrt PorousFlow variables.

Definition at line 68 of file PorousFlowDispersiveFlux.h.

Referenced by computeQpJac().

◆ _fluid_component

const unsigned int PorousFlowDispersiveFlux::_fluid_component

protected

Index of the fluid component that this kernel acts on.

Definition at line 80 of file PorousFlowDispersiveFlux.h.

Referenced by computeQpJac(), and computeQpResidual().

◆ _fluid_density_qp

const MaterialProperty<std::vector<Real> >& PorousFlowDispersiveFlux::_fluid_density_qp

protected

Fluid density for each phase (at the qp)

Definition at line 47 of file PorousFlowDispersiveFlux.h.

Referenced by computeQpJac(), and computeQpResidual().

◆ _fluid_viscosity

const MaterialProperty<std::vector<Real> >& PorousFlowDispersiveFlux::_fluid_viscosity

protected

Viscosity of each component in each phase.

Definition at line 95 of file PorousFlowDispersiveFlux.h.

Referenced by computeQpJac(), and computeQpResidual().

◆ _grad_mass_frac

const MaterialProperty<std::vector<std::vector<RealGradient> > >& PorousFlowDispersiveFlux::_grad_mass_frac

protected

Gradient of mass fraction of each component in each phase.

Definition at line 53 of file PorousFlowDispersiveFlux.h.

Referenced by computeQpJac(), and computeQpResidual().

◆ _grad_p

const MaterialProperty<std::vector<RealGradient> >& PorousFlowDispersiveFlux::_grad_p

protected

Gradient of the pore pressure in each phase.

Definition at line 110 of file PorousFlowDispersiveFlux.h.

Referenced by computeQpJac(), and computeQpResidual().

◆ _gravity

const RealVectorValue PorousFlowDispersiveFlux::_gravity

protected

Gravitational acceleration.

Definition at line 119 of file PorousFlowDispersiveFlux.h.

Referenced by computeQpJac(), and computeQpResidual().

◆ _identity_tensor

const RankTwoTensor PorousFlowDispersiveFlux::_identity_tensor

protected

Identity tensor.

Definition at line 86 of file PorousFlowDispersiveFlux.h.

Referenced by computeQpJac(), and computeQpResidual().

◆ _num_phases

const unsigned int PorousFlowDispersiveFlux::_num_phases

protected

The number of fluid phases.

Definition at line 83 of file PorousFlowDispersiveFlux.h.

Referenced by computeQpJac(), computeQpResidual(), and PorousFlowDispersiveFlux().

◆ _perm_derivs

const bool PorousFlowDispersiveFlux::_perm_derivs

protected

Flag to check whether permeabiity derivatives are non-zero.

Definition at line 128 of file PorousFlowDispersiveFlux.h.

Referenced by computeQpJac().

◆ _permeability

const MaterialProperty<RealTensorValue>& PorousFlowDispersiveFlux::_permeability

protected

Permeability of porous material.

Definition at line 101 of file PorousFlowDispersiveFlux.h.

Referenced by computeQpJac(), and computeQpResidual().

◆ _porosity_qp

const MaterialProperty<Real>& PorousFlowDispersiveFlux::_porosity_qp

protected

Porosity at the qps.

Definition at line 59 of file PorousFlowDispersiveFlux.h.

Referenced by computeQpJac(), and computeQpResidual().

◆ _relative_permeability

const MaterialProperty<std::vector<Real> >& PorousFlowDispersiveFlux::_relative_permeability

protected

Relative permeability of each phase.

Definition at line 89 of file PorousFlowDispersiveFlux.h.

Referenced by computeQpJac(), and computeQpResidual().

◆ _tortuosity

const MaterialProperty<std::vector<Real> >& PorousFlowDispersiveFlux::_tortuosity

protected

Tortuosity tau_0 * tau_{alpha} for fluid phase alpha.

Definition at line 65 of file PorousFlowDispersiveFlux.h.

Referenced by computeQpJac(), and computeQpResidual().

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

porous_flow/include/kernels/PorousFlowDispersiveFlux.h
porous_flow/src/kernels/PorousFlowDispersiveFlux.C

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ PorousFlowDispersiveFlux()

Member Function Documentation

◆ computeQpJac()

◆ computeQpJacobian()

◆ computeQpOffDiagJacobian()

◆ computeQpResidual()

Member Data Documentation

◆ _ddiffusion_coeff_dvar

◆ _dfluid_density_qp_dvar

◆ _dfluid_viscosity_dvar

◆ _dgrad_p_dgrad_var

◆ _dgrad_p_dvar

◆ _dictator

◆ _diffusion_coeff

◆ _disp_long

◆ _disp_trans

◆ _dmass_frac_dvar

◆ _dpermeability_dgradvar

◆ _dpermeability_dvar

◆ _dporosity_qp_dvar

◆ _drelative_permeability_dvar

◆ _dtortuosity_dvar

◆ _fluid_component

◆ _fluid_density_qp

◆ _fluid_viscosity

◆ _grad_mass_frac

◆ _grad_p

◆ _gravity

◆ _identity_tensor

◆ _num_phases

◆ _perm_derivs

◆ _permeability

◆ _porosity_qp

◆ _relative_permeability

◆ _tortuosity