diff --git a/src/lagrangian/intermediate/IntegrationScheme/Analytical/Analytical.C b/src/lagrangian/intermediate/IntegrationScheme/Analytical/Analytical.C index 7075f2152a..fa96bdec0a 100644 --- a/src/lagrangian/intermediate/IntegrationScheme/Analytical/Analytical.C +++ b/src/lagrangian/intermediate/IntegrationScheme/Analytical/Analytical.C @@ -69,7 +69,8 @@ Foam::Analytical::integrate const scalar expTerm = exp(min(50, -beta*dt)); - retValue.average() = alpha + (phi - alpha)*(1 - expTerm)/(beta*dt); + retValue.average() = + alpha + (phi - alpha)*(1 - expTerm)/(beta*dt + ROOTVSMALL); retValue.value() = alpha + (phi - alpha)*expTerm; return retValue;