Mesh dependency of gradient scheme
 

Hello;

I have been writing a solver based on stressedFoam, which needs to evaluate the gradient of field U for the governing equations. I have noticed that if I transform a 2d mesh from xy plane to xz or yz, I get a "Floating point exception" as the solver is trying to calculate the gradient (fvc::grad(U)) when the leastSquares method is set. However if Gauss linear method is chosen the solver works fine. Hence my question is, is there a dependence of the least squares method on the mesh orientation?

Best;

Bercan