diff --git a/dumux/boxmodels/1p2c/1p2clocalresidual.hh b/dumux/boxmodels/1p2c/1p2clocalresidual.hh index b2b38b584d64ea08bcacef6049cf5ecd52dc13a3..66e1d0c275e3763a2c614f2f110c0aeb213d135b 100644 --- a/dumux/boxmodels/1p2c/1p2clocalresidual.hh +++ b/dumux/boxmodels/1p2c/1p2clocalresidual.hh @@ -139,10 +139,10 @@ public: { // storage term of continuity equation - massfractions result[contiEqIdx] += - volVars.density()*volVars.porosity(); + volVars.fluidState().density(/*phaseIdx*/0)*volVars.porosity(); //storage term of the transport equation - massfractions result[transEqIdx] += - volVars.density() * volVars.massFrac(comp1Idx) * volVars.porosity(); + volVars.fluidState().density(/*phaseIdx*/0) * volVars.fluidState().massFraction(/*phaseIdx*/0, comp1Idx) * volVars.porosity(); } else { @@ -151,7 +151,7 @@ public: result[contiEqIdx] += volVars.molarDensity()*volVars.porosity(); // storage term of the transport equation - molefractions result[transEqIdx] += - volVars.molarDensity()*volVars.moleFrac(comp1Idx) * + volVars.fluidState().molarDensity(/*phaseIdx*/0)*volVars.fluidState().moleFraction(/*phaseIdx*/0, comp1Idx) * volVars.porosity(); } @@ -210,9 +210,9 @@ public: // advective flux of the second component - massfraction flux[transEqIdx] += fluxVars.KmvpNormal() * - (( upwindWeight_)*up.density() * up.massFrac(comp1Idx)/up.viscosity() + (( upwindWeight_)*up.fluidState().density(/*phaseIdx=*/0) * up.fluidState().massFraction(/*phaseIdx=*/0, comp1Idx)/up.viscosity() + - (1 - upwindWeight_)*dn.density()*dn.massFrac(comp1Idx)/dn.viscosity()); + (1 - upwindWeight_)*dn.fluidState().density(/*phaseIdx=*/0)*dn.fluidState().massFraction(/*phaseIdx=*/0, comp1Idx)/dn.viscosity()); } else { @@ -227,9 +227,9 @@ public: // advective flux of the second component -molefraction flux[transEqIdx] += fluxVars.KmvpNormal() * - (( upwindWeight_)*up.molarDensity() * up.moleFrac(comp1Idx)/up.viscosity() - + - (1 - upwindWeight_)*dn.molarDensity() * dn.moleFrac(comp1Idx)/dn.viscosity()); + (( upwindWeight_)*up.molarDensity() * up.fluidState().moleFraction(/*phaseIdx=*/0, comp1Idx)/up.viscosity() + + + (1 - upwindWeight_)*dn.molarDensity() * dn.fluidState().moleFraction(/*phaseIdx=*/0, comp1Idx)/dn.viscosity()); } } @@ -411,8 +411,9 @@ protected: if(!useMoles)//use massfractions { // advective flux - values[transEqIdx]+= boundaryVars.KmvpNormal()*vertVars.density()/vertVars.viscosity() - *vertVars.fluidState().massFrac(phaseIdx, comp1Idx); + values[transEqIdx]+= + boundaryVars.KmvpNormal()*vertVars.density()/vertVars.viscosity() + *vertVars.fluidState().massFraction(phaseIdx, comp1Idx); // diffusive flux of comp1 component in phase0 Scalar tmp = 0; @@ -426,8 +427,9 @@ protected: else //use molefractions { // advective flux - values[transEqIdx]+= boundaryVars.KmvpNormal()*vertVars.molarDensity()/vertVars.viscosity() - *vertVars.fluidState().moleFrac(phaseIdx, comp1Idx); + values[transEqIdx]+= + boundaryVars.KmvpNormal()*vertVars.molarDensity()/vertVars.viscosity() + *vertVars.fluidState().moleFraction(phaseIdx, comp1Idx); // diffusive flux of comp1 component in phase0 Scalar tmp = 0; diff --git a/dumux/boxmodels/MpNc/MpNcvolumevariables.hh b/dumux/boxmodels/MpNc/MpNcvolumevariables.hh index de56494757e76b45caf50f1b60829e5b6166f538..c98ae3a50100af564c8a649c4df85697c80a68da 100644 --- a/dumux/boxmodels/MpNc/MpNcvolumevariables.hh +++ b/dumux/boxmodels/MpNc/MpNcvolumevariables.hh @@ -151,7 +151,6 @@ public: elemGeom, scvIdx, problem); - EnergyVolumeVariables::checkDefinedTemp(); ///////////// // set the phase pressures