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

curvature correction term, material derivative of a tensor

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 1, 2010, 11:08
Default curvature correction term, material derivative of a tensor
  #1
New Member
 
Volker Tritschler
Join Date: Jan 2010
Posts: 20
Rep Power: 16
volker is on a distinguished road
Hi,

I'm planning to implement a curvature correction term to an explicit algebraic reynolds stress model based on the existing k-e-model of Launder and Sharma in order to sensitize the Launder-Sharma-model to streamline curvature.
The curvature correction term is to be implemented according to P.E. Smirnov and F.R. Menter "Sensitization of the SST Turbulence Model to Rotation and Curvature by Applying the Spalart-Shur Correction Term".
I ran into problems quite fast. There is a term, which is the material derivative of the shear strain tensor. The implementation is not straight-forward, especially in tensor notation. Either one has somehow to implement a summing over all faces of a control volume or to evaluate the gradient of the shear strain tensor, which leads to a tensor of third order. And right now, I do not see a way to get rid of this problem.
Does anyone has already tried to implement a curvature correction or has any experience on that? Or even know how to implement a material derivative of a tensor?

I'm glad for any hint.

Greets, volker
volker is offline   Reply With Quote

Old   April 1, 2010, 11:46
Default
  #2
Senior Member
 
wayne.zhang
Join Date: Mar 2009
Location: Shanghai, Shanghai, P.R.China
Posts: 309
Rep Power: 18
waynezw0618 is on a distinguished road
Send a message via MSN to waynezw0618 Send a message via Skype™ to waynezw0618
I have try RC correction as you mentioned in ref paper. you can send me a mail. if interesting.

my email address: waynezw0618@163.com
waynezw0618 is offline   Reply With Quote

Old   May 6, 2010, 16:55
Default Gradient of a second order tensor
  #3
New Member
 
Michael B Martell Jr
Join Date: Feb 2010
Location: Amherst, MA
Posts: 18
Rep Power: 16
theory37 is on a distinguished road
I am attempting to add a term to a Reynolds stress transport model which involves taking the gradient of the Reynolds stress tensor, which should result in a tensor of rank three.

Even a simple attempt fails:

Code:
volTensorField test = fvc::grad(R);
(If I attempt to take the gradient of a non-symmetric tensor, I still cannot compile.) In both cases, the result is the same:

Code:
error: no type named ‘type’ in ‘class Foam::typeOfRank<double, 3>’
error: no matching function for call to 'grad(Foam::GeometricField<Foam::Tensor<double>, Foam::fvPatchField, Foam::volMesh>&)’
Is it possible to take the gradient of a second order tensor in FOAM?

Any ideas or help would be appreciated!
theory37 is offline   Reply With Quote

Old   May 6, 2010, 17:55
Default
  #4
New Member
 
Volker Tritschler
Join Date: Jan 2010
Posts: 20
Rep Power: 16
volker is on a distinguished road
hi michael,

I had exactly the same problem, though I found it kind of surprising because OpenFOAM is supposed to handle tensors of rank three.
But the gradient seems to be restricted to tensors of rank one. I hope that there might be a workaround.

hope someone can help us.

cheers, volker
volker is offline   Reply With Quote

Old   May 6, 2010, 18:03
Default
  #5
New Member
 
Michael B Martell Jr
Join Date: Feb 2010
Location: Amherst, MA
Posts: 18
Rep Power: 16
theory37 is on a distinguished road
Hi Volker. I have discussed this problem with some of my colleagues here and we agree that this would be non-trivial to implement. I am going to try, and I'll keep you posted on my progress. Feel free to contact me through http://www.ecs.umass.edu/~mmartell if you'd like.

FYI, I have tried not assigning the gradient to a particular type by simple putting

Code:
fvc::grad(R);
to see if I could figure out what type FOAM was looking for. No luck.

Cheers,

Mike
theory37 is offline   Reply With Quote

Old   May 7, 2010, 08:29
Default
  #6
New Member
 
Volker Tritschler
Join Date: Jan 2010
Posts: 20
Rep Power: 16
volker is on a distinguished road
would be great, if you would come up with something new.

keep me posted!
volker is offline   Reply With Quote

Old   May 27, 2010, 04:50
Default
  #7
New Member
 
Volker Tritschler
Join Date: Jan 2010
Posts: 20
Rep Power: 16
volker is on a distinguished road
hi michael,

did you manage to implement the gradient of a tensor?

cheers, volker
volker is offline   Reply With Quote

Old   June 3, 2010, 09:08
Default
  #8
New Member
 
Michael B Martell Jr
Join Date: Feb 2010
Location: Amherst, MA
Posts: 18
Rep Power: 16
theory37 is on a distinguished road
Hi Volker,

I managed to do what I needed "by hand" but it is sloppy and not generic. I'll attach the source for it. It has no external dependencies (you can ignore the inclusion of OEC.H) and it doesn't need to be cast in the incompressible or RASModels name space - I simply took it right from my existing code.

In my case, I needed something like fvc::grad(kINT) & fvc::grad(RINT) with kINT being a volScalarField and RINT being a volSymmTensorField. Note that this inner product (of grad(kINT), a vector, and grad(RINT), a rank 3 tensor) results in a rank 2 symmetric tensor. I know this isn't exactly what you were looking for, but the structure is there.

Let me know if you need any help

Mike
Attached Files
File Type: c calcGradKgradR.C (7.7 KB, 121 views)
theory37 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
mass flow in is not equal to mass flow out saii CFX 12 March 19, 2018 06:21
Curvature Correction & Reattachment Mod. Felix CFX 5 December 19, 2010 17:04
How to use SST Curvature Correction mike wilson CFX 11 December 19, 2010 17:04
ATTENTION! Reliability problems in CFX 5.7 Joseph CFX 14 April 20, 2010 16:45
Two-Phase Buoyant Flow Issue Miguel Baritto CFX 4 August 31, 2006 13:02


All times are GMT -4. The time now is 12:44.