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 bhp as a parameter of computeWellFlux in StandardWellsDense

Browse Source 
To make the function usable for computing well potentionals.
This commit is contained in:
Kai Bao
2017-01-06 14:42:37 +01:00
parent 5955cb07ba
commit 92a4d8687f
1 changed files with 4 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

8
opm/autodiff/StandardWellsDense.hpp
View File

@@ -235,7 +235,8 @@ enum WellVariablePositions {
const int cell_idx = wells().well_cells[perf]; const int cell_idx = wells().well_cells[perf];
const auto& intQuants = *(ebosSimulator.model().cachedIntensiveQuantities(cell_idx, /*timeIdx=*/0)); const auto& intQuants = *(ebosSimulator.model().cachedIntensiveQuantities(cell_idx, /*timeIdx=*/0));
std::vector<EvalWell> cq_s(np,0.0); std::vector<EvalWell> cq_s(np,0.0);
computeWellFlux(w, wells().WI[perf], intQuants, wellPerforationPressureDiffs()[perf], allow_cf, cq_s); const EvalWell bhp = getBhp(w);
computeWellFlux(w, wells().WI[perf], intQuants, bhp, wellPerforationPressureDiffs()[perf], allow_cf, cq_s);
for (int p1 = 0; p1 < np; ++p1) { for (int p1 = 0; p1 < np; ++p1) {
@@ -531,10 +532,9 @@ enum WellVariablePositions {
template<typename intensiveQuants> template<typename intensiveQuants>
void void
computeWellFlux(const int& w, const double& Tw, const intensiveQuants& intQuants, const double& cdp, const bool& allow_cf, std::vector<EvalWell>& cq_s) const computeWellFlux(const int& w, const double& Tw, const intensiveQuants& intQuants, const EvalWell& bhp, const double& cdp, const bool& allow_cf, std::vector<EvalWell>& cq_s) const
{ {
const Opm::PhaseUsage& pu = phase_usage_; const Opm::PhaseUsage& pu = fluid_->phaseUsage();
EvalWell bhp = getBhp(w);
const int np = wells().number_of_phases; const int np = wells().number_of_phases;
std::vector<EvalWell> cmix_s(np,0.0); std::vector<EvalWell> cmix_s(np,0.0);
for (int phase = 0; phase < np; ++phase) { for (int phase = 0; phase < np; ++phase) {