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>
86 unsigned upIdx =
static_cast<unsigned>(
extQuants.upstreamIndex(liquidPhaseIdx));
92 const Evaluation& rho =
up.fluidState().density(liquidPhaseIdx);
96 flux[contiEqIdx] =
extQuants.volumeFlux(liquidPhaseIdx)*Toolbox::value(rho);
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.