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 source term (measured at reservoir condition).

Browse Source
This commit is contained in:
Xavier Raynaud 2012-09-05 13:38:30 +02:00
parent 938fc0f0d4
commit 55486d3c03
1 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
opm/polymer/TransportModelCompressiblePolymer.cpp
View File

@ -364,8 +364,8 @@ namespace Opm
B_cell0 = 1.0/tm.A0_[np*np*cell + 0];
B_cell = 1.0/tm.A_[np*np*cell + 0];
// Not clear why we multiply by B_cell source terms.
influx = src_is_inflow ? B_cell*src_flux : 0.0;
outflux = !src_is_inflow ? B_cell*src_flux : 0.0;
influx = src_is_inflow ? src_flux : 0.0;
outflux = !src_is_inflow ? src_flux : 0.0;
porevolume0 = tm.porevolume0_[cell];
porevolume = tm.porevolume_[cell];
const double vol_cell = tm.grid_.cell_volumes[cell];
@ -377,7 +377,7 @@ namespace Opm
tm.polyprops_.adsorption(c0, cmax0, ads0);
double mc;
tm.computeMc(tm.inflow_c_, mc);
influx_polymer = src_is_inflow ? B_cell*src_flux*mc : 0.0;
influx_polymer = src_is_inflow ? src_flux*mc : 0.0;
for (int i = tm.grid_.cell_facepos[cell]; i < tm.grid_.cell_facepos[cell+1]; ++i) {
int f = tm.grid_.cell_faces[i];
double flux;