[{"data":1,"prerenderedAt":161},["ShallowReactive",2],{"$f3c2i3qdvokjkk":3},{"slug":4,"publishedAt":5,"category":6,"title":7,"summary":8,"image":9,"imageAlt":14,"readingMinutes":15,"serviceId":16,"keywords":17,"imageCaption":22,"sections":23,"sources":151},"k-omega-sst-vs-k-epsilon","2026-10-12","methods","k-omega SST vs k-epsilon: Which Model and What y+","Choosing between k-omega SST and k-epsilon depends on how you treat the wall. Compare the two models, see which y+ each needs, and calculate the first cell height.",{"src":10,"thumbnail":11,"width":12,"height":13},"\u002Fassets\u002Fcfd-cases\u002Faero\u002Fimages\u002F01.webp","\u002Fbrand\u002Fstudies\u002Faero.svg",1672,941,"DEEB-FLOW illustration of velocity streamlines in an aerodynamics case",7,"flow",[18,19,20,21],"k-omega SST vs k-epsilon","y+ wall functions","first cell height calculation","turbulence model selection","Concept illustration from the DEEB-FLOW case library. It is not a result from this article.",[24,30,41,75,91,128,135,146],{"id":25,"title":26,"paragraphs":27},"choice-about-the-wall","The model choice is a choice about the wall",[28,29],"Ask a group of CFD engineers which turbulence model to use and most will say k-ω SST. Ask what y+ their last mesh had at the wall and the answers come more slowly.","The second question decides whether the first answer was any good. Choosing between k-ω SST and k-ε is largely a choice about how you treat the wall, and a good model on the wrong near-wall mesh gives a worse answer than a simpler model used properly. This article explains what each model does well, which y+ each needs, and how to calculate the first cell height before you mesh.",{"id":31,"title":32,"paragraphs":33},"two-models","What the two models are",[34,35,40],"Both are two-equation models. Both add a turbulent viscosity to the flow equations and calculate it from two transported quantities. They differ in the second quantity.",{"list":36},[37,38,39],"**k-ε** transports the turbulent kinetic energy k and its dissipation rate ε. It dates from the early 1970s. It is well behaved in the bulk of a flow, insensitive to what you assume about turbulence at the inlet, and cheap. Its weakness is the wall. The standard form cannot be integrated through the viscous sublayer, so it relies on wall functions to bridge the gap. It also tends to predict separation too late and too small under adverse pressure gradients, and it overproduces turbulence at stagnation points.","**k-ω** transports k and the specific dissipation rate ω. It can be integrated down to the wall and handles adverse pressure gradients better. Its weakness is the opposite one. In its original form, the solution depends on the free-stream value of ω you happen to set.","**k-ω SST**, published by Menter in 1994, uses k-ω near the wall and blends to a k-ε formulation away from it, taking the stronger half of each. It also limits the turbulent shear stress in adverse pressure gradients, which improves its prediction of separation.","That history explains why SST became the default, and a default still has to be checked against the case in front of you.",{"id":42,"title":43,"paragraphs":44},"what-is-y-plus","What is y+ and why does it matter?",[45,46,49,50,69,70,74],"y+ is a dimensionless distance from the wall.",{"formula":47},[48],"y+ = y × u_τ \u002F ν","Here y is the distance from the wall to the centre of the first cell, u_τ is the friction velocity, which is the square root of the wall shear stress divided by the density, and ν is the kinematic viscosity. The value tells you which part of the boundary layer your first cell sits in.",{"table":51},{"head":52,"rows":56},[53,54,55],"Region","y+","What happens there",[57,61,65],[58,59,60],"Viscous sublayer","y+ \u003C 5","Viscosity dominates and velocity is linear in distance",[62,63,64],"Buffer layer","5 \u003C y+ \u003C 30","Neither description fits well",[66,67,68],"Log-law region","30 \u003C y+ \u003C 300","Turbulence dominates and velocity follows the logarithmic law","That leaves two strategies.",{"list":71},[72,73],"**Resolve the boundary layer.** Put the first cell in the viscous sublayer, at y+ of about 1, with a dozen or more cells across the boundary layer. The model is integrated to the wall. This is what k-ω SST was built for.","**Use wall functions.** Put the first cell in the log-law region, at y+ of 30 or above, and let an algebraic formula supply the wall shear stress. This is how standard k-ε is meant to be used.","The place to avoid, with standard wall functions, is in between. A first cell at y+ of 10 sits in the buffer layer, where the assumptions of the formula do not hold. Many current solvers offer blended or automatic wall treatments that reduce this penalty without removing it.",{"id":76,"title":77,"paragraphs":78},"first-cell-height","Working out the first cell height",[79,80,89,90],"Do this before meshing. It takes five minutes and saves a remesh. As a worked example, take air at 20 m\u002Fs over a surface 1 m long, with a density of 1.2 kg\u002Fm³ and a kinematic viscosity of 1.5 × 10⁻⁵ m²\u002Fs. The skin friction comes from a flat-plate correlation.",{"formula":81},[82,83,84,85,86,87,88],"Re = U × L \u002F ν = 20 × 1 \u002F (1.5 × 10⁻⁵) = 1.33 × 10⁶","C_f ≈ 0.058 × Re^(-0.2) = 0.00346","τ_w = ½ × ρ × U² × C_f = 0.5 × 1.2 × 400 × 0.00346 = 0.83 Pa","u_τ = √(τ_w \u002F ρ) = √(0.83 \u002F 1.2) = 0.83 m\u002Fs","y = y+ × ν \u002F u_τ","y+ = 1:  y = 1.8 × 10⁻⁵ m = 0.018 mm","y+ = 30:  y = 5.4 × 10⁻⁴ m = 0.54 mm","That y is the distance to the cell centre. In a cell-centred solver the first cell is twice as tall, about 0.036 mm to resolve the wall and about 1.1 mm for wall functions. A factor of thirty in first cell height is a large difference in cell count and in the aspect ratio of those cells. It is the real cost of the choice.","Two cautions apply. The flat-plate correlation is an estimate, so check the y+ the solver reports after the first run and adjust. And y+ varies over a surface. It falls towards zero at separation and reattachment points, which is exactly where wall functions fail.",{"id":92,"title":93,"paragraphs":94},"which-model","Which model for which job?",[95,127],{"table":96},{"head":97,"rows":101},[98,99,100],"Situation","I would start with","Wall treatment",[102,106,109,112,116,120,124],[103,104,105],"External aerodynamics, lift and drag","k-ω SST","Resolved, y+ about 1",[107,104,108],"Separation or adverse pressure gradients, as in diffusers, bends and bluff bodies","Resolved",[110,104,111],"Wall heat transfer","Resolved, because heat transfer is decided in the sublayer",[113,114,115],"Large internal flow, mostly attached, limited cell budget","Realizable k-ε","Wall functions, y+ above 30",[117,118,119],"Free jets and mixing layers far from walls","k-ε family","Not critical",[121,122,123],"Strong swirl or strongly curved streamlines","Consider a Reynolds stress model","Depends on the model",[125,126,105],"Transition from laminar to turbulent flow matters","A transition model based on SST","If you cannot afford to resolve the wall on a large industrial geometry, a k-ε model with wall functions used correctly is a defensible choice. An SST model on a mesh with y+ of 60 everywhere is running on its wall functions and has lost much of the reason you chose it.",{"id":129,"title":130,"paragraphs":131},"limits","What neither model can do",[132,133,134],"Both assume a single turbulent viscosity that is the same in every direction. That fails in strongly swirling flow, in secondary flows in non-circular ducts and in flows with strong streamline curvature.","Both average the turbulence away. If your question is about unsteady loads, noise or mixing driven by large eddies, a RANS model answers a different question from the one you asked, and scale-resolving methods are the tool.","Both were calibrated on simple flows. Far from those flows, treat the result as an estimate and check it against data, or at least against a second model.",{"id":136,"title":137,"paragraphs":138},"common-mistakes","Common mistakes",[139],{"list":140},[141,142,143,144,145],"**Choosing SST by habit and never checking y+.** The model then runs on a wall treatment you did not choose.","**Refining the wall mesh while using standard wall functions.** The first cell moves into the buffer layer and the result gets worse.","**Judging y+ by its average.** Look at the distribution on the surfaces that matter.","**Comparing two models on two different meshes.** The difference then belongs to the mesh as much as to the model.","**Expecting the turbulence model to fix a mesh problem.** No model compensates for cells that are too coarse or too distorted.",{"id":147,"title":148,"paragraphs":149},"short-version","The short version",[150],"Decide how you will treat the wall, which is a decision about cell budget and about whether separation or heat transfer matters. That decision mostly picks the model for you. Then compute the first cell height, mesh to it, and check the y+ you actually got.",[152,155,158],{"href":153,"label":154},"https:\u002F\u002Fdoi.org\u002F10.2514\u002F3.12149","Menter (1994), AIAA Journal: two-equation eddy-viscosity turbulence models for engineering applications",{"href":156,"label":157},"https:\u002F\u002Fdoi.org\u002F10.1016\u002F0045-7825(74)90029-2","Launder and Spalding (1974), Computer Methods in Applied Mechanics and Engineering: the numerical computation of turbulent flows",{"href":159,"label":160},"https:\u002F\u002Fturbmodels.larc.nasa.gov\u002F","NASA Langley Research Center: Turbulence Modeling Resource",1791759706843]