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

Transformation of volSymmTensorField to volTensorField

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

Reply
 
LinkBack Thread Tools Display Modes
Old   August 20, 2009, 22:20
Default Transformation of volSymmTensorField to volTensorField
  #1
Member
 
Sven Schweikert
Join Date: Jun 2009
Posts: 38
Rep Power: 7
svens is on a distinguished road
Hello together

Is there a way to transform a existing volSymmTensorField to a volTensorField?

I want to implement a new turbulence model which is based on the use of the symmetric part of the mean velocity gradient Sij as well as the skew-symmetric ones Wij.

In my eyes - defining the symmetric part Sij is no big deal:
Code:
 volSymmTensorField S = symm(fvc::grad(U_));
but the skew-symmetric part Wij has to be a tensor with all 9 componets (I guess...?). I defined it like this:
Code:
 volTensorField W = skew(fvc::grad(U_));
Bringing then all together leads to an error while compiling. The operator '==' in the REqn (in line: 431 - I used the LaunderGibsonRSTM as basis for my new model - switched from tmp<fvSymmTensorMatrix> to tmp<fvTensorMatrix> REqn) couldn't been recognized anymore.
So far I guess it's because of the different number of components I used while defining with volSymmTensorField and volTensorField.

Therefore my quetion: Is there a way to transform a existing volSymmTensorField to a volTensorField?
But maybe there is a better way to deal with my problem? I'm grateful about every input!

Thanks a lot & greetings
Sven

Last edited by svens; August 21, 2009 at 02:10.
svens is offline   Reply With Quote

Old   September 4, 2009, 01:28
Default
  #2
Member
 
Sven Schweikert
Join Date: Jun 2009
Posts: 38
Rep Power: 7
svens is on a distinguished road
Hello together

It doesn't really answer my question but I found a way to solve my problem:

I had a further look at my eqautions and recognized that I'm multiply my new skew-symmetric tensor W with a symmetric ones b. After that I'm summing this product like (bW + Wb) which creates a symmetric tensor again. While implementing it like symm(bW+Wb) to the REqn all my problems with tmp<fvSymmTensorMatrix> dissapeared.

Maybe someone of you can help this.

Cheers - Sven
svens is offline   Reply With Quote

Old   November 19, 2009, 05:47
Default Q&D fix
  #3
New Member
 
Iyer Arvind Sundaram
Join Date: May 2009
Posts: 12
Rep Power: 7
iyer.arvind is on a distinguished road
A quick and dirty fix.
Just multiply with tensor(1,0,0,0,1,0,0,0,1)

too bad but if performance is not a concern go ahead, kill the CPU
iyer.arvind is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Transformation of a cartesian velocity gradient volTensorField to cylindrical coord. alexey7783 OpenFOAM 5 April 26, 2012 18:29
grid transformation lavan Main CFD Forum 1 October 22, 2005 04:57
elliptic transformation in 3D Dinesh Godavarty Main CFD Forum 1 August 21, 2003 14:17
Roberts transformation Eqn. Yang Main CFD Forum 2 January 10, 2002 09:20
transformation between Cartesian and spherical coordinates? Xiaofeng Yang Main CFD Forum 6 February 23, 2000 10:19


All times are GMT -4. The time now is 19:51.