OilfieldToolsNet/opm-core
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Rewriting the d_mu_d_p for deadGas.

Browse Source
This commit is contained in:
Kai Bao 2014-09-14 16:24:37 +02:00
parent d9eb3ea30f
commit d9acd0f3b8
1 changed files with 4 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
opm/core/props/pvt/PvtDead.cpp
View File

@ -156,7 +156,10 @@ namespace Opm
double tempInvBV = inverseBV_[regionIdx](p[i]); double tempInvBV = inverseBV_[regionIdx](p[i]);
// output_mu[i] = viscosity_[regionIdx](p[i]); // output_mu[i] = viscosity_[regionIdx](p[i]);
output_mu[i] = tempInvB / tempInvBV; output_mu[i] = tempInvB / tempInvBV;
output_dmudp[i] = viscosity_[regionIdx].derivative(p[i]); // output_dmudp[i] = viscosity_[regionIdx].derivative(p[i]);
// output_dmudp[i] = tempInvB / (tempInvBV * tempInvBV) * inverseBV_[regionIdx].derivative(p[i]);
output_dmudp[i] = (tempInvBV * b_[regionIdx].derivative(p[i])
- tempInvB * inverseBV_[regionIdx].derivative(p[i])) / (tempInvBV * tempInvBV);
} }
std::fill(output_dmudr, output_dmudr + n, 0.0); std::fill(output_dmudr, output_dmudr + n, 0.0);