400 using ConstEvaluation =
typename std::remove_reference<typename FluidState::Scalar>::type;
401 using FsEvaluation =
typename std::remove_const<ConstEvaluation>::type;
404 bool gasPresent = FluidSystem::phaseIsActive(gasPhaseIdx)?(fluidState.saturation(gasPhaseIdx) > 0.0):
false;
405 bool oilPresent = FluidSystem::phaseIsActive(oilPhaseIdx)?(fluidState.saturation(oilPhaseIdx) > 0.0):
false;
406 bool waterPresent = FluidSystem::phaseIsActive(waterPhaseIdx)?(fluidState.saturation(waterPhaseIdx) > 0.0):
false;
407 const auto& saltSaturation = BlackOil::getSaltSaturation_<FluidSystem, FluidState, Scalar>(fluidState, pvtRegionIdx_);
417 primaryVarsMeaningPressure_ = PressureMeaning::Pg;
418 }
else if (FluidSystem::phaseIsActive(oilPhaseIdx)) {
419 primaryVarsMeaningPressure_ = PressureMeaning::Po;
421 primaryVarsMeaningPressure_ = PressureMeaning::Pw;
422 }
else if (FluidSystem::phaseIsActive(gasPhaseIdx)) {
423 primaryVarsMeaningPressure_ = PressureMeaning::Pg;
425 assert(FluidSystem::phaseIsActive(waterPhaseIdx));
426 primaryVarsMeaningPressure_ = PressureMeaning::Pw;
434 primaryVarsMeaningWater_ = WaterMeaning::Sw;
435 }
else if (
gasPresent && FluidSystem::enableVaporizedWater()) {
436 primaryVarsMeaningWater_ = WaterMeaning::Rvw;
437 }
else if (
waterPresent && FluidSystem::enableDissolvedGasInWater()) {
438 primaryVarsMeaningWater_ = WaterMeaning::Rsw;
439 }
else if (FluidSystem::phaseIsActive(waterPhaseIdx) && !
oneActivePhases) {
440 primaryVarsMeaningWater_ = WaterMeaning::Sw;
442 primaryVarsMeaningWater_ = WaterMeaning::Disabled;
451 primaryVarsMeaningGas_ = GasMeaning::Sg;
452 }
else if (
oilPresent && FluidSystem::enableDissolvedGas()) {
453 primaryVarsMeaningGas_ = GasMeaning::Rs;
454 }
else if (
gasPresent && FluidSystem::enableVaporizedOil()){
455 primaryVarsMeaningGas_ = GasMeaning::Rv;
456 }
else if (FluidSystem::phaseIsActive(gasPhaseIdx) && FluidSystem::phaseIsActive(oilPhaseIdx)) {
457 primaryVarsMeaningGas_ = GasMeaning::Sg;
459 primaryVarsMeaningGas_ = GasMeaning::Disabled;
463 if constexpr (enableSaltPrecipitation){
465 primaryVarsMeaningBrine_ = BrineMeaning::Sp;
467 primaryVarsMeaningBrine_ = BrineMeaning::Cs;
469 primaryVarsMeaningBrine_ = BrineMeaning::Disabled;
474 case PressureMeaning::Po:
475 this->setScaledPressure_(FsToolbox::value(fluidState.pressure(oilPhaseIdx)));
477 case PressureMeaning::Pg:
478 this->setScaledPressure_(FsToolbox::value(fluidState.pressure(gasPhaseIdx)));
480 case PressureMeaning::Pw:
481 this->setScaledPressure_(FsToolbox::value(fluidState.pressure(waterPhaseIdx)));
484 throw std::logic_error(
"No valid primary variable selected for pressure");
487 case WaterMeaning::Sw:
489 (*this)[waterSwitchIdx] = FsToolbox::value(fluidState.saturation(waterPhaseIdx));
492 case WaterMeaning::Rvw:
494 const auto& rvw = BlackOil::getRvw_<FluidSystem, FluidState, Scalar>(fluidState, pvtRegionIdx_);
495 (*this)[waterSwitchIdx] = rvw;
498 case WaterMeaning::Rsw:
500 const auto& Rsw = BlackOil::getRsw_<FluidSystem, FluidState, Scalar>(fluidState, pvtRegionIdx_);
501 (*this)[waterSwitchIdx] = Rsw;
504 case WaterMeaning::Disabled:
509 throw std::logic_error(
"No valid primary variable selected for water");
514 (*this)[compositionSwitchIdx] = FsToolbox::value(fluidState.saturation(gasPhaseIdx));
519 const auto& rs = BlackOil::getRs_<FluidSystem, FluidState, Scalar>(fluidState, pvtRegionIdx_);
520 (*this)[compositionSwitchIdx] = rs;
525 const auto& rv = BlackOil::getRv_<FluidSystem, FluidState, Scalar>(fluidState, pvtRegionIdx_);
526 (*this)[compositionSwitchIdx] = rv;
529 case GasMeaning::Disabled:
534 throw std::logic_error(
"No valid primary variable selected for composision");
569 Scalar saltConcentration = 0.0;
570 const Scalar& T = asImp_().temperature_(problem,
globalDofIdx);
572 sw = (*this)[waterSwitchIdx];
574 sg = (*this)[compositionSwitchIdx];
583 if constexpr (enableSaltPrecipitation) {
585 if (primaryVarsMeaningBrine() == BrineMeaning::Sp) {
586 saltConcentration = saltSolubility;
587 Scalar
saltSat = (*this)[saltConcentrationIdx];
590 (*this)[saltConcentrationIdx] = saltSolubility;
593 else if (primaryVarsMeaningBrine() == BrineMeaning::Cs) {
594 saltConcentration = (*this)[saltConcentrationIdx];
595 if (saltConcentration > saltSolubility +
eps){
597 (*this)[saltConcentrationIdx] = 0.0;
605 if constexpr (enableSolvent) {
606 if (SolventModule::isSolubleInWater()) {
607 Scalar
p = (*this)[pressureSwitchIdx];
609 if (primaryVarsMeaningSolvent() == SolventMeaning::Ss) {
610 Scalar
solSat = (*this)[solventSaturationIdx];
613 (*this)[solventSaturationIdx] =
solLimit;
617 else if (primaryVarsMeaningSolvent() == SolventMeaning::Rsolw) {
618 Scalar
rsolw = (*this)[solventSaturationIdx];
621 (*this)[solventSaturationIdx] = 0.0;
639 if constexpr (waterEnabled) {
640 (*this)[Indices::waterSwitchIdx] = std::min(
swMaximum,
sw);
644 if constexpr (compositionSwitchEnabled)
645 (*
this)[Indices::compositionSwitchIdx] = 0.0;
649 if constexpr (compositionSwitchEnabled)
657 if (BrineModule::hasPcfactTables() && primaryVarsMeaningBrine() == BrineMeaning::Sp) {
659 Scalar Sp = saltConcentration_();
669 case WaterMeaning::Sw:
673 Scalar
p = this->pressure_();
675 std::array<Scalar, numPhases>
pC = { 0.0 };
677 Scalar
so = 1.0 -
sg - solventSaturation_();
678 computeCapillaryPressures_(
pC,
so,
sg + solventSaturation_(), 0.0,
matParams);
679 p += pcFactor_ * (
pC[gasPhaseIdx] -
pC[oilPhaseIdx]);
681 Scalar
rvwSat = FluidSystem::gasPvt().saturatedWaterVaporizationFactor(pvtRegionIdx_,
686 (*this)[Indices::waterSwitchIdx] =
rvwSat;
693 const Scalar pg = this->pressure_();
695 std::array<Scalar, numPhases>
pC = { 0.0 };
697 Scalar
so = 1.0 -
sw - solventSaturation_();
699 Scalar pw = pg + pcFactor_ * (
pC[waterPhaseIdx] -
pC[gasPhaseIdx]);
700 Scalar
rswSat = FluidSystem::waterPvt().saturatedGasDissolutionFactor(pvtRegionIdx_,
708 this->setScaledPressure_(pw);
714 case WaterMeaning::Rvw:
716 const Scalar& rvw = (*this)[waterSwitchIdx];
717 Scalar
p = this->pressure_();
719 std::array<Scalar, numPhases>
pC = { 0.0 };
721 Scalar
so = 1.0 -
sg - solventSaturation_();
722 computeCapillaryPressures_(
pC,
so,
sg + solventSaturation_(), 0.0,
matParams);
723 p += pcFactor_ * (
pC[gasPhaseIdx] -
pC[oilPhaseIdx]);
725 Scalar
rvwSat = FluidSystem::gasPvt().saturatedWaterVaporizationFactor(pvtRegionIdx_,
732 (*this)[Indices::waterSwitchIdx] = 0.0;
737 case WaterMeaning::Rsw:
742 const Scalar& pw = this->pressure_();
744 Scalar
rswSat = FluidSystem::waterPvt().saturatedGasDissolutionFactor(pvtRegionIdx_,
749 Scalar rsw = (*this)[Indices::waterSwitchIdx];
754 (*this)[Indices::waterSwitchIdx] = 1.0;
756 std::array<Scalar, numPhases>
pC = { 0.0 };
758 computeCapillaryPressures_(
pC, 0.0, 0.0, 1.0,
matParams);
759 Scalar pg = pw + pcFactor_ * (
pC[gasPhaseIdx] -
pC[waterPhaseIdx]);
760 this->setScaledPressure_(pg);
765 case WaterMeaning::Disabled:
770 throw std::logic_error(
"No valid primary variable selected for water");
784 Scalar s = 1.0 -
sw - solventSaturation_();
786 const Scalar po = this->pressure_();
793 : FluidSystem::oilPvt().saturatedGasDissolutionFactor(pvtRegionIdx_,
798 (*this)[Indices::compositionSwitchIdx] = std::min(
rsMax,
rsSat);
801 Scalar
so = 1.0 -
sw - solventSaturation_() -
sg;
807 const Scalar po = this->pressure_();
808 std::array<Scalar, numPhases>
pC = { 0.0 };
810 computeCapillaryPressures_(
pC, 0.0,
sg + solventSaturation_(),
sw,
matParams);
811 Scalar pg = po + pcFactor_ * (
pC[gasPhaseIdx] -
pC[oilPhaseIdx]);
816 this->setScaledPressure_(pg);
822 : FluidSystem::gasPvt().saturatedOilVaporizationFactor(pvtRegionIdx_,
828 (*this)[Indices::compositionSwitchIdx] = std::min(
rvMax,
rvSat);
838 const Scalar po = this->pressure_();
839 Scalar
so = 1.0 -
sw - solventSaturation_();
845 : FluidSystem::oilPvt().saturatedGasDissolutionFactor(pvtRegionIdx_,
851 Scalar rs = (*this)[Indices::compositionSwitchIdx];
855 (*this)[Indices::compositionSwitchIdx] = 0.0;
866 const Scalar pg = this->pressure_();
872 : FluidSystem::gasPvt().saturatedOilVaporizationFactor(pvtRegionIdx_,
878 Scalar rv = (*this)[Indices::compositionSwitchIdx];
883 Scalar
sg2 = 1.0 -
sw - solventSaturation_();
884 std::array<Scalar, numPhases>
pC = { 0.0 };
886 computeCapillaryPressures_(
pC,
888 sg2 + solventSaturation_(),
891 Scalar po = pg + pcFactor_ * (
pC[oilPhaseIdx] -
pC[gasPhaseIdx]);
895 this->setScaledPressure_(po);
896 (*this)[Indices::compositionSwitchIdx] =
sg2;
901 case GasMeaning::Disabled:
906 throw std::logic_error(
"No valid primary variable selected for water");