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

Also fix gas-oil initialization with transision zone

Browse Source
This commit is contained in:
Tor Harald Sandve
2024-04-19 10:12:14 +02:00
committed by Lisa Julia Nebel
parent 5abf487710
commit 4171bb0961
1 changed files with 8 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

11
opm/simulators/flow/equil/EquilibrationHelpers_impl.hpp
View File

@@ -446,10 +446,15 @@ operator()(double s) const
fluidState.setSaturation(FluidSystem::waterPhaseIdx, 0.0);
fluidState.setSaturation(FluidSystem::oilPhaseIdx, 0.0);
fluidState.setSaturation(FluidSystem::gasPhaseIdx, 0.0);
fluidState.setSaturation(phase_, s);
if (phase_ == FluidSystem::gasPhaseIdx)
fluidState.setSaturation(FluidSystem::waterPhaseIdx, 1.0 - s);
fluidState.setSaturation(phase_, s);
// This code should only be called for water or gas phase
assert(phase_ != FluidSystem::oilPhaseIdx);
// The capillaryPressure code only uses the wetting phase saturation
if (phase_ == FluidSystem::gasPhaseIdx) {
fluidState.setSaturation(FluidSystem::waterPhaseIdx, 1.0 - s);
fluidState.setSaturation(FluidSystem::oilPhaseIdx, 1.0 - s);
}
std::array<double, FluidSystem::numPhases> pc{0.0};
using MaterialLaw = typename MaterialLawManager::MaterialLaw;
MaterialLaw::capillaryPressures(pc, matParams, fluidState);