27#ifndef DUMUX_RICHARDS_LENSPROBLEM_HH
28#define DUMUX_RICHARDS_LENSPROBLEM_HH
30#include <dune/grid/yaspgrid.hh>
32#include <dune/alugrid/grid.hh>
35#include <dune/grid/uggrid.hh>
50template <
class TypeTag>
51class RichardsLensProblem;
64template<
class TypeTag>
65struct Grid<TypeTag, TTag::RichardsLens> {
using type = Dune::YaspGrid<2>; };
68template<
class TypeTag>
73template<
class TypeTag>
77template<
class TypeTag>
111template <
class TypeTag>
124 pressureIdx = Indices::pressureIdx,
125 conti0EqIdx = Indices::conti0EqIdx,
126 bothPhases = Indices::bothPhases,
129 dimWorld = GridView::dimensionworld
131 using Element =
typename GridView::template Codim<0>::Entity;
133 using GlobalPosition =
typename Element::Geometry::GlobalCoordinate;
139 name_ = getParam<std::string>(
"Problem.Name");
152 const std::string&
name()
const
161 {
return 273.15 + 10; };
187 BoundaryTypes bcTypes;
188 if (onLeftBoundary_(globalPos) || onRightBoundary_(globalPos))
189 bcTypes.setAllDirichlet();
191 bcTypes.setAllNeumann();
204 return initial_(globalPos);
217 NumEqVector values(0.0);
218 if (onInlet_(globalPos))
219 values[conti0EqIdx] = -0.04;
237 {
return initial_(globalPos); };
242 PrimaryVariables initial_(
const GlobalPosition &globalPos)
const
244 PrimaryVariables values(0.0);
245 const Scalar sw = 0.0;
246 using MaterialLaw =
typename ParentType::SpatialParams::MaterialLaw;
247 const Scalar pc = MaterialLaw::pc(this->
spatialParams().materialLawParamsAtPos(globalPos), sw);
249 values.setState(bothPhases);
253 bool onLeftBoundary_(
const GlobalPosition &globalPos)
const
255 return globalPos[0] < this->
gridGeometry().bBoxMin()[0] + eps_;
258 bool onRightBoundary_(
const GlobalPosition &globalPos)
const
260 return globalPos[0] > this->
gridGeometry().bBoxMax()[0] - eps_;
263 bool onLowerBoundary_(
const GlobalPosition &globalPos)
const
265 return globalPos[1] < this->
gridGeometry().bBoxMin()[1] + eps_;
268 bool onUpperBoundary_(
const GlobalPosition &globalPos)
const
270 return globalPos[1] > this->
gridGeometry().bBoxMax()[1] - eps_;
273 bool onInlet_(
const GlobalPosition &globalPos)
const
276 Scalar lambda = (this->
gridGeometry().bBoxMax()[0] - globalPos[0])/width;
277 return onUpperBoundary_(globalPos) && 0.5 < lambda + eps_ && lambda < 2.0/3.0 + eps_;
280 static constexpr Scalar eps_ = 1.5e-7;
282 GlobalPosition lensLowerLeft_;
283 GlobalPosition lensUpperRight_;
Defines a type tag and some properties for models using the box scheme.
Properties for all models using cell-centered finite volume scheme with TPFA.
A much simpler (and thus potentially less buggy) version of pure water.
A liquid phase consisting of a single component.
make the local view function available whenever we use the grid geometry
Definition: adapt.hh:29
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:149
Base class for all finite-volume problems.
Definition: common/fvproblem.hh:50
const GridGeometry & gridGeometry() const
The finite volume grid geometry.
Definition: common/fvproblem.hh:588
The DUNE grid type.
Definition: common/properties.hh:57
UndefinedProperty type
Definition: common/properties.hh:57
Property to specify the type of a problem which has to be solved.
Definition: common/properties.hh:69
The type of the spatial parameters object.
Definition: common/properties.hh:221
Base class for all fully implicit porous media problems.
Definition: dumux/porousmediumflow/problem.hh:39
SpatialParams & spatialParams()
Returns the spatial parameters object.
Definition: dumux/porousmediumflow/problem.hh:146
A water infiltration problem with a low-permeability lens embedded into a high-permeability domain wh...
Definition: test/porousmediumflow/richards/implicit/lens/problem.hh:113
Scalar nonWettingReferencePressure() const
Returns the reference pressure [Pa] of the non-wetting fluid phase within a finite volume.
Definition: test/porousmediumflow/richards/implicit/lens/problem.hh:169
BoundaryTypes boundaryTypesAtPos(const GlobalPosition &globalPos) const
Specifies which kind of boundary condition should be used for which equation on a given boundary segm...
Definition: test/porousmediumflow/richards/implicit/lens/problem.hh:185
NumEqVector neumannAtPos(const GlobalPosition &globalPos) const
Evaluates the boundary conditions for a Neumann boundary segment.
Definition: test/porousmediumflow/richards/implicit/lens/problem.hh:215
PrimaryVariables initialAtPos(const GlobalPosition &globalPos) const
Evaluates the initial values for a control volume.
Definition: test/porousmediumflow/richards/implicit/lens/problem.hh:236
Scalar temperature() const
Returns the temperature [K] within a finite volume.
Definition: test/porousmediumflow/richards/implicit/lens/problem.hh:160
PrimaryVariables dirichletAtPos(const GlobalPosition &globalPos) const
Evaluates the boundary conditions for a Dirichlet boundary segment.
Definition: test/porousmediumflow/richards/implicit/lens/problem.hh:202
const std::string & name() const
The problem name.
Definition: test/porousmediumflow/richards/implicit/lens/problem.hh:152
RichardsLensProblem(std::shared_ptr< const GridGeometry > gridGeometry)
Definition: test/porousmediumflow/richards/implicit/lens/problem.hh:136
Definition: test/porousmediumflow/richards/implicit/lens/problem.hh:57
std::tuple< Richards > InheritsFrom
Definition: test/porousmediumflow/richards/implicit/lens/problem.hh:57
Definition: test/porousmediumflow/richards/implicit/lens/problem.hh:58
std::tuple< RichardsLens, BoxModel > InheritsFrom
Definition: test/porousmediumflow/richards/implicit/lens/problem.hh:58
Definition: test/porousmediumflow/richards/implicit/lens/problem.hh:59
std::tuple< RichardsLens, CCTpfaModel > InheritsFrom
Definition: test/porousmediumflow/richards/implicit/lens/problem.hh:59
Dune::YaspGrid< 2 > type
Definition: test/porousmediumflow/richards/implicit/lens/problem.hh:65
The spatial parameters for the RichardsLensProblem.
Definition: porousmediumflow/richards/implicit/lens/spatialparams.hh:43
This model implements a variant of the Richards' equation for quasi-twophase flow.
Base class for all porous media problems.
#define GRIDTYPE
Definition: test/porousmediumflow/1p/implicit/incompressible/problem.hh:48
Definition of the spatial parameters for the MaxwellStefan problem.