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 computation of derivatives for polymer.

Browse Source
This commit is contained in:
Xavier Raynaud
2012-06-14 09:35:45 +02:00
parent ee89824e0d
commit 9df87db47f
1 changed files with 1 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files

2
opm/polymer/PolymerProperties.cpp
View File

@@ -211,7 +211,7 @@ namespace Opm
double rk = 1 + (res_factor_ - 1)*c_ads/c_max_ads_; double rk = 1 + (res_factor_ - 1)*c_ads/c_max_ads_;
eff_relperm_wat = relperm[0]/rk; eff_relperm_wat = relperm[0]/rk;
if (if_with_der) { if (if_with_der) {
deff_relperm_wat_ds = (drelperm_ds[0]-drelperm_ds[3])/rk; //derivative with respect to sw deff_relperm_wat_ds = (drelperm_ds[0]-drelperm_ds[2])/rk; //derivative with respect to sw
deff_relperm_wat_dc = dc_ads_dc*relperm[0]/(rk*rk*c_max_ads_); deff_relperm_wat_dc = dc_ads_dc*relperm[0]/(rk*rk*c_max_ads_);
} else { } else {
deff_relperm_wat_ds = -1.0; deff_relperm_wat_ds = -1.0;