Search results

Hybrid Reynolds-averaged Navier–Stokes large eddy simulation (RANS-LES) methods have become popular for simulation of massively separated flows at high Reynolds numbers due to their reduced computational cost and good accuracy. The current study aims to examine the performance of LES and hybrid RANS-LES model for a given grid resolution.

For better assessment and contrast of model performance, both mean and instantaneous flow fields have been investigated. For studying instantaneous flow, proper orthogonal decomposition has been used.

Current analysis shows that hybrid RANS-LES is capable of achieving similar accuracy in prediction of both mean and instantaneous flow fields at a very coarse grid as compared to LES.

Focusing mostly on the practical applications of computation, most of the attention has been given to the prediction of one-point flow statistics and little consideration has been put to two-point statistics. Here, two-point statistics has been considered using POD to investigate unsteady turbulent flow.

Applicability of two k‐ω hybrid RANS/LES and a k‐ω RANS models is studied for simulation of round impinging jets at nozzle‐plate distance H/D=2 with Reynolds number 70000, H/D=2 with Reynolds number 5000 and H/D=10 with Reynolds number 5000 (D is the nozzle exit diameter). The aim is to verify two concepts of unified hybrid RANS/LES formulations, one of DES (Detached Eddy Simulation) type and one of LNS (Limited Number Scales) type in analysis of impinging jet flow and heat transfer. The grid resolution requirements are also discussed.

The simulations are performed with two k‐ω based hybrid RANS/LES models of very different nature, one of DES type and one of LNS type, and the RANS k‐ω model. For the lower Reynolds number (5000), also dynamic Smagorinsky LES is done. Both hybrid model formulations converge to the same RANS k‐ω model in the near‐wall region and have the same Smagorinsky limit on fine isotropic grids in the LES mode of the hybrid models.

With the hybrid RANS/LES models, improved fluid flow and heat transfer results are obtained compared to RANS, in the impact region and in the developing wall‐jet region. For accurate predictions at low nozzle‐plate distance, where the impact region is in the core of the jet, it is necessary to sufficiently resolve the formation and breakup of the near‐wall vortices in the jet impingement region and the developing wall‐jet region, as these determine largely the level of fluctuating velocity and the heat transfer. This requires high grid resolution for high Reynolds number, while the grid resolution requirements stay modest for low Reynolds number.

The paper demonstrates that two formulations of hybrid RANS/LES models of different nature, one of DES type and one of LES type, lead to equivalent results. Consistency has been guaranteed in the sense that the RANS limit of both models is the same and that the LES limit on fine, isotropic, grids is the same. In the intermediate range, however, the repartition into resolved and modelled fluctuations may differ considerably.

This study aims to investigate the effect of gusts on aircraft wake vortices. Aircraft wake vortices present a potential risk to following aircraft, particularly during final approach and landing, as wake vortices may remain in the flight corridor for a long time. Wind and turbulence are key factors that influence the wake vortex evolution and the wake vortex generation in the aircraft. Flying through a gust influences the wake vortex roll-up process and its evolution. Note that vertical and lateral gusts may affect counter-rotating wake vortices differently. Both vortices influence each other by inducing a downward velocity. Disturbances may therefore lead to local vortex tilting and later to a complex three-dimensional deformation. This work uses two different hybrid Reynolds-averaged Navier–Stokes/large-eddy simulation (RANS-LES) approaches to investigate the effect of gusts on wake vortex evolution. In a one-way coupling, a pre-calculated RANS velocity field of the aircraft’s near-field is being swept through an LES domain. The effect of a sine gust on the turbulent wake is modeled by manipulating the RANS-field accordingly. As a more sophisticated approach, the concept of a two-way coupling is being presented. Here an LES solver is bi-directionally coupled with an unsteady RANS (URANS) solver, exchanging values at every physical time step of the simulation.

A one-way coupling approach of the LES code MGLET and the RANS code TAU is presented to simulate the gust effect on aircraft wake vortices. Additionally, the concept of the two-way coupling of these two codes incorporating a coupling module.

The gust effect of wake vortices subjected to a crosswind can be simulated. The vortex physics is analyzed. Unexpected behavior like fast upwind vortex decay is revealed.

The understanding of the aircraft wake vortex physics during landing provides valuable information for wake vortex advisory systems.

The effect of gust on wake vortices during and after landing has not been studied so far. The hybrid one-way coupling approach, as well as the concept of the two-way coupling, are relatively new.

This paper aims to conduct, a detailed investigation of various Reynolds averaged Navier–Stokes (RANS) models to study their performance in attached and separated flows. The turbulent flow over two airfoils, namely, National Advisory Committee for Aeronautics (NACA)-0012 and National Aeronautics and Space Administration (NASA) MS(1)-0317 with a static stall setup at a Reynolds number of 6 million, is chosen to investigate these models. The pre-stall and post-stall regions, which are in the range of angles of attack 0°–20°, are simulated.

RANS turbulence models with the Boussinesq approximation are the most commonly used cost-effective models for engineering flows. Four RANS models are considered to predict the static stall of two airfoils: Spalart–Allmaras (SA), Menter’s k – ω shear stress transport (SST), k – kL and SA-Bas Cakmakcioglu modified (BCM) transition model. All the simulations are performed on an in-house unstructured-grid compressible flow solver.

All the turbulence models considered predicted the lift and drag coefficients in good agreement with experimental data for both airfoils in the attached pre-stall region. For the NACA-0012 airfoil, all models except the SA-BCM over-predicted the stall angle by 2°, whereas SA-BCM failed to predict stall. For the NASA MS(1)-0317 airfoil, all models predicted the lift and drag coefficients accurately for attached flow. But the first three models showed even further delayed stall, whereas SA-BCM again did not predict stall.

The numerical results at high Re obtained from this work, especially that of the NASA MS(1)-0317, are new to the literature in the knowledge of the authors. This paper highlights the inability of RANS models to predict the stall phenomenon and suggests a need for improvement in modeling flow physics in near- and post-stall flows.

The present paper aims at evaluating the lattice Boltzmann method (LBM) on a high-subsonic high-pressure compressor stage at nominal regime.

The studied configuration corresponds to the H25 compressor operated in a closed-loop test rig at the von Karman Institute. Several operating points are simulated with LBM for two grids of successive refinements. A detailed analysis is performed on the time-averaged flow predicted by LBM, using a comparison with experimental and existing RANS data.

The finest grid is found to correctly predict the mean flow across the machine, as well as the influence of the rotor tip gap size. Going beyond time-averaged data, some flow analysis is performed to show the relevance of such a high-fidelity method applied to a compressor configuration. In particular, vortical structures and their evolution with the operating points are clearly highlighted. Spectral analyses finally hint at a proper prediction of tonal and broadband contents by LBM.

The application of LBM to high-speed turbomachinery flows is very recent. This paper validates one of the first LBM simulations of a high-subsonic high-pressure compressor stage.

This paper aims to study passive control techniques for transonic flow over a backward-facing step (BFS) using square-lobed trailing edges. The study investigates the efficacy of upward and downward lobe patterns, different lobe widths and deflection angles on flow separation, aiming for a deeper understanding of the flow physics behind the passive flow control system.

Large Eddy Simulation and Reynolds-averaged Navier–Stokes were used to evaluate the results of the study. The research explores the impact of upward and downward patterns of lobes on flow separation through the effects of different lobe widths and deflection angles. Numerical methods are used to analyse the behaviour of transonic flow over BFS and compared it to existing experimental results.

The square-lobed trailing edges significantly enhance the reduction of mean reattachment length by up to 80%. At Ma = 0.8, the up-downward configuration demonstrates increased effectiveness in reducing the root mean square of pressure fluctuations at a proximity of 5-step height in the wake region, with a reduction of 50%, while the flat-downward configuration proves to be more efficient in reducing the root mean square of pressure fluctuations at a proximity of 1-step height in the near wake region, achieving a reduction of 71%. Furthermore, the study shows that the up-downward configuration triggers early spanwise velocity fluctuations, whereas the standalone flat-downward configuration displays less intense crosswise velocity fluctuations within the wake region.

The findings demonstrate the effectiveness of square-lobed trailing edges as passive control techniques, showing significant implications for improving efficiency, performance and safety of the design in aerospace and industrial systems.

This paper demonstrates that the square-lobed trailing edges are effective in reducing the mean reattachment length and pressure fluctuations in transonic conditions. The study evaluates the efficacy of different configurations, deflection angles and lobe widths on flow and provides insights into the flow physics of passive flow control systems.

The flow phenomenon of particle image velocimetry has revealed the transition process of the complex multi-scale vortex between the boundary layer and mainstream region. Nonetheless, present computational fluid dynamics methods inadequately distinguish the discernable flows in detail. A multi-physical field coupling model, which was applied in rotor-stator fluid machinery (Umavathi, 2015; Syawitri et al., 2020), was put forward to ensure the identification of multi-scale vortexes and the improvement of performance prediction in torque converter.

A newly-developed multi-physical field simulation framework that coupled the scale-resolving simulation method with a dynamic modified viscosity coefficient was proposed to comparatively investigate the influence of energy exchange on thermal and flow characteristics and the description of the flow field in detail.

Regardless of whether quantitative or qualitative, its description ability on turbulence statistics, pressure-streamline, vortex structure and eddy viscosity ratio were visually experimentally and numerically analyzed. The results revealed that the modification of transmission medium viscous can identify flows more exactly between the viscous sublayer and outer boundary layer. Compared with RANS and large eddy simulation, a stress-blended eddy simulation model with a dynamic modified viscosity coefficient, which was further used to achieve blending on the stress level, can effectively solve the calculating problem of the transition region between the near-wall boundary layer and mainstream region.

This indeed provides an excellent description of the transient flow field and vortex structure in different physical flow states. Furthermore, the experimental data has proven that the maximum error of the external performance prediction was less than 4%.

An improved model was applied to simulate and analyze the flow mechanism through the evolution of vortex structures in a working chamber, to deepen the designer with a fundamental understanding on how to reduce flow losses and flow non-uniformity in manufacturing.

This paper aims to investigate numerical solution of transonic flow around NACA0012 airfoil under sinusoidal pitch oscillation. Accordingly, effects of the amplitude and frequency of oscillations on aerodynamic coefficients are evaluated and the efficiency of the turbulent models, K-ω shear-stress transport (SST), scale adaptive simulation (SAS) and delayed detached eddy simulation (DDES), in simulation of the nonlinear phenomena – i.e. the interaction between shock and boundary layer and the shock oscillations – is studied.

K-ω SST, SAS and DDES models are used as turbulence approaches. The numerical results are compared with available experimental and numerical information.

According to the results inside the buffet boundaries, the DDES turbulent model expresses results that are more appropriate; however, SAS and SST models are not efficient enough in evaluating the characteristics of nonlinear flow.

In this research study, hybrid RANS-LES turbulence model is engaged to simulate transonic flow around pitching NACA0012 airfoil, and results are compared to the SAS and Reynolds Average Navier–Stocks simulations as well as available numerical and experimental data. In addition, effects of the amplitude and frequency of oscillations on aerodynamic coefficients are evaluated in buffet regions.

The design of a space launcher requires some considerations about the unsteady loads and heat transfer occurring at the base of the structure. In particular, these phenomena are predominant during the early stage of the flight. This paper aims to evaluate the ability of the unstructured, high order finite-volume CFD solver FLUSEPA, developed by Airbus Safran Launchers, to accurately describe these phenomena.

This paper first performs a steady simulation on a base flow around a four-clustered rocket configuration. Results are compared with NASA experiments and Loci-CHEM simulations. Then, unsteady simulations of supersonic H2/air reacting mixing layer based on the experiment of Miller, Bowman and Mungal are performed. Three meshes with different cells number are used to study the impact of spatial resolution. Instantaneous and time-averaged concentrations are compared with the combined OH/acetone planar laser-induced fluorescence imaging from the experiment.

FLUSEPA satisfactorily predicts the base heat flux at the base of a four-clustered rocket configuration. NASA Loci-CHEM reactive simulations indicate that afterburning plays an important role and should not be neglected. The unsteady reactive computation of a supersonic mixing layer shows that FLUSEPA is also able to accurately predict flow structures and interactions. However, the complexity of the experiment and the lack of details concerning the facility prevents from obtaining satisfactory converged results.

This study is the first step on the development of a cost-effective method aiming at predicting unsteady loads and heat transfer on space launchers using an unsteady and reactive model for the CDF calculations. It uses original techniques such as conservative CHIMERA-like overset grids, local re-centering of fluxes and local adaptive time-stepping to reduce computational cost while being robust and accurate.

This study aims to methodologically investigate heat transfer effects on reacting flow inside a liquid-fueled, swirl-stabilized burner. Furthermore, particular attention is paid to turbulence modeling and the results of Reynolds-averaged Navier–Stokes and large eddy simulation approaches are compared in terms of velocity field and flame temperature.

Simulations consist liquid fuel distribution using Eulerian–Lagrangian approach. Flamelet-Generated Manifold combustion model, which is a mixture fraction-progress variable formulation, is used to obtain reacting flow field. Discrete ordinates method is also added for modeling radiation heat transfer effect inside the burner. As a parametric study, different thermal boundary conditions namely: adiabatic wall, constant temperature and heat transfer coefficient are applied. Because of the fact that the burner is designed for operating with different materials, the effects of burner material on heat transfer and combustion processes are investigated. Additionally, material temperatures have been calculated using 1 D method. Finally, soot particles, which are source of luminous radiation in gas turbine combustors, are calculated using Moss-Brookes model.

The results show that the flow behavior is obviously different in recirculation region for both turbulence modeling approach, and this difference causes change on flame temperature distribution, particularly in the outer recirculation zone and region close to swirler. In thermal assessment of the burner, it is predicted that material of the burner walls and the applied thermal boundary conditions have significant influence on flame temperature, wall temperature and flow field. The radiation heat transfer also makes a strong impact on combustion inside the burner; however, luminous radiation arising from soot particles is negligible for the current case.

These types of burners are widely used in research of gas turbine combustion, and it can be seen that the heat transfer effects are generally neglected or oversimplified in the literature. This parametric study provides a basic understanding and methodology of the heat transfer effects on combustion to the researchers.