blackoilonephaseindices.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 This file is part of the Open Porous Media project (OPM).
5
6 OPM 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 2 of the License, or
9 (at your option) any later version.
10
11 OPM 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 OPM. If not, see <http://www.gnu.org/licenses/>.
18
19 Consult the COPYING file in the top-level source directory of this
20 module for the precise wording of the license and the list of
21 copyright holders.
22*/
28#ifndef EWOMS_BLACK_OIL_ONE_PHASE_INDICES_HH
29#define EWOMS_BLACK_OIL_ONE_PHASE_INDICES_HH
30
31#include <cassert>
32
33#include <opm/common/utility/ConstexprAssert.hpp>
34
35namespace Opm {
36
42template<unsigned numSolventsV,
43 unsigned numExtbosV,
44 unsigned numPolymersV,
45 unsigned numEnergyV,
46 bool enableFoam,
47 bool enableBrine,
48 unsigned PVOffset,
49 unsigned canonicalCompIdx,
50 unsigned numBioCompV>
52{
54 static constexpr bool oilEnabled = canonicalCompIdx == 0;
55 static constexpr bool waterEnabled = canonicalCompIdx == 1;
56 static constexpr bool gasEnabled = canonicalCompIdx == 2;
57
59 static constexpr bool enableSolvent = numSolventsV > 0;
60
62 static constexpr bool enableExtbo = numExtbosV > 0;
63
65 static constexpr bool enablePolymer = numPolymersV > 0;
66
68 static constexpr bool enableEnergy = numEnergyV > 0;
69
71 static constexpr bool enableMICP = numBioCompV == 5;
72
74 static constexpr bool enableBiofilm = false;
75
77 static constexpr int numSolvents = enableSolvent ? numSolventsV : 0;
78
80 static constexpr int numExtbos = enableExtbo ? numExtbosV : 0;
81
83 static constexpr int numPolymers = enablePolymer ? numPolymersV : 0;
84
86 static constexpr int numEnergy = enableEnergy ? numEnergyV : 0;
87
89 static constexpr int numFoam = enableFoam? 1 : 0;
90
92 static constexpr int numBrine = enableBrine? 1 : 0;
93
95 static constexpr int numPhases = 1;
96
98 static constexpr int numMICPs = enableMICP ? numBioCompV : 0;
99
101 static constexpr int numBioInWat = enableMICP ? 3 : 0;
102
104 static constexpr int numEq = numPhases + numSolvents + numExtbos + numPolymers +
106
108 // Primary variable indices
110
119 static constexpr int waterSwitchIdx = -10000;
120
127 static constexpr int pressureSwitchIdx = PVOffset + 0;
128
135 static constexpr int compositionSwitchIdx = -10000;
136
138 static constexpr int solventSaturationIdx =
139 enableSolvent ? PVOffset + numPhases : -1000;
140
142 static constexpr int zFractionIdx =
143 enableExtbo ? PVOffset + numPhases + numSolvents : -1000;
144
146 static constexpr int polymerConcentrationIdx =
147 enablePolymer ? PVOffset + numPhases + numSolvents : -1000;
148
150 static constexpr int polymerMoleWeightIdx =
151 numPolymers > 1 ? polymerConcentrationIdx + 1 : -1000;
152
154 static constexpr int microbialConcentrationIdx =
155 enableMICP ? PVOffset + numPhases + numSolvents : -1000;
156
158 static constexpr int oxygenConcentrationIdx =
159 numMICPs > 1 ? microbialConcentrationIdx + 1 : -1000;
160
162 static constexpr int ureaConcentrationIdx =
163 numMICPs > 2 ? oxygenConcentrationIdx + 1 : -1000;
164
166 static constexpr int biofilmVolumeFractionIdx =
167 numMICPs > 3 ? ureaConcentrationIdx + 1 : -1000;
168
170 static constexpr int calciteVolumeFractionIdx =
171 numMICPs > 4 ? biofilmVolumeFractionIdx + 1 : -1000;
172
174 static constexpr int foamConcentrationIdx =
175 enableFoam ? PVOffset + numPhases + numSolvents + numPolymers + numMICPs : -1000;
176
178 static constexpr int saltConcentrationIdx =
179 enableBrine ? PVOffset + numPhases + numSolvents + numExtbos + numPolymers + numMICPs + numFoam : -1000;
180
182 static constexpr int temperatureIdx =
184
186 // Equation indices
188
190 static constexpr int conti0EqIdx = PVOffset + 0;
191
193 static constexpr int contiSolventEqIdx =
194 enableSolvent ? PVOffset + numPhases : -1000;
195
197 static constexpr int contiZfracEqIdx =
198 enableExtbo ? PVOffset + numPhases + numSolvents : -1000;
199
201 static constexpr int contiPolymerEqIdx =
202 enablePolymer ? PVOffset + numPhases + numSolvents : -1000;
203
205 static constexpr int contiPolymerMWEqIdx =
206 numPolymers > 1 ? contiPolymerEqIdx + 1 : -1000;
207
209 static constexpr int contiMicrobialEqIdx =
210 enableMICP ? PVOffset + numPhases + numSolvents : -1000;
211
213 static constexpr int contiOxygenEqIdx =
214 numMICPs > 1 ? contiMicrobialEqIdx + 1 : -1000;
215
217 static constexpr int contiUreaEqIdx =
218 numMICPs > 2 ? contiOxygenEqIdx + 1 : -1000;
219
221 static constexpr int contiBiofilmEqIdx =
222 numMICPs > 3 ? contiUreaEqIdx + 1 : -1000;
223
225 static constexpr int contiCalciteEqIdx =
226 numMICPs > 4 ? contiBiofilmEqIdx + 1 : -1000;
227
229 static constexpr int contiFoamEqIdx =
230 enableFoam ? PVOffset + numPhases + numSolvents + numPolymers + numMICPs : -1000;
231
233 static constexpr int contiBrineEqIdx =
234 enableBrine ? PVOffset + numPhases + numSolvents + numExtbos + numPolymers + numMICPs + numFoam : -1000;
235
237 static constexpr int contiEnergyEqIdx =
239};
240
241} // namespace Opm
242
243#endif
Definition: blackoilbioeffectsmodules.hh:43
The primary variable and equation indices for the black-oil model.
Definition: blackoilonephaseindices.hh:52
static constexpr int numBioInWat
Number of biocomponents in the water phase.
Definition: blackoilonephaseindices.hh:101
static constexpr int numBrine
Number of salt equations to be considered.
Definition: blackoilonephaseindices.hh:92
static constexpr int ureaConcentrationIdx
Index of the primary variable for the third MICP component.
Definition: blackoilonephaseindices.hh:162
static constexpr int contiUreaEqIdx
Index of the continuity equation for the third MICP component.
Definition: blackoilonephaseindices.hh:217
static constexpr bool gasEnabled
Definition: blackoilonephaseindices.hh:56
static constexpr int contiEnergyEqIdx
Index of the continuity equation for energy.
Definition: blackoilonephaseindices.hh:237
static constexpr int numPolymers
Number of polymer components to be considered.
Definition: blackoilonephaseindices.hh:83
static constexpr int foamConcentrationIdx
Index of the primary variable for the foam.
Definition: blackoilonephaseindices.hh:174
static constexpr int numEq
The number of equations.
Definition: blackoilonephaseindices.hh:104
static constexpr int polymerConcentrationIdx
Index of the primary variable for the first polymer.
Definition: blackoilonephaseindices.hh:146
static constexpr int contiCalciteEqIdx
Index of the continuity equation for the fifth MICP component.
Definition: blackoilonephaseindices.hh:225
static constexpr int numPhases
The number of fluid phases.
Definition: blackoilonephaseindices.hh:95
static constexpr bool enableSolvent
Are solvents involved?
Definition: blackoilonephaseindices.hh:59
static constexpr int contiZfracEqIdx
Index of the continuity equation for the first extbo component.
Definition: blackoilonephaseindices.hh:197
static constexpr bool enablePolymer
Are polymers involved?
Definition: blackoilonephaseindices.hh:65
static constexpr bool enableExtbo
Is extbo invoked?
Definition: blackoilonephaseindices.hh:62
static constexpr int pressureSwitchIdx
Index of the switching variable which determines the pressure.
Definition: blackoilonephaseindices.hh:127
static constexpr int waterSwitchIdx
Index of the switching variable which determines the composistion of the water phase.
Definition: blackoilonephaseindices.hh:119
static constexpr int compositionSwitchIdx
Index of the switching variable which determines the composition of the hydrocarbon phases.
Definition: blackoilonephaseindices.hh:135
static constexpr int conti0EqIdx
Index of the continuity equation of the first (and only) phase.
Definition: blackoilonephaseindices.hh:190
static constexpr bool enableEnergy
Shall energy be conserved?
Definition: blackoilonephaseindices.hh:68
static constexpr int contiPolymerEqIdx
Index of the continuity equation for the first polymer component.
Definition: blackoilonephaseindices.hh:201
static constexpr int numSolvents
Number of solvent components to be considered.
Definition: blackoilonephaseindices.hh:77
static constexpr int microbialConcentrationIdx
Index of the primary variable for the first MICP component.
Definition: blackoilonephaseindices.hh:154
static constexpr int contiFoamEqIdx
Index of the continuity equation for the foam component.
Definition: blackoilonephaseindices.hh:229
static constexpr int contiBiofilmEqIdx
Index of the continuity equation for the fourth MICP component.
Definition: blackoilonephaseindices.hh:221
static constexpr int temperatureIdx
Index of the primary variable for temperature.
Definition: blackoilonephaseindices.hh:182
static constexpr int calciteVolumeFractionIdx
Index of the primary variable for the fifth MICP component.
Definition: blackoilonephaseindices.hh:170
static constexpr int contiMicrobialEqIdx
Index of the continuity equation for the first MICP component.
Definition: blackoilonephaseindices.hh:209
static constexpr int numMICPs
Number of MICP components to be considered.
Definition: blackoilonephaseindices.hh:98
static constexpr int numFoam
Number of foam equations to be considered.
Definition: blackoilonephaseindices.hh:89
static constexpr int zFractionIdx
Index of the primary variable for the first extbo component.
Definition: blackoilonephaseindices.hh:142
static constexpr int solventSaturationIdx
Index of the primary variable for the first solvent.
Definition: blackoilonephaseindices.hh:138
static constexpr bool enableMICP
Is MICP involved? (microbes, oxygen, urea, biofilm, and calcite)
Definition: blackoilonephaseindices.hh:71
static constexpr bool oilEnabled
Is phase enabled or not.
Definition: blackoilonephaseindices.hh:54
static constexpr int oxygenConcentrationIdx
Index of the primary variable for the second MICP component.
Definition: blackoilonephaseindices.hh:158
static constexpr bool enableBiofilm
Biofilm effects on co2/h2store only for two phase indices.
Definition: blackoilonephaseindices.hh:74
static constexpr int contiPolymerMWEqIdx
Index of the continuity equation for the second polymer component (molecular weight)
Definition: blackoilonephaseindices.hh:205
static constexpr int polymerMoleWeightIdx
Index of the primary variable for the second polymer primary variable (molecular weight)
Definition: blackoilonephaseindices.hh:150
static constexpr int biofilmVolumeFractionIdx
Index of the primary variable for the fourth MICP component.
Definition: blackoilonephaseindices.hh:166
static constexpr int contiOxygenEqIdx
Index of the continuity equation for the second MICP component.
Definition: blackoilonephaseindices.hh:213
static constexpr int saltConcentrationIdx
Index of the primary variable for the salt.
Definition: blackoilonephaseindices.hh:178
static constexpr int contiBrineEqIdx
Index of the continuity equation for the salt component.
Definition: blackoilonephaseindices.hh:233
static constexpr int contiSolventEqIdx
Index of the continuity equation for the first solvent component.
Definition: blackoilonephaseindices.hh:193
static constexpr bool waterEnabled
Definition: blackoilonephaseindices.hh:55
static constexpr int numEnergy
Number of energy equations to be considered.
Definition: blackoilonephaseindices.hh:86
static constexpr int numExtbos
Number of components to be considered for extbo.
Definition: blackoilonephaseindices.hh:80