|
[Sponsors] |
September 29, 2016, 09:56 |
One more time: wall functions - SST
|
#1 | ||||
Member
Join Date: Apr 2016
Posts: 91
Rep Power: 10 |
Hi everyone,
I struggle choosing the right wall functions. I am using SST in my simulation, since it combines the advantages of k-epsilon and k-omega. As far as I searched for wall functions and SST in this forum I found: http://www.cfd-online.com/Forums/ope...ion-usage.html Quote:
Quote:
Quote:
I am wondering now: one source says SST uses wall functions, one says: no . So one is wrong or do I confuse anything? The range of my y+ values are for a really fine grid on some walls 0.01 and on other walls y+=45. So when I use nutSpalartAllmarasWallFunction , what range should my y+ be? Member vkrastev says: Quote:
Could someone bring light into the dark (dark for me). By the way: currently I am using mutkWallFunction which makes trouble for very fine meshes with SST. So that's why I have doubts about it. Sorry I am a rookie in CFD :-/ Thank you guys! |
|||||
October 5, 2016, 10:58 |
|
#2 | |
Senior Member
Join Date: Mar 2015
Posts: 250
Rep Power: 12 |
Quote:
Best regards, Kate Last edited by KateEisenhower; October 5, 2016 at 11:02. Reason: typo |
||
October 7, 2016, 08:52 |
|
#3 |
Member
Join Date: Apr 2016
Posts: 91
Rep Power: 10 |
I think yes. If the grid is too fine, the wall laws can lead to wrong results. Does anyone have some suggestions?
|
|
February 3, 2017, 07:07 |
|
#4 |
New Member
Martin
Join Date: Oct 2013
Location: Newcastle
Posts: 21
Rep Power: 13 |
Hi There,
did you manage to clarify the above? I'm looking for some more information regarding yplus levels in combination with KomegaSST and the nutUSpaldingWallFunction. I'm running simulations at the moment with yplus ranging from <1 to <100. If someone can clarify or suggest limits on yplus in this approach would be very helpful as i'm trying to work with automatically generated meshes in Openfoam which makes mesh manipulation quite hard to find good yplus. Thanks |
|
February 3, 2017, 11:34 |
|
#5 |
Senior Member
Join Date: Mar 2014
Posts: 112
Rep Power: 12 |
I have done a test to understand the same issue about 2 years ago with 2.3.1.
The set up was a 2D external case with naca0012 profile. The kwSST model in combination with nutUSpladingWallFunction run on 11 different meshes (structured O-grid / Re 1,8e06) which have following average yPlus values: 1, 5, 15, 30, 50, 70, 100, 150, 300, 500 and 1000. The first two have diverged and/or produced unphysical results and all others worked and produced acceptable results. If I remember correctly average yPlus 30 was the best against the wind tunnel data. The same setup worked well with SA model on all mesh range since the SA implementation on OpenFOAM can work on low and high Re modes but kwSST can not. There are also some lowRe implementations of kwSST for OpenFOAM around though and I guess you can find it in this forum... |
|
February 3, 2017, 13:34 |
|
#6 |
Senior Member
Mikko
Join Date: Jul 2014
Location: The Hague, The Netherlands
Posts: 243
Rep Power: 13 |
I have simulated NACA0012 with recent foundation and foam-extend versions and my results are close to NASA's. I used k-omega SST turbulence model without wall functions (y+<1) and following boundary condition for the wall:
nut=0, k=0 and omega=omegaWallFunction. |
|
February 8, 2017, 07:29 |
|
#7 | |
Member
jey
Join Date: Nov 2016
Location: Greece
Posts: 30
Rep Power: 9 |
Quote:
Dear Flowkersma I can understand that since y+~1, you don't have to use wall functions (as you already did for nut and k). But why did you still use wall functions for omega? Is there a specific reason for that? Thanks in advance for your help! Kind regards |
||
February 9, 2017, 09:00 |
|
#9 |
Senior Member
Join Date: Mar 2014
Posts: 112
Rep Power: 12 |
I was not aware that the standart kOmegaSST model in OpneFOAM is capable to work on on low-Re grids and had given up after the failed results with the continious wall function (nutUSplading) on this model as decribed my above post.
I have just tested the standart kOmegaSST without wall functions on y+ ~1 meshed cases and got good results with that approcach so thanks for the information to Flowkersma. |
|
February 9, 2017, 09:50 |
|
#10 |
Member
jey
Join Date: Nov 2016
Location: Greece
Posts: 30
Rep Power: 9 |
Dear Flowkersma and mzzmrt
I am using the conditions that Flowkersma suggested (post#6) for a y+~1 case, but I can not make it work (it crashes after some time) But.. I tried using a fixedValue for omega at wall (as suggested by Menter - see attached picture) and the simulations works fine until now. omegaAtWall.PNG I am really frustrated that I can not make it work with omegaWallFunction while you people can do it. I am probably missing something? Please find here my /0 and /system folder, in case you can help with any suggestions. p.s.1: The value I am calculating at y=0 according to Menter, is of the same order of magnitude as the one that omegaWallFunction sets at wall. p.s.2: This is the way I implement the omegaWallFunction type. I noticed that the function is completely independent of the "value" entry. So, the "value" entry is just a placeholder. Is it correct? Code:
<patchName> { type omegaWallFunction; value $internalField; } Thank you in advance! Kind regards |
|
February 9, 2017, 17:42 |
|
#11 |
Senior Member
Mikko
Join Date: Jul 2014
Location: The Hague, The Netherlands
Posts: 243
Rep Power: 13 |
The value field is the initial condition. Generally, a good initial value results in faster convergence. Bad initial value may result in divergence. In this case it probably does not play a big/any role.
I use first order upwind scheme for divergence for the turbulence quantities. Give it a try. |
|
February 14, 2017, 03:18 |
|
#12 |
Super Moderator
Tobias Holzmann
Join Date: Oct 2010
Location: Bad Wörishofen
Posts: 2,711
Blog Entries: 6
Rep Power: 51 |
Dear everybody,
I just want to give you some information about the usage of the wallFunctions. It might be, that I am not correct in all details but as I got it out of a few literatures it is like the following; I refer to [1], given at my webpage. The flow close to a wall has a viscose sublayer which depends on different parameters like the roughness and the velocity components. You can imagine -- and it is easy to understand -- that with higher velocity (or Re-Number), the viscose sublayer decrease. It means that the layer thickness gets smaller. Okay, so far so good, right? (Check). Now, there occur several problems which I will not talk about in detail (I think turbulence guys are much more advanced and experienced in that topic and they might say that we do everything wrong). You should keep in mind that the standard models that we use are based on isotropy which is not really valid in a lot of cases but however, ... The main disadvantage now is related to these viscose sublayer because it is obvious that the damping of the velocity normal to the wall is much higher than tangential to the wall, which makes the isotropic assumption even worse (check). Okay, now you can imagine that if you have high Re flows, the sublayer decreases extremely. In order to resolve the velocity change in the sublayer you have two options:
Keep in mind that the sublayer is important for mass / heat-transfer at the boundaries and will influence the result dramatically. However, we also should keep in mind that RANS is an averaging method. Using compressible -- density based -- flows it's getting even worse because of the FAVRE averaging method which is only correct in a mathematical point of view. In addition, if you ever derived the Reynolds-Stress equation and out of that the turbulence equations like k, epsilon etc. you will see that there are so many terms we just model and it is just an approximation which is good or bad based on the case we are looking at etc. In addition: We will never get a mesh where the y+ is always correct. Just think about flow separation. At last, I am not an expert and the stuff above might be wrong or not complete. It is just my opinion about the problem based on the literature I read. As I already stated above, the people who really deal with turbulence flow might tell us that we are doing everything wrong :P CFD in RANS is just modeling - keep that in mind. If you want to see the modeled terms in the k and epsilon equations, check out my book.
__________________
Keep foaming, Tobias Holzmann |
|
December 16, 2017, 06:39 |
no wall function
|
#13 | |
Senior Member
Guilherme
Join Date: Apr 2017
Posts: 232
Rep Power: 10 |
Quote:
|
||
December 17, 2017, 06:01 |
|
#14 |
Super Moderator
Tobias Holzmann
Join Date: Oct 2010
Location: Bad Wörishofen
Posts: 2,711
Blog Entries: 6
Rep Power: 51 |
Hi,
actually I am not familiar with turbulence modeling. As far as I understand the problem --- which should be stated in my last post --- in a numerical/mathematical way, if you have y+~1, then you should not apply wall functions (sure, the subsequent cells should resolve the viscous sub layer). Therefore, I do not know which value one has to set. ZeroGradient? I read somewhere that you fix k to very small value such as 0.001. But please don't ask me about that. There are much more advanced users here who should give a fundamental and clear statement.
__________________
Keep foaming, Tobias Holzmann |
|
April 19, 2018, 17:42 |
|
#15 |
Senior Member
Join Date: Mar 2018
Posts: 115
Rep Power: 8 |
Hello Tobias
How can I disable wall functions in OpenFOAM if I have y+ << 1 ? |
|
June 29, 2018, 00:00 |
|
#16 |
New Member
Juan Felipe Monsalvo Salazar
Join Date: Apr 2013
Location: Medellín / Colombia
Posts: 4
Rep Power: 13 |
Dear Tobias
Just to be sure, reviewing the following link https://en.wikipedia.org/wiki/Law_of_the_wall I saw that the log law is only valid in the range between [ 30 < y+ < 200 or 300] and not at [ 0 < y+ < 16.8] as you state. Neverthless I'm not sure if you was referring to the wall of law when talking about log law Sincerely, Juan Felipe Monsalvo |
|
June 29, 2018, 16:44 |
|
#17 | |
New Member
Juan Felipe Monsalvo Salazar
Join Date: Apr 2013
Location: Medellín / Colombia
Posts: 4
Rep Power: 13 |
Quote:
Just to be sure, reviewing the following link https://en.wikipedia.org/wiki/Law_of_the_wall I saw that the log law is only valid in the range between [ 30 < y+ < 200 or 300] and not at [ 0 < y+ < 16.8] as you state. Neverthless I'm not sure if you was referring to the wall of law when talking about log law Sincerely, Juan Felipe Monsalvo |
||
September 21, 2020, 07:25 |
|
#19 |
Super Moderator
Tobias Holzmann
Join Date: Oct 2010
Location: Bad Wörishofen
Posts: 2,711
Blog Entries: 6
Rep Power: 51 |
Nice summary but a small review would increase its qulaity significantly in terms of readablity, typo and others.
__________________
Keep foaming, Tobias Holzmann |
|
March 30, 2021, 19:53 |
|
#20 |
Senior Member
Join Date: Jun 2011
Posts: 206
Rep Power: 16 |
also it would be very helpful (and practical) if one knows and lists the Y+ values of the grid needed (i.e. providing the highest theoretical accuracy) when using different turbulent models (the most popular ones).
The information in the links above is too difficult to digest to a simple conclusion especially by not CFD professionals. As Tobias said a short summary on this important subject would be very helpful to a lot of users. |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Natural convection in a closed domain STILL NEEDING help! | Yr0gErG | FLUENT | 4 | December 2, 2019 00:04 |
Hardware-Configuration for Fluent HPC-Pack (8x) | JohHaas | Hardware | 9 | March 3, 2015 13:25 |
Low Mixing time Problem | Mavier | CFX | 5 | April 29, 2013 00:00 |
Boundary Field and Wall Functions, time 0 & time >0 | NJG | OpenFOAM Running, Solving & CFD | 1 | April 5, 2013 10:48 |
Multicomponent fluid | Andrea | CFX | 2 | October 11, 2004 05:12 |