PorousFlowBasicAdvection Class Reference

Kernel = grad(test) * darcy_velocity * u. More...

#include <PorousFlowBasicAdvection.h>

Inheritance diagram for PorousFlowBasicAdvection:

Public Member Functions
	PorousFlowBasicAdvection (const InputParameters &parameters)

Protected Member Functions
virtual Real	computeQpResidual () override
virtual Real	computeQpJacobian () override
virtual Real	computeQpOffDiagJacobian (unsigned int jvar) override

Protected Attributes
const PorousFlowDictator &	_dictator
	Holds info on the Porous Flow variables. More...
const unsigned	_ph
	Phase of Darcy velocity. More...
const MaterialProperty< std::vector< RealVectorValue > > &	_darcy_velocity
	_darcy_velocity[_qp][ph](j) = j^th component of the Darcy velocity of phase ph More...
const MaterialProperty< std::vector< std::vector< RealVectorValue > > > &	_ddarcy_velocity_dvar
	_ddarcy_velocity_dvar[_qp][ph][v](j) = d(j^th component of the Darcy velocity of phase ph)/d(PorousFlow variable v) More...
const MaterialProperty< std::vector< std::vector< std::vector< RealVectorValue > > > > &	_ddarcy_velocity_dgradvar
	_ddarcy_velocity_dgradvar[_qp][ph][j][v](k) = d(k^th component of the Darcy velocity of phase ph)/d(j^th component of grad(PorousFlow variable v)) More...

Detailed Description

Kernel = grad(test) * darcy_velocity * u.

Definition at line 23 of file PorousFlowBasicAdvection.h.

Constructor & Destructor Documentation

PorousFlowBasicAdvection()

PorousFlowBasicAdvection::PorousFlowBasicAdvection

(

const InputParameters &

parameters

)

Definition at line 27 of file PorousFlowBasicAdvection.C.

  : Kernel(parameters),
    _dictator(getUserObject<PorousFlowDictator>("PorousFlowDictator")),
    _ph(getParam<unsigned int>("phase")),
    _darcy_velocity(
        getMaterialProperty<std::vector<RealVectorValue>>("PorousFlow_darcy_velocity_qp")),
    _ddarcy_velocity_dvar(getMaterialProperty<std::vector<std::vector<RealVectorValue>>>(
        "dPorousFlow_darcy_velocity_qp_dvar")),
    _ddarcy_velocity_dgradvar(
        getMaterialProperty<std::vector<std::vector<std::vector<RealVectorValue>>>>(
            "dPorousFlow_darcy_velocity_qp_dgradvar"))
{
  if (_ph >= _dictator.numPhases())
    paramError("phase",
               "The Dictator proclaims that the maximum phase index in this simulation is ",
               _dictator.numPhases() - 1,
               " whereas you have used ",
               _ph,
               ". Remember that indexing starts at 0. The Dictator is watching you, to "
               "ensure your wellbeing.");
}

Member Function Documentation

computeQpJacobian()

Real PorousFlowBasicAdvection::computeQpJacobian ( )

overrideprotectedvirtual

Definition at line 56 of file PorousFlowBasicAdvection.C.

{
  const Real result = -_grad_test[_i][_qp] * _darcy_velocity[_qp][_ph] * _phi[_j][_qp];
  return result + computeQpOffDiagJacobian(_var.number());
}

computeQpOffDiagJacobian()

Real PorousFlowBasicAdvection::computeQpOffDiagJacobian ( unsigned int  jvar )

overrideprotectedvirtual

Definition at line 63 of file PorousFlowBasicAdvection.C.

{
  if (_dictator.notPorousFlowVariable(jvar))
    return 0.0;
 
  const unsigned pvar = _dictator.porousFlowVariableNum(jvar);
  Real result =
      -_grad_test[_i][_qp] * _ddarcy_velocity_dvar[_qp][_ph][pvar] * _phi[_j][_qp] * _u[_qp];
  for (unsigned j = 0; j < LIBMESH_DIM; ++j)
    result -= _grad_test[_i][_qp] *
              (_ddarcy_velocity_dgradvar[_qp][_ph][j][pvar] * _grad_phi[_j][_qp](j)) * _u[_qp];
 
  return result;
}

Referenced by computeQpJacobian().

computeQpResidual()

Real PorousFlowBasicAdvection::computeQpResidual ( )

overrideprotectedvirtual

Definition at line 50 of file PorousFlowBasicAdvection.C.

{
  return -_grad_test[_i][_qp] * _darcy_velocity[_qp][_ph] * _u[_qp];
}

Member Data Documentation

_darcy_velocity

const MaterialProperty<std::vector<RealVectorValue> >& PorousFlowBasicAdvection::_darcy_velocity

protected

_darcy_velocity[_qp][ph](j) = j^th component of the Darcy velocity of phase ph

Definition at line 40 of file PorousFlowBasicAdvection.h.

Referenced by computeQpJacobian(), and computeQpResidual().

_ddarcy_velocity_dgradvar

const MaterialProperty<std::vector<std::vector<std::vector<RealVectorValue> > > >& PorousFlowBasicAdvection::_ddarcy_velocity_dgradvar

protected

_ddarcy_velocity_dgradvar[_qp][ph][j][v](k) = d(k^th component of the Darcy velocity of phase ph)/d(j^th component of grad(PorousFlow variable v))

Definition at line 54 of file PorousFlowBasicAdvection.h.

Referenced by computeQpOffDiagJacobian().

_ddarcy_velocity_dvar

const MaterialProperty<std::vector<std::vector<RealVectorValue> > >& PorousFlowBasicAdvection::_ddarcy_velocity_dvar

protected

_ddarcy_velocity_dvar[_qp][ph][v](j) = d(j^th component of the Darcy velocity of phase ph)/d(PorousFlow variable v)

Definition at line 46 of file PorousFlowBasicAdvection.h.

Referenced by computeQpOffDiagJacobian().

_dictator

const PorousFlowDictator& PorousFlowBasicAdvection::_dictator

protected

Holds info on the Porous Flow variables.

Definition at line 34 of file PorousFlowBasicAdvection.h.

Referenced by computeQpOffDiagJacobian(), and PorousFlowBasicAdvection().

_ph

const unsigned PorousFlowBasicAdvection::_ph

protected

Phase of Darcy velocity.

Definition at line 37 of file PorousFlowBasicAdvection.h.

Referenced by computeQpJacobian(), computeQpOffDiagJacobian(), computeQpResidual(), and PorousFlowBasicAdvection().

---

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

• porous_flow/include/kernels/PorousFlowBasicAdvection.h
• porous_flow/src/kernels/PorousFlowBasicAdvection.C