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

Output support for RSWSAT and RVWSAT

Browse Source
This commit is contained in:
David Landa Marban 2024-07-17 12:57:40 +02:00
parent 6e74fea547
commit 8320a9abdf
3 changed files with 25 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
opm/simulators/flow/GenericOutputBlackoilModule.cpp
View File

@ -556,8 +556,10 @@ assignToSolution(data::Solution& sol)
DataEntry{"PRESPOTF", UnitSystem::measure::pressure, mechPotentialPressForce_}, DataEntry{"PRESPOTF", UnitSystem::measure::pressure, mechPotentialPressForce_},
DataEntry{"PRES_OVB", UnitSystem::measure::pressure, overburdenPressure_}, DataEntry{"PRES_OVB", UnitSystem::measure::pressure, overburdenPressure_},
DataEntry{"RSW", UnitSystem::measure::gas_oil_ratio, rsw_}, DataEntry{"RSW", UnitSystem::measure::gas_oil_ratio, rsw_},
DataEntry{"RSWSAT", UnitSystem::measure::gas_oil_ratio, gasDissolutionFactorInWater_},
DataEntry{"RSWSOL", UnitSystem::measure::gas_oil_ratio, rswSol_}, DataEntry{"RSWSOL", UnitSystem::measure::gas_oil_ratio, rswSol_},
DataEntry{"RVW", UnitSystem::measure::oil_gas_ratio, rvw_}, DataEntry{"RVW", UnitSystem::measure::oil_gas_ratio, rvw_},
DataEntry{"RVWSAT", UnitSystem::measure::oil_gas_ratio, waterVaporizationFactor_},
DataEntry{"SALTP", UnitSystem::measure::identity, pSalt_}, DataEntry{"SALTP", UnitSystem::measure::identity, pSalt_},
DataEntry{"SS_X", UnitSystem::measure::identity, extboX_}, DataEntry{"SS_X", UnitSystem::measure::identity, extboX_},
DataEntry{"SS_Y", UnitSystem::measure::identity, extboY_}, DataEntry{"SS_Y", UnitSystem::measure::identity, extboY_},
@ -1186,6 +1188,14 @@ doAllocBuffers(const unsigned bufferSize,
rstKeywords["RVSAT"] = 0; rstKeywords["RVSAT"] = 0;
oilVaporizationFactor_.resize(bufferSize, 0.0); oilVaporizationFactor_.resize(bufferSize, 0.0);
} }
if (FluidSystem::enableDissolvedGasInWater() && rstKeywords["RSWSAT"] > 0) {
rstKeywords["RSWSAT"] = 0;
gasDissolutionFactorInWater_.resize(bufferSize, 0.0);
}
if (FluidSystem::enableVaporizedWater() && rstKeywords["RVWSAT"] > 0) {
rstKeywords["RVWSAT"] = 0;
waterVaporizationFactor_.resize(bufferSize, 0.0);
}
if (FluidSystem::phaseIsActive(waterPhaseIdx) && rstKeywords["BW"] > 0) { if (FluidSystem::phaseIsActive(waterPhaseIdx) && rstKeywords["BW"] > 0) {
rstKeywords["BW"] = 0; rstKeywords["BW"] = 0;

2
opm/simulators/flow/GenericOutputBlackoilModule.hpp
View File

@ -459,6 +459,8 @@ protected:
ScalarBuffer ppcw_; ScalarBuffer ppcw_;
ScalarBuffer gasDissolutionFactor_; ScalarBuffer gasDissolutionFactor_;
ScalarBuffer oilVaporizationFactor_; ScalarBuffer oilVaporizationFactor_;
ScalarBuffer gasDissolutionFactorInWater_;
ScalarBuffer waterVaporizationFactor_;
ScalarBuffer bubblePointPressure_; ScalarBuffer bubblePointPressure_;
ScalarBuffer dewPointPressure_; ScalarBuffer dewPointPressure_;
ScalarBuffer rockCompPorvMultiplier_; ScalarBuffer rockCompPorvMultiplier_;

13
opm/simulators/flow/OutputBlackoilModule.hpp
View File

@ -362,6 +362,19 @@ public:
fs, gasPhaseIdx, pvtRegionIdx, SoMax); fs, gasPhaseIdx, pvtRegionIdx, SoMax);
Valgrind::CheckDefined(this->oilVaporizationFactor_[globalDofIdx]); Valgrind::CheckDefined(this->oilVaporizationFactor_[globalDofIdx]);
} }
if (!this->gasDissolutionFactorInWater_.empty()) {
Scalar SwMax = elemCtx.problem().maxWaterSaturation(globalDofIdx);
this->gasDissolutionFactorInWater_[globalDofIdx]
= FluidSystem::template saturatedDissolutionFactor<FluidState, Scalar>(
fs, waterPhaseIdx, pvtRegionIdx, SwMax);
Valgrind::CheckDefined(this->gasDissolutionFactorInWater_[globalDofIdx]);
}
if (!this->waterVaporizationFactor_.empty()) {
this->waterVaporizationFactor_[globalDofIdx]
= FluidSystem::template saturatedVaporizationFactor<FluidState, Scalar>(
fs, gasPhaseIdx, pvtRegionIdx);
Valgrind::CheckDefined(this->waterVaporizationFactor_[globalDofIdx]);
}
if (!this->gasFormationVolumeFactor_.empty()) { if (!this->gasFormationVolumeFactor_.empty()) {
this->gasFormationVolumeFactor_[globalDofIdx] = 1.0 this->gasFormationVolumeFactor_[globalDofIdx] = 1.0
/ FluidSystem::template inverseFormationVolumeFactor<FluidState, Scalar>( / FluidSystem::template inverseFormationVolumeFactor<FluidState, Scalar>(