“Turbulence is the most important unsolved problem of classical physics.” (Richard Feynman, American theoretical physicist, Nobel Prize in Physics in 1965)
“I am an old man now, and when I die and go to heaven, there are two matters on which I hope for enlightenment. One is quantum electrodynamics, and the other is the turbulent motion of fluids. And about the former, I am rather optimistic.” (Horace Lamb, English applied mathematician, Advisor: George Gabriel Stokes)
Even a few decades after these great scientists expressed these remarks, turbulence modeling is still no easy thing.
The idea of using fluid to extract or deliver energy is not recent. In the 17th century, Leonardo Da Vinci conducted various experiments to visualize fluid flow, talking about vortex flow, vorticity, swirls, and eddies.
Fluid flow is classified into two main categories—laminar or turbulence—regarding the driven forces (inertial, viscous, etc.), and also a transition regime between them. Notably, being laminar or turbulent is a property of fluid flow under dynamic conditions, not a property of being fluid.
Laminar Flow: A fluid flows through a smooth path with no disruption between adjacent paths. It is quite compatible to examine laminar flow both numerically and experimentally.
Turbulent Flow: A fluid flows through a chaotic path that comprises eddies, swirls, and flow instabilities. It is uneasy, almost impossible in some cases, to examine turbulent flow both numerically and experimentally.
Earlier, the determination of the type of fluid flow numerically was hard to realize. Irish scientist Osborne Reynolds (1883) discovered the dimensionless number that predicts fluid flow based on static and dynamic properties such as velocity, density, dynamic viscosity, and length:
Re=(inertial force)/(viscous force)=ρVL/μ
Where ρ (kg/m^{3}) is the density of the fluid, V (m/s^{2}) is the characteristic velocity of the flow, L (m) is the characteristic length scale of flow, and μ (Pa*s) is the dynamic viscosity of the fluid.
Type 
Flow type 
Reynolds Number 
Internal 
Laminar regime 
up to Re=2300 
Transition regime 
2300<Re<4000 

Turbulent regime 
Re>4000 

External  Laminar to Turbulence 
Re>3×10^{5} 
Table 1: Flow types according to the Reynolds number
If the inertial forces—which resist a change in velocity of an object and cause of the fluid movement—are dominant, the flow is turbulent. On the contrary, if the viscous forces, defined as the resistance to flow, are dominant, the flow is laminar. It is common to generate turbulence for a fluid with low viscosity, though it is rare for fluids with high viscosity. A detailed description of the Reynolds number can be obtained from the SimWiki: What is the Reynolds number?
Aside from laminar behavior, turbulence encompasses several hurdles, and thus requires rigorous effort during experimental and numerical examinations. Turbulence is a type of fluid flow which is unsteady, enormously irregular in space and time, threedimensional, rotational, dissipative (in terms of energy), and diffusive (transport phenomenon) at high Reynolds numbers. Due to those divergences in turbulent flow, extremely smallscale fluctuations emerge in velocity, pressure, and temperature. Despite the fact that direct implementation of fluctuated values into the NavierStokes equation is possible, called a Direct Numerical Solution (DNS), it requires an extreme amount of resources in terms of hardware, software, and human effort. Therefore, an appropriate numerical model should be implemented when modeling turbulent flow.
In this case study, the SimScale cloudbased CAE platform was used to investigate a ducting system and optimize its performance. Download it for free to learn how.
Turbulent Flow Turbulence Modeling
Which turbulence model is convenient for your CFD analysis is a troublesome question. To select an appropriate model and simulate physical incident as accurately as possible, you must:
 Scrutinize the physical incident to understand the phenomenon
 Research the literature in detail to define a suitable model
 If literature is poor, try some models concurrently to get an accurate prediction
 In the validation process, take care of the model that you would apply
Implementation of the turbulence model into the numerical scheme is substantial and makes a big difference to the simulation results. At the first step, a quick examination must be carried out—which pertains to the Reynolds number—to detect the type of fluid flow. For instance, as you keep on with the laminar model (no turbulence) for a fluid flow over a cylinder which is turbulent in the reality, the effect of the driven forces, eddies, vorticities and so forth are destructively negated. The numerical study hereby diverges, as demonstrated in Figure 3.
Although there is a number of miscellaneous turbulence models that investigate the motion of the fluid, these rely on turbulent viscosity, and no universal turbulence model exists yet. Generally, turbulence models are classified regarding governing equation and numerical method used to calculate turbulent viscosity, for which a solution is sought for turbulence. Reynoldsaveraged NavierStokes equations (RANS) and large eddy simulation equations (LES) are the common ones that require a compatible amount of resources during examination against DNS. Beyond that, Unsteady Reynoldsaveraged NavierStokes (URANS), in which motion of the solid body or flow separation causes unsteady flow, has been broadly implemented.
The main purpose of turbulence modeling is to prompt equations to anticipate the timeaveraged velocity, pressure, and temperature fields, without calculating the complete turbulent flow pattern as a function of time as in RANS and LES. It is unnecessary to solve the NavierStokes equations for every value of fluctuation since most engineering problems do not require such a comprehensive solution. The turbulence models can be summarized as follows:
DNS: Direct implementation of fluctuated values into the NavierStokes equation without any turbulence model.
LES: An average turbulence model between DNS and RANS in which filtered NavierStokes equations are used for largescale eddies. An appropriate model is preferred to solve smallscale eddies.
RANS: A mathematical model based on average values of variables for both steadystate and dynamic flows (unsteady for URANS). The numerical simulation is driven by a turbulence model which is arbitrarily selected to find out the effect of turbulence fluctuation on the mean fluid flow.
Commonly Applied Turbulence Models in CFD Analysis
Requiring a modest amount of hardware, computational time, and human effort, RANS/URANS methods, and submodels are highly applied for various computational fluid dynamics problems. The implementation of LES is rare but possible in some cases which specifically need much more computational facilities against URANS/RANS. Here are the most prominent turbulence models you can run in the SimScale environment [3]:
SpalartAllmaras
 Oneequation model
 No wall functions
 Stable with good convergence
 Convenient: Aerodynamics flows, transonic flows over airfoils
 Limitations: Solving shear flows, separated flow, decaying turbulence
kepsilon (on SimScale)
 Twoequation model (turbulent kinetic energy and dissipation)
 Uses wall functions
 Good convergence and low memory requirements
 Convenient: Compressible/incompressible, external flow interactions with complex geometry
 Limitations: Not accurate for noslip walls, adverse pressure gradients, strong curvature into flow, and jet flows
komega
 Twoequation model (turbulent kinetic energy and dissipation)
 Omega used as it is easier to solve than epsilon
 Uses wall functions
 Good convergence and low memory requirements
 Convenient: Similar to kepsilon improved accuracy for internal flows, curvatures, separated flows and jets
 Limitations: Hard to converge and sensitive to initial conditions
komega SST (on SimScale)
 Twoequation model (turbulent kinetic energy and dissipation)
 Wall functions
 Combination of kepsilon (in the outer region of and outside of the boundary layer) and komega (in the inner boundary layer)
 Convenient: Separated flows and jets
 Limitations: Difficult to converge
LES Smagorinsky & LES Spart Allmaras
 Selfsimilarity theory
 Solves small eddies with a subgrid scale model and large ones based on geometry calculation
 Separates velocity field into resolved (large eddies) and subgrid parts (small eddies)
 Convenient: Thermal fatigue, vibration, buoyant flows (ship design)
 Limitations: Difficulties in nearwall regions
Turbulence Models Conclusion
I hope this article helped you better understand turbulence. If you’d like put your knowledge into practice and set up your own CFD analysis, SimScale offers the possibility to carry out simulations in the web browser. To discover all the features provided by the SimScale cloudbased simulation platform, download this document. In addition, here are a few resources that you might find interesting:
 Validated turbulence in the SimScale environment: Turbulent Pipe Flow
 Various online public CFD analysis projects that comprise turbulent flow: SimScale Public Projects
 Turbulent flow tutorials: Spoiler Flow, Front Wing, Incompressible Cyclone
References
 Reynolds, Osborne. “An experimental investigation of the circumstances which determine whether the motion of water shall be direct or sinuous, and of the law of resistance in parallel channels”. Philosophical Transactions of the Royal Society. 174 (0), 1883, P. 935–982
 Bird, R.B., Stewart, W.E. and Lightfoot, E.N. “Transport Phenomena”, 2ndth edition, John Wiley & Sons, 2001, ISBN 0471410772
 https://www.bakker.org/dartmouth06/engs150/10rans.pdf
 https://cfd.mace.manchester.ac.uk/desider/workplan.html
 https://cecas.clemson.edu/~rm/quote.html