|
[Sponsors] |
February 13, 2021, 17:11 |
uPlus vs yPlus ?
|
#1 |
New Member
Kermit
Join Date: Jan 2021
Location: France
Posts: 19
Rep Power: 5 |
Hello everyone
I'm currently computing the turbulent flow over a flate plate in order to prove that wall functions are useful to reproduce the theoretical laws of the wall. My simulation settings are quickly summed up below: - solver: simpleFoam. - boundary layer mesh spec: around 120 cells in 0.1m with refinement near the wall. - turbulence model: k-epsilon. - wall functions: kqRWallFunction, epsilonWallFunction, nutkWallFunction. Here is how I evaluate y+ and u+: 1) First, I get the value of the wallShearStress t_w (simpleFoam -postProcessing -func wallShearStress): t_w = sqrt(t_w(x)² + t_w(y)² + t_w(z)²) 2) Then, I evaluate the friction velocity value: u_tau = sqrt(t_w) NB: Units of wallShearStress values in openFoam are m²/s² instead of kg.m^-1/s². So we don't need do divide by rho (density). 3) Third, I evaluate y+ and u+ for each y: y+ = y * u_tau / nu, where nu is the kinematic viscosity. u+ = u(y) / u_tau, where u(y) is the x-component of the velocity field Is there anything wrong in my methodology ? Someone know how yPlus is calculated by openFoam with the utility: simpleFoam -postProcessing -func yPlus ? I hope someone will be able to help me. Thanks in advance! |
|
February 13, 2021, 17:17 |
|
#2 |
Senior Member
Herpes Free Engineer
Join Date: Sep 2019
Location: The Home Under The Ground with the Lost Boys
Posts: 931
Rep Power: 13 |
The line of computations is OK to my knowledge. Yet didn't quite understand the main question.
Two points: - For the zero-gradient flat plate case, consider t_w(x) only - and ignore the other two components. - Make sure using mag(t_w(x)) to avoid the negative sign - otherwise sqrt will throw a domain error.
__________________
The OpenFOAM community is the biggest contributor to OpenFOAM: User guide/Wiki-1/Wiki-2/Code guide/Code Wiki/Journal Nilsson/Guerrero/Holzinger/Holzmann/Nagy/Santos/Nozaki/Jasak/Primer Governance Bugs/Features: OpenFOAM (ESI-OpenCFD-Trademark) Bugs/Features: FOAM-Extend (Wikki-FSB) Bugs: OpenFOAM.org How to create a MWE New: Forkable OpenFOAM mirror |
|
February 13, 2021, 17:38 |
|
#3 |
New Member
Kermit
Join Date: Jan 2021
Location: France
Posts: 19
Rep Power: 5 |
Thanks for your quick reply and advice !
I just dont understand why my y+ and u+ dont match with theoretical values (Spalding law). Is that because my mesh is too refined near the wall (openfoam gives me a yPlus of 0.32) for a kEpsilon model ? Otherwise, I noticed that the pressure isnt really converged (residual=10^-3). NB: My case geometry is based on the following one: https://www.openfoam.com/documentati...plate-zpg.html |
|
February 13, 2021, 18:10 |
|
#4 |
Senior Member
Herpes Free Engineer
Join Date: Sep 2019
Location: The Home Under The Ground with the Lost Boys
Posts: 931
Rep Power: 13 |
`kEpsilon` is tricky for y+<30 meshes since the derivation assumption of this turbulence closure model is y+>30 (i.e. it cannot model below the inertial sublayer). Having said that you are using wall functions - which should allow you to model below y+ 30, so I presume the issue is relevant to the settings used for `kEpsilon` or `k`-`epsilon` fields.
Note that in the link you have shared, the turbulence model being used is `kOmegaSST`. See if the ZPG flat plate tutorial helps you: https://develop.openfoam.com/Develop...ulentFlatPlate where you can find some settings for `kEpsilon`. This tutorial is more or less the replica of the ZPG case available in the NASA turbulence website. Hope these help.
__________________
The OpenFOAM community is the biggest contributor to OpenFOAM: User guide/Wiki-1/Wiki-2/Code guide/Code Wiki/Journal Nilsson/Guerrero/Holzinger/Holzmann/Nagy/Santos/Nozaki/Jasak/Primer Governance Bugs/Features: OpenFOAM (ESI-OpenCFD-Trademark) Bugs/Features: FOAM-Extend (Wikki-FSB) Bugs: OpenFOAM.org How to create a MWE New: Forkable OpenFOAM mirror |
|
February 14, 2021, 06:01 |
|
#5 |
New Member
Kermit
Join Date: Jan 2021
Location: France
Posts: 19
Rep Power: 5 |
Thanks, you help a lot !
I computed the case from the link you shared. It seems that with kOmegaSST it's easier to reproduce the log-law and I'm not surprised with that results because this model is way better for modelling the boundary layer. However, for the kEpsilon model it still not good... It looks like wall functions aren't working well in my case. |
|
Tags |
simplefoam, uplus, wallfunction, yplus |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Uplus vs Yplus with K-Epsilon Turbulent Model in Fluent | sandpants | FLUENT | 2 | December 17, 2022 13:08 |
[Other] Contribution a new utility: refine wall layer mesh based on yPlus field | lakeat | OpenFOAM Community Contributions | 58 | December 23, 2021 03:36 |
Uplus vs Yplus Non-Newtonian Fluids | C.C | FLUENT | 1 | November 19, 2012 13:09 |
Uplus vs Yplus | C.C | FLUENT | 7 | November 19, 2012 10:00 |
Compressible yPlus and uPlus for LES | owayz | OpenFOAM Running, Solving & CFD | 0 | August 13, 2012 20:03 |