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

a few bug fixing for MultisegmentWell

Browse Source
This commit is contained in:
Kai Bao
2017-09-15 16:28:10 +02:00
parent 5ffd2bf026
commit f784cf97eb
2 changed files with 11 additions and 8 deletions
Download Patch File Download Diff File Expand all files Collapse all files

2
opm/autodiff/MultisegmentWell.hpp
View File

@@ -332,7 +332,7 @@ namespace Opm
// handling the overshooting and undershooting of the fractions // handling the overshooting and undershooting of the fractions
void processFractions(const int seg) const; void processFractions(const int seg) const;
void updateWellStateFromPrimaryVariabls(WellState& well_state) const; void updateWellStateFromPrimaryVariables(WellState& well_state) const;
}; };
// obtain y = D^-1 * x // obtain y = D^-1 * x

17
opm/autodiff/MultisegmentWell_impl.hpp
View File

@@ -787,24 +787,27 @@ namespace Opm
primary_variables_[seg][GFrac] = old_primary_variables[seg][GFrac] - dx_limited; primary_variables_[seg][GFrac] = old_primary_variables[seg][GFrac] - dx_limited;
} }
// handling the overshooting or undershooting of the fractions
processFractions(seg);
// update the segment pressure // update the segment pressure
{ {
const int sign = dwells[seg][SPres] > 0.? 1 : -1; const int sign = dwells[seg][SPres] > 0.? 1 : -1;
const double dx_limited = sign * std::min(std::abs(dwells[seg][SPres]), dBHPLimit); const double current_pressure = old_primary_variables[seg][SPres];
const double dx_limited = sign * std::min(std::abs(dwells[seg][SPres]), dBHPLimit * current_pressure);
primary_variables_[seg][SPres] = old_primary_variables[seg][SPres] - dx_limited; primary_variables_[seg][SPres] = old_primary_variables[seg][SPres] - dx_limited;
} }
// update the total rate // TODO: should we have a limitation of the total rate change
// update the total rate // TODO: should we have a limitation of the total rate change?
{ {
primary_variables_[seg][SPres] = old_primary_variables[seg][SPres] - dwells[seg][SPres]; primary_variables_[seg][GTotal] = old_primary_variables[seg][GTotal] - dwells[seg][GTotal];
} }
// TODO: not handling solvent related for now // TODO: not handling solvent related for now
// handling the overshooting or undershooting of the fraction first
processFractions(seg);
} }
updateWellStateFromPrimaryVariabls(well_state); updateWellStateFromPrimaryVariables(well_state);
} }
@@ -1583,7 +1586,7 @@ namespace Opm
template<typename TypeTag> template<typename TypeTag>
void void
MultisegmentWell<TypeTag>:: MultisegmentWell<TypeTag>::
updateWellStateFromPrimaryVariabls(WellState& well_state) const updateWellStateFromPrimaryVariables(WellState& well_state) const
{ {
for (int seg = 0; seg < numberOfSegments(); ++seg) { for (int seg = 0; seg < numberOfSegments(); ++seg) {
std::vector<double> fractions(number_of_phases_, 0.0); std::vector<double> fractions(number_of_phases_, 0.0);