26#ifndef DUMUX_1PTEST_PROBLEM_HH
27#define DUMUX_1PTEST_PROBLEM_HH
29#include <dune/grid/yaspgrid.hh>
47template <
class TypeTag>
61template<
class TypeTag>
69template<
class TypeTag>
70struct Grid<TypeTag, TTag::OnePTest>
71{
using type = Dune::YaspGrid<2>; };
74template<
class TypeTag>
79template<
class TypeTag>
89template<
class TypeTag>
117template <
class TypeTag>
118class OnePTestProblem :
public PorousMediumFlowProblem<TypeTag>
120 using ParentType = PorousMediumFlowProblem<TypeTag>;
121 using GridView = GetPropType<TypeTag, Properties::GridView>;
122 using Element =
typename GridView::template Codim<0>::Entity;
123 using Scalar = GetPropType<TypeTag, Properties::Scalar>;
130 pressureIdx = Indices::pressureIdx
133 using PrimaryVariables = GetPropType<TypeTag, Properties::PrimaryVariables>;
134 using BoundaryTypes = GetPropType<TypeTag, Properties::BoundaryTypes>;
135 using GridGeometry = GetPropType<TypeTag, Properties::GridGeometry>;
137 static constexpr int dimWorld = GridView::dimensionworld;
139 using GlobalPosition =
typename Element::Geometry::GlobalCoordinate;
145 name_ = getParam<std::string>(
"Problem.Name");
169 {
return 273.15 + 10; }
185 BoundaryTypes values;
187 if (globalPos[dimWorld-1] < eps_ || globalPos[dimWorld-1] > this->
gridGeometry().bBoxMax()[dimWorld-1] - eps_)
188 values.setAllDirichlet();
190 values.setAllNeumann();
204 PrimaryVariables values(0);
205 values[pressureIdx] = 1.0e+5*(2.0 - globalPos[dimWorld-1]);
224 PrimaryVariables priVars(0.0);
225 priVars[pressureIdx] = 1.0e+5;
233 static constexpr Scalar eps_ = 1.0e-6;
Defines a type tag and some properties for models using the box scheme.
Properties for all models using cell-centered finite volume scheme with mpfa.
Properties for all models using cell-centered finite volume scheme with TPFA.
Forchheimer's law specialized for different discretization schemes This file contains the data which ...
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
Property tag Indices
Definition: porousmediumflow/sequential/properties.hh:59
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
Property to specify the type of a problem which has to be solved.
Definition: common/properties.hh:69
The type for the calculation the advective fluxes.
Definition: common/properties.hh:208
UndefinedProperty type
Definition: common/properties.hh:208
The type of the spatial parameters object.
Definition: common/properties.hh:221
The type of the fluid system to use.
Definition: common/properties.hh:223
forward declare
Definition: forchheimerslaw.hh:38
A liquid phase consisting of a single component.
Definition: 1pliquid.hh:46
Definition: test_propertysystem.cc:60
std::tuple< OnePModel, BoxDisc > InheritsFrom
Definition: test_propertysystem.cc:60
PrimaryVariables dirichletAtPos(const GlobalPosition &globalPos) const
Evaluates the boundary conditions for a Dirichlet control volume.
Definition: test/porousmediumflow/1p/implicit/isothermal/problem.hh:202
OnePTestProblem(std::shared_ptr< const GridGeometry > gridGeometry)
Definition: test/porousmediumflow/1p/implicit/isothermal/problem.hh:142
PrimaryVariables initialAtPos(const GlobalPosition &globalPos) const
Evaluates the initial value for a control volume.
Definition: test/porousmediumflow/1p/implicit/isothermal/problem.hh:222
BoundaryTypes boundaryTypesAtPos(const GlobalPosition &globalPos) const
Specifies which kind of boundary condition should be used for which equation on a given boundary cont...
Definition: test/porousmediumflow/1p/implicit/isothermal/problem.hh:183
Scalar temperature() const
Returns the temperature within the domain in [K].
Definition: test/porousmediumflow/1p/implicit/isothermal/problem.hh:168
std::string name() const
The problem name.
Definition: test/porousmediumflow/1p/implicit/isothermal/problem.hh:158
The spatial parameters class for the test problem using the incompressible 1p model.
Definition: multidomain/boundary/darcydarcy/1p_1p/spatialparams.hh:62
Definition: test/porousmediumflow/1p/implicit/isothermal/problem.hh:55
std::tuple< BoxModel, OnePTest > InheritsFrom
Definition: test/porousmediumflow/1p/implicit/isothermal/problem.hh:55
Definition: test/porousmediumflow/1p/implicit/isothermal/problem.hh:56
std::tuple< CCTpfaModel, OnePTest > InheritsFrom
Definition: test/porousmediumflow/1p/implicit/isothermal/problem.hh:56
Definition: test/porousmediumflow/1p/implicit/isothermal/problem.hh:57
std::tuple< CCMpfaModel, OnePTest > InheritsFrom
Definition: test/porousmediumflow/1p/implicit/isothermal/problem.hh:57
GetPropType< TypeTag, Scalar > Scalar
Definition: test/porousmediumflow/1p/implicit/isothermal/problem.hh:64
Dune::YaspGrid< 2 > type
Definition: test/porousmediumflow/1p/implicit/isothermal/problem.hh:71
GetPropType< TypeTag, Scalar > Scalar
Definition: test/porousmediumflow/1p/implicit/isothermal/problem.hh:83
GetPropType< TypeTag, GridGeometry > GridGeometry
Definition: test/porousmediumflow/1p/implicit/isothermal/problem.hh:82
A single-phase, isothermal flow model using the fully implicit scheme.
Base class for all porous media problems.
Definition of the spatial parameters for the MaxwellStefan problem.