13#ifndef DUMUX_POROMECHANICS_PRORELASTIC_LOCAL_RESIDUAL_HH 
   14#define DUMUX_POROMECHANICS_PRORELASTIC_LOCAL_RESIDUAL_HH 
   27template<
class TypeTag>
 
   33    using Element = 
typename GridView::template Codim<0>::Entity;
 
   38    using SubControlVolume = 
typename FVElementGeometry::SubControlVolume;
 
   39    using SubControlVolumeFace = 
typename FVElementGeometry::SubControlVolumeFace;
 
   43    using VolumeVariables = 
typename ElementVolumeVariables::VolumeVariables;
 
   46    using ParentType::ParentType;
 
   62                              const Element& element,
 
   63                              const FVElementGeometry& fvGeometry,
 
   64                              const ElementVolumeVariables& elemVolVars,
 
   65                              const SubControlVolume &scv)
 const 
   67        NumEqVector source(0.0);
 
   70        source += problem.source(element, fvGeometry, elemVolVars, scv);
 
   73        source += problem.scvPointSources(element, fvGeometry, elemVolVars, scv);
 
   80            const auto& vv = elemVolVars[scv];
 
   81            const auto phi = vv.porosity();
 
   82            const auto rhoFluid = problem.spatialParams().effectiveFluidDensity(element, scv);
 
   83            const auto rhoAverage = phi*rhoFluid + (1.0 - phi)*vv.solidDensity();
 
   86            const auto& g = problem.spatialParams().gravity(scv.center());
 
   87            for (
int dir = 0; dir < GridView::dimensionworld; ++dir)
 
   88                source[ Indices::momentum(dir) ] += rhoAverage*g[dir];
 
 
 
Element-wise calculation of the local residual for problems using the elastic model considering linea...
Definition solidmechanics/elastic/localresidual.hh:32
Element-wise calculation of the local residual for problems using the poroelastic model.
Definition poromechanics/poroelastic/localresidual.hh:29
NumEqVector computeSource(const Problem &problem, const Element &element, const FVElementGeometry &fvGeometry, const ElementVolumeVariables &elemVolVars, const SubControlVolume &scv) const
Calculate the source term of the equation.
Definition poromechanics/poroelastic/localresidual.hh:61
Defines all properties used in Dumux.
typename NumEqVectorTraits< PrimaryVariables >::type NumEqVector
A vector with the same size as numbers of equations This is the default implementation and has to be ...
Definition numeqvector.hh:34
T getParamFromGroup(Args &&... args)
A free function to get a parameter from the parameter tree singleton with a model group.
Definition parameters.hh:149
typename GetProp< TypeTag, Property >::type GetPropType
get the type alias defined in the property
Definition propertysystem.hh:296
A helper to deduce a vector with the same size as numbers of equations.
Element-wise calculation of the local residual for problems using the elastic model considering linea...