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

calculating the surface volume of component in segment

Browse Source 
For convinience, we calculate and store the surface volume by dt
instead.
This commit is contained in:
Kai Bao 2015-10-15 15:18:22 +02:00
parent e82828527c
commit 95eb44e530
2 changed files with 26 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
opm/autodiff/BlackoilMultiSegmentModel.hpp
View File

@ -199,6 +199,12 @@ namespace Opm {
// for the top segment, they should always be zero for the moment.
ADB well_segment_pressures_delta_;
// the surface volume of components in the segments
// initial one at the beginning of the time step
std::vector<V> segment_comp_surf_volume_initial_;
// the current one at the current iteration.
std::vector<ADB> segment_comp_surf_volume_current_;
const std::vector<WellMultiSegmentConstPtr> wells_multisegment_;
std::vector<int> top_well_segments_;
@ -270,8 +276,10 @@ namespace Opm {
std::vector<ADB>& vars,
SolutionState& state) const;
// Calculate the density of the mixture in the segments
// And the surface volume of the components in the segments by dt
void
computeSegmentDensities(const SolutionState& state);
computeSegmentDensitiesAndCompVolumeDt(const SolutionState& state);
void
computeSegmentPressuresDelta(const SolutionState& state);

19
opm/autodiff/BlackoilMultiSegmentModel_impl.hpp
View File

@ -84,6 +84,8 @@ namespace Opm {
, well_perforation_cell_densities_adb_(ADB::null())
, well_segment_densities_(ADB::null())
, well_segment_pressures_delta_(ADB::null())
, segment_comp_surf_volume_initial_(fluid.numPhases())
, segment_comp_surf_volume_current_(fluid.numPhases(), ADB::null())
, wells_multisegment_(wells_multisegment)
{
// Modify the wops_.well_cell member, since the
@ -507,6 +509,12 @@ namespace Opm {
// Compute initial accumulation contributions
// and well connection pressures.
asImpl().computeAccum(state0, 0);
asImpl().computeSegmentDensitiesAndCompVolumeDt(state0);
const int np = numPhases();
assert(np == int(segment_comp_surf_volume_initial_.size()));
for (int phase = 0; phase < np; ++phase) {
segment_comp_surf_volume_initial_[phase] = segment_comp_surf_volume_current_[phase].value();
}
asImpl().computeWellConnectionPressures(state0, well_state);
}
@ -528,7 +536,7 @@ namespace Opm {
return;
}
asImpl().computeSegmentDensities(state);
asImpl().computeSegmentDensitiesAndCompVolumeDt(state);
asImpl().computeSegmentPressuresDelta(state);
std::vector<ADB> mob_perfcells;
@ -1649,7 +1657,7 @@ namespace Opm {
template <class Grid>
void
BlackoilMultiSegmentModel<Grid>::computeSegmentDensities(const SolutionState& state)
BlackoilMultiSegmentModel<Grid>::computeSegmentDensitiesAndCompVolumeDt(const SolutionState& state)
{
const int nw = wellsMultiSegment().size();
const int nseg_total = state.segp.size();
@ -1829,6 +1837,13 @@ namespace Opm {
}
well_segment_densities_ = dens / volrat;
// Calculating the surface volume of each component in the segment
assert(np == int(segment_comp_surf_volume_current_.size()));
const ADB segment_surface_volume = segvdt_ / volrat;
for (int phase = 0; phase < np; ++phase) {
segment_comp_surf_volume_current_[phase] = segment_surface_volume * mix[phase];
}
#if 0
std::cout << " output the well_segment_densities_ " << std::endl;
std::cout << well_segment_densities_.value() << std::endl;