Large Eddy Simulation Framework Research Articles

Impact of helical grooves on drag force and flow-induced noise of a cylinder under subcritical Reynolds numbers

The drag force and flow-induced noise of underwater vehicles significantly affect their hydrodynamic and stealth performance. This paper investigates the impact of helical grooves on the drag force and flow-induced noise of underwater vehicles through numerical simulations of the flow around cylinders with two types of helical grooves under various subcritical Reynolds numbers. The simulation scheme employs the large-eddy simulation framework combined with the Lighthill acoustic analogy method. The results show that the helical-groove structure can achieve reductions of up to 30% in drag and 5 dB in noise. These helical grooves have a significant effect in terms of suppressing the formation of a Karman vortex street downstream of the cylinder. Under subcritical Reynolds numbers, the drag-reduction effect of the helically grooved cylinder decreases as the number of helical grooves increases, while the noise-reduction effect increases with increasing number of helical grooves.

Large-Eddy Simulations with remeshed Vortex methods: An assessment and calibration of subgrid-scale models

This study assesses various subgrid-scale models within the framework of Large Eddy Simulation (LES) using a remeshed Vortex method (RVM). RVM is a semi-Lagrangian method discretizing the vorticity-velocity Navier–Stokes equations that has proven to be a stable and less dissipative alternative to more classical Eulerian methods. The subgrid-scale models are first tested on the well-known Taylor–Green Vortex case at Re=5000. Notably, the Variational Multiscale (VMS) variant of the Smagorinsky model and the Spectral Vanishing Viscosity (SVV) approaches emerge as the best-suited to the RVM, as they add diffusion to only the smallest resolved vorticity scales. Then, a stochastic uncertainty quantification analysis is conducted for both selected models, and the model coefficients are calibrated against direct numerical simulation. These coefficients are then applied to additional cases (different regimes, grid resolutions and test cases), showing the robustness of the calibration within the RVM-LES framework.

Liquid sheet formation and spray characterization of N-heptane spray jet from a swirl atomizer: Numerical analysis and validation

This paper analyzes the liquid atomization mechanisms of an N-heptane spray jet emerging from a swirl simplex atomizer using numerical and experimental techniques. In particular, a direct seamless coupled Volume of Fluid and Lagrangian Particle Tracking approach together with adaptive mesh refinement within the Large Eddy Simulation framework offers a suitable way to accurately simulate the complex behavior of spray atomization, spray evolution, and droplet dispersion as a whole while using manageable computational cost. The achieved simulation results are first presented in terms of qualitative properties, characteristics of liquid sheet, air core generation, flow recirculation zones, and vortex patterns. For validation purposes, the numerical results are then compared with detailed experimental data obtained by a two-component Phase Doppler Anemometry technique. The assessment includes especially droplet statistics which strongly determine subsequent possible spray combustion process and related product and species emissions. The overall reported agreement demonstrates the capability of the adopted methodology in predicting and comprehensively investigating the complex phenomena associated with a pressure swirl fuel atomizer. In particular, it is found out that the normalized number-based probability density function of droplet size fits well with a lognormal distribution.

Large eddy simulation-bivariate sectional model for non-spherical TiO2 nanoparticle synthesis in flame

Numerical simulation of nanoparticle synthesis in flame necessitates a high level of fidelity for turbulent flame, as well as delivers more detail insight into the evolution of structure and size distribution of particles undergoing various dynamic events. In the present work, a new simulation approach coupling large-eddy simulation (LES) and bivariate sectional model is proposed to explore the synthesis of nanoparticles in turbulent flame. More specifically, the nonlinear LES-partially stirred reactor model (NLES-PaSR) is utilized to simulate the turbulent flame, wherein the gradient-type structural subgrid-scale (SGS) model for turbulent flow and partially stirred reactor (PaSR) model for turbulent combustion are considered. In addition, the volume-based bivariate sectional model is thoroughly coupled with the LES model for the first time to describe the spatiotemporally resolved formation and growth of nanoparticles using particle number concentration and surface area concentration as two sets of variables. Compared directly with the LES-monodisperse (LES-Mono) and LES-univariate sectional (LES-UniSe) models, the LES-bivariate sectional (LES-BiSe) model has significant advantages in its ability to analyze and predict the size, morphology, as well as polydispersed particle size distribution (PSD) evolution of nanoparticles at a moderate computational cost simultaneously. In the LES framework, the proposed model successfully simulates the spatiotemporally resolved formation and growth of non-spherical TiO2 nanoparticle in turbulent diffusion burner. The predicted results derived from the LES-BiSe model are in good agreement with the experimental measurements. Also, the aggregate and primary particle size distribution shows highly temporal fluctuation, which mainly results from the continuous variations of particle path induced by the unsteady turbulent combustion. Furthermore, the present work demonstrates the profound dependence of particle morphological on both spatiotemporal evolution and their intrinsic size distribution. In particular, the large aggregates in the higher sections are commonly found as irregular structure with fractal-like dimensions.

Numerical investigation of supercritical combustion dynamics in a multi-element LOx–methane combustor using flamelet-generated manifold approach

The article investigates liquid oxygen (LOx)–methane supercritical combustion dynamics in a multi-element rocket-scale combustor using large eddy simulation (LES). A complex framework of real gas thermodynamics and flamelet-generated manifold (FGM) combustion model is invoked to simulate transcritical oxygen injection and supercritical methane combustion. A benchmark Mascotte chamber, rocket combustion Modeling (RCM) test case, i.e., RCM-3 (V04)/G2 test case, is used to validate the real gas FGM model in the LES framework. The validation study accurately reproduces experimental flame structure and OH concentration, demonstrating the FGM model's importance in incorporating finite rate kinetics in LOx–methane combustion. Subsequently, the numerical framework investigates a specially designed multi-element combustion chamber featuring seven bidirectional swirl coaxial injectors. The analyses capture the complex hydrodynamics and combustion dynamics associated with multiple swirl injectors operating at supercritical pressure, effectively demonstrating the initiation of transverse acoustic waves and examining the effect of local sound speed on the evolution of acoustic modes in the combustor. The dominant frequency modes shed light on understanding the role of injectors in enhanced combustor dynamics. Spectral analysis reveals the interplay of the upstream injector and chamber acoustics due to possible frequency coupling. The results also highlight the effect of fuel injection temperature on the stability of the combustor, revealing a violent dynamic activity for lower fuel injection temperature associated with the longitudinal acoustic mode of the combustor. The investigation appropriately reproduces self-sustained limit-cycle oscillations at lower fuel injection temperatures and corroborates the conventional understanding of combustor instability.

Large-Eddy Simulations of the Tropical Cyclone Boundary Layer at Landfall in an Idealized Urban Environment

Abstract Populated urban areas along many coastal regions are vulnerable to landfalling tropical cyclones (TCs). To the detriment of surface parameterizations in mesoscale models, the complexities of turbulence at high TC wind speeds in urban canopies are presently poorly understood. Thus, this study explores the impacts of urban morphology on TC-strength winds and boundary layer turbulence in landfalling TCs. To better quantify how urban structures interact with TC winds, large-eddy simulations (LESs) are conducted with the Cloud Model 1 (CM1). This implementation of CM1 includes immersed boundary conditions (IBCs) to represent buildings and eddy recycling to maintain realistic turbulent flow perturbations. Within the IBCs, an idealized coastal city with varying scales is introduced. TC winds impinge perpendicularly to the urbanized coastline. Numerical experiments show that buildings generate distinct, intricate flow patterns that vary significantly as the city structure is varied. Urban IBCs produce much stronger turbulent kinetic energy than is produced by conventional surface parameterizations. Strong effective eddy viscosity due to resolved eddy mixing is displayed in the wake of buildings within the urban canopy, while deep and enhanced effective eddy viscosity is present downstream. Such effects are not seen in a comparison LES using a simple surface parameterization with high roughness values. Wind tunneling effects in streamwise canyons enhance pedestrian-level winds well beyond what is possible without buildings. In the arena of regional mesoscale modeling, this type of LES framework with IBCs can be used to improve parameters in surface and boundary layer schemes to more accurately represent the drag coefficient and the eddy viscosity in landfalling TC boundary layers. Significance Statement This is among the first large-eddy simulation model studies to examine the impacts of tropical cyclone–like winds around explicitly resolved buildings. This work is a step forward in bridging the gap between engineering studies that use computational fluid dynamics models or laboratory experiments for flow through cities and mesoscale model simulations of landfalling tropical cyclones that use surface parameterizations specialized for urban land use.

LES investigation of soot formation in a turbulent non-premixed jet flame with sectional method and FGM chemistry

This study proposes a detailed soot modeling framework for large-eddy simulation (LES) to accurately predict soot formation and particle size distributions (PSD) in turbulent reacting flows. The framework incorporates Flamelet Generated Manifold (FGM) chemistry and a soot model based on the discrete sectional method (DSM) to predict both qualitative and quantitative sooting behavior while keeping the computational cost affordable. Two elementary modeling strategies are considered in the LES formalism for describing soot formation rates. These strategies rely on an a-priori tabulation of soot formation rates and their run-time computation. The LES formalism is applied to the simulations of a well-characterized, non-premixed, turbulent jet flame. A comparative analysis of strategies employed for filtered soot source term treatment is conducted to investigate their impact on the prediction of soot quantities and the evolution of soot PSDs. The LES results for the gas phase and soot phase are compared against the available experimental data. A good prediction of soot evolution is achieved with the two methodologies. The tabulation of soot formation rates leads to a significant reduction in computational cost compared to the model based on their explicit runtime computation. The LES results reveal that the modeling of filtered soot source terms has a significant impact on the quantitative prediction of soot formation. The possible reasons for the observed differences in the soot prediction are discussed. The run-time computation-based model provides a more consistent treatment of the non-linear interactions between the gas and soot phases in soot source terms compared to the tabulated soot chemistry approach. On the other hand, the tabulated soot chemistry model is an interesting and efficient modeling approach for predicting soot formation in turbulent conditions. Overall, both approaches have their strengths and limitations, and the choice of approach may depend on the specific needs of the application.

A priori test of perfectly stirred reactor approach to evaluating mean fuel consumption and heat release rates in highly turbulent premixed flames

Unsteady three-dimensional Direct Numerical Simulation (DNS) data obtained earlier by Dave et al. ( J Fluid Mech 2020; 884: A46) from a statistically planar and one-dimensional, highly turbulent, moderately lean hydrogen-air flame propagating in a box are processed to perform a priori test of perfectly stirred reactor model. The test aims in particularly at estimating mesh resolution (or filter width within large eddy simulation framework) required to neglect variations in the temperature and mixture composition within a computational cell when evaluating mean (or filtered) fuel consumption and heat release rates. For this purpose, fuel consumption and heat release rates sampled directly from the DNS data and averaged over a cube of width [Formula: see text] are compared with fuel consumption and heat release rates calculated using the temperature and species concentrations averaged over the same cube. Moreover, turbulent burning velocities computed by integrating the former and latter rates are compared with one another. A ratio of [Formula: see text] to a laminar flame thickness [Formula: see text] is varied from 0.44 to 1.8. The obtained results indicate that the tested simple approach performs reasonably well (poor) if [Formula: see text] ([Formula: see text], respectively). This result is further supported by directly filtering fuel consumption rate in a laminar premixed flame. The values of the thickness [Formula: see text], calculated using detail chemical mechanisms for different fuels under elevated temperatures and pressures associated with combustion in piston engines, indicate that it is difficult to satisfy the constraint of [Formula: see text] in contemporary unsteady multidimensional numerical simulations of turbulent burning in such engines.

Large-Eddy Simulation of Low-Reynolds-Number Flow Around Partially Porous Airfoils

This paper uses the large-eddy simulation framework in OpenFOAM to investigate the flow characteristics around two-dimensional partially porous airfoils with symmetrical flat geometries, which were used in a previous experiment. The effects of a porous medium near the trailing edge are investigated at a chord Reynolds number of 350,000. Two porous media are modeled to investigate how the porous medium influences the flow around the airfoil and aerodynamic performance: One is a normal porous medium, defined as a homogeneous and isotropic material with identical properties in all directions. The other medium consists of straight holes perpendicular to the airfoil surface, allowing fluid to pass through unidirectionally. It is observed that drag and lift forces vary, depending on the internal structure of the porous medium. The size of the laminar separation bubbles also varies, depending on the parameters of the porous medium, such as the internal structure, the length of the porous section, and the porosity. It is shown that an airfoil with a low-porosity unidirectional porous medium increases the lift-to-drag ratio at high angles of attack as compared with the nonporous airfoil.

Toward the use of LES for industrial complex geometries. Part I: automatic mesh definition

With the constant increase of computational power for the past years, Computational Fluid Dynamics (CFD) has become an essential part of the design in complex industrial processes. In this context, among the scale resolving numerical methods, Large-Eddy Simulation (LES) has become a valuable tool for the simulation of complex unsteady flows. To generalise the industrial use of LES, two main limitations are identified. First, the generation of a proper mesh can be a difficult task, which often relies on user-experience. Secondly, the ‘time-to-solution’ associated with the LES approach can be prohibitive in an industrial context. In this work, these two challenges are addressed in two parts. In this Part I, an automatic procedure for mesh definition is proposed, whereas the Part II is devoted to numerical technique to reduce the LES ‘time-to-solution’. The main goal of these works is then to develop an accurate LES strategy at an optimised computational cost. Concerning the mesh definition, because LES is based on separation between resolved and modelled subgrid-scales, the quality of the computed solution is then directly linked to the quality of the mesh. However, the definition of an adequate mesh is still an issue when LES is used to predict the flow in an industrial complex geometry without a priori knowledge of the flow dynamics. This first part presents a user-independent approach for both the generation of an initial mesh and the convergence of the mesh in the LES framework. An automatic mesh convergence strategy is proposed to ensure LES accuracy. This strategy is built to guarantee a mesh-independent mean field kinetic energy budget. The mean field kinetic energy is indeed expected to be mesh independent since only turbulent scales should be unresolved in LES. The approach is validated on canonical cases, a turbulent round jet and a turbulent pipe flow. Finally, the PRECCINSTA swirl burner is considered as a representative case of complex geometry. First, an algorithm for the generation of an unstructured mesh from a STL file is proposed to generate a coarse initial mesh, before applying the mesh convergence procedure. The overall strategy including automatic first mesh generation and its automatic adaptation paves the way to use LES approach as a decision support tool for various applications, provided that the ‘time-to-solution’ is compatible with the applications constraint. A second paper, referred as Part II, is devoted to the reduction of this time.

Application of Vortex Identification Methods in Vertical Slit Fishways

The reproduction and survival of fish are often negatively affected by the construction of dams and other hydroelectric projects, which cut off their migratory routes. Building effective fish passage facilities that allow fish to pass through dams smoothly alleviates the negative impact of hydroelectric projects on the ecological environment, thus protecting the diversity of aquatic species and preventing the extinction of indigenous fish. Vertical slit fishways are highly effective, but turbulence inside the fishway pools directly affects fish passage. In this study, the large-eddy simulation framework is used to capture the vortex characteristics in the interior of vertical slit fishway pools, and the volume of fluid method is applied to simulate the free surface. The independence of the grid is assessed by the large-eddy simulation quality index, and the simulation results are compared with experimental acoustic Doppler velocimetry data. This work characterizes the vortex flow field inside the vertical slit fishway using the Q-criterion, Omega method, and Liutex vortex identification method. The results show that the vortex structure inside the fishway pool has obvious three-dimensional characteristics and vortex structure varies within the different fishway pool chambers. The analysis and comparison of the three different vortex identification methods show that the vortex structure captured by the Liutex method is more consistent with the actual motion pattern of the fishway water flow.

Large eddy simulation of turbulent wake flow around a marine propeller under the influence of incident waves

In practical applications, propellers often operate beneath incident waves. To investigate the effects of waves on turbulent propeller wakes, a comparative study of the instantaneous flow fields and turbulent statistics of propeller wakes under open-water condition and with incident waves is conducted through the large-eddy simulation framework. The simulations are performed on a Cartesian grid, with the air–water interface captured using a coupled level-set and volume-of-fluid method. The complex geometry of the propeller is captured using an immersed boundary method. The results show that, near the water surface, the existence of incident waves accelerates the wake flow in the axial direction. Both the mean axial velocity and turbulent kinetic energy exhibit phase-lag behavior in the presence of incident waves. Waves increase the instability of the flow in the far-wake region, resulting in the rising of turbulent kinetic energy. Further analysis of the power spectral density shows that the velocity fluctuations gain energy from the waves at low frequencies, and nonlinear interactions between wave-induced motion and turbulent fluctuations transport energy from low to high frequencies as the wake flows downstream. The present study shows the potential of the large-eddy simulation framework to provide engineering guidance and a theoretical basis for the design and operation of propellers in wave environments.

How Fast is Fast Enough? Industry Perspectives on the Use of Large-eddy Simulation in Wind Energy

The use of graphics processing units (GPUs) has facilitated unprecedented performance gains for computational fluids dynamics in recent years. In many industries this has enabled the integration of large-eddy simulation (LES) in the engineering practice. Flow modelling in the wind industry though still primarily relies on models with significantly lower fidelity. This paper seeks to investigate the reasons why wind energy applications of LES are still an exception in the industrial practice. On that account, we present a survey among industry experts on the matter. The survey shows that the large runtimes and computational costs of LES are still seen as a main obstacle. However, other reasons such as a lack of expertise and user experience, the need for more validation, and lacking trust in the potential benefits of LES reveal that computational efficiency is not the only concern. Lastly, we present an exemplary simulation of a generic offshore wind farm using a GPU-resident Lattice Boltzmann LES framework. The example shows that the runtime requirements stated by a large part of the respondents can already now be fulfilled with reasonable hardware effort.

Experimental and mechanistic study of dispersed micrometer-sized particle resuspension in a square straight duct with rough walls

The resuspension of graphite dust is an important phenomenon in the release of radioactivity and the safety of nuclear reactors during severe accidents. In this study, a visualization experimental platform is constructed to study effects of particle size, flow velocity, and wall roughness on the resuspension characteristics of graphite particles. A statistical model of particle resuspension applicable to monolayer dispersed particles is developed based on the moment equilibrium of the particles and the flow field characteristics, as calculated by the large-eddy simulation framework. The results show that particle resuspension can be divided into short- and long-term resuspension stages. Most particle resuspension occurs during the short-term stage. With increases in flow velocity and particle diameter, the aerodynamic or adhesion force acting on the particles increases, and corresponding particle resuspension fraction increases. The influence of rough walls on particle resuspension is related to both the force on the particles and the arm ratio between the wall morphology and the particle diameter. A comparison with the experimental results demonstrates that the particle resuspension model developed in this study accurately predicts the impact of flow velocity, particle size, and wall roughness on particle resuspension.

CABARET on rotating meshes

The Compact Accurately Boundary-Adjusting high-REsolution Technique (CABARET) method for 3D gas dynamics equations is generalised to rotating hexahedral mesh zones using the sliding interface approach, which extends applications of the original algorithm to a wide range of unsteady rotor flow problems. Important properties of the suggested extension of the CABARET scheme include the preservation of good linear dispersion and dissipation properties of the original CABARET method on fixed grids. To strictly preserve flux conservation at the sliding interface, the so-called supermesh approach is implemented that uses the projection of the contact faces on both sides of the sliding interface onto an intermediate virtual surface. The implementation is based on applying the sub-cell and sub-face solution reconstruction on the sliding interface of hexahedral meshes rather than deforming unstructured meshes. The sub-faces are defined by constructing projection overlap polygons using the Sutherland-Hodgman clipping algorithm. A local characteristic flux splitting is applied at each overlap polygon face, while updating the local stencil connectivity table at each time step. To further improve the dispersion properties of the CABARET grid at small Courant–Friedrichs–Lewy numbers, the dispersion-improved version of the CABARET method previously developed on rectangular Cartesian meshes is extended to non-uniform curvilinear meshes typical of open-rotor applications. Numerical examples are provided for the analytical problem of normal acoustic wave propagation through a rotating cylindrical grid zone embedded in an outer stationary grid. The error decay rates are calculated and compared with the theoretical order of CABARET convergence. In addition, the same algorithm accelerated on Graphics Processing Units (GPUs) is applied for simulations the flow around a generic two-bladed rotor in hover in the Wall Model Large Eddy Simulation framework. For the latter case, flow solutions are analysed, and their numerical grid sensitivity is discussed.

Impacts of an aerosol layer on a midlatitude continental system of cumulus clouds: how do these impacts depend on the vertical location of the aerosol layer?

Abstract. Effects of an aerosol layer on warm cumulus clouds in the Korean Peninsula when the layer is above or around the cloud tops in the free atmosphere are compared to effects when the layer is around or below the cloud bases in the planetary boundary layer (PBL). For this comparison, simulations are performed using the large-eddy simulation framework. When the aerosol layer is in the PBL, aerosols absorb solar radiation and radiatively heat up air enough to induce greater instability, stronger updrafts and more cloud mass than when the layer is in the free atmosphere. Hence, there is a variation of cloud mass with the location (or altitude) of the aerosol layer. It is found that this variation of cloud mass is reduced as aerosol concentrations in the layer decrease or aerosol impacts on radiation are absent. The transportation of aerosols by updrafts reduces aerosol concentrations in the PBL. This in turn reduces the aerosol radiative heating, updraft intensity and cloud mass.

Coupling turbulent flow with blade aeroelastics and control modules in large-eddy simulation of utility-scale wind turbines

We present a large-eddy simulation framework capable of control co-design of large wind turbines, coupling the turbulent flow environment with blade aeroelastics and turbine controllers. The geometry and aerodynamics of the rotor blades and the turbine nacelle are parameterized using an actuator surface model. The baseline collective pitch control and individual pitch control (IPC) algorithms, consisting of a single-input, single-output proportional–integral controller and two integral controllers, respectively, are incorporated into the simulation framework. Furthermore, a second-order model based on the Euler–Bernoulli beam theory is implemented to describe the blade deformation. Simulations are carried out to investigate the impact of collective and individual pitch control strategies on the deflection of turbine blades. Our results show that the IPC reduces the blade tip deflection fluctuations in the out-of-plane direction, while the fluctuations of the blade tip deflection along the in-plane direction are barely affected by the IPC. Furthermore, the blade out-of-plane deformation fluctuation is underestimated by the one-way coupling approach compared to the two-way coupling approach. The findings of this study reveal the importance of advanced control systems in reducing the dynamic loads on wind turbine blades and underscore the potential of control co-design to reduce the levelized cost of wind energy.

Toward control co-design of utility-scale wind turbines: Collective vs. individual blade pitch control

A large-eddy simulation framework has been coupled with controller modules to systematically investigate the impacts of collective (CPC) and individual (IPC) pitch control strategies on utility-scale wind turbine energy production and fatigue loads. Wind turbine components were parameterized using an actuator surface model to simulate the rotor blades and the turbine nacelle. The baseline CPC and IPC algorithms, consisting of single-input single-output proportional–integral controllers and two integral controllers, respectively, were incorporated into the numerical framework. A series of simulations were carried out to investigate the relative performance of the two controllers under various turbulent inflow conditions, spanning hub-height velocities of 7 to 14 m/s. The numerical simulation results of this study showed that, in comparison to the CPC, the IPC controller could successfully reduce the damage equivalent loads of utility-scale turbines at regions 2 and 3 of turbine operation by about 3% and 40%, respectively, without any penalty on the power production of the turbine. It was also shown that, despite its minor impact on the turbulence kinetic energy of the wake, the IPC controller did not influence the recovery of the turbine wake.

A turbulent combustion model for soot formation at the LES subgrid-scale using virtual chemistry approach

Designing efficient and clean combustion devices to reduce soot formation in practical applications is a major concern nowadays. This could be achieved by performing predictive simulations in turbulent sooting combustors. Large Eddy Simulation (LES) is an attractive approach to predict flame/turbulence interactions at small scales including soot phenomena. However, because of the peculiarly complex nature of soot production phenomena, turbulent sooting flames remain very difficult to simulate and are computationally expensive. Thus, reduced gas-phase kinetics and soot models should be included in numerical simulations. The reduced virtual chemistry approach is a great alternative to address this problem since it introduces an optimized and flexible strategy to tackle multi-mode turbulent combustion and pollutants prediction. Also, the virtual chemistry approach has recently demonstrated its capability in reproducing soot phenomena in laminar flames including radiative heat transfer. In the present work, the virtual chemistry methodology is extended to resolve turbulent sooting flames using LES. Based on analysis of flame and soot chemical time scales along with turbulent time scale, an original hybrid turbulent combustion model formulation is proposed, which degenerates to the flamelet regime in the flame front region and to a PSR formalism in the post-flame zone, characterized by slow soot chemistry. The subgrid-scale model for turbulence-chemistry-soot interactions is then implemented in an LES framework. The model is challenged in the turbulent non-premixed sooting flame from Sandia Laboratory. Results indicate that the LES subgrid-scale model has a significant impact on soot formation and confirms previous observations from the literature in LES of non-premixed turbulent sooting flames. Additionally, results show that the virtual chemistry can reproduce flame-soot characteristics of a turbulent sooting jet flame and are comparable to other state-of-the-art works.

Large-eddy simulation of Cambridge-Sandia stratified flames under high swirl

The simulation of turbulent stratified swirling flame has received much attention recently owing to its importance for research on power generation and transportation devices. In this study, we have simulated the Cambridge-Sandia flames under the high swirl flow, i.e., premixed flame (SwB3), moderately stratified flame (SwB7), and highly stratified flame (SwB11) using a large-eddy simulation (LES) framework that integrates the mixed gradient-type dynamic structure model, the modulated gradient model (MGM), and the partially-stirred-reactor combustion approach coupled with a reduced (from GRI-Mech) methane oxidation mechanism. The current LES results are assessed through a comparison with experimental data and reference LES results and show good agreement with the experimental data and an evident improvement compared to the reference results. The adopted LES framework successfully captures the secondary recirculation zones in SwB7 and SwB11. Moreover, the temperature above the central bluff body is predicted well in the present work by setting the temperature profiles on the bluff body according to the experimental data. The MGM predictions indicate that the SGS flux is not directly related to the strain rate in the highly swirling and stratified flame. The stratification promotes the reactions in the reaction zones at the height of 10 mm and shows an evident back-supported combustion mode in the reaction zone. At the height of 30 mm, the effects of the stratification are beginning to fade. At the height of 50 mm, the effects of stratification are no longer prominent.