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

Merge pull request #5435 from steink/updates_in_preparation_for_5157

Browse Source 
Minor updates in preparation for New defaults #5157.
This commit is contained in:
Kai Bao 2024-06-26 00:24:29 +02:00 committed by GitHub
commit eb83fff9d5
No known key found for this signature in database
GPG Key ID: B5690EEEBB952194
1 changed files with 18 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

21
opm/simulators/wells/StandardWell_impl.hpp
View File

@ -1624,6 +1624,11 @@ namespace Opm
// prepare/modify well state and control
well_copy.prepareForPotentialCalculations(summary_state, well_state_copy, inj_controls, prod_controls);
// update connection pressures relative to updated bhp to get better estimate of connection dp
const int num_perf = ws.perf_data.size();
for (int perf = 0; perf < num_perf; ++perf) {
ws.perf_data.pressure[perf] = ws.bhp + well_copy.connections_.pressure_diff(perf);
}
// initialize rates from previous potentials
const int np = this->number_of_phases_;
bool trivial = true;
@ -1650,9 +1655,19 @@ namespace Opm
// fetch potentials (sign is updated on the outside).
well_potentials.clear();
well_potentials.resize(np, 0.0);
for (int compIdx = 0; compIdx < this->num_components_; ++compIdx) {
const EvalWell rate = well_copy.primary_variables_.getQs(compIdx);
well_potentials[this->modelCompIdxToFlowCompIdx(compIdx)] = rate.value();
for (int comp_idx = 0; comp_idx < this->num_components_; ++comp_idx) {
if (has_solvent && comp_idx == Indices::contiSolventEqIdx) continue; // we do not store the solvent in the well_potentials
const EvalWell rate = well_copy.primary_variables_.getQs(comp_idx);
well_potentials[this->modelCompIdxToFlowCompIdx(comp_idx)] = rate.value();
}
// the solvent contribution is added to the gas potentials
if constexpr (has_solvent) {
const auto& pu = this->phaseUsage();
assert(pu.phase_used[Gas]);
const int gas_pos = pu.phase_pos[Gas];
const EvalWell rate = well_copy.primary_variables_.getQs(Indices::contiSolventEqIdx);
well_potentials[gas_pos] += rate.value();
}
return converged;
}