Validation of a novel very large eddy simulation method for simulation of turbulent separated flow

The paper describes the validation of a newly developed very LES (VLES) method for the simulation of turbulent separated flow. The new VLES method is a unified simulation approach that can change seamlessly from Reynolds‐averaged Navier–Stokes to DNS depending on the numerical resolution. Four complex test cases are selected to validate the performance of the new method, that is, the flow past a square cylinder at Re = 3000 confined in a channel (with a blockage ratio of 20%), the turbulent flow over a circular cylinder at Re = 3900 as well as Re = 140,000, and a turbulent backward‐facing step flow with a thick incoming boundary layer at Re = 40,000. The simulation results are compared with available experimental, LES, and detached eddy simulation‐type results. The new VLES model performs well overall, and the predictions are satisfactory compared with previous experimental and numerical results. It is observed that the new VLES method is quite efficient for the turbulent flow simulations; that is, good predictions can be obtained using a quite coarse mesh compared with the previous LES method. Discussions of the implementation of the present VLES modeling are also conducted on the basis of the simulations of turbulent channel flow up to high Reynolds number of Re_τ = 4000. The efficiency of the present VLES modeling is also observed in the channel flow simulation. From a practical point of view, this new method has considerable potential for more complex turbulent flow simulations at relative high Reynolds numbers. Copyright © 2013 John Wiley & Sons, Ltd.     

DNS of a spatially developing turbulent boundary layer with passive scalar transport

A direct numerical simulation (DNS) of a spatially developing turbulent boundary layer over a flat plate under zero pressure gradient (ZPG) has been carried out. The evolution of several passive scalars with both isoscalar and isoflux wall boundary condition are computed during the simulation. The Navier–Stokes equations as well as the scalar transport equation are solved using a fully spectral method. The highest Reynolds number based on the free-stream velocity U_∞ and momentum thickness θ is Re_θ=830, and the molecular Prandtl numbers are 0.2, 0.71 and 2. To the authors’ knowledge, this Reynolds number is to date the highest with such a variety of scalars. A large number of turbulence statistics for both flow and scalar fields are obtained and compared when possible to existing experimental and numerical simulations at comparable Reynolds number. The main focus of the present paper is on the statistical behaviour of the scalars in the outer region of the boundary layer, distinctly different from the channel-flow simulations. Agreements as well as discrepancies are discussed while the influence of the molecular Prandtl number and wall boundary conditions is also highlighted. A Pr scaling for various quantities is proposed in outer scalings. In addition, spanwise two-point correlation and instantaneous fields are employed to investigate the near-wall streak spacing and the coherence between the velocity and the scalar fields. Probability density functions (PDF) and joint probability density functions (JPDF) are shown to identify the intermittency both near the wall and in the outer region of the boundary layer. The present simulation data will be available online for the research community.     

LES of temporally evolving mixing layers by an eighth‐order filter scheme

An eighth‐order filter method for a wide range of compressible flow speeds (H. C. Yee and B. Sjogreen, Proceedings of ICOSAHOM09, June 22–26, 2009, Trondheim, Norway) is employed for large eddy simulations (LES) of temporally evolving mixing layers (TML) for different convective Mach numbers (M_c) and Reynolds numbers. The high‐order filter method is designed for accurate and efficient simulations of shock‐free compressible turbulence, turbulence with shocklets, and turbulence with strong shocks with minimum tuning of scheme parameters. The value of the M_c considered is for the TML range from the quasi‐incompressible regime to the highly compressible supersonic regime. The three main characteristics of compressible TML (the self‐similarity property, compressibility effects, and the presence of large‐scale structures with shocklets for high M_c) are considered for the LES study. The LES results that used the same scheme parameters for all studied cases agree well with experimental results and published direct numerical simulations (DNS). Published 2012. This article is a US Government work and is in the public domain in the USA.     

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.     

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.     

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.     

Turbulent flow past a bubble and an ellipsoid using shadow-image and PIV techniques

An experimental investigation on flow around an oscillating bubble and solid ellipsoid with a flat bottom was conducted. A single air bubble (equivalent diameter D_e=9.12 mm) was attached to a small disk (1 mm) at the end of a needle and suspended across a vertical square channel (100 mm) by wire wherein water flowed downward at a constant flowrate. The solid ellipsoid (D_e9.1 mm) was suspended across the square channel in the same manner. The equivalent diameter-based Reynolds and Eotvos number range, 1950<Re<2250 and 11<Eo<11.5, placed the bubble in the ‘wobbly’ regime while the flow in its wake was turbulent. A constant flowrate and one bubble size was used such that flow in the wake was turbulent. Velocity measurements of the flow field around the bubble or solid were made using a one CCD camera Digital Particle Image Velocimetry (DPIV) system enhanced by Laser Induced Fluorescence (LIF). The shape of the bubble or solid was simultaneously recorded along with the velocity using a second CCD camera and an Infrared Shadow Technique (IST). In this way both the flow-field and the boundary of the bubble (solid) were measured. The velocity vector plots of flow around and in the wake of a bubble/solid, supplemented by profiles and contours of the average and root-mean-square velocities, vorticity, Reynolds stress and turbulent kinetic energy, revealed differences in the wake flow structure behind a bubble and solid. One of the significant differences was in the inherent, oscillatory motion of the bubble which not only produced vorticity in the near-wake, but as a result of apparent vorticity stretching distributed the turbulent kinetic energy associated with this flow more uniformly on its wake, in contrast to the solid.     

Progress in subgrid-scale transport modelling for continuous hybrid non-zonal RANS/LES simulations

The partially integrated transport modelling (PITM) method can be viewed as a continuous approach for hybrid RANS/LES modelling allowing seamless coupling between the RANS and the LES regions. The subgrid turbulence quantities are thus calculated from spectral equations depending on the varying spectral cutoff location [Schiestel, R., Dejoan, A., 2005. Towards a new partially integrated transport model for coarse grid and unsteady turbulent flow simulations. Theoretical and Computational Fluid Dynamics 18, 443–468; Chaouat, B., Schiestel, R., 2005. A new partially integrated transport model for subgrid-scale stresses and dissipation rate for turbulent developing flows. Physics of Fluids, 17 (6)] The PITM method can be applied to almost all statistical models to derive its hybrid LES counterpart. In the present work, the PITM version based on the transport equations for the turbulent Reynolds stresses together with the dissipation transport rate equation is now developed in a general formulation based on a new accurate energy spectrum function E(κ) valid in both large and small eddy ranges that allows to calibrate more precisely the c_sgs2 function involved in the subgrid dissipation rate _sgs transport equation. The model is also proposed here in an extended form which remains valid in low Reynolds number turbulent flows. This is achieved by considering a characteristic turbulence length-scale based on the total turbulent energy and the total dissipation rate taking into account the subgrid and resolved parts of the dissipation rate. These improvements allow to consider a large range of flows including various free flows as well as bounded flows. The present model is first tested on the decay of homogeneous isotropic turbulence by referring to the well known experiment of Comte-Bellot and Corrsin. Then, initial perturbed spectra E(κ) with a peak or a defect of energy are considered for analysing the model capabilities in strong non-equilibrium flow situations. The second test case is the classical fully turbulent channel flow that allows to assess the performance of the model in non-homogeneous flows characterised by important anisotropy effects. Different simulations are performed on coarse and refined meshes for checking the grid independence of solutions as well as the consistency of the subgrid-scale model when the filter width is changed. A special attention is devoted to the sharing out of the energy between the subgrid-scales and the resolved scales. Both the mean velocity and the turbulent stress computations are compared with data from direct numerical simulations.     

Numerical prediction of turbulent boundary layer noise from a sharp-edged flat plate

An efficient hybrid uncorrelated wall plane waves–boundary element method (UWPW-BEM) technique is proposed to predict the flow-induced noise from a structure in low Mach number turbulent flow. Reynolds-averaged Navier-Stokes equations are used to estimate the turbulent boundary layer parameters such as convective velocity, boundary layer thickness, and wall shear stress over the surface of the structure. The spectrum of the wall pressure fluctuations is evaluated from the turbulent boundary layer parameters and by using semi-empirical models from literature. The wall pressure field underneath the turbulent boundary layer is synthesized by realizations of uncorrelated wall plane waves (UWPW). An acoustic BEM solver is then employed to compute the acoustic pressure scattered by the structure from the synthesized wall pressure field. Finally, the acoustic response of the structure in turbulent flow is obtained as an ensemble average of the acoustic pressures due to all realizations of uncorrelated plane waves. To demonstrate the hybrid UWPW-BEM approach, the self-noise generated by a flat plate in turbulent flow with Reynolds number based on chord Re_c = 4.9 × 10⁵ is predicted. The results are compared with those obtained from a large eddy simulation (LES)-BEM technique as well as with experimental data from literature.     

Direct numerical simulation of turbulent heat transfer in annuli: Effect of heat flux ratio

Fully developed turbulent flow and heat transfer in a concentric annular duct is investigated for the first time by using a direct numerical simulation (DNS) with isoflux conditions imposed at both walls. The Reynolds number based on the half-width between inner and outer walls, δ=(r₂-r₁)/2, and the laminar maximum velocity is Re_δ=3500. A Prandtl number Pr=0.71 and a radius ratio r^*=0.1 were retained. The main objective of this work is to examine the effect of the heat flux density ratio, q^*=q₁/q₂, on different thermal statistics (mean temperature profiles, root mean square (rms) of temperature fluctuations, turbulent heat fluxes, heat transfer, etc.). To validate the present DNS calculations, predictions of the flow and thermal fields with q^*=1 are compared to results recently reported in the archival literature. A good agreement with available DNS data is shown. The effect of heat flux ratio q^* on turbulent thermal statistics in annular duct with arbitrarily prescribed heat flux is discussed then. This investigation highlights that heat flux ratio has a marked influence on the thermal field. When q^* varies from 0 to 0.01, the rms of temperature fluctuations and the turbulent heat fluxes are more intense near the outer wall while changes in q^* from 1 to 100, lead to opposite trends.     

Performance of a 12-sensor vorticity probe in the near field of a rectangular turbulent jet

The design and operational characteristics of a 12-sensor hot wire probe for three-dimensional velocity–vorticity measurements in turbulent flow fields is described and discussed. The performance of the probe is investigated in comparison with X-sensor probe measurements in the near field of a rectangular turbulent jet with aspect ratio 6. Measurements have been conducted at Reynolds number Re _D = 21,000 at nozzle distances of x/D = 1, 3, 6 and 11, where D is the width of the nozzle. The results obtained with the 12-sensor probe compare well to the results of the X-sensor probe. Distributions of mean and fluctuating velocity–vorticity fields are presented and discussed. Among the results the most prominent is the experimental confirmation of the high levels of fluctuating vorticity in the shear layers.     

Effect of particles on turbulent thermal field of channel flow with different Prandtl numbers

The direct numerical simulation(DNS) of heat transfer in a fully developed non-isothermal particle-laden turbulent channel flow is performed.The focus of this paper is on the modulation of the particles on turbulent thermal statistics in the particle-laden flow with three Prandtl numbers(P r = 0.71,1.5,and 3.0) and a shear Reynolds number(Reτ = 180).Some typical thermal statistics,including normalized mean temperature and their fluctuations,turbulent heat fluxes,Nusselt number and so on,are analyzed.The results show that the particles have less effects on turbulent thermal fields with the increase of Prandtl number.Two reasons can explain this.First,the correlation between fluid thermal field and velocity field decreases as the Prandtl number increases,and the modulation of turbulent velocity field induced by the particles has less influence on the turbulent thermal field.Second,the heat exchange between turbulence and particles decreases for the particle-laden flow with the larger Prandtl number,and the thermal feedback of the particles to turbulence becomes weak.     

The stability of low Reynolds number round jets

Multigrid cross-correlation digital particle image velocimetry (MCCDPIV) is used to investigate the stability and structure of low Reynolds number axisymmetric jets. The in-plane velocities, out-of-plane vorticity and some of the components of the Reynolds stress tensor are measured. Two Reynolds numbers based on the orifice outlet diameter are examined (680 and 1,030) at two different positions: one close to the orifice, ranging from 2D ₀ to 5D ₀ (D ₀ is the orifice diameter); and the other further from the orifice, ranging from 10D ₀ to 14.4D ₀. The results show that the lower Reynolds number jet (Re=680) is marginally unstable in the near-orifice region and is best described as laminar. Further downstream some intermittent structures are observed in the jet, and the growth in integrated turbulent kinetic energy with axial position indicates that the jet is also unstable in this region. For the higher Reynolds number jet (Re=1,030) the increasing size and intensity of vortical structures in the jet in the near-orifice region observed from the MCCDPIV data and the growth in integrated turbulent kinetic energy indicate that the jet is unstable. Further downstream this jet is best described as transitional or turbulent. From flow visualisation images in the near-orifice region it seems that, for both Reynolds numbers, shear layer roll-up occurs when the jet exits the orifice and enters the quiescent fluid in the tank, resulting in vortical structures that appear to grow as the jet proceeds. This is indicative of instability in both cases and is consistent with previous flow visualisation studies of low Reynolds number round jets. Discrepancies observed between the flow visualisation results and the MCCDPIV data is addressed. On the basis of flow visualisation results it is generally assumed that round jets are unstable at very low Reynolds number, however the present work shows that this assertion may be incorrect.     

DNS of backward‐facing step flow with fully turbulent inflow

Direct numerical simulation (DNS) has been performed to study the channel flow over a backward‐facing step at a Reynolds number Re_b=5600 based on the step height h and the inflow bulk velocity U_b. A dynamic method has been used in order to generate realistic turbulent inflow conditions. The results upstream of the step compared well with the fully developed channel flow. Downstream of the step our results show excellent agreement with experimental data. Copyright © 2009 John Wiley & Sons, Ltd.     

DNS and LES of turbulent flow in a closed channel featuring a pattern of hemispherical roughness elements

Direct Numerical Simulations (DNS) and Large Eddy Simulations (LES) were performed for fully-developed turbulent flow in channels with smooth walls and walls featuring hemispherical roughness elements at shear Reynolds numbers Re_τ = 180 and 400, with the goal of studying the effect of these roughness elements on the wall-layer structure and on the friction factor. The LES and DNS approaches were verified first by comparison with existing DNS databases for smooth walls. Then, a parametric study for the hemispherical roughness elements was conducted, including the effects of shear Reynolds number, normalized roughness height (k⁺ = 10–20) and relative roughness spacing (s⁺/k⁺ = 2–6). The sensitivity study also included the effect of distribution pattern (regular square lattice vs. random pattern) of the roughness elements on the walls. The hemispherical roughness elements generate turbulence, thus increasing the friction factor with respect to the smooth-wall case, and causing a downward shift in the mean velocity profiles. The simulations revealed that the friction factor decreases with increasing Reynolds number and roughness spacing, and increases strongly with increasing roughness height. The effect of random element distribution on friction factor and mean velocities is however weak. In all cases, there is a clear cut between the inner layer near the wall, which is affected by the presence of the roughness elements, and the outer layer, which remains relatively unaffected. The study reveals that the presence of roughness elements of this shape promotes locally the instantaneous flow motion in the lateral direction in the wall layer, causing a transfer of energy from the streamwise Reynolds stress to the lateral component. The study indicates also that the coherent structures developing in the wall layer are rather similar to the smooth case but are lifted up by almost a constant wall-unit shift y⁺ (∼10–15), which, interestingly, corresponds to the relative roughness k⁺ = 10.     

DNS and scaling laws from new symmetry groups of ZPG turbulent boundary layer flow

The Lie group, or symmetry approach, developed by Oberlack (see e.g. Oberlack [26] and references therein) is used to derive new scaling laws for various quantities of a zero pressure gradient turbulent boundary layer flow. The approach unifies and extends the work done by Oberlack for the mean velocity of stationary parallel turbulent shear flows. From the two-point correlation (TPC) equations the knowledge of the symmetries allows us to derive a variety of invariant solutions (scaling laws) for turbulent flows, one of which is the new exponential mean velocity profile that is found in the mid-wake region of flat-plate boundary layers. Further, a third scaling group was found in the TPC equations for the one-dimensional turbulent boundary layer. This is in contrast to the Navier–Stokes and Euler equations, which have one and two scaling groups, respectively. The present focus is on the exponential law in the outer region of turbulent boundary layer corresponding new scaling laws for one- and two-point correlation functions. A direct numerical simulation (DNS) of a flat plate turbulent boundary layer with zero pressure gradient was performed at two different Reynolds numbers Re=750,2240. The Navier–Stokes equations were numerically solved using a spectral method with up to 140 million grid points. The results of the numerical simulations are compared with the new scaling laws. TPC functions are presented. The numerical simulation shows good agreement with the theoretical results, however only for a limited range of applicability. PACS 02.20.-a, 47.11.+j, 47.27.Nz, 47.27.Eq     

High‐order discontinuous Galerkin spectral element methods for transitional and turbulent flow simulations

In this paper, we investigate the accuracy and efficiency of discontinuous Galerkin spectral method simulations of under‐resolved transitional and turbulent flows at moderate Reynolds numbers, where the accurate prediction of closely coupled laminar regions, transition and developed turbulence presents a great challenge to large eddy simulation modelling. We take full advantage of the low numerical errors and associated superior scale resolving capabilities of high‐order spectral methods by using high‐order ansatz functions up to 12th order. We employ polynomial de‐aliasing techniques to prevent instabilities arising from inexact quadrature of nonlinearities. Without the need for any additional filtering, explicit or implicit modelling, or artificial dissipation, our high‐order schemes capture the turbulent flow at the considered Reynolds number range very well. Three classical large eddy simulation benchmark problems are considered: a circular cylinder flow at Re_D=3900, a confined periodic hill flow at Re_h=2800 and the transitional flow over a SD7003 airfoil at Re_c=60,000. For all computations, the total number of degrees of freedom used for the discontinuous Galerkin spectral method simulations is chosen to be equal or considerably less than the reported data in literature. In all three cases, we achieve an equal or better match to direct numerical simulation results, compared with other schemes of lower order with explicitly or implicitly added subgrid scale models. Copyright © 2014 John Wiley & Sons, Ltd.     

Finite element study of free surface turbulent flow in rotating channels

Under certain conditions of liquid flow through rotating channels, the Coriolis force can induce a free surface to be formed. This problem is of practical importance in a Coriolis wear tester, which is used for determining the sliding wear coefficient of wear materials in slurry handling equipment. A deforming Galerkin finite element method is presented for predicting two‐dimensional turbulent free surface mean flow in rotating channels. Reynolds‐averaged Navier–Stokes (RANS) equations are cast into weak(algebraic) form using primitive variables (velocity and pressure). Eddy viscosity is determined via a mixing length model. Velocity is interpolated biquadratically, while pressure is interpolated bilinearly. The kinematic condition is used to form the Galerkin residual for the free surface. The free surface is represented by Hermite polynomials of zeroeth order for continuity of position and slope. Combined Newton's iteration is used to simultaneously solve for the free surface and the field variables. Results of velocity and pressure fields, as well as the free surface are shown to converge with mesh‐size refinement. There is excellent respect for mass conservation. Results are presented for various values of Rossby number (Ro) and height‐based Reynolds number (Re_H). Parameter continuation in Ro and Re_H space is used to compute solutions at higher values of flow rate and angular velocity. Copyright © 2002 John Wiley & Sons, Ltd.     

Turbulent channel flow with 2D wedges of random height on one wall

Direct numerical simulations (DNS) of a turbulent channel flow with 2D wedges of random height on the bottom wall have been performed. In addition, two other simulations have been carried out to assess the effect of the geometry on the overlying flow. In the first simulation, the four smallest elements were removed while in the other, a uniform distribution of wedges with the same area was used. Two Reynolds numbers were studied, Re_b=2500 and Re_b=5000 which correspond in case of smooth walls to Re_τ=180 and 300, respectively. Roughness on the wall induces separated regions, the reattachment occurring on the walls of the wedges or on the bottom wall. The pressure gradients on the walls increase the ejections and inrushes towards the wall. As a consequence the flow is more isotropic. The mechanism inducing an improved isotropy has been explained in term of the spectra and budgets of Reynolds stress. The comparison of the 3 surfaces has shown that near the wall, the uniformly distributed roughness represents only a poor approximation of the surface with wedges of random height. The Reynolds stresses, pressure distribution and spectra on the modified wall agree well with those on the random surface. Energy spectra show the pitch to height ratio of the largest elements to be the more appropriate geometrical parameter to describe the geometry.     

On the role of the smallest scales of a passive scalar field in a near-wall turbulent flow

Role of the smallest diffusive scales of a passive scalar field in the near-wall turbulent flow was examined with pseudo-spectral numerical simulations. Temperature fields were analyzed at friction Reynolds number Re _τ=171 and at Prandtl numbers, Pr=1 and Pr=5.4. Results of direct numerical simulations (DNS) were compared with the under-resolved simulations where the velocity field was still resolved with the DNS accuracy, while a coarser grid was used to describe the temperature fields. Since the smallest temperature scales remained unresolved in these simulations, an appropriate spectral turbulent thermal diffusivity was applied to avoid pile-up at the higher wave numbers. In spite of coarser numerical grids, the temperature fields are still highly correlated with the DNS results, including instantaneous temperature fields. Results point to practically negligible role of the diffusive temperature scales on the macroscopic behavior of the turbulent heat transfer.