[david_h]

The double dot product operator seems to calculate the incorrect result for when one argument is of class "symmTensor" and the other class "tensor"

I have attached two examples:

checkTensor - calculates several operation with symmTensor and tensor's

dnsFoamCheckTensor - modification of dnsFoam which calculates S:S, S:D using S of class symmTensor and S of class tensor.

Below is some output of dnsFoamCheckTensor
S:S

S:D

S2:S2

S2:D


D = grad(U)
S = symm(D)
S2 = 0.5 (D + D^T)


checkTensors.tgz
dnsFoamCheckTensor.tgz

Dave

[henry]

Yes you are quite right there is a bug in the && operators between symmetric and non-symmetric tensors, here is the fix TensorI.H to replace the TensorI.H in OpenFOAM/OpenFOAM-1.4.1/src/OpenFOAM/primitives/Tensor