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

Rename cell_idx to perf, since it is the perforation index actually

Browse Source
This commit is contained in:
Lisa Julia Nebel 2024-10-11 16:04:47 +02:00
parent 449847da9e
commit cde38bc23f
2 changed files with 4 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
opm/simulators/wells/MultisegmentWellAssemble.cpp
View File

@ -352,7 +352,7 @@ assembleInflowTerm(const int seg,
template<class FluidSystem, class Indices> template<class FluidSystem, class Indices>
void MultisegmentWellAssemble<FluidSystem,Indices>:: void MultisegmentWellAssemble<FluidSystem,Indices>::
assemblePerforationEq(const int seg, assemblePerforationEq(const int seg,
const int cell_idx, const int local_perf_index,
const int comp_idx, const int comp_idx,
const EvalWell& cq_s_effective, const EvalWell& cq_s_effective,
Equations& eqns1) const Equations& eqns1) const
@ -364,7 +364,7 @@ assemblePerforationEq(const int seg,
// assemble the jacobians // assemble the jacobians
for (int pv_idx = 0; pv_idx < numWellEq; ++pv_idx) { for (int pv_idx = 0; pv_idx < numWellEq; ++pv_idx) {
// also need to consider the efficiency factor when manipulating the jacobians. // also need to consider the efficiency factor when manipulating the jacobians.
eqns.C()[seg][cell_idx][pv_idx][comp_idx] -= cq_s_effective.derivative(pv_idx + Indices::numEq); // input in transformed matrix eqns.C()[seg][local_perf_index][pv_idx][comp_idx] -= cq_s_effective.derivative(pv_idx + Indices::numEq); // input in transformed matrix
// the index name for the D should be eq_idx / pv_idx // the index name for the D should be eq_idx / pv_idx
eqns.D()[seg][seg][comp_idx][pv_idx] += cq_s_effective.derivative(pv_idx + Indices::numEq); eqns.D()[seg][seg][comp_idx][pv_idx] += cq_s_effective.derivative(pv_idx + Indices::numEq);
@ -372,7 +372,7 @@ assemblePerforationEq(const int seg,
for (int pv_idx = 0; pv_idx < Indices::numEq; ++pv_idx) { for (int pv_idx = 0; pv_idx < Indices::numEq; ++pv_idx) {
// also need to consider the efficiency factor when manipulating the jacobians. // also need to consider the efficiency factor when manipulating the jacobians.
eqns.B()[seg][cell_idx][comp_idx][pv_idx] += cq_s_effective.derivative(pv_idx); eqns.B()[seg][local_perf_index][comp_idx][pv_idx] += cq_s_effective.derivative(pv_idx);
} }
} }

2
opm/simulators/wells/MultisegmentWellAssemble.hpp
View File

@ -139,7 +139,7 @@ public:
//! \brief Assemble equation for a perforation. //! \brief Assemble equation for a perforation.
void assemblePerforationEq(const int seg, void assemblePerforationEq(const int seg,
const int cell_idx, const int local_perf_index,
const int comp_idx, const int comp_idx,
const EvalWell& cq_s_effective, const EvalWell& cq_s_effective,
Equations& eqns) const; Equations& eqns) const;