mirror of
https://github.com/OPM/opm-simulators.git
synced 2025-02-25 18:55:30 -06:00
rearranged to make other grids easier to integrate with or without fem
This commit is contained in:
hnil
Arne Morten Kvarving
parent
66ff026008
commit
030c06afb2
@ -20,874 +20,10 @@
|
|||||||
module for the precise wording of the license and the list of
|
module for the precise wording of the license and the list of
|
||||||
copyright holders.
|
copyright holders.
|
||||||
*/
|
*/
|
||||||
|
|
||||||
#include <config.h>
|
#include <config.h>
|
||||||
#include <ebos/eclgenericproblem.hh>
|
#include <ebos/eclgenericproblem.hh>
|
||||||
|
#include <ebos/eclgenericproblem_impl.hh>
|
||||||
|
|
||||||
#include <opm/material/fluidsystems/BlackOilFluidSystem.hpp>
|
template class Opm::EclGenericProblem<Dune::GridView<Dune::DefaultLeafGridViewTraits<Dune::CpGrid>>,
|
||||||
|
Opm::BlackOilFluidSystem<double,Opm::BlackOilDefaultIndexTraits>,
|
||||||
#include <opm/input/eclipse/EclipseState/EclipseState.hpp>
|
double>;
|
||||||
#include <opm/input/eclipse/EclipseState/Tables/OverburdTable.hpp>
|
|
||||||
#include <opm/input/eclipse/EclipseState/Tables/RockwnodTable.hpp>
|
|
||||||
|
|
||||||
#include <opm/input/eclipse/Schedule/Schedule.hpp>
|
|
||||||
|
|
||||||
#include <opm/input/eclipse/Units/Units.hpp>
|
|
||||||
|
|
||||||
#include <opm/grid/CpGrid.hpp>
|
|
||||||
#include <opm/grid/polyhedralgrid.hh>
|
|
||||||
|
|
||||||
#include <dune/common/parametertree.hh>
|
|
||||||
|
|
||||||
#if HAVE_DUNE_ALUGRID
|
|
||||||
#include <dune/alugrid/grid.hh>
|
|
||||||
#include <dune/alugrid/3d/gridview.hh>
|
|
||||||
#include "alucartesianindexmapper.hh"
|
|
||||||
#endif // HAVE_DUNE_ALUGRID
|
|
||||||
|
|
||||||
#if HAVE_DUNE_FEM
|
|
||||||
#include <dune/fem/gridpart/adaptiveleafgridpart.hh>
|
|
||||||
#include <dune/fem/gridpart/common/gridpart2gridview.hh>
|
|
||||||
#include <ebos/femcpgridcompat.hh>
|
|
||||||
#endif // HAVE_DUNE_FEM
|
|
||||||
|
|
||||||
#include <boost/date_time.hpp>
|
|
||||||
|
|
||||||
#include <algorithm>
|
|
||||||
#include <limits>
|
|
||||||
#include <stdexcept>
|
|
||||||
#include <iostream>
|
|
||||||
|
|
||||||
namespace Opm {
|
|
||||||
|
|
||||||
int eclPositionalParameter(Dune::ParameterTree& tree,
|
|
||||||
std::set<std::string>& seenParams,
|
|
||||||
std::string& errorMsg,
|
|
||||||
const char** argv,
|
|
||||||
int paramIdx)
|
|
||||||
{
|
|
||||||
std::string param = argv[paramIdx];
|
|
||||||
size_t i = param.find('=');
|
|
||||||
if (i != std::string::npos) {
|
|
||||||
std::string oldParamName = param.substr(0, i);
|
|
||||||
std::string oldParamValue = param.substr(i+1);
|
|
||||||
std::string newParamName = "--" + oldParamName;
|
|
||||||
std::replace(newParamName.begin(),
|
|
||||||
newParamName.end(), '_' , '-');
|
|
||||||
errorMsg =
|
|
||||||
"The old syntax to specify parameters on the command line is no longer supported: "
|
|
||||||
"Try replacing '" + oldParamName + "=" + oldParamValue + "' with "+
|
|
||||||
"'" + newParamName + "=" + oldParamValue + "'!";
|
|
||||||
return 0;
|
|
||||||
}
|
|
||||||
|
|
||||||
if (seenParams.count("EclDeckFileName") > 0) {
|
|
||||||
errorMsg =
|
|
||||||
"Parameter 'EclDeckFileName' specified multiple times"
|
|
||||||
" as a command line parameter";
|
|
||||||
return 0;
|
|
||||||
}
|
|
||||||
|
|
||||||
tree["EclDeckFileName"] = argv[paramIdx];
|
|
||||||
seenParams.insert("EclDeckFileName");
|
|
||||||
return 1;
|
|
||||||
}
|
|
||||||
|
|
||||||
template<class GridView, class FluidSystem, class Scalar>
|
|
||||||
EclGenericProblem<GridView,FluidSystem,Scalar>::
|
|
||||||
EclGenericProblem(const EclipseState& eclState,
|
|
||||||
const Schedule& schedule,
|
|
||||||
const GridView& gridView)
|
|
||||||
: eclState_(eclState)
|
|
||||||
, schedule_(schedule)
|
|
||||||
, gridView_(gridView)
|
|
||||||
{
|
|
||||||
}
|
|
||||||
|
|
||||||
template<class GridView, class FluidSystem, class Scalar>
|
|
||||||
EclGenericProblem<GridView,FluidSystem,Scalar>
|
|
||||||
EclGenericProblem<GridView,FluidSystem,Scalar>::
|
|
||||||
serializationTestObject(const EclipseState& eclState,
|
|
||||||
const Schedule& schedule,
|
|
||||||
const GridView& gridView)
|
|
||||||
{
|
|
||||||
EclGenericProblem result(eclState, schedule, gridView);
|
|
||||||
result.maxOilSaturation_ = {1.0, 2.0};
|
|
||||||
result.maxWaterSaturation_ = {6.0};
|
|
||||||
result.minOilPressure_ = {7.0, 8.0, 9.0, 10.0};
|
|
||||||
result.overburdenPressure_ = {11.0};
|
|
||||||
result.solventSaturation_ = {15.0};
|
|
||||||
result.polymer_ = PolymerSolutionContainer<Scalar>::serializationTestObject();
|
|
||||||
result.micp_ = MICPSolutionContainer<Scalar>::serializationTestObject();
|
|
||||||
result.lastRv_ = {21.0};
|
|
||||||
result.maxDRv_ = {22.0, 23.0};
|
|
||||||
result.convectiveDrs_ = {24.0, 25.0, 26.0};
|
|
||||||
result.lastRs_ = {27.0};
|
|
||||||
result.maxDRs_ = {28.0};
|
|
||||||
result.dRsDtOnlyFreeGas_ = {false, true};
|
|
||||||
|
|
||||||
return result;
|
|
||||||
}
|
|
||||||
|
|
||||||
template<class GridView, class FluidSystem, class Scalar>
|
|
||||||
std::string
|
|
||||||
EclGenericProblem<GridView,FluidSystem,Scalar>::
|
|
||||||
helpPreamble(int,
|
|
||||||
const char **argv)
|
|
||||||
{
|
|
||||||
std::string desc = EclGenericProblem::briefDescription();
|
|
||||||
if (!desc.empty())
|
|
||||||
desc = desc + "\n";
|
|
||||||
|
|
||||||
return
|
|
||||||
"Usage: "+std::string(argv[0]) + " [OPTIONS] [ECL_DECK_FILENAME]\n"
|
|
||||||
+ desc;
|
|
||||||
}
|
|
||||||
|
|
||||||
template<class GridView, class FluidSystem, class Scalar>
|
|
||||||
std::string
|
|
||||||
EclGenericProblem<GridView,FluidSystem,Scalar>::
|
|
||||||
briefDescription()
|
|
||||||
{
|
|
||||||
if (briefDescription_.empty())
|
|
||||||
return
|
|
||||||
"The Ecl-deck Black-Oil reservoir Simulator (ebos); a hydrocarbon "
|
|
||||||
"reservoir simulation program that processes ECL-formatted input "
|
|
||||||
"files that is part of the Open Porous Media project "
|
|
||||||
"(https://opm-project.org).\n"
|
|
||||||
"\n"
|
|
||||||
"THE GOAL OF THE `ebos` SIMULATOR IS TO CATER FOR THE NEEDS OF "
|
|
||||||
"DEVELOPMENT AND RESEARCH. No guarantees are made for production use!";
|
|
||||||
else
|
|
||||||
return briefDescription_;
|
|
||||||
}
|
|
||||||
|
|
||||||
template<class GridView, class FluidSystem, class Scalar>
|
|
||||||
void EclGenericProblem<GridView,FluidSystem,Scalar>::
|
|
||||||
readRockParameters_(const std::vector<Scalar>& cellCenterDepths)
|
|
||||||
{
|
|
||||||
const auto& rock_config = eclState_.getSimulationConfig().rock_config();
|
|
||||||
|
|
||||||
// read the rock compressibility parameters
|
|
||||||
{
|
|
||||||
const auto& comp = rock_config.comp();
|
|
||||||
rockParams_.clear();
|
|
||||||
for (const auto& c : comp)
|
|
||||||
rockParams_.push_back( { c.pref, c.compressibility } );
|
|
||||||
}
|
|
||||||
|
|
||||||
// read the parameters for water-induced rock compaction
|
|
||||||
readRockCompactionParameters_();
|
|
||||||
|
|
||||||
unsigned numElem = gridView_.size(0);
|
|
||||||
if (eclState_.fieldProps().has_int(rock_config.rocknum_property())) {
|
|
||||||
rockTableIdx_.resize(numElem);
|
|
||||||
const auto& num = eclState_.fieldProps().get_int(rock_config.rocknum_property());
|
|
||||||
for (size_t elemIdx = 0; elemIdx < numElem; ++ elemIdx) {
|
|
||||||
rockTableIdx_[elemIdx] = num[elemIdx] - 1;
|
|
||||||
}
|
|
||||||
}
|
|
||||||
|
|
||||||
// Store overburden pressure pr element
|
|
||||||
const auto& overburdTables = eclState_.getTableManager().getOverburdTables();
|
|
||||||
if (!overburdTables.empty()) {
|
|
||||||
overburdenPressure_.resize(numElem,0.0);
|
|
||||||
size_t numRocktabTables = rock_config.num_rock_tables();
|
|
||||||
|
|
||||||
if (overburdTables.size() != numRocktabTables)
|
|
||||||
throw std::runtime_error(std::to_string(numRocktabTables) +" OVERBURD tables is expected, but " + std::to_string(overburdTables.size()) +" is provided");
|
|
||||||
|
|
||||||
std::vector<Tabulated1DFunction<Scalar>> overburdenTables(numRocktabTables);
|
|
||||||
for (size_t regionIdx = 0; regionIdx < numRocktabTables; ++regionIdx) {
|
|
||||||
const OverburdTable& overburdTable = overburdTables.template getTable<OverburdTable>(regionIdx);
|
|
||||||
overburdenTables[regionIdx].setXYContainers(overburdTable.getDepthColumn(),overburdTable.getOverburdenPressureColumn());
|
|
||||||
}
|
|
||||||
|
|
||||||
for (size_t elemIdx = 0; elemIdx < numElem; ++ elemIdx) {
|
|
||||||
unsigned tableIdx = 0;
|
|
||||||
if (!rockTableIdx_.empty()) {
|
|
||||||
tableIdx = rockTableIdx_[elemIdx];
|
|
||||||
}
|
|
||||||
overburdenPressure_[elemIdx] = overburdenTables[tableIdx].eval(cellCenterDepths[elemIdx], /*extrapolation=*/true);
|
|
||||||
}
|
|
||||||
}
|
|
||||||
}
|
|
||||||
|
|
||||||
template<class GridView, class FluidSystem, class Scalar>
|
|
||||||
void EclGenericProblem<GridView,FluidSystem,Scalar>::
|
|
||||||
readRockCompactionParameters_()
|
|
||||||
{
|
|
||||||
const auto& rock_config = eclState_.getSimulationConfig().rock_config();
|
|
||||||
|
|
||||||
if (!rock_config.active())
|
|
||||||
return; // deck does not enable rock compaction
|
|
||||||
|
|
||||||
unsigned numElem = gridView_.size(0);
|
|
||||||
switch (rock_config.hysteresis_mode()) {
|
|
||||||
case RockConfig::Hysteresis::REVERS:
|
|
||||||
break;
|
|
||||||
case RockConfig::Hysteresis::IRREVERS:
|
|
||||||
// interpolate the porv volume multiplier using the minimum pressure in the cell
|
|
||||||
// i.e. don't allow re-inflation.
|
|
||||||
minOilPressure_.resize(numElem, 1e99);
|
|
||||||
break;
|
|
||||||
default:
|
|
||||||
throw std::runtime_error("Not support ROCKOMP hysteresis option ");
|
|
||||||
}
|
|
||||||
|
|
||||||
size_t numRocktabTables = rock_config.num_rock_tables();
|
|
||||||
bool waterCompaction = rock_config.water_compaction();
|
|
||||||
|
|
||||||
if (!waterCompaction) {
|
|
||||||
const auto& rocktabTables = eclState_.getTableManager().getRocktabTables();
|
|
||||||
if (rocktabTables.size() != numRocktabTables)
|
|
||||||
throw std::runtime_error("ROCKCOMP is activated." + std::to_string(numRocktabTables)
|
|
||||||
+" ROCKTAB tables is expected, but " + std::to_string(rocktabTables.size()) +" is provided");
|
|
||||||
|
|
||||||
rockCompPoroMult_.resize(numRocktabTables);
|
|
||||||
rockCompTransMult_.resize(numRocktabTables);
|
|
||||||
for (size_t regionIdx = 0; regionIdx < numRocktabTables; ++regionIdx) {
|
|
||||||
const auto& rocktabTable = rocktabTables.template getTable<RocktabTable>(regionIdx);
|
|
||||||
const auto& pressureColumn = rocktabTable.getPressureColumn();
|
|
||||||
const auto& poroColumn = rocktabTable.getPoreVolumeMultiplierColumn();
|
|
||||||
const auto& transColumn = rocktabTable.getTransmissibilityMultiplierColumn();
|
|
||||||
rockCompPoroMult_[regionIdx].setXYContainers(pressureColumn, poroColumn);
|
|
||||||
rockCompTransMult_[regionIdx].setXYContainers(pressureColumn, transColumn);
|
|
||||||
}
|
|
||||||
} else {
|
|
||||||
const auto& rock2dTables = eclState_.getTableManager().getRock2dTables();
|
|
||||||
const auto& rock2dtrTables = eclState_.getTableManager().getRock2dtrTables();
|
|
||||||
const auto& rockwnodTables = eclState_.getTableManager().getRockwnodTables();
|
|
||||||
maxWaterSaturation_.resize(numElem, 0.0);
|
|
||||||
|
|
||||||
if (rock2dTables.size() != numRocktabTables)
|
|
||||||
throw std::runtime_error("Water compation option is selected in ROCKCOMP." + std::to_string(numRocktabTables)
|
|
||||||
+" ROCK2D tables is expected, but " + std::to_string(rock2dTables.size()) +" is provided");
|
|
||||||
|
|
||||||
if (rockwnodTables.size() != numRocktabTables)
|
|
||||||
throw std::runtime_error("Water compation option is selected in ROCKCOMP." + std::to_string(numRocktabTables)
|
|
||||||
+" ROCKWNOD tables is expected, but " + std::to_string(rockwnodTables.size()) +" is provided");
|
|
||||||
//TODO check size match
|
|
||||||
rockCompPoroMultWc_.resize(numRocktabTables, TabulatedTwoDFunction(TabulatedTwoDFunction::InterpolationPolicy::Vertical));
|
|
||||||
for (size_t regionIdx = 0; regionIdx < numRocktabTables; ++regionIdx) {
|
|
||||||
const RockwnodTable& rockwnodTable = rockwnodTables.template getTable<RockwnodTable>(regionIdx);
|
|
||||||
const auto& rock2dTable = rock2dTables[regionIdx];
|
|
||||||
|
|
||||||
if (rockwnodTable.getSaturationColumn().size() != rock2dTable.sizeMultValues())
|
|
||||||
throw std::runtime_error("Number of entries in ROCKWNOD and ROCK2D needs to match.");
|
|
||||||
|
|
||||||
for (size_t xIdx = 0; xIdx < rock2dTable.size(); ++xIdx) {
|
|
||||||
rockCompPoroMultWc_[regionIdx].appendXPos(rock2dTable.getPressureValue(xIdx));
|
|
||||||
for (size_t yIdx = 0; yIdx < rockwnodTable.getSaturationColumn().size(); ++yIdx)
|
|
||||||
rockCompPoroMultWc_[regionIdx].appendSamplePoint(xIdx,
|
|
||||||
rockwnodTable.getSaturationColumn()[yIdx],
|
|
||||||
rock2dTable.getPvmultValue(xIdx, yIdx));
|
|
||||||
}
|
|
||||||
}
|
|
||||||
|
|
||||||
if (!rock2dtrTables.empty()) {
|
|
||||||
rockCompTransMultWc_.resize(numRocktabTables, TabulatedTwoDFunction(TabulatedTwoDFunction::InterpolationPolicy::Vertical));
|
|
||||||
for (size_t regionIdx = 0; regionIdx < numRocktabTables; ++regionIdx) {
|
|
||||||
const RockwnodTable& rockwnodTable = rockwnodTables.template getTable<RockwnodTable>(regionIdx);
|
|
||||||
const auto& rock2dtrTable = rock2dtrTables[regionIdx];
|
|
||||||
|
|
||||||
if (rockwnodTable.getSaturationColumn().size() != rock2dtrTable.sizeMultValues())
|
|
||||||
throw std::runtime_error("Number of entries in ROCKWNOD and ROCK2DTR needs to match.");
|
|
||||||
|
|
||||||
for (size_t xIdx = 0; xIdx < rock2dtrTable.size(); ++xIdx) {
|
|
||||||
rockCompTransMultWc_[regionIdx].appendXPos(rock2dtrTable.getPressureValue(xIdx));
|
|
||||||
for (size_t yIdx = 0; yIdx < rockwnodTable.getSaturationColumn().size(); ++yIdx)
|
|
||||||
rockCompTransMultWc_[regionIdx].appendSamplePoint(xIdx,
|
|
||||||
rockwnodTable.getSaturationColumn()[yIdx],
|
|
||||||
rock2dtrTable.getTransMultValue(xIdx, yIdx));
|
|
||||||
}
|
|
||||||
}
|
|
||||||
}
|
|
||||||
}
|
|
||||||
}
|
|
||||||
|
|
||||||
template<class GridView, class FluidSystem, class Scalar>
|
|
||||||
Scalar EclGenericProblem<GridView,FluidSystem,Scalar>::
|
|
||||||
rockCompressibility(unsigned globalSpaceIdx) const
|
|
||||||
{
|
|
||||||
if (this->rockParams_.empty())
|
|
||||||
return 0.0;
|
|
||||||
|
|
||||||
unsigned tableIdx = 0;
|
|
||||||
if (!this->rockTableIdx_.empty()) {
|
|
||||||
tableIdx = this->rockTableIdx_[globalSpaceIdx];
|
|
||||||
}
|
|
||||||
return this->rockParams_[tableIdx].compressibility;
|
|
||||||
}
|
|
||||||
|
|
||||||
template<class GridView, class FluidSystem, class Scalar>
|
|
||||||
Scalar EclGenericProblem<GridView,FluidSystem,Scalar>::
|
|
||||||
rockReferencePressure(unsigned globalSpaceIdx) const
|
|
||||||
{
|
|
||||||
if (this->rockParams_.empty())
|
|
||||||
return 1e5;
|
|
||||||
|
|
||||||
unsigned tableIdx = 0;
|
|
||||||
if (!this->rockTableIdx_.empty()) {
|
|
||||||
tableIdx = this->rockTableIdx_[globalSpaceIdx];
|
|
||||||
}
|
|
||||||
return this->rockParams_[tableIdx].referencePressure;
|
|
||||||
}
|
|
||||||
|
|
||||||
template<class GridView, class FluidSystem, class Scalar>
|
|
||||||
Scalar EclGenericProblem<GridView,FluidSystem,Scalar>::
|
|
||||||
porosity(unsigned globalSpaceIdx, unsigned timeIdx) const
|
|
||||||
{
|
|
||||||
return this->referencePorosity_[timeIdx][globalSpaceIdx];
|
|
||||||
}
|
|
||||||
|
|
||||||
template<class GridView, class FluidSystem, class Scalar>
|
|
||||||
Scalar EclGenericProblem<GridView,FluidSystem,Scalar>::
|
|
||||||
rockFraction(unsigned elementIdx, unsigned timeIdx) const
|
|
||||||
{
|
|
||||||
const auto& fp = eclState_.fieldProps();
|
|
||||||
const std::vector<double>& poro = fp.get_double("PORO");
|
|
||||||
// the reference porosity is defined as the accumulated pore volume divided by the
|
|
||||||
// geometric volume of the element. Note that it can
|
|
||||||
// be larger than 1.0 if porevolume multipliers are used
|
|
||||||
// to for instance implement larger boundary cells
|
|
||||||
Scalar porosity = poro[elementIdx];
|
|
||||||
return referencePorosity(elementIdx, timeIdx) / porosity * (1 - porosity);
|
|
||||||
}
|
|
||||||
|
|
||||||
template<class GridView, class FluidSystem, class Scalar>
|
|
||||||
template<class T>
|
|
||||||
void EclGenericProblem<GridView,FluidSystem,Scalar>::
|
|
||||||
updateNum(const std::string& name, std::vector<T>& numbers, size_t num_regions)
|
|
||||||
{
|
|
||||||
if (!eclState_.fieldProps().has_int(name))
|
|
||||||
return;
|
|
||||||
|
|
||||||
const auto& numData = eclState_.fieldProps().get_int(name);
|
|
||||||
|
|
||||||
unsigned numElems = gridView_.size(/*codim=*/0);
|
|
||||||
numbers.resize(numElems);
|
|
||||||
for (unsigned elemIdx = 0; elemIdx < numElems; ++elemIdx) {
|
|
||||||
if (numData[elemIdx] > (int)num_regions) {
|
|
||||||
throw std::runtime_error("Values larger than maximum number of regions " + std::to_string(num_regions) + " provided in " + name);
|
|
||||||
} else if (numData[elemIdx] > 0) {
|
|
||||||
numbers[elemIdx] = static_cast<T>(numData[elemIdx]) - 1;
|
|
||||||
} else {
|
|
||||||
throw std::runtime_error("zero or negative values provided for region array: " + name);
|
|
||||||
}
|
|
||||||
}
|
|
||||||
}
|
|
||||||
|
|
||||||
template<class GridView, class FluidSystem, class Scalar>
|
|
||||||
void EclGenericProblem<GridView,FluidSystem,Scalar>::
|
|
||||||
updatePvtnum_()
|
|
||||||
{
|
|
||||||
const auto num_regions = eclState_.getTableManager().getTabdims().getNumPVTTables();
|
|
||||||
updateNum("PVTNUM", pvtnum_, num_regions);
|
|
||||||
}
|
|
||||||
|
|
||||||
template<class GridView, class FluidSystem, class Scalar>
|
|
||||||
void EclGenericProblem<GridView,FluidSystem,Scalar>::
|
|
||||||
updateSatnum_()
|
|
||||||
{
|
|
||||||
const auto num_regions = eclState_.getTableManager().getTabdims().getNumSatTables();
|
|
||||||
updateNum("SATNUM", satnum_, num_regions);
|
|
||||||
}
|
|
||||||
|
|
||||||
template<class GridView, class FluidSystem, class Scalar>
|
|
||||||
void EclGenericProblem<GridView,FluidSystem,Scalar>::
|
|
||||||
updateMiscnum_()
|
|
||||||
{
|
|
||||||
const auto num_regions = 1; // we only support single region
|
|
||||||
updateNum("MISCNUM", miscnum_, num_regions);
|
|
||||||
}
|
|
||||||
|
|
||||||
template<class GridView, class FluidSystem, class Scalar>
|
|
||||||
void EclGenericProblem<GridView,FluidSystem,Scalar>::
|
|
||||||
updatePlmixnum_()
|
|
||||||
{
|
|
||||||
const auto num_regions = 1; // we only support single region
|
|
||||||
updateNum("PLMIXNUM", plmixnum_, num_regions);
|
|
||||||
}
|
|
||||||
|
|
||||||
template<class GridView, class FluidSystem, class Scalar>
|
|
||||||
void EclGenericProblem<GridView,FluidSystem,Scalar>::
|
|
||||||
updateKrnum_()
|
|
||||||
{
|
|
||||||
const auto num_regions = eclState_.getTableManager().getTabdims().getNumSatTables();
|
|
||||||
updateNum("KRNUMX", krnumx_, num_regions);
|
|
||||||
updateNum("KRNUMY", krnumy_, num_regions);
|
|
||||||
updateNum("KRNUMZ", krnumz_, num_regions);
|
|
||||||
}
|
|
||||||
|
|
||||||
template<class GridView, class FluidSystem, class Scalar>
|
|
||||||
void EclGenericProblem<GridView,FluidSystem,Scalar>::
|
|
||||||
updateImbnum_()
|
|
||||||
{
|
|
||||||
const auto num_regions = eclState_.getTableManager().getTabdims().getNumSatTables();
|
|
||||||
updateNum("IMBNUMX", imbnumx_, num_regions);
|
|
||||||
updateNum("IMBNUMY", imbnumy_, num_regions);
|
|
||||||
updateNum("IMBNUMZ", imbnumz_, num_regions);
|
|
||||||
}
|
|
||||||
|
|
||||||
template<class GridView, class FluidSystem, class Scalar>
|
|
||||||
bool EclGenericProblem<GridView,FluidSystem,Scalar>::
|
|
||||||
vapparsActive(int episodeIdx) const
|
|
||||||
{
|
|
||||||
const auto& oilVaporizationControl = schedule_[episodeIdx].oilvap();
|
|
||||||
return (oilVaporizationControl.getType() == OilVaporizationProperties::OilVaporization::VAPPARS);
|
|
||||||
}
|
|
||||||
|
|
||||||
template<class GridView, class FluidSystem, class Scalar>
|
|
||||||
bool EclGenericProblem<GridView,FluidSystem,Scalar>::
|
|
||||||
drsdtActive_(int episodeIdx) const
|
|
||||||
{
|
|
||||||
const auto& oilVaporizationControl = schedule_[episodeIdx].oilvap();
|
|
||||||
const bool bothOilGasActive = FluidSystem::phaseIsActive(FluidSystem::oilPhaseIdx) &&
|
|
||||||
FluidSystem::phaseIsActive(FluidSystem::gasPhaseIdx);
|
|
||||||
return (oilVaporizationControl.drsdtActive() && bothOilGasActive);
|
|
||||||
}
|
|
||||||
|
|
||||||
template<class GridView, class FluidSystem, class Scalar>
|
|
||||||
bool EclGenericProblem<GridView,FluidSystem,Scalar>::
|
|
||||||
drvdtActive_(int episodeIdx) const
|
|
||||||
{
|
|
||||||
const auto& oilVaporizationControl = schedule_[episodeIdx].oilvap();
|
|
||||||
const bool bothOilGasActive = FluidSystem::phaseIsActive(FluidSystem::oilPhaseIdx) &&
|
|
||||||
FluidSystem::phaseIsActive(FluidSystem::gasPhaseIdx);
|
|
||||||
return (oilVaporizationControl.drvdtActive() && bothOilGasActive);
|
|
||||||
}
|
|
||||||
|
|
||||||
template<class GridView, class FluidSystem, class Scalar>
|
|
||||||
bool EclGenericProblem<GridView,FluidSystem,Scalar>::
|
|
||||||
drsdtConvective_(int episodeIdx) const
|
|
||||||
{
|
|
||||||
const auto& oilVaporizationControl = schedule_[episodeIdx].oilvap();
|
|
||||||
const bool bothOilGasActive = FluidSystem::phaseIsActive(FluidSystem::oilPhaseIdx) &&
|
|
||||||
FluidSystem::phaseIsActive(FluidSystem::gasPhaseIdx);
|
|
||||||
return (oilVaporizationControl.drsdtConvective() && bothOilGasActive);
|
|
||||||
}
|
|
||||||
|
|
||||||
template<class GridView, class FluidSystem, class Scalar>
|
|
||||||
bool EclGenericProblem<GridView,FluidSystem,Scalar>::
|
|
||||||
beginEpisode_(bool enableExperiments,
|
|
||||||
int episodeIdx)
|
|
||||||
{
|
|
||||||
if (enableExperiments && gridView_.comm().rank() == 0 && episodeIdx >= 0) {
|
|
||||||
// print some useful information in experimental mode. (the production
|
|
||||||
// simulator does this externally.)
|
|
||||||
std::ostringstream ss;
|
|
||||||
boost::posix_time::time_facet* facet = new boost::posix_time::time_facet("%d-%b-%Y");
|
|
||||||
boost::posix_time::ptime curDateTime =
|
|
||||||
boost::posix_time::from_time_t(schedule_.simTime(episodeIdx));
|
|
||||||
ss.imbue(std::locale(std::locale::classic(), facet));
|
|
||||||
ss << "Report step " << episodeIdx + 1
|
|
||||||
<< "/" << schedule_.size() - 1
|
|
||||||
<< " at day " << schedule_.seconds(episodeIdx)/(24*3600)
|
|
||||||
<< "/" << schedule_.seconds(schedule_.size() - 1)/(24*3600)
|
|
||||||
<< ", date = " << curDateTime.date()
|
|
||||||
<< "\n ";
|
|
||||||
OpmLog::info(ss.str());
|
|
||||||
}
|
|
||||||
|
|
||||||
const auto& events = schedule_[episodeIdx].events();
|
|
||||||
|
|
||||||
// react to TUNING changes
|
|
||||||
if (episodeIdx > 0 && enableTuning_ && events.hasEvent(ScheduleEvents::TUNING_CHANGE))
|
|
||||||
{
|
|
||||||
const auto& sched_state = schedule_[episodeIdx];
|
|
||||||
const auto& tuning = sched_state.tuning();
|
|
||||||
initialTimeStepSize_ = sched_state.max_next_tstep();
|
|
||||||
maxTimeStepAfterWellEvent_ = tuning.TMAXWC;
|
|
||||||
maxTimeStepSize_ = tuning.TSMAXZ;
|
|
||||||
restartShrinkFactor_ = 1./tuning.TSFCNV;
|
|
||||||
minTimeStepSize_ = tuning.TSMINZ;
|
|
||||||
return true;
|
|
||||||
}
|
|
||||||
|
|
||||||
return false;
|
|
||||||
}
|
|
||||||
|
|
||||||
|
|
||||||
template<class GridView, class FluidSystem, class Scalar>
|
|
||||||
void EclGenericProblem<GridView,FluidSystem,Scalar>::
|
|
||||||
beginTimeStep_(bool enableExperiments,
|
|
||||||
int episodeIdx,
|
|
||||||
int timeStepIndex,
|
|
||||||
Scalar startTime,
|
|
||||||
Scalar time,
|
|
||||||
Scalar timeStepSize,
|
|
||||||
Scalar endTime)
|
|
||||||
{
|
|
||||||
if (enableExperiments && gridView_.comm().rank() == 0 && episodeIdx >= 0) {
|
|
||||||
std::ostringstream ss;
|
|
||||||
boost::posix_time::time_facet* facet = new boost::posix_time::time_facet("%d-%b-%Y");
|
|
||||||
boost::posix_time::ptime date = boost::posix_time::from_time_t(startTime) + boost::posix_time::milliseconds(static_cast<long long>(time / prefix::milli));
|
|
||||||
ss.imbue(std::locale(std::locale::classic(), facet));
|
|
||||||
ss <<"\nTime step " << timeStepIndex << ", stepsize "
|
|
||||||
<< unit::convert::to(timeStepSize, unit::day) << " days,"
|
|
||||||
<< " at day " << (double)unit::convert::to(time, unit::day)
|
|
||||||
<< "/" << (double)unit::convert::to(endTime, unit::day)
|
|
||||||
<< ", date = " << date;
|
|
||||||
OpmLog::info(ss.str());
|
|
||||||
}
|
|
||||||
|
|
||||||
// update explicit quantities between timesteps.
|
|
||||||
const auto& oilVaporizationControl = schedule_[episodeIdx].oilvap();
|
|
||||||
if (drsdtActive_(episodeIdx))
|
|
||||||
// DRSDT is enabled
|
|
||||||
for (size_t pvtRegionIdx = 0; pvtRegionIdx < maxDRs_.size(); ++pvtRegionIdx)
|
|
||||||
maxDRs_[pvtRegionIdx] = oilVaporizationControl.getMaxDRSDT(pvtRegionIdx)*timeStepSize;
|
|
||||||
|
|
||||||
if (drvdtActive_(episodeIdx))
|
|
||||||
// DRVDT is enabled
|
|
||||||
for (size_t pvtRegionIdx = 0; pvtRegionIdx < maxDRv_.size(); ++pvtRegionIdx)
|
|
||||||
maxDRv_[pvtRegionIdx] = oilVaporizationControl.getMaxDRVDT(pvtRegionIdx)*timeStepSize;
|
|
||||||
}
|
|
||||||
|
|
||||||
template<class GridView, class FluidSystem, class Scalar>
|
|
||||||
void EclGenericProblem<GridView,FluidSystem,Scalar>::
|
|
||||||
initFluidSystem_()
|
|
||||||
{
|
|
||||||
FluidSystem::initFromState(eclState_, schedule_);
|
|
||||||
}
|
|
||||||
|
|
||||||
template<class GridView, class FluidSystem, class Scalar>
|
|
||||||
void EclGenericProblem<GridView,FluidSystem,Scalar>::
|
|
||||||
readBlackoilExtentionsInitialConditions_(size_t numDof,
|
|
||||||
bool enableSolvent,
|
|
||||||
bool enablePolymer,
|
|
||||||
bool enablePolymerMolarWeight,
|
|
||||||
bool enableMICP)
|
|
||||||
{
|
|
||||||
if (enableSolvent) {
|
|
||||||
if (eclState_.fieldProps().has_double("SSOL"))
|
|
||||||
solventSaturation_ = eclState_.fieldProps().get_double("SSOL");
|
|
||||||
else
|
|
||||||
solventSaturation_.resize(numDof, 0.0);
|
|
||||||
}
|
|
||||||
|
|
||||||
if (enablePolymer) {
|
|
||||||
if (eclState_.fieldProps().has_double("SPOLY"))
|
|
||||||
polymer_.concentration = eclState_.fieldProps().get_double("SPOLY");
|
|
||||||
else
|
|
||||||
polymer_.concentration.resize(numDof, 0.0);
|
|
||||||
}
|
|
||||||
|
|
||||||
if (enablePolymerMolarWeight) {
|
|
||||||
if (eclState_.fieldProps().has_double("SPOLYMW"))
|
|
||||||
polymer_.moleWeight = eclState_.fieldProps().get_double("SPOLYMW");
|
|
||||||
else
|
|
||||||
polymer_.moleWeight.resize(numDof, 0.0);
|
|
||||||
}
|
|
||||||
|
|
||||||
if (enableMICP) {
|
|
||||||
if (eclState_.fieldProps().has_double("SMICR"))
|
|
||||||
micp_.microbialConcentration = eclState_.fieldProps().get_double("SMICR");
|
|
||||||
else
|
|
||||||
micp_.microbialConcentration.resize(numDof, 0.0);
|
|
||||||
if (eclState_.fieldProps().has_double("SOXYG"))
|
|
||||||
micp_.oxygenConcentration = eclState_.fieldProps().get_double("SOXYG");
|
|
||||||
else
|
|
||||||
micp_.oxygenConcentration.resize(numDof, 0.0);
|
|
||||||
if (eclState_.fieldProps().has_double("SUREA"))
|
|
||||||
micp_.ureaConcentration = eclState_.fieldProps().get_double("SUREA");
|
|
||||||
else
|
|
||||||
micp_.ureaConcentration.resize(numDof, 0.0);
|
|
||||||
if (eclState_.fieldProps().has_double("SBIOF"))
|
|
||||||
micp_.biofilmConcentration = eclState_.fieldProps().get_double("SBIOF");
|
|
||||||
else
|
|
||||||
micp_.biofilmConcentration.resize(numDof, 0.0);
|
|
||||||
if (eclState_.fieldProps().has_double("SCALC"))
|
|
||||||
micp_.calciteConcentration = eclState_.fieldProps().get_double("SCALC");
|
|
||||||
else
|
|
||||||
micp_.calciteConcentration.resize(numDof, 0.0);
|
|
||||||
}
|
|
||||||
}
|
|
||||||
|
|
||||||
|
|
||||||
template<class GridView, class FluidSystem, class Scalar>
|
|
||||||
Scalar EclGenericProblem<GridView,FluidSystem,Scalar>::
|
|
||||||
maxWaterSaturation(unsigned globalDofIdx) const
|
|
||||||
{
|
|
||||||
if (maxWaterSaturation_.empty())
|
|
||||||
return 0.0;
|
|
||||||
|
|
||||||
return maxWaterSaturation_[globalDofIdx];
|
|
||||||
}
|
|
||||||
|
|
||||||
template<class GridView, class FluidSystem, class Scalar>
|
|
||||||
Scalar EclGenericProblem<GridView,FluidSystem,Scalar>::
|
|
||||||
minOilPressure(unsigned globalDofIdx) const
|
|
||||||
{
|
|
||||||
if (minOilPressure_.empty())
|
|
||||||
return 0.0;
|
|
||||||
|
|
||||||
return minOilPressure_[globalDofIdx];
|
|
||||||
}
|
|
||||||
|
|
||||||
template<class GridView, class FluidSystem, class Scalar>
|
|
||||||
Scalar EclGenericProblem<GridView,FluidSystem,Scalar>::
|
|
||||||
overburdenPressure(unsigned elementIdx) const
|
|
||||||
{
|
|
||||||
if (overburdenPressure_.empty())
|
|
||||||
return 0.0;
|
|
||||||
|
|
||||||
return overburdenPressure_[elementIdx];
|
|
||||||
}
|
|
||||||
|
|
||||||
template<class GridView, class FluidSystem, class Scalar>
|
|
||||||
Scalar EclGenericProblem<GridView,FluidSystem,Scalar>::
|
|
||||||
solventSaturation(unsigned elemIdx) const
|
|
||||||
{
|
|
||||||
if (solventSaturation_.empty())
|
|
||||||
return 0;
|
|
||||||
|
|
||||||
return solventSaturation_[elemIdx];
|
|
||||||
}
|
|
||||||
|
|
||||||
template<class GridView, class FluidSystem, class Scalar>
|
|
||||||
Scalar EclGenericProblem<GridView,FluidSystem,Scalar>::
|
|
||||||
drsdtcon(unsigned elemIdx, int episodeIdx) const
|
|
||||||
{
|
|
||||||
if (convectiveDrs_.empty())
|
|
||||||
return 0;
|
|
||||||
|
|
||||||
// The episode index is set to -1 in the initialization phase.
|
|
||||||
// Output drsdt value for index 0
|
|
||||||
episodeIdx = std::max(episodeIdx, 0);
|
|
||||||
const auto& oilVaporizationControl = schedule_[episodeIdx].oilvap();
|
|
||||||
int pvtRegionIdx = pvtRegionIndex(elemIdx);
|
|
||||||
return oilVaporizationControl.getMaxDRSDT(pvtRegionIdx)*convectiveDrs_[elemIdx];
|
|
||||||
}
|
|
||||||
|
|
||||||
template<class GridView, class FluidSystem, class Scalar>
|
|
||||||
Scalar EclGenericProblem<GridView,FluidSystem,Scalar>::
|
|
||||||
polymerConcentration(unsigned elemIdx) const
|
|
||||||
{
|
|
||||||
if (polymer_.concentration.empty()) {
|
|
||||||
return 0;
|
|
||||||
}
|
|
||||||
|
|
||||||
return polymer_.concentration[elemIdx];
|
|
||||||
}
|
|
||||||
|
|
||||||
template<class GridView, class FluidSystem, class Scalar>
|
|
||||||
Scalar EclGenericProblem<GridView,FluidSystem,Scalar>::
|
|
||||||
polymerMolecularWeight(const unsigned elemIdx) const
|
|
||||||
{
|
|
||||||
if (polymer_.moleWeight.empty()) {
|
|
||||||
return 0.0;
|
|
||||||
}
|
|
||||||
|
|
||||||
return polymer_.moleWeight[elemIdx];
|
|
||||||
}
|
|
||||||
|
|
||||||
template<class GridView, class FluidSystem, class Scalar>
|
|
||||||
Scalar EclGenericProblem<GridView,FluidSystem,Scalar>::
|
|
||||||
microbialConcentration(unsigned elemIdx) const
|
|
||||||
{
|
|
||||||
if (micp_.microbialConcentration.empty()) {
|
|
||||||
return 0;
|
|
||||||
}
|
|
||||||
|
|
||||||
return micp_.microbialConcentration[elemIdx];
|
|
||||||
}
|
|
||||||
|
|
||||||
template<class GridView, class FluidSystem, class Scalar>
|
|
||||||
Scalar EclGenericProblem<GridView,FluidSystem,Scalar>::
|
|
||||||
oxygenConcentration(unsigned elemIdx) const
|
|
||||||
{
|
|
||||||
if (micp_.oxygenConcentration.empty()) {
|
|
||||||
return 0;
|
|
||||||
}
|
|
||||||
|
|
||||||
return micp_.oxygenConcentration[elemIdx];
|
|
||||||
}
|
|
||||||
|
|
||||||
template<class GridView, class FluidSystem, class Scalar>
|
|
||||||
Scalar EclGenericProblem<GridView,FluidSystem,Scalar>::
|
|
||||||
ureaConcentration(unsigned elemIdx) const
|
|
||||||
{
|
|
||||||
if (micp_.ureaConcentration.empty()) {
|
|
||||||
return 0;
|
|
||||||
}
|
|
||||||
|
|
||||||
return micp_.ureaConcentration[elemIdx];
|
|
||||||
}
|
|
||||||
|
|
||||||
template<class GridView, class FluidSystem, class Scalar>
|
|
||||||
Scalar EclGenericProblem<GridView,FluidSystem,Scalar>::
|
|
||||||
biofilmConcentration(unsigned elemIdx) const
|
|
||||||
{
|
|
||||||
if (micp_.biofilmConcentration.empty()) {
|
|
||||||
return 0;
|
|
||||||
}
|
|
||||||
|
|
||||||
return micp_.biofilmConcentration[elemIdx];
|
|
||||||
}
|
|
||||||
|
|
||||||
template<class GridView, class FluidSystem, class Scalar>
|
|
||||||
Scalar EclGenericProblem<GridView,FluidSystem,Scalar>::
|
|
||||||
calciteConcentration(unsigned elemIdx) const
|
|
||||||
{
|
|
||||||
if (micp_.calciteConcentration.empty()) {
|
|
||||||
return 0;
|
|
||||||
}
|
|
||||||
|
|
||||||
return micp_.calciteConcentration[elemIdx];
|
|
||||||
}
|
|
||||||
|
|
||||||
template<class GridView, class FluidSystem, class Scalar>
|
|
||||||
unsigned EclGenericProblem<GridView,FluidSystem,Scalar>::
|
|
||||||
pvtRegionIndex(unsigned elemIdx) const
|
|
||||||
{
|
|
||||||
if (pvtnum_.empty())
|
|
||||||
return 0;
|
|
||||||
|
|
||||||
return pvtnum_[elemIdx];
|
|
||||||
}
|
|
||||||
|
|
||||||
template<class GridView, class FluidSystem, class Scalar>
|
|
||||||
unsigned EclGenericProblem<GridView,FluidSystem,Scalar>::
|
|
||||||
satnumRegionIndex(unsigned elemIdx) const
|
|
||||||
{
|
|
||||||
if (satnum_.empty())
|
|
||||||
return 0;
|
|
||||||
|
|
||||||
return satnum_[elemIdx];
|
|
||||||
}
|
|
||||||
|
|
||||||
template<class GridView, class FluidSystem, class Scalar>
|
|
||||||
unsigned EclGenericProblem<GridView,FluidSystem,Scalar>::
|
|
||||||
miscnumRegionIndex(unsigned elemIdx) const
|
|
||||||
{
|
|
||||||
if (miscnum_.empty())
|
|
||||||
return 0;
|
|
||||||
|
|
||||||
return miscnum_[elemIdx];
|
|
||||||
}
|
|
||||||
|
|
||||||
template<class GridView, class FluidSystem, class Scalar>
|
|
||||||
unsigned EclGenericProblem<GridView,FluidSystem,Scalar>::
|
|
||||||
plmixnumRegionIndex(unsigned elemIdx) const
|
|
||||||
{
|
|
||||||
if (plmixnum_.empty())
|
|
||||||
return 0;
|
|
||||||
|
|
||||||
return plmixnum_[elemIdx];
|
|
||||||
}
|
|
||||||
|
|
||||||
template<class GridView, class FluidSystem, class Scalar>
|
|
||||||
Scalar EclGenericProblem<GridView,FluidSystem,Scalar>::
|
|
||||||
maxPolymerAdsorption(unsigned elemIdx) const
|
|
||||||
{
|
|
||||||
if (polymer_.maxAdsorption.empty()) {
|
|
||||||
return 0;
|
|
||||||
}
|
|
||||||
|
|
||||||
return polymer_.maxAdsorption[elemIdx];
|
|
||||||
}
|
|
||||||
|
|
||||||
template<class GridView, class FluidSystem, class Scalar>
|
|
||||||
void EclGenericProblem<GridView,FluidSystem,Scalar>::
|
|
||||||
initDRSDT_(size_t numDof,
|
|
||||||
int episodeIdx)
|
|
||||||
{
|
|
||||||
// deal with DRSDT
|
|
||||||
unsigned ntpvt = eclState_.runspec().tabdims().getNumPVTTables();
|
|
||||||
//TODO We may want to only allocate these properties only if active.
|
|
||||||
//But since they may be activated at later time we need some more
|
|
||||||
//intrastructure to handle it
|
|
||||||
if (FluidSystem::phaseIsActive(FluidSystem::oilPhaseIdx) && FluidSystem::phaseIsActive(FluidSystem::gasPhaseIdx)) {
|
|
||||||
maxDRv_.resize(ntpvt, 1e30);
|
|
||||||
lastRv_.resize(numDof, 0.0);
|
|
||||||
maxDRs_.resize(ntpvt, 1e30);
|
|
||||||
dRsDtOnlyFreeGas_.resize(ntpvt, false);
|
|
||||||
lastRs_.resize(numDof, 0.0);
|
|
||||||
maxDRv_.resize(ntpvt, 1e30);
|
|
||||||
lastRv_.resize(numDof, 0.0);
|
|
||||||
maxOilSaturation_.resize(numDof, 0.0);
|
|
||||||
if (drsdtConvective_(episodeIdx)) {
|
|
||||||
convectiveDrs_.resize(numDof, 1.0);
|
|
||||||
}
|
|
||||||
}
|
|
||||||
}
|
|
||||||
|
|
||||||
template<class GridView, class FluidSystem, class Scalar>
|
|
||||||
bool EclGenericProblem<GridView,FluidSystem,Scalar>::
|
|
||||||
operator==(const EclGenericProblem& rhs) const
|
|
||||||
{
|
|
||||||
return this->maxOilSaturation_ == rhs.maxOilSaturation_ &&
|
|
||||||
this->maxWaterSaturation_ == rhs.maxWaterSaturation_ &&
|
|
||||||
this->minOilPressure_ == rhs.minOilPressure_ &&
|
|
||||||
this->overburdenPressure_ == rhs.overburdenPressure_ &&
|
|
||||||
this->polymer_ == rhs.polymer_ &&
|
|
||||||
this->solventSaturation_ == rhs.solventSaturation_ &&
|
|
||||||
this->micp_ == rhs.micp_ &&
|
|
||||||
this->lastRv_ == rhs.lastRv_ &&
|
|
||||||
this->maxDRv_ == rhs.maxDRv_ &&
|
|
||||||
this->convectiveDrs_ == rhs.convectiveDrs_ &&
|
|
||||||
this->lastRs_ == rhs.lastRs_ &&
|
|
||||||
this->maxDRs_ == rhs.maxDRs_ &&
|
|
||||||
this->dRsDtOnlyFreeGas_ == rhs.dRsDtOnlyFreeGas_;
|
|
||||||
}
|
|
||||||
|
|
||||||
#if HAVE_DUNE_FEM
|
|
||||||
template class EclGenericProblem<Dune::GridView<Dune::Fem::GridPart2GridViewTraits<Dune::Fem::AdaptiveLeafGridPart<Dune::CpGrid, Dune::PartitionIteratorType(4), false>>>,
|
|
||||||
BlackOilFluidSystem<double,BlackOilDefaultIndexTraits>,
|
|
||||||
double>;
|
|
||||||
template class EclGenericProblem<Dune::Fem::GridPart2GridViewImpl<
|
|
||||||
Dune::Fem::AdaptiveLeafGridPart<
|
|
||||||
Dune::CpGrid,
|
|
||||||
Dune::PartitionIteratorType(4),
|
|
||||||
false>>,
|
|
||||||
Opm::BlackOilFluidSystem<
|
|
||||||
double,
|
|
||||||
Opm::BlackOilDefaultIndexTraits>,
|
|
||||||
double>;
|
|
||||||
#if HAVE_DUNE_ALUGRID
|
|
||||||
|
|
||||||
#if HAVE_MPI
|
|
||||||
using ALUGrid3CN = Dune::ALUGrid<3, 3, Dune::cube, Dune::nonconforming, Dune::ALUGridMPIComm>;
|
|
||||||
#else
|
|
||||||
using ALUGrid3CN = Dune::ALUGrid<3, 3, Dune::cube, Dune::nonconforming, Dune::ALUGridNoComm>;
|
|
||||||
#endif //HAVE_MPI
|
|
||||||
|
|
||||||
template class EclGenericProblem<Dune::GridView<Dune::Fem::GridPart2GridViewTraits<Dune::Fem::AdaptiveLeafGridPart<ALUGrid3CN, Dune::PartitionIteratorType(4), false>>>,
|
|
||||||
BlackOilFluidSystem<double,BlackOilDefaultIndexTraits>,
|
|
||||||
double>;
|
|
||||||
template class EclGenericProblem<Dune::Fem::GridPart2GridViewImpl<
|
|
||||||
Dune::Fem::AdaptiveLeafGridPart<
|
|
||||||
ALUGrid3CN,
|
|
||||||
Dune::PartitionIteratorType(4),
|
|
||||||
false>>,
|
|
||||||
Opm::BlackOilFluidSystem<
|
|
||||||
double,
|
|
||||||
Opm::BlackOilDefaultIndexTraits>,
|
|
||||||
double>;
|
|
||||||
|
|
||||||
#endif //HAVE_DUNE_ALUGRID
|
|
||||||
#else
|
|
||||||
template class EclGenericProblem<Dune::GridView<Dune::DefaultLeafGridViewTraits<Dune::CpGrid>>,
|
|
||||||
BlackOilFluidSystem<double,BlackOilDefaultIndexTraits>,
|
|
||||||
double>;
|
|
||||||
#if HAVE_DUNE_ALUGRID
|
|
||||||
#if HAVE_MPI
|
|
||||||
using ALUGrid3CN = const Dune::ALUGrid<3, 3, Dune::cube, Dune::nonconforming, Dune::ALUGridMPIComm>;
|
|
||||||
#else
|
|
||||||
using ALUGrid3CN = const Dune::ALUGrid<3, 3, Dune::cube, Dune::nonconforming, Dune::ALUGridNoComm>;
|
|
||||||
#endif //HAVE_MPI
|
|
||||||
|
|
||||||
template class EclGenericProblem<Dune::GridView<Dune::ALU3dLeafGridViewTraits<ALUGrid3CN, Dune::PartitionIteratorType(4)>>,
|
|
||||||
BlackOilFluidSystem<double,BlackOilDefaultIndexTraits>,
|
|
||||||
double>;
|
|
||||||
|
|
||||||
#endif //HAVE_DUNE_ALUGRID
|
|
||||||
#endif //HAVE_DUNE_FEM
|
|
||||||
|
|
||||||
template class EclGenericProblem<Dune::GridView<Dune::PolyhedralGridViewTraits<3,3,double,Dune::PartitionIteratorType(4)>>,
|
|
||||||
BlackOilFluidSystem<double,BlackOilDefaultIndexTraits>,
|
|
||||||
double>;
|
|
||||||
|
|
||||||
} // namespace Opm
|
|
||||||
|
|||||||
826
ebos/eclgenericproblem_impl.hh
Normal file
826
ebos/eclgenericproblem_impl.hh
Normal file
@ -0,0 +1,826 @@
|
|||||||
|
// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
|
||||||
|
// vi: set et ts=4 sw=4 sts=4:
|
||||||
|
/*
|
||||||
|
This file is part of the Open Porous Media project (OPM).
|
||||||
|
|
||||||
|
OPM is free software: you can redistribute it and/or modify
|
||||||
|
it under the terms of the GNU General Public License as published by
|
||||||
|
the Free Software Foundation, either version 2 of the License, or
|
||||||
|
(at your option) any later version.
|
||||||
|
|
||||||
|
OPM is distributed in the hope that it will be useful,
|
||||||
|
but WITHOUT ANY WARRANTY; without even the implied warranty of
|
||||||
|
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
||||||
|
GNU General Public License for more details.
|
||||||
|
|
||||||
|
You should have received a copy of the GNU General Public License
|
||||||
|
along with OPM. If not, see <http://www.gnu.org/licenses/>.
|
||||||
|
|
||||||
|
Consult the COPYING file in the top-level source directory of this
|
||||||
|
module for the precise wording of the license and the list of
|
||||||
|
copyright holders.
|
||||||
|
*/
|
||||||
|
#ifndef EWOMS_GENERIC_ECL_PROBLEM_IMPL_HH
|
||||||
|
#define EWOMS_GENERIC_ECL_PROBLEM_IPML_HH
|
||||||
|
|
||||||
|
#include <ebos/eclgenericproblem.hh>
|
||||||
|
|
||||||
|
#include <dune/common/parametertree.hh>
|
||||||
|
|
||||||
|
#include <opm/input/eclipse/EclipseState/EclipseState.hpp>
|
||||||
|
#include <opm/input/eclipse/EclipseState/Tables/OverburdTable.hpp>
|
||||||
|
#include <opm/input/eclipse/EclipseState/Tables/RockwnodTable.hpp>
|
||||||
|
#include <opm/input/eclipse/Schedule/Schedule.hpp>
|
||||||
|
#include <opm/input/eclipse/Units/Units.hpp>
|
||||||
|
|
||||||
|
#include <opm/material/fluidsystems/BlackOilFluidSystem.hpp>
|
||||||
|
|
||||||
|
#include <boost/date_time.hpp>
|
||||||
|
|
||||||
|
#include <algorithm>
|
||||||
|
#include <limits>
|
||||||
|
#include <stdexcept>
|
||||||
|
#include <iostream>
|
||||||
|
|
||||||
|
namespace Opm {
|
||||||
|
|
||||||
|
int eclPositionalParameter(Dune::ParameterTree& tree,
|
||||||
|
std::set<std::string>& seenParams,
|
||||||
|
std::string& errorMsg,
|
||||||
|
const char** argv,
|
||||||
|
int paramIdx)
|
||||||
|
{
|
||||||
|
std::string param = argv[paramIdx];
|
||||||
|
size_t i = param.find('=');
|
||||||
|
if (i != std::string::npos) {
|
||||||
|
std::string oldParamName = param.substr(0, i);
|
||||||
|
std::string oldParamValue = param.substr(i+1);
|
||||||
|
std::string newParamName = "--" + oldParamName;
|
||||||
|
std::replace(newParamName.begin(),
|
||||||
|
newParamName.end(), '_' , '-');
|
||||||
|
errorMsg =
|
||||||
|
"The old syntax to specify parameters on the command line is no longer supported: "
|
||||||
|
"Try replacing '" + oldParamName + "=" + oldParamValue + "' with "+
|
||||||
|
"'" + newParamName + "=" + oldParamValue + "'!";
|
||||||
|
return 0;
|
||||||
|
}
|
||||||
|
|
||||||
|
if (seenParams.count("EclDeckFileName") > 0) {
|
||||||
|
errorMsg =
|
||||||
|
"Parameter 'EclDeckFileName' specified multiple times"
|
||||||
|
" as a command line parameter";
|
||||||
|
return 0;
|
||||||
|
}
|
||||||
|
|
||||||
|
tree["EclDeckFileName"] = argv[paramIdx];
|
||||||
|
seenParams.insert("EclDeckFileName");
|
||||||
|
return 1;
|
||||||
|
}
|
||||||
|
|
||||||
|
template<class GridView, class FluidSystem, class Scalar>
|
||||||
|
EclGenericProblem<GridView,FluidSystem,Scalar>::
|
||||||
|
EclGenericProblem(const EclipseState& eclState,
|
||||||
|
const Schedule& schedule,
|
||||||
|
const GridView& gridView)
|
||||||
|
: eclState_(eclState)
|
||||||
|
, schedule_(schedule)
|
||||||
|
, gridView_(gridView)
|
||||||
|
{
|
||||||
|
}
|
||||||
|
|
||||||
|
template<class GridView, class FluidSystem, class Scalar>
|
||||||
|
EclGenericProblem<GridView,FluidSystem,Scalar>
|
||||||
|
EclGenericProblem<GridView,FluidSystem,Scalar>::
|
||||||
|
serializationTestObject(const EclipseState& eclState,
|
||||||
|
const Schedule& schedule,
|
||||||
|
const GridView& gridView)
|
||||||
|
{
|
||||||
|
EclGenericProblem result(eclState, schedule, gridView);
|
||||||
|
result.maxOilSaturation_ = {1.0, 2.0};
|
||||||
|
result.maxPolymerAdsorption_ = {3.0, 4.0, 5.0};
|
||||||
|
result.maxWaterSaturation_ = {6.0};
|
||||||
|
result.minOilPressure_ = {7.0, 8.0, 9.0, 10.0};
|
||||||
|
result.overburdenPressure_ = {11.0};
|
||||||
|
result.polymerConcentration_ = {12.0};
|
||||||
|
result.polymerMoleWeight_ = {13.0, 14.0};
|
||||||
|
result.solventSaturation_ = {15.0};
|
||||||
|
result.microbialConcentration_ = {16.0};
|
||||||
|
result.oxygenConcentration_ = {17.0};
|
||||||
|
result.ureaConcentration_ = {18.0};
|
||||||
|
result.biofilmConcentration_ = {19.0};
|
||||||
|
result.calciteConcentration_ = {20.0};
|
||||||
|
result.lastRv_ = {21.0};
|
||||||
|
result.maxDRv_ = {22.0, 23.0};
|
||||||
|
result.convectiveDrs_ = {24.0, 25.0, 26.0};
|
||||||
|
result.lastRs_ = {27.0};
|
||||||
|
result.maxDRs_ = {28.0};
|
||||||
|
result.dRsDtOnlyFreeGas_ = {false, true};
|
||||||
|
|
||||||
|
return result;
|
||||||
|
}
|
||||||
|
|
||||||
|
template<class GridView, class FluidSystem, class Scalar>
|
||||||
|
std::string
|
||||||
|
EclGenericProblem<GridView,FluidSystem,Scalar>::
|
||||||
|
helpPreamble(int,
|
||||||
|
const char **argv)
|
||||||
|
{
|
||||||
|
std::string desc = EclGenericProblem::briefDescription();
|
||||||
|
if (!desc.empty())
|
||||||
|
desc = desc + "\n";
|
||||||
|
|
||||||
|
return
|
||||||
|
"Usage: "+std::string(argv[0]) + " [OPTIONS] [ECL_DECK_FILENAME]\n"
|
||||||
|
+ desc;
|
||||||
|
}
|
||||||
|
|
||||||
|
template<class GridView, class FluidSystem, class Scalar>
|
||||||
|
std::string
|
||||||
|
EclGenericProblem<GridView,FluidSystem,Scalar>::
|
||||||
|
briefDescription()
|
||||||
|
{
|
||||||
|
if (briefDescription_.empty())
|
||||||
|
return
|
||||||
|
"The Ecl-deck Black-Oil reservoir Simulator (ebos); a hydrocarbon "
|
||||||
|
"reservoir simulation program that processes ECL-formatted input "
|
||||||
|
"files that is part of the Open Porous Media project "
|
||||||
|
"(https://opm-project.org).\n"
|
||||||
|
"\n"
|
||||||
|
"THE GOAL OF THE `ebos` SIMULATOR IS TO CATER FOR THE NEEDS OF "
|
||||||
|
"DEVELOPMENT AND RESEARCH. No guarantees are made for production use!";
|
||||||
|
else
|
||||||
|
return briefDescription_;
|
||||||
|
}
|
||||||
|
|
||||||
|
template<class GridView, class FluidSystem, class Scalar>
|
||||||
|
void EclGenericProblem<GridView,FluidSystem,Scalar>::
|
||||||
|
readRockParameters_(const std::vector<Scalar>& cellCenterDepths)
|
||||||
|
{
|
||||||
|
const auto& rock_config = eclState_.getSimulationConfig().rock_config();
|
||||||
|
|
||||||
|
// read the rock compressibility parameters
|
||||||
|
{
|
||||||
|
const auto& comp = rock_config.comp();
|
||||||
|
rockParams_.clear();
|
||||||
|
for (const auto& c : comp)
|
||||||
|
rockParams_.push_back( { c.pref, c.compressibility } );
|
||||||
|
}
|
||||||
|
|
||||||
|
// read the parameters for water-induced rock compaction
|
||||||
|
readRockCompactionParameters_();
|
||||||
|
|
||||||
|
unsigned numElem = gridView_.size(0);
|
||||||
|
if (eclState_.fieldProps().has_int(rock_config.rocknum_property())) {
|
||||||
|
rockTableIdx_.resize(numElem);
|
||||||
|
const auto& num = eclState_.fieldProps().get_int(rock_config.rocknum_property());
|
||||||
|
for (size_t elemIdx = 0; elemIdx < numElem; ++ elemIdx) {
|
||||||
|
rockTableIdx_[elemIdx] = num[elemIdx] - 1;
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
// Store overburden pressure pr element
|
||||||
|
const auto& overburdTables = eclState_.getTableManager().getOverburdTables();
|
||||||
|
if (!overburdTables.empty()) {
|
||||||
|
overburdenPressure_.resize(numElem,0.0);
|
||||||
|
size_t numRocktabTables = rock_config.num_rock_tables();
|
||||||
|
|
||||||
|
if (overburdTables.size() != numRocktabTables)
|
||||||
|
throw std::runtime_error(std::to_string(numRocktabTables) +" OVERBURD tables is expected, but " + std::to_string(overburdTables.size()) +" is provided");
|
||||||
|
|
||||||
|
std::vector<Tabulated1DFunction<Scalar>> overburdenTables(numRocktabTables);
|
||||||
|
for (size_t regionIdx = 0; regionIdx < numRocktabTables; ++regionIdx) {
|
||||||
|
const OverburdTable& overburdTable = overburdTables.template getTable<OverburdTable>(regionIdx);
|
||||||
|
overburdenTables[regionIdx].setXYContainers(overburdTable.getDepthColumn(),overburdTable.getOverburdenPressureColumn());
|
||||||
|
}
|
||||||
|
|
||||||
|
for (size_t elemIdx = 0; elemIdx < numElem; ++ elemIdx) {
|
||||||
|
unsigned tableIdx = 0;
|
||||||
|
if (!rockTableIdx_.empty()) {
|
||||||
|
tableIdx = rockTableIdx_[elemIdx];
|
||||||
|
}
|
||||||
|
overburdenPressure_[elemIdx] = overburdenTables[tableIdx].eval(cellCenterDepths[elemIdx], /*extrapolation=*/true);
|
||||||
|
}
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
template<class GridView, class FluidSystem, class Scalar>
|
||||||
|
void EclGenericProblem<GridView,FluidSystem,Scalar>::
|
||||||
|
readRockCompactionParameters_()
|
||||||
|
{
|
||||||
|
const auto& rock_config = eclState_.getSimulationConfig().rock_config();
|
||||||
|
|
||||||
|
if (!rock_config.active())
|
||||||
|
return; // deck does not enable rock compaction
|
||||||
|
|
||||||
|
unsigned numElem = gridView_.size(0);
|
||||||
|
switch (rock_config.hysteresis_mode()) {
|
||||||
|
case RockConfig::Hysteresis::REVERS:
|
||||||
|
break;
|
||||||
|
case RockConfig::Hysteresis::IRREVERS:
|
||||||
|
// interpolate the porv volume multiplier using the minimum pressure in the cell
|
||||||
|
// i.e. don't allow re-inflation.
|
||||||
|
minOilPressure_.resize(numElem, 1e99);
|
||||||
|
break;
|
||||||
|
default:
|
||||||
|
throw std::runtime_error("Not support ROCKOMP hysteresis option ");
|
||||||
|
}
|
||||||
|
|
||||||
|
size_t numRocktabTables = rock_config.num_rock_tables();
|
||||||
|
bool waterCompaction = rock_config.water_compaction();
|
||||||
|
|
||||||
|
if (!waterCompaction) {
|
||||||
|
const auto& rocktabTables = eclState_.getTableManager().getRocktabTables();
|
||||||
|
if (rocktabTables.size() != numRocktabTables)
|
||||||
|
throw std::runtime_error("ROCKCOMP is activated." + std::to_string(numRocktabTables)
|
||||||
|
+" ROCKTAB tables is expected, but " + std::to_string(rocktabTables.size()) +" is provided");
|
||||||
|
|
||||||
|
rockCompPoroMult_.resize(numRocktabTables);
|
||||||
|
rockCompTransMult_.resize(numRocktabTables);
|
||||||
|
for (size_t regionIdx = 0; regionIdx < numRocktabTables; ++regionIdx) {
|
||||||
|
const auto& rocktabTable = rocktabTables.template getTable<RocktabTable>(regionIdx);
|
||||||
|
const auto& pressureColumn = rocktabTable.getPressureColumn();
|
||||||
|
const auto& poroColumn = rocktabTable.getPoreVolumeMultiplierColumn();
|
||||||
|
const auto& transColumn = rocktabTable.getTransmissibilityMultiplierColumn();
|
||||||
|
rockCompPoroMult_[regionIdx].setXYContainers(pressureColumn, poroColumn);
|
||||||
|
rockCompTransMult_[regionIdx].setXYContainers(pressureColumn, transColumn);
|
||||||
|
}
|
||||||
|
} else {
|
||||||
|
const auto& rock2dTables = eclState_.getTableManager().getRock2dTables();
|
||||||
|
const auto& rock2dtrTables = eclState_.getTableManager().getRock2dtrTables();
|
||||||
|
const auto& rockwnodTables = eclState_.getTableManager().getRockwnodTables();
|
||||||
|
maxWaterSaturation_.resize(numElem, 0.0);
|
||||||
|
|
||||||
|
if (rock2dTables.size() != numRocktabTables)
|
||||||
|
throw std::runtime_error("Water compation option is selected in ROCKCOMP." + std::to_string(numRocktabTables)
|
||||||
|
+" ROCK2D tables is expected, but " + std::to_string(rock2dTables.size()) +" is provided");
|
||||||
|
|
||||||
|
if (rockwnodTables.size() != numRocktabTables)
|
||||||
|
throw std::runtime_error("Water compation option is selected in ROCKCOMP." + std::to_string(numRocktabTables)
|
||||||
|
+" ROCKWNOD tables is expected, but " + std::to_string(rockwnodTables.size()) +" is provided");
|
||||||
|
//TODO check size match
|
||||||
|
rockCompPoroMultWc_.resize(numRocktabTables, TabulatedTwoDFunction(TabulatedTwoDFunction::InterpolationPolicy::Vertical));
|
||||||
|
for (size_t regionIdx = 0; regionIdx < numRocktabTables; ++regionIdx) {
|
||||||
|
const RockwnodTable& rockwnodTable = rockwnodTables.template getTable<RockwnodTable>(regionIdx);
|
||||||
|
const auto& rock2dTable = rock2dTables[regionIdx];
|
||||||
|
|
||||||
|
if (rockwnodTable.getSaturationColumn().size() != rock2dTable.sizeMultValues())
|
||||||
|
throw std::runtime_error("Number of entries in ROCKWNOD and ROCK2D needs to match.");
|
||||||
|
|
||||||
|
for (size_t xIdx = 0; xIdx < rock2dTable.size(); ++xIdx) {
|
||||||
|
rockCompPoroMultWc_[regionIdx].appendXPos(rock2dTable.getPressureValue(xIdx));
|
||||||
|
for (size_t yIdx = 0; yIdx < rockwnodTable.getSaturationColumn().size(); ++yIdx)
|
||||||
|
rockCompPoroMultWc_[regionIdx].appendSamplePoint(xIdx,
|
||||||
|
rockwnodTable.getSaturationColumn()[yIdx],
|
||||||
|
rock2dTable.getPvmultValue(xIdx, yIdx));
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
if (!rock2dtrTables.empty()) {
|
||||||
|
rockCompTransMultWc_.resize(numRocktabTables, TabulatedTwoDFunction(TabulatedTwoDFunction::InterpolationPolicy::Vertical));
|
||||||
|
for (size_t regionIdx = 0; regionIdx < numRocktabTables; ++regionIdx) {
|
||||||
|
const RockwnodTable& rockwnodTable = rockwnodTables.template getTable<RockwnodTable>(regionIdx);
|
||||||
|
const auto& rock2dtrTable = rock2dtrTables[regionIdx];
|
||||||
|
|
||||||
|
if (rockwnodTable.getSaturationColumn().size() != rock2dtrTable.sizeMultValues())
|
||||||
|
throw std::runtime_error("Number of entries in ROCKWNOD and ROCK2DTR needs to match.");
|
||||||
|
|
||||||
|
for (size_t xIdx = 0; xIdx < rock2dtrTable.size(); ++xIdx) {
|
||||||
|
rockCompTransMultWc_[regionIdx].appendXPos(rock2dtrTable.getPressureValue(xIdx));
|
||||||
|
for (size_t yIdx = 0; yIdx < rockwnodTable.getSaturationColumn().size(); ++yIdx)
|
||||||
|
rockCompTransMultWc_[regionIdx].appendSamplePoint(xIdx,
|
||||||
|
rockwnodTable.getSaturationColumn()[yIdx],
|
||||||
|
rock2dtrTable.getTransMultValue(xIdx, yIdx));
|
||||||
|
}
|
||||||
|
}
|
||||||
|
}
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
template<class GridView, class FluidSystem, class Scalar>
|
||||||
|
Scalar EclGenericProblem<GridView,FluidSystem,Scalar>::
|
||||||
|
rockCompressibility(unsigned globalSpaceIdx) const
|
||||||
|
{
|
||||||
|
if (this->rockParams_.empty())
|
||||||
|
return 0.0;
|
||||||
|
|
||||||
|
unsigned tableIdx = 0;
|
||||||
|
if (!this->rockTableIdx_.empty()) {
|
||||||
|
tableIdx = this->rockTableIdx_[globalSpaceIdx];
|
||||||
|
}
|
||||||
|
return this->rockParams_[tableIdx].compressibility;
|
||||||
|
}
|
||||||
|
|
||||||
|
template<class GridView, class FluidSystem, class Scalar>
|
||||||
|
Scalar EclGenericProblem<GridView,FluidSystem,Scalar>::
|
||||||
|
rockReferencePressure(unsigned globalSpaceIdx) const
|
||||||
|
{
|
||||||
|
if (this->rockParams_.empty())
|
||||||
|
return 1e5;
|
||||||
|
|
||||||
|
unsigned tableIdx = 0;
|
||||||
|
if (!this->rockTableIdx_.empty()) {
|
||||||
|
tableIdx = this->rockTableIdx_[globalSpaceIdx];
|
||||||
|
}
|
||||||
|
return this->rockParams_[tableIdx].referencePressure;
|
||||||
|
}
|
||||||
|
|
||||||
|
template<class GridView, class FluidSystem, class Scalar>
|
||||||
|
Scalar EclGenericProblem<GridView,FluidSystem,Scalar>::
|
||||||
|
porosity(unsigned globalSpaceIdx, unsigned timeIdx) const
|
||||||
|
{
|
||||||
|
return this->referencePorosity_[timeIdx][globalSpaceIdx];
|
||||||
|
}
|
||||||
|
|
||||||
|
template<class GridView, class FluidSystem, class Scalar>
|
||||||
|
Scalar EclGenericProblem<GridView,FluidSystem,Scalar>::
|
||||||
|
rockFraction(unsigned elementIdx, unsigned timeIdx) const
|
||||||
|
{
|
||||||
|
const auto& fp = eclState_.fieldProps();
|
||||||
|
const std::vector<double>& poro = fp.get_double("PORO");
|
||||||
|
// the reference porosity is defined as the accumulated pore volume divided by the
|
||||||
|
// geometric volume of the element. Note that it can
|
||||||
|
// be larger than 1.0 if porevolume multipliers are used
|
||||||
|
// to for instance implement larger boundary cells
|
||||||
|
Scalar porosity = poro[elementIdx];
|
||||||
|
return referencePorosity(elementIdx, timeIdx) / porosity * (1 - porosity);
|
||||||
|
}
|
||||||
|
|
||||||
|
template<class GridView, class FluidSystem, class Scalar>
|
||||||
|
template<class T>
|
||||||
|
void EclGenericProblem<GridView,FluidSystem,Scalar>::
|
||||||
|
updateNum(const std::string& name, std::vector<T>& numbers, size_t num_regions)
|
||||||
|
{
|
||||||
|
if (!eclState_.fieldProps().has_int(name))
|
||||||
|
return;
|
||||||
|
|
||||||
|
const auto& numData = eclState_.fieldProps().get_int(name);
|
||||||
|
|
||||||
|
unsigned numElems = gridView_.size(/*codim=*/0);
|
||||||
|
numbers.resize(numElems);
|
||||||
|
for (unsigned elemIdx = 0; elemIdx < numElems; ++elemIdx) {
|
||||||
|
if (numData[elemIdx] > (int)num_regions) {
|
||||||
|
throw std::runtime_error("Values larger than maximum number of regions " + std::to_string(num_regions) + " provided in " + name);
|
||||||
|
} else if (numData[elemIdx] > 0) {
|
||||||
|
numbers[elemIdx] = static_cast<T>(numData[elemIdx]) - 1;
|
||||||
|
} else {
|
||||||
|
throw std::runtime_error("zero or negative values provided for region array: " + name);
|
||||||
|
}
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
template<class GridView, class FluidSystem, class Scalar>
|
||||||
|
void EclGenericProblem<GridView,FluidSystem,Scalar>::
|
||||||
|
updatePvtnum_()
|
||||||
|
{
|
||||||
|
const auto num_regions = eclState_.getTableManager().getTabdims().getNumPVTTables();
|
||||||
|
updateNum("PVTNUM", pvtnum_, num_regions);
|
||||||
|
}
|
||||||
|
|
||||||
|
template<class GridView, class FluidSystem, class Scalar>
|
||||||
|
void EclGenericProblem<GridView,FluidSystem,Scalar>::
|
||||||
|
updateSatnum_()
|
||||||
|
{
|
||||||
|
const auto num_regions = eclState_.getTableManager().getTabdims().getNumSatTables();
|
||||||
|
updateNum("SATNUM", satnum_, num_regions);
|
||||||
|
}
|
||||||
|
|
||||||
|
template<class GridView, class FluidSystem, class Scalar>
|
||||||
|
void EclGenericProblem<GridView,FluidSystem,Scalar>::
|
||||||
|
updateMiscnum_()
|
||||||
|
{
|
||||||
|
const auto num_regions = 1; // we only support single region
|
||||||
|
updateNum("MISCNUM", miscnum_, num_regions);
|
||||||
|
}
|
||||||
|
|
||||||
|
template<class GridView, class FluidSystem, class Scalar>
|
||||||
|
void EclGenericProblem<GridView,FluidSystem,Scalar>::
|
||||||
|
updatePlmixnum_()
|
||||||
|
{
|
||||||
|
const auto num_regions = 1; // we only support single region
|
||||||
|
updateNum("PLMIXNUM", plmixnum_, num_regions);
|
||||||
|
}
|
||||||
|
|
||||||
|
template<class GridView, class FluidSystem, class Scalar>
|
||||||
|
void EclGenericProblem<GridView,FluidSystem,Scalar>::
|
||||||
|
updateKrnum_()
|
||||||
|
{
|
||||||
|
const auto num_regions = eclState_.getTableManager().getTabdims().getNumSatTables();
|
||||||
|
updateNum("KRNUMX", krnumx_, num_regions);
|
||||||
|
updateNum("KRNUMY", krnumy_, num_regions);
|
||||||
|
updateNum("KRNUMZ", krnumz_, num_regions);
|
||||||
|
}
|
||||||
|
|
||||||
|
template<class GridView, class FluidSystem, class Scalar>
|
||||||
|
void EclGenericProblem<GridView,FluidSystem,Scalar>::
|
||||||
|
updateImbnum_()
|
||||||
|
{
|
||||||
|
const auto num_regions = eclState_.getTableManager().getTabdims().getNumSatTables();
|
||||||
|
updateNum("IMBNUMX", imbnumx_, num_regions);
|
||||||
|
updateNum("IMBNUMY", imbnumy_, num_regions);
|
||||||
|
updateNum("IMBNUMZ", imbnumz_, num_regions);
|
||||||
|
}
|
||||||
|
|
||||||
|
template<class GridView, class FluidSystem, class Scalar>
|
||||||
|
bool EclGenericProblem<GridView,FluidSystem,Scalar>::
|
||||||
|
vapparsActive(int episodeIdx) const
|
||||||
|
{
|
||||||
|
const auto& oilVaporizationControl = schedule_[episodeIdx].oilvap();
|
||||||
|
return (oilVaporizationControl.getType() == OilVaporizationProperties::OilVaporization::VAPPARS);
|
||||||
|
}
|
||||||
|
|
||||||
|
template<class GridView, class FluidSystem, class Scalar>
|
||||||
|
bool EclGenericProblem<GridView,FluidSystem,Scalar>::
|
||||||
|
drsdtActive_(int episodeIdx) const
|
||||||
|
{
|
||||||
|
const auto& oilVaporizationControl = schedule_[episodeIdx].oilvap();
|
||||||
|
const bool bothOilGasActive = FluidSystem::phaseIsActive(FluidSystem::oilPhaseIdx) &&
|
||||||
|
FluidSystem::phaseIsActive(FluidSystem::gasPhaseIdx);
|
||||||
|
return (oilVaporizationControl.drsdtActive() && bothOilGasActive);
|
||||||
|
}
|
||||||
|
|
||||||
|
template<class GridView, class FluidSystem, class Scalar>
|
||||||
|
bool EclGenericProblem<GridView,FluidSystem,Scalar>::
|
||||||
|
drvdtActive_(int episodeIdx) const
|
||||||
|
{
|
||||||
|
const auto& oilVaporizationControl = schedule_[episodeIdx].oilvap();
|
||||||
|
const bool bothOilGasActive = FluidSystem::phaseIsActive(FluidSystem::oilPhaseIdx) &&
|
||||||
|
FluidSystem::phaseIsActive(FluidSystem::gasPhaseIdx);
|
||||||
|
return (oilVaporizationControl.drvdtActive() && bothOilGasActive);
|
||||||
|
}
|
||||||
|
|
||||||
|
template<class GridView, class FluidSystem, class Scalar>
|
||||||
|
bool EclGenericProblem<GridView,FluidSystem,Scalar>::
|
||||||
|
drsdtConvective_(int episodeIdx) const
|
||||||
|
{
|
||||||
|
const auto& oilVaporizationControl = schedule_[episodeIdx].oilvap();
|
||||||
|
const bool bothOilGasActive = FluidSystem::phaseIsActive(FluidSystem::oilPhaseIdx) &&
|
||||||
|
FluidSystem::phaseIsActive(FluidSystem::gasPhaseIdx);
|
||||||
|
return (oilVaporizationControl.drsdtConvective() && bothOilGasActive);
|
||||||
|
}
|
||||||
|
|
||||||
|
template<class GridView, class FluidSystem, class Scalar>
|
||||||
|
bool EclGenericProblem<GridView,FluidSystem,Scalar>::
|
||||||
|
beginEpisode_(bool enableExperiments,
|
||||||
|
int episodeIdx)
|
||||||
|
{
|
||||||
|
if (enableExperiments && gridView_.comm().rank() == 0 && episodeIdx >= 0) {
|
||||||
|
// print some useful information in experimental mode. (the production
|
||||||
|
// simulator does this externally.)
|
||||||
|
std::ostringstream ss;
|
||||||
|
boost::posix_time::time_facet* facet = new boost::posix_time::time_facet("%d-%b-%Y");
|
||||||
|
boost::posix_time::ptime curDateTime =
|
||||||
|
boost::posix_time::from_time_t(schedule_.simTime(episodeIdx));
|
||||||
|
ss.imbue(std::locale(std::locale::classic(), facet));
|
||||||
|
ss << "Report step " << episodeIdx + 1
|
||||||
|
<< "/" << schedule_.size() - 1
|
||||||
|
<< " at day " << schedule_.seconds(episodeIdx)/(24*3600)
|
||||||
|
<< "/" << schedule_.seconds(schedule_.size() - 1)/(24*3600)
|
||||||
|
<< ", date = " << curDateTime.date()
|
||||||
|
<< "\n ";
|
||||||
|
OpmLog::info(ss.str());
|
||||||
|
}
|
||||||
|
|
||||||
|
const auto& events = schedule_[episodeIdx].events();
|
||||||
|
|
||||||
|
// react to TUNING changes
|
||||||
|
if (episodeIdx > 0 && enableTuning_ && events.hasEvent(ScheduleEvents::TUNING_CHANGE))
|
||||||
|
{
|
||||||
|
const auto& sched_state = schedule_[episodeIdx];
|
||||||
|
const auto& tuning = sched_state.tuning();
|
||||||
|
initialTimeStepSize_ = sched_state.max_next_tstep();
|
||||||
|
maxTimeStepAfterWellEvent_ = tuning.TMAXWC;
|
||||||
|
maxTimeStepSize_ = tuning.TSMAXZ;
|
||||||
|
restartShrinkFactor_ = 1./tuning.TSFCNV;
|
||||||
|
minTimeStepSize_ = tuning.TSMINZ;
|
||||||
|
return true;
|
||||||
|
}
|
||||||
|
|
||||||
|
return false;
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
template<class GridView, class FluidSystem, class Scalar>
|
||||||
|
void EclGenericProblem<GridView,FluidSystem,Scalar>::
|
||||||
|
beginTimeStep_(bool enableExperiments,
|
||||||
|
int episodeIdx,
|
||||||
|
int timeStepIndex,
|
||||||
|
Scalar startTime,
|
||||||
|
Scalar time,
|
||||||
|
Scalar timeStepSize,
|
||||||
|
Scalar endTime)
|
||||||
|
{
|
||||||
|
if (enableExperiments && gridView_.comm().rank() == 0 && episodeIdx >= 0) {
|
||||||
|
std::ostringstream ss;
|
||||||
|
boost::posix_time::time_facet* facet = new boost::posix_time::time_facet("%d-%b-%Y");
|
||||||
|
boost::posix_time::ptime date = boost::posix_time::from_time_t(startTime) + boost::posix_time::milliseconds(static_cast<long long>(time / prefix::milli));
|
||||||
|
ss.imbue(std::locale(std::locale::classic(), facet));
|
||||||
|
ss <<"\nTime step " << timeStepIndex << ", stepsize "
|
||||||
|
<< unit::convert::to(timeStepSize, unit::day) << " days,"
|
||||||
|
<< " at day " << (double)unit::convert::to(time, unit::day)
|
||||||
|
<< "/" << (double)unit::convert::to(endTime, unit::day)
|
||||||
|
<< ", date = " << date;
|
||||||
|
OpmLog::info(ss.str());
|
||||||
|
}
|
||||||
|
|
||||||
|
// update explicit quantities between timesteps.
|
||||||
|
const auto& oilVaporizationControl = schedule_[episodeIdx].oilvap();
|
||||||
|
if (drsdtActive_(episodeIdx))
|
||||||
|
// DRSDT is enabled
|
||||||
|
for (size_t pvtRegionIdx = 0; pvtRegionIdx < maxDRs_.size(); ++pvtRegionIdx)
|
||||||
|
maxDRs_[pvtRegionIdx] = oilVaporizationControl.getMaxDRSDT(pvtRegionIdx)*timeStepSize;
|
||||||
|
|
||||||
|
if (drvdtActive_(episodeIdx))
|
||||||
|
// DRVDT is enabled
|
||||||
|
for (size_t pvtRegionIdx = 0; pvtRegionIdx < maxDRv_.size(); ++pvtRegionIdx)
|
||||||
|
maxDRv_[pvtRegionIdx] = oilVaporizationControl.getMaxDRVDT(pvtRegionIdx)*timeStepSize;
|
||||||
|
}
|
||||||
|
|
||||||
|
template<class GridView, class FluidSystem, class Scalar>
|
||||||
|
void EclGenericProblem<GridView,FluidSystem,Scalar>::
|
||||||
|
initFluidSystem_()
|
||||||
|
{
|
||||||
|
FluidSystem::initFromState(eclState_, schedule_);
|
||||||
|
}
|
||||||
|
|
||||||
|
template<class GridView, class FluidSystem, class Scalar>
|
||||||
|
void EclGenericProblem<GridView,FluidSystem,Scalar>::
|
||||||
|
readBlackoilExtentionsInitialConditions_(size_t numDof,
|
||||||
|
bool enableSolvent,
|
||||||
|
bool enablePolymer,
|
||||||
|
bool enablePolymerMolarWeight,
|
||||||
|
bool enableMICP)
|
||||||
|
{
|
||||||
|
if (enableSolvent) {
|
||||||
|
if (eclState_.fieldProps().has_double("SSOL"))
|
||||||
|
solventSaturation_ = eclState_.fieldProps().get_double("SSOL");
|
||||||
|
else
|
||||||
|
solventSaturation_.resize(numDof, 0.0);
|
||||||
|
}
|
||||||
|
|
||||||
|
if (enablePolymer) {
|
||||||
|
if (eclState_.fieldProps().has_double("SPOLY"))
|
||||||
|
polymerConcentration_ = eclState_.fieldProps().get_double("SPOLY");
|
||||||
|
else
|
||||||
|
polymerConcentration_.resize(numDof, 0.0);
|
||||||
|
}
|
||||||
|
|
||||||
|
if (enablePolymerMolarWeight) {
|
||||||
|
if (eclState_.fieldProps().has_double("SPOLYMW"))
|
||||||
|
polymerMoleWeight_ = eclState_.fieldProps().get_double("SPOLYMW");
|
||||||
|
else
|
||||||
|
polymerMoleWeight_.resize(numDof, 0.0);
|
||||||
|
}
|
||||||
|
|
||||||
|
if (enableMICP) {
|
||||||
|
if (eclState_.fieldProps().has_double("SMICR"))
|
||||||
|
microbialConcentration_ = eclState_.fieldProps().get_double("SMICR");
|
||||||
|
else
|
||||||
|
microbialConcentration_.resize(numDof, 0.0);
|
||||||
|
if (eclState_.fieldProps().has_double("SOXYG"))
|
||||||
|
oxygenConcentration_ = eclState_.fieldProps().get_double("SOXYG");
|
||||||
|
else
|
||||||
|
oxygenConcentration_.resize(numDof, 0.0);
|
||||||
|
if (eclState_.fieldProps().has_double("SUREA"))
|
||||||
|
ureaConcentration_ = eclState_.fieldProps().get_double("SUREA");
|
||||||
|
else
|
||||||
|
ureaConcentration_.resize(numDof, 0.0);
|
||||||
|
if (eclState_.fieldProps().has_double("SBIOF"))
|
||||||
|
biofilmConcentration_ = eclState_.fieldProps().get_double("SBIOF");
|
||||||
|
else
|
||||||
|
biofilmConcentration_.resize(numDof, 0.0);
|
||||||
|
if (eclState_.fieldProps().has_double("SCALC"))
|
||||||
|
calciteConcentration_ = eclState_.fieldProps().get_double("SCALC");
|
||||||
|
else
|
||||||
|
calciteConcentration_.resize(numDof, 0.0);
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
template<class GridView, class FluidSystem, class Scalar>
|
||||||
|
Scalar EclGenericProblem<GridView,FluidSystem,Scalar>::
|
||||||
|
maxWaterSaturation(unsigned globalDofIdx) const
|
||||||
|
{
|
||||||
|
if (maxWaterSaturation_.empty())
|
||||||
|
return 0.0;
|
||||||
|
|
||||||
|
return maxWaterSaturation_[globalDofIdx];
|
||||||
|
}
|
||||||
|
|
||||||
|
template<class GridView, class FluidSystem, class Scalar>
|
||||||
|
Scalar EclGenericProblem<GridView,FluidSystem,Scalar>::
|
||||||
|
minOilPressure(unsigned globalDofIdx) const
|
||||||
|
{
|
||||||
|
if (minOilPressure_.empty())
|
||||||
|
return 0.0;
|
||||||
|
|
||||||
|
return minOilPressure_[globalDofIdx];
|
||||||
|
}
|
||||||
|
|
||||||
|
template<class GridView, class FluidSystem, class Scalar>
|
||||||
|
Scalar EclGenericProblem<GridView,FluidSystem,Scalar>::
|
||||||
|
overburdenPressure(unsigned elementIdx) const
|
||||||
|
{
|
||||||
|
if (overburdenPressure_.empty())
|
||||||
|
return 0.0;
|
||||||
|
|
||||||
|
return overburdenPressure_[elementIdx];
|
||||||
|
}
|
||||||
|
|
||||||
|
template<class GridView, class FluidSystem, class Scalar>
|
||||||
|
Scalar EclGenericProblem<GridView,FluidSystem,Scalar>::
|
||||||
|
solventSaturation(unsigned elemIdx) const
|
||||||
|
{
|
||||||
|
if (solventSaturation_.empty())
|
||||||
|
return 0;
|
||||||
|
|
||||||
|
return solventSaturation_[elemIdx];
|
||||||
|
}
|
||||||
|
|
||||||
|
template<class GridView, class FluidSystem, class Scalar>
|
||||||
|
Scalar EclGenericProblem<GridView,FluidSystem,Scalar>::
|
||||||
|
drsdtcon(unsigned elemIdx, int episodeIdx) const
|
||||||
|
{
|
||||||
|
if (convectiveDrs_.empty())
|
||||||
|
return 0;
|
||||||
|
|
||||||
|
// The episode index is set to -1 in the initialization phase.
|
||||||
|
// Output drsdt value for index 0
|
||||||
|
episodeIdx = std::max(episodeIdx, 0);
|
||||||
|
const auto& oilVaporizationControl = schedule_[episodeIdx].oilvap();
|
||||||
|
int pvtRegionIdx = pvtRegionIndex(elemIdx);
|
||||||
|
return oilVaporizationControl.getMaxDRSDT(pvtRegionIdx)*convectiveDrs_[elemIdx];
|
||||||
|
}
|
||||||
|
|
||||||
|
template<class GridView, class FluidSystem, class Scalar>
|
||||||
|
Scalar EclGenericProblem<GridView,FluidSystem,Scalar>::
|
||||||
|
polymerConcentration(unsigned elemIdx) const
|
||||||
|
{
|
||||||
|
if (polymerConcentration_.empty())
|
||||||
|
return 0;
|
||||||
|
|
||||||
|
return polymerConcentration_[elemIdx];
|
||||||
|
}
|
||||||
|
|
||||||
|
template<class GridView, class FluidSystem, class Scalar>
|
||||||
|
Scalar EclGenericProblem<GridView,FluidSystem,Scalar>::
|
||||||
|
polymerMolecularWeight(const unsigned elemIdx) const
|
||||||
|
{
|
||||||
|
if (polymerMoleWeight_.empty())
|
||||||
|
return 0.0;
|
||||||
|
|
||||||
|
return polymerMoleWeight_[elemIdx];
|
||||||
|
}
|
||||||
|
|
||||||
|
template<class GridView, class FluidSystem, class Scalar>
|
||||||
|
Scalar EclGenericProblem<GridView,FluidSystem,Scalar>::
|
||||||
|
microbialConcentration(unsigned elemIdx) const
|
||||||
|
{
|
||||||
|
if (microbialConcentration_.empty())
|
||||||
|
return 0;
|
||||||
|
|
||||||
|
return microbialConcentration_[elemIdx];
|
||||||
|
}
|
||||||
|
|
||||||
|
template<class GridView, class FluidSystem, class Scalar>
|
||||||
|
Scalar EclGenericProblem<GridView,FluidSystem,Scalar>::
|
||||||
|
oxygenConcentration(unsigned elemIdx) const
|
||||||
|
{
|
||||||
|
if (oxygenConcentration_.empty())
|
||||||
|
return 0;
|
||||||
|
|
||||||
|
return oxygenConcentration_[elemIdx];
|
||||||
|
}
|
||||||
|
|
||||||
|
template<class GridView, class FluidSystem, class Scalar>
|
||||||
|
Scalar EclGenericProblem<GridView,FluidSystem,Scalar>::
|
||||||
|
ureaConcentration(unsigned elemIdx) const
|
||||||
|
{
|
||||||
|
if (ureaConcentration_.empty())
|
||||||
|
return 0;
|
||||||
|
|
||||||
|
return ureaConcentration_[elemIdx];
|
||||||
|
}
|
||||||
|
|
||||||
|
template<class GridView, class FluidSystem, class Scalar>
|
||||||
|
Scalar EclGenericProblem<GridView,FluidSystem,Scalar>::
|
||||||
|
biofilmConcentration(unsigned elemIdx) const
|
||||||
|
{
|
||||||
|
if (biofilmConcentration_.empty())
|
||||||
|
return 0;
|
||||||
|
|
||||||
|
return biofilmConcentration_[elemIdx];
|
||||||
|
}
|
||||||
|
|
||||||
|
template<class GridView, class FluidSystem, class Scalar>
|
||||||
|
Scalar EclGenericProblem<GridView,FluidSystem,Scalar>::
|
||||||
|
calciteConcentration(unsigned elemIdx) const
|
||||||
|
{
|
||||||
|
if (calciteConcentration_.empty())
|
||||||
|
return 0;
|
||||||
|
|
||||||
|
return calciteConcentration_[elemIdx];
|
||||||
|
}
|
||||||
|
|
||||||
|
template<class GridView, class FluidSystem, class Scalar>
|
||||||
|
unsigned EclGenericProblem<GridView,FluidSystem,Scalar>::
|
||||||
|
pvtRegionIndex(unsigned elemIdx) const
|
||||||
|
{
|
||||||
|
if (pvtnum_.empty())
|
||||||
|
return 0;
|
||||||
|
|
||||||
|
return pvtnum_[elemIdx];
|
||||||
|
}
|
||||||
|
|
||||||
|
template<class GridView, class FluidSystem, class Scalar>
|
||||||
|
unsigned EclGenericProblem<GridView,FluidSystem,Scalar>::
|
||||||
|
satnumRegionIndex(unsigned elemIdx) const
|
||||||
|
{
|
||||||
|
if (satnum_.empty())
|
||||||
|
return 0;
|
||||||
|
|
||||||
|
return satnum_[elemIdx];
|
||||||
|
}
|
||||||
|
|
||||||
|
template<class GridView, class FluidSystem, class Scalar>
|
||||||
|
unsigned EclGenericProblem<GridView,FluidSystem,Scalar>::
|
||||||
|
miscnumRegionIndex(unsigned elemIdx) const
|
||||||
|
{
|
||||||
|
if (miscnum_.empty())
|
||||||
|
return 0;
|
||||||
|
|
||||||
|
return miscnum_[elemIdx];
|
||||||
|
}
|
||||||
|
|
||||||
|
template<class GridView, class FluidSystem, class Scalar>
|
||||||
|
unsigned EclGenericProblem<GridView,FluidSystem,Scalar>::
|
||||||
|
plmixnumRegionIndex(unsigned elemIdx) const
|
||||||
|
{
|
||||||
|
if (plmixnum_.empty())
|
||||||
|
return 0;
|
||||||
|
|
||||||
|
return plmixnum_[elemIdx];
|
||||||
|
}
|
||||||
|
|
||||||
|
template<class GridView, class FluidSystem, class Scalar>
|
||||||
|
Scalar EclGenericProblem<GridView,FluidSystem,Scalar>::
|
||||||
|
maxPolymerAdsorption(unsigned elemIdx) const
|
||||||
|
{
|
||||||
|
if (maxPolymerAdsorption_.empty())
|
||||||
|
return 0;
|
||||||
|
|
||||||
|
return maxPolymerAdsorption_[elemIdx];
|
||||||
|
}
|
||||||
|
|
||||||
|
template<class GridView, class FluidSystem, class Scalar>
|
||||||
|
void EclGenericProblem<GridView,FluidSystem,Scalar>::
|
||||||
|
initDRSDT_(size_t numDof,
|
||||||
|
int episodeIdx)
|
||||||
|
{
|
||||||
|
// deal with DRSDT
|
||||||
|
unsigned ntpvt = eclState_.runspec().tabdims().getNumPVTTables();
|
||||||
|
//TODO We may want to only allocate these properties only if active.
|
||||||
|
//But since they may be activated at later time we need some more
|
||||||
|
//intrastructure to handle it
|
||||||
|
if (FluidSystem::phaseIsActive(FluidSystem::oilPhaseIdx) && FluidSystem::phaseIsActive(FluidSystem::gasPhaseIdx)) {
|
||||||
|
maxDRv_.resize(ntpvt, 1e30);
|
||||||
|
lastRv_.resize(numDof, 0.0);
|
||||||
|
maxDRs_.resize(ntpvt, 1e30);
|
||||||
|
dRsDtOnlyFreeGas_.resize(ntpvt, false);
|
||||||
|
lastRs_.resize(numDof, 0.0);
|
||||||
|
maxDRv_.resize(ntpvt, 1e30);
|
||||||
|
lastRv_.resize(numDof, 0.0);
|
||||||
|
maxOilSaturation_.resize(numDof, 0.0);
|
||||||
|
if (drsdtConvective_(episodeIdx)) {
|
||||||
|
convectiveDrs_.resize(numDof, 1.0);
|
||||||
|
}
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
template<class GridView, class FluidSystem, class Scalar>
|
||||||
|
bool EclGenericProblem<GridView,FluidSystem,Scalar>::
|
||||||
|
operator==(const EclGenericProblem& rhs) const
|
||||||
|
{
|
||||||
|
return this->maxOilSaturation_ == rhs.maxOilSaturation_ &&
|
||||||
|
this->maxPolymerAdsorption_ == rhs.maxPolymerAdsorption_ &&
|
||||||
|
this->maxWaterSaturation_ == rhs.maxWaterSaturation_ &&
|
||||||
|
this->minOilPressure_ == rhs.minOilPressure_ &&
|
||||||
|
this->overburdenPressure_ == rhs.overburdenPressure_ &&
|
||||||
|
this->polymerConcentration_ == rhs.polymerConcentration_ &&
|
||||||
|
this->polymerMoleWeight_ == rhs.polymerMoleWeight_ &&
|
||||||
|
this->solventSaturation_ == rhs.solventSaturation_ &&
|
||||||
|
this->microbialConcentration_ == rhs.microbialConcentration_ &&
|
||||||
|
this->oxygenConcentration_ == rhs.oxygenConcentration_ &&
|
||||||
|
this->ureaConcentration_ == rhs.ureaConcentration_ &&
|
||||||
|
this->biofilmConcentration_ == rhs.biofilmConcentration_ &&
|
||||||
|
this->calciteConcentration_ == rhs.calciteConcentration_ &&
|
||||||
|
this->lastRv_ == rhs.lastRv_ &&
|
||||||
|
this->maxDRv_ == rhs.maxDRv_ &&
|
||||||
|
this->convectiveDrs_ == rhs.convectiveDrs_ &&
|
||||||
|
this->lastRs_ == rhs.lastRs_ &&
|
||||||
|
this->maxDRs_ == rhs.maxDRs_ &&
|
||||||
|
this->dRsDtOnlyFreeGas_ == rhs.dRsDtOnlyFreeGas_;
|
||||||
|
}
|
||||||
|
|
||||||
|
} // namespace Opm
|
||||||
|
|
||||||
|
#endif // EWOMS_GENERIC_ECL_PROBLEM_IMPL_HH
|
||||||
@ -20,6 +20,7 @@
|
|||||||
#include <opm/simulators/flow/Main.hpp>
|
#include <opm/simulators/flow/Main.hpp>
|
||||||
#include <dune/alugrid/grid.hh>
|
#include <dune/alugrid/grid.hh>
|
||||||
#include <ebos/eclalugridvanguard.hh>
|
#include <ebos/eclalugridvanguard.hh>
|
||||||
|
#include <ebos/eclgenericproblem_impl.hh>
|
||||||
namespace Opm {
|
namespace Opm {
|
||||||
namespace Properties {
|
namespace Properties {
|
||||||
namespace TTag {
|
namespace TTag {
|
||||||
|
|||||||
@ -23,7 +23,7 @@
|
|||||||
#include <ebos/equil/initstateequil_impl.hh>
|
#include <ebos/equil/initstateequil_impl.hh>
|
||||||
#include <opm/models/blackoil/blackoillocalresidualtpfa.hh>
|
#include <opm/models/blackoil/blackoillocalresidualtpfa.hh>
|
||||||
#include <opm/models/discretization/common/tpfalinearizer.hh>
|
#include <opm/models/discretization/common/tpfalinearizer.hh>
|
||||||
|
#include <ebos/eclgenericproblem_impl.hh>
|
||||||
namespace Opm {
|
namespace Opm {
|
||||||
namespace Properties {
|
namespace Properties {
|
||||||
namespace TTag {
|
namespace TTag {
|
||||||
|
|||||||
Loading…
Reference in New Issue
Block a user