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

Fixed bug in SinglePointUpwindTwoPhasePolymer.

Browse Source
This commit is contained in:
Xavier Raynaud
2012-04-24 11:20:06 +02:00
parent 81d90b28e2
commit 6eb68c9652
2 changed files with 5 additions and 5 deletions
Download Patch File Download Diff File Expand all files Collapse all files

6
opm/polymer/PolymerProperties.cpp
View File

@@ -238,7 +238,7 @@ namespace Opm
double& dmobwatdc) const
{
effectiveMobilitiesBoth(c, cmax, visc,
relperm, drelpermds, mob,dmobds,
relperm, drelpermds, mob, dmobds,
dmobwatdc, true);
}
@@ -271,8 +271,8 @@ namespace Opm
dmobwat_dc = - eff_relperm_wat*dmu_w_eff_dc/(mu_w_eff*mu_w_eff)
+ deff_relperm_wat_dc/mu_w_eff;
dmob_ds[0*2 + 0] = deff_relperm_wat_ds/visc_eff[0];
dmob_ds[0*2 + 1] = -dmob_ds[0*2 + 0];
dmob_ds[1*2 + 0] = drelperm_ds[1*2 + 0]/visc_eff[1];
dmob_ds[0*2 + 1] = drelperm_ds[0*2 + 1]/visc_eff[1];
dmob_ds[1*2 + 0] = drelperm_ds[1*2 + 0]/visc_eff[0];
dmob_ds[1*2 + 1] = drelperm_ds[1*2 + 1]/visc_eff[1];
} else {
dmob_ds.clear();

4
opm/polymer/SinglePointUpwindTwoPhasePolymer.hpp
View File

@@ -519,14 +519,14 @@ namespace Opm {
store_.ds(cell) = x[2*cell + 0];
s[0] = sat[cell*2 + 0] + x[2*cell + 0];
c = cpoly[cell] + x[2*cell + 1];
cmax = std::max(cpoly[cell] + x[2*cell + 1], cmaxpoly[cell]);
store_.sw(cell) = s[0];
store_.c(cell) = c;
cmax = std::max(c, cmaxpoly[cell]);
store_.cmax(cell) = cmax;
store_.dsc(cell) = s[0]*c - sat[cell*2 + 0]*cpoly[cell];
double dcadsdc;
double cads;
fluid_.adsorption(cpoly[cell], cmax, cads, dcadsdc);
fluid_.adsorption(cpoly[cell], cmaxpoly[cell], cads, dcadsdc);
store_.dcads(cell) = -cads;
fluid_.adsorption(c, cmax, cads, dcadsdc);
store_.dcads(cell) += cads;