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

Corrected fluidDensity() function to account for dissolved gas.

Browse Source
This commit is contained in:
Atgeirr Flø Rasmussen 2013-05-31 14:59:56 +02:00
parent 59e36f72b3
commit d49be41cae
1 changed files with 5 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
opm/autodiff/FullyImplicitBlackoilSolver.cpp
View File

@ -954,8 +954,12 @@ namespace Opm {
const std::vector<int>& cells) const const std::vector<int>& cells) const
{ {
const double* rhos = fluid_.surfaceDensity(); const double* rhos = fluid_.surfaceDensity();
ADB b = fluidReciprocFVF(phase, p, rs, cells); ADB b = fluidReciprocFVF(phase, p, rs, cells);
ADB rho = V::Constant(p.size(), 1, rhos[phase]) * b; ADB rho = V::Constant(p.size(), 1, rhos[phase]) * b;
if (phase == Oil && active_[Gas]) {
// It is correct to index into rhos with canonical phase indices.
rho += V::Constant(p.size(), 1, rhos[Gas]) * rs * b;
}
return rho; return rho;
} }