28#ifndef OPM_FLOW_GENERIC_PROBLEM_HPP
29#define OPM_FLOW_GENERIC_PROBLEM_HPP
31#include <opm/material/common/UniformXTabulated2DFunction.hpp>
32#include <opm/material/common/Tabulated1DFunction.hpp>
51template<
typename Gr
id,
typename Gr
idView>
class LookUpData;
59template<
class Gr
idView,
class Flu
idSystem>
63 using Scalar =
typename FluidSystem::Scalar;
68 Scalar referencePressure;
69 Scalar compressibility;
74 const GridView& gridView);
79 const GridView& gridView);
98 { briefDescription_ =
msg; }
253 int numPressurePointsEquil()
const
254 {
return numPressurePointsEquil_; }
256 bool operator==(
const FlowGenericProblem& rhs)
const;
258 template<
class Serializer>
272 void initFluidSystem_();
289 bool beginEpisode_(
bool enableExperiments,
291 void beginTimeStep_(
bool enableExperiments,
299 void readRockParameters_(
const std::vector<Scalar>& cellCenterDepths,
300 std::function<std::array<int,3>(
const unsigned)>
ijkIndex);
301 void readRockCompactionParameters_();
303 void readBlackoilExtentionsInitialConditions_(std::size_t numDof,
306 bool enablePolymerMolarWeight,
309 void updatePvtnum_();
310 void updateSatnum_();
311 void updateMiscnum_();
312 void updatePlmixnum_();
314 const EclipseState& eclState_;
316 const GridView& gridView_;
318 static inline std::string briefDescription_;
319 std::array<std::vector<Scalar>, 2> referencePorosity_;
321 std::vector<int> pvtnum_;
322 std::vector<unsigned short> satnum_;
323 std::vector<unsigned short> miscnum_;
324 std::vector<unsigned short> plmixnum_;
326 std::vector<RockParams> rockParams_;
327 std::vector<unsigned short> rockTableIdx_;
328 std::vector<TabulatedTwoDFunction> rockCompPoroMultWc_;
329 std::vector<TabulatedTwoDFunction> rockCompTransMultWc_;
330 std::vector<TabulatedFunction> rockCompPoroMult_;
331 std::vector<TabulatedFunction> rockCompTransMult_;
332 std::vector<Scalar> rockCompTransMultVal_;
335 std::vector<Scalar> maxOilSaturation_;
336 std::vector<Scalar> maxWaterSaturation_;
337 std::vector<Scalar> minRefPressure_;
338 std::vector<Scalar> overburdenPressure_;
339 std::vector<Scalar> solventSaturation_;
340 std::vector<Scalar> solventRsw_;
345 Scalar initialTimeStepSize_;
346 Scalar maxTimeStepAfterWellEvent_;
349 int numPressurePointsEquil_;
351 bool enableDriftCompensation_;
352 bool explicitRockCompaction_;
355 using Grid = std::remove_cv_t<
typename std::remove_reference<
decltype(gridView_.grid())>::type>;
363 std::function<
unsigned(
unsigned)> lookupIdxOnLevelZeroAssigner_()
370 void updateNum(
const std::string& name, std::vector<T>&
numbers, std::size_t
num_regions);
This problem simulates an input file given in the data format used by the commercial ECLiPSE simulato...
This problem simulates an input file given in the data format used by the commercial ECLiPSE simulato...
Definition FlowGenericProblem.hpp:61
void setPorosity(Scalar poro, unsigned elementIdx, unsigned timeIdx=0)
Sets the porosity of an element.
Definition FlowGenericProblem.hpp:155
Scalar maxPolymerAdsorption(unsigned elemIdx) const
Returns the max polymer adsorption value.
Definition FlowGenericProblem_impl.hpp:731
unsigned pvtRegionIndex(unsigned elemIdx) const
Returns the index the relevant PVT region given a cell index.
Definition FlowGenericProblem_impl.hpp:690
Scalar oxygenConcentration(unsigned elemIdx) const
Returns the initial oxygen concentration for a given a cell index.
Definition FlowGenericProblem_impl.hpp:643
Scalar microbialConcentration(unsigned elemIdx) const
Returns the initial microbial concentration for a given a cell index.
Definition FlowGenericProblem_impl.hpp:631
static std::string briefDescription()
Returns a human readable description of the problem for the help message.
Definition FlowGenericProblem_impl.hpp:126
Scalar solventRsw(unsigned elemIdx) const
Returns the initial solvent dissolved in water for a given a cell index.
Definition FlowGenericProblem_impl.hpp:594
Scalar overburdenPressure(unsigned elementIdx) const
Get the pressure of the overburden.
Definition FlowGenericProblem_impl.hpp:572
Scalar porosity(unsigned globalSpaceIdx, unsigned timeIdx) const
Direct indexed access to the porosity.
Definition FlowGenericProblem_impl.hpp:327
Scalar rockCompressibility(unsigned globalSpaceIdx) const
Direct access to rock compressibility.
Definition FlowGenericProblem_impl.hpp:312
unsigned miscnumRegionIndex(unsigned elemIdx) const
Returns the index the relevant MISC region given a cell index.
Definition FlowGenericProblem_impl.hpp:710
Scalar referencePorosity(unsigned elementIdx, unsigned timeIdx) const
Returns the porosity of an element.
Definition FlowGenericProblem.hpp:137
unsigned satnumRegionIndex(unsigned elemIdx) const
Returns the index the relevant saturation function region given a cell index.
Definition FlowGenericProblem_impl.hpp:700
Scalar rockFraction(unsigned elementIdx, unsigned timeIdx) const
Returns the rockFraction of an element.
Definition FlowGenericProblem_impl.hpp:335
Scalar maxWaterSaturation(unsigned globalDofIdx) const
Returns an element's historic maximum water phase saturation that was observed during the simulation.
Definition FlowGenericProblem_impl.hpp:550
unsigned plmixnumRegionIndex(unsigned elemIdx) const
Returns the index the relevant PLMIXNUM (for polymer module) region given a cell index.
Definition FlowGenericProblem_impl.hpp:720
Scalar solventSaturation(unsigned elemIdx) const
Returns the initial solvent saturation for a given a cell index.
Definition FlowGenericProblem_impl.hpp:583
Scalar polymerMolecularWeight(const unsigned elemIdx) const
Returns the polymer molecule weight for a given cell index.
Definition FlowGenericProblem_impl.hpp:619
Scalar biofilmConcentration(unsigned elemIdx) const
Returns the initial biofilm concentration for a given a cell index.
Definition FlowGenericProblem_impl.hpp:667
Scalar minOilPressure(unsigned globalDofIdx) const
Returns an element's historic minimum pressure of the oil phase that was observed during the simulati...
Definition FlowGenericProblem_impl.hpp:561
Scalar calciteConcentration(unsigned elemIdx) const
Returns the initial calcite concentration for a given a cell index.
Definition FlowGenericProblem_impl.hpp:679
Scalar ureaConcentration(unsigned elemIdx) const
Returns the initial urea concentration for a given a cell index.
Definition FlowGenericProblem_impl.hpp:655
bool shouldWriteOutput() const
Always returns true.
Definition FlowGenericProblem.hpp:277
static void setBriefDescription(const std::string &msg)
Specifies the string returned by briefDescription()
Definition FlowGenericProblem.hpp:97
static std::string helpPreamble(int, const char **argv)
Returns the string that is printed before the list of command line parameters in the help message.
Definition FlowGenericProblem_impl.hpp:111
bool shouldWriteRestartFile() const
Returns true if an eWoms restart file should be written to disk.
Definition FlowGenericProblem.hpp:286
Scalar polymerConcentration(unsigned elemIdx) const
Returns the initial polymer concentration for a given a cell index.
Definition FlowGenericProblem_impl.hpp:607
Definition FlowGenericProblem.hpp:51
This file contains a set of helper functions used by VFPProd / VFPInj.
Definition blackoilboundaryratevector.hh:37
constexpr auto getPropValue()
get the value data member of a property
Definition propertysystem.hh:242
Definition FlowGenericProblem.hpp:67
Struct holding MICP extension data.
Definition SolutionContainers.hpp:57
Struct holding polymer extension data.
Definition SolutionContainers.hpp:37