The double dot product operato
 

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
http://www.cfd-online.com/OpenFOAM_D...s/126/5400.png
S:D
http://www.cfd-online.com/OpenFOAM_D...s/126/5401.png
S2:S2
http://www.cfd-online.com/OpenFOAM_D...s/126/5402.png
S2:D
http://www.cfd-online.com/OpenFOAM_D...s/126/5403.png

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


http://www.cfd-online.com/OpenFOAM_D...hment_icon.gif checkTensors.tgz
http://www.cfd-online.com/OpenFOAM_D...hment_icon.gif dnsFoamCheckTensor.tgz

Dave

Yes you are quite right there
 

Yes you are quite right there is a bug in the && operators between symmetric and non-symmetric tensors, here is the fix http://www.cfd-online.com/OpenFOAM_D...hment_icon.gif TensorI.H to replace the TensorI.H in OpenFOAM/OpenFOAM-1.4.1/src/OpenFOAM/primitives/Tensor