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

minimum formulation in function calcUTau (S-A wall function)

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

LinkBack Thread Tools Display Modes
Old   May 8, 2010, 19:21
Default minimum formulation in function calcUTau (S-A wall function)
New Member
Join Date: Jan 2010
Posts: 12
Rep Power: 9
Luise is on a distinguished road

I took a closer look at nutSpalartAllmarasWallFunction.C. In function calcUTau the calculation of kUu reads

kUu = min(50, kappa*Up/ut)

Why is this minimum formulation introduced here and why is it restricting to the particular value of 50?
Maybe it should keep the argument of the taylor expansion in a converging range, but 50 seems much too large for this. Does anybody have an idea?

Thanks in advance,
Luise is offline   Reply With Quote

Old   May 11, 2010, 03:38
New Member
Join Date: Jan 2010
Posts: 12
Rep Power: 9
Luise is on a distinguished road
Nobody knows anything?

In my calculations (a simple airfoil at 5deg incidence) the velocity profiles (sampled on the suction side at x/chord={0.07, 0.17, 0.3}) in the boundary layer show a strange behaviour, which seems to be a rather systematic problem:

For the cell closest to the wall the yPlus/uPlus mark is glued to the viscous sublayer, i.e. uPlus=yPlus.

Now, the Newtons method for finding the zero of f(uTau) is a gradient based method, while the function f shows a hyperbolic behaviour. So it might be better to engage other methods, which do not use gradients.
Otherwise, if all this has its reasons in some other part of the wallfunction or whatever (maybe the restriction of kUu), working on this would be a waste of time.

Thats what I figured!
Maybe someone knows better? I'm looking forward for suggestions!

Attached Images
File Type: gif Gs.07.gif (4.9 KB, 9 views)
File Type: gif Gs.17.gif (4.9 KB, 5 views)
File Type: gif Gs.30.gif (4.9 KB, 6 views)
Luise is offline   Reply With Quote


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
BlockMesh FOAM warning gaottino OpenFOAM Native Meshers: blockMesh 7 July 19, 2010 14:11
Problem with compile the setParabolicInlet ivanyao OpenFOAM Running, Solving & CFD 6 September 5, 2008 20:50
Wall function problem in Fluent mefpz FLUENT 1 October 10, 2007 13:43
wall function Alex CFX 9 May 19, 2005 10:05
Wall function in adverse pressure gradients stephane baralon Main CFD Forum 11 September 2, 1999 04:05

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