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

Set eq weights externally in the blackoil model

Browse Source
This commit is contained in:
Tor Harald Sandve 2022-06-30 11:25:44 +02:00
parent 6e0aef4860
commit 6405b2b6e0
1 changed files with 9 additions and 35 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

42
opm/models/blackoil/blackoilmodel.hh
View File

@ -300,7 +300,9 @@ class BlackOilModel
public: public:
BlackOilModel(Simulator& simulator) BlackOilModel(Simulator& simulator)
: ParentType(simulator) : ParentType(simulator)
{} {
eqWeights_.resize(numEq, 1.0);
}
/*! /*!
* \brief Register all run-time parameters for the immiscible model. * \brief Register all run-time parameters for the immiscible model.
@ -440,41 +442,11 @@ public:
if (globalDofIdx >= this->numGridDof()) if (globalDofIdx >= this->numGridDof())
return 1.0; return 1.0;
// we do not care much about water, so it gets de-prioritized by a factor of 100 return eqWeights_[eqIdx];
static constexpr Scalar waterPriority = 1e-2;
if constexpr (getPropValue<TypeTag, Properties::BlackoilConserveSurfaceVolume>()) {
// Roughly convert the surface volume of the fluids from m^3 to kg. (in this
// context, it does not really matter if the actual densities are off by a
// factor of two or three.)
switch (eqIdx) {
case Indices::conti0EqIdx + FluidSystem::waterCompIdx:
return 1000.0*waterPriority;
case Indices::conti0EqIdx + FluidSystem::gasCompIdx:
return 1.0;
case Indices::conti0EqIdx + FluidSystem::oilCompIdx:
return 650.0;
}
} }
if (SolventModule::eqApplies(eqIdx)) void setEqWeight(unsigned eqIdx, Scalar value) {
return SolventModule::eqWeight(eqIdx); eqWeights_[eqIdx] = value;
if (ExtboModule::eqApplies(eqIdx))
return ExtboModule::eqWeight(eqIdx);
else if (PolymerModule::eqApplies(eqIdx))
return PolymerModule::eqWeight(eqIdx);
else if (EnergyModule::eqApplies(eqIdx))
return EnergyModule::eqWeight(eqIdx);
// it is said that all kilograms are born equal (except water)!
if (eqIdx == Indices::conti0EqIdx + FluidSystem::waterCompIdx)
return waterPriority;
return 1.0;
} }
/*! /*!
@ -615,6 +587,8 @@ protected:
} }
private: private:
std::vector<Scalar> eqWeights_;
Implementation& asImp_() Implementation& asImp_()
{ return *static_cast<Implementation*>(this); } { return *static_cast<Implementation*>(this); }
const Implementation& asImp_() const const Implementation& asImp_() const