OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-25 18:55:30 -06:00
Code Issues Packages Projects Releases Wiki Activity

Add support for region DR[SV]DTR

Browse Source
This commit is contained in:
Tor Harald Sandve 2018-10-15 12:20:39 +02:00
parent 3185e8662b
commit fad7c68446
1 changed files with 63 additions and 33 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

96
ebos/eclproblem.hh
View File

@ -538,30 +538,56 @@ public:
} }
// deal with DRSDT // deal with DRSDT
maxDRsDt_ = 0.0; const auto& eclState = simulator.vanguard().eclState();
maxDRs_ = -1.0; unsigned ntpvt = eclState.runspec().tabdims().getNumPVTTables();
if (!vapparsActive_ && deck.hasKeyword("DRSDT")) { maxDRsDt_.resize(ntpvt, 0.0);
drsdtActive_ = !vapparsActive_; maxDRs_.resize(ntpvt, -1.0);
const auto& drsdtKeyword = deck.getKeyword("DRSDT"); dRsDtOnlyFreeGas_.resize(ntpvt, false);
maxDRsDt_ = drsdtKeyword.getRecord(0).getItem("DRSDT_MAX").getSIDouble(0); if (!vapparsActive_) {
size_t numDof = this->model().numGridDof(); if (deck.hasKeyword("DRSDT")) {
lastRs_.resize(numDof, 0.0); drsdtActive_ = !vapparsActive_;
const auto& drsdtKeyword = deck.getKeyword("DRSDT");
std::fill(maxDRsDt_.begin(), maxDRsDt_.end(), drsdtKeyword.getRecord(0).getItem("DRSDT_MAX").getSIDouble(0));
size_t numDof = this->model().numGridDof();
lastRs_.resize(numDof, 0.0);
std::string drsdtFlag =
drsdtKeyword.getRecord(0).getItem("Option").getTrimmedString(0);
std::transform(drsdtFlag.begin(), drsdtFlag.end(), drsdtFlag.begin(), ::toupper);
std::string drsdtFlag = std::fill(dRsDtOnlyFreeGas_.begin(), dRsDtOnlyFreeGas_.end(), (drsdtFlag == "FREE"));
drsdtKeyword.getRecord(0).getItem("Option").getTrimmedString(0); } else if (deck.hasKeyword("DRSDTR")) {
std::transform(drsdtFlag.begin(), drsdtFlag.end(), drsdtFlag.begin(), ::toupper); drsdtActive_ = !vapparsActive_;
const auto& drsdtrKeyword = deck.getKeyword("DRSDTR");
dRsDtOnlyFreeGas_ = (drsdtFlag == "FREE"); size_t numDof = this->model().numGridDof();
lastRs_.resize(numDof, 0.0);
for (size_t recordIdx = 0; recordIdx < drsdtrKeyword.size(); ++recordIdx) {
const Opm::DeckRecord& record = drsdtrKeyword.getRecord(recordIdx);
maxDRsDt_[recordIdx] = record.getItem("DRSDT_MAX").getSIDouble(0);
std::string drsdtFlag = record.getItem("Option").getTrimmedString(0);
std::transform(drsdtFlag.begin(), drsdtFlag.end(), drsdtFlag.begin(), ::toupper);
dRsDtOnlyFreeGas_[recordIdx] = (drsdtFlag == "FREE");
}
}
} }
// deal with DRVDT // deal with DRVDT
maxDRvDt_ = 0.0; maxDRvDt_.resize(ntpvt, 0.0);
maxDRv_ = -1.0; maxDRv_.resize(ntpvt, -1.0);
if (!vapparsActive_ && deck.hasKeyword("DRVDT")) { if (!vapparsActive_) {
const auto& drvdtKeyword = deck.getKeyword("DVSDT"); if (deck.hasKeyword("DRVDT")) {
maxDRvDt_ = drvdtKeyword.getRecord(0).getItem("DRVDT_MAX").getSIDouble(0); const auto& drvdtKeyword = deck.getKeyword("DRVDT");
size_t numDof = this->model().numGridDof(); std::fill(maxDRvDt_.begin(), maxDRvDt_.end(), drvdtKeyword.getRecord(0).getItem("DRVDT_MAX").getSIDouble(0));
lastRv_.resize(numDof, 0.0); size_t numDof = this->model().numGridDof();
lastRv_.resize(numDof, 0.0);
} else if (deck.hasKeyword("DRVDTR")) {
const auto& drvdtrKeyword = deck.getKeyword("DRVDTR");
size_t numDof = this->model().numGridDof();
lastRv_.resize(numDof, 0.0);
for (size_t recordIdx = 0; recordIdx < drvdtrKeyword.size(); ++recordIdx) {
const Opm::DeckRecord& record = drvdtrKeyword.getRecord(recordIdx);
maxDRvDt_[recordIdx] = record.getItem("DRVDT_MAX").getSIDouble(0);
}
}
} }
initFluidSystem_(); initFluidSystem_();
@ -571,7 +597,6 @@ public:
readThermalParameters_(); readThermalParameters_();
transmissibilities_.finishInit(); transmissibilities_.finishInit();
const auto& eclState = simulator.vanguard().eclState();
const auto& initconfig = eclState.getInitConfig(); const auto& initconfig = eclState.getInitConfig();
const auto& timeMap = simulator.vanguard().schedule().getTimeMap(); const auto& timeMap = simulator.vanguard().schedule().getTimeMap();
if(initconfig.restartRequested()) { if(initconfig.restartRequested()) {
@ -737,11 +762,13 @@ public:
{ {
if (drsdtActive_) if (drsdtActive_)
// DRSDT is enabled // DRSDT is enabled
maxDRs_ = maxDRsDt_*this->simulator().timeStepSize(); for (size_t pvtRegionIdx = 0; pvtRegionIdx < maxDRs_.size(); ++pvtRegionIdx )
maxDRs_[pvtRegionIdx] = maxDRsDt_[pvtRegionIdx]*this->simulator().timeStepSize();
if (drvdtActive_) if (drvdtActive_)
// DRVDT is enabled // DRVDT is enabled
maxDRv_ = maxDRvDt_*this->simulator().timeStepSize(); for (size_t pvtRegionIdx = 0; pvtRegionIdx < maxDRv_.size(); ++pvtRegionIdx )
maxDRv_[pvtRegionIdx] = maxDRvDt_[pvtRegionIdx]*this->simulator().timeStepSize();
if (!GET_PROP_VALUE(TypeTag, DisableWells)) if (!GET_PROP_VALUE(TypeTag, DisableWells))
wellModel_.beginTimeStep(); wellModel_.beginTimeStep();
@ -1354,10 +1381,11 @@ public:
*/ */
Scalar maxGasDissolutionFactor(unsigned globalDofIdx) const Scalar maxGasDissolutionFactor(unsigned globalDofIdx) const
{ {
if (!drsdtActive_ || maxDRs_ < 0.0) int pvtRegionIdx = pvtRegionIndex(globalDofIdx);
if (!drsdtActive_ || maxDRs_[pvtRegionIdx] < 0.0)
return std::numeric_limits<Scalar>::max()/2; return std::numeric_limits<Scalar>::max()/2;
return lastRs_[globalDofIdx] + maxDRs_; return lastRs_[globalDofIdx] + maxDRs_[pvtRegionIdx];
} }
/*! /*!
@ -1366,10 +1394,11 @@ public:
*/ */
Scalar maxOilVaporizationFactor(unsigned globalDofIdx) const Scalar maxOilVaporizationFactor(unsigned globalDofIdx) const
{ {
if (!drvdtActive_ || maxDRv_ < 0.0) int pvtRegionIdx = pvtRegionIndex(globalDofIdx);
if (!drvdtActive_ || maxDRv_[pvtRegionIdx] < 0.0)
return std::numeric_limits<Scalar>::max()/2; return std::numeric_limits<Scalar>::max()/2;
return lastRv_[globalDofIdx] + maxDRv_; return lastRv_[globalDofIdx] + maxDRv_[pvtRegionIdx];
} }
/*! /*!
@ -1483,7 +1512,8 @@ private:
typedef typename std::decay<decltype(fs) >::type FluidState; typedef typename std::decay<decltype(fs) >::type FluidState;
if (!dRsDtOnlyFreeGas_ || fs.saturation(gasPhaseIdx) > freeGasMinSaturation_) int pvtRegionIdx = pvtRegionIndex(compressedDofIdx);
if (!dRsDtOnlyFreeGas_[pvtRegionIdx] || fs.saturation(gasPhaseIdx) > freeGasMinSaturation_)
lastRs_[compressedDofIdx] = lastRs_[compressedDofIdx] =
Opm::BlackOil::template getRs_<FluidSystem, Opm::BlackOil::template getRs_<FluidSystem,
FluidState, FluidState,
@ -2123,14 +2153,14 @@ private:
std::vector<Scalar> solventSaturation_; std::vector<Scalar> solventSaturation_;
bool drsdtActive_; // if no, VAPPARS *might* be active bool drsdtActive_; // if no, VAPPARS *might* be active
bool dRsDtOnlyFreeGas_; // apply the DRSDT rate limit only to cells that exhibit free gas std::vector<bool> dRsDtOnlyFreeGas_; // apply the DRSDT rate limit only to cells that exhibit free gas
std::vector<Scalar> lastRs_; std::vector<Scalar> lastRs_;
Scalar maxDRsDt_; std::vector<Scalar> maxDRsDt_;
Scalar maxDRs_; std::vector<Scalar> maxDRs_;
bool drvdtActive_; // if no, VAPPARS *might* be active bool drvdtActive_; // if no, VAPPARS *might* be active
std::vector<Scalar> lastRv_; std::vector<Scalar> lastRv_;
Scalar maxDRvDt_; std::vector<Scalar> maxDRvDt_;
Scalar maxDRv_; std::vector<Scalar> maxDRv_;
constexpr static Scalar freeGasMinSaturation_ = 1e-7; constexpr static Scalar freeGasMinSaturation_ = 1e-7;
bool vapparsActive_; // if no, DRSDT and/or DRVDT *might* be active bool vapparsActive_; // if no, DRSDT and/or DRVDT *might* be active