trivial clean up

This commit is contained in:
Tor Harald Sandve 2020-06-16 15:01:29 +02:00
parent e9d040a284
commit 43d434d42e
3 changed files with 2 additions and 6 deletions

View File

@ -206,6 +206,8 @@ namespace mswellhelpers
template <typename ValueType>
ValueType velocityHead(const double area, const ValueType& mass_rate, const ValueType& density)
{
// \Note: a factor of 2 is added to the formulation in order to match results from the
// reference simulator. This is inline with what is done for the friction loss.
return (mass_rate * mass_rate / (area * area * density));
}

View File

@ -278,9 +278,6 @@ namespace Opm
mutable int debug_cost_counter_ = 0;
// TODO: this is the old implementation, it is possible the new value does not need it anymore
std::vector<EvalWell> segment_reservoir_volume_rates_;
std::vector<std::vector<EvalWell>> segment_phase_fractions_;
std::vector<std::vector<EvalWell>> segment_phase_viscosities_;

View File

@ -50,7 +50,6 @@ namespace Opm
, segment_mass_rates_(numberOfSegments(), 0.0)
, segment_depth_diffs_(numberOfSegments(), 0.0)
, upwinding_segments_(numberOfSegments(), 0)
, segment_reservoir_volume_rates_(numberOfSegments(), 0.0)
, segment_phase_fractions_(numberOfSegments(), std::vector<EvalWell>(num_components_, 0.0)) // number of phase here?
, segment_phase_viscosities_(numberOfSegments(), std::vector<EvalWell>(num_components_, 0.0)) // number of phase here?
{
@ -1656,8 +1655,6 @@ namespace Opm
const EvalWell rate = getSegmentRateUpwinding(seg, comp_idx);
segment_mass_rates_[seg] += rate * surf_dens[comp_idx];
}
segment_reservoir_volume_rates_[seg] = segment_mass_rates_[seg] / segment_densities_[seg];
}
}