[hrushi.397]

Dear all,

I have developed a variant of interFOAM solver that solves multiphase flows. It uses PIMPLE algorithm to solve pressure-velocity coupling of Navier-Stokes solver. While my divergence of surface flux is very low (fvc::div(phi) is of the order of 10^-12) and my global continuity as spat out by continuityErr.H is of the order of 10^-20-10^-18; my divergence in velocity is very high (fvc::div(U) is of the order of 10^2).

I was wondering if there are any settings available in the OpenFOAM that can reduce the divergence of velocity U. I do not use div(U) at any point in my solution of the Navier-Stokes solver, but I would still like to reduce my divergence in velocity. All the current settings can only affect calculation of phi and not the calculation of div(U) directly, so I am really confused where to start modifying the solver settings.

Regards,

Hrushi