OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2024-12-26 17:20:59 -06:00
Code Issues Packages Projects Releases Wiki Activity

using global index when accessing pore volume in numerical aquifer

Browse Source
This commit is contained in:
Kai Bao 2020-11-26 20:39:50 +01:00
parent 65bde5d65b
commit b6fdb7cc8e
2 changed files with 13 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
opm/simulators/aquifers/AquiferNumerical.hpp
View File

@ -51,11 +51,13 @@ public:
// Constructor
AquiferNumerical(const SingleNumericalAquifer& aquifer,
const std::unordered_map<int, int>& cartesian_to_compressed,
const Simulator& ebos_simulator)
const Simulator& ebos_simulator,
const int* global_cell)
: aquifer_(aquifer)
, ebos_simulator_(ebos_simulator)
, flux_rate_(0.)
, cumulative_flux_(0.)
, global_cell_(global_cell)
{
this->cell_to_aquifer_cell_idx_.resize(this->ebos_simulator_.gridView().size(/*codim=*/0), -1);
@ -106,6 +108,7 @@ private:
const Simulator& ebos_simulator_;
double flux_rate_; // aquifer influx rate
double cumulative_flux_; // cumulative aquifer influx
const int* global_cell_; // mapping to global index
double init_pressure_;
double pressure_; // aquifer pressure
@ -137,7 +140,7 @@ private:
const double water_pressure_reservoir = fs.pressure(waterPhaseIdx).value();
const auto& pvs = this->ebos_simulator_.vanguard().eclState().fieldProps().porv(true);
// TODO: should get this PV, how to consider the rock compressiblity
const double pv = pvs[cell_index];
const double pv = pvs[global_cell_[cell_index]];
sum_pv_pressure += pv * water_pressure_reservoir;
sum_pv += pv;
}

9
opm/simulators/aquifers/BlackoilAquiferModel_impl.hpp
View File

@ -183,6 +183,13 @@ BlackoilAquiferModel<TypeTag>::init()
}
// Get all the carter tracy aquifer properties data and put it in aquifers vector
const auto& ugrid = simulator_.vanguard().grid();
const int number_of_cells = simulator_.gridView().size(0);
const int* global_cell = Opm::UgGridHelpers::globalCell(ugrid);
cartesian_to_compressed_ = cartesianToCompressed(number_of_cells,
global_cell);
const auto& connections = aquifer.connections();
for (const auto& aq : aquifer.ct()) {
aquifers_CarterTracy.emplace_back(connections[aq.aquiferID],
@ -197,7 +204,7 @@ BlackoilAquiferModel<TypeTag>::init()
if (aquifer.hasNumericalAquifer()) {
for (const auto& elem : aquifer.numericalAquifers().aquifers()) {
this->aquifers_numerical.emplace_back(elem.second,
this->cartesian_to_compressed_, this->simulator_);
this->cartesian_to_compressed_, this->simulator_, global_cell);
}
}
}