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

help with wall functions

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

Reply
 
LinkBack Thread Tools Display Modes
Old   February 10, 2000, 16:51
Default help with wall functions
  #1
Nick Georgiadis
Guest
 
Posts: n/a
I am trying to use the wall-function approach of Ota and Goldberg, in which none of the flow variables are reset, and in fact, a no-slip surface is used. The essence of the model is that you use a law of the wall (they propose using the White-Cristoph compressible law of the wall) to solve for the wall shear stress, and then simply use this shear stress in the solution scheme. Any turbulence model (even an algebraic model like I am using, Cebeci-Smith) may then be employed with this method.

My questions to anyone familiar with this (or a similar) approach is:

(1) When solving the law of the wall expression, the iterated quatity is friction velocity. This quantity is then squared and multiplied by the density to get the wall shear stress. How can this (especially through the squaring process) be expected to get an accurate wall shear stress?

(2) Apart from solving for the wall shear stress, the (solved) shear stress throughout the rest of the bdy. layer should be constant, particularly near the wall. Can this happen in practice when using such a wall function? especially with a higher order scheme that would use the flow variables at the wall point and one off the wall?

This wall function approach has been used with success by several authors, so I am sure my doubts are due to my ignornace on the topic. As a result, any explanations to my 2 questions will be very appreciated.
  Reply With Quote

Old   February 10, 2000, 17:56
Default Re: help with wall functions
  #2
John C. Chien
Guest
 
Posts: n/a
(1). Why not send the e-mail directly to the author(s) of the method? (2). Unless you have the name of the reference paper, it is somewhat difficult to know what you are talking about. (3). In general, the law of the wall is a one parameter family of curves, see Schlichting's Theory of Boundary Layer book. The unknown parameter is the v*, the friction velocity. It is the dimensionless wall shear stress. (4). This log-law of the wall can be derived from the assumption of constant shear stress throughout the boundary layer. I think, the book also covers this part. (5). So, the constant shear stress assumption, and the universal law of the wall are well-known. (6). Exactly how this should be used in a turbulence model varies greatly depending upon the model and the numerical method used. So, the best place is the author(s) of the model, I think.
  Reply With Quote

Old   February 19, 2000, 03:23
Default Re: help with wall functions
  #3
Bob Anderson
Guest
 
Posts: n/a
Your first question confuses me. u_\tau is by definition equal to sqrt(\tau_w/\rho). So if you can find u_\tau and you know \rho, you know \tau_w. Its that simple.

The answer to your second question is that if your pressure gradients are small, constant stress is a good assumption through the log layer, in fact the law-of-the-wall is based upon this assumption. If your pressure gradients are not small, then you have to do something else, like add a first order correction to \tau_w.

Bob
  Reply With Quote

Old   February 20, 2000, 10:28
Default Re: help with wall functions
  #4
Nick Georgiadis
Guest
 
Posts: n/a
Bob,

Thank you for your response. Regarding my first question, I am not challenging whether the formulation is physically correct, I know u-tau = sqrt(tauw / rho).....the problem is when you rewrite this, tauw = rho * (u-tau)**2 - I am wondering as to how accurate this can be expected to be in a computational sense, since you are indirectly solving for u-tau, then squaring this value (and as a result, any innacuracies of u-tau are greatly magnified) Then you use this value for the wall-shear where it very much determines the bdy. layer characteristics.

I should mention that at the time I posted the first message, my bdy. layers (Mach 1.36 and Mach 1.91 flat plates with zero pressure gradient) were not correct - but then I took the suggestion of a co-worker to completely turn off the artificial dissipation at the points near the wall, and this substantially improved the results. The behavior of the first point off the wall seems critical in getting the wall function method to work - you have to be careful with your solution scheme (how you do the differencing) and any artificial dissipation.
  Reply With Quote

Old   February 20, 2000, 18:07
Default Re: help with wall functions
  #5
Bob Anderson
Guest
 
Posts: n/a
If you converge your iteration for u_\tau to machine accuracy, then it should be no less accurate than any other squared quanitity in your calculation, of which there are likely dozens or hundreds, so I would not worry about it.

Artificial dissipation is an important issue for boundary layer calculations with or without wall functions, so yes, you do have to be careful there. Seems like you have it under control now.

Out of curiousity, how do you handle the v component of velocity?

Bob
  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
Modelling Turbulent Wall Functions victorwfreire CFX 11 January 14, 2015 16:48
LES with wall functions, wall boundary, Re 10000 subba CFX 4 July 10, 2007 05:04
Wall Functions pierre OpenFOAM Running, Solving & CFD 0 October 1, 2005 13:13
Wall functions - iterative approach agg Main CFD Forum 1 May 13, 2005 12:25
LES wall functions barrat Main CFD Forum 10 September 1, 2004 08:29


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