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.

Browse Source
This commit is contained in:
Xavier Raynaud 2012-09-03 16:46:02 +02:00
parent 27dea09311
commit 612a1ffc63
1 changed files with 3 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
opm/polymer/CompressibleTpfaPolymer.cpp
View File

@ -102,6 +102,7 @@ namespace Opm
const int nc = grid_.number_of_cells;
const int np = props_.numPhases();
cell_relperm_.resize(nc*np);
cell_eff_viscosity_.resize(nc*np);
const double* cell_s = &state.saturation()[0];
props_.relperm(nc, cell_s, &allcells_[0], &cell_relperm_[0], 0);
computeWellPotentials(state);
@ -137,8 +138,8 @@ namespace Opm
props_.viscosity(nc, cell_p, cell_z, &allcells_[0], &cell_viscosity_[0], 0);
cell_phasemob_.resize(nc*np);
for (int cell = 0; cell < nc; ++cell) {
poly_props_.effectiveVisc((*c_)[cell], &cell_viscosity_[nc + 0], cell_eff_viscosity_[nc + 0]);
poly_props_.effectiveMobilities((*c_)[cell], (*cmax_)[cell], &cell_viscosity_[nc + 0], &cell_relperm_[nc + 0], &cell_phasemob_[nc + 0]);
poly_props_.effectiveVisc((*c_)[cell], &cell_viscosity_[cell + 0], cell_eff_viscosity_[cell + 0]);
poly_props_.effectiveMobilities((*c_)[cell], (*cmax_)[cell], &cell_viscosity_[cell + 0], &cell_relperm_[cell + 0], &cell_phasemob_[cell + 0]);
}
// Volume discrepancy: we have that