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

adding function computeInitialComposition to MultisegmentWell

Browse Source
This commit is contained in:
Kai Bao 2017-09-05 15:34:31 +02:00
parent eb119f245d
commit 1e32b40837
2 changed files with 8 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
opm/autodiff/MultisegmentWell.hpp
View File

@ -236,7 +236,7 @@ namespace Opm
std::vector<double> perforation_cell_pressure_diffs_;
// the intial component compistion of segments
std::vector<double> segment_comp_initial_;
std::vector<std::vector<double> > segment_comp_initial_;
void initMatrixAndVectors(const int num_cells) const;

10
opm/autodiff/MultisegmentWell_impl.hpp
View File

@ -32,7 +32,7 @@ namespace Opm
, segment_perforations_(numberOfSegments())
, segment_inlets_(numberOfSegments())
, perforation_cell_pressure_diffs_(number_of_perforations_, 0.0)
, segment_comp_initial_(numberOfSegments(), 0.0)
, segment_comp_initial_(numberOfSegments(), std::vector<double>(numWellEq, 0.0))
{
// TODO: to see what information we need to process here later.
// const auto& completion_set = well->getCompletions(time_step);
@ -511,8 +511,11 @@ namespace Opm
MultisegmentWell<TypeTag>::
computeInitialComposition()
{
// TODO this, we need to implement the wellSurfaceVolumeFraction() function
// should we provide member variables to store all these values to avoid calculate again and again?
for (int seg = 0; seg < numberOfSegments(); ++seg) {
for (int eq_idx = 0; eq_idx < numWellEq; ++eq_idx) {
segment_comp_initial_[seg][eq_idx] = surfaceVolumeFraction(seg, eq_idx).value();
}
}
}
@ -526,6 +529,7 @@ namespace Opm
const WellState& /* well_state */)
{
computePerfCellPressDiffs(ebosSimulator);
computeInitialComposition();
}