|
[Sponsors] |
what is correct boundary condition for K-w sst model at y+<1? |
![]() |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
![]() |
![]() |
#1 |
Senior Member
Join Date: Jun 2011
Posts: 163
Rep Power: 15 ![]() |
Hi
I found some thread about k-w sst boundary conditions unfortunately they are not clear in low y+ please let me know what is correct boundary conditions for k, w and nut for y+<5 Best Regards |
|
![]() |
![]() |
![]() |
![]() |
#2 |
New Member
Andre Weiner
Join Date: Aug 2012
Posts: 29
Rep Power: 14 ![]() |
Hello!
omega: Menter's automatic wall funcion works for all parts of the momentum boundary layer, in openfoam it's the omegaWallFunction http://www.cfd-online.com/Wiki/Near-...k-omega_models or you use the solution for the viscos sublayer to set the boundary value http://turbmodels.larc.nasa.gov/sst.html (there is an equation for omega_wall) for y^+<5 both give the same results k: fixed value 0 or zero gradient if you look at same results for k near the wall (from other simulations for example), you will see that both are valid (again for y^+<5) nut: calculated or fixed value 0 or nutLowReWallFunction (sets the value to zero and provides a function for calculating y^+ which is needed for the yPlusRAS tool) This is what i tried so far and it worked for mi :-) Best regards, Andre |
|
![]() |
![]() |
![]() |
![]() |
#3 | |
Senior Member
Join Date: Jun 2011
Posts: 163
Rep Power: 15 ![]() |
Quote:
what function should be set for calculating y+ ? and how set it ? can you attached an example for k,w,nut (initial file at folder 0)? Best Regards |
||
![]() |
![]() |
![]() |
![]() |
#4 |
Senior Member
Joachim
Join Date: Mar 2012
Location: Paris, France
Posts: 145
Rep Power: 15 ![]() |
Actually, I think that nutLowReWallFunction is a wall function for low-Re turbulence models. If you use LaunderSharmaKE and have y+ ~ 60, you can use it for example.
Also, try using 1e-9 (or any small value) rather than 0 when imposing the BC. ![]() Good luck! |
|
![]() |
![]() |
![]() |
![]() |
#5 |
New Member
Andre Weiner
Join Date: Aug 2012
Posts: 29
Rep Power: 14 ![]() |
for k:
Code:
myWall { type zeroGradient; } // or myWall { type fixedValue; value uniform 0; } Code:
myWall { type omegaWallFunction; } Code:
myWall { type fixedValue; value uniform 0; } // or myWall { type calculated; value uniform 0; // optional value entry } // what i prefer myPatch { type nutLowReWallFunction; } it sets nut to zero: Code:
tmp<scalarField> nutLowReWallFunctionFvPatchScalarField::calcNut() const { return tmp<scalarField>(new scalarField(patch().size(), 0.0)); } Code:
tmp<scalarField> nutLowReWallFunctionFvPatchScalarField::yPlus() const { .... return y*sqrt(nuw*mag(Uw.snGrad()))/nuw; } |
|
![]() |
![]() |
![]() |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
Epsilon boundary condition for walls in LamBremhorst LowRe model | maruthamuthu_venkatraman | OpenFOAM | 1 | April 21, 2019 06:25 |
Error finding variable "THERMX" | sunilpatil | CFX | 8 | April 26, 2013 07:00 |
CFX13 Post Periodic interface | EtaEta | CFX | 7 | December 8, 2011 17:15 |
SST model and boundary layers | Rodrigo Basniak | CFX | 7 | March 27, 2008 14:19 |