GitLab

Using OpenFOAM 1706, with pimpleFoam, and the turbulentDFSEMInlet boundary condition, defined for velocity as below.

inlet
{
    type            turbulentDFSEMInlet;
    value           uniform (0 18 0);
	delta           1;
	R               uniform (3.2 0.1 0.1 3.2 0.1 3.2);
	L               uniform 10;
	U		    uniform (0 10 0);
	mapMethod       nearestCell;
}

Results in warning: --> FOAM Warning : From function Foam::vector Foam::eigenValues(const tensor&) in file primitives/Tensor/tensor/tensor.C at line 100 Complex eigenvalues detected for tensor: (3.2 0.1 0.1 0.1 3.2 0.1 0.1 0.1 3.2)

However, eigenvalues of this tensor are real (3.1, 3.1, 3.4).

The code in tensor.C appears correct, so I suspect the problem may lie in the cubicEqn code, which is significantly more opaque.