CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Bugs (http://www.cfd-online.com/Forums/openfoam-bugs/)
-   -   Possible Bug in double dot product ampamp operator (http://www.cfd-online.com/Forums/openfoam-bugs/62544-possible-bug-double-dot-product-ampamp-operator.html)

david_h September 13, 2007 14:23

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

henry September 14, 2007 15:58

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


All times are GMT -4. The time now is 15:18.