25#ifndef DUMUX_RICHARDS_LENS_SPATIAL_PARAMETERS_HH
26#define DUMUX_RICHARDS_LENS_SPATIAL_PARAMETERS_HH
40template<
class Gr
idGeometry,
class Scalar>
42:
public FVSpatialParams<GridGeometry, Scalar, RichardsLensSpatialParams<GridGeometry, Scalar>>
46 using GridView =
typename GridGeometry::GridView;
47 using FVElementGeometry =
typename GridGeometry::LocalView;
48 using SubControlVolume =
typename FVElementGeometry::SubControlVolume;
49 using Element =
typename GridView::template Codim<0>::Entity;
50 using GlobalPosition =
typename Element::Geometry::GlobalCoordinate;
52 enum { dimWorld = GridView::dimensionworld };
64 lensLowerLeft_ = {1.0, 2.0};
65 lensUpperRight_ = {4.0, 3.0};
68 lensMaterialParams_.setSwr(0.18);
69 lensMaterialParams_.setSnr(0.0);
70 outerMaterialParams_.setSwr(0.05);
71 outerMaterialParams_.setSnr(0.0);
75 lensMaterialParams_.setVgAlpha(0.00045);
76 lensMaterialParams_.setVgn(7.3);
77 outerMaterialParams_.setVgAlpha(0.0037);
78 outerMaterialParams_.setVgn(4.7);
91 if (isInLens_(globalPos))
114 template<
class ElementSolution>
116 const SubControlVolume& scv,
117 const ElementSolution& elemSol)
const
119 const auto& globalPos = scv.dofPosition();
134 if (isInLens_(globalPos))
135 return lensMaterialParams_;
136 return outerMaterialParams_;
140 bool isInLens_(
const GlobalPosition &globalPos)
const
142 for (
int i = 0; i < dimWorld; ++i)
143 if (globalPos[i] < lensLowerLeft_[i] - eps_ || globalPos[i] > lensUpperRight_[i] + eps_)
149 static constexpr Scalar eps_ = 1e-6;
151 GlobalPosition lensLowerLeft_;
152 GlobalPosition lensUpperRight_;
Implementation of the regularized version of the van Genuchten's capillary pressure / relative permea...
The base class for spatial parameters of multi-phase problems using a fully implicit discretization m...
make the local view function available whenever we use the grid geometry
Definition: adapt.hh:29
This material law takes a material law defined for effective saturations and converts it to a materia...
Definition: 2p/efftoabslaw.hh:60
AbsParamsT Params
Definition: 2p/efftoabslaw.hh:64
The base class for spatial parameters of multi-phase problems using a fully implicit discretization m...
Definition: fv.hh:57
const GridGeometry & gridGeometry() const
The finite volume grid geometry.
Definition: fv1p.hh:334
The spatial parameters for the RichardsLensProblem.
Definition: porousmediumflow/richards/implicit/lens/spatialparams.hh:43
RichardsLensSpatialParams(std::shared_ptr< const GridGeometry > gridGeometry)
Definition: porousmediumflow/richards/implicit/lens/spatialparams.hh:60
const MaterialLawParams & materialLawParams(const Element &element, const SubControlVolume &scv, const ElementSolution &elemSol) const
Returns the parameters for the material law for the sub-control volume.
Definition: porousmediumflow/richards/implicit/lens/spatialparams.hh:115
Scalar porosityAtPos(const GlobalPosition &globalPos) const
Returns the porosity [] at a given location.
Definition: porousmediumflow/richards/implicit/lens/spatialparams.hh:101
PermeabilityType permeabilityAtPos(const GlobalPosition &globalPos) const
Returns the intrinsic permeability tensor [m^2] at a given location.
Definition: porousmediumflow/richards/implicit/lens/spatialparams.hh:89
Scalar PermeabilityType
Definition: porousmediumflow/richards/implicit/lens/spatialparams.hh:58
const MaterialLawParams & materialLawParamsAtPos(const GlobalPosition &globalPos) const
Returns the parameters for the material law at a given location.
Definition: porousmediumflow/richards/implicit/lens/spatialparams.hh:132
typename MaterialLaw::Params MaterialLawParams
Definition: porousmediumflow/richards/implicit/lens/spatialparams.hh:56
This model implements a variant of the Richards' equation for quasi-twophase flow.
This material law takes a material law defined for effective saturations and converts it to a materia...