3.4
DUNE for Multi-{Phase, Component, Scale, Physics, ...} flow and transport in porous media
porousmediumflow/2p/iofields.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 *****************************************************************************/
25#ifndef DUMUX_TWOP_IO_FIELDS_HH
26#define DUMUX_TWOP_IO_FIELDS_HH
27
28#include <dumux/io/name.hh>
29
30#include "model.hh"
31
32namespace Dumux {
33
39{
40public:
41 template <class OutputModule>
42 static void initOutputModule(OutputModule& out)
43 {
44 using VolumeVariables = typename OutputModule::VolumeVariables;
45 using FS = typename VolumeVariables::FluidSystem;
46
47 for (int phaseIdx = 0; phaseIdx < FS::numPhases; ++phaseIdx)
48 {
49 out.addVolumeVariable([phaseIdx](const VolumeVariables& v){ return v.saturation(phaseIdx); },
50 IOName::saturation<FS>(phaseIdx));
51 out.addVolumeVariable([phaseIdx](const VolumeVariables& v){ return v.pressure(phaseIdx); },
52 IOName::pressure<FS>(phaseIdx));
53 out.addVolumeVariable([phaseIdx](const auto& v){ return v.density(phaseIdx); },
54 IOName::density<FS>(phaseIdx));
55 out.addVolumeVariable([phaseIdx](const auto& v){ return v.mobility(phaseIdx); },
56 IOName::mobility<FS>(phaseIdx));
57 }
58
59 out.addVolumeVariable([](const auto& v){ return v.capillaryPressure(); },
61 out.addVolumeVariable([](const auto& v){ return v.porosity(); },
63 }
64
65 template <class ModelTraits, class FluidSystem, class SolidSystem = void>
66 static std::string primaryVariableName(int pvIdx, int state = 0)
67 {
68 if (ModelTraits::priVarFormulation() == TwoPFormulation::p0s1)
69 return pvIdx == 0 ? IOName::pressure<FluidSystem>(FluidSystem::phase0Idx)
70 : IOName::saturation<FluidSystem>(FluidSystem::phase1Idx);
71 else
72 return pvIdx == 0 ? IOName::pressure<FluidSystem>(FluidSystem::phase1Idx)
73 : IOName::saturation<FluidSystem>(FluidSystem::phase0Idx);
74 }
75};
76
77} // end namespace Dumux
78
79#endif
A collection of input/output field names for common physical quantities.
@ p0s1
first phase pressure and second phase saturation as primary variables
Definition: adapt.hh:29
std::string capillaryPressure() noexcept
I/O name of capillary pressure.
Definition: name.hh:135
std::string porosity() noexcept
I/O name of porosity.
Definition: name.hh:139
Adds I/O fields specific to the two-phase model.
Definition: porousmediumflow/2p/iofields.hh:39
static std::string primaryVariableName(int pvIdx, int state=0)
Definition: porousmediumflow/2p/iofields.hh:66
static void initOutputModule(OutputModule &out)
Definition: porousmediumflow/2p/iofields.hh:42
Adaption of the fully implicit scheme to the tracer transport model.