3.1-git
DUNE for Multi-{Phase, Component, Scale, Physics, ...} flow and transport in porous media
gravity/problem_bulk.hh
Go to the documentation of this file.
1// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
2// vi: set et ts=4 sw=4 sts=4:
3/*****************************************************************************
4 * See the file COPYING for full copying permissions. *
5 * *
6 * This program is free software: you can redistribute it and/or modify *
7 * it under the terms of the GNU General Public License as published by *
8 * the Free Software Foundation, either version 3 of the License, or *
9 * (at your option) any later version. *
10 * *
11 * This program is distributed in the hope that it will be useful, *
12 * but WITHOUT ANY WARRANTY; without even the implied warranty of *
13 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
14 * GNU General Public License for more details. *
15 * *
16 * You should have received a copy of the GNU General Public License *
17 * along with this program. If not, see <http://www.gnu.org/licenses/>. *
18 *****************************************************************************/
24#ifndef DUMUX_TEST_TPFAFACETCOUPLING_ONEP_BULKPROBLEM_HH
25#define DUMUX_TEST_TPFAFACETCOUPLING_ONEP_BULKPROBLEM_HH
26
27#include <dune/alugrid/grid.hh>
28
31
35
38
39#include "spatialparams.hh"
40
41// default for the bulk grid type
42#ifndef BULKGRIDTYPE
43#define BULKGRIDTYPE Dune::ALUGrid<2, 2, Dune::cube, Dune::nonconforming>
44#endif
45
46namespace Dumux {
47// forward declarations
48template<class TypeTag> class OnePBulkProblem;
49
50namespace Properties {
51// create the type tag nodes
52namespace TTag {
53struct OnePBulk { using InheritsFrom = std::tuple<OneP>; };
54struct OnePBulkTpfa { using InheritsFrom = std::tuple<CCTpfaFacetCouplingModel, OnePBulk>; };
55struct OnePBulkMpfa { using InheritsFrom = std::tuple<CCMpfaFacetCouplingModel, OnePBulk>; };
56} // end namespace TTag
57
58// Set the grid type
59template<class TypeTag>
60struct Grid<TypeTag, TTag::OnePBulk> { using type = BULKGRIDTYPE; };
61// Set the problem type
62template<class TypeTag>
63struct Problem<TypeTag, TTag::OnePBulk> { using type = OnePBulkProblem<TypeTag>; };
64// set the spatial params
65template<class TypeTag>
66struct SpatialParams<TypeTag, TTag::OnePBulk>
67{
70};
71
72// the fluid system
73template<class TypeTag>
74struct FluidSystem<TypeTag, TTag::OnePBulk>
75{
76private:
78public:
80};
81
82} // end namespace Properties
83
89template<class TypeTag>
90class OnePBulkProblem : public PorousMediumFlowProblem<TypeTag>
91{
92 using ParentType = PorousMediumFlowProblem<TypeTag>;
93
95 using PrimaryVariables = typename GridVariables::PrimaryVariables;
96 using Scalar = typename GridVariables::Scalar;
97
98 using GridGeometry = typename GridVariables::GridGeometry;
99 using FVElementGeometry = typename GridGeometry::LocalView;
100 using SubControlVolume = typename GridGeometry::SubControlVolume;
101 using SubControlVolumeFace = typename GridGeometry::SubControlVolumeFace;
102 using GridView = typename GridGeometry::GridView;
103 using Element = typename GridView::template Codim<0>::Entity;
104 using GlobalPosition = typename Element::Geometry::GlobalCoordinate;
105
109
110 static constexpr int dimWorld = GridView::dimensionworld;
111
112public:
113 OnePBulkProblem(std::shared_ptr<const GridGeometry> gridGeometry,
114 std::shared_ptr<typename ParentType::SpatialParams> spatialParams,
115 std::shared_ptr<CouplingManager> couplingManager,
116 const std::string& paramGroup = "")
118 , couplingManagerPtr_(couplingManager)
119 {
120 problemName_ = getParam<std::string>("Vtk.OutputName") + "_" +
121 getParamFromGroup<std::string>(this->paramGroup(), "Problem.Name");
122 }
123
125 const std::string& name() const
126 { return problemName_; }
127
129 BoundaryTypes boundaryTypesAtPos(const GlobalPosition& globalPos) const
130 {
131 BoundaryTypes values;
132 values.setAllNeumann();
133 if (globalPos[1] < 1e-6 || globalPos[1] > this->gridGeometry().bBoxMax()[1] - 1e-6)
134 values.setAllDirichlet();
135 return values;
136 }
137
139 BoundaryTypes interiorBoundaryTypes(const Element& element, const SubControlVolumeFace& scvf) const
140 {
141 BoundaryTypes values;
142 values.setAllNeumann();
143 return values;
144 }
145
147 NumEqVector sourceAtPos(const GlobalPosition& globalPos) const
148 { return NumEqVector(0.0); }
149
151 PrimaryVariables dirichletAtPos(const GlobalPosition& globalPos) const
152 { return initialAtPos(globalPos); }
153
155 NumEqVector neumannAtPos(const GlobalPosition& globalPos) const
156 { return NumEqVector(0.0); }
157
159 PrimaryVariables initialAtPos(const GlobalPosition& globalPos) const
160 {
161 const auto g = this->spatialParams().gravity(globalPos)[dimWorld-1];
162 const auto h = 1.0 - (3.0-globalPos[dimWorld-1])*g;
163 return PrimaryVariables(h);
164 }
165
167 Scalar temperature() const
168 { return 283.15; /*10°*/ }
169
171 const CouplingManager& couplingManager() const
172 { return *couplingManagerPtr_; }
173
174private:
175 std::shared_ptr<CouplingManager> couplingManagerPtr_;
176 std::string problemName_;
177};
178
179} // end namespace Dumux
180
181#endif
Setting constant fluid properties via the input file.
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
Class to specify the type of a boundary.
Definition: common/boundarytypes.hh:38
Base class for all finite-volume problems.
Definition: common/fvproblem.hh:50
const std::string & paramGroup() const
The parameter group in which to retrieve runtime parameters.
Definition: common/fvproblem.hh:592
const GridGeometry & gridGeometry() const
The finite volume grid geometry.
Definition: common/fvproblem.hh:588
Property to specify the type of scalar values.
Definition: common/properties.hh:53
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 of the spatial parameters object.
Definition: common/properties.hh:221
The type of the fluid system to use.
Definition: common/properties.hh:223
A liquid phase consisting of a single component.
Definition: 1pliquid.hh:46
Definition: multidomain/couplingmanager.hh:46
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
The spatial parameters class for the test problem using the 1p cc model.
Definition: multidomain/boundary/stokesdarcy/1p2c_1p2c/spatialparams.hh:41
Test problem for the incompressible one-phase model with coupling across the bulk grid facets.
Definition: problem_1p_bulk.hh:90
PrimaryVariables dirichletAtPos(const GlobalPosition &globalPos) const
Evaluates the Dirichlet boundary condition for a given position.
Definition: gravity/problem_bulk.hh:151
NumEqVector sourceAtPos(const GlobalPosition &globalPos) const
Evaluates the source term at a given position.
Definition: gravity/problem_bulk.hh:147
const std::string & name() const
The problem name.
Definition: gravity/problem_bulk.hh:125
BoundaryTypes boundaryTypesAtPos(const GlobalPosition &globalPos) const
Specifies the type of boundary condition at a given position.
Definition: gravity/problem_bulk.hh:129
BoundaryTypes interiorBoundaryTypes(const Element &element, const SubControlVolumeFace &scvf) const
Specifies the type of interior boundary condition at a given position.
Definition: gravity/problem_bulk.hh:139
NumEqVector neumannAtPos(const GlobalPosition &globalPos) const
Evaluates the Neumann boundary condition for a boundary segment.
Definition: gravity/problem_bulk.hh:155
OnePBulkProblem(std::shared_ptr< const GridGeometry > gridGeometry, std::shared_ptr< typename ParentType::SpatialParams > spatialParams, std::shared_ptr< CouplingManager > couplingManager, const std::string &paramGroup="")
Definition: gravity/problem_bulk.hh:113
Scalar temperature() const
Returns the temperature in in the domain.
Definition: gravity/problem_bulk.hh:167
const CouplingManager & couplingManager() const
Returns reference to the coupling manager.
Definition: analytical/problem_bulk.hh:207
PrimaryVariables initialAtPos(const GlobalPosition &globalPos) const
Evaluates the initial conditions.
Definition: analytical/problem_bulk.hh:199
Definition: analytical/problem_bulk.hh:49
std::tuple< OneP > InheritsFrom
Definition: analytical/problem_bulk.hh:49
Definition: analytical/problem_bulk.hh:50
std::tuple< CCTpfaFacetCouplingModel, OnePBulk > InheritsFrom
Definition: analytical/problem_bulk.hh:50
Dune::ALUGrid< 2, 2, Dune::cube, Dune::nonconforming > type
Definition: analytical/problem_bulk.hh:56
Definition: gravity/problem_bulk.hh:55
std::tuple< CCMpfaFacetCouplingModel, OnePBulk > InheritsFrom
Definition: gravity/problem_bulk.hh:55
Properties (and default properties) for all models using the box scheme together with coupling across...
Properties (and default properties) for all models using cell-centered finite volume scheme with MPFA...
Properties (and default properties) for all models using cell-centered finite volume scheme with TPFA...
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.
#define BULKGRIDTYPE
Definition: gravity/problem_bulk.hh:43