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

Merge pull request #44 from qilicun/fix_eff_relperm

Browse Source 
fix derivative of effective relperm respect to concentration.
This commit is contained in:
Atgeirr Flø Rasmussen 2013-12-16 01:55:59 -08:00
commit d3453c4ded
1 changed files with 2 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
opm/polymer/PolymerProperties.cpp
View File

@ -212,7 +212,8 @@ namespace Opm
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[2])/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_); //\frac{\partial k_{rw_eff}}{\parital c} = -\frac{krw}{rk^2}\frac{(RRF-1)}{c^a_{max}}\frac{\partial c^a}{\partial c}.
deff_relperm_wat_dc = -(res_factor - 1)*dc_ads_dc*relperm[0]/(rk*rk*c_max_ads_);
} else { } else {
deff_relperm_wat_ds = -1.0; deff_relperm_wat_ds = -1.0;
deff_relperm_wat_dc = -1.0; deff_relperm_wat_dc = -1.0;