48 static constexpr int numEq = BlackoilIndices::numEq;
49 using Eval = DenseAd::Evaluation<Scalar, numEq>;
52 const Simulator& simulator,
57 , connection_flux_ (connections_.size(), Eval{0})
59 this->total_face_area_ = this->initializeConnections();
65 result.cumulative_flux_ = 1.0;
70 void computeFaceAreaFraction(
const std::vector<Scalar>&
total_face_area)
override
74 this->area_fraction_ = this->totalFaceArea()
78 Scalar totalFaceArea()
const override
80 return this->total_face_area_;
88 void initFromRestart(
const data::Aquifers&
aquiferSoln)
override
95 this->cumulative_flux_ = this->area_fraction_ *
xaqPos->second.volume;
98 void initialSolutionApplied()
override
101 void beginTimeStep()
override
104 void endTimeStep()
override
106 this->flux_rate_ = this->totalFluxRate();
107 this->cumulative_flux_ +=
108 this->flux_rate_ * this->simulator_.timeStepSize();
111 data::AquiferData aquiferData()
const override
113 data::AquiferData data;
115 data.aquiferID = this->aquifer_data_.id;
119 data.fluxRate = this->totalFluxRate();
121 data.volume = this->cumulative_flux_;
124 data.initPressure = 0.0;
129 void addToSource(RateVector& rates,
133 const int idx = this->cellToConnectionIdx_[
cellIdx];
138 const auto& model = this->simulator_.model();
140 const auto fw = this->aquifer_data_.flux;
142 this->connection_flux_[
idx] =
fw * this->connections_[
idx].effective_facearea;
144 rates[BlackoilIndices::conti0EqIdx + compIdx_()]
145 += this->connection_flux_[
idx] / model.dofTotalVolume(
cellIdx);
148 template<
class Serializer>
156 return this->cumulative_flux_ == rhs.cumulative_flux_;
160 const std::vector<Aquancon::AquancCell>& connections_;
163 std::vector<Eval> connection_flux_{};
164 std::vector<int> cellToConnectionIdx_{};
166 Scalar cumulative_flux_{};
167 Scalar total_face_area_{0.0};
168 Scalar area_fraction_{1.0};
170 Scalar initializeConnections()
174 this->cellToConnectionIdx_
175 .resize(this->simulator_.gridView().size(0), -1);
177 for (std::size_t
idx = 0;
idx < this->connections_.size(); ++
idx) {
179 const int cell_index = this->simulator_.vanguard()
182 if (cell_index < 0) {
186 this->cellToConnectionIdx_[cell_index] =
idx;
211 if (this->co2store_or_h2store_())
212 return FluidSystem::oilCompIdx;
214 return FluidSystem::waterCompIdx;
217 Scalar totalFluxRate()
const
219 return std::accumulate(this->connection_flux_.begin(),
220 this->connection_flux_.end(), 0.0,
221 [](
const Scalar rate,
const auto&
q)
223 return rate + getValue(q);
constexpr auto getPropValue()
get the value data member of a property
Definition propertysystem.hh:242
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