diff --git a/src/lagrangian/intermediate/IntegrationScheme/Analytical/Analytical.C b/src/lagrangian/intermediate/IntegrationScheme/Analytical/Analytical.C index 55881be097..cb73887cb4 100644 --- a/src/lagrangian/intermediate/IntegrationScheme/Analytical/Analytical.C +++ b/src/lagrangian/intermediate/IntegrationScheme/Analytical/Analytical.C @@ -60,7 +60,7 @@ Foam::Analytical::integrate ) const { typename IntegrationScheme::IntegrationResult retValue; - retValue.average() = alpha + (phi - alpha)*(1-exp(-beta*dt))/(beta*dt); + retValue.average() = alpha + (phi - alpha)*(1 - exp(-beta*dt))/(beta*dt); retValue.value() = alpha + (phi - alpha)*exp(-beta*dt); return retValue; diff --git a/src/lagrangian/intermediate/IntegrationScheme/Euler/Euler.C b/src/lagrangian/intermediate/IntegrationScheme/Euler/Euler.C index 3897cd44f3..e56eb6e8d3 100644 --- a/src/lagrangian/intermediate/IntegrationScheme/Euler/Euler.C +++ b/src/lagrangian/intermediate/IntegrationScheme/Euler/Euler.C @@ -60,7 +60,7 @@ Foam::Euler::integrate ) const { typename IntegrationScheme::IntegrationResult retValue; - retValue.value() = (phi + dt*alpha)/(1.0 + dt/beta); + retValue.value() = (phi + beta*dt*alpha)/(1.0 + beta*dt); retValue.average() = 0.5*(phi + retValue.value()); return retValue;