An efficient very large eddy simulation model for simulation of turbulent flow

Among the various hybrid methodologies, Speziale's very large eddy simulation (VLES) is one that was proposed very early. It is a unified simulation approach that can change seamlessly from Reynolds Averaged Navier–Stokes (RANS) to direct numerical simulation (DNS) depending on the numerical resolution. The present study proposes a new improved variant of the original VLES model. The advantages are achieved in two ways: (i) RANS simulation can be recovered near the wall which is similar to the detached eddy simulation concept; (ii) a LES subgrid scale model can be reached by the introduction of a third length scale, that is, the integral turbulence length scale. Thus, the new model can provide a proper LES mode between the RANS and DNS limits. This new methodology is implemented in the standard k ? ? model. Applications are conducted for the turbulent channel flow at Reynolds number of Re_τ = 395, periodic hill flow at Re = 10,595, and turbulent flow past a square cylinder at Re = 22,000. In comparison with the available experimental data, DNS or LES, the new VLES model produces better predictions than the original VLES model. Furthermore, it is demonstrated that the new method is quite efficient in resolving the large flow structures and can give satisfactory predictions on a coarse mesh. Copyright © 2012 John Wiley & Sons, Ltd.     

Simple lattice Boltzmann approach for turbulent buoyant flow simulation

In the present work, a simple large eddy simulation （LES）-based lattice Boltz- mann model （LBM） is developed for thermal turbulence research. This model is validated by some benchmark tests. The numerical results demonstrate the good performance of the present model for turbulent buoyant flow simulation.     

Large eddy simulation of bluff body wake in planar shear flow

Large eddy simulation of planar shear flow past a square cylinder has been investigated. Dynamic Smagorinsky model has been used to model subgrid scale stress. The shear parameter, K, namely the nondimensional streamwise velocity gradient in the lateral direction, is 0.0, 0.1 and 0.2. Reynolds number based on the centerline velocity is fixed at Re=21400. The time and span‐averaged velocity components, pressure coefficient, Reynolds stresses for uniform are in good agreement with the literature. In shear flow the calculated flow structure and mean velocity components are shown to be markedly different from those of the uniform flow. With increasing shear parameter, the cylinder wake is dominated by clockwise vortices. Both the velocity components in shear flow are compared with respective components in uniform flow. Comparison of normal and shear stresses between shear and no shear case have also been presented. Copyright © 2009 John Wiley & Sons, Ltd.     

Large eddy simulation of turbulent channel flow using an algebraic model

In this paper an algebraic model from the constitutive equations of the subgrid stresses has been developed. This model has an additional term in comparison with the mixed model, which represents the backscatter of energy explicitly. The proposed model thus provides independent modelling of the different energy transfer mechanisms, thereby capturing the effect of subgrid scales more accurately. The model is also found to depict the flow anisotropy better than the linear and mixed models. The energy transfer capability of the model is analysed for the isotropic decay and the forced isotropic turbulence. The turbulent plane channel flow simulation is performed over three Reynolds numbers, Re_τ=180, 395 and 590, and the results are compared with that of the dynamic model, Smagorinsky model, and the DNS data. Both the algebraic and dynamic models are in good agreement with the DNS data for the mean flow quantities. However, the algebraic model is found to be more accurate for the turbulence intensities and the higher‐order statistics. The capability of the algebraic model to represent backscatter is also demonstrated. Copyright © 2005 John Wiley & Sons, Ltd.     

Adaptive mesh refinement method-based large eddy simulation for the flow over circular cylinder at ReD = 3900

With the development of computational power, large eddy simulation (LES) method is increasingly used in simulating complex flow. However, there still exist many factors affecting the LES quality and appropriate mesh resolution is among one of them. This work aims to develop an automatic procedure to refine the LES mesh by combining adaptive mesh refinement (AMR) and LES quality criteria. An LES refinement criterion is developed by estimating the proper grid length scale which meets the accuracy requirement of LES method. With this criterion, the baseline mesh is automatically refined with the AMR method. In this work, an efficient one-shot refinement strategy is also proposed to reduce the overall simulation time. Current AMR-based LES method is verified with the typical LES test case about the flow past circular cylinder at Re _D = 3900. Results show that the automatically refined mesh provides systematically better agreement with experimental results and with current method the balance between accuracy and computational expense for LES can be obtained.     

A large eddy simulation of the near wake of a circular cylinder

Viscous flow around a circular cylinder at a subcritical Reynolds number is investigated using a large eddy simulation (LES) coupled with the Smagorinsky subgrid-scale (SGS) model. A fractional-step method with a second-order in time and a combined finite-difference/spectral approximations are used to solve the filtered three-dimensional incompressible Navier-Stokes equations. Calculations have been performed with and without the SGS model. Turbulence statistical behaviors and flow structures in the near wake of the cylinder are studied. Some calculated results, including the lift and drag coefficients, shedding frequency, peak Reynolds stresses, and time-average velocity profile, are in good agreement with the experimental and computational data, which shows that the Smagorinsky model can reasonably predict the global features of the flow and some turbulent statistical behaviors. The project supported by the National Science Fund for Distinguished Scholars (10125210), the Special Funds for Major State Basic Research Project (G1999032801) and the National Natural Science Foundation of China (19772062)     

Large eddy simulation and a simple wall model for turbulent flow calculation by a particle method

A turbulent channel flow and the flow around a cubic obstacle are calculated by the moving particle semi‐implicit method with the subparticle‐scale turbulent model and a wall model, which is based on the zero equation RANS (Reynolds Averaged Navier‐Stokes). The wall model is useful in practical problems that often involve high Reynolds numbers and wall turbulence, because it is difficult to keep high resolution in the near‐wall region in particle simulation. A turbulent channel flow is calculated by the present method to validate our wall model. The mean velocity distribution agrees with the log‐law velocity profile near the wall. Statistical values are also the same order and tendency as experimental results with emulating viscous layer by the wall model. We also investigated the influence of numerical oscillations on turbulence analysis in using the moving particle semi‐implicit method. Finally, the turbulent flow around a cubic obstacle is calculated by the present method to demonstrate capability of calculating practical turbulent flows. Three characteristic eddies appear in front of, over, and in the back of the cube both in our calculation and the experimental result that was obtained by Martinuzzi and Tropea. Mean velocity and turbulent intensity profiles are predicted in the same order and have similar tendency as the experimental result. Copyright © 2012 John Wiley & Sons, Ltd.     

An investigation of thermally stratified turbulent channel flow with temperature oscillation on the bottom wall by large eddy simulation

Thermally stratified shear turbulent channel flow with temperature oscillation on the bottom wall of the channel is calculated to investigate the behavior of turbulent flow and heat transfer by use of large eddy simulation (LES) approach coupled with dynamic subgrid-scale (SGS) models. The objective of this study is to deal with the effect of the temperature oscillation on turbulent behavior of thermally stratified turbulent channel flow and to examine the effectiveness of the LES technique for predicting statistically unsteady turbulent flow driven by time-varying buoyancy force. To validate the present calculation, thermally stratified shear turbulent channel flow is computed and compared with available data obtained by direct numerical simulation (DNS), which confirm that the present approach can be used to predict thermally stratified turbulent channel flow satisfactorily. Further, to illustrate the effect of the temperature oscillation with different Richardson numbers and periods of the oscillation on turbulence characteristics, the phase-averaged mean value and fluctuation of the resolved velocities and temperature, and instantaneous velocity fluctuation structures are analyzed.     

Passive heat transfer in a turbulent channel flow simulation using large eddy simulation based on the lattice Boltzmann method framework

In this paper, a large eddy simulation based on the lattice Boltzmann framework is carried out to simulate the heat transfer in a turbulent channel flow, in which the temperature can be regarded as a passive scalar. A double multiple relaxation time (DMRT) thermal lattice Boltzmann model is employed. While applying DMRT, a multiple relaxation time D3Q19 model is used to simulate the flow field, and a multiple relaxation time D3Q7 model is used to simulate the temperature field. The dynamic subgrid stress model, in which the turbulent eddy viscosity and the turbulent Prandtl number are dynamically computed, is integrated to describe the subgrid effect. Not only the strain rate but also the temperature gradient is calculated locally by the non-equilibrium moments. The Reynolds number based on the shear velocity and channel half height is 180. The molecular Prandtl numbers are set to be 0.025 and 0.71. Statistical quantities, such as the average velocity, average temperature, Reynolds stress, root mean square (RMS) velocity fluctuations, RMS temperature and turbulent heat flux are obtained and compared with the available data. The results demonstrate great reliability of DMRT–LES in studying turbulence.     

Large eddy simulation of turbulent flows via domain decomposition techniques. Part 2: applications

The present paper discusses the application of large eddy simulation to incompressible turbulent flows in complex geometries. Algorithmic developments concerning the flow solver were provided in the companion paper (Int. J. Numer. Meth. Fluids, 2003; submitted), which addressed the development and validation of a multi‐domain kernel suitable for the integration of the elliptic partial differential equations arising from the fractional step procedure applied to the incompressible Navier–Stokes equations. Numerical results for several test problems are compared to reference experimental and numerical data to demonstrate the potential of the method. Copyright © 2005 John Wiley & Sons, Ltd.     

Large eddy simulation in a turbulent channel flow using exit boundary conditions

The influence of the exit boundary conditions on the vanishing first derivative of the velocity components and constant pressure on the large eddy simulation of the fully developed turbulent channel flow has been investigated for equidistant and stretched grids at the channel exit. Results show that the chosen exit boundary conditions introduce some small disturbances that are mostly damped by the grid stretching. The difference of rms values between the fully developed turbulent channel flow with periodicity conditions and the fully developed channel flow using inlet and the exit boundary conditions is less than 10% for the equidistant grids and less than 5% for the stretched grids. The chosen boundary conditions are of interest because they may be used in complex problems with back flow at the exit. Copyright © 1999 John Wiley & Sons, Ltd.     

Numerical study of an oscillatory turbulent flow over a flat plate

Oscillatory turbulent flow over a flat plate is studied using large eddy simulation (LES) and Reynolds-average Navier-Stokes (RANS) methods. A dynamic subgrid-scale model is employed in LES and Saffman's turbulence model is used in RANS. The flow behaviors are discussed for the accelerating and decelerating phases during the oscillating cycle. The friction force on the wall and its phase shift from laminar to turbulent regime are also investigated for different Reynolds numbers. The project supported by the Youngster Funding of Academia Sinica and by the National Natural Science Foundation of China     

Large eddy simulation of droplet dispersion for inhomogeneous turbulent wall flow

A large eddy simulation (LES) coupled with a Lagrangian stochastic model has been applied to the study of droplet dispersion in a turbulent boundary layer. Droplets are tracked in a Lagrangian way. The velocity of the fluid particle along the droplet trajectory is considered to have a large-scale part and a small-scale part given by a modified three-dimensional Langevin model using the filtered subgrid scale (SGS) statistics. An appropriate Lagrangian correlation timescale is considered in order to include the influences of gravity and inertia. Two-way coupling is also taken into account. The inter-droplet collision has been introduced as the main mechanism of secondary breakup. A stochastic model for breakup has been generalized for coalescence simulation, thereby two phenomena, coalescence and breakup are simulated in the framework of a single stochastic model. The parameters of this model, selectively for coalescence and for breakup, are computed dynamically by relating them to the local resolved properties of the dispersed phase compared to the main fluid. The model is validated by comparison with an agglomeration model and with experimental results on secondary breakup. The LES coupled with Lagrangian particle tracking and the model for droplet coalescence and breakup is applied to the study of the atmospheric dispersion of wet cooling tower plumes. The simulations are done for different droplet size distributions and volume fractions. We focused on the influence of these parameters on mean concentration, concentration variance and mass flux profiles.     

Large eddy simulations of flow interference between two unequal sized square cylinders

A uniform flow past two unequal sized square cylinders arranged in a side-by-side pattern and at a Reynolds number of 50,000 has been investigated using large eddy simulation (LES) technique. The modelling of sub-grid scales of turbulence is done using the Smagorinsky model. The effect of the transverse gap ratio (T/D) on the flow characteristics has been studied. Numerical simulations are carried out for five different transverse gap ratios (T/D), namely 1.120, 1.250, 1.375, 1.750 and 2.500. Results in terms of the aerodynamic forces, Strouhal number, mean base pressure coefficient, streamlines, vorticity, surface pressure distribution, normal and shear stresses are presented. A shift in the stagnation point for the small square cylinder from the centre of its front face towards its gap side is seen at smaller T/D ratios. The presence of a jet-like flow seen in the gap side is more pronounced at T/D = 1.12. A biased gap side flow towards the near wake of the small square cylinder is seen at smaller T/D ratios. No interference effect is seen at T/D = 2.5. The flow behaviour is similar to that of the isolated square cylinder at this gap ratio.     

方形钝体受限绕流的三维数值模拟

采用一种具有二阶精度的分裂步有限元方法作为大涡模拟的空间离散格式,经过标准算例的验证后,对Re=1.0×104条件下的方形钝体三维受限湍流绕流流场进行了数值模拟.计算中,为消除初始效应,略去初始段的计算结果.数值分析表明在均匀来流条件下,该湍流场沿槽道轴面对称,并呈现出一定的拟周期特性.在流场特性分析的基础上,进行了湍流能耗场的分析,结果表明,方形钝体受限绕流的能耗主要集中在大涡丰富的流动区段内.计算过程反映出,采用该空间离散格式的大涡模拟方法,能够捕捉到非常丰富的涡系及涡动的时变过程,适用于方形钝体受限绕流的数值模拟.     

Extension of the local domain-free discretization method to large eddy simulation of turbulent flows

In this work, an immersed boundary method, called the local domain-free discretization (DFD) method, is extended to large eddy simulation (LES) of turbulent flows. The discrete form of partial differential equations at an interior node may involve some nodes outside the solution domain. The flow variables at these exterior dependent nodes are evaluated via linear extrapolation along the direction normal to the wall. To alleviate the requirement of mesh resolution in the near-wall region, a wall model based on the turbulence boundary layer equations is introduced. The wall shear stress yielded by the wall model and the no-penetration condition are enforced at the immersed boundary to evaluate the velocity components at an exterior dependent node. For turbulence closure, a dynamic subgrid scale (SGS) model is adopted and the Lagrangian averaging procedure is used to compute the model coefficient. The SGS eddy viscosity at an exterior dependent node is set to be equal to that at the outer layer. To maintain the mass conservation near the immersed boundary, a mass source/sink term is added into the continuity equation. Numerical experiments on relatively coarse meshes with stationary or moving solid boundaries have been conducted to verify the ability of the present LES-DFD method. The predicted results agree well with the published experimental or numerical data.