diff --git a/src/lagrangian/intermediate/submodels/Reacting/PhaseChangeModel/LiquidEvaporation/LiquidEvaporation.C b/src/lagrangian/intermediate/submodels/Reacting/PhaseChangeModel/LiquidEvaporation/LiquidEvaporation.C index 61dfb6c6ee..b387a99ea7 100644 --- a/src/lagrangian/intermediate/submodels/Reacting/PhaseChangeModel/LiquidEvaporation/LiquidEvaporation.C +++ b/src/lagrangian/intermediate/submodels/Reacting/PhaseChangeModel/LiquidEvaporation/LiquidEvaporation.C @@ -155,7 +155,7 @@ void Foam::LiquidEvaporation::calculate const scalarField X(liquids_.X(Yl)); // immediately evaporate mass that has reached critical condition - if (mag(T - liquids_.Tc(X)) < SMALL) + if ((liquids_.Tc(X) - T) < SMALL) { if (debug) { diff --git a/src/lagrangian/intermediate/submodels/Reacting/PhaseChangeModel/LiquidEvaporationBoil/LiquidEvaporationBoil.C b/src/lagrangian/intermediate/submodels/Reacting/PhaseChangeModel/LiquidEvaporationBoil/LiquidEvaporationBoil.C index d78d6b110b..adde0740a0 100644 --- a/src/lagrangian/intermediate/submodels/Reacting/PhaseChangeModel/LiquidEvaporationBoil/LiquidEvaporationBoil.C +++ b/src/lagrangian/intermediate/submodels/Reacting/PhaseChangeModel/LiquidEvaporationBoil/LiquidEvaporationBoil.C @@ -155,7 +155,7 @@ void Foam::LiquidEvaporationBoil::calculate const scalarField X(liquids_.X(Yl)); // immediately evaporate mass that has reached critical condition - if (mag(T - liquids_.Tc(X)) < SMALL) + if ((liquids_.Tc(X) - T) < SMALL) { if (debug) {