|
[Sponsors] |
in k-epsilon wall function approach high Re turbulence models: question of velocity |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
June 27, 2013, 05:04 |
in k-epsilon wall function approach high Re turbulence models: question of velocity
|
#1 |
Senior Member
Roman Thiele
Join Date: Aug 2009
Location: Eindhoven, NL
Posts: 374
Rep Power: 21 |
Hej,
I looked through the different wall function for k-epsilon turbulence models, however I feel I am missing something. In Wilcox (Turbulence modeling for CFD 2006) as well as many other publications, in wall functions the velocity seems to be set in the first grid point. In OpenFOAM, the following properties are changed with wall functions in the example of the k-epsilon turbulence model,
__________________
~roman |
|
July 2, 2013, 06:18 |
took a while but here is a solution to the problem
|
#2 |
Senior Member
Roman Thiele
Join Date: Aug 2009
Location: Eindhoven, NL
Posts: 374
Rep Power: 21 |
Hej,
I think I found the solution to this. I found it in Pope, Turbulent Flow, 2000 under the way that wall functions are created. The original question was why the velocity is not set by any of the wall function. The answer is, that wall functions are apparently more robust by setting the shear stress and not the velocity, which makes sense when looking at separation points and other nasty flow phenomena, in which the velocity can become essentially 0. The way this is for the k-epsilon model done is nutkWallFunction
__________________
~roman |
|
July 2, 2013, 14:22 |
|
#3 |
Senior Member
Mohsen KiaMansouri
Join Date: Jan 2010
Location: CFD Lab
Posts: 118
Rep Power: 16 |
Thanks roman
it was very helpful.
__________________
“If you have an apple and I have an apple and we exchange these apples then you and I will still each have one apple. But if you have an idea and I have an idea and we exchange these ideas, then each of us will have two ideas.” |
|
July 2, 2013, 20:02 |
|
#4 |
Senior Member
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,208
Rep Power: 27 |
hi
Whould you explain more and Why stage 4 is done?
__________________
Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King. To Be or Not To Be,Thats the Question! The Only Stupid Question Is the One that Goes Unasked. |
|
July 3, 2013, 05:26 |
|
#5 |
Senior Member
Roman Thiele
Join Date: Aug 2009
Location: Eindhoven, NL
Posts: 374
Rep Power: 21 |
Step 4 is introduced in order to have a more stable solver. Instead of setting the velocity in the first point away from the wall, a force (the wall shear force) is introduced in the equation. This way, the solver won't have any 0 velocities introduced by accident (think of stagnation points), this could destabilize the solver as you would be introducing 0's into the matrix.
The easiest way to introduce this source term into the equation is by modifying the turbulent viscosity, see Pope, Turbulent Flows, 2000, pp 442-443 and Versteeg and Malalasekera, An introduction to computational fluid dynamics: The finite volume method, pp. 276-277
__________________
~roman |
|
March 15, 2015, 22:12 |
|
#6 | |
Senior Member
Huang Xianbei
Join Date: Sep 2013
Location: Yangzhou,China
Posts: 302
Rep Power: 14 |
Quote:
Thank you for your explain. I have a question about the wall function for non-linear models. As the non-linear model can't be depicted by a simple turbulent viscosity, how can we set the shear stress? say we obtained the shear stress as in your form: u_shear * U_p,real / U_p,nominal, so we should set the Reynolds shear stress by u_shear * U_p,real / U_p,nominal - nu*dU/dy, here, y represents the wall-normal direction. However, the non-linear Reynolds stress is not symmetric, it has 6 shear stress components, which one should be modified? I'm confused by this problem for some days. I'm looking forward for your reply. Thank you Xianbei |
||
May 26, 2016, 10:14 |
|
#7 | |
New Member
karar
Join Date: Feb 2013
Posts: 13
Rep Power: 13 |
Quote:
Did you solve your wall function problem. I have tried to implement a non-linear model and I have stuck the same problem. What kind of wall function should I use. thank you |
||
Tags |
turbulence model, velocity, wall function |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[Other] mesh airfoil NACA0012 | anand_30 | OpenFOAM Meshing & Mesh Conversion | 13 | March 7, 2022 18:22 |
Discussion: Reason of Turbulence!! | Wen Long | Main CFD Forum | 3 | May 15, 2009 10:52 |
Variables Definition in CFX Solver 5.6 | R P | CFX | 2 | October 26, 2004 03:13 |
Quick Question - Wall Function | D.Tandra | Main CFD Forum | 2 | March 16, 2004 05:29 |