OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-25 18:55:30 -06:00
Code Issues Packages Projects Releases Wiki Activity

Fix handling of multiple pvt-regions in BlackoilPropsAdFromDeck.

Browse Source
This commit is contained in:
Atgeirr Flø Rasmussen 2015-02-23 14:05:33 +01:00
parent d53cf1cf54
commit ba86dc191c
2 changed files with 62 additions and 20 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

75
opm/autodiff/BlackoilPropsAdFromDeck.cpp
View File

@ -333,13 +333,14 @@ BlackoilPropsAdFromDeck::BlackoilPropsAdFromDeck(const BlackoilPropsAdFromDeck&
OPM_THROW(std::runtime_error, "Cannot call muWat(): water phase not present."); OPM_THROW(std::runtime_error, "Cannot call muWat(): water phase not present.");
} }
const int n = cells.size(); const int n = cells.size();
mapPvtRegions(cells);
assert(pw.size() == n); assert(pw.size() == n);
V mu(n); V mu(n);
V dmudp(n); V dmudp(n);
V dmudr(n); V dmudr(n);
const double* rs = 0; const double* rs = 0;
props_[phase_usage_.phase_pos[Water]]->mu(n, &pvtTableIdx_[0], pw.data(), T.data(), rs, props_[phase_usage_.phase_pos[Water]]->mu(n, pvt_region_.data(), pw.data(), T.data(), rs,
mu.data(), dmudp.data(), dmudr.data()); mu.data(), dmudp.data(), dmudr.data());
return mu; return mu;
} }
@ -361,12 +362,13 @@ BlackoilPropsAdFromDeck::BlackoilPropsAdFromDeck(const BlackoilPropsAdFromDeck&
OPM_THROW(std::runtime_error, "Cannot call muOil(): oil phase not present."); OPM_THROW(std::runtime_error, "Cannot call muOil(): oil phase not present.");
} }
const int n = cells.size(); const int n = cells.size();
mapPvtRegions(cells);
assert(po.size() == n); assert(po.size() == n);
V mu(n); V mu(n);
V dmudp(n); V dmudp(n);
V dmudr(n); V dmudr(n);
props_[phase_usage_.phase_pos[Oil]]->mu(n, &pvtTableIdx_[0], po.data(), T.data(), rs.data(), &cond[0], props_[phase_usage_.phase_pos[Oil]]->mu(n, pvt_region_.data(), po.data(), T.data(), rs.data(), &cond[0],
mu.data(), dmudp.data(), dmudr.data()); mu.data(), dmudp.data(), dmudr.data());
return mu; return mu;
} }
@ -384,13 +386,14 @@ BlackoilPropsAdFromDeck::BlackoilPropsAdFromDeck(const BlackoilPropsAdFromDeck&
OPM_THROW(std::runtime_error, "Cannot call muGas(): gas phase not present."); OPM_THROW(std::runtime_error, "Cannot call muGas(): gas phase not present.");
} }
const int n = cells.size(); const int n = cells.size();
mapPvtRegions(cells);
assert(pg.size() == n); assert(pg.size() == n);
V mu(n); V mu(n);
V dmudp(n); V dmudp(n);
V dmudr(n); V dmudr(n);
const double* rs = 0; const double* rs = 0;
props_[phase_usage_.phase_pos[Gas]]->mu(n, &pvtTableIdx_[0], pg.data(), T.data(), rs, props_[phase_usage_.phase_pos[Gas]]->mu(n, pvt_region_.data(), pg.data(), T.data(), rs,
mu.data(), dmudp.data(), dmudr.data()); mu.data(), dmudp.data(), dmudr.data());
return mu; return mu;
} }
@ -410,12 +413,13 @@ BlackoilPropsAdFromDeck::BlackoilPropsAdFromDeck(const BlackoilPropsAdFromDeck&
OPM_THROW(std::runtime_error, "Cannot call muGas(): gas phase not present."); OPM_THROW(std::runtime_error, "Cannot call muGas(): gas phase not present.");
} }
const int n = cells.size(); const int n = cells.size();
mapPvtRegions(cells);
assert(pg.size() == n); assert(pg.size() == n);
V mu(n); V mu(n);
V dmudp(n); V dmudp(n);
V dmudr(n); V dmudr(n);
props_[phase_usage_.phase_pos[Gas]]->mu(n, &pvtTableIdx_[0], pg.data(), T.data(), rv.data(),&cond[0], props_[phase_usage_.phase_pos[Gas]]->mu(n, pvt_region_.data(), pg.data(), T.data(), rv.data(),&cond[0],
mu.data(), dmudp.data(), dmudr.data()); mu.data(), dmudp.data(), dmudr.data());
return mu; return mu;
} }
@ -433,13 +437,14 @@ BlackoilPropsAdFromDeck::BlackoilPropsAdFromDeck(const BlackoilPropsAdFromDeck&
OPM_THROW(std::runtime_error, "Cannot call muWat(): water phase not present."); OPM_THROW(std::runtime_error, "Cannot call muWat(): water phase not present.");
} }
const int n = cells.size(); const int n = cells.size();
mapPvtRegions(cells);
assert(pw.size() == n); assert(pw.size() == n);
V mu(n); V mu(n);
V dmudp(n); V dmudp(n);
V dmudr(n); V dmudr(n);
const double* rs = 0; const double* rs = 0;
props_[phase_usage_.phase_pos[Water]]->mu(n, &pvtTableIdx_[0], pw.value().data(), T.value().data(), rs, props_[phase_usage_.phase_pos[Water]]->mu(n, pvt_region_.data(), pw.value().data(), T.value().data(), rs,
mu.data(), dmudp.data(), dmudr.data()); mu.data(), dmudp.data(), dmudr.data());
ADB::M dmudp_diag = spdiag(dmudp); ADB::M dmudp_diag = spdiag(dmudp);
const int num_blocks = pw.numBlocks(); const int num_blocks = pw.numBlocks();
@ -467,12 +472,13 @@ BlackoilPropsAdFromDeck::BlackoilPropsAdFromDeck(const BlackoilPropsAdFromDeck&
OPM_THROW(std::runtime_error, "Cannot call muOil(): oil phase not present."); OPM_THROW(std::runtime_error, "Cannot call muOil(): oil phase not present.");
} }
const int n = cells.size(); const int n = cells.size();
mapPvtRegions(cells);
assert(po.size() == n); assert(po.size() == n);
V mu(n); V mu(n);
V dmudp(n); V dmudp(n);
V dmudr(n); V dmudr(n);
props_[phase_usage_.phase_pos[Oil]]->mu(n, &pvtTableIdx_[0], po.value().data(), T.value().data(), rs.value().data(), props_[phase_usage_.phase_pos[Oil]]->mu(n, pvt_region_.data(), po.value().data(), T.value().data(), rs.value().data(),
&cond[0], mu.data(), dmudp.data(), dmudr.data()); &cond[0], mu.data(), dmudp.data(), dmudr.data());
ADB::M dmudp_diag = spdiag(dmudp); ADB::M dmudp_diag = spdiag(dmudp);
@ -501,13 +507,14 @@ BlackoilPropsAdFromDeck::BlackoilPropsAdFromDeck(const BlackoilPropsAdFromDeck&
OPM_THROW(std::runtime_error, "Cannot call muGas(): gas phase not present."); OPM_THROW(std::runtime_error, "Cannot call muGas(): gas phase not present.");
} }
const int n = cells.size(); const int n = cells.size();
mapPvtRegions(cells);
assert(pg.value().size() == n); assert(pg.value().size() == n);
V mu(n); V mu(n);
V dmudp(n); V dmudp(n);
V dmudr(n); V dmudr(n);
const double* rv = 0; const double* rv = 0;
props_[phase_usage_.phase_pos[Gas]]->mu(n, &pvtTableIdx_[0], pg.value().data(), T.value().data(), rv, props_[phase_usage_.phase_pos[Gas]]->mu(n, pvt_region_.data(), pg.value().data(), T.value().data(), rv,
mu.data(), dmudp.data(), dmudr.data()); mu.data(), dmudp.data(), dmudr.data());
ADB::M dmudp_diag = spdiag(dmudp); ADB::M dmudp_diag = spdiag(dmudp);
@ -536,12 +543,13 @@ BlackoilPropsAdFromDeck::BlackoilPropsAdFromDeck(const BlackoilPropsAdFromDeck&
OPM_THROW(std::runtime_error, "Cannot call muGas(): gas phase not present."); OPM_THROW(std::runtime_error, "Cannot call muGas(): gas phase not present.");
} }
const int n = cells.size(); const int n = cells.size();
mapPvtRegions(cells);
assert(pg.value().size() == n); assert(pg.value().size() == n);
V mu(n); V mu(n);
V dmudp(n); V dmudp(n);
V dmudr(n); V dmudr(n);
props_[phase_usage_.phase_pos[Gas]]->mu(n, &pvtTableIdx_[0], pg.value().data(), T.value().data(), rv.value().data(),&cond[0], props_[phase_usage_.phase_pos[Gas]]->mu(n, pvt_region_.data(), pg.value().data(), T.value().data(), rv.value().data(),&cond[0],
mu.data(), dmudp.data(), dmudr.data()); mu.data(), dmudp.data(), dmudr.data());
ADB::M dmudp_diag = spdiag(dmudp); ADB::M dmudp_diag = spdiag(dmudp);
@ -588,6 +596,7 @@ BlackoilPropsAdFromDeck::BlackoilPropsAdFromDeck(const BlackoilPropsAdFromDeck&
OPM_THROW(std::runtime_error, "Cannot call bWat(): water phase not present."); OPM_THROW(std::runtime_error, "Cannot call bWat(): water phase not present.");
} }
const int n = cells.size(); const int n = cells.size();
mapPvtRegions(cells);
assert(pw.size() == n); assert(pw.size() == n);
V b(n); V b(n);
@ -595,7 +604,7 @@ BlackoilPropsAdFromDeck::BlackoilPropsAdFromDeck(const BlackoilPropsAdFromDeck&
V dbdr(n); V dbdr(n);
const double* rs = 0; const double* rs = 0;
props_[phase_usage_.phase_pos[Water]]->b(n, &pvtTableIdx_[0], pw.data(), T.data(), rs, props_[phase_usage_.phase_pos[Water]]->b(n, pvt_region_.data(), pw.data(), T.data(), rs,
b.data(), dbdp.data(), dbdr.data()); b.data(), dbdp.data(), dbdr.data());
return b; return b;
@ -618,13 +627,14 @@ BlackoilPropsAdFromDeck::BlackoilPropsAdFromDeck(const BlackoilPropsAdFromDeck&
OPM_THROW(std::runtime_error, "Cannot call bOil(): oil phase not present."); OPM_THROW(std::runtime_error, "Cannot call bOil(): oil phase not present.");
} }
const int n = cells.size(); const int n = cells.size();
mapPvtRegions(cells);
assert(po.size() == n); assert(po.size() == n);
V b(n); V b(n);
V dbdp(n); V dbdp(n);
V dbdr(n); V dbdr(n);
props_[phase_usage_.phase_pos[Oil]]->b(n, &pvtTableIdx_[0], po.data(), T.data(), rs.data(), &cond[0], props_[phase_usage_.phase_pos[Oil]]->b(n, pvt_region_.data(), po.data(), T.data(), rs.data(), &cond[0],
b.data(), dbdp.data(), dbdr.data()); b.data(), dbdp.data(), dbdr.data());
return b; return b;
@ -643,6 +653,7 @@ BlackoilPropsAdFromDeck::BlackoilPropsAdFromDeck(const BlackoilPropsAdFromDeck&
OPM_THROW(std::runtime_error, "Cannot call bGas(): gas phase not present."); OPM_THROW(std::runtime_error, "Cannot call bGas(): gas phase not present.");
} }
const int n = cells.size(); const int n = cells.size();
mapPvtRegions(cells);
assert(pg.size() == n); assert(pg.size() == n);
V b(n); V b(n);
@ -650,7 +661,7 @@ BlackoilPropsAdFromDeck::BlackoilPropsAdFromDeck(const BlackoilPropsAdFromDeck&
V dbdr(n); V dbdr(n);
const double* rs = 0; const double* rs = 0;
props_[phase_usage_.phase_pos[Gas]]->b(n, &pvtTableIdx_[0], pg.data(), T.data(), rs, props_[phase_usage_.phase_pos[Gas]]->b(n, pvt_region_.data(), pg.data(), T.data(), rs,
b.data(), dbdp.data(), dbdr.data()); b.data(), dbdp.data(), dbdr.data());
return b; return b;
@ -673,13 +684,14 @@ BlackoilPropsAdFromDeck::BlackoilPropsAdFromDeck(const BlackoilPropsAdFromDeck&
OPM_THROW(std::runtime_error, "Cannot call muGas(): gas phase not present."); OPM_THROW(std::runtime_error, "Cannot call muGas(): gas phase not present.");
} }
const int n = cells.size(); const int n = cells.size();
mapPvtRegions(cells);
assert(pg.size() == n); assert(pg.size() == n);
V b(n); V b(n);
V dbdp(n); V dbdp(n);
V dbdr(n); V dbdr(n);
props_[phase_usage_.phase_pos[Gas]]->b(n, &pvtTableIdx_[0], pg.data(), T.data(), rv.data(), &cond[0], props_[phase_usage_.phase_pos[Gas]]->b(n, pvt_region_.data(), pg.data(), T.data(), rv.data(), &cond[0],
b.data(), dbdp.data(), dbdr.data()); b.data(), dbdp.data(), dbdr.data());
return b; return b;
@ -698,6 +710,7 @@ BlackoilPropsAdFromDeck::BlackoilPropsAdFromDeck(const BlackoilPropsAdFromDeck&
OPM_THROW(std::runtime_error, "Cannot call muWat(): water phase not present."); OPM_THROW(std::runtime_error, "Cannot call muWat(): water phase not present.");
} }
const int n = cells.size(); const int n = cells.size();
mapPvtRegions(cells);
assert(pw.size() == n); assert(pw.size() == n);
V b(n); V b(n);
@ -705,7 +718,7 @@ BlackoilPropsAdFromDeck::BlackoilPropsAdFromDeck(const BlackoilPropsAdFromDeck&
V dbdr(n); V dbdr(n);
const double* rs = 0; const double* rs = 0;
props_[phase_usage_.phase_pos[Water]]->b(n, &pvtTableIdx_[0], pw.value().data(), T.value().data(), rs, props_[phase_usage_.phase_pos[Water]]->b(n, pvt_region_.data(), pw.value().data(), T.value().data(), rs,
b.data(), dbdp.data(), dbdr.data()); b.data(), dbdp.data(), dbdr.data());
ADB::M dbdp_diag = spdiag(dbdp); ADB::M dbdp_diag = spdiag(dbdp);
@ -734,13 +747,14 @@ BlackoilPropsAdFromDeck::BlackoilPropsAdFromDeck(const BlackoilPropsAdFromDeck&
OPM_THROW(std::runtime_error, "Cannot call muOil(): oil phase not present."); OPM_THROW(std::runtime_error, "Cannot call muOil(): oil phase not present.");
} }
const int n = cells.size(); const int n = cells.size();
mapPvtRegions(cells);
assert(po.size() == n); assert(po.size() == n);
V b(n); V b(n);
V dbdp(n); V dbdp(n);
V dbdr(n); V dbdr(n);
props_[phase_usage_.phase_pos[Oil]]->b(n, &pvtTableIdx_[0], po.value().data(), T.value().data(), rs.value().data(), props_[phase_usage_.phase_pos[Oil]]->b(n, pvt_region_.data(), po.value().data(), T.value().data(), rs.value().data(),
&cond[0], b.data(), dbdp.data(), dbdr.data()); &cond[0], b.data(), dbdp.data(), dbdr.data());
ADB::M dbdp_diag = spdiag(dbdp); ADB::M dbdp_diag = spdiag(dbdp);
@ -769,6 +783,7 @@ BlackoilPropsAdFromDeck::BlackoilPropsAdFromDeck(const BlackoilPropsAdFromDeck&
OPM_THROW(std::runtime_error, "Cannot call muGas(): gas phase not present."); OPM_THROW(std::runtime_error, "Cannot call muGas(): gas phase not present.");
} }
const int n = cells.size(); const int n = cells.size();
mapPvtRegions(cells);
assert(pg.size() == n); assert(pg.size() == n);
V b(n); V b(n);
@ -776,7 +791,7 @@ BlackoilPropsAdFromDeck::BlackoilPropsAdFromDeck(const BlackoilPropsAdFromDeck&
V dbdr(n); V dbdr(n);
const double* rv = 0; const double* rv = 0;
props_[phase_usage_.phase_pos[Gas]]->b(n, &pvtTableIdx_[0], pg.value().data(), T.value().data(), rv, props_[phase_usage_.phase_pos[Gas]]->b(n, pvt_region_.data(), pg.value().data(), T.value().data(), rv,
b.data(), dbdp.data(), dbdr.data()); b.data(), dbdp.data(), dbdr.data());
ADB::M dbdp_diag = spdiag(dbdp); ADB::M dbdp_diag = spdiag(dbdp);
@ -805,13 +820,14 @@ BlackoilPropsAdFromDeck::BlackoilPropsAdFromDeck(const BlackoilPropsAdFromDeck&
OPM_THROW(std::runtime_error, "Cannot call muGas(): gas phase not present."); OPM_THROW(std::runtime_error, "Cannot call muGas(): gas phase not present.");
} }
const int n = cells.size(); const int n = cells.size();
mapPvtRegions(cells);
assert(pg.size() == n); assert(pg.size() == n);
V b(n); V b(n);
V dbdp(n); V dbdp(n);
V dbdr(n); V dbdr(n);
props_[phase_usage_.phase_pos[Gas]]->b(n, &pvtTableIdx_[0], pg.value().data(), T.value().data(), rv.value().data(), &cond[0], props_[phase_usage_.phase_pos[Gas]]->b(n, pvt_region_.data(), pg.value().data(), T.value().data(), rv.value().data(), &cond[0],
b.data(), dbdp.data(), dbdr.data()); b.data(), dbdp.data(), dbdr.data());
ADB::M dbdp_diag = spdiag(dbdp); ADB::M dbdp_diag = spdiag(dbdp);
@ -842,10 +858,11 @@ BlackoilPropsAdFromDeck::BlackoilPropsAdFromDeck(const BlackoilPropsAdFromDeck&
OPM_THROW(std::runtime_error, "Cannot call rsMax(): oil phase not present."); OPM_THROW(std::runtime_error, "Cannot call rsMax(): oil phase not present.");
} }
const int n = cells.size(); const int n = cells.size();
mapPvtRegions(cells);
assert(po.size() == n); assert(po.size() == n);
V rbub(n); V rbub(n);
V drbubdp(n); V drbubdp(n);
props_[phase_usage_.phase_pos[Oil]]->rsSat(n, &pvtTableIdx_[0], po.data(), rbub.data(), drbubdp.data()); props_[phase_usage_.phase_pos[Oil]]->rsSat(n, pvt_region_.data(), po.data(), rbub.data(), drbubdp.data());
return rbub; return rbub;
} }
@ -874,10 +891,11 @@ BlackoilPropsAdFromDeck::BlackoilPropsAdFromDeck(const BlackoilPropsAdFromDeck&
OPM_THROW(std::runtime_error, "Cannot call rsMax(): oil phase not present."); OPM_THROW(std::runtime_error, "Cannot call rsMax(): oil phase not present.");
} }
const int n = cells.size(); const int n = cells.size();
mapPvtRegions(cells);
assert(po.size() == n); assert(po.size() == n);
V rbub(n); V rbub(n);
V drbubdp(n); V drbubdp(n);
props_[phase_usage_.phase_pos[Oil]]->rsSat(n, &pvtTableIdx_[0], po.value().data(), rbub.data(), drbubdp.data()); props_[phase_usage_.phase_pos[Oil]]->rsSat(n, pvt_region_.data(), po.value().data(), rbub.data(), drbubdp.data());
ADB::M drbubdp_diag = spdiag(drbubdp); ADB::M drbubdp_diag = spdiag(drbubdp);
const int num_blocks = po.numBlocks(); const int num_blocks = po.numBlocks();
std::vector<ADB::M> jacs(num_blocks); std::vector<ADB::M> jacs(num_blocks);
@ -914,10 +932,11 @@ BlackoilPropsAdFromDeck::BlackoilPropsAdFromDeck(const BlackoilPropsAdFromDeck&
OPM_THROW(std::runtime_error, "Cannot call rvMax(): gas phase not present."); OPM_THROW(std::runtime_error, "Cannot call rvMax(): gas phase not present.");
} }
const int n = cells.size(); const int n = cells.size();
mapPvtRegions(cells);
assert(po.size() == n); assert(po.size() == n);
V rv(n); V rv(n);
V drvdp(n); V drvdp(n);
props_[phase_usage_.phase_pos[Gas]]->rvSat(n, &pvtTableIdx_[0], po.data(), rv.data(), drvdp.data()); props_[phase_usage_.phase_pos[Gas]]->rvSat(n, pvt_region_.data(), po.data(), rv.data(), drvdp.data());
return rv; return rv;
} }
@ -946,10 +965,11 @@ BlackoilPropsAdFromDeck::BlackoilPropsAdFromDeck(const BlackoilPropsAdFromDeck&
OPM_THROW(std::runtime_error, "Cannot call rvMax(): gas phase not present."); OPM_THROW(std::runtime_error, "Cannot call rvMax(): gas phase not present.");
} }
const int n = cells.size(); const int n = cells.size();
mapPvtRegions(cells);
assert(po.size() == n); assert(po.size() == n);
V rv(n); V rv(n);
V drvdp(n); V drvdp(n);
props_[phase_usage_.phase_pos[Gas]]->rvSat(n, &pvtTableIdx_[0], po.value().data(), rv.data(), drvdp.data()); props_[phase_usage_.phase_pos[Gas]]->rvSat(n, pvt_region_.data(), po.value().data(), rv.data(), drvdp.data());
ADB::M drvdp_diag = spdiag(drvdp); ADB::M drvdp_diag = spdiag(drvdp);
const int num_blocks = po.numBlocks(); const int num_blocks = po.numBlocks();
std::vector<ADB::M> jacs(num_blocks); std::vector<ADB::M> jacs(num_blocks);
@ -1241,5 +1261,20 @@ BlackoilPropsAdFromDeck::BlackoilPropsAdFromDeck(const BlackoilPropsAdFromDeck&
} }
} }
// Fills pvt_region_ with cellPvtRegionIdx_[cells].
void BlackoilPropsAdFromDeck::mapPvtRegions(const std::vector<int>& cells) const
{
const int n = cells.size();
pvt_region_.resize(n);
for (int ii = 0; ii < n; ++ii) {
pvt_region_[ii] = cellPvtRegionIdx_[cells[ii]];
}
}
} // namespace Opm } // namespace Opm

7
opm/autodiff/BlackoilPropsAdFromDeck.hpp
View File

@ -455,6 +455,9 @@ namespace Opm
const std::vector<int>& cells, const std::vector<int>& cells,
const double vap) const; const double vap) const;
// Fills pvt_region_ with cellPvtRegionIdx_[cells].
void mapPvtRegions(const std::vector<int>& cells) const;
RockFromDeck rock_; RockFromDeck rock_;
// This has to be a shared pointer as we must // This has to be a shared pointer as we must
// be able to make a copy of *this in the parallel case. // be able to make a copy of *this in the parallel case.
@ -467,6 +470,10 @@ namespace Opm
// The PVT region which is to be used for each cell // The PVT region which is to be used for each cell
std::vector<int> cellPvtRegionIdx_; std::vector<int> cellPvtRegionIdx_;
// Used for storing the region-per-cell array computed in calls
// to pvt functions.
mutable std::vector<int> pvt_region_;
// The PVT properties. One object per active fluid phase. // The PVT properties. One object per active fluid phase.
std::vector<std::shared_ptr<Opm::PvtInterface> > props_; std::vector<std::shared_ptr<Opm::PvtInterface> > props_;