mirror of
https://github.com/OPM/opm-simulators.git
synced 2026-09-05 04:40:19 -05:00
Merge pull request #866 from totto82/fixSolventRs
BUGFIX allow for solvent flux when disgasw
This commit is contained in:
@@ -485,12 +485,12 @@ public:
|
||||
|
||||
if (isSolubleInWater()) {
|
||||
if (upIdx == inIdx)
|
||||
flux[contiSolventEqIdx] =
|
||||
flux[contiSolventEqIdx] +=
|
||||
extQuants.volumeFlux(waterPhaseIdx)
|
||||
* up.fluidState().invB(waterPhaseIdx)
|
||||
* up.rsSolw();
|
||||
else
|
||||
flux[contiSolventEqIdx] =
|
||||
flux[contiSolventEqIdx] +=
|
||||
extQuants.volumeFlux(waterPhaseIdx)
|
||||
*decay<Scalar>(up.fluidState().invB(waterPhaseIdx))
|
||||
*decay<Scalar>(up.rsSolw());
|
||||
@@ -509,12 +509,12 @@ public:
|
||||
|
||||
if (isSolubleInWater()) {
|
||||
if (upIdx == inIdx)
|
||||
flux[contiSolventEqIdx] =
|
||||
flux[contiSolventEqIdx] +=
|
||||
extQuants.volumeFlux(waterPhaseIdx)
|
||||
* up.fluidState().density(waterPhaseIdx)
|
||||
* up.rsSolw();
|
||||
else
|
||||
flux[contiSolventEqIdx] =
|
||||
flux[contiSolventEqIdx] +=
|
||||
extQuants.volumeFlux(waterPhaseIdx)
|
||||
*decay<Scalar>(up.fluidState().density(waterPhaseIdx))
|
||||
*decay<Scalar>(up.rsSolw());
|
||||
|
||||
Reference in New Issue
Block a user