diff --git a/applications/solvers/multiphase/interFoam/LTSInterFoam/setrDeltaT.H b/applications/solvers/multiphase/interFoam/LTSInterFoam/setrDeltaT.H index f419153506c1ef5962aad3693546bf00ee6a8ae1..3f8ba982cf9c644a203541499d90fd9aa7f51520 100644 --- a/applications/solvers/multiphase/interFoam/LTSInterFoam/setrDeltaT.H +++ b/applications/solvers/multiphase/interFoam/LTSInterFoam/setrDeltaT.H @@ -121,9 +121,11 @@ && runTime.timeIndex() > runTime.startTimeIndex() + 1 ) { - rDeltaT = - rDeltaT0 - *max(rDeltaT/rDeltaT0, scalar(1.0) - rDeltaTDampingCoeff); + rDeltaT = max + ( + rDeltaT, + (scalar(1.0) - rDeltaTDampingCoeff)*rDeltaT0 + ); Info<< "Damped flow time scale min/max = " << gMin(1/rDeltaT.internalField())