Search results

Aircraft noise is dominant for residents near airports when planes fly at low altitudes such as during departure and landing. Flaps, wings, landing gear contribute significantly to the total sound emission. This paper aims to present a passive flow control (in the sense that there is no power input) to reduce the noise radiation induced by the flow over the cavity of the landing gear during take-off and landing.

The understanding of the noise source mechanism is normally caused by the unsteady interactions between the cavity surface and the turbulent flows as well as some studies that have shown tonal noise because of cavity resonances; this tonal noise is dependent on cavity geometry and incoming flow that lead us to use of a sinusoidal surface modification application upstream of a cavity as a passive acoustics control device in approach conditions.

It is demonstrated that the proposed surface waviness showed a potential reduction in cavity resonance and in the overall sound pressure level at the majority of the points investigated in the low Mach number. Furthermore, optimum sinusoidal amplitude and frequency were determined by the means of a two-dimensional computational fluid dynamics analysis for a cavity with a length to depth ratio of four.

The noise control by surface waviness has not implemented in real flight test yet, as all the tests are conducted in the credible numerical simulation.

The application of passive control method on the cavity requires a global aerodynamic study of the air frame is a matter of ongoing debate between aerodynamicists and acousticians. The latter is aimed at the reduction of the noise, whereas the former fears a corruption of flow conditions. To balance aerodynamic performance and acoustics, the use of the surface waviness in cavity leading edge is the most optimal solution.

The proposed leading-edge modification it has important theoretical basis and reference value for engineering application it can meet the demands of engineering practice. Particularly, to contribute to the reduce the aircraft noise adopted by the “European Visions 2020”.

The investigate cavity noise with and without surface waviness generation and propagation by using a hybrid approach, the computation of flow based on the large-eddy simulation method, is decoupled from the computation of sound, which can be performed during a post-processing based on Curle’s acoustic analogy as implemented in OpenFOAM.

The simulations of grid-resolved rod vortex generators (RVGs) require high computational cost and time. Additionally, the computational mesh topology must be adjusted to rods geometries. The purpose of this study is to propose the new source term model for RVG.

The model was proposed by modification of Bender, Anderson, Yagle (BAY) model used to predict flows around different type of vortex generators (VGs) – vanes. Original BAY model was built on lifting line theory. The proposed model was implemented in ANSYS Fluent by means of the user-defined function technique. Additional momentum and energy sources are imposed to transport equations.

The computational results of source term model were validated against experimental data and numerical simulation results for grid-resolved rod. It was shown that modified BAY model can be successfully used for RVG in complex cases. An example of BAY model application for RVG on transonic V2C airfoil with strongly oscillating shock waves is presented. Aerodynamic performance predicted numerically by means of both approaches (grid resolved RVG and modeled) is in good agreement, what indicates application opportunity of the proposed model to complex cases.

Modified BAY model can be used to simulate the influence of RVGs in complex real cases. It allows for time/cost reduction if the location or distribution of RVG has to be optimized on a profile, wing or in the channel.

In the paper, the new modification of BAY model was proposed to simulate RVGs. The presented results are innovative because of original approach to model RVGs.

The purpose of this study is to investigate the accuracy and applicability of the Flowfield Dependent Variation (FDV) method for large-eddy simulations (LES) of decaying isotropic turbulence.

In an earlier paper, the FDV method was successfully demonstrated for simulations of laminar flows with speeds varying from low subsonic to high supersonic Mach numbers. In the current study, the FDV method, implemented in a finite element framework, is used to perform LESs of decaying isotropic turbulence. The FDV method is fundamentally derived from the Lax–Wendroff Scheme (LWS) by replacing the explicit time derivatives in LWS with a weighted combination of explicit and implicit time derivatives. The increased implicitness and the inherent numerical dissipation of FDV contribute to the scheme’s numerical stability and monotonicity. Understanding the role of numerical dissipation that is inherent to the FDV method is essential for the maturation of FDV into a robust scheme for LES of turbulent flows. Accordingly, three types of LES of decaying isotropic turbulence were performed. The first two types of LES utilized explicit subgrid scale (SGS) models, namely, the constant-coefficient Smagorinsky and dynamic Smagorinsky models. In the third, no explicit SGS model was employed; instead, the numerical dissipation inherent to FDV was used to emulate the role played by explicit SGS models. Such an approach is commonly known as Implicit LES (ILES). A new formulation was also developed for quantifying the FDV numerical viscosity that principally arises from the convective terms of the filtered Navier–Stokes equations.

The temporal variation of the turbulent kinetic energy and enstrophy and the energy spectra are presented and analyzed. At all grid resolutions, the temporal profiles of kinetic energy showed good agreement with t^(−1.43) theoretical scaling in the fully developed turbulent flow regime, where t represents time. The energy spectra also showed reasonable agreement with the Kolmogorov’s k^(−5/3) power law in the inertial subrange, with the spectra moving closer to the Kolmogorov scaling at higher-grid resolutions. The intrinsic numerical viscosity and the dissipation rate of the FDV scheme are quantified, both in physical and spectral spaces, and compared with those of the two SGS LES runs. Furthermore, at a finite number of flow realizations, the numerical viscosities of FDV and of the Streamline Upwind/Petrov–Galerkin (SUPG) finite element method are compared. In the initial stages of turbulence development, all three LES cases have similar viscosities. But, once the turbulence is fully developed, implicit LES is less dissipative compared to the two SGS LES runs. It was also observed that the SUPG method is significantly more dissipative than the three LES approaches.

Just as any computational method, the limitations are based on the available computational resources.

Solving problems involving turbulent flows is by far the biggest challenge facing engineers and scientists in the twenty-first century, this is the road that the authors have embarked upon in this paper and the road ahead of is very long.

Understanding turbulence is a very lofty goal and a challenging one as well; however, if the authors succeed, the rewards are limitless.

The derivation of an explicit expression for the numerical viscosity tensor of FDV is an important contribution of this study, and is a crucial step forward in elucidating the fundamental properties of the FDV method. The comparison of viscosities for the three LES cases and the SUPG method has important implications for the application of ILES approach for turbulent flow simulations.

The purpose of this paper is to develop a subgrid-scale (SGS) model for large eddy simulation (LES) of buoyancy- and thermally driven transitional and turbulent flows and further examine its performance.

Favre-filtered, non-dimensional LES equations are solved using non-dissipative, fully implicit, kinetic energy conserving, finite-volume algorithm which uses an iterative predictor-corrector approach based on pressure correction. Also, to develop a new SGS model which accounts for buoyancy, turbulent generation term in SGS viscosity is properly modified and enhanced by buoyancy production.

The proposed model has been successfully applied to turbulent Rayleigh–Bénard convection. The results show that the model is able to reproduce the complex physics of turbulent thermal convection. In comparison with the original wall-adapting local eddy-viscosity (WALE) and buoyancy-modified (BM) Smagorinsky models, turbulent diagnostics predicted by the new model are in better agreement with direct numerical simulation.

A BM variant of the WALE SGS model is newly developed and analyzed.

The purpose of this paper is to numerically investigate the mildly separated flow phenomena on a near-stall NACA0015 airfoil, by using Detached-Eddy Simulation (DES) type methods. It includes a comparison of different choices of underlying Reynolds-averaged Navier–Stokes model as well as subgrid-scale stress model in Large-Eddy simulation mode.

The unsteady flow phenomena are simulated by using delayed DES (DDES) and improved DDES (IDDES) methods, with an in-house computational fluid dynamics solver. Characteristic frequencies in different flow regions are extracted using fast Fourier transform. Dynamic mode decomposition (DMD) method is applied to uncover the critical dynamic modes.

Among all the DES type methods investigated in this paper, only the Spalart–Allmaras-based IDDES captures the separation point as measured in the experiments. The classical vortex-shedding and the shear-layer flapping modes for airfoil flows with shallow separation are also found from the IDDES results by using DMD.

The value of this paper lies in the assessment of five different DES-type models through the detailed investigation of the Reynolds stresses as well as the separation and reattachment.

This study aims to investigate effects of sonic jet injection into supersonic cross-flow (JISC) numerically in different dynamic pressure ratio values and free stream Mach numbers.

Large Eddy simulation (LES) with dynamic Smagorinsky model is used as the turbulence approach. The numerical results are compared with the experimental data, and the comparison shows acceptable validation.

According to the results, the dynamic pressure ratio has critical effects on the zone related to barrel shock. Despite free stream Mach number, increasing dynamic pressure ratio leads to expansion of barrel shock zone. Consequently, expanded barrel shock zone would bring about more obstruction effect. In addition, the height of counter-rotating vortex pair increases, and the high-pressure area before jet and low-pressure area after jet will rise. The results show that the position of barrel shock is deviated by increasing free stream Mach number, and the Bow shock zone becomes stronger and close to barrel shock. Moreover, high pressure zone, which is located before the jet, decreases by high free stream Mach number.

In this study, LES with a dynamic Smagorinsky model is used as the turbulence approach. Effects of sonic JISC are investigated numerically in different dynamic pressure ratio values and free stream Mach numbers.

As summary, the following are the contribution of this paper in the field of JISC subjects: several case studies of jet condition have been performed. In all the cases, the flow at the nozzle exit is sonic, and the free stream static pressure is constant. To generate proper grid, a cut cell method is used for domain modelling. Boundary condition effect on the wall pressure distribution around the jet and velocity profiles, especially S shape profiles, is investigated. The results show that the relation between representing the location of Mach disk centre and at transonic regime is a function of second-order polynomial, whereas at supersonic regime, the relationship is modelled as a first-order polynomial. In addition, the numerical results are compared with the experimental data demonstrating acceptable validation.

The purpose of the paper is to propose some modifications to the SIMPLE (semi-implicit method for pressure-linked equations) algorithm. These modifications can ensure the numerical robustness and optimize computational efficiency. They remarkably promote the ability of the SIMPLE algorithm for incompressible DNS (direct numerical simulation) of multiscale problems, such as transitional flows and turbulent flows, by improving the properties of dispersion and dissipation.

The MDCD (minimized dispersion and controllable dissipation) scheme and MMIM (modified momentum interpolation method) are introduced. Six typical test cases are used to validate the modified algorithm, including the linear convective flow, lid-driven cavity flow, laminar boundary layer, Taylor vortex and DHIT (decaying homogenous isotropic turbulence). Particularly, a highly unsteady DNS of separated-flow transition in turbomachinery is precisely predicted by the modified algorithm.

The numerical examples show the distinct superiority of the modified algorithm in both internal flows and external flows. The advantages of the MDCD scheme and MMIM make the SIMPLE algorithm a promising method for DNS.

Some effective modifications to the SIMPLE algorithm are addressed. It is the first attempt to introduce the MDCD approach into the SIMPLE-type algorithms. The new algorithm is especially suitable for the incompressible DNS of convection-dominated flows.

This paper aims to show results of numerical simulations of transonic flow around a supercritical airfoil at chord Reynolds number Re_c = 3 × 10⁶, with the aim of elucidating the mechanisms responsible for large-scale shock oscillations, namely, transonic buffet.

Unsteady Reynolds-averaged Navier–Stokes simulations and detached-eddy simulations provide a preliminary buffet map, while a high fidelity implicit large-eddy simulation with an upstream laminar boundary layer is used to ascertain the physical feasibility of the various buffet mechanisms. Numerical experiments with unsteady RANS highlight the role of waves travelling on pressure side in the buffet mechanism. Estimates of the propagation velocities of coherent disturbances and of acoustic waves are obtained, to check the validity of popular mechanisms based on acoustic feedback from the trailing edge.

Unsteady RANS numerical experiments demonstrate that the pressure side of the airfoil plays a marginal role in the buffet mechanism. Implicit LES data show that the only plausible self-sustaining mechanism involves waves scattered from the trailing edge and penetrating the sonic region from above the suction side shock. An interesting side result of this study is that buffet appears to be more intense in the case that the boundary layer state upstream of the shock is turbulent, rather than laminar.

The results of the study will be of interest to any researcher involved with transonic buffet.

The paper aims to expand the scope of application of the lattice Boltzmann method (LBM), especially in the field of aircraft engineering. The traditional LBM is usually applied to incompressible flows at a low Reynolds number, which is not sufficient to satisfy the needs of aircraft engineering. Devoted to tackling the defect, the paper proposes a developed LBM combining the subgrid model and the multiple relaxation time (MRT) approach. A multilayer adaptive Cartesian grid method to improve the computing efficiency of the traditional LBM is also employed.

The subgrid model and the multilayer adaptive Cartesian grid are introduced into MRT-LBM for simulations of incompressible flows at a high Reynolds number. Validated by several typical flow simulations, the numerical methods in this paper can efficiently study the flows under high Reynolds numbers.

Some numerical simulations for the lid-driven flow of cavity, flow around iced GLC305, LB606b and ONERA-M6 are completed. The paper presents the investigation results, indicating that the methods are accurate and effective for the separated flow after icing.

LBM is developed with the addition of the subgrid model and the MRT method. A numerical strategy is proposed using a multilayer adaptive Cartesian grid method and its treatment of boundary conditions. The paper refers to innovative algorithm developments and applications to the aircraft engineering, especially for iced wing simulations with flow separations.

The purpose of this paper devoted to the application of modal analysis to analyze the flow structure of trailing edge cutback film cooling and the effects of vortex structure on the film cooling effectiveness of the cutback surface.

Large eddy simulation (LES) is used to simulate the trailing edge cutback film cooling. The results of LES are analyzed by proper orthogonal decomposition (POD) method and dynamic mode decomposition (DMD) method. The POD method is used to determine the dominated vortex structure and the energy level of these structures. The DMD method is used to analyze the relationship between vortex structures and wall temperature.

The POD method shows that the flow field consists of three main vortices – streamwise vortex, lip vortex and coolant vortex. The DMD results show that the lip vortex mainly acts on the middle section of the cutback surface, while the streamwise vortex mainly acts on the back section of the cutback surface.

The modal analysis is only based on numerical simulation but the modal analysis of experimental results will be further studied in the future.

This paper presents the powerful ability of the modal analysis method to study complex flows in trailing edge cutback film cooling. Establishing the relationship between vortex and wall temperature by modal analysis method can provide a new idea for studying convective heat transfer problems.

The role of streamwise vortex in the flow of the trailing edge cutback cooling and its effect on the cooling effectiveness of the cutback surface is found.