52 static constexpr int numEq = BlackoilIndices::numEq;
53 using Eval = DenseAd::Evaluation<
double, numEq>;
61 , connection_flux_ (connections_.size(), Eval{0})
63 this->total_face_area_ = this->initializeConnections();
69 result.cumulative_flux_ = 1.0;
76 void computeFaceAreaFraction(
const std::vector<double>&
total_face_area)
override
80 this->area_fraction_ = this->totalFaceArea()
84 double totalFaceArea()
const override
86 return this->total_face_area_;
94 void initFromRestart(
const data::Aquifers&
aquiferSoln)
override
101 this->cumulative_flux_ = this->area_fraction_ *
xaqPos->second.volume;
104 void initialSolutionApplied()
override
107 void beginTimeStep()
override
110 void endTimeStep()
override
112 this->flux_rate_ = this->totalFluxRate();
113 this->cumulative_flux_ +=
114 this->flux_rate_ * this->ebos_simulator_.timeStepSize();
117 data::AquiferData aquiferData()
const override
119 data::AquiferData data;
121 data.aquiferID = this->aquifer_data_.id;
125 data.fluxRate = this->totalFluxRate();
127 data.volume = this->cumulative_flux_;
130 data.initPressure = 0.0;
135 void addToSource(RateVector& rates,
139 const int idx = this->cellToConnectionIdx_[
cellIdx];
144 const auto& model = this->ebos_simulator_.model();
146 const auto fw = this->aquifer_data_.flux;
148 this->connection_flux_[
idx] =
fw * this->connections_[
idx].effective_facearea;
150 rates[BlackoilIndices::conti0EqIdx + compIdx_()]
151 += this->connection_flux_[
idx] / model.dofTotalVolume(
cellIdx);
154 template<
class Serializer>
162 return this->cumulative_flux_ == rhs.cumulative_flux_;
166 const std::vector<Aquancon::AquancCell>& connections_;
169 std::vector<Eval> connection_flux_{};
170 std::vector<int> cellToConnectionIdx_{};
172 double cumulative_flux_{};
173 double total_face_area_{0.0};
174 double area_fraction_{1.0};
176 double initializeConnections()
180 this->cellToConnectionIdx_
181 .resize(this->ebos_simulator_.gridView().size(0), -1);
183 for (std::size_t
idx = 0;
idx < this->connections_.size(); ++
idx) {
185 const int cell_index = this->ebos_simulator_.vanguard()
188 if (cell_index < 0) {
192 this->cellToConnectionIdx_[cell_index] =
idx;
217 if (this->co2store_or_h2store_())
218 return FluidSystem::oilCompIdx;
220 return FluidSystem::waterCompIdx;
223 double totalFluxRate()
const
225 return std::accumulate(this->connection_flux_.begin(),
226 this->connection_flux_.end(), 0.0,
227 [](
const double rate,
const auto&
q)
229 return rate + getValue(q);