Large-eddy simulation of circular cylinder flow at subcritical Reynolds number: Turbulent wake and sound radiation

The flows past a circular cylinder at Reynolds number 3900 are simulated using large-eddy simulation(LES) and the far-field sound is calculated from the LES results. A low dissipation energy-conserving finite volume scheme is used to discretize the incompressible Navier–Stokes equations. The dynamic global coefficient version of the Vreman's subgrid scale(SGS) model is used to compute the sub-grid stresses. Curle's integral of Lighthill's acoustic analogy is used to extract the sound radiated from the cylinder. The profiles of mean velocity and turbulent fluctuations obtained are consistent with the previous experimental and computational results. The sound radiation at far field exhibits the characteristic of a dipole and directivity. The sound spectra display the-5/3 power law. It is shown that Vreman's SGS model in company with dynamic procedure is suitable for LES of turbulence generated noise.     

Large eddy simulation of hot and cold fluids mixing in a T-junction for predicting thermal fluctuations

Temperature fluctuations in a mixing T-junction have been simulated on the FLUENT platform using the large eddy simulation （LES） turbulent flow model and a sub-grid scale Smagorinsky-Lilly model. The normalized mean and root mean square temperatures for describing time-averaged temperature and temperature fluctuation intensity, and the velocity are obtained. The power spectrum densities of temperature fluctuations, which are key parameters for thermal fatigue analysis and lifetime evaluation, are analyzed. Simulation results are consistent with experimental data published in the literature, showing that the LES is reliable. Several mixing processes under different conditions are simulated in order to analyze the effects of varying Reynolds number and Richardson number on the mixing course and thermal fluctuations.     

Multi-relaxation-time lattice Boltzmann simulations of lid driven flows using graphics processing unit

Large eddy simulation(LES) using the Smagorinsky eddy viscosity model is added to the two-dimensional nine velocity components(D2Q9) lattice Boltzmann equation(LBE) with multi-relaxation-time(MRT) to simulate incompressible turbulent cavity flows with the Reynolds numbers up to 1×10~7. To improve the computation efficiency of LBM on the numerical simulations of turbulent flows, the massively parallel computing power from a graphic processing unit(GPU) with a computing unified device architecture(CUDA) is introduced into the MRT-LBE-LES model. The model performs well, compared with the results from others, with an increase of 76 times in computation efficiency.It appears that the higher the Reynolds numbers is, the smaller the Smagorinsky constant should be, if the lattice number is fixed. Also, for a selected high Reynolds number and a selected proper Smagorinsky constant, there is a minimum requirement for the lattice number so that the Smagorinsky eddy viscosity will not be excessively large.     

Large eddy simulation of fuel injection and mixing process in a diesel engine

The large eddy simulation(LES) approach implemented in the KIVA-3V code and based on one-equation sub-grid turbulent kinetic energy model are employed for numerical computation of diesel sprays in a constant volume vessel and in a Caterpillar 3400 series diesel engine.Computational results are compared with those obtained by an RANS(RNG k-ε) model as well as with experimental data.The sensitivity of the LES results to mesh resolution is also discussed.The results show that LES generally provides flow and spray characteristics in better agreement with experimental data than RANS;and that small-scale random vortical structures of the in-cylinder turbulent spray field can be captured by LES.Furthermore,the penetrations of fuel droplets and vapors calculated by LES are larger than the RANS result,and the sub-grid turbulent kinetic energy and sub-grid turbulent viscosity provided by the LES model are evidently less than those calculated by the RANS model.Finally,it is found that the initial swirl significantly affects the spray penetration and the distribution of fuel vapor within the combustion chamber.     

Applications of URANS on predicting unsteady turbulent separated flows

Accurate prediction of unsteady separated turbulent flows remains one of the toughest tasks and a practi cal challenge for turbulence modeling. In this paper, a 2D flow past a circular cylinder at Reynolds number 3,900 is numerically investigated by using the technique of unsteady RANS （URANS）. Some typical linear and nonlinear eddy viscosity turbulence models （LEVM and NLEVM） and a quadratic explicit algebraic stress model （EASM） are evaluated. Numerical results have shown that a high-performance cubic NLEVM, such as CLS, are superior to the others in simulating turbulent separated flows with unsteady vortex shedding.     

Nonlinear turbulence models for predicting strong curvature effects

Prediction of the characteristics of turbulent flows with strong streamline curvature, such as flows in turbomachines, curved channel flows, flows around airfoils and buildings, is of great importance in engineering applications and poses a very practical challenge for turbulence modeling. In this paper, we analyze qualitatively the curvature effects on the structure of turbulence and conduct numerical simulations of a turbulent Uduct flow with a number of turbulence models in order to assess their overall performance. The models evaluated in this work are some typical linear eddy viscosity turbulence models, nonlinear eddy viscosity turbulence models （NLEVM） （quadratic and cubic）, a quadratic explicit algebraic stress model （EASM） and a Reynolds stress model （RSM） developed based on the second-moment closure. Our numerical results show that a cubic NLEVM that performs considerably well in other benchmark turbulent flows, such as the Craft, Launder and Suga model and the Huang and Ma model, is able to capture the major features of the highly curved turbulent U-duct flow, including the damping of turbulence near the convex wall, the enhancement of turbulence near the concave wall, and the subsequent turbulent flow separation. The predictions of the cubic models are quite close to that of the RSM, in relatively good agreement with the experimental data, which suggests that these models may be employed to simulate the turbulent curved flows in engineering applications.     

LES prediction of space-time correlations in turbulent shear flows

We compare the space-time correlations calculated from direct numerical simulation(DNS) and large-eddy simulation(LES) of turbulent channel flows.It is found from the comparisons that the LES with an eddy-viscosity subgrid scale(SGS) model over-predicts the space-time correlations than the DNS.The overpredictions are further quantified by the integral scales of directional correlations and convection velocities.A physical argument for the overprediction is provided that the eddy-viscosity SGS model alone does not includes the backscatter effects although it correctly represents the energy dissipations of SGS motions.This argument is confirmed by the recently developed elliptic model for space-time correlations in turbulent shear flows.It suggests that enstrophy is crucial to the LES prediction of spacetime correlations.The random forcing models and stochastic SGS models are proposed to overcome the overpredictions on space-time correlations.     

AN EXTENDED k-ε MODEL FOR NUMERICAL SIMULATION OF WIND FLOW AROUND BUILDINGS

It is assumed in this paper that for a high Reynolds number nearly homogeneouswind flow, the Reynolds stresses are uniquely related to the mean velocity gradientsand the two independent turbulent scaling parameters k and E. By applying dimensionalanalysis and owing to the Cayley-Hamilton theorem for tensors, a new turbulenceenclosure model so-called the axtended k-ε model has been developed. The coefficientsof the model expression were detemined by the wind tunnel experimental data ofhomogeneous shear turbulent flow. The model was compared with the standard k-εmodel in in composition and the prediction of the Reynold’s normal Stresses. Using thenew model the numerical simulation of wind flow around a square cross-section tallbuilding was performed. The results show that the extended k-ε model improves theprediction of wind velocities around the building the building and wind pressures on the buildingenvelope.     

Studies on gas turbulence and particle fluctuation in dense gas-particle flows

Particle fluctuation and gas turbulence in dense gas-particle flows are less studied due to complexity of the phenomena. In the present study, simulations of gas turbulent flows passing over a single particle are carried out first by using RANS modeling with a Reynolds stress equation turbulence model and sufficiently fine grids, and then by using LES. The turbulence enhancement by the particle wake effect is studied under various particle sizes and relative gas velocities, and the turbulence enhancement is found proportional to the particle diameter and the square of velocity. Based on the above results, a turbulence enhancement model for the particle-wake effect is proposed and is incorporated as a sub-model into a comprehensive two-phase flow model, which is then used to simulate dilute gas-particle flows in a horizontal channel. The simulation results show that the predicted gas turbulence by using the present model accounting for the particle wake effect is obviously in better agreement with the experimental results than the prediction given by the model not accounting for the wake effect. Finally, the proposed model is incorporated into another two-phase flow model to simulate dense gasparticle flows in a downer. The results show that the particle wake effect not only enhances the gas turbulence, but also amplifies the particle fluctuation.     

INVESTIGATION OF THE TURBULENT MODEL FOR PRESSURE FLUCTUATIONS WITH SPECTRAL THEORY

The pressure fluctuations in turbulent shear flows are investigated with the theory of spectral analysis.An expression for pressure spectra is analytically derived in terms of velocity spectra.This derivation is based on a formal solution of the Navier-Stokes equation and quasi-normal assumption to express the third and fourth order velocity correlations in terms of double velocity correlation.Then,a turbulent model for the computation of pressure fluctuation intensity with Renolds stress and mean flow velocity gradients is established.The turbulent constants in the model are calculated from the assumptions about the general behaviour of velocity spectra in high Renolds number flows.Comparison with direct simulation of turbulent boundary layer is made.It is found that the turbulent-turbulent,cross correlation,and turbulent-shear source terms for mean square value of pressure fluctuation are about the same magnitude.     

LES of incompressible heat and fluid flows past a square cylinder at high Reynolds numbers

This paper presents large eddy simulation (LES) results of incompressible heat and fluid flows around a square cylinder (SC) at zero incident angle at high Reynolds numbers (Re) in the range from 1.25×10⁵ to 3.5×10⁵. LES results are obtained on the basis of swirling strength based sub-grid model, and a higher order upwind scheme developed with respect to the Taylor expansion. It was found that, for the zero incident SC wake flows at a Reynolds number in the range {Re₅ = Re/10⁵ ∈ [1.25, 3.5]}, the Strouhal number equals to 0.1079, completely independent of the Reynolds number; the coefficient of drag is around 1.835 with an uncertainty of about 1.9%, almost non-sensitive to the Re. When Re is beyond 3.0×10⁵, the time-averaged peak value of sub-grid viscosity is over 340, implying that the role of sub-grid model is crucial in some regions where vortex motion is active and vortex interaction is intense. The time–spanwise (t-z) averaged sub-grid viscosity ratio profiles and the profiles of fluctuations of the sub-grid viscosity ratio and velocity components at four locations downstream of the SC are presented. The fields of the t-z averaged sub-grid viscosity ratio, and the instantaneous fields of streamwise and spanwise vorticities are also reported and discussed. The predicted mean Nusselt number is compared with empirical correlations, revealing that swirling strength based LES has its potential in predicting natural and industrial flows.     

Wall modeling for implicit large-eddy simulation and immersed-interface methods

We propose and analyze a wall model based on the turbulent boundary layer equations (TBLE) for implicit large-eddy simulation (LES) of high Reynolds number wall-bounded flows in conjunction with a conservative immersed-interface method for mapping complex boundaries onto Cartesian meshes. Both implicit subgrid-scale model and immersed-interface treatment of boundaries offer high computational efficiency for complex flow configurations. The wall model operates directly on the Cartesian computational mesh without the need for a dual boundary-conforming mesh. The combination of wall model and implicit LES is investigated in detail for turbulent channel flow at friction Reynolds numbers from Re _τ  = 395 up to Re _τ =100,000 on very coarse meshes. The TBLE wall model with implicit LES gives results of better quality than current explicit LES based on eddy viscosity subgrid-scale models with similar wall models. A straightforward formulation of the wall model performs well at moderately large Reynolds numbers. A logarithmic-layer mismatch, observed only at very large Reynolds numbers, is removed by introducing a new structure-based damping function. The performance of the overall approach is assessed for two generic configurations with flow separation: the backward-facing step at Re _h = 5,000 and the periodic hill at Re _H = 10,595 and Re _H = 37,000 on very coarse meshes. The results confirm the observations made for the channel flow with respect to the good prediction quality and indicate that the combination of implicit LES, immersed-interface method, and TBLE-based wall modeling is a viable approach for simulating complex aerodynamic flows at high Reynolds numbers. They also reflect the limitations of TBLE-based wall models.     

LES and RANS for Turbulent Flow over Arrays of Wall-Mounted Obstacles

Large-eddy simulation (LES) has been applied to calculate the turbulent flow over staggered wall-mounted cubes and staggered random arrays of obstacles with area density 25%, at Reynolds numbers between 5 × 10³ and 5 10⁶, based on the free stream velocity and the obstacle height. Re = 5 × 10³ data were intensively validated against direct numerical simulation (DNS) results at the same Re and experimental data obtained in a boundary layer developing over an identical roughness and at a rather higher Re. The results collectively confirm that Reynolds number dependency is very weak, principally because the surface drag is predominantly form drag and the turbulence production process is at scales comparable to the roughness element sizes. LES is thus able to simulate turbulent flow over the urban-like obstacles at high Re with grids that would be far too coarse for adequate computation of corresponding smooth-wall flows. Comparison between LES and steady Reynolds-averaged Navier-Stokes (RANS) results are included, emphasising that the latter are inadequate, especially within the canopy region.     

Vortices for computing: the engines of turbulence simulation

Vortices have been described as the “sinews of turbulence”. They are also, increasingly, the computational engines driving numerical simulations of turbulence. In this paper, I review some recent advances in vortex-based numerical methods for simulating high Reynolds number turbulent flows. I focus on coherent vortex simulation, where nonlinear wavelet filtering is used to identify and track the few high energy multiscale vortices that dominate the flow dynamics. This filtering drastically reduces the computational complexity for high Reynolds number simulations, e.g. by a factor of 1000 for fluid–structure interaction calculations (Kevlahan and Vasilyevvon in SIAM J Sci Comput 26(6):1894–1915, 2005). It also has the advantage of decomposing the flow into two physically important components: coherent vortices and background noise. In addition to its computational efficiency, this decomposition provides a way of directly estimating how space and space–time intermittency scales with Reynolds number, Re ^α . Comparing α to its non-intermittent values gives a realistic Reynolds number upper bound for adaptive direct numerical simulation of turbulent flows. This direct measure of intermittency also guides the development of new mathematical theories for the structure of high Reynolds number turbulence.     

Large Eddy Simulations Using an Unstructured Grid Compressible Navier - Stokes Algorithm

Large Eddy Simulation (LES) of the decay of isotropic turbulence and of channel flow has been performed using an explicit second-order unstructured grid algorithm for tetrahedral cells. The algorithm solves for cell-averaged values using the finite volume form of the unsteady compressible Jittered Navier-Stokes equations. The inviscid fluxes are obtained from Godunov's exact Riemann solver. Reconstruction of the flow variables to the left and right sides of each face is performed using least squares or Frink's method. The viscous fluxes and heat transfer are obtained by application of Gauss' theorem. LES of the decay of nearly incompressible isotropic turbulence has been performed using two models for the SGS stresses: the Monotone Integrated Large Eddy Simulation (MILES) approach, wherein the inherent numerical dissipation models the sub-grid scale (SGS) dissipation, and the Smagorinsky SGS model. The results using the MILES approach with least squares reconstruction show good agreement with incompressible experimental data. The contribution of the Smagorinsky SGS model is negligible. LES of turbulent channel flow was performed at a Reynolds number (based on channel height and bulk velocity) of 5600 and Mach number of 0.5 (at which compressibility effects are minimal) using Smagorinsky's SGS model with van Driest damping. The results show good agreement with experimental data and direct numerical simulations for incompressible channel flow. The SGS eddy viscosity is less than 10% of the molecular viscosity, and therefore the LES is effectively MILES with molecular viscosity.     

An improved anisotropy-resolving subgrid-scale model for flows in laminar–turbulent transition region

Some types of mixed subgrid-scale (SGS) models combining an isotropic eddy-viscosity model and a scale-similarity model can be used to effectively improve the accuracy of large eddy simulation (LES) in predicting wall turbulence. Abe (2013) has recently proposed a stabilized mixed model that maintains its computational stability through a unique procedure that prevents the energy transfer between the grid-scale (GS) and SGS components induced by the scale-similarity term. At the same time, since this model can successfully predict the anisotropy of the SGS stress, the predictive performance, particularly at coarse grid resolutions, is remarkably improved in comparison with other mixed models. However, since the stabilized anisotropy-resolving SGS model includes a transport equation of the SGS turbulence energy, k_SGS, containing a production term proportional to the square root of k_SGS, its applicability to flows with both laminar and turbulent regions is not so high. This is because such a production term causes k_SGS to self-reproduce. Consequently, the laminar–turbulent transition region predicted by this model depends on the inflow or initial condition of k_SGS. To resolve these issues, in the present study, the mixed-timescale (MTS) SGS model proposed by Inagaki et al. (2005) is introduced into the stabilized mixed model as the isotropic eddy-viscosity part and the production term in the k_SGS transport equation. In the MTS model, the SGS turbulence energy, k_es, estimated by filtering the instantaneous flow field is used. Since the k_es approaches zero by itself in the laminar flow region, the self-reproduction property brought about by using the conventional k_SGS transport equation model is eliminated in this modified model. Therefore, this modification is expected to enhance the applicability of the model to flows with both laminar and turbulent regions. The model performance is tested in plane channel flows with different Reynolds numbers and in a backward-facing step flow. The results demonstrate that the proposed model successfully predicts a parabolic velocity profile under laminar flow conditions and reduces the dependence on the grid resolution to the same degree as the unmodified model by Abe (2013) for turbulent flow conditions. Moreover, it is shown that the present model is effective at transitional Reynolds numbers. Furthermore, the present model successfully provides accurate results for the backward-facing step flow with various grid resolutions. Thus, the proposed model is considered to be a refined anisotropy-resolving SGS model applicable to laminar, transitional, and turbulent flows.     

Implicit large-eddy simulation applied to turbulent channel flow with periodic constrictions

The subgrid-scale (SGS) model in a large-eddy simulation (LES) operates on a range of scales which is marginally resolved by discretization schemes. Accordingly, the discretization scheme and the subgrid-scale model are linked. One can exploit this link by developing discretization methods from subgrid-scale models, or the converse. Approaches where SGS models and numerical discretizations are fully merged are called implicit LES (ILES). Recently, we have proposed a systematic framework for the design, analysis, and optimization of nonlinear discretization schemes for implicit LES. In this framework parameters inherent to the discretization scheme are determined in such a way that the numerical truncation error acts as a physically motivated SGS model. The resulting so-called adaptive local deconvolution method (ALDM) for implicit LES allows for reliable predictions of isotropic forced and decaying turbulence and of unbounded transitional flows for a wide range of Reynolds numbers. In the present paper, ALDM is evaluated for the separated flow through a channel with streamwise-periodic constrictions at two Reynolds numbers Re = 2,808 and Re = 10,595. We demonstrate that, although model parameters of ALDM have been determined for isotropic turbulence at infinite Reynolds number, it successfully predicts mean flow and turbulence statistics in the considered physically complex, anisotropic, and inhomogeneous flow regime. It is shown that the implicit model performs at least as well as an established explicit model.      

Parallel FEM LES with one-equation subgrid-scale model for incompressible flows

This article develops a parallel large-eddy simulation (LES) with a one-equation subgrid-scale (SGS) model based on the Galerkin finite element method and three-dimensional (3D) brick elements. The governing filtered Navier–Stokes equations were solved by a second-order accurate fractional-step method, which decomposed the implicit velocity–pressure coupling in incompressible flow and segregated the solution to the advection and diffusion terms. The transport equation for the SGS turbulent kinetic energy was solved to calculate the SGS processes. This FEM LES model was applied to study the turbulence of the benchmark open channel flow at a Reynolds number Re_τ = 180 (based on the friction velocity and channel height) using different model constants and grid resolutions. By comparing the turbulence statistics calculated by the current model with those obtained from direct numerical simulation (DNS) and experiments in literature, an optimum set of model constants for the current FEM LES model was established. The budgets of turbulent kinetic energy and vertical Reynolds stress were then analysed for the open channel flow. Finally, the flow structures were visualised to further reveal some important characteristics. It was demonstrated that the current model with the optimum model constants can predict well the organised structure near the wall and free surface, and can be further applied to other fundamental and engineering applications.