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

Merge pull request #4647 from totto82/initRVW

Browse Source 
initialize trivial rvw if no table is given
This commit is contained in:
Tor Harald Sandve 2023-05-31 15:13:08 +02:00 committed by GitHub
commit c70dea53a8
No known key found for this signature in database
GPG Key ID: 4AEE18F83AFDEB23
1 changed files with 9 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

16
ebos/equil/initstateequil_impl.hh
View File

@ -1436,14 +1436,16 @@ InitialStateComputer(MaterialLawManager& materialLawManager,
} }
else { else {
if (rec[i].gasOilContactDepth() != rec[i].datumDepth()) { if (rec[i].gasOilContactDepth() != rec[i].datumDepth()) {
throw std::runtime_error( rvwFunc_.push_back(std::make_shared<Miscibility::NoMixing>());
"Cannot initialise: when no explicit RVWVD table is given, \n" const auto msg = "No explicit RVWVD table is given for EQUIL region " + std::to_string(i + 1) +". \n"
"datum depth must be at the gas-oil-contact. " "and datum depth is not at the gas-water-contact. \n"
"In EQUIL region "+std::to_string(i + 1)+" (counting from 1), this does not hold."); "Rvw is set to 0.0 in all cells. \n";
OpmLog::warning(msg);
} else {
const double pContact = rec[i].datumDepthPressure() + rec[i].gasOilContactCapillaryPressure();
const double TContact = 273.15 + 20; // standard temperature for now
rvwFunc_.push_back(std::make_shared<Miscibility::RvwSatAtContact<FluidSystem>>(pvtIdx,pContact, TContact));
} }
const double pContact = rec[i].datumDepthPressure() + rec[i].gasOilContactCapillaryPressure();
const double TContact = 273.15 + 20; // standard temperature for now
rvwFunc_.push_back(std::make_shared<Miscibility::RvwSatAtContact<FluidSystem>>(pvtIdx,pContact, TContact));
} }
} }
} }