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 processFractions()

Browse Source 
to reduce the length of the function updatePrimaryVariablesNewton()
This commit is contained in:
Kai Bao
2018-05-11 14:58:44 +02:00
parent e2e24ec355
commit b2bace449e
2 changed files with 27 additions and 15 deletions
Download Patch File Download Diff File Expand all files Collapse all files

3
opm/autodiff/StandardWell.hpp
View File

@@ -328,6 +328,9 @@ namespace Opm
void updateWellStateFromPrimaryVariables(WellState& well_state) const; void updateWellStateFromPrimaryVariables(WellState& well_state) const;
void assembleControlEq(); void assembleControlEq();
// handle the non reasonable fractions due to numerical overshoot
void processFractions() const;
}; };
} }

39
opm/autodiff/StandardWell_impl.hpp
View File

@@ -881,9 +881,31 @@ namespace Opm
primary_variables_[SFrac] = old_primary_variables[SFrac] - dx4_limited; primary_variables_[SFrac] = old_primary_variables[SFrac] - dx4_limited;
} }
assert(FluidSystem::phaseIsActive(FluidSystem::oilPhaseIdx)); processFractions();
// TODO: the following processing of the fractions should go to a different function // updating the total rates G_t
primary_variables_[GTotal] = old_primary_variables[GTotal] - dwells[0][GTotal];
// updating the bottom hole pressure
{
const double dBHPLimit = param_.dbhp_max_rel_;
const int sign1 = dwells[0][Bhp] > 0 ? 1: -1;
const double dx1_limited = sign1 * std::min(std::abs(dwells[0][Bhp]), std::abs(old_primary_variables[Bhp]) * dBHPLimit);
// 1e5 to make sure bhp will not be below 1bar
primary_variables_[Bhp] = std::max(old_primary_variables[Bhp] - dx1_limited, 1e5);
}
}
template<typename TypeTag>
void
StandardWell<TypeTag>::
processFractions() const
{
assert(FluidSystem::phaseIsActive(FluidSystem::oilPhaseIdx));
const auto pu = phaseUsage(); const auto pu = phaseUsage();
std::vector<double> F(number_of_phases_, 0.0); std::vector<double> F(number_of_phases_, 0.0);
F[pu.phase_pos[Oil]] = 1.0; F[pu.phase_pos[Oil]] = 1.0;
@@ -952,19 +974,6 @@ namespace Opm
if(has_solvent) { if(has_solvent) {
primary_variables_[SFrac] = F_solvent; primary_variables_[SFrac] = F_solvent;
} }
// updating the total rates G_t
primary_variables_[GTotal] = old_primary_variables[GTotal] - dwells[0][GTotal];
// updating the bottom hole pressure
{
const double dBHPLimit = param_.dbhp_max_rel_;
const int sign1 = dwells[0][Bhp] > 0 ? 1: -1;
const double dx1_limited = sign1 * std::min(std::abs(dwells[0][Bhp]), std::abs(old_primary_variables[Bhp]) * dBHPLimit);
// 1e5 to make sure bhp will not be below 1bar
primary_variables_[Bhp] = std::max(old_primary_variables[Bhp] - dx1_limited, 1e5);
}
} }