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

Merge pull request #745 from goncalvesmachadoc/salt-density

Browse Source 
add Salt Density
This commit is contained in:
Bård Skaflestad 2022-11-24 16:32:16 +01:00 committed by GitHub
commit 6cd200a254
2 changed files with 22 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

22
opm/models/blackoil/blackoilbrinemodules.hh
View File

@ -136,10 +136,12 @@ public:
const TableContainer& saltsolTables = tableManager.getSaltsolTables();
if (!saltsolTables.empty()) {
params_.saltsolTable_.resize(numPvtRegions);
params_.saltdenTable_.resize(numPvtRegions);
assert(numPvtRegions == saltsolTables.size());
for (unsigned pvtRegionIdx = 0; pvtRegionIdx < numPvtRegions; ++ pvtRegionIdx) {
const SaltsolTable& saltsolTable = saltsolTables.getTable<SaltsolTable>(pvtRegionIdx );
params_.saltsolTable_[pvtRegionIdx] = saltsolTable.getSaltsolColumn().front();
params_.saltdenTable_[pvtRegionIdx] = saltsolTable.getSaltdenColumn().front();
}
}
}
@ -231,7 +233,7 @@ public:
* Toolbox::template decay<LhsEval>(fs.saltConcentration());
if (enableSaltPrecipitation){
double saltDensity = 2170; // Solid salt density kg/m3
double saltDensity = intQuants.saltDensity(); // Solid salt density kg/m3
const LhsEval solidSalt =
Toolbox::template decay<LhsEval>(intQuants.porosity())
/ (1.0 - Toolbox::template decay<LhsEval>(intQuants.saltSaturation()) + 1.e-8)
@ -342,6 +344,14 @@ public:
return params_.saltsolTable_[pvtnumRegionIdx];
}
static const Scalar saltdenTable(const ElementContext& elemCtx,
unsigned scvIdx,
unsigned timeIdx)
{
unsigned pvtnumRegionIdx = elemCtx.problem().pvtRegionIndex(elemCtx, scvIdx, timeIdx);
return params_.saltdenTable_[pvtnumRegionIdx];
}
static bool hasBDensityTables()
{
return !params_.bdensityTable_.empty();
@ -415,6 +425,9 @@ public:
const auto& saltsolTable = BrineModule::saltsolTable(elemCtx, dofIdx, timeIdx);
saltSolubility_ = saltsolTable;
const auto& saltdenTable = BrineModule::saltdenTable(elemCtx, dofIdx, timeIdx);
saltDensity_ = saltdenTable;
if (priVars.primaryVarsMeaningBrine() == PrimaryVariables::Sp) {
saltSaturation_ = priVars.makeEvaluation(saltConcentrationIdx, timeIdx);
fs.setSaltConcentration(saltSolubility_);
@ -462,6 +475,9 @@ public:
Scalar saltSolubility() const
{ return saltSolubility_; }
Scalar saltDensity() const
{ return saltDensity_; }
const Evaluation& permFactor() const
{ return permFactor_; }
@ -474,6 +490,7 @@ protected:
Evaluation saltSaturation_;
Evaluation permFactor_;
Scalar saltSolubility_;
Scalar saltDensity_;
};
@ -507,6 +524,9 @@ public:
const Scalar saltSolubility() const
{ throw std::logic_error("saltSolubility() called but salt precipitation is disabled"); }
const Scalar saltDensity() const
{ throw std::logic_error("saltDensity() called but salt precipitation is disabled"); }
const Evaluation& permFactor() const
{ throw std::logic_error("permFactor() called but salt precipitation is disabled"); }

1
opm/models/blackoil/blackoilbrineparams.hh
View File

@ -42,6 +42,7 @@ struct BlackOilBrineParams {
std::vector<TabulatedFunction> bdensityTable_;
std::vector<TabulatedFunction> permfactTable_;
std::vector<Scalar> saltsolTable_;
std::vector<Scalar> saltdenTable_;
std::vector<Scalar> referencePressure_;
};