26#ifndef DUMUX_INJECTION_SPATIAL_PARAMS_HH
27#define DUMUX_INJECTION_SPATIAL_PARAMS_HH
41template<
class Gr
idGeometry,
class Scalar>
44 InjectionSpatialParams<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;
52 static constexpr int dimWorld = GridView::dimensionworld;
56 using GlobalPosition =
typename Element::Geometry::GlobalCoordinate;
76 coarsePorosity_ = 0.3;
79 fineMaterialParams_.setSwr(0.2);
80 fineMaterialParams_.setSnr(0.0);
81 coarseMaterialParams_.setSwr(0.2);
82 coarseMaterialParams_.setSnr(0.0);
85 fineMaterialParams_.setPe(1e4);
86 coarseMaterialParams_.setPe(1e4);
87 fineMaterialParams_.setLambda(2.0);
88 coarseMaterialParams_.setLambda(2.0);
98 if (isFineMaterial_(globalPos))
110 if (isFineMaterial_(globalPos))
111 return finePorosity_;
112 return coarsePorosity_;
124 if (isFineMaterial_(globalPos))
125 return fineMaterialParams_;
126 return coarseMaterialParams_;
135 template<
class Flu
idSystem>
137 {
return FluidSystem::H2OIdx; }
140 bool isFineMaterial_(
const GlobalPosition &globalPos)
const
141 {
return globalPos[dimWorld-1] > layerBottom_; }
147 Scalar finePorosity_;
148 Scalar coarsePorosity_;
Implementation of a regularized version of the Brooks-Corey capillary pressure / relative permeabilit...
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
Implementation of the regularized Brooks-Corey capillary pressure / relative permeability <-> saturat...
Definition: regularizedbrookscorey.hh:62
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
Definition of the spatial parameters for the injection problem which uses the isothermal two-phase tw...
Definition: porousmediumflow/2p2c/implicit/injection/spatialparams.hh:45
int wettingPhaseAtPos(const GlobalPosition &globalPos) const
Function for defining which phase is to be considered as the wetting phase.
Definition: porousmediumflow/2p2c/implicit/injection/spatialparams.hh:136
InjectionSpatialParams(std::shared_ptr< const GridGeometry > gridGeometry)
Definition: porousmediumflow/2p2c/implicit/injection/spatialparams.hh:65
Scalar PermeabilityType
Export the type used for the permeability.
Definition: porousmediumflow/2p2c/implicit/injection/spatialparams.hh:60
const MaterialLawParams & materialLawParamsAtPos(const GlobalPosition &globalPos) const
Returns the parameter object for the capillary-pressure/ saturation material law.
Definition: porousmediumflow/2p2c/implicit/injection/spatialparams.hh:122
Scalar porosityAtPos(const GlobalPosition &globalPos) const
Returns the porosity .
Definition: porousmediumflow/2p2c/implicit/injection/spatialparams.hh:108
typename MaterialLaw::Params MaterialLawParams
Definition: porousmediumflow/2p2c/implicit/injection/spatialparams.hh:63
PermeabilityType permeabilityAtPos(const GlobalPosition &globalPos) const
Returns the intrinsic permeability tensor .
Definition: porousmediumflow/2p2c/implicit/injection/spatialparams.hh:96
Defines a type tag and some properties for models using the box scheme.
This material law takes a material law defined for effective saturations and converts it to a materia...