[xuegy]

Dear FOAMers,

I have a very frustrating bug in my code. I need to calculate 3 eigenvalues (saved as volVectorField) and 3 eigenvectors (saved as volTensorField) in my code. Obviously there's no direction in an eigenvalue vector. It's just 3 numbers. But OpenFOAM automatically cuts my eigenvalues and eigenvectors along the normal direction of the symmetryPlane. This is really really stupid. All I need is zeroGradient (OpenFOAM doesn't allow me to use zeroGradient on symmetryPlane boundaries).
One of the possible solution is to store 3+9=12 scalarFields, and convert it back to a vector/tensor while running. But before doing that, is there any elegant way to resolve this?

Thanks.

[HPE]

have you solved your problem? out of curiosity, how do you use eigenvalues and eigenvectors within a simulation - as a part of the flow field (or I misunderstood)?

maybe keeping them as vectorField/tensorField instead of vol*Fields could solve the issue?

[xuegy]

Originally Posted by HPE

have you solved your problem? out of curiosity, how do you use eigenvalues and eigenvectors within a simulation - as a part of the flow field (or I misunderstood)?

maybe keeping them as vectorField/tensorField instead of vol*Fields could solve the issue?

Hello.

Yes I solved this problem by not reading/saving these two fields.

The eigenvalue/vector is used in the diagonalization of the stress tensor, so the tensor can be positive-definite and more stable.

Best

[HPE]

just curiosity: would be using ::Internal help for the object to be read/written? e.g. volVectorField::Internal eigenvalues?