57 using Simulator = GetPropType<TypeTag, Properties::Simulator>;
58 using ElementContext = GetPropType<TypeTag, Properties::ElementContext>;
59 using FluidSystem = GetPropType<TypeTag, Properties::FluidSystem>;
60 using BlackoilIndices = GetPropType<TypeTag, Properties::Indices>;
61 using RateVector = GetPropType<TypeTag, Properties::RateVector>;
62 using IntensiveQuantities = GetPropType<TypeTag, Properties::IntensiveQuantities>;
63 using ElementMapper = GetPropType<TypeTag, Properties::ElementMapper>;
65 enum { enableTemperature = getPropValue<TypeTag, Properties::EnableTemperature>() };
66 enum { enableEnergy = getPropValue<TypeTag, Properties::EnableEnergy>() };
67 enum { enableBrine = getPropValue<TypeTag, Properties::EnableBrine>() };
68 enum { enableEvaporation = getPropValue<TypeTag, Properties::EnableEvaporation>() };
69 enum { has_disgas_in_water = getPropValue<TypeTag, Properties::EnableDisgasInWater>() };
71 enum { enableSaltPrecipitation = getPropValue<TypeTag, Properties::EnableSaltPrecipitation>() };
73 static constexpr int numEq = BlackoilIndices::numEq;
74 using Scalar = double;
76 using Eval = DenseAd::Evaluation<double, numEq>;
78 using FluidState = BlackOilFluidState<Eval,
82 BlackoilIndices::gasEnabled,
85 enableSaltPrecipitation,
87 BlackoilIndices::numPhases>;
91 const std::vector<Aquancon::AquancCell>& connections,
92 const Simulator& ebosSimulator)
103 void initFromRestart(
const data::Aquifers&
aquiferSoln)
override
109 this->assignRestartData(
xaqPos->second);
111 this->W_flux_ =
xaqPos->second.volume;
112 this->pa0_ =
xaqPos->second.initPressure;
113 this->solution_set_from_restart_ =
true;
116 void initialSolutionApplied()
override
121 void beginTimeStep()
override
123 ElementContext
elemCtx(this->ebos_simulator_);
124 OPM_BEGIN_PARALLEL_TRY_CATCH();
126 for (
const auto&
elem :
elements(this->ebos_simulator_.gridView())) {
130 const int idx = cellToConnectionIdx_[
cellIdx];
134 elemCtx.updateIntensiveQuantities(0);
135 const auto&
iq =
elemCtx.intensiveQuantities(0, 0);
139 OPM_END_PARALLEL_TRY_CATCH(
"AquiferAnalytical::beginTimeStep() failed: ",
140 this->ebos_simulator_.vanguard().grid().comm());
143 void addToSource(RateVector& rates,
145 const unsigned timeIdx)
override
147 const auto& model = this->ebos_simulator_.model();
149 const int idx = this->cellToConnectionIdx_[
cellIdx];
156 this->updateCellPressure(this->pressure_current_,
idx,
intQuants);
157 this->calculateInflowRate(
idx, this->ebos_simulator_);
159 rates[BlackoilIndices::conti0EqIdx + compIdx_()]
160 += this->Qai_[
idx] / model.dofTotalVolume(
cellIdx);
162 if constexpr (enableEnergy) {
164 if (this->Ta0_.has_value() && this->Qai_[
idx] > 0)
166 fs.setTemperature(this->Ta0_.value());
167 typedef typename std::decay<
decltype(
fs)>::type::Scalar
FsScalar;
169 const unsigned pvtRegionIdx =
intQuants.pvtRegionIndex();
171 paramCache.setMaxOilSat(this->ebos_simulator_.problem().maxOilSaturation(
cellIdx));
173 const auto&
h = FluidSystem::enthalpy(
fs,
paramCache, this->phaseIdx_());
174 fs.setEnthalpy(this->phaseIdx_(),
h);
176 rates[BlackoilIndices::contiEnergyEqIdx]
177 += this->Qai_[
idx] *
fs.enthalpy(this->phaseIdx_()) * FluidSystem::referenceDensity( this->phaseIdx_(),
intQuants.pvtRegionIndex()) / model.dofTotalVolume(
cellIdx);
182 std::size_t size()
const
184 return this->connections_.size();
187 template<
class Serializer>
199 return this->pressure_previous_ == rhs.pressure_previous_ &&
200 this->pressure_current_ == rhs.pressure_current_ &&
201 this->Qai_ == rhs.Qai_ &&
202 this->rhow_ == rhs.rhow_ &&
203 this->W_flux_ == rhs.W_flux_;
207 virtual void assignRestartData(
const data::AquiferData&
xaq) = 0;
208 virtual void calculateInflowRate(
int idx,
const Simulator&
simulator) = 0;
209 virtual void calculateAquiferCondition() = 0;
210 virtual void calculateAquiferConstants() = 0;
211 virtual Scalar aquiferDepth()
const = 0;
213 Scalar gravity_()
const
215 return this->ebos_simulator_.problem().gravity()[2];
220 if (this->co2store_())
221 return FluidSystem::oilCompIdx;
223 return FluidSystem::waterCompIdx;
226 void initQuantities()
229 if (!this->solution_set_from_restart_) {
235 initializeConnections();
236 calculateAquiferCondition();
237 calculateAquiferConstants();
239 pressure_previous_.resize(this->connections_.size(), Scalar{0});
240 pressure_current_.resize(this->connections_.size(), Scalar{0});
241 Qai_.resize(this->connections_.size(), Scalar{0});
260 void initializeConnections()
262 this->cell_depth_.resize(this->size(), this->aquiferDepth());
263 this->alphai_.resize(this->size(), 1.0);
264 this->faceArea_connected_.resize(this->size(), Scalar{0});
273 this->cellToConnectionIdx_.resize(this->ebos_simulator_.gridView().size(0), -1);
274 const auto&
gridView = this->ebos_simulator_.vanguard().gridView();
275 for (std::size_t
idx = 0;
idx < this->size(); ++
idx) {
277 const int cell_index = this->ebos_simulator_.vanguard().compressedIndex(
global_index);
280 std::advance(
elemIt, cell_index);
288 this->cellToConnectionIdx_[cell_index] =
idx;
289 this->cell_depth_.at(
idx) = this->ebos_simulator_.vanguard().cellCenterDepth(cell_index);
295 int idx = this->cellToConnectionIdx_[cell_index];
328 "Internal error in initialization of aquifer.");
333 this->faceArea_connected_[
idx] = this->connections_[
idx].influx_coeff;
340 const auto& comm = this->ebos_simulator_.vanguard().grid().comm();
342 const double eps_sqrt = std::sqrt(std::numeric_limits<double>::epsilon());
343 for (std::size_t
idx = 0;
idx < this->size(); ++
idx) {
350 if (this->solution_set_from_restart_) {
364 ? std::accumulate(this->alphai_.begin(),
369 const auto tot_area = this->ebos_simulator_.vanguard()
376 Scalar calculateReservoirEquilibrium()
382 ElementContext
elemCtx(this->ebos_simulator_);
383 const auto&
gridView = this->ebos_simulator_.gridView();
388 const auto idx = this->cellToConnectionIdx_[
cellIdx];
392 elemCtx.updatePrimaryIntensiveQuantities(0);
393 const auto&
iq0 =
elemCtx.intensiveQuantities(0, 0);
394 const auto&
fs =
iq0.fluidState();
400 this->gravity_() * (this->cell_depth_[
idx] - this->aquiferDepth());
407 const auto& comm = this->ebos_simulator_.vanguard().grid().comm();
410 vals[0] = std::accumulate(this->alphai_.begin(),
this->alphai_.end(), Scalar{0});
418 const std::vector<Aquancon::AquancCell> connections_;
421 std::vector<Scalar> faceArea_connected_;
422 std::vector<int> cellToConnectionIdx_;
425 std::vector<Scalar> cell_depth_;
426 std::vector<Scalar> pressure_previous_;
427 std::vector<Eval> pressure_current_;
428 std::vector<Eval> Qai_;
429 std::vector<Scalar> alphai_;
433 std::optional<Scalar> Ta0_{};
438 bool solution_set_from_restart_ {
false};