|
[Sponsors] |
Wall function Spalart-Allmaras Turbulence Model |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
June 18, 2020, 10:46 |
Wall function Spalart-Allmaras Turbulence Model
|
#1 |
Senior Member
Join Date: Mar 2016
Location: Bergamo
Posts: 157
Rep Power: 10 |
Hello Foamers,
I'm going deeper to the implementaton of the Spalart-Alamars turbulence model and I have a big doubt about where wall function is called. In the kEpsilon and kOmega turbulence models the wall function is applied through Code:
epsilon_.boundaryFieldRef().updateCoeffs(); Code:
omega_.boundaryFieldRef().updateCoeffs(); On contrary, for he SA turbulence model there is no updateCoeffs() function. Any idea where wall functon is called? Best regards |
|
June 18, 2020, 17:02 |
|
#2 |
Senior Member
Herpes Free Engineer
Join Date: Sep 2019
Location: The Home Under The Ground with the Lost Boys
Posts: 932
Rep Power: 12 |
Hi,
- `nuTilda` BCs are updated by nuTilda_.correctBoundaryConditions(); - There is no `nuTilda` wall function in OpenFOAM - at least in the standard variants. - One of the derivation assumptions of the standard SA closure model is that y1+ is around 1. Therefore, by definition, wall functions are not required, and should not be used. - Having said that some researchers have been providing various approaches to modify the standard SA model, so that it can behave as a "high-Re" number type model to reduce computational costs, e.g. A generalised wall function including compressibility and pressure-gradient terms for the Spalart–Allmaras turbulence model. Hope it helps.
__________________
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 |
|
June 19, 2020, 05:04 |
|
#3 |
Senior Member
Join Date: Mar 2016
Location: Bergamo
Posts: 157
Rep Power: 10 |
Thank you for the reply. It makes a lot of sense. Anyway it rises a new question because the wall function is applied somewhere to nut and it also affects the value of nuTilda.
In the tutorial located in tutorial tutorials/incompressible/simpleFoam/airFoil2D the nuTilda file doesn't show any wall function (as you said) while the nut variable has a nutUSpaldingWallFunction wall function at the wall. I run two simulations with the nutUSpaldingWallFunction activated and not activated and the results for the nuTilda variable are shown in the attached files. So my question now is where the nutUSpaldingWallFunction is called. |
|
June 19, 2020, 06:45 |
|
#4 |
Senior Member
Herpes Free Engineer
Join Date: Sep 2019
Location: The Home Under The Ground with the Lost Boys
Posts: 932
Rep Power: 12 |
Hi,
- `nutUSpaldingWallFunction` has been called by this line: this->nut_.correctBoundaryConditions(); - As you know, `nut` and `nuTilda` are different quantities.
__________________
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 |
|
June 21, 2020, 04:47 |
|
#5 |
Senior Member
Join Date: Mar 2016
Location: Bergamo
Posts: 157
Rep Power: 10 |
Thank you, now it's more clear. I was just confused by the fact that the name of the function which updates the value at the wall for nut is different than the name of the functions that update the value of the k and omega.
In summary the wall function is applied to epsilon and omega through Code:
epsilon_.boundaryFieldRef().updateCoeffs(); Code:
omega_.boundaryFieldRef().updateCoeffs(); Code:
nut_.correctBoundaryConditions(); As last thing I wonder where the wall function is applied to the turbulent kinetic energy (k) for the models that calculate it. Best regards and thanks for the contribution |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Does Baldwin-Lomax turbulence model need wall function? | Abolfazl_cfd | Main CFD Forum | 0 | May 8, 2020 11:53 |
What's the problem with turbulence models near the wall region? | Jaydi_21 | Main CFD Forum | 6 | July 7, 2017 02:39 |
Ansys Fluent Turbulence Model wall function | sr71bb | ANSYS | 0 | August 13, 2015 15:57 |
Spalart Allmaras turbulence model | tareqkh | OpenFOAM Pre-Processing | 0 | December 9, 2014 21:30 |
Water subcooled boiling | Attesz | CFX | 7 | January 5, 2013 03:32 |