Near wall treatment in komega 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 komega SST model. Still a few questions remain, esp. concerning the implementation of the model in OF1.7.x, so hopefully someone can help me a bit out. Question in short form: Is the current komega SST implementation (in OF1.7.x) capable of automatically switching between lowRe and highRe model? Question in long form: If I understood it right, the komega SST model is a combination of the komega and kepsilon 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 highRe model. LowRe models don't use any wall functions as they are gridresolved. 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.58? 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 timevarying different y+ and u+ values. Therefore I need some kind of near wall treatment, automatically switching between lowRe and highRe formulation, depending on the location and flowsituation. So my question is: Is the current komega SST implementation (in OF1.7.x) capable of automatically switching between lowRe and highRe 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 lowRe model? I already ran some test cases using omegaWallFuntions, nutRoughWallFunction, kqRWallFunction and a y+ between 0.58. Nevertheless, the results seemed to be quite griddepending... Best, Arne [1] Menter: TwoEquation EddyViscosity 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 komegaSST 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 lowRe 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.cfdonline.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.cfdonline.com/Forums/ope...megasst3.html). Am I right? Thanks for your time and attention. Regards, Robert 
Wall functions should only be used for 30 < y+ < 100300, 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 
kOmega 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 komega 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/OpenFOAM2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #1 Foam::sigFpe::sigHandler(int) in "/home/cerecam/OpenFOAM/OpenFOAM2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #2 in "/lib/x86_64linuxgnu/libc.so.6" #3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/home/cerecam/OpenFOAM/OpenFOAM2.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/OpenFOAM2.1.1/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so" #5 at kOmegaSST.C:0 #6 Foam::incompressible::RASModels::kOmegaSST::F2() const in "/home/cerecam/OpenFOAM/OpenFOAM2.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/OpenFOAM2.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/OpenFOAM2.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/OpenFOAM2.1.1/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so" #10 in "/home/cerecam/OpenFOAM/OpenFOAM2.1.1/platforms/linux64GccDPOpt/bin/MRFSimpleFoam" #11 __libc_start_main in "/lib/x86_64linuxgnu/libc.so.6" #12 in "/home/cerecam/OpenFOAM/OpenFOAM2.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 kvalue to 1e11 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 kw SST and I uncheck LowReynoldsNumber 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?

All times are GMT 4. The time now is 22:51. 