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

applying the volume effects to the well flux equation.

Browse Source 
The effects on the result are really tiny, which makes the real
verfication not easy. It did not bring anything bad for the current
working case though.
This commit is contained in:
Kai Bao 2015-10-15 17:12:02 +02:00
parent 9dea3c6b57
commit 8bde14d7ce
1 changed files with 11 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
opm/autodiff/BlackoilMultiSegmentModel_impl.hpp
View File

@ -870,6 +870,13 @@ namespace Opm {
ADB segqs = state.segqs;
// Gain of the surface volume of each component in the segment by dt
std::vector<ADB> segment_volume_change_dt(np, ADB::null());
for (int phase = 0; phase < np; ++phase) {
segment_volume_change_dt[phase] = segment_comp_surf_volume_current_[phase] -
segment_comp_surf_volume_initial_[phase];
}
for (int phase = 0; phase < np; ++phase) {
int start_segment = 0;
int start_perforation = 0;
@ -906,7 +913,10 @@ namespace Opm {
// the segment rates of this well
const ADB& segqs_well = subset(segqs, Span(nseg, 1, start_position));
segqs -= superset(cq_s_seg + well->wellOps().s2s_inlets * segqs_well, Span(nseg, 1, start_position), np * nseg_total);
const ADB segment_volume_change_dt_well = subset(segment_volume_change_dt[phase],
Span(nseg, 1, start_segment));
segqs -= superset(cq_s_seg + well->wellOps().s2s_inlets * segqs_well + segment_volume_change_dt_well,
Span(nseg, 1, start_position), np * nseg_total);
// another form for non-segment wells, keep here for future reference
// segqs -= superset(cq_s_seg, Span(1, 1, start_position), np * nseg_total);
start_segment += nseg;