Skip to content
Snippets Groups Projects
Commit f4b66927 authored by Dennis Gläser's avatar Dennis Gläser
Browse files

Merge branch 'cleanup/FS211-UseConvFunc' into 'master'

Cleanup/fs211 use conv func



See merge request !57
parents 129978c0 fc2ab1f1
No related branches found
No related tags found
No related merge requests found
...@@ -227,10 +227,10 @@ private: ...@@ -227,10 +227,10 @@ private:
// calculate the phase density at the integration point. we // calculate the phase density at the integration point. we
// only do this if the wetting phase is present in both cells // only do this if the wetting phase is present in both cells
Scalar SI = elemVolVars[this->face().i].fluidState().saturation(phaseIdx); Scalar SI = elemVolVars[this->face().i].saturation(phaseIdx);
Scalar SJ = elemVolVars[this->face().j].fluidState().saturation(phaseIdx); Scalar SJ = elemVolVars[this->face().j].saturation(phaseIdx);
Scalar rhoI = elemVolVars[this->face().i].fluidState().density(phaseIdx); Scalar rhoI = elemVolVars[this->face().i].density(phaseIdx);
Scalar rhoJ = elemVolVars[this->face().j].fluidState().density(phaseIdx); Scalar rhoJ = elemVolVars[this->face().j].density(phaseIdx);
Scalar fI = std::max(0.0, std::min(SI/1e-5, 0.5)); Scalar fI = std::max(0.0, std::min(SI/1e-5, 0.5));
Scalar fJ = std::max(0.0, std::min(SJ/1e-5, 0.5)); Scalar fJ = std::max(0.0, std::min(SJ/1e-5, 0.5));
if (Dune::FloatCmp::eq<Scalar, Dune::FloatCmp::absolute>(fI + fJ, 0.0, 1.0e-30)) if (Dune::FloatCmp::eq<Scalar, Dune::FloatCmp::absolute>(fI + fJ, 0.0, 1.0e-30))
......
...@@ -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
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment