28#ifndef EWOMS_RICHARDS_LOCAL_RESIDUAL_HH
29#define EWOMS_RICHARDS_LOCAL_RESIDUAL_HH
41template <
class TypeTag>
51 enum { contiEqIdx = Indices::contiEqIdx };
60 template <
class LhsEval>
62 const ElementContext& elemCtx,
66 const IntensiveQuantities& intQuants = elemCtx.intensiveQuantities(dofIdx,
timeIdx);
70 Toolbox::template
decay<LhsEval>(intQuants.fluidState().density(liquidPhaseIdx))
71 *Toolbox::template
decay<LhsEval>(intQuants.fluidState().saturation(liquidPhaseIdx))
79 const ElementContext& elemCtx,
86 unsigned upIdx =
static_cast<unsigned>(
extQuants.upstreamIndex(liquidPhaseIdx));
88 const IntensiveQuantities&
up = elemCtx.intensiveQuantities(
upIdx,
timeIdx);
92 const Evaluation& rho =
up.fluidState().density(liquidPhaseIdx);
96 flux[contiEqIdx] =
extQuants.volumeFlux(liquidPhaseIdx)*Toolbox::value(rho);
103 const ElementContext& elemCtx,
106 { elemCtx.problem().source(source, elemCtx, dofIdx,
timeIdx); }
Element-wise calculation of the residual for the Richards model.
Definition richardslocalresidual.hh:43
void computeFlux(RateVector &flux, const ElementContext &elemCtx, unsigned scvfIdx, unsigned timeIdx) const
Evaluates the total mass flux of all conservation quantities over a face of a sub-control volume.
Definition richardslocalresidual.hh:78
void computeStorage(Dune::FieldVector< LhsEval, numEq > &storage, const ElementContext &elemCtx, unsigned dofIdx, unsigned timeIdx) const
Evaluate the amount all conservation quantities (e.g.
Definition richardslocalresidual.hh:61
void computeSource(RateVector &source, const ElementContext &elemCtx, unsigned dofIdx, unsigned timeIdx) const
Calculate the source term of the equation.
Definition richardslocalresidual.hh:102
This file contains a set of helper functions used by VFPProd / VFPInj.
Definition blackoilboundaryratevector.hh:37
constexpr auto getPropValue()
get the value data member of a property
Definition propertysystem.hh:242
typename Properties::Detail::GetPropImpl< TypeTag, Property >::type::type GetPropType
get the type alias defined in the property (equivalent to old macro GET_PROP_TYPE(....
Definition propertysystem.hh:235
Calculates and stores the data which is required to calculate the flux of fluid over a face of a fini...
Intensive quantities required by the Richards model.