diff --git a/test/porousmediumflow/2pncmin/implicit/problem.hh b/test/porousmediumflow/2pncmin/implicit/problem.hh index 1fc9655b1591ed552a9f96544dadfac0d6abaee7..2b9941d2e4f79b7384c60270b4cab8363e908447 100644 --- a/test/porousmediumflow/2pncmin/implicit/problem.hh +++ b/test/porousmediumflow/2pncmin/implicit/problem.hh @@ -368,10 +368,8 @@ public: const auto& volVars = elemVolVars[scv]; Scalar moleFracNaCl_wPhase = volVars.moleFraction(liquidPhaseIdx, NaClIdx); - Scalar moleFracNaCl_nPhase = volVars.moleFraction(gasPhaseIdx, NaClIdx); Scalar massFracNaCl_Max_wPhase = this->spatialParams().solubilityLimit(); Scalar moleFracNaCl_Max_wPhase = massToMoleFrac_(massFracNaCl_Max_wPhase); - Scalar moleFracNaCl_Max_nPhase = moleFracNaCl_Max_wPhase / volVars.pressure(gasPhaseIdx); Scalar saltPorosity = this->spatialParams().minimalPorosity(element, scv); // liquid phase @@ -382,11 +380,6 @@ public: if (moleFracNaCl_wPhase < moleFracNaCl_Max_wPhase) precipSalt *= -1; - // gas phase - precipSalt += volVars.porosity() * volVars.molarDensity(gasPhaseIdx) - * volVars.saturation(gasPhaseIdx) - * abs(moleFracNaCl_nPhase - moleFracNaCl_Max_nPhase); - // make sure we don't dissolve more salt than previously precipitated if (precipSalt*timeStepSize_ + volVars.solidVolumeFraction(sPhaseIdx)* volVars.solidComponentMolarDensity(sPhaseIdx)< 0) precipSalt = -volVars.solidVolumeFraction(sPhaseIdx)* volVars.solidComponentMolarDensity(sPhaseIdx)/timeStepSize_;