Turbulent flow around a wall-mounted cube: A direct numerical simulation

An immersed-boundary method was employed to perform a direct numerical simulation (DNS) of flow around a wall-mounted cube in a fully developed turbulent channel for a Reynolds number Re = 5610, based on the bulk velocity and the channel height. Instantaneous results of the DNS of a plain channel flow were used as a fully developed inflow condition for the main channel. The results confirm the unsteadiness of the considered flow caused by the unstable interaction of a horseshoe vortex formed in front of the cube and on both its sides with an arch-type vortex behind the cube. The time-averaged data of the turbulence mean-square intensities, Reynolds shear stresses, kinetic energy and dissipation rate are presented. The negative turbulence production is predicted in the region in front of the cube where the main horseshoe vortex originates.     

Large eddy simulation of flow around a rectangular cylinder

A large eddy simulation (LES) of turbulent flow around a stationary rectangular cylinder at high Reynolds number of 2.2×10⁴ is conducted as the first step to prove the applicability of LES to practical engineering problems. Time-averaged and phase-averaged velocities and turbulent stresses are obtained and they are compared with the experimental data. To investigate mesh dependence on computational results of the LES, two kinds of grid resolution are used. In addition, the effect of a second-order upwind scheme QUICK for convection terms is also investigated due to its dependence on grid resolution. The drag coefficients, the base pressure coefficients and Strouhal numbers are in fairly good agreement with the experimental results, while the computational results show that the artificial dissipative and dispersive effect of QUICK is large in the vicinity of the cylinder in our computation. Thus, it is necessary to use higher-order upwind schemes to reduce the numerical errors, since it is effective in applying LES to practical engineering problems with complicated geometry.     

Detached eddy simulation of cyclic large scale fluctuations in a simplified engine setup

Computational Fluid Dynamics using RANS-based modelling approaches have become an important tool in the internal combustion engine development and optimization process. However, these models cannot resolve cycle to cycle variations, which are an important aspect in the design of new combustion systems. In this study the feasibility of using a Detached Eddy Simulation (DES) SST model, which is a hybrid RANS/LES model, to predict cycle to cycle variations is investigated. In the near wall region or in regions where the grid resolution is not sufficiently fine to resolve smaller structures, the two-equation RANS SST model is used. In the other regions with higher grid resolution an LES model is applied. The case considered is a geometrically simplified engine, for which detailed experimental data for the ensemble averaged and single cycle velocity field are available from Boreé et al. [Boreé, J., Maurel, S., Bazile, R., 2002. Disruption of a compressed vortex, Physics of Fluids 14 (7), 2543–2556]. The fluid flow shows a strong tumbling motion, which is a major characteristic for modern turbo-charged, direct-injection gasoline engines. The general flow structure is analyzed first and the extent of the LES region and the amount of resolved fluctuations are discussed. Multiple consecutive cycles are computed and turbulent statistics of DES SST, URANS and the measured velocity field are compared for different piston positions. Cycle to cycle variations of the velocity field are analyzed for both computation and experiment with a special emphasis on the useability of the DES SST model to predict cyclic variations.     

Large eddy simulation of compressible channel flow

Astract The present study is a contribution to the analysis of wall-bounded compressible flows, including a special focus on wall modeling for compressible turbulent boundary layer in a plane channel. large eddy simulation (LES) of fully developed isothermal channel flows at Re = 3,000 and Re = 4,880 with a sufficient mesh refinement at the wall are carried out in the Mach number range 0.3 ≤ M ≤ 3 for two different source term formulations: first the classical extension of the incompressible configuration by Coleman et al. (J. Fluid Mech. 305:159–183, 1995), second a formulation presently derived to model both streamwise pressure drop and streamwise internal energy loss in a spatially developed compressible channel flow. It is shown that the second formulation is consistent with the spatial problem and yields a much stronger cooling effect at the wall than the classical formulation. Based on the present LES data bank, compressibility and low Reynolds number effects are analysed in terms of coherent structure and statistics. A study of the universality of the structure of the turbulence in non-hypersonic compressible boundary layers (M≤5) is performed in reference to Bradshaw (Annu. Rev. Fluid. Mech. 9:33–54, 1977). An improvement of the van Driest transformation is proposed; it accounts for both density and viscosity changes in the wall layer. Consistently, a new integral wall scaling (y ^c+) which accounts for strong temperature gradients at the wall is developed for the present non-adiabatic compressible flow. The modification of the strong Reynolds analogy proposed by Huang et al. (J. Fluid Mech. 305:185–218, 1995) to model the correlation between velocity and temperature for non-adiabatic wall layers is assessed on the basis of a Crocco–Busemann relation specific to channel flow. The key role of the mixing turbulent Prandtl number Pr _m is pointed out. Results show very good agreement for both source formulations although each of them involve a very different amount of energy transfer at the wall. The present work was performed within the framework of the French–German research initiative “large eddy simulation of complex flows’ (UR 507). The computing resources were provided by IDRIS-France. The authors gratefully acknowledge the financial support from the Centre National de la Recherche Scientifique (CNRS), the Centre d’été Mathématique de Recherche Avancée en Calcul Scientifique (CEMRACS) and the Direction Générale de l’Armement (DGA/D4S).     

Experimental study and large eddy simulation for the turbulent flow around four cylinders in an in-line square configuration

The turbulent flows around four cylinders in an in-line square configuration with different spacing ratios of 1.5, 2.5, 3.5 and 5.0 have been investigated experimentally at subcritical Reynolds numbers from 11,000 to 20,000. The mean and fluctuating velocity distributions were obtained using the laser Doppler anemometry (LDA) measurement. The digital particle image velocimetry (DPIV) was employed to characterize the full field vorticity and velocity distributions as well as other turbulent quantities. The experimental study indicated that several distinct flow patterns exist depending on the spacing ratio and subcritical Reynolds number for turbulent flow. The three-dimensional numerical simulations were also carried out using the large eddy simulation (LES) at Reynolds number of 15,000 with the spacing ratio of 1.5 and 3.5. The results show that the LES numerical predictions are in good agreement with the experimental measurements. Therefore, the three-dimensional vortex structures and the full field instantaneous and mean quantities of the flow field such as velocity field, vorticity field, etc., which are very difficult to obtain experimentally, can be extracted from the simulation results for the deepening of our understanding on the complex flow phenomena around four cylinders in in-line configuration.     

圆柱绕流流场结构的大涡模拟研究

为进一步揭示绕流现象的形成机理,本文分别对处于层流稳态区、尾流过渡区、剪切层转换区Re分别为26、200、1.4×105的三种典型流态下的单圆柱绕流进行了二维数值模拟研究。Re为26时应用层流模型直接求解N-S方程,而Re分别为200、1.4×105时使用大涡模拟的方法进行计算。数值模拟很好地再现了稳定的涡旋结构、周期性交替脱落的卡门涡街结构、不规则的涡旋结构,在此基础上分析了尾流结构的基本特征及其压强分布规律、平均的流场特性、积分参数(如升力系数、阻力系数、斯特劳哈尔等),并与有关研究成果进行了对比。研究发现,采用不同流动介质时流场特性有所差异,空气为介质时的计算结果更符合实验的成果,而水为介质时计算结果偏差较大,这主要是由尾流涡旋产生的不合理负压造成的。     

Numerical simulation of flow around a square cylinder in uniform-shear flow

This paper presents results obtained from a numerical simulation of a two-dimensional (2-D) incompressible linear shear flow over a square cylinder. Numerical simulations are performed, using the lattice Boltzmann method, in the ranges of 50⩽Re⩽200 and 0⩽K⩽0.5, where Re and K are the Reynolds number and the shear rate, respectively. The effect of the shear rate on the frequency of vortex shedding from the cylinder, and the lift and drag forces exerted on the cylinder are quantified together with the flow patterns around the cylinder. The present results show that vortex structure behind the cylinder is strongly dependant on both the shear rate and Reynolds number. When Re=50, a small K can disturb the steady state and cause an alternative vortex shedding with uneven intensity. In contrast, a large value of K will suppress the vortex shedding from the cylinder. When Re>50, the differences in the strength and size of vortices shed from the upper and lower sides of the cylinder become more pronounced as K increases. Vortex shedding disappears when K is larger than a critical value, which depends on Re. The flow patterns around the cylinder for different Re tend towards self-similarity with increasing K. The lift and drag forces exerted on the cylinder, in general, decrease with increasing K. Unlike a shear flow past a circular cylinder, the vortex shedding frequency past a square cylinder decreases with increasing the shear rate. A significant reduction of the drag force occurs in the range 0.15<K<0.3.     

Large eddy simulation of boundary layer flow under cnoidal waves

Water waves in coastal areas are generally nonlinear, exhibiting asymmetric velocity profiles with different amplitudes of crest and trough. The behaviors of the boundary layer under asymmetric waves are of great significance for sediment transport in natural circumstances. While previous studies have mainly focused on linear or symmetric waves, asymmetric wave-induced flows remain unclear, particularly in the flow regime with high Reynolds numbers.Taking cnoidal wave as a typical example of asymmetric waves, we propose to use an infinite immersed plate oscillating cnoidally in its own plane in quiescent water to simulate asymmetric wave boundary layer. A large eddy simulation approach with Smagorinsky subgrid model is adopted to investigate the flow characteristics of the boundary layer. It is verified that the model well reproduces experimental and theoretical results. Then a series of numerical experiments are carried out to study the boundary layer beneath cnoidal waves from laminar to fully developed turbulent regimes at high Reynolds numbers, larger than ever studied before.Results of velocity profile, wall shear stress, friction coefficient, phase lead between velocity and wall shear stress, and the boundary layer thickness are obtained. The dependencies of these boundary layer properties on the asymmetric degree and Reynolds number are discussed in detail.     

Large eddy simulation of fully developed turbulent flow in a curved channel

In this paper large eddy simulation of the fully developed turbulent flow in a curved channel is carried out. The computational results are presented and compared with the experimental results of Eskinazi and Yeh^[1]. It is shown that the numerical results of the present LES are reliable and the influence of the curvature on the turbulence feature is correctly revealed.     

Detached Eddy Simulation of Flow around A-Airfoil

In this study the flow around the Aérospatiale A-airfoil at maximum lift (α = 13.3°) and Re = 2 × 10⁶ is investigated by RANS and DES to quantify the influence of transient flow patterns on the quality of the flow prediction. This flow features a highly unsteady pressure-induced trailing-edge separation. The fairly high Reynolds number makes LES rather unattractive from a practical point of view as the numerical costs scale with ≈ Re^1.95 [9] and lead to an unreasonable effort to carry out these simulations. Like LES, DES is designed to capture the unsteady flow features and consequently is supposed to give results superior to RANS. The RANS solution captures the separation, but is notable to predict transient vortex shedding. The application of DES to a 2D domain leads to results comparable to those obtained by 2D-RANS. Only the 3D-DES is capable of predicting both the unsteady flow features and the development of the boundary layer on the airfoil. As expected, the spanwise resolution has a strong impact on the vortex shedding and on the size of the separation bubble. Attention must be paid to both accommodating the full extent of the flow structures and providing the adequate spatial resolution in order to accomplish satisfactory results – that compare favourably with LES – at only a fraction of computational cost.     

Large eddy simulation of free surface shallow‐water flow

Shallow‐water flow with free surface frequently occurs in ambient water bodies, in which the horizontal scale of motion is generally two orders of magnitude greater than the water depth. To accurately predict this flow phenomenon in more detail, a three‐dimensional numerical model incorporating the method of large eddy simulation (LES) has been developed and assessed. The governing equations are split into three parts in the finite difference solution: advection, dispersion and propagation. The advection part is solved by the QUICKEST scheme. The dispersion part is solved by the central difference method and the propagation part is solved implicitly using the Gauss–Seidel iteration method. The model has been applied to free surface channel flow for which ample experimental data are available for verification. The inflow boundary condition for turbulence is generated by a spectral line processor. The computed results compare favourably with the experimental data and those results obtained by using a periodic boundary condition. The performance of the model is also assessed for the case in which anisotropic grids and filters with horizontal grid size of the order of the water depth are used for computational efficiency. The coarse horizontal grid was found to cause a significant reduction in the large‐scale turbulent motion generated by the bottom turbulence, and the turbulent motion is predominately described by the sub‐grid scale (SGS) terms. The use of the Smagorinsky model for SGS turbulence in this situation is found inappropriate. A parabolic mixing length model, which accounts for the filtered turbulence, is then proposed. The new model can reproduce more accurately the flow quantities. Copyright © 2000 John Wiley & Sons, Ltd.     

A three-dimensional simulation of wind flow around buildings

The steady Navier–Stokes equation is solved to simulate the wind-flow environment of three-dimensional configurations of buildings. The method assumes an incident wind described by a power-law velocity profile. A new method for controlling the two-part nested solution iteration is introduced. The simulation is compared to some published wind-tunnel measurements.     

Large eddy simulation in a fully developed turbulent flow in a channel and comparison of subgrid eddy viscosity models

The accuracy and computational efficiency are compared for a number of models of subgrid eddy viscosity (Smagorinsky model, renormalization group model, and dynamic and one-parameter models). Space-filtered Navier-Stokes equations are solved numerically by the control-volume approach on a nonuniform grid with the use of high-resolution schemes in time and space. The numerical data are compared with the results of a physical experiment and direct numerical simulation. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 47, No. 3, pp. 31–42, May–June, 2006.     

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.     

Numerical simulation of viscous flow around unrestrained cylinders

A 2-D analysis is made for the dynamic interactions between viscous flow and one or more circular cylinders. The cylinder is free to respond to the fluid excitation and its motions are part of the solution. The numerical procedure is based on the finite volume discretization of the Navier–Stokes equations on adaptive tri-tree grids which are unstructured and nonorthogonal. Both a fully implicit scheme and a semi-implicit scheme in the time domain have been used for the momentum equations, while the pressure correction method based on the SIMPLE technique is adopted to satisfy the continuity equation. A new upwind method is developed for the triangular and unstructured mesh, which requires information only from two neighbouring cells but is of order of accuracy higher than linear. A new procedure is also introduced to deal with the nonorthogonal term. The pressure on the body surface required in solving the momentum equation is obtained through the Poisson equation in the local cell. Results including flow field, pressure distribution and force are provided for fixed single and multiple cylinders and for an unrestrained cylinder in steady incoming flow with Reynolds numbers at 200 and 500 and in unsteady flow with Keulegan–Carpenter numbers at 5 and 10.     

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.     

Local heat transfer around a wall-mounted cube at 45° to flow in a turbulent boundary layer

The flow and local heat transfer around a wall-mounted cube oriented 45° to the flow is investigated experimentally in the range of Reynolds number 4.2 × 10³–3.3 × 10⁴ based on the cube height. The distribution of local heat transfer on the cube and its base wall are examined, and it is clarified that the heat transfer distribution under the angled condition differs markedly to that for cube oriented perpendicular to the flow, particularly on the top face of the cube. The surface pressure distribution is also investigated, revealing a well-formed pair of leading-edge vortices extending from the front corner of the top face downstream along both front edges for Re>(1−2)×10⁴. Regions of high heat transfer and low pressure are formed along the flow reattachment and separation lines caused by these vortices. In particular, near the front corner of the top face, pressure suction and heat transfer enhancement are pronounced. The average heat transfer on the top face is enhanced at Re>(1−2)×10⁴ over that of a cube aligned perpendicular to the flow.     

Flow around two tandem square cylinders near a plane wall

An experimental study has been conducted to investigate the flow around two identical square cylinders in tandem arrangement and placed near a plane wall at a Reynolds number of 6,300. The inter-cylinder spacing ratio was varied from S ^* = 0.5 to 6, and the cylinder-to-wall gap ratio from G ^* = 0.25 to 2. Totally, 42 cases were considered to systematically examine the effects of wall proximity and the mutual interference between the two cylinders in the normalized gap–spacing (G ^*–S ^*) plane. The flow fields were captured using digital particle image velocimetry, in conjunction with measurements of the fluid forces (drag and lift) acting on the downstream cylinder using a piezoelectric load cell. The results show that the flow is highly dependent on the combined values of G ^* and S ^*. Categories relating to G ^* could be broadly classified as small-gap regime (G ^* < 0.5) at which periodic vortex shedding from the cylinders is suppressed, intermediate-gap regime (0.5 < G ^* < 1) where vortex shedding occurs but is under the influence of the wall proximity, and large-gap regime (G ^* > 1) where the wall effects become negligible. Similarly, the flow interference between the two cylinders can be divided into three basic categories as a function of S ^*, namely, shielding regime at S ^* < 1, reattachment regime at 1 < S ^* < 3, and impinging regime at S ^* > 3. Variations of force coefficients, amplitude spectra, Strouhal numbers, and Reynolds shear stress with G ^* and S ^* are presented to characterize the different flow regimes.     

Numerical simulation of hypersonic flow around a sharp cone

A hypersonic viscous flow in the shock layer near sharp cones is considered. The profiles of density and velocity, the slopes of the shock wave, and the pressures and heat fluxes calculated on the basis of the full viscous shock layer equations are compared with available experimental and theoretical data.