Commit fb332ca5 authored by Andrew Heather's avatar Andrew Heather
Browse files

ENH: Added stabilisation to parcel analaytical integration scheme

parent c4db7d05
......@@ -69,7 +69,8 @@ Foam::Analytical<Type>::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;
......
Supports Markdown
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