Skip to content
Snippets Groups Projects
Commit 839ea2c7 authored by Andrew Heather's avatar Andrew Heather
Browse files

corrections

parent 531e22bc
Branches
Tags
No related merge requests found
......@@ -60,7 +60,7 @@ Foam::Analytical<Type>::integrate
) const
{
typename IntegrationScheme<Type>::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;
......
......@@ -60,7 +60,7 @@ Foam::Euler<Type>::integrate
) const
{
typename IntegrationScheme<Type>::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;
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment