[uyelken]

Hi guys,

While im solving on ddtfoam(custom solver) with openfoam 2.1.1. div(Ut) returns inf value. How can i solve it?

volVectorField Ut = U + sL*Xi*n;
volScalarField Xisigmat = (n & n)*fvc::div(Ut) - (n & fvc::grad(Ut) & n);