Commit 583e1bf7 authored by Kilian Weishaupt's avatar Kilian Weishaupt
Browse files

[cleanup] Use conv functions in 2pncmin dissolutionproblem

Use volVars.moleFraction(wPhaseIdx, NaClIdx) instead
of volVars.fluidState().moleFraction(wPhaseIdx, NaClIdx)
parent 496eb859
...@@ -343,11 +343,11 @@ public: ...@@ -343,11 +343,11 @@ public:
{ {
source = 0; source = 0;
const VolumeVariables &volVars = elemVolVars[scvIdx]; const VolumeVariables &volVars = elemVolVars[scvIdx];
Scalar moleFracNaCl_lPhase = volVars.fluidState().moleFraction(wPhaseIdx, NaClIdx); Scalar moleFracNaCl_lPhase = volVars.moleFraction(wPhaseIdx, NaClIdx);
Scalar moleFracNaCl_gPhase = volVars.fluidState().moleFraction(nPhaseIdx, NaClIdx); Scalar moleFracNaCl_gPhase = volVars.moleFraction(nPhaseIdx, NaClIdx);
Scalar massFracNaCl_Max_lPhase = this->spatialParams().SolubilityLimit(); Scalar massFracNaCl_Max_lPhase = this->spatialParams().SolubilityLimit();
Scalar moleFracNaCl_Max_lPhase = massTomoleFrac_(massFracNaCl_Max_lPhase); Scalar moleFracNaCl_Max_lPhase = massTomoleFrac_(massFracNaCl_Max_lPhase);
Scalar moleFracNaCl_Max_gPhase = moleFracNaCl_Max_lPhase / volVars.fluidState().pressure(nPhaseIdx); Scalar moleFracNaCl_Max_gPhase = moleFracNaCl_Max_lPhase / volVars.pressure(nPhaseIdx);
Scalar saltPorosity = this->spatialParams().porosityMin(element, fvGeometry, scvIdx); Scalar saltPorosity = this->spatialParams().porosityMin(element, fvGeometry, scvIdx);
// liquid phase // liquid phase
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment