Large eddy simulations are performed to investigate the impact of a solid obstacle on the flow around a multiperforated plate typical of aeronautical combustion chambers. The reference configuration is a perforated plate with approximately 200 holes immersed between a cold vein and hot vein at a typical operating point of helicopter combustors. The micro-jet Reynolds number is of the order of 4000, while the blowing and momentum ratios are close to 4 and 8, respectively. A variant configuration is considered that features an additional cylindrical obstacle located in the cold vein and mimicking a spark plug. The study reveals that, downstream of the obstacle, the cooling effectiveness of the plate is reduced by approximately 40 % compared with the reference case, mainly due to the absence of perforation at the obstacle location. The mass flow rate within the holes in the wake produced by the obstacle is reduced by 7 %, which is likely to locally influence the plate cooling. The reduction is attributed to the wake's pressure loss and its impact on the discharge coefficient. Additionally, the cooling effectiveness outside the wake shows a 5 % increase that can be linked to the mass flow rate increase within corresponding holes.

The performance of horizontal dilution tubes is investigated by Reynolds-averaged Navier–Stokes and large-eddy simulations. The flame gas enters the dilution tube through a pinhole. The orifice flow and the dilution process inside the tube are studied. The volume flow through the orifice is shown to be proportional to the square root of the pressure drop. The discharge coefficient is 0.9 ± 0.3 in the cold air (calibration) case and drops to 0.35 under hot (flame) conditions. The resulting dilution ratio is roughly a factor of five below typical literature data. The gas sample remains in the wall boundary layer and the mixing process is not complete at the end of the dilution tube. Turbulence decays rapidly behind the tube inlet, which shifts the flow into the laminar to turbulent transition regime. Turbulence increases significantly in the outlet section which has much smaller pipe cross-sections. Despite its relatively low Reynolds number, the outlet flow to the particle sizer (or to the gas analyzer) is clearly turbulent, and interactions with the wall are probable. The results are in agreement with previous findings from laminar jets in cross-flow. Guidelines for optimization of the sampling conditions are suggested.

While there have been numerous applications of large eddy simulations (LES) to complex flows, their application to practical engineering configurations, such as full aircraft models, have been limited to date. Recently, however, advances in rapid, high quality mesh generation, low-dissipation numerical schemes and physics-based subgrid-scale and wall models have led to, for the first time, accurate simulations of a realistic aircraft in landing configuration (the Japanese Aerospace Exploration Agency Standard Model) in less than a day of turnaround time with modest resource requirements. In this paper, a systematic study of the predictive capability of LES across a range of angles of attack (including maximum lift and post-stall regimes), the robustness of the predictions to grid resolution and the incorporation of wind tunnel effects is carried out. Integrated engineering quantities of interest, such as lift, drag and pitching moment will be compared with experimental data, while sectional pressure forces will be used to corroborate the accuracy of the integrated quantities. Good agreement with experimental $C_L$ data is obtained across the lift curve with the coefficient of lift at maximum lift, $C_{L,max}$, consistently being predicted to within five lift counts of the experimental value. The grid point requirements to achieve this level of accuracy are reduced compared with recent estimates (even for wall modelled LES), with the solutions showing systematic improvement upon grid refinement, with the exception of the solution at the lowest angles of attack, which will be discussed later in the text. Simulations that include the wind tunnel walls and aircraft body mounting system are able to replicate important features of the flow field noted in the experiment that are absent from free air calculations of the same geometry, namely, the onset of inboard flow separation in the post-stall regime. Turnaround times of the order of a day are made possible in part by algorithmic advances made to leverage graphical processing units. The results presented herein suggest that this combined approach (meshing, numerical algorithms, modelling, efficient computer implementation) is on the threshold of readiness for industrial use in aeronautical design.

Prediction of aeroengine exhaust plume near-field development requires knowledge of velocity and turbulence distributions at nozzle exit. The high Reynolds number nozzle inlet boundary layers of engineering practice are fully turbulent, but acceleration can induce re-laminarisation. Thus, to reproduce nozzle exit conditions accurately, large eddy simulation (LES) for plume prediction must be capable of capturing re-laminarisation and any subsequent boundary layer recovery. Validation is essential to establish a credible LES methodology, but previous studies have suffered from lack of nozzle inlet/exit measurements in the test cases selected. Validation data were here taken from an experiment on a convergent round nozzle with a parallel exit extension to allow boundary layer recovery. LES inlet condition generation applied a rescaling/recycling method (R2M), whose performance was validated against measurements of first and second moment statistics as well as the turbulence integral length scale. Simulations employed two sub-grid-scale (SGS) models; these produced similar predictions up to the end of the nozzle convergent section, but marked differences appeared for the nozzle exit turbulence field. The Smagorinsky model predicted much lower turbulence levels than measured, whereas the Piomelli and Geurts model revealed the presence of a small separation region at the convergence/parallel section corner, which led to higher exit turbulence and much better agreement with measured data.

Numerical investigation of the underexpanded sonic coaxial jets is carried out using large eddy simulation for three typical inner nozzle lip-thicknesses. Various fundamental mechanisms dictating the flow phenomena including shock structure, shear layer evolution and sound production are investigated. It is found that the inner nozzle lip induces a recirculation zone between inner and outer jets, which significantly influences the behaviors of shock structures and shear layers. The sound properties of the coaxial jets are further analyzed in detail. As the inner lip-thickness increases, the helical screech mode switches to an axisymmetric one and high-frequency screech also occurs with an oscillation frequency of recirculation zone. Based on the temporal Fourier transform and correlation analysis, the primary sources of low- and high-frequency screeches are associated with the downstream shock cells in the jet column and the secondary shock structures in the outer annular jet, respectively. The proper orthogonal decomposition analysis reveals that the dominant structures constructed by the most energetic modes shift from the downstream shock cells region to the upstream secondary shock region as the lip-thickness increases. The results obtained in this study provide physical insight into the understanding of the mechanisms relevant to the coherent structures and sound properties in sonic coaxial jets.

Jet impingement induced heat transfer is an important issue in engineering science. This paper presents results of large eddy simulation (LES) of normally impinging elliptic air-jet heat transfer at a Reynolds number of 4400, with orifice-to-plate distance fixed to be 5 in the unit of jet nozzle effective diameter . The elliptic aspect ratio (a/b) is 3/2. While the target wall is heated under some condition of constant heat flux. The LES are carried out using dynamic subgrid model and Open-FOAM. The distributions ofmean velocity components, velocity fluctuations, and subgrid stresses in vertical and radial directions, and the Nusselt numbers involving heat transfer through the target wall are discussed. The comparison with existing experimental and numerical results shows good agreement.

Scale similarity is found in many natural phenomena in the universe, from fluid dynamics to astrophysics. In large eddy simulations of turbulent flows, some sub-grid scale (SGS) models are based on scale similarity. The earliest scale similarity SGS model was developed by Bardina et al., which produced SGS stresses with good correlation to the true stresses. In the present study, we perform a mathematical analysis of scale similarity. The analysis has revealed that the ratio of the resolved stress to the SGS stress is γ2, where γ is the ratio of the second filter width to the first filter width, under the assumption of small filter width. The implications of this analysis are discussed in the context of large eddy simulation.

Dynamics and instability of a vortex ring impinging on a wall were investigated by means of large eddy simulation for two vortex core thicknesses corresponding to thin and thick vortex rings. Various fundamental mechanisms dictating the flow behaviors, such as evolution of vortical structures, formation of vortices wrapping around vortex rings, instability and breakdown of vortex rings, and transition from laminar to turbulent state, have been studied systematically. The evolution of vortical structures is elucidated and the formation of the loop-like and hair-pin vortices wrapping around the vortex rings (called briefly wrapping vortices) is clarified. Analysis of the enstrophy of wrapping vortices and turbulent kinetic energy (TKE) in flow field indicates that the formation and evolution of wrapping vortices are closely associated with the flow transition to turbulent state. It is found that the temporal development of wrapping vortices and the growth rate of axial flow generated around the circumference of the core region for the thin ring are faster than those for the thick ring. The azimuthal instabilities of primary and secondary vortex rings are analyzed and the development of modal energies is investigated to reveal the flow transition to turbulent state. The modal energy decay follows a characteristic –5/3 power law, indicating that the vortical flow has become turbulence. Moreover, it is identified that the TKE with a major contribution of the azimuthal component is mainly distributed in the core region of vortex rings. The results obtained in this study provide physical insight of the mechanisms relevant to the vortical flow evolution from laminar to turbulent state.

A vortex ring impacting a three-dimensional bump is studied using large eddy simulation for a Reynolds number Re = 4 × 104 based on the initial diameter and translational speed of the vortex ring. The effects of bump height and vortex core thickness for thin and thick vortex rings on the vortical flow phenomena and the underlying physical mechanisms are investigated. Based on the analysis of the evolution of vortical structures, two typical kinds of vortical structures, i.e., the wrapping vortices and the hair-pin vortices, are identified and play an important role in the flow state evolution. The boundary vorticity flux is analyzed to reveal the mechanism of the vorticity generation on the bump surface. The circulation of the primary vortex ring reasonably elucidates some typical phases of flow evolution. Further, the analysis of turbulent kinetic energy reveals the transition from laminar to turbulent state. The results obtained in this study provide physical insight into the understanding of the mechanisms relevant to the flow evolution and the flow transition to turbulent state.

A submerged turbulent plane jet in shallow water impinging vertically onto the free surface will produce a large-scale flapping motion when the jet exit velocity is larger than a critical one. The flapping phenomenon is verified in this paper through a large eddy simulation where the free surface is modeled by volume of fluid approach. The quantitative results for flapping jet are found to be in good agreement with available experimental data in terms of mean velocity, flapping-induced velocity and turbulence intensity. Results show that the flapping motion is a new flow pattern with characteristic flapping frequency for submerged turbulent plane jets, the mean centerline velocity decay is considerably faster than that of the stable impinging jet without flapping motion, and the flapping-induced velocities are as important as the turbulent fluctuations.

We consider in this paper the problem of finding appropriate models for Large Eddy Simulations of turbulent incompressible flows from a mathematical point of view. The Smagorinsky model is analyzed and the vorticity formulation of the Navier–Stokes equations is used to explore more efficient subgrid-scale models as minimal regularizations of these equations. Two classes of variants of the Smagorinsky model emerge from this approach: a model based on anisotropic turbulent viscosity and a selective model based on vorticity angles. The efficiency of these models is demonstrated by comparisons with reference results on decaying turbulence experiments.