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

help with wall functions

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

Reply
 
LinkBack Thread Tools Search this Thread 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

Old   January 10, 2017, 04:52
Default
  #6
Senior Member
 
Ashwani
Join Date: Sep 2013
Location: Hyderabad
Posts: 154
Rep Power: 12
AshwaniAssam is on a distinguished road
I hope the problem was solved. Can you please suggest how was pressure gradient calculated? Did you tried to solve this way for compression ramp problem?
AshwaniAssam is offline   Reply With Quote

Old   January 11, 2017, 23:09
Default
  #7
Member
 
Jingchang.Shi
Join Date: Aug 2012
Location: Hang Zhou, China
Posts: 78
Rep Power: 13
aerosjc is on a distinguished road
Are you working on wall model for hybrid LES/RANS? And try to account for separation in flow over ramp? Do you try just using pressure gradient calculated in your program?
aerosjc is offline   Reply With Quote

Old   January 16, 2017, 02:15
Default
  #8
Senior Member
 
Ashwani
Join Date: Sep 2013
Location: Hyderabad
Posts: 154
Rep Power: 12
AshwaniAssam is on a distinguished road
Yes I am working on a wall model with RANS and I am trying to account for separation over ramp. Do you can provide any help.
AshwaniAssam is offline   Reply With Quote

Old   January 16, 2017, 17:13
Default
  #9
Member
 
Jingchang.Shi
Join Date: Aug 2012
Location: Hang Zhou, China
Posts: 78
Rep Power: 13
aerosjc is on a distinguished road
I'm working on a wall model with LES. And my final goal is to simulate flow over ramp. Yes, there will be separation. And I want to study the unsteadiness of shock wave/boundary layer interaction. So RANS does not work for my case, especially the unsteadiness of the separation. Why do you choose RANS to capture separation? Does that work?
aerosjc is offline   Reply With Quote

Old   January 17, 2017, 01:52
Default
  #10
Senior Member
 
Ashwani
Join Date: Sep 2013
Location: Hyderabad
Posts: 154
Rep Power: 12
AshwaniAssam is on a distinguished road
Yes, it works mostly link, I more interested for getting steady results. RANS would need to be modified to get something of sort URANS for unsteadiness (which is presently not in my scope). So, presently I am trying to make wall function work with our RANS code. LES is something we have decided in our group to look on later. DO, you have results for flow over ramp steady case using RANS?

Thanks
AshwaniAssam is offline   Reply With Quote

Old   January 17, 2017, 11:03
Default
  #11
Member
 
Jingchang.Shi
Join Date: Aug 2012
Location: Hang Zhou, China
Posts: 78
Rep Power: 13
aerosjc is on a distinguished road
Sorry, I've no results for your case with RANS.
aerosjc 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
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 06:04
Wall Functions pierre OpenFOAM Running, Solving & CFD 0 October 1, 2005 14:13
Wall functions - iterative approach agg Main CFD Forum 1 May 13, 2005 13:25
LES wall functions barrat Main CFD Forum 10 September 1, 2004 09:29


All times are GMT -4. The time now is 06:11.