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

Use proper transport source term to compute injected and produced fluid.

Browse Source
This commit is contained in:
Atgeirr Flø Rasmussen 2012-05-08 15:52:07 +02:00
parent 0f25f9b13e
commit 1e61783cd3
1 changed files with 2 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
examples/spu_2p.cpp
View File

@ -594,7 +594,7 @@ main(int argc, char** argv)
reorder_model.solve(&state.faceflux()[0], &porevol[0], &reorder_src[0], reorder_model.solve(&state.faceflux()[0], &porevol[0], &reorder_src[0],
stepsize, &reorder_sat[0]); stepsize, &reorder_sat[0]);
Opm::toBothSat(reorder_sat, state.saturation()); Opm::toBothSat(reorder_sat, state.saturation());
Opm::computeInjectedProduced(*props, state.saturation(), src, stepsize, injected, produced); Opm::computeInjectedProduced(*props, state.saturation(), reorder_src, stepsize, injected, produced);
if (use_segregation_split) { if (use_segregation_split) {
if (use_column_solver) { if (use_column_solver) {
if (use_gauss_seidel_gravity) { if (use_gauss_seidel_gravity) {
@ -614,7 +614,7 @@ main(int argc, char** argv)
} else { } else {
tsolver.solve(*grid->c_grid(), tsrc, stepsize, ctrl, state, linsolve, rpt); tsolver.solve(*grid->c_grid(), tsrc, stepsize, ctrl, state, linsolve, rpt);
std::cout << rpt; std::cout << rpt;
Opm::computeInjectedProduced(*props, state.saturation(), src, stepsize, injected, produced); Opm::computeInjectedProduced(*props, state.saturation(), reorder_src, stepsize, injected, produced);
} }
} }
transport_timer.stop(); transport_timer.stop();
@ -624,7 +624,6 @@ main(int argc, char** argv)
// Report volume balances. // Report volume balances.
Opm::computeSaturatedVol(porevol, state.saturation(), satvol); Opm::computeSaturatedVol(porevol, state.saturation(), satvol);
// Opm::computeInjectedProduced(*props, state.saturation(), src, simtimer.currentStepLength(), injected, produced);
tot_injected[0] += injected[0]; tot_injected[0] += injected[0];
tot_injected[1] += injected[1]; tot_injected[1] += injected[1];
tot_produced[0] += produced[0]; tot_produced[0] += produced[0];