[andersson.j]

Hello everybody

I'm using interFoam to generate waves and having some troubles with periodic intervals of high velocities yielding very small time steps. These high velocities are not present during post processing. However the time step size is calculated using the flux field phi where these artificial high velocities appear.


Could someone please explain why do we have the extra term added to phiU (in bold)?

from pEqn.H:

 
phi = phiU +
    (
        fvc::interpolate(interface.sigmaK())*fvc::snGrad(gamma)
      - ghf*fvc::snGrad(rho)
    )*rUAf*mesh.magSf();

Shouldn't phi be "reset" (equal to U*faceArea) before include CourantNo.H?

I'm using OF1.5dev.

Thanks in advance,
Jonas