CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Possible turbulence modelling bug in SRF solvers

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

LinkBack Thread Tools Search this Thread Display Modes
Old   April 27, 2012, 09:14
Default Possible turbulence modelling bug in SRF solvers
New Member
Join Date: Jun 2009
Posts: 22
Rep Power: 17
otm is on a distinguished road

I'm working a lot with development of machines that rotates. Hence, I need to use the SRF solvers provided in OF (have not had time to use the MRF solvers yet).

My issue is that I'm not sure that the effects of system rotation, on the turbulence are included properly in the SRFSimpleFoam and SRFPimpleFoam solvers. As far as I can see the turbulence model is using the "Urel" field to calculate the tensors used in the modelling (but I should admit that I often have a hard time tracking down what the code does since there are so many different source files involved). In a linear eddy viscosity model this doesn't matter because only the strain rate tensor ((dUi/dxj+dUj/dxi)/2), which is invariant to rotation is used in the modelling of the Reynolds stresses. However in e.g. the LRR model the rotation rate tensor ((dUi/dxj-dUj/dxi)/2), which is not invariant to rotation is used.

A good test case for turbulence models in rotating frames is the fully developed flow in a pipe rotating around its axis. The flow is three dimensional but still very simple. The azimuthal velocity component (in the stationary frame) should have a parabolic shaped profile in the radial direction if the Reynolds stresses are predicted correctly. With a linear eddy viscosity model the azimuthal velocity will correspond to solid body rotation (i.e. a straight line).

I have tested the LRR model with the SRFSimpleFoam solver in the rotating pipe and did not get the correct profile. So it seems as if there is a bug (or at least an inconsistent method) in the SRF solvers.

Next i will run the pipe case in a stationary frame (with a moving wall). If I can get the right profiles in this case with the LRR model, I think it can be concluded that the implementation in the SRF solvers is wrong.

I welcome input from other users on this issue, if you agree with me or not. And if I have misunderstood something, please inform me.

I'm running OF-2.1.x.

otm is offline   Reply With Quote

Old   April 28, 2012, 02:12
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,912
Rep Power: 36
alberto will become famous soon enoughalberto will become famous soon enough
As a suggestion, if the problem is confirmed, please report it on mantis:

Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats (current status:
OpenQBMM - An open-source implementation of quadrature-based moment methods.

To obtain more accurate answers, please specify the version of OpenFOAM you are using.
alberto is offline   Reply With Quote

Old   April 30, 2012, 05:52
Senior Member
Vesselin Krastev
Join Date: Jan 2010
Location: University of Tor Vergata, Rome
Posts: 368
Rep Power: 20
vkrastev is on a distinguished road
That's quite interesting...Any news about this issue? From a rapid look on Mantis' page it seems it hasn't been reported yet


vkrastev is offline   Reply With Quote

Old   May 29, 2012, 04:03
New Member
Join Date: Jun 2009
Posts: 22
Rep Power: 17
otm is on a distinguished road

I still think that it could be an inadequate implementation of the turbulence models in the SRF-solvers. However, the way I wanted to show this did not work out. The reason is that in order to correctly predict the parabolic profile of the azimuthal velocity in the rotating pipe one needs to have a model for the inter-component transfer between the Reynolds stresses i.e. use a RSTM or EARSM. Unfortunately it turns out that the Reynolds stress models implemented in OF (LRR and Launder Gibson) are based on the simplified "isotropization of production" model (often referred to as IP-LRR) for the inter-component transfer, and this model is known to be unable of predicting the azimuthal velocity profile correctly.

So, the possible bug may not have a big influence until there are more advanced models available in OF.

otm is offline   Reply With Quote


srfpimplefoam, srfsimplefoam, turbulence

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
modelling zero gradient of turbulence amirbahador FLUENT 0 September 25, 2010 06:37
Turbulence modelling kluni Main CFD Forum 6 December 1, 2009 16:46
k-ep turbulence modelling allan thomson Main CFD Forum 5 February 20, 2000 20:39
Question about Low Re turbulence modelling. ghlee Main CFD Forum 1 May 28, 1999 19:29
Multigrid Method in Turbulence Modelling Bipin Lokhande Main CFD Forum 1 November 11, 1998 19:20

All times are GMT -4. The time now is 10:59.