Near wall treatment in k-omega SST
Hi all,
the last days I tried to get a bit deeper in the general understanding of turbulence modeling, and esp. in Menter's k-omega SST model. Still a few questions remain, esp. concerning the implementation of the model in OF-1.7.x, so hopefully someone can help me a bit out. Question in short form: Is the current k-omega SST implementation (in OF-1.7.x) capable of automatically switching between low-Re and high-Re model? Question in long form: If I understood it right, the k-omega SST model is a combination of the k-omega and k-epsilon model plus a shear stress transport model, using blending functions to choose the appropriate model depending on the region (boundary layer or free stream). So the basic implementation from [1] is a high-Re model. Low-Re models don't use any wall functions as they are grid-resolved. Some paper say that the viscous sublayer is between 0 < y+ < 5. Others say, that a good solution can only be achieved when y+ < 1, sometimes even < 0.1. So what is 'right' here, how big is the error if I'm in a range of y+ 0.5-8? I have a (interFoam) model using different grid resolutions (esp. first layer cell heights) in different regions at the walls. Additionally, different flow situations, i.e. different water velocities and also water waves are simulated, which lead to local- and time-varying different y+ and u+ values. Therefore I need some kind of near wall treatment, automatically switching between low-Re and high-Re formulation, depending on the location and flow-situation. So my question is: Is the current k-omega SST implementation (in OF-1.7.x) capable of automatically switching between low-Re and high-Re model? If this is hopefully true, do I then select omegaWallFunction? Looking at the source code of omegaWallFunction and comparing it to [2] I found that the blending function for y+ is a bit different. E.g. in the implementation, Code:
scalar omegaLog = sqrt(k[faceCellI])/(Cmu25*kappa_*y[faceI]) Code:
Cmu25 = pow(Cmu_, 0.25) Code:
omegaLog = 1/0.3*kappa *ut/y In the wall functions implementation, a wall roughness can be specified. Does this still work if switched to a low-Re model? I already ran some test cases using omegaWallFuntions, nutRoughWallFunction, kqRWallFunction and a y+ between 0.5-8. Nevertheless, the results seemed to be quite grid-depending... Best, Arne [1] Menter: Two-Equation Eddy-Viscosity Turbulence Models for Engineering Applications, 1994 [2] Menter: The SST Turbulence Model with Improved Wall Treatment for Heat Transfer Predictions in Gas Turbines, 2003 |
Hello Arne,
Nice thread you started here. I'm also interested in wall functions those days, and I took a closer loop to what it is actually implemented. Regarding the omegawallfunction first. What disturbs me is that the omega which is actually computed at the wall cell is a combination of omega_vis and omega_log. Omega=sqrt(omega_vis^2+omega_log^2). If you make a parallel with the epsilonwallfunction, with epsilon=c_mu*k*omega, you find that what is prescribed for epsilon is equivalent to what is prescribed for omega_log. I'm not yet very familiar with the komega-SST model, so maybe this combination of the viscous value and the log value is linked to equation of omega, I don’t know. What I would have except is that for y+<y_lam the code uses omega_vis, and for y+>y_lam, the code uses omega_log. It's what you call switching from low to high Reynolds. I'm simulating a train with a yaw angle. For low yaw angle the forces calculated match well the values found in wind tunnel, for higher yaw angles the forces are far too low. I suspect the wall functions. But again I don't know. Regarding the roughness, as far as I have understood, it is implemented via nut. You calculate a new value of nut=nu(y+*kappa/(ln(E_star+)-1). E_star is the only parameter modified by the roughness. Then I suppose that nut is used to calculate the value of the velocity in the wall cell using the wall shear stress tau_w/rho=u_star^2=(nu_t+nu)*du/dy which yields to u(y)=u_star*log(y+*E_star) And I think that is all for the roughness, means there is no direct connection I'm aware of between the values of omega (or epsilon) in the wall cell and the surface roughness. I'm also trying to compute atmospheric boundary layers, I couldn't' succeed in maintain the theorical profile in OpenFOAM adjusting the correct parameters of Ks and Cs, Fluent does. I think I'm missed something, but I don't know what. |
Hi Arne,
I found that post this morning, might be helpfull. http://www.openfoam.com/mantisbt/view.php?id=179 I quote Henry Quote:
|
Hi Sylvain,
thanks for giving me a feedback to my topic and sorry for not having answered yet. I'm quite busy with other things in OF right now, regarding implementing new things (but still at the very beginning...) I already got the information from somebody else and other threads that using komegaSST works for high- and low-Re grids (mixture) when using it together with nutSpalartAllmarasWallFunction and omegaWallFunction. If Henry says so as well, then it well be right... Concerning omega, this is what I got from another member (thanks for that!): "The menters version uses y+ for wall distance, but OF uses y* (compare with http://my.fit.edu/itresources/manual...ug/node512.htm). This seems to be standard procedure in OpenFOAM, for example if you use yPlusRAS you don't get y+ but y*". I will have a detailed look on this again and check it with some simple test when there is some more time... Greetings, Arne In the meanwhile, please let me know when you got new insights ;-) |
Hello Arne and Sylvain,
Not sure if this helps but there's a handy utility created by Florian, you may find it here http://www.cfd-online.com/Forums/ope...tml#post306730 On the same note, if I have a fine mesh at y* (or y+) less than 1 everywhere, I can opt not to use wall functions (according to Dr. Weller at post#54 in http://www.cfd-online.com/Forums/ope...megasst-3.html). Am I right? Thanks for your time and attention. Regards, Robert |
Wall functions should only be used for 30 < y+ < 100-300, which means that the first cell is in the log layer and not in the sublayer. If you want to use kOmega_SST without a wall function, y+ should be around 1.
|
Quote:
Arne |
k-Omega SST low Re considerations
hej,
as I have noticed now, the change in the "wall function" as suggested by henry is not to reliable. I found however, that when using k-omega SST in a low Re mode, one should have a very smooth transition of the cell size around the y+ < 200. Then the results become quite reliable. |
Quote:
|
not checked
hej,
i haven't checked at which time the model breaks down, but in the end I used 1.1 to 1.15 as my cell growth ratio. |
Low Re BC for K, Omega and nut
Dear all,
I am trying to simulate fan with low re number turb model, I have gone through so many post regarding wall treatment, still it is not clear for me. My yPlus vale is 3 to 6, I tried for k --> top0 { type fixedValue; value uniform 0; } top1 { type fixedValue; value uniform 0; } For Omega; top0 { type omegaWallFunction; value uniform 0; } top1 { type omegaWallFunction; value uniform 0; } and for nut; top0 { type calculated; value uniform 0; } top1 { type calculated; value uniform 0; } my questions are; is it right setting? for low re number turb model. when I try this settings, it is reporting the following error Create time Create mesh for time = 0 Reading field p Reading field U Reading/calculating face flux field phi Selecting incompressible transport model Newtonian Selecting RAS turbulence model kOmegaSST bounding k, min: 0 max: 0 average: 0 bounding omega, min: 0 max: 0 average: 0 #0 Foam::error::printStack(Foam::Ostream&) in "/home/cerecam/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #1 Foam::sigFpe::sigHandler(int) in "/home/cerecam/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #2 in "/lib/x86_64-linux-gnu/libc.so.6" #3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/home/cerecam/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #4 void Foam::divide<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/home/cerecam/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so" #5 at kOmegaSST.C:0 #6 Foam::incompressible::RASModels::kOmegaSST::F2() const in "/home/cerecam/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so" #7 Foam::incompressible::RASModels::kOmegaSST::kOmega SST(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&, Foam::word const&) in "/home/cerecam/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so" #8 Foam::incompressible::RASModel::adddictionaryConst ructorToTable<Foam::incompressible::RASModels::kOm egaSST>::New(Foam::GeometricField<Foam::Vector<dou ble>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&) in "/home/cerecam/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so" #9 Foam::incompressible::RASModel::New(Foam::Geometri cField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&) in "/home/cerecam/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so" #10 in "/home/cerecam/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/bin/MRFSimpleFoam" #11 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #12 in "/home/cerecam/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/bin/MRFSimpleFoam" Floating point exception (core dumped) please give me some idea what is wrong, Tanks for your time and help, Sivakumar |
You've got a floating point exception. Mostly it's caused by some division by 0.
Try setting your k-value to 1e-11 instead to zero. |
OMEGA on the wall surface in SST model
Dear all,
I am working with Fluent and I have simulated a flat plate and the mesh is enogh fine (Y+<=1). The turbulence model that I use is k-w SST and I uncheck Low-Reynolds-Number in the turbulence panel. I export the omega value on the plate surface and now I am trying to find what equation fluent uses for it. It would be appriciated if you help me about the equation. I used this equation but the answer was not the same with Fluent. OMG=sqrt(OMG_vis^2+OMG_log^2) OMG_vis=6*ro*(U_tau^2)/(0.075*Mu*Yplus^2) OMG_log=ro*(U_tau^2)/(0.3*k*Mu*Yplus) ro=dinsity U_tau=friction velosity Mu=viscosity k=0.41 |
Ali, this may be just a matter of cell and face value. Do you actually know, whether Fluent shows you the value of the wall or the value of the near wall cell?
Anyway, I think Fluent boundary condition of omega is not a boundary condition but setting the frist cell value to a fixed value. So there is actually no wall value of omega and if you plot the wall values this is some extrapolation of the (fixed) first cell value. I may be wrong with this. |
Thank you very much for replying.
To be honest, Fluent gives us value of the wall. My case study is a 2D flat plate and all values are the face one. To get the OMGa value of the wall, I did it in Display/plot/xy plot and then I selected Specific Dissipation Rate (Omega) and chose plate in the Surfaces. The Fluent gives us the Omega value of the plate surface. I checked the value of the cell adjecent (next cell after the wall boundary) and that value was not the same that I obtained with mentioned equation. In Fluent, we need a boundary condition for Omega near the wall to solve turbulence equation. Actually, Omega value of the wall surface is obtained by a formule (my problem is to find this formule in SST). Moreover, I do not think that the Omega value of the wall is extrapolation of the first cell value because Omega value has an increasing behaviour when closing to the wall and finally, it is logically accepted that dissipation should be sharply increased on the wall surface. Quote:
|
Let's assume you y+<1, then you can just calculate the OMG_vis value you mentioned above. OMG_log << OMG_vis.
What value do you get for OMG_vis and what does Fluent show you in the plot? |
3 Attachment(s)
When y+<1, then OMG_vis should give us the omega value of the wall. The problem is that the value of calculating the OMG_vis is not match to the Fluent one. I attach the Omega value of the wall, obtained by Fluent and equation OMG_vis. You can find that two results are not the same. It seems that Fluent uses different equation or another damping function or ....
Quote:
|
What is "x"?
|
"x" is a direction along the plate, x=0 is the leading edge of the plate and x=1.7 is the end of it.
Quote:
|
So your y+ changes along the plate?
|
2 Attachment(s)
Yes, it changes. I put all information in the excel file, enclosed here
Quote:
|
This looks strange. Is "wall_dist" the first cell distance in "m"? Is this exported from Fluent or from your meshing tool?
I am asking because, did you write down the OMG_vis formula without wall units? It's OMG_vis = 6 * nu / (0.072 y^2) So for constant y (in m) you should have a constant omega wall value... |
Yes, "wall_dist" is the adjcent cell to the wall in "m" and it is exported from Fluent.
The OMG_vis formula is based on y+ not y and y+ is a nondimensional parameter. Moreover, Omega depents on wall shear, too; so that it has not constant value along the plate. The main problem is that what formula Fluent uses for the Omega value of the wall. Quote:
|
Quote:
You will get OMG_vis = 6 * nu / (0.072 y^2) No y+, no wall shear. |
Yes, you are right. Then, how is Fluent calculate the omega of the wall?
Actually, the same formula is written in Fluent theory guide, Quote:
|
Are you 100% sure about the grid? Is the wall distance really constant?
|
2 Attachment(s)
Yes, I am sure, because my case study is a flat plate and I use hexahedral elements (structural mesh) and cells have the same wall distance.
I enclosed two figures of the mesh Quote:
|
I remember when I checked these values for one of my cases everything was alright. There must be some mistake in the data. Can you upload the whole fluent case?
|
Yes, sure, you can find the cas and data file in the Dropbox, linked here:
https://www.dropbox.com/sh/rg87mkr95...MyjuKLP-a?dl=0 Quote:
|
Ali, you should try to understand my post #14 in this thread (very first post of me). I think you can not plot omega on the surface, because Fluent doesn't let you. You are currently plotting interpolated values. Try this: plot "k" along the plate. You will see, that it isn't zero everywhere. This doesn't make sense, because k-boundary condition is zero. This can only be explained by interpolation errors during post-processing.
I also found in the userguide: In ANSYS FLUENT you can choose to display the computed cell-center values or values that have been interpolated to the nodes. By default, the Node Values option is turned on, and the interpolated values are displayed. For line contours or profiles, node values are always used. If you are displaying filled contours or profiles and you prefer to display the cell values, turn the Node Values option off. Filled contours/profiles of node values will show a smooth gradation of color, while filled contours/profiles of cell values may show sharp changes in color from one cell to the next. For face-only functions (e.g., Wall Shear Stress), the cell values that are displayed for boundary zone surfaces will actually be the face values. This is only true in the case of boundary zone surfaces created for postprocessing, where the actual cell values are used for the part of the surface which lies in the interior. These face values are more accurate, as face-only functions are computed on the faces and not on the cells. |
Hi!
I don't understand this: Quote:
Code:
/*---------------------------------------------------------------------------*\ |
yPlus value is calculated for each cell next to the wall, each cell has it's y+ value, that's why you have min/max and avg values.
|
I know..
But they are not in the supposed interval, that's because I posted that |
Quote:
|
Infinite thanks, Roman! :)
If I use the same wall functions of the motorbike tutorial (nutKwallfunction, kqRwallfunction and omegawallfunction) and standard k omega model what values of y+ I should have? And for standard k epsilon with epsilonwallfunction? |
Quote:
|
Quote:
More information on wall functions and near wall behavior can be found in 1. Wilcox, D. C. (2006). Turbulence modeling for CFD (3rd ed.). La Canada Calif.: DCW Industries. 2. Pope, S. B. (2000). Turbulent flows. Cambridge [u.a.]: Cambridge Univ. Press. 3. Schlichting, H., & Gersten, K. (2000). Boundary Layer Theory (8th Revise). Darmstadt: Springer Berlin Heidelberg. 4. Bredberg, J. (2000). On the Wall Boundary Condition for Turbulence Models. Göteborg, Sweden. (gives explanation of wall functions in OpenFOAM) 5. Schiestel, R. (2008). Modeling and Simulation of Turbulent Flows. London, UK: ISTE. http://doi.org/10.1002/9780470610848 |
Quote:
It's the first time I hear about y+ values below 0.2 leading to wrong results. Do you have more information about this? Best regards, Kate |
Helllo,
So by having all in hand, now in modelling high-speed flow with K-w sst, if your first grid is in Y+=150 (log layer). and you are interested in studying the vorticity of separated flow downstream of a backstep (Backward Facing Step), does this putting first grid point on log layer make any problem? :confused: what is the difference between Fluent and OpenFoam in this? Can anyone introduce me some papers in this regard? Thanks |
All times are GMT -4. The time now is 22:33. |