Hybrid RANS/LES computations of plane impinging jets with DES and PANS models

The qualities of a DES (Detached Eddy Simulation) and a PANS (Partially-Averaged Navier–Stokes) hybrid RANS/LES model, both based on the k–ω RANS turbulence model of Wilcox (2008, “Formulation of the k–ω turbulence model revisited” AIAA J., 46: 2823–2838), are analysed for simulation of plane impinging jets at a high nozzle-plate distance (H/B = 10, Re = 13,500; H is nozzle-plate distance, B is slot width; Reynolds number based on slot width and maximum velocity at nozzle exit) and a low nozzle-plate distance (H/B = 4, Re = 20,000). The mean velocity field, fluctuating velocity components, Reynolds stresses and skin friction at the impingement plate are compared with experimental data and LES (Large Eddy Simulation) results. The k–ω DES model is a double substitution type, following Davidson and Peng (2003, “Hybrid LES–RANS modelling: a one-equation SGS model combined with a k–ω model for predicting recirculating flows” Int. J. Numer. Meth. Fluids, 43: 1003–1018). This means that the turbulent length scale is replaced by the grid size in the destruction term of the k-equation and in the eddy viscosity formula. The k–ω PANS model is derived following Girimaji (2006, “Partially-Averaged Navier–Stokes model for turbulence: a Reynolds-Averaged Navier–Stokes to Direct Numerical Simulation bridging method” J. Appl. Mech., 73: 413–421). The turbulent length scale in the PANS model is constructed from the total turbulent kinetic energy and the sub-filter dissipation rate. Both hybrid models change between RANS (Reynolds-Averaged Navier–Stokes) and LES based on the cube root of the cell volume. The hybrid techniques, in contrast to RANS, are able to reproduce the turbulent flow dynamics in the shear layers of the impacting jet. The change from RANS to LES is much slower however for the PANS model than for the DES model on fine enough grids. This delays the break-up process of the vortices generated in the shear layers with as a consequence that the DES model produces better results than the PANS model.     

Superiority of PANS compared to LES in predicting a rudimentary landing gear flow with affordable meshes

Results of simulations of the flow around a rudimentary landing gear are presented in the paper. A newly proposed improved Partially-Averaged Navier–Stokes (PANS) method using k − ε − ζ − f turbulence model is used for prediction of the flow. The results are compared with the experimental data but also with the results of two LES simulations performed using the PANS computational grids. PANS simulations predicted the flow in good agreement with the experimental data. LES predicted a non-physical creation of separation over the front wheels that does not exist in the PANS prediction and was not observed in the experimental oil film. PANS simulations showed low sensitivity to the grid refinement. They show clear advantage compared with the LES simulations when the computational grid is inadequate for resolution of the near-wall flow structures.     

A low Reynolds number variant of partially-averaged Navier-Stokes model for turbulence

A low Reynolds number (LRN) formulation based on the Partially Averaged Navier-Stokes (PANS) modelling method is presented, which incorporates improved asymptotic representation in near-wall turbulence modelling. The effect of near-wall viscous damping can thus be better accounted for in simulations of wall-bounded turbulent flows. The proposed LRN PANS model uses an LRN k-ε model as the base model and introduces directly its model functions into the PANS formulation. As a result, the inappropriate wall-limiting behavior inherent in the original PANS model is corrected. An interesting feature of the PANS model is that the turbulent Prandtl numbers in the k and ε equations are modified compared to the base model. It is found that this modification has a significant effect on the modelled turbulence. The proposed LRN PANS model is scrutinized in computations of decaying grid turbulence, turbulent channel flow and periodic hill flow, of which the latter has been computed at two different Reynolds numbers of Re = 10,600 and 37,000. In comparison with available DNS, LES or experimental data, the LRN PANS model produces improved predictions over the standard PANS model, particularly in the near-wall region and for resolved turbulence statistics. Furthermore, the LRN PANS model gives similar or better results - at a reduced CPU time - as compared to the Dynamic Smagorinsky model.     

LES of turbulent jet in cross-flow: Part 1 – A numerical validation study

The paper presents results of a LES based numerical simulation of the turbulent jet-in-cross-flow (JICF) flowfield, with Reynolds number based on cross-flow velocity and jet diameter Re = 2400 and jet-to-cross-flow velocity ratio of R = 3.3. The JICF flow case has been investigated in great detail, involving conduction of two independent precursor simulations, prior to the main JICF simulation, as the considered case has turbulent inflow conditions on both jet and cross-stream side. The LES results are directly compared to pointwise Laser Doppler Anemometry (LDA) measurements, showing a very good agreement on the level of various statistical quantities in all flow regions but the immediate jet-to-cross-flow exhaustion zone. Several LES computations involving grids of up to 15 million grid points have been conducted, showing no improvement in the agreement between numerical results and measurements, possibly indicating a LDA measurement problem in this particular region.     

Hybrid RANS/LES of flow and heat transfer in round impinging jets

Fluid flow and convective heat transfer predictions are presented of round impinging jets for several combinations of nozzle-plate distances H/D = 2, 6 and 13.5 (where D is the nozzle diameter) and Reynolds numbers Re = 5000, 23,000 and 70,000 with the newest version of the k-ω model of Wilcox (2008) and three hybrid RANS/LES models. In the RANS mode of the hybrid RANS/LES models, the k-ω model is recovered. Three formulations are considered to activate the LES mode. The first model is similar to the hybrid models of Davidson and Peng (2003) and Kok et al. (2004). The turbulent length scale is replaced by the grid size in the destruction term of the k-equation and in the definition of the RANS eddy viscosity. As grid size, a maximum measure of the hexahedral grid cell is used. The second model has the same k-equation, but the eddy viscosity is the minimum of the k-ω eddy viscosity and the Smagorinsky eddy viscosity, following a proposal by Batten et al. (2004). The Smagorinsky eddy viscosity is formed with the cube root of the cell volume. The third model has, again, the same k-equation, but has an eddy viscosity which is an intermediate between the eddy viscosities of the first and second models. This is reached by using the cube root of the cell volume in the eddy viscosity formula of the first model.The simulation results are compared with experimental data for the high Reynolds number cases Re = 23,000 and Re = 70,000 and LES data for the low-Reynolds number case Re = 5000. The Reynolds numbers are defined with the nozzle diameter and the bulk velocity at nozzle outlet. At low nozzle-plate distance (the impingement plate is in the core of the jet), turbulent kinetic energy is overpredicted by RANS in the stagnation flow region. This leads to overprediction of the heat transfer rate along the impingement plate in the impact zone. At high nozzle-plate distance (the impingement plate is in the mixed-out region of the jet), the turbulence mixing is underpredicted by RANS in the shear layer of the jet which gives a too high length of the jet core. This also results in overprediction of the heat transfer rate in the impingement zone caused by too big temperature gradients at impingement.All hybrid RANS/LES models are able to correct the heat transfer overprediction of the RANS model. For good predictions at low nozzle-plate distance, it is necessary to sufficiently resolve the formation and development of the near-wall vortices in the jet impingement region. At high nozzle-plate distance, the essence is to capture the evolution and breakup of the flow unsteadiness in the shear layer of the jet, so that accurate mean and fluctuating velocity profiles are obtained in the impingement region. Although the models have a quite different theoretical justification and generate a quite different eddy viscosity in some flow regions, their overall results are very comparable. The reason is that in zones that are crucial for the results, the models behave similarly.     

Euler–Euler large eddy simulations for dispersed turbulent bubbly flows

In this paper we present detailed Euler–Euler Large Eddy Simulations (LES) of dispersed bubbly flow in a rectangular bubble column. The motivation of this study is to investigate the potential of this approach for the prediction of bubbly flows, in terms of mean quantities. The physical models describing the momentum exchange between the phases including drag, lift and wall force were chosen according to previous experiences of the authors. Experimental data, Euler–Lagrange LES and unsteady Euler–Euler Reynolds-Averaged Navier–Stokes results are used for comparison. It is found that the present model combination provides good agreement with experimental data for the mean flow and liquid velocity fluctuations. The energy spectrum obtained from the resolved velocity of the Euler–Euler LES is presented as well.     

A reformulated synthetic turbulence generation method for a zonal RANS–LES method and its application to zero-pressure gradient boundary layers

A synthetic turbulence generation (STG) method for subsonic and supersonic flows at low and moderate Reynolds numbers to provide inflow distributions of zonal Reynolds-averaged Navier–Stokes (RANS) – large-eddy simulation (LES) methods is presented. The STG method splits the LES inflow region into three planes where a local velocity signal is decomposed from the turbulent flow properties of the upstream RANS solution. Based on the wall-normal position and the local flow Reynolds number, specific length and velocity scales with different vorticity content are imposed at the inlet plane of the boundary layer. The quality of the STG method for incompressible and compressible zero-pressure gradient boundary layers is shown by comparing the zonal RANS–LES data with pure LES, pure RANS, and direct numerical simulation (DNS) solutions. The distributions of the time and spanwise wall-shear stress, Reynolds stress distributions, and two point correlations of the zonal RANS–LES simulations are smooth in the transition region and in good agreement with the pure LES and reference DNS findings. The STG approach reduces the RANS-to-LES transition length to less than four boundary-layer thicknesses.     

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.     

Reducing the pressure drag of a D-shaped bluff body using linear feedback control

The pressure drag of blunt bluff bodies is highly relevant in many practical applications, including to the aerodynamic drag of road vehicles. This paper presents theory revealing that a mean drag reduction can be achieved by manipulating wake flow fluctuations. A linear feedback control strategy then exploits this idea, targeting attenuation of the spatially integrated base (back face) pressure fluctuations. Large-eddy simulations of the flow over a D-shaped blunt bluff body are used as a test-bed for this control strategy. The flow response to synthetic jet actuation is characterised using system identification, and controller design is via shaping of the frequency response to achieve fluctuation attenuation. The designed controller successfully attenuates integrated base pressure fluctuations, increasing the time-averaged pressure on the body base by 38%. The effect on the flow field is to push the roll-up of vortices further downstream and increase the extent of the recirculation bubble. This control approach uses only body-mounted sensing/actuation and input–output model identification, meaning that it could be applied experimentally.     

Wall jet similarity of impinging planar underexpanded jets

Velocity profiles and wall shear stress values in the wall jet region of planar underexpanded impinging jets are parameterized based on nozzle parameters (stand-off height, jet hydraulic diameter, and nozzle pressure ratio). Computational fluid dynamics is used to calculate the velocity fields of impinging jets with height-to-diameter ratios in the range of 15–30 and nozzle pressure ratio in the range of 1.2–3.0. The wall jet has an incomplete self-similar profile with a typical triple-layer structure as in traditional wall jets. The effects of compressibility are found to be insignificant for wall jets with Ma < 0.8. Wall jet analysis yielded power-law relationships with source dependent coefficients describing maximum velocity, friction velocity, and wall distances for maximum and half-maximum velocities. Source dependency is determined using the conjugate gradient method. These power-law relationships can be used for mapping wall shear stress as a function of nozzle parameters.     

The effect of a systematic change in surface roughness skewness on turbulence and drag

Previous work by the authors (Flack and Schultz, 2010) has identified the root-mean-square roughness height, k_rms, and the skewness, Sk, of the surface elevation distribution as important parameters in scaling the skin-friction drag on rough surfaces. In this study, three surfaces are tested in turbulent boundary layer flow at a friction Reynolds number, Re_τ = 1600–2200. All the surfaces have similar root-mean-square roughness height, while the skewness is varied. Measurements are presented using both two-component LDV and PIV. The results show the anticipated trend of increasing skin-friction drag with increasing skewness. The largest increase in drag occurs going from negative skewness to zero skewness with a more modest increase going from zero to positive skewness. Some differences in the mean velocity and Reynolds stress profiles are observed for the three surfaces. However, these differences are confined to a region close to the rough surface, and the mean velocity and Reynolds stress profiles collapse away from the wall when scaled in outer variables. The turbulence structure as documented through two-point spatial correlations of velocity is also observed to be very similar over the three surfaces. These results support Townsend’s (1976) concept of outer-layer similarity that the wall boundary condition exerts no direct influence on the turbulence structure away from the wall except in setting the velocity and length scales for the outer layer.     

A hybrid RANS/LES simulation of high-Reynolds-number channel flow using additional filtering at the interface

A hybrid method combining large eddy simulation (LES) with the Reynolds-averaged Navier-Stokes (RANS) equation is used to simulate a turbulent channel flow at high Reynolds number. It is known that the mean velocity profile has a mismatch between the RANS and LES regions in hybrid simulations of a channel flow. The velocity mismatch is reproduced and its dependence on the location of the RANS/LES interface and on the type of RANS model is examined in order to better understand its properties. To remove the mismatch and to obtain better velocity profiles, additional filtering is applied to the velocity components in the wall-parallel planes near the interface. The additional filtering was previously introduced to simulate a channel flow at low Reynolds number. It is shown that the filtering is effective in reducing the mismatch even at high Reynolds number. Profiles of the velocity fluctuations of runs with and without the additional filtering are examined to help understand the reason for the mismatch. Due to the additional filtering, the wall-normal velocity fluctuation increases at the bottom of the LES region. The resulting velocity field creates the grid-scale shear stress more efficiently, and an overestimate of the velocity gradient is removed. The dependence of the velocity profile on the grid point number is also investigated. It is found that the velocity gradient in the core region is underestimated in the case of a coarse grid. Attention should be paid not only to the velocity mismatch near the interface but also to the velocity profile in the core region in hybrid simulations of a channel flow at high Reynolds number. PACS47.27.Eq; 47.27.Nz; 47.60.+i     

Direct noise computation of adaptive control applied to a cavity flowSimulation directe de l'atténuation par contrôle adaptatif du rayonnement d'une cavité affleurée par un écoulement

The Large Eddy Simulation of closed-loop active flow control applied to a 3D cavity excited by a compressible airflow with a Mach number of 0.6 is presented. The control actuator is an idealized synthetic jet located at the upstream cavity edge, and the control function is supplied by a feedback LMS-type algorithm whose input is a pressure signal measured inside the cavity. The radiated sound, provided directly by the LES simulation, was shown to decrease substantially when active control was applied. A simultaneous reduction of the vertical velocity fluctuations in the shear layer was observed. The intensity of vortical structures inside the cavity was also reduced, although the general aspect of the recirculation zone was not modified. The direct noise computation technique, which supplies the pressure field by solving the fluid mechanics equations, is shown to constitute a powerful tool for studying active aeroacoustic noise control. To cite this article: O. Marsden et al., C. R. Mecanique 331 (2003).     

A Hybrid-Filter Approach to Turbulence Simulation

Germano (Theor Comput Fluid Dyn 17:225–331, 2004) proposed a hybrid-filter approach, which additively combines an LES-like filter operator (F) and a RANS-like statistical operator (E) using a blending function k: H?=?kF?+?(1???k)E. Using turbulent channel flow as an example, we first conducted a priori tests in order to gain some insights into this hybrid-filter approach, and then performed full simulations to further assess the approach in actual simulations. For a priori tests, two separate simulations, RANS (E) and LES (F), were performed using the same grid in order to construct a hybrid-filtered field (H). It was shown that the extra terms arising out of the hybrid-filtered Navier–Stokes (HFNS) equations provided additional energy transfer from the RANS region to the LES region, thus alleviating the need for the ad hoc forcing term that has been used by some investigators. The complexity of the governing equations necessitated several modifications in order to render it suitable for a full numerical simulation. Despite some issues associated with the numerical implementation, good results were obtained for the mean velocity and skin friction coefficient. The mean velocity profile did not have an overshoot in the logarithmic region for most blending functions, confirming that proper energy transfer from the RANS to the LES region was a key to successful hybrid models. It is shown that Germano’s hybrid-filter approach is a viable and mathematically more appealing approach to simulate high Reynolds number turbulent flows.     

The interaction of an asymmetrical localised synthetic jet on a side-supported sphere

A localised synthetic jet offers promise of an optimum and cost-effective practical method of delaying separation and promoting reattachment in fluids with solid body interactions. The asymmetric flow that may result from its use may also be beneficial in improving the aerodynamic performance of a lifting body. There are insufficient studies of synthetic jets, particularly on three-dimensional bluff bodies that are more representative of complex flows in real situations. A comprehensive study on an 80 mm diameter sphere designed with localised synthetic jet orifices was, therefore, conducted in an 18 in×18 in open circuit closed test-section wind tunnel at a Reynolds number of 5×10⁴. The coefficient of pressure distribution was measured by continuously varying the location of the synthetic jet and compared with the no synthetic jet condition. The three-dimensional effects on the flow over the sphere body are particularly made apparent through the growth and the effects of the boundary layer and the deviation from potential flow. Overall, the synthetic jet had the effect of delaying the separation point and extending it further downstream on the sphere surface concomitantly producing a significant reduction in drag, providing solid support to the viability of strategically located synthetic jet when higher lift or lower drag is desired. A surprising discovery was the ability of the synthetic jet to improve the flow at the junction of the sting support and sphere. This has promising implications in devising methods to reduce interference drag that are common in many practical applications such as near junctions between wing and the fuselage.     

合成射流对失速状态下翼型大分离流动控制的试验研究

为研究低速状态合成射流在抑制翼型气流分离和推迟失速方面的控制机理, 开展了NACA0021 翼型失速特性射流控制的风洞试验研究. 通过系统性的模型测力、翼型瞬态流场粒子图像测速和边界层速度测定的对比试验, 深入探索了合成射流各参数对翼型失速特性控制效果的影响规律. 试验结果表明射流偏角在翼型升力和失速迎角控制方面的效果对射流动量系数较为敏感: 当动量系数较大时, 近切向射流的控制效果更好. 射流动量系数为0.033 时, 偏角30°的射流使得翼型升力系数峰值提高23.56%, 失速迎角增大5°; 而动量系数较小时, 偏角较大的射流能够获得最佳控制效果. 射流动量系数为0.0026 时, 法向射流对翼型最大升力系数控制效果最好(提高9.2%).      

LES of the flow around several cuboids in a row

This paper presents Large Eddy Simulations (LES) of flow around a four-vehicle platoon when one of the platoon members was forced to undergo in-line oscillations. The LES were made at the Reynolds number of 10⁵ based on the height of the vehicles. Combinations of two different frequencies corresponding to non-dimensional frequencies at the Strouhal numbers St₁ = 0.025 and St₂ = 0.013 and two oscillation amplitudes were used in this study. The methodology was validated by comparisons with data from previous experimental investigations. In order to highlight the dynamic effects, comparisons were made with steady results on a single vehicle and on a four-vehicle platoon. Large differences were found in the flow structures between quasi-steady and dynamic results. Furthermore, the behavior of the drag coefficient of the upstream neighbor of the oscillating model was investigated.     

Active flow control on the 25°Ahmed body using a new unsteady jet

A numerical simulation based on the Large eddy simulation method is carried out on the near wake flow behind a 25° slant angle Ahmed body to analyze and establish a new method to control the near wake flow. An active flow control using a new unsteady jet derived from the traditional synthetic jet is applied to reduce the aerodynamic drag. The control devices are distributed along the separation edges on the rear part of the body. Their effects on the near wake and the rear body by influencing the flow topology and the static pressure distribution are examined respectively. The control frequency of the jet as the key forcing parameter is taken into consideration as well. The different actuation set-ups lead to a max drag reduction of up to 13.6%, which demonstrates a good correlation with the static pressure distribution at the rear end of the body.     

Numerical simulation of the flow around a simplified vehicle model with active flow control

Large-eddy simulation (LES) was used to study the influence and the resulting flow mechanisms of active flow control applied to a two-dimensional vehicle geometry. The LES results were validated against existing Particle Image Velocimetry (PIV) and force measurement data. This was followed by an exploration of the influence of flow actuation on the near-wake flow and resulting aerodynamic forces. Not only was good agreement found with the previous experimental study, but new knowledge was gained in the form of a complex interaction of the actuation with the coherent flow structures. The resulting time-averaged flow shows a strong influence of the extension of the actuation slots and the lateral solid walls on the near-wake flow structures and thereby on the resulting drag.     

Turbulent Scalar Mixing in a Skewed Jet in Crossflow: Experiments and Modeling

Turbulent mixing of an inclined, skewed jet injected into a crossflow is investigated using MRI-based experiments and a high-fidelity LES of the same configuration. The MRI technique provides three-dimensional fields of mean velocity and mean jet concentration. The 30° skew of the jet relative to the crossflow produces a single dominant vortex which introduces spanwise asymmetries to the velocity and concentration fields. The turbulent scalar transport of the skewed jet is investigated in further detail using the LES, which is validated against the experimental measurements. Mixing is found to be highly anisotropic throughout the jet region. Isotropic turbulent diffusivity and viscosity are used to calculate an optimal value of the turbulent Schmidt number, which varies widely over the jet region and lies mostly outside of the typically accepted range 0.7 ≤ S c _t ≤ 0.9. Finally, three common scalar flux models of increasing complexity are evaluated based on their ability to capture the anisotropy and predict the scalar concentration field of the present configuration. The higher order models are shown to better represent the turbulent scalar flux vector, leading to more accurate calculations of the concentration field. While more complex models are better able to capture the turbulent mixing, optimization of model constants is shown to significantly affect the results.