70 GetPropType<TypeTag, Properties::Indices>>
85 using GLiftOptWells =
typename BlackoilWellModel<TypeTag>::GLiftOptWells;
86 using GLiftProdWells =
typename BlackoilWellModel<TypeTag>::GLiftProdWells;
87 using GLiftWellStateMap =
88 typename BlackoilWellModel<TypeTag>::GLiftWellStateMap;
89 using GLiftSyncGroups =
typename GasLiftSingleWellGeneric<Scalar>::GLiftSyncGroups;
91 using VectorBlockType = Dune::FieldVector<Scalar, Indices::numEq>;
92 using MatrixBlockType = Dune::FieldMatrix<Scalar, Indices::numEq, Indices::numEq>;
93 using Eval =
typename Base::Eval;
94 using BVector = Dune::BlockVector<VectorBlockType>;
95 using PressureMatrix = Dune::BCRSMatrix<Opm::MatrixBlock<Scalar, 1, 1>>;
97 using RateConverterType =
104 using ModelParameters =
typename Base::ModelParameters;
125 Indices::compositionSwitchIdx >= 0,
130 Indices::numPhases >;
135 const ModelParameters& param,
137 const int pvtRegionIdx,
139 const int num_phases,
149 const std::vector<Scalar>&
B_avg,
152 virtual void initPrimaryVariablesEvaluation() = 0;
156 const std::vector<Scalar>&
B_avg,
160 virtual void solveEqAndUpdateWellState(
const Simulator& simulator,
164 void assembleWellEq(
const Simulator& simulator,
170 void assembleWellEqWithoutIteration(
const Simulator& simulator,
177 void prepareWellBeforeAssembling(
const Simulator& simulator,
184 virtual void computeWellRatesWithBhp(
const Simulator&
ebosSimulator,
189 virtual std::optional<Scalar>
203 virtual void apply(
const BVector& x, BVector&
Ax)
const = 0;
206 virtual void apply(BVector&
r)
const = 0;
209 virtual void computeWellPotentials(
const Simulator& simulator,
211 std::vector<Scalar>& well_potentials,
214 virtual void updateWellStateWithTarget(
const Simulator& simulator,
219 virtual void computeWellRatesWithBhpIterations(
const Simulator& simulator,
224 bool wellUnderZeroRateTarget(
const Simulator& simulator,
228 bool wellUnderZeroGroupRateTarget(
const Simulator& simulator,
233 bool stoppedOrZeroRateTarget(
const Simulator& simulator,
237 bool updateWellStateWithTHPTargetProd(
const Simulator& simulator,
241 enum class IndividualOrGroup { Individual, Group, Both };
242 bool updateWellControl(
const Simulator& simulator,
243 const IndividualOrGroup
iog,
248 bool updateWellControlAndStatusLocalIteration(
const Simulator& simulator,
253 const Scalar WQTotal,
258 virtual void updatePrimaryVariables(
const Simulator& simulator,
262 virtual void calculateExplicitQuantities(
const Simulator& simulator,
266 virtual void updateProductivityIndex(
const Simulator& simulator,
281 virtual void addWellContributions(SparseMatrixAdapter&)
const = 0;
283 virtual void addWellPressureEquations(PressureMatrix& mat,
285 const int pressureVarIndex,
289 void addCellRates(RateVector& rates,
int cellIdx)
const;
291 Scalar volumetricSurfaceRateForConnection(
int cellIdx,
int phaseIdx)
const;
295 void wellTesting(
const Simulator& simulator,
302 void checkWellOperability(
const Simulator& simulator,
312 void gliftBeginTimeStepWellTestUpdateALQ(
const Simulator& simulator,
318 void updateWellOperability(
const Simulator& simulator,
322 bool updateWellOperabilityFromWellEq(
const Simulator& simulator,
327 virtual void updateWaterThroughput(
const double dt,
332 virtual std::vector<Scalar>
343 void solveWellEquation(
const Simulator& simulator,
348 const std::vector<RateVector>& connectionRates()
const
350 return connectionRates_;
353 virtual std::vector<Scalar> getPrimaryVars()
const
358 virtual int setPrimaryVars(
typename std::vector<Scalar>::const_iterator)
363 std::vector<Scalar> wellIndex(
const int perf,
368 void updateConnectionDFactor(
const Simulator& simulator,
371 void updateConnectionTransmissibilityFactor(
const Simulator& simulator,
374 virtual bool iterateWellEqWithSwitching(
const Simulator& simulator,
385 std::vector<RateVector> connectionRates_;
386 std::vector<Scalar> B_avg_;
387 bool changed_to_stopped_this_step_ =
false;
388 bool thp_update_iterations =
false;
390 Scalar wpolymer()
const;
391 Scalar wfoam()
const;
392 Scalar wsalt()
const;
393 Scalar wmicrobes()
const;
394 Scalar woxygen()
const;
395 Scalar wurea()
const;
397 virtual Scalar getRefDensity()
const = 0;
400 const std::vector<Scalar>& compFrac()
const;
408 const Simulator& simulator,
412 virtual void checkOperabilityUnderTHPLimit(
const Simulator& simulator,
416 virtual void updateIPR(
const Simulator& simulator,
419 virtual void assembleWellEqWithoutIteration(
const Simulator& simulator,
428 virtual bool iterateWellEqWithControl(
const Simulator& simulator,
436 virtual void updateIPRImplicit(
const Simulator& simulator,
440 bool iterateWellEquations(
const Simulator& simulator,
446 bool solveWellWithTHPConstraint(
const Simulator& simulator,
454 std::optional<Scalar>
461 bool solveWellWithBhp(
const Simulator& simulator,
467 bool solveWellWithZeroRate(
const Simulator& simulator,
472 bool solveWellForTesting(
const Simulator& simulator,
477 Eval getPerfCellPressure(
const FluidState&
fs)
const;
480 template<
class Value,
class Callback>
481 void getMobility(
const Simulator& simulator,
483 std::vector<Value>&
mob,
487 void computeConnLevelProdInd(
const FluidState&
fs,
488 const std::function<Scalar(
const Scalar)>&
connPICalc,
489 const std::vector<Scalar>& mobility,
492 void computeConnLevelInjInd(
const FluidState&
fs,
494 const std::function<Scalar(
const Scalar)>&
connIICalc,
495 const std::vector<Scalar>& mobility,
499 Scalar computeConnectionDFactor(
const int perf,
virtual void recoverWellSolutionAndUpdateWellState(const Simulator &simulator, const BVector &x, WellState< Scalar > &well_state, DeferredLogger &deferred_logger)=0
using the solution x to recover the solution xw for wells and applying xw to update Well State
void updateWellStateRates(const Simulator &simulator, WellState< Scalar > &well_state, DeferredLogger &deferred_logger) const
Modify the well_state's rates if there is only one nonzero rate.
Definition WellInterface_impl.hpp:1572
typename Properties::Detail::GetPropImpl< TypeTag, Property >::type::type GetPropType
get the type alias defined in the property (equivalent to old macro GET_PROP_TYPE(....
Definition propertysystem.hh:235