OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2024-11-29 12:33:49 -06:00
Code Issues Packages Projects Releases Wiki Activity

fix initialization EQUIL + water and gas

Browse Source
This commit is contained in:
Tor Harald Sandve 2022-03-11 10:12:14 +01:00
parent ae32b44124
commit 82099705ac
1 changed files with 16 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

18
ebos/equil/initstateequil.hh
View File

@ -626,13 +626,20 @@ equil_WOG(const Region& reg, const VSpan& span)
this->makeOilPressure(ic, reg, span); this->makeOilPressure(ic, reg, span);
} }
if (this->gasActive()) { if (this->gasActive() && this->oilActive()) {
// Pcgo = Pg - Po => Pg = Po + Pcgo // Pcgo = Pg - Po => Pg = Po + Pcgo
const auto ic = typename GPress::InitCond { const auto ic = typename GPress::InitCond {
reg.zgoc(), reg.zgoc(),
this->oil(reg.zgoc()) + reg.pcgoGoc() this->oil(reg.zgoc()) + reg.pcgoGoc()
}; };
this->makeGasPressure(ic, reg, span);
} else if (this->gasActive() && !this->oilActive()) {
// No oil phase set Pg = Pw + Pcgw
const auto ic = typename GPress::InitCond {
reg.zwoc(), // The WOC is really the GWC for gas/water cases
this->water(reg.zwoc()) + reg.pcowWoc() // Pcow(WOC) is really Pcgw(GWC) for gas/water cases
};
this->makeGasPressure(ic, reg, span); this->makeGasPressure(ic, reg, span);
} }
} }
@ -669,13 +676,20 @@ equil_GOW(const Region& reg, const VSpan& span)
this->makeOilPressure(ic, reg, span); this->makeOilPressure(ic, reg, span);
} }
if (this->waterActive()) { if (this->waterActive() && this->oilActive()) {
// Pcow = Po - Pw => Pw = Po - Pcow // Pcow = Po - Pw => Pw = Po - Pcow
const auto ic = typename WPress::InitCond { const auto ic = typename WPress::InitCond {
reg.zwoc(), reg.zwoc(),
this->oil(reg.zwoc()) - reg.pcowWoc() this->oil(reg.zwoc()) - reg.pcowWoc()
}; };
this->makeWatPressure(ic, reg, span);
} else if (this->waterActive() && !this->oilActive()) {
// No oil phase set Pw = Pg - Pcgw
const auto ic = typename WPress::InitCond {
reg.zwoc(), // The WOC is really the GWC for gas/water cases
this->gas(reg.zwoc()) - reg.pcowWoc() // Pcow(WOC) is really Pcgw(GWC) for gas/water cases
};
this->makeWatPressure(ic, reg, span); this->makeWatPressure(ic, reg, span);
} }
} }