- Please check and comment entries here.

# Turbulence Simulation Approaches

## Definition

Turbulent flow can be numerically resolved with different levels of accuracy. Many numerical approaches for solving turbulence have been proposed, such as the Reynolds-Averaged Navier–Stokes (RANS), the Large Eddy Simulation (LES), and Direct Numerical Simulation (DNS) approaches. Among these numerical methods, the RANS approach, specifically the Eddy Viscosity Model (EVM), is widely used for calculating turbulent flows thanks to its relatively high accuracy in predicting the mean flow features and its more limited computational demands. However, this approach suffers from several weaknesses, e.g., compromised accuracy and uncertainties due to assumptions in the model construction and insufficient incorporation of the fluid physics. In the LES approach, the whole eddy range is separated into two parts, namely, the large-scale eddy and subgrid-scale (SGS) eddy. The former can be directly resolved, while the latter is computed using the SGS model. As the computing power rapidly increases, this approach is extensively used to study turbulence physics and to resolve low-to-medium Reynolds number flows.

## 1. RANS

In the RANS approach, instantaneous solution variables in the governing equations are decomposed into the mean and fluctuating components, as expressed in Equation (1).

Substituting this variable expression into the instantaneous continuity and Navier-Stokes (N-S) equations yields the ensemble- or time-averaged forms for single-phase Newtonian flow, as shown in Equations (2) and (3). Henceforth, repeated-suffix summation convention is used in the formulae.

The fluctuating quantities are included in the Reynolds stress tensor (−ρu′iu′j¯¯¯¯¯¯) with six components. In order to close the equation set, the Reynolds stress tensor needs to be appropriately solved. One of the solutions for this closure problem employs the Boussinesq assumption ^{[1]}, which relates the Reynolds stresses to the mean velocity gradients. The advantage of the Boussinesq assumption is the relatively low computational cost due to its simplicity. This works well for the engineering flows, which are dominated by only one turbulent shear stress such as the jet flow, wall boundary layer flow, and mixing layers flow. However, this approach is insensitive to the streamline curvature, rotation and body forces, and it exhibits a poor performance in the flows with a strong anisotropy or stress transport effect ^{[2]}^{[3]}. It also has difficulty in predicting transitional flows.

It is worth noting that the turbulent viscosity used in solving the Reynolds stress terms is a function of the space and flow features, rather than a physical parameter such as the fluid viscosity, which is dependent on the molecular structure of the fluid. Obviously, the turbulent viscosity needs to be solved before computing the Reynolds stress terms. In this paper, two-equation models, which include two additional transport equations, are reviewed. Usually the turbulence kinetic energy (k) is adopted as one equation and the turbulent kinetic energy dissipation rate (ε) or the specific dissipation rate (ω) as another one. The modified versions of the ε/ω model will also be presented here. Due to length restrictions, the zero- and one-equation models are not included in the article, but they may be found elsewhere ^{[4]}^{[5]}^{[6]}.

### 1.1. The k-ε Model

The standard k-ε (SKE) model was originally proposed by Launder and Spalding ^{[7]}. The model has been widely applied for resolving turbulent flows without a severe pressure gradient or strong swirling effect (e.g., plane jet, mixing layer, and boundary layer flows) because of its relatively high robustness, low computational cost, and reasonable accuracy. Equations (4) and (5) show the general form of k and ε. By solving these two transport equations, the turbulent viscosity can be calculated as expressed in Equation (6).

The terms from left to right in Equations (6) and (7) are the respective local time derivative, convection, diffusion, production, sink, and source terms. The SKE model is derived from a fully turbulent flow with high Reynolds numbers. The viscous effect is ignored in the model. However, this cannot be applied in the vicinity of the wall, where the viscous force dominates the flow characteristics. In order to deal with this problem, either a wall function is adopted with the SKE model, or a low Reynolds number model is used. The former confuses the users’ judgement whether the weakness of this method lies in the basic SKE model itself or in the wall function. The latter requires additional functions to modify the standard transport equations. With respect to the low Reynolds number models ^{[8]}^{[9]}^{[10]}^{[11]}^{[12]}^{[13]}^{[14]}, an example proposed by Lam and Bremhorst (LB model) ^{[13]} is presented in Equations (7) and (8).

where:

Compared with the SKE model, different formulations of functions f1, f2 and fμ are developed in the LB model to describe the near-wall behavior better. It has been confirmed that the function fμ has a predominant influence on the model performance, and functions f1, f2 play a secondary role in the performance ^{[15]}.

There are also other modified versions of the SKE model, amongst which the widely used Renormalization Group (RNG) k-ε model ^{[16]} and the Realizable k-ε (RKE) model ^{[17]} are introduced in the following section. The main differences between the RNG k-ε model and the SKE model are the modifications of the turbulent viscosity and ε sink term. A differential equation is analytically derived for effective viscosity μeff to account for the low Reynolds number effect. This feature can improve the predictive ability of the RNG k-ε model for low Reynolds number flows or near-wall flows. Additionally, a new ε destruction term is used to account for the rapid strain by modifying the constant of this term. The RKE model is modified mainly with regard to the turbulent viscosity and the ε equation. By defining a variable Cμ ^{[17]}^{[18]}^{[19]} instead of a constant value in turbulent viscosity formulation, the RKE model satisfies the realizability constraints, i.e., positive values for the normal stresses and the Schwartz inequality for the shear stresses. In order to increase the robustness of the model, a new ε equation is employed based on a dynamic equation for fluctuating vorticity. The new ε equation describes turbulent vortex stretching and turbulent dissipation more appropriately compared to the ε equation in the SKE model. With the modified ε equation, the well-known round-jet anomaly that is a poor prediction of the spreading rate of round/axisymmetric jet may be solved ^{[17]}.

### 1.2. The k-ω Model

The k-ω model is widely used for turbulence modelling ^{[20]}. Different versions of this model have been developed in the last decades ^{[21]}^{[22]}^{[23]}^{[24]}^{[25]}^{[26]}^{[27]}. In this paper, the most well-known k-ω model proposed by Wilcox ^{[21]} is reviewed. The transport equations of this model are presented in Equations (10) and (11), where the calculation of β_{1} refers to the work of Wilcox ^{[21]}.

Unlike the k-ε model, the k-ω model can be integrated through the viscous sublayer without any damping function to account for the low Reynolds number effect with high numerical stability. Therefore, it is well applied in aerodynamic flows ^{[22]}^{[20]}. However, the k-ω model is highly sensitive to the empirical value of ω at the free edge of the turbulent shear layer, which can lead to a large prediction error. In order to solve this problem, a modified model was proposed with combination of the original k-ω model and the SKE model by adding a blending function ^{[22]}. This new model is termed the Baseline k-ω model, which applies the original k-ω model in the near-wall region and switches to the SKE model in the outer region. The Baseline k-ω model has a similar performance to the original k-ω model in boundary layer flows, but the former one avoids the strong freestream dependence. However, both k-ω models fail to predict the onset and amount of separation in adverse pressure gradient flows. Based on the Baseline k-ω model, further modification to eddy viscosity is proposed to account for the transport effects of the principal turbulent shear stress, leading to a significant improvement in predicting the adverse pressure gradient flows ^{[22]}. However, the introduced blending function depends on empiricism (e.g., the distance to wall), limiting its application to flows in complex geometries.

### 1.3. Advanced Eddy Viscosity Model

All the aforementioned methods belong to the category of the Eddy Viscosity Model (EVM). Some other advanced EVMs were developed ^{[28]}^{[29]}^{[30]}^{[31]}^{[32]}. In order to account for the strong anisotropy in the near-wall region, Durbin ^{[28]} adopted a new turbulent viscosity term defined in Equation (12), which is considered to be more appropriate than that defined in Equation (6) in a near-wall region. A separate transport equation for a wall-normal turbulent stress υ2¯¯¯ was proposed and solved with the aid of the elliptic relaxation concept. This model is termed the υ2-f model, where f represents the elliptic relaxation function. Subsequently, several modified versions were proposed with respect to the velocity scale ^{[31]}, the characteristic length ^{[29]}, the function f ^{[32]}, and the variable υ2¯¯¯ ^{[30]}. This model category performs well for pressure-induced separating flow, buoyancy impairing turbulent flow, and backstep flow ^{[31]}^{[32]}^{[33]}.

where:

In addition to the modified versions of the linear EVM, the idea of non-linear EVM has been substantially used ^{[34]}^{[35]}^{[36]}^{[37]}^{[38]}^{[39]}^{[40]}. Even though these modified models demonstrated certain improvements over linear EVMs in predicting flows with a strong streamline curvature or turbulent stresses in the near-wall sublayer, they are still inferior to the more advanced model, e.g.

### 1.4. Reynolds Stress Model

In order to overcome the limitations of the EVM, Second-Moment Closure (SMC) models abandoning the Boussinesq assumption have been developed. The SMC model directly solves the transport equation for each of the Reynolds stress terms. Since the SMC approach considers the effects of streamline curvature, rotation, and rapid change of strain rate in a more rigorous manner, it is long expected to replace the currently widely applied two-equation models. The SMC model class consists of the Algebraic Stress Model (ASM) and the Differential Stress Model (DSM). ASM is derived from differential stress transport equations by invoking the weak-equilibrium assumption ^{[41]}^{[42]}^{[43]}. It ignores the transport terms of the anisotropy by assuming that the transport of the Reynolds stress is proportional to that of turbulent kinetic energy ^{[44]}^{[45]}. In general, the ASM is considered an intermediate tool between the LEVM and the DSM. Due to the space limitation, only the DSM is presented in this paper. A symbolic representation of the stress transport equation is expressed in Equation (14). In addition, a scale-determining equation, i.e., the ε equation, is needed to complete the DSM.

The terms from left to right represent the local time derivative of Reynolds stress, convection, turbulent diffusion, molecular diffusion, stress production, buoyancy production, pressure strain, dissipation and source, respectively. The sij term is user-defined for a specific stress transport source. If there is no source, this term becomes zero. It is required to model DT,ij, Gij, ϕij and εij to close the equation, while it is not necessary to model Lij, Cij, DL,ij and Pij.

The DT,ij term includes the velocity transport and the pressure transport. The velocity triple moments can be measured, whereas the pressure transport is intractable. Usually, the pressure transport is considered to be negligible ^{[46]}. Therefore, the model is mainly designed for the velocity triple moments. The most popular model is the generalized gradient-diffusion model proposed by Daly and Harlow (DH) ^{[47]}. The DH model has a symmetry problem in the indices, leading to dependence on the coordinate frame. Subsequently, some variants of this model were developed by Hanjalic and Launder ^{[48]}, Shir ^{[49]}, Mellor and Herring ^{[50]}. More complex models were also put forward by Nagano and Tagawa ^{[51]} and Magnaudet ^{[52]}. Due to uncertainties in the modelling equations, the complex models may not necessarily outperform the simplified models. Given more computing resources consumed by the complex models and a poor convergence, application of these models in engineering has been doubted ^{[53]}^{[54]}.

Compared to the εij and Gij terms, it is necessary to pay extra attention to model ϕij. Usually, the ϕij term is decomposed into three parts, namely, the slow pressure-strain term ϕij,1, rapid pressure-strain term ϕij,2, and wall-reflection term ϕij,w. Not all of the models introduced below include the third term. Rotta ^{[55]} proposed a linear model for ϕij,1, which considers that ϕij,1 is proportional to the stress anisotropy tensor. However, this linear model is unable to satisfy the realizability constraints. A general quadratic model ^{[56]} is proposed to solve this problem. Linear ^{[48]}^{[57]}^{[58]}^{[59]} and nonlinear models ^{[60]}^{[61]}^{[62]} were proposed to model the ϕij,2 term. Even though the nonlinear model is considered to be theoretically advanced, the complexity of the formulation prohibits its application for engineering computation. The turbulence anisotropy is enhanced due to the damping effect of the normal stress by the wall. This damping affects both the pressure-strain terms. In order to account for the damping effect, a commonly used model ^{[49]}^{[63]} is presented to model the wall-reflection term ϕij,w. However, this model involves a variable, i.e., the normal distance to the wall. This is believed to be a major weakness. For purpose of overcoming this weakness, the elliptic relaxation concept and elliptic-blending method were proposed to account for the near-wall inhomogeneity, and more information on that can be found in ^{[64]}^{[65]}.

The DSM is the most elaborate model in the RANS approach, which has an indisputable superiority over the rudimentary two-equation models in predicting complex flows, e.g., highly swirling and rotating flow, separating flow, and secondary flow. However, its application is limited by (1) a high degree of uncertainty in modelling the high-order correlation terms (e.g., pressure-strain and dissipative correlation) due to an insufficient knowledge of physics; (2) a high demand for computational resource (approximately 50–150% more computing time than a two-equation EVM ^{[2]}). Fortunately, due to the use of more advanced models (e.g., the DNS), the accuracy and robustness of the DSM have been improved. A rapid development in computer science (e.g., parallel processing and improved performance) satisfies the high computational need for the use of DSM model. The DSM has received more attention recently because of the failure of the EVM in predicting complex turbulent flows.

## 2. LES Approach

### 2.1. Formulation and Subgrid-Scale Model

Turbulent flow features a wide range of eddy scales from the Kolmogorov length scale to the size comparable to the characteristic length of the mean flow. The large eddies contain most of the turbulent energy and are mainly responsible for the momentum and energy transfer. They are strongly affected by boundary conditions. The small eddies tend to be more isotropic and homogeneous, and their dissipation process is linked to fluid viscosity. For this reason, it is very difficult for the RANS approach to model all the eddies in a single model. The LES approach separates the large eddies from the small ones by employing a spatial filtering method ^{[66]} for the instantaneous governing equations. After that, the large eddies are directly resolved by the filtered equation, and the small ones (i.e., the Subgrid-Scale (SGS) eddies) are modelled by the SGS model. The filtered variable (donated by an overbar) is defined by Equation (15). The resulting continuity and momentum equations are expressed in Equations (16) and (17), showing similar forms but a different physical meaning for those in the RANS approach.

where σij refers to the stress tensor due to molecular viscosity, and τij represents the SGS stress term. Most of current SGS models adopt the Boussinesq assumption (termed the Eddy-Viscosity model), which relates the SGS stress to the large-scale strain-rate tensor, as shown in Equation (18). Based on the definition of the eddy viscosity, various SGS models ^{[67]}^{[68]}^{[69]}^{[70]}^{[71]} have been proposed. Smagorinsky ^{[67]} developed the first SGS model by assuming a local energy equilibrium between the large scale and the subgrid scale. The eddy viscosity in this model is defined in Equation (20). This model becomes very popular to date due to its simplicity, numerical robustness, and stability. However, it has several drawbacks: (1) the model constant varies with different flows; (2) the model cannot predict the inverse energy transfer (i.e., backscatter) due to its purely dissipative nature; (3) the model has difficulty in reproducing the correct mean quantities (e.g., SGS dissipation) as the grid scale approaches the integral scale; (4) the model does not yield a zero eddy viscosity in near-wall regions.

where:

In Equation (20), where:

where Δ represents the local grid scale. In order to solve the model constant problem, a Dynamic Smagorinsky–Lilly (DSL) model was proposed ^{[68]}^{[69]}, where the model constant is dynamically calculated by using the resolved eddies with the scale size between the grid filter and test filter. The main advantage of this model is that it is not necessary to prescribe and/or tune the model constant. However, the DSL model is subjected to a numerical instability and a variable model constant. Germano ^{[68]} proposed an averaging method to overcome this weakness. A good performance was achieved in a channel flow simulation ^{[72]}. Another variant of Smagorinsky–Lilly model is the Dynamic Kinetic Energy (DKE) model ^{[73]}^{[74]}^{[75]}. Unlike the algebraic form in Smagorinsky–Lilly and DSL models, the DKE model solves an additional transport equation for the SGS turbulent kinetic energy instead of adopting the local equilibrium assumption. This model can better account for the energy transfer from the large-scale eddy at the cost of computational expenses. Some other variants of the Smagorinsky–Lilly model were formulated to solve the low Reynolds number effect in the near-wall region, one of which was based on the square of the velocity gradient tensor named the Wall-Adapting Local Eddy-viscosity (WALE) model ^{[70]}. Compared to the original Smagorinsky–Lilly model, the WALE model can produce a zero eddy viscosity in the vicinity of the wall or in a pure shear flow. Hence, this model does not need a damping function. In addition to the WALE model, a hybrid model ^{[76]} was proposed by combining the Smagorinsky–Lilly model with a damping function ^{[77]} to improve the predictive capability for wall-bounded flows. This hybrid model demonstrated a good performance in plane channel flow with different Reynolds numbers. However, this model involves a variable, i.e., the wall normal distance. The determination of this wall normal distance requires an empirical approach for specific flow.

In the framework of the eddy-viscosity SGS model, there are several alternatives to the Smagorinsky-type SGS model, such as Vreman’s model ^{[78]}, the QR model ^{[79]}^{[80]}, the σ-model ^{[81]}, and the S3PQR model ^{[82]}. Compared to the Smagorinsky-type SGS model, Vreman’s model can predict zero eddy viscosity in near-wall regions or in transitional flows without explicit filtering, averaging or clipping procedures. However, it was found that the model coefficient in Vreman’s model is far from universal. To solve this problem, two procedures were proposed to dynamically determine the model coefficient, i.e., the one based on the global equilibrium between the subgrid-scale dissipation and the viscous dissipation ^{[83]}^{[84]} and the other one based on the Germano identity ^{[85]}. It was reported that the latter is better suited for transient flows ^{[85]}. The QR model, which is a minimum-dissipation eddy-viscosity model, gives the minimum eddy dissipation required to dissipate the energy of sub-filter scales. The advantages of this model lie in appropriately switching off for laminar and transitional flows, the low computational complexity, and consistency with the exact sub-filter tensor on isotropic grids. The disadvantage of this model is the insufficient eddy dissipation, which can be corrected by increasing the model constant. Moreover, the QR model requires an approximation of the filter width to be consistent with the exact sub-filter tensor on anisotropic grids. It was noted that the accuracy of the model result for anisotropic grids is highly dependent on the used filter width approximation. By modifying the Poincaré inequality used in the QR model, the dependence can be removed, leading to the construction of an anisotropic minimum-dissipation model that generalizes the desirable properties of the QR model to anisotropic grids ^{[86]}. For the purpose of meeting a set of properties based on the practical/physical considerations, the σ-model based on the singular values of the velocity gradient tensor was developed ^{[81]}. Owing to its unique properties, ease of implementation, and low computational cost, the σ-model is considered to be suitable for complex flow configurations. Subsequently, through comparison between static and dynamic σ-models, it was found that the local dynamic procedure is not suited for the σ-model, and a global dynamic procedure is suggested ^{[87]}. Trias et al. ^{[82]} built a general framework for LES eddy-viscosity models, which is based on the 5D phase space of invariants. By imposing appropriate restrictions in this space, a new eddy-viscosity model, i.e., the S3PQR model, was developed. In addition to meeting a set of desirable properties such as positiveness, locality, Galilean invariance, proper near-wall behavior, and automatic switch-off for laminar, 2D and axisymmetric flows, this new model is well-conditioned and has a low computational cost, with no intrinsic limitations for statistically inhomogeneous flows. Despite of all the merits of this model, special attention should be given to the calculation of the characteristic length scale and the determination of the model constant before engineering applications.

Alternatives to the eddy-viscosity SGS model are the similarity model ^{[88]}^{[89]}^{[90]}^{[91]}, the velocity estimation model ^{[92]}^{[93]}, the Approximate Deconvolution model (ADM) ^{[94]}^{[95]}^{[96]}, and the regularization model ^{[97]}^{[98]}^{[99]}. The similarity model class adopts the idea that an accurate approximation for a SGS model can be reconstructed from the information contained in the resolved field. Therefore, the similarity models ^{[88]}^{[89]}^{[90]}^{[91]} approximate the SGS stress tensor by a stress tensor calculated from the resolved scales. Due to this nature, the similarity model can naturally account for the inverse energy transfer (i.e., backscatter). This is different from the eddy viscosity model, which only considers the global SGS dissipation (i.e., the net energy flux from the resolved scales to the subgrid scales). It is worth mentioning that, due to the importance of accurate prediction of the inverse energy transfer, a dynamic two-component SGS model was proposed to include the non-local and local interactions between the resolved scales and subgrid scales ^{[91]}. The model correctly predicted the breakdown of the net transfer into forward and inverse contributions in a priori tests. In some cases, however, the similarity model underestimated the SGS dissipation. An extra dissipative term was added to solve this issue. This new model formulated is also referred to as the mixed model. Furthermore, the similarity model and the mixed model need additional computational resources due to the implementation of the second filtering. Special attention should be paid to choose an appropriate filter type and size. Following the same idea, Domaradzki et al. ^{[92]} improved the SGS stress approximation by replacing the unknown unfiltered variables by approximately deconvolved filtered variables. Subsequently, this SGS model based on the estimation of the unfiltered velocity, which was originally formulated in spectral space, was extended to the physical space ^{[93]}. It was found that both versions of this velocity estimation model perform better than or are comparable to classical eddy viscosity models for most physical quantities. This model can account for backscatter without any adverse effects on the numerical stability. Several questions for improving the model need to be addressed, such as the modelling of nonequilibrium and high Reynolds number turbulence in three Cartesian directions. Stolz et al. ^{[94]}^{[95]}^{[96]} proposed a formulation of the ADM for LES, in which an approximation of the unfiltered solution is obtained from the filtered solution by a series expansion involving repeated filtering. Given a good approximation of the unfiltered solution at a time instant, the nonlinear flux terms of the filtered N-S equations can be computed directly, avoiding the explicit computation of the SGS closures. The effect of the non–represented scales is modelled by a relaxation regularization involving a second filtering and a dynamically estimated relaxation parameter. The ADM is evaluated for incompressible wall-bounded flow ^{[95]} and compressible flows ^{[94]}^{[96]}. The results showed that the ADM can have a significant improvement over the standard and dynamic Smagorinsky models, while at less computational cost compared with that of the dynamic models or the velocity estimation model. The high-Reynolds-number supersonic flow ^{[100]} and transitional flow ^{[101]} were investigated by the ADM. Agreement was observed between the ADM and experiments or DNS. For the former flow, a rescaling and recycling method was used to have a better control on the desired inflow data. Recently, the ADM was extended for a two-phase flow simulation ^{[102]}. By comparing the macroscopic flow characteristics, the ADM showed a better performance than the conventional LES model. However, further investigations should be performed on the relaxation term model for a two-phase simulation and microscopic characteristics of the dispersed phase in a 3D simulation. Another important class of the SGS model is the regularization model, which combines a regularization principle with an explicit filter and its inversion. The regularization model includes many versions, such as the Leray model ^{[97]}, the Leray-α model ^{[98]}^{[103]}, the Clark-α model ^{[99]}, the Navier–Stokes-α (NS-α) model ^{[104]}, etc. For the last one, several variants are proposed, including the NS-α deconvolution model ^{[105]}^{[106]}^{[107]}^{[108]} and the reduced order NS-α model ^{[109]}^{[110]}. It was found that the NS-α deconvolution model can significantly improve the prediction accuracy by carefully choosing the filtering radius and by correctly selecting the approximate deconvolution order ^{[108]}. Given the difficulties in efficient and stable simulation of the NS-α model for incompressible flows on coarse grids, the reduced order NS-α model is introduced by using deconvolution as an approximation to the filter inverse, reducing the fourth-order NS-α formulation to a second-order model. In spite of the success of the reduced order NS-α model, future work needs to be conducted on locally and dynamically choosing α and numerical testing on different benchmark flows, to name but a few. In addition, comparative studies have been performed between different regularization models ^{[111]}^{[112]}^{[113]}, in which the capability of the regularization model has been demonstrated for a specific flow.

The LES is considered a compromise between the RANS and the DNS. It is more accurate than the RANS and it needs less computational resources than the DNS. However, the LES model has not reached the maturity stage as a numerical tool for the design or the parametric study of complex engineering flows, due to not only a high computational requirement, but also many unresolved issues such as ill-defined boundary conditions, wall-resolved flow, turbulent flow with chemical reactions, and compressible flow. Nevertheless, the LES model has been successfully applied in transitional flow ^{[114]}^{[115]}^{[116]}^{[117]}, separated flow ^{[118]}^{[119]}, and bubbly flow ^{[120]}^{[121]}. **Figure 1** shows the calculation results in a separated boundary layer transition on a flat plate with a semi-circular leading edge of radius of 5 mm under elevated free-stream turbulence. A periodic boundary condition was adopted in the spanwise direction. Free-slip and no-slip conditions were used at lateral boundaries and the plate surface, respectively. The simulation agrees well with experimental data on mean and fluctuating streamwise velocities for an Enhanced-Turbulence-Level (ETL) case, demonstrating a good performance of the LES model for the transitional flow. **Figure 2** shows a comparison between a modified k-ε model and the LES model in a bubbly flow. Compared with the experimental data, the LES model is superior in predicting the turbulence.

**Figure 1.** (**a**) Mean streamwise velocity, (**b**) rms streamwise velocity at different streamwise stations (ETL-case) (LES: solid lines, experimental data: symbols, x/l: normalized streamwise, y/l: wall-normal direction, U/U_{0}: normalized mean streamwise velocity, u’/U_{0}: normalized rms streamwise velocity) ^{[117]}.

**Figure 2.** Comparison of model predictions with experimental data for (**a**) axial and (**b**) lateral liquid velocity fluctuations at a height of 0.25 m in a square cross-sectioned bubble column with a superficial gas velocity of 4.9 × 10^{−3} m/s ^{[120]}.

### 2.2. Inlet Boundary Condition

The fluid behavior in the domain is largely determined by the inflow condition ^{[122]}. The treatment of the inflow condition is of significant importance for LES modelling. Currently, there are two main categories for generating the inflow data, namely, precursor simulation and a synthetic method. The former involves a separate simulation, where the periodic boundary condition or the recycling method can be used. The flow data are stored at each time step in this simulation and then introduced to the inlet boundary for modelling the flow of interest. The main advantage of this method is to obtain more realistic inflow conditions, which represent the required flow characteristics (e.g., velocity profile, turbulence intensity, shear stress, power spectrum, and turbulent structures). The inlet boundary, however, needs to be placed in an equilibrium region for scaling arguments in the precursor simulation, which may even not exist in some flows. This method may lead to a spurious periodicity for the time series ^{[123]}. In addition, running a separate simulation requires high computational costs especially for a high-Reynolds number flow. This restrains its application to complex engineering flows. The limitation can be reduced by an internal mapping method. This method integrates the precursor simulation into the main domain, mapping the data downstream of the inlet back to the inlet boundary ^{[124]}^{[125]}.

The synthetic method as an alternative to the precursor simulation is expected to construct the inflow conditions for practical flows. The simplest way is to impose a white-noise random component on the inlet velocity. The magnitude of this random component is determined by the turbulent intensity. Since the turbulence-like component is rapidly dissipated due to the lack of spatial and temporal correlation, this white-noise method is inappropriate to generate the inflow data ^{[126]}. In order to impose realistic inflow data on the inlet boundary, other advanced synthetic methods have been developed. These advanced methods consist of the Fourier technique ^{[127]}^{[128]}, principal orthogonal decomposition (POD) method ^{[129]}^{[130]}, digital filtering technique ^{[131]}^{[132]}, and synthetic eddy method (SEM) ^{[133]}^{[134]}^{[135]}. Several comparative studies have been performed on different synthetic methods ^{[134]}^{[135]}^{[136]}^{[137]}. Jarrin et al. ^{[134]} used the SEM in the hybrid RANS/LES simulations for turbulent flows from simple channel and square duct flows to the flow over an airfoil trailing edge. Compared to other synthetic methods (i.e., Batten’s method (Fourier method) ^{[127]} and random method), the SEM can substantially reduce the inlet section, leading to a large decrease in the CPU time. **Figure 3** shows a better performance of the SEM in the inlet velocity vector compared to that of the others. The SEM realistically reproduced the magnitude and length scale of the fluctuations. The fluctuations in the Batten’s model decorrelated in space in the near-wall region due to the decomposition in Fourier’s mode. Recently, Skillen et al. ^{[135]} improved the SEM of Jarrin et al. (Original SEM) ^{[133]} with respect to the normalization algorithm and the eddy placement. The former leads to an improvement over the original SEM model. The latter saves a cost of around 1–2 orders of magnitude. **Figure 4** and **Figure 5** show the turbulent shear stresses and skin-friction coefficients from the original SEM ^{[133]}, the improved SEM of Skillen et al. ^{[135]}, and the precursor LES of Kaltenbach et al. ^{[138]}. In comparison with the original SEM, the improved SEM shows a better agreement with the precursor LES results both for the turbulent shear stress and the skin-friction coefficient. Although there is a rapid development of synthetic methods, the available synthetic methods are limited to construct the inflow condition with all required turbulence characteristics mentioned above. Further development is needed. It is unnecessary to claim which method is the best; however, the most appropriate method can be selected by considering the accuracy and the computational cost.

**Figure 3.** Velocity vectors of the LES inlet condition for hybrid simulation of channel flow: (**a**) precursor LES, (**b**) SEM, (**c**) Batten’s method, (**d**) random method (y/δ: normalized y distance, z/δ: normalized z distance, δ: initial boundary layer thickness) ^{[134]}.

**Figure 4.** Profiles of turbulent shear stress (the number at the top of the figure indicates profile locations (the distance from the start of the incline, normalized by the initial boundary layer thickness)) ^{[135]}.

**Figure 5.** Wall shear stress for the inclined wall (Cf: wall shear stress, x/δ: normalized distance from the start of the incline) ^{[135]}.

The entry is from 10.3390/met11081297

## References

- Hinze, J.O. Turbulence; McGraw-Hill Publishing Co.: New York, NY, USA, 1975.
- Launder, B.E. Second moment closure and its use in modelling turbulent industrial flows. Int. J. Numer. Meth. Fluids 1989, 9, 963–985.
- Hanjalić, K. Introduction to Turbulence Modelling; Von Kaman Institute for Fluid Dynamics: Sint-Genesius-Rode, Belgium, 2004; pp. 1–75.
- Prandtl, L. Bericht fiber untersuchungen zur ausgebildeten turbulent. ZAMM 1925, 5, 136–139.
- Spalart, P.R.; Allmaras, S.R. A one-equation turbulence model for aerodynamic flows. In Proceedings of the 30th Aerospace Sciences Meeting and Exhibit, Reno, NV, USA, 6–9 January 1992.
- Nallasamy, M. Turbulence models and their applications to the prediction of internal flows: A review. Comput. Fluids 1987, 15, 151–194.
- Launder, B.E.; Spalding, D.B. The numerical computation of turbulent flows. Meth. Appl. Mech. Eng. 1974, 3, 269–289.
- Jones, W.P.; Launder, B.E. The prediction of laminarization with a two-equation model of turbulence. Int. J. Heat Mass Trans. 1972, 15, 301–314.
- Jones, W.P.; Launder, B.E. The calculation of low-Reynolds-number phenomena with a two-equation model of turbulence. Int. J. Heat Mass Trans. 1973, 16, 1119–1130.
- Launder, B.E.; Sharma, B.I. Application of the energy-dissipation model of turbulence to the calculation of flow near a spinning disc. Lett. Heat Mass Trans. 1974, 1, 131–138.
- Hoffmann, G.H. Improved form of the low Reynolds number k−ε turbulence model. Phys. Fluids 1975, 18, 309–312.
- Hassid, S.; Poreh, M. A turbulent energy dissipation model for flows with drag reduction. J. Fluids Eng. 1978, 100, 107–112.
- Lam, C.K.G.; Bremhorst, K. A modified form of the k-ε model for predicting wall turbulence. Tans. ASME 1981, 103, 456–460.
- Chien, K.Y. Predictions of channel and boundary-layer flows with a low-Reynolds-number turbulence model. AIAA J. 1982, 20, 33–38.
- Patel, V.C.; Rodi, W.; Scheuerer, G. Turbulence models for near-wall and low Reynolds number flows—A review. AIAA J. 1985, 23, 1308–1319.
- Yakhot, V.; Orszag, S.A. Renormalization group analysis of turbulence: Basic theory. J. Sci. Comput. 1986, 1, 3–51.
- Shih, T.H.; Liou, W.W.; Shabbir, A.; Yang, Z.; Zhu, J. A new k-ϵ eddy viscosity model for high Reynolds number turbulent flows. Comput. Fluids 1995, 24, 227–238.
- Reynolds, W.C. Fundamentals of Turbulence for Turbulence Modeling and Simulation in Lecture Notes for Von Karman Institute; AGARD Lecture Note Series: New York, NY, USA, 1987.
- Shih, T.H.; Zhu, J.; Lumley, J.L. A new Reynolds stress algebraic model. NASA TM 1995, 125, 1–4.
- Menter, F.R. Review of the shear-stress transport turbulence model experience from an industrial perspective. Int. J. Comput. Fluid Dyn. 2009, 23, 305–316.
- Wilcox, D.C. Reassessment of the scale-determining equation for advanced turbulence models. AIAA J. 1988, 26, 1299–1310.
- Menter, F.R. Two-equation eddy-viscosity turbulence models for engineering applications. AIAA J. 1994, 32, 1598–1605.
- Wilcox, D.C. Comparison of two-equation turbulence models for boundary layers with pressure gradient. AIAA J. 1993, 31, 1414–1421.
- Wilcox, D.C. Simulation of transition with a two-equation turbulence model. AIAA J. 1994, 32, 247–255.
- Peng, S.H.; Davidson, L.; Holmberg, S. A modified low-Reynolds-number k-ω model for recirculating flows. J. Fluid Eng. 1997, 119, 867–875.
- Bredberg, J.; Peng, S.H.; Davidson, L. An improved k-ω turbulence model applied to recirculating flows. Int. J. Heat Fluid Flow 2002, 23, 731–743.
- Resende, P.R.; Pinho, F.T.; Younis, B.A.; Kim, K.; Sureshkumar, R. Development of a low-Reynolds-number k-ω model for FENE-P fluids. Flow Turbul. Combust. 2013, 90, 69–94.
- Durbin, P.A. Near-wall turbulence closure modeling without “damping functions”. Theor. Comput. Fluid Dyn. 1991, 3, 1–13.
- Durbin, P.A. On the k-ε stagnation point anomaly. Int. J. Heat Fluid Flow 1996, 17, 89–90.
- Laurence, D.R.; Uribe, J.C.; Utyuzhnikov, S.V. A robust formulation of the v2−f model. Flow Turb. Combust. 2004, 73, 169–185.
- Hanjalić, K.; Popovac, M.; Hadžiabdić, M. A robust near-wall elliptic-relaxation eddy-viscosity turbulence model for CFD. Int. J. Heat Fluid Flow 2004, 25, 1047–1051.
- Billard, F.; Laurence, D. A robust k-ε v ¯ 2k elliptic blending turbulence model applied to near-wall, separated and buoyant flows. Int. J. Heat Fluid Flow 2012, 33, 45–58.
- Billard, F.; Revell, A.; Craft, T. Application of recently developed elliptic blending based models to separated flows. Int. J. Heat Fluid Flow 2012, 35, 141–151.
- Speziale, C.G. On nonlinear kl and k-ε models of turbulence. J. Fluid Mech. 1987, 178, 459–475.
- Lien, F.S.; Leschziner, M.A. Assessment of turbulence-transport models including non-linear RNG eddy-viscosity formulation and second-moment closure for flow over a backward-facing step. Comput. Fluids 1994, 23, 983–1004.
- Craft, T.J.; Launder, B.E.; Suga, K. Development and application of a cubic eddy-viscosity model of turbulence. Int. J. Heat Fluid Flow 1996, 17, 08–115.
- Craft, T.J.; Launder, B.E.; Suga, K. Prediction of turbulent transitional phenomena with a nonlinear eddy-viscosity model. Int. J. Heat Fluid Flow 1997, 18, 15–28.
- Apsley, D.D.; Leschziner, M.A. A new low-Reynolds-number nonlinear two-equation turbulence model for complex flows. Int. J. Heat Fluid Flow 1998, 19, 209–222.
- Jongen, T.; Machiels, L.; Gatski, T.B. Predicting noninertial effects with linear and nonlinear eddy-viscosity, and algebraic stress models. Flow Turb. Combust. 1998, 60, 215–234.
- Hellsten, A.; Wallin, S. Explicit algebraic Reynolds stress and non-linear eddy-viscosity models. Int. J. Comput. Fluid Dyn. 2009, 23, 349–361.
- Gatski, T.B.; Speziale, C.G. On explicit algebraic stress models for complex turbulent flows. J. Fluid Mech. 1993, 254, 59–78.
- Gatski, T.B.; Jongen, T. Nonlinear eddy viscosity and algebraic stress models for solving complex turbulent flows. Prog. Aerosp. Sci. 2000, 36, 655–682.
- Gatski, T.B.; Wallin, S. Extending the weak-equilibrium condition for algebraic Reynolds stress models to rotating and curved flows. J. Fluid Mech. 2004, 518, 147–155.
- Rodi, W. The Prediction of Free Turbulent Shear Flows by Use of a 2-Equation Model of Turbulence. Ph.D. Thesis, University of London, London, UK, 1972.
- Rodi, W. A new algebraic relation for calculating the Reynolds stresses. ZAMM 1976, 56, 219–221.
- Kim, J.; Moin, P.; Moser, R. Turbulence statistics in fully developed channel flow at low Reynolds number. J. Fluid Mech. 1987, 177, 133–166.
- Daly, B.J.; Harlow, F.H. Transport equations in turbulence. Phys. Fluids 1970, 13, 2634–2649.
- Hanjalić, K.; Launder, B.E. A Reynolds stress model of turbulence and its application to thin shear flows. J. Fluid Mech. 1972, 52, 609–638.
- Shir, C.C. A preliminary numerical study of atmospheric turbulent flows in the idealized planetary boundary layer. J. Atmospher. Sci. 1973, 30, 1327–1339.
- Mellor, G.L.; Herring, H.J. A survey of the mean turbulent field closure models. AIAA J. 1973, 11, 590–599.
- Nagano, Y.; Tagawa, M. Turbulence model for triple velocity and scalar correlations. In Turbulent Shear Flows 7; Springer: Berlin, Germany, 1991; Volume 7, pp. 47–62.
- Magnaudet, J. The modelling of inhomogeneous turbulence in the absence of mean velocity gradients. In Proceedings of the 4th European Turbulence Conference, Delft, The Netherlands, 30 June–3 July 1992.
- Hanjalić, K. Advanced turbulence closure models: A view of current status and future prospects. Int. J. Heat Fluid Flow 1994, 15, 178–203.
- Launder, B.E. Phenomenological modelling: Present and future? In Proceedings of the Workshop Held at Cornell University, Ithaca, NY, USA, 22–24 March 1989.
- Rotta, J.C. Statistische Theorie nichthomogener Turbulenz. Z. Phys. 1951, 129, 547–572.
- Lumley, J.L. Computational modeling of turbulent flows. Adv. Appl. Mech. 1979, 18, 123–176.
- Naot, D. Interactions between components of the turbulent velocity correlation tensor. Isr. J. Technol. 1970, 8, 259–269.
- Launder, B.E.; Reece, G.J.; Rodi, W. Progress in the development of a Reynolds-stress turbulence closure. J. Fluid Mech. 1975, 68, 537–566.
- Speziale, C.G.; Sarkar, S.; Gatski, T.B. Modelling the pressure-strain correlation of turbulence: An invariant dynamical systems approach. J. Fluid Mech. 1991, 227, 245–272.
- Shih, T.H.; Lumley, J.L. Modeling of Pressure Correlation Terms in Reynolds Stress and Scalar Flux Equations; Cornell University: Ithaca, NY, USA, 1985.
- Craft, T.J.; Launder, B.E. Computation of impinging flows using second-moment closures. In Proceedings of the 8th Symposium on Turbulent Shear Flows, Munich, Germany, 9–11 September 1991; pp. 851–856.
- Launder, B.E.; Tselepidakis, D.P. Progress and paradoxes in modelling near-wall turbulence. In Proceedings of the 8th Symposium on Turbulent Shear Flows, Munich, Germany, 9–11 September 1991; pp. 2911–2916.
- Gibson, M.M.; Launder, B.E. Ground effects on pressure fluctuations in the atmospheric boundary layer. J. Fluid Mech. 1978, 86, 491–511.
- Durbin, P.A. A Reynolds stress model for near-wall turbulence. J. Fluid Mech. 1993, 249, 465–498.
- Manceau, R.; Hanjalić, K. Elliptic blending model: A new near-wall Reynolds-stress turbulence closure. Phys. Fluids 2002, 14, 744–754.
- Leonard, A. Energy cascade in large-eddy simulations of turbulent fluid flows. Adv. Geophys. 1975, 18, 237–248.
- Smagorinsky, J. General circulation experiments with the primitive equations. Month. Wea. Rev. 1963, 91, 99–164.
- Germano, M.; Piomelli, U.; Moin, P.; Cabot, W.H. A dynamic subgrid-scale eddy viscosity model. Phys. Fluids 1991, 3, 1760–1765.
- Lilly, D.K. A proposed modification of the Germano subgrid-scale closure method. Phys. Fluids 1992, 4, 633–635.
- Nicoud, F.; Ducros, F. Subgrid-scale stress modelling based on the square of the velocity gradient tensor. Flow Turb. Combust. 1999, 62, 183–200.
- Schumann, U. Subgrid scale model for finite difference simulations of turbulent flows in plane channels and annuli. J. Comput. Phys. 1975, 18, 376–404.
- Piomelli, U. High Reynolds number calculations using the dynamic subgrid-scale stress model. Phys. Fluids 1993, 5, 1484–1490.
- Kim, W.W.; Menon, S. A new dynamic one-equation subgrid-scale model for large eddy simulations. In Proceedings of the AIAA 34th Aerospace Sciences Meeting and Exhibit, Reno, NV, USA, 15–18 January 1996.
- Kajishima, T.; Nomachi, T. One-equation subgrid scale model using dynamic procedure for the energy production. J. Appl. Mech. 2006, 73, 368–373.
- Huang, S.; Li, Q.S. A new dynamic one equation subgrid scale model for large eddy simulations. Int. J. Numer. Meth. Eng. 2010, 81, 835–865.
- Shur, M.L.; Spalart, P.R.; Strelets, M.K.; Travin, A.K. A hybrid RANS-LES approach with delayed-DES and wall-modelled LES capabilities. Int. J. Heat Fluid Flow 2008, 29, 1638–1649.
- Piomelli, U.; Moin, P.; Ferziger, J.H. Model consistency in large eddy simulation of turbulent channel flows. Phys. Fluids 1988, 31, 884–1891.
- Vreman, A.W. An eddy-viscosity subgrid-scale model for turbulent shear flow: Algebraic theory and applications. Phys. Fluids 2004, 16, 3670–3681.
- Verstappen, R.W.C.P.; Bose, S.T.; Lee, J.; Choi, H.; Moin, P. A dynamic eddy-viscosity model based on the invariants of the rate-of-strain. In Proceedings of the Summer Program of the Center for Turbulence Research, Stanford, CA, USA, 2010; pp. 183–192.
- Verstappen, R.W.C.P. When does eddy viscosity damp subfilter scales sufficiently? J. Sci. Comput. 2011, 49, 94–110.
- Nicoud, F.; Toda, H.B.; Cabrit, O.; Bose, S.; Lee, J. Using singular values to build a subgrid-scale model for large eddy simulations. Phys. Fluids 2011, 23, 085106.
- Trias, F.X.; Folch, D.; Gorobets, A.; Oliva, A. Building proper invariants for eddy-viscosity subgrid-scale models. Phys. Fluids 2015, 27, 065103.
- Park, N.; Lee, S.; Lee, J.; Choi, H. A dynamic subgrid-scale eddy viscosity model with a global model coefficient. Phys. Fluids 2006, 18, 125109.
- You, D.; Moin, P. A dynamic global-coefficient subgrid-scale eddy-viscosity model for large-eddy simulation in complex geometries. Phys. Fluids 2007, 19, 065110.
- Lee, J.; Choi, H.; Park, N. Dynamic global model for large eddy simulation of transient flow. Phys. Fluids 2010, 22, 075106.
- Rozema, W.; Bae, H.J.; Moin, P.; Verstappen, R. Minimum-dissipation models for large-eddy simulation. Phys. Fluids 2015, 27, 085107.
- Rieth, M.; Proch, F.; Stein, O.T.; Pettit, M.W.A.; Kempf, A.M. Comparison of the Sigma and Smagorinsky LES models for grid generated turbulence and a channel flow. Comput. Fluids 2014, 99, 172–181.
- Bardina, J.; Ferziger, J.H.; Reynolds, W.C. Improved subgrid-scale models for large-eddy simulation. In Proceedings of the AlAA 13th Fluid and Plasma Dynamics Coference, Colorado, CO, USA, 14–16 July 1980.
- Liu, S.; Meneveau, C.; Katz, J. On the properties of similarity subgrid-scale models as deduced from measurements in a turbulent jet. J. Fluid Mech. 1994, 275, 83–119.
- Liu, S.; Meneveau, C.; Katz, J. Experimental study of similarity subgrid-scale models of turbulence in the far-field of a jet. Appl. Sci. Res. 1995, 54, 177–190.
- Akhavan, R.; Ansari, A.; Kang, S.; Mangiavacchi, N. Subgrid-scale interactions in a numerically simulated planar turbulent jet and implications for modelling. J. Fluid Mech. 2000, 408, 83–120.
- Domaradzki, J.A.; Saiki, E.M. A subgrid-scale model based on the estimation of unresolved scales of turbulence. Phys. Fluids. 1997, 9, 2148–2164.
- Domaradzki, J.A.; Loh, K.C. The subgrid-scale estimation model in the physical space representation. Phys. Fluids 1999, 11, 2330–2342.
- Stolz, S.; Adams, N.A. An approximate deconvolution procedure for large-eddy simulation. Phys. Fluids 1999, 11, 1699–1701.
- Stolz, S.; Adams, N.A.; Kleiser, L. An approximate deconvolution model for large-eddy simulation with application to incompressible wall-bounded flows. Phys. Fluids 2001, 13, 997–1015.
- Stolz, S.; Adams, N.A.; Kleiser, L. The approximate deconvolution model for large-eddy simulations of compressible flows and its application to shock-turbulent-boundary-layer interaction. Phys. Fluids 2001, 13, 2985–3001.
- Geurts, B.J.; Holm, D.D. Regularization modeling for large-eddy simulation. Phys. Fluids 2003, 15, 13–16.
- Cheskidov, A.; Holm, D.D.; Olson, E.; Titi, E.S. On a Leray-α model of turbulence. Proc. R. Soc. 2005, 461, 629–649.
- Cao, C.; Holm, D.D.; Titi, E.S. On the Clark-α model of turbulence: Global regularity and long-time dynamics. J. Turb. 2005, 6, 1–11.
- Stolz, S.; Adams, N.A. Large-eddy simulation of high-Reynolds-number supersonic boundary layers using the approximate deconvolution model and a rescaling and recycling technique. Phys. Fluids 2003, 15, 2398–2412.
- Schlatter, P.; Stolz, S.; Kleiser, L. LES of transitional flows using the approximate deconvolution model. Int. J. Heat Fluid Flow 2004, 25, 549–558.
- Saeedipour, M.; Vincent, S.; Pirker, S. Large eddy simulation of turbulent interfacial flows using approximate deconvolution model. Int. J. Multiph. Flow 2019, 112, 286–299.
- Ilyin, A.A.; Lunasin, E.M.; Titi, E.S. A modified-Leray-α subgrid scale model of turbulence. Nonlinearity 2006, 19, 879–897.
- Foias, C.; Holm, D.D.; Titi, E.S. The Navier-stokes-alpha model of fluid turbulence. Phys. D 2001, 152–153, 505–519.
- Rebholz, L.G. A family of new, high order NS-α models arising from helicity correction in Leray turbulence models. J. Math. Anal. Appl. 2008, 342, 246–254.
- Layton, W.; Rebholz, L.; Sussman, M. Energy and helicity dissipation rates of the NS-alpha and NS-alpha-deconvolution models. J. Appl. Math. 2010, 75, 56–74.
- Rebholz, L.G.; Sussman, M.M. On the high accuracy NS-alpha-deconvolution turbulence model. Math. Models Methods Appl. Sci. 2010, 20, 611–633.
- Breckling, S.; Neda, M. Numerical study of the Navier–Stokes-α deconvolution model with pointwise mass conservation. Int. J. Comput. Math. 2018, 95, 1727–1760.
- Cuff, V.M.; Dunca, A.A.; Monica, C.C.; Rebholz, L.G. The reduced order NS-α model for incompressible flow: Theory, numerical analysis and benchmark testing. Math. Modell. Numer. Anal. 2015, 49, 641–662.
- Bowers, A.L.; Rebholz, L.G. The reduced NS-α model for incompressible flow: A review of recent progress. Fluids 2017, 38, 38.
- Geurts, B.J.; Kuczaj, A.K.; Titi, E.S. Regularization modeling for large-eddy simulation of homogeneous isotropic decaying turbulence. J. Phys. Math. Theor. 2008, 41, 344008.
- Graham, J.P.; Holm, D.D.; Mininni, P.D.; Pouquet, A. Three regularization models of the navier-stokes equations. Phys. Fluids 2008, 20, 035107.
- Sheu, T.W.H.; Lin, Y.X.; Yu, C.H. Numerical study of two regularization models for simulating the turbulent flows. Comput. Fluids 2013, 74, 13–31.
- Huai, X.; Joslin, R.D.; Piomelli, U. Large-eddy simulation of transition to turbulence in boundary layers. Theor. Comput. Fluid Dyn. 1997, 9, 149–163.
- Yang, Z.; Voke, P.R. Large-eddy simulation of boundary-layer separation and transition at a change of surface curvature. J. Fluid Mech. 2001, 439, 305–333.
- Sayadi, T.; Moin, P. Large eddy simulation of controlled transition to turbulence. Phys. Fluids 2012, 24, 114103.
- Langari, M.; Yang, Z. Numerical study of the primary instability in a separated boundary layer transition under elevated free-stream turbulence. Phys. Fluids 2013, 25, 074106.
- Mittal, R.; Moin, P. Suitability of upwind-biased finite difference schemes for large-eddy simulation of turbulent flows. AIAA J. 1997, 35, 1415–1417.
- Kravchenko, A.G.; Moin, P. Numerical studies of flow over a circular cylinder at ReD = 3900. Phys. Fluids 2000, 12, 403–417.
- Dhotre, M.T.; Niceno, B.; Smith, B.L. Large eddy simulation of a bubble column using dynamic sub-grid scale model. Chem. Eng. J. 2008, 136, 337–348.
- Wang, Y.; Vanierschot, M.; Cao, L.; Cheng, Z.; Blanpain, B.; Guo, M. Hydrodynamics study of bubbly flow in a top-submerged lance vessel. Chem. Eng. Sci. 2018, 192, 1091–1104.
- Nilsen, K.M.; Kong, B.; Fox, R.O.; Hill, J.C.; Olsen, M.G. Effect of inlet conditions on the accuracy of large eddy simulations of a turbulent rectangular wake. Chem. Eng. J. 2014, 250, 175–189.
- Spille-Kohoff, A.; Kaltenbach, H.J. Generation of turbulent inflow data with a prescribed shear-stress profile. In Proceedings of the 3rd AFOSR International conference on DNS/LES, Arlington, TX, USA, 5–9 August 2001.
- Baba-Ahmadi, M.H.; Tabor, G.R. Inlet conditions for LES using mapping and feedback control. Comput. Fluids 2009, 38, 1299–1311.
- Baba-Ahmadi, M.H.; Tabor, G.R. Inlet conditions for large eddy simulation of gas-turbine swirl injectors. AIAA J. 2008, 46, 1782–1790.
- Aider, J.L.; Danet, A.; Lesieur, M. Large-eddy simulation applied to study the influence of upstream conditions on the time-dependant and averaged characteristics of a backward-facing step flow. J. Turb. 2007, 8, 1–30.
- Batten, P.; Goldberg, U.; Chakravarthy, S. Interfacing statistical turbulence closures with large-eddy simulation. AIAA J. 2004, 42, 485–492.
- Andersson, N.; Eriksson, L.E.; Davidson, L. LES prediction of flow and acoustic field of a coaxial jet. In Proceedings of the 11th AIAA/CES Aeroacoustics Conference, Monterey, CA, USA, 23–25 May 2005.
- Druault, P.; Lardeau, S.; Bonnet, J.P.; Coiffet, F.; Delville, J.; Lamballais, E.; Largeau, J.F.; Perret, L. Generation of three-dimensional turbulent inlet conditions for large-eddy simulation. AIAA J. 2004, 42, 447–456.
- Perret, L.; Delville, J.; Manceau, R.; Bonnet, J.P. Turbulent inflow conditions for large-eddy simulation based on low-order empirical model. Phys. Fluids 2008, 20, 075107.
- Klein, M.; Sadiki, A.; Janicka, J. A digital filter based generation of inflow data for spatially developing direct numerical or large eddy simulations. J. Comput. Phys. 2003, 186, 652–665.
- Di Mare, L.; Klein, M.; Jones, W.P.; Janicka, J. Synthetic turbulence inflow conditions for large-eddy simulation. Phys. Fluids 2006, 18, 025107.
- Jarrin, N.; Benhamadouche, S.; Laurence, D.; Prosser, R. A synthetic-eddy-method for generating inflow conditions for large-eddy simulations. Int. J. Heat Fluid Flow 2006, 27, 585–593.
- Jarrin, N.; Prosser, R.; Uribe, J.C.; Benhamadouche, S.; Laurence, D. Reconstruction of turbulent fluctuations for hybrid RANS/LES simulations using a synthetic-eddy method. Int. J. Heat Fluid Flow 2009, 30, 435–442.
- Skillen, A.; Revell, A.; Craft, T. Accuracy and efficiency improvements in synthetic eddy methods. Int. J. Heat Fluid Flow 2016, 62, 386–394.
- Keating, A.; Piomelli, U.; Balaras, E.; Kaltenbach, H.J. A priori and a posteriori tests of inflow conditions for large-eddy simulation. Phys. Fluids 2004, 16, 4696–4712.
- Benhamadouche, S.; Jarrin, N.; Addad, Y.; Laurence, D. Synthetic turbulent inflow conditions based on a vortex method for large-eddy simulation. Prog. Comput. Fluid Dyn. 2006, 6, 50–57.
- Kaltenbach, H.J.; Fatica, M.; Mittal, R.; Lund, T.S.; Moin, P. Study of flow in a planar asymmetric diffuser using large-eddy simulation. J. Fluid Mech. 1999, 390, 151–185.