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

Added output of dew and bubble point pressures

Browse Source
This commit is contained in:
babrodtk 2017-02-15 14:27:32 +01:00
parent dda6f58de0
commit 4a46451d95
2 changed files with 17 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
opm/autodiff/BlackoilModelEbos.hpp
View File

@ -1271,6 +1271,13 @@ namespace Opm {
VectorType& RsSat = simData.getCellData( "RSSAT" ); VectorType& RsSat = simData.getCellData( "RSSAT" );
VectorType& RvSat = simData.getCellData( "RVSAT" ); VectorType& RvSat = simData.getCellData( "RVSAT" );
simData.registerCellData( "PBUB", 1 );
simData.registerCellData( "PDEW", 1 );
VectorType& Pb = simData.getCellData( "PBUB" );
VectorType& Pd = simData.getCellData( "PDEW" );
for (int cellIdx = 0; cellIdx < numCells; ++cellIdx) { for (int cellIdx = 0; cellIdx < numCells; ++cellIdx) {
const auto& intQuants = *ebosModel.cachedIntensiveQuantities(cellIdx, /*timeIdx=*/0); const auto& intQuants = *ebosModel.cachedIntensiveQuantities(cellIdx, /*timeIdx=*/0);
const auto& fs = intQuants.fluidState(); const auto& fs = intQuants.fluidState();
@ -1304,6 +1311,8 @@ namespace Opm {
FluidSystem::gasPhaseIdx, FluidSystem::gasPhaseIdx,
intQuants.pvtRegionIndex(), intQuants.pvtRegionIndex(),
/*maxOilSaturation=*/1.0).value(); /*maxOilSaturation=*/1.0).value();
Pb[cellIdx] = FluidSystem::bubblePointPressure(fs, intQuants.pvtRegionIndex()).value();
Pd[cellIdx] = FluidSystem::dewPointPressure(fs, intQuants.pvtRegionIndex()).value();
} }
if( liquid_active ) if( liquid_active )
{ {

11
opm/autodiff/SimulatorFullyImplicitBlackoilOutput.hpp
View File

@ -723,9 +723,14 @@ namespace Opm
if (log && vapour_active && if (log && vapour_active &&
liquid_active && rstKeywords["PBPD"] > 0) { liquid_active && rstKeywords["PBPD"] > 0) {
rstKeywords["PBPD"] = 0; rstKeywords["PBPD"] = 0;
Opm::OpmLog::warning("Bubble/dew point pressure output unsupported", output.insert("PBUB",
"Writing bubble points and dew points (PBPD) to file is unsupported, " Opm::UnitSystem::measure::pressure,
"as the simulator does not use these internally."); std::move( sd.getCellData("PBUB") ),
data::TargetType::RESTART_AUXILLARY);
output.insert("PDEW",
Opm::UnitSystem::measure::pressure,
std::move( sd.getCellData("PDEW") ),
data::TargetType::RESTART_AUXILLARY);
} }
//Warn for any unhandled keyword //Warn for any unhandled keyword