CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Bugs

Possible Bug in double dot product ampamp operator

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   September 13, 2007, 15:23
Default The double dot product operato
  #1
Member
 
E. David Huckaby
Join Date: Mar 2009
Posts: 57
Rep Power: 17
david_h is on a distinguished road
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
david_h is offline   Reply With Quote

Old   September 14, 2007, 16:58
Default Yes you are quite right there
  #2
Senior Member
 
Join Date: Mar 2009
Posts: 854
Rep Power: 22
henry is on a distinguished road
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
henry is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Operator Splitting. Maria. Main CFD Forum 5 September 17, 2005 23:10
Operator precedence hemph OpenFOAM 1 September 13, 2005 13:40
a math operator in UDF lichun Dong FLUENT 7 June 18, 2005 23:04
double fabs(double ) gts FLUENT 0 October 27, 2004 10:51
Double the # of cells->double the computing time? JB Main CFD Forum 13 April 5, 2002 19:13


All times are GMT -4. The time now is 05:16.