In this study, direct numerical simulation of the particle dispersion and turbulence modulation in a sonic transverse jet injected into a supersonic cross-flow with a Mach number of 2 was carried out with the Eulerian–Lagrangian point-particle method. One single-phase case and two particle-laden cases with different particle diameters were simulated. The jet and particle trajectories, the dispersion characteristics of particles, and the modulation effect of particles on the flow were investigated systematically. It was found that large particles primarily accumulate around shear layer structures situated on the windward side of the jet trajectory. In contrast, small particles exhibit radial transport, accessing both upstream and downstream recirculation zones. Moreover, small particles disperse extensively within the boundary layer and large-scale shear layers, evidently influenced by the streamwise vortices. The particles increase the mean wall-normal velocity near the wall in the wake region of the transverse jet, while reducing the mean streamwise and wall-normal velocities in outer regions. Particles significantly alter the flow velocity adjacent to shock fronts. In particular, the turbulent fluctuations near the windward barrel shock and bow shock are reduced, while those around the leeward barrel shock are increased. An upward displacement of the bow shock in the wall-normal direction is also observed due to particles. In the regions away from the shocks, small particles tend to amplify the Reynolds stress, while large particles attenuate the turbulent kinetic energy.

Lean hydrogen/air premixed flame flashback in a turbulent boundary layer over a flat plate is investigated using three-dimensional direct numerical simulation with detailed chemical kinetics. The upstream propagation of the flame takes place in near-wall turbulence and the interaction between the flame and the approaching reactant flow is studied. It is found that backflow regions are always present immediately upstream of flame bulges that are convex towards the reactants, confirming earlier observations. A flashback speed, including the effects of flame displacement speed and flow velocity, is introduced to quantify the flame flashback behaviour. This analysis indicates that the flashback speed is overall positive and it is considerably affected by the presence of the backflow regions. A budget analysis of the pressure transport equation is performed to explain the presence of the backflow regions. It is suggested that the positive dilatation and thermal diffusion terms near the leading edge of flame bulges are the main reasons for the pressure increase, leading to an adverse pressure gradient. The effects of the flame-induced adverse pressure gradient on the structures of the turbulent boundary layer are also investigated. It is revealed that the near-wall mean velocity and skin-friction coefficient are reduced due to the adverse pressure gradient. The coherent vortical structures of the boundary layer turbulence are lifted by the adverse pressure gradient. The analysis of the Reynolds stress component showed that the ejection event is augmented by combustion while the sweep event is attenuated, which facilitates the occurrence of flame flashback.

In the present work, direct numerical simulation of a laboratory-scale lean premixed reacting jet in cross-flow was performed to explore the flow–flame structures and turbulence–flame interactions. A jet of lean premixed ethylene–air mixtures (equivalence ratio $\phi = 0.6$) was injected into a hot vitiated cross-flow. Both non-reacting and reacting cases were simulated. It was found that the reacting jet penetrates deeper in the cross-flow with a weaker shear layer compared with the non-reacting one. The wake of the non-reacting and reacting jet is characterized by vertical vortices and recirculation zones, respectively. As for the flame structure of the reacting case, the reaction intensity varies considerably in different flame zones. The heat release rate on the leeward side is higher than that on the windward side, but lower than that of the corresponding laminar flame. The analysis of the turbulence–flame interactions of the reacting case showed that the large local Damköhler number ($Da$) related to reaction-induced dilatations results in an increased tendency of the scalar gradient to align with the most extensive strain rate, which is more evident in the regions with high heat release rate on the leeward side. Negative dilatation regions with positive tangential strain rate and negative normal strain rate are observed on the windward side. High positive dilatations appear on the flame front of the leeward side. The tangential strain rate is negatively correlated with the normal strain rate and curvature. Regions with a high local $Da$ on the windward side correspond with high positive curvature regions.

In the present study, turbulent particle-laden boundary layer combustion over a flat plate is investigated using direct numerical simulation (DNS). A two-way coupled Eulerian–Lagrangian point particle method is used for the solid phase. The effects of particle Stokes numbers, mass loadings and chemical reactions on the interactions between particles and boundary layer turbulence in the near-wall region are explored. It was found that particle heat transfer is dominant over wall heat transfer in the reacting case with heavy particles and large mass loadings, resulting in a lower fluid temperature. Particle accumulation due to the turbophoresis effect in the near-wall region is observed, which is more prominent in the cases with a large Stokes number. The turbophoresis effect is examined via the magnitude of streamwise vorticity $\varOmega _x$. It is shown that $\varOmega _x$ is attenuated by heavy particles, and the attenuation increases with increasing mass loadings. Therefore, particle wall–accumulation is less prominent in the cases with large mass loadings. Compared with the non-reacting cases, the distribution of particles is more inhomogeneous for the reacting cases, where the particles move faster due to intense reactions with increasing wall-normal distance. Finally, the flow topologies and the Reynolds stress anisotropy are examined to understand turbulence modulation by combustion and particles. It was suggested that light particles augment the vortex-dominant topologies, whereas heavy particles have an opposite effect. The anisotropy mapping of the Reynolds stress shows that the flows become more one-dimensional in the near-wall region for the cases with combustion and/or large mass loadings.

The aim of this study is to establish a fundamental understanding of the flame structure and autoignition characteristics of supercritical hydrothermal flames in three-dimensional shear-driven turbulence. The study involves direct numerical simulation of a non-premixed flame (with fuel comprising a mixture of 10 % $\textrm {H}_2$ and 90 % $\textrm {H}_2\textrm {O}$ in terms of mole fraction) at 25.0 MPa in a slot jet; detailed reaction mechanism and multispecies real-fluid properties are considered in the simulation. Qualitative transient inspection revealed that the flame undergoes a three-stage development process in the streamwise direction: sparse autoignition kernels in the upstream region, intense ignitions and establishment of a continuous flame surface in the middle region, and massive flamelets in the downstream region. Ignition kernels primarily form in the interior of large-scale shear-driven vortices featuring a low scalar dissipation rate. Probability density function (p.d.f.) analysis further confirmed that these kernels mainly form in the premixed combustion mode and on the fuel-lean side, in contrast to the authors’ previous findings concerning autoignition in a two-dimensional mixing layer. Analysis of the preignition chemistry indicator (i.e. $\textrm {H}_2\textrm {O}_{2}$ radicals) revealed that although the fuel-rich condition has a shorter homogeneous autoignition delay time, it does not exhibit any remarkable preignition chemistry or intense heat release in the upstream or middle regions because of its large-scale flow structure. A volume rendering of the dimensionless Damköhler number ($Da$) reveals the distribution of autoignition spots and propagating flames. The joint p.d.f. of the mixture fraction and $Da$ reveals the transition from sparse ignition to intense ignition and, finally, to flame propagation.

The present paper focuses on the structures and dynamics of flame edges in planar turbulent non-premixed flames bounded with a wall using direct numerical simulation (DNS). The global quenching behaviour was first examined and the flame edges were identified based on the intersections of mixture fraction and OH mass fraction iso-surfaces. For the upper branch of the planar jet flame, it is observed that the structures of flame edges change from tribrachial to monobrachial with increasing scalar dissipation rate. The flame edge speed is negatively correlated with the scalar dissipation rate in regions away from the wall, highlighting the role of turbulent mixing on the flame edge dynamics. During flame–wall interactions, the propagation speed of flame edges is mainly affected by the projection of edge flame normal in the wall-normal direction, i.e. $\boldsymbol {N}_{Z}\boldsymbol {\cdot }\boldsymbol {N}_{wall}$. In particular, the propagation speed increases with increasing $\boldsymbol {N}_{Z}\boldsymbol {\cdot }\boldsymbol {N}_{wall}$ in the near-wall region. The interactions of flame edges and turbulence bounded with a wall are characterized by the alignment between edge flame normal and principal strain rates. The normal of quenching edges has a tendency to align with the most extensive strain rate $\boldsymbol {e}_{1}$ in regions where the heat-release-induced dilatation is dominant over turbulent strain. In contrast, when the heat loss by cold wall effect is large enough to counteract the heat release induced by chemical reactions, turbulent strain is prevalent and the edge flame normal of the quenching edges preferentially aligns with the most compressive strain rate $\boldsymbol {e}_{3}$.

Three-dimensional (3-D) measurements of flame stretch are experimentally challenging. In this paper, two-dimensional (2-D) and 3-D measurements of flame stretch and turbulence–flame interactions were examined using direct numerical simulation (DNS) data of turbulent premixed flames, and models to estimate 3-D statistics of flame stretch-related quantities by correcting 2-D measurements were developed. A variety of DNS cases were simulated, including three freely propagating planar flames without a mean shear and a slot-jet flame with a mean shear. The main findings are summarized as follows. First, the mean shear mainly influences the flame orientations. However, it does not change the flame stretch and turbulence–flame interactions qualitatively. The distributions of out-of-plane angle of all cases are nearly isotropic. Second, models were proposed to approximate the 3-D statistics of flame stretch-related quantities using 2-D measurements, the performance of which was verified by comparing modelled and actual 3-D surface averages and probability density functions of tangential strain rate, curvature and displacement velocity. Third, 2-D measurements of flame stretch capture properly the trends of the 3-D results, with flame surface area being produced in low curvature regions and destroyed in highly curved regions. However, the magnitude of flame stretch was under-estimated in 2-D measurements. Finally, 2-D and 3-D turbulence–flame interactions were examined. The flame normal vector is aligned with the most compressive strain rate in both 2-D and 3-D measurements. Meanwhile, the flame normal vector is misaligned (weakly aligned) with the most extensive strain rate in 3-D (2-D) measurements, highlighting the difference in 2-D and 3-D results of turbulence–flame interactions.

In this study, direct numerical simulation of the dispersion and motion of inertial particles in a spatially developing compressible turbulent boundary layer at a Mach number of 2 is performed with the Eulerian–Lagrangian point particle method. Two cases are simulated with different particle diameters (Stokes number) but identical inflow particle numbers. Statistical characteristics and preferential accumulation of particles in the very-near-wall and wake regions are systematically investigated through conditional sampling and mechanism analysis. The results reveal that particle streaks are formed in low-speed regions near the wall because of the influence of dominating ejection events. After normalization with the local minimum particle number density, the particle number density profile reveals a self-similar feature at different streamwise positions. Compared with small particles, large particles are more significantly influenced by turbophoresis and demonstrate stronger preferential accumulation; thus, more large particles are clustered in the near-wall regions and the deviation between the mean velocities of the particle and the fluid increases. With the wall effect, both large and small particles are selectively accumulated in high-vorticity regions in the buffer layer in contrast to turbulence without walls. In comparison with incompressible wall-bounded turbulence, a new mechanism for particle preferential accumulation based on local fluid density is discovered. Large particles are located in low-density regions in the inner zones and high-density regions in the outer zones. Nevertheless, small particles remain located in regions with low fluid density, as illustrated by the mechanism analysis of particle dilatation.

Direct numerical simulations of particle-laden flows in a spatially developing turbulent thermal boundary layer over an isothermally heated wall have been performed with realistic fully developed turbulent inflow boundary conditions. To the authors’ best knowledge, this is the first time the effects of inertial solid particles on turbulent flow and heat transfer in a flat-plate turbulent boundary layer have been investigated, using a two-way coupled Eulerian–Lagrangian method. Results indicate that the presence of particles increases the mean streamwise velocity and temperature gradients of the fluid in the near-wall region. As a result, the skin-friction drag and heat transfer are significantly enhanced in the particle-laden flows with respect to the single-phase flow. The near-wall sweep and ejection motions are suppressed by the particles and hence the Reynolds shear stress and wall-normal turbulent heat flux are attenuated, which leads to reductions in the production of the turbulent kinetic energy and temperature fluctuations. In addition, the coherence and spacing of the near-wall velocity and temperature streaky structures are distinctly increased, while the turbulent vortical structures appear to be disorganized under the effect of the particles. Moreover, the intensity of the streamwise vortices decreases monotonically with increasing particle inertia.

In this paper, a systematic investigation of turbulence modulation by particles and its underlying physical mechanisms in decaying compressible isotropic turbulence is performed by using direct numerical simulations with the Eulerian–Lagrangian point-source approach. Particles interact with turbulence through two-way coupling and the initial turbulent Mach number is 1.2. Five simulations with different particle diameters (or initial Stokes numbers, $St_{0}$) are conducted while fixing both their volume fraction and particle densities. The underlying physical mechanisms responsible for turbulence modulation are analysed through investigating the particle motion in the different cases and the transport equations of turbulent kinetic energy, vorticity and dilatation, especially the two-way coupling terms. Our results show that microparticles ($St_{0}\leqslant 0.5$) augment turbulent kinetic energy and the rotational motion of fluid, critical particles ($St_{0}\approx 1.0$) enhance the rotational motion of fluid, and large particles ($St_{0}\geqslant 5.0$) attenuate turbulent kinetic energy and the rotational motion of fluid. The compressibility of the turbulence field is suppressed for all the cases, and the suppression is more significant if the Stokes number of particles is close to 1. The modifications of turbulent kinetic energy, the rotational motion and the compressibility are all related with the particle inertia and distributions, and the suppression of the compressibility is attributed to the preferential concentration and the inertia of particles.

The Brownian motion of a nanoparticle in fluid depends on the molecular forces acting on it. Because of the small size and the high frequency, it is difficult to make experimental measurements of these forces. In the present work, Brownian forces acting on a nanoparticle are numerically investigated with the molecular dynamics method. Some new phenomena are disclosed. (i) The probability distribution shows that the Brownian forces conform to the Gaussian distribution and self-similarity of the probability distribution is also found for different $1/Kn$ numbers which are characterized with the particle radius and the mean path $\unicode[STIX]{x1D706}$ of the gas molecule $(1/Kn=R/\unicode[STIX]{x1D706})$. (ii) The frequency spectrum distribution of the Brownian force is not a white noise spectrum, which is different from the assumption commonly used in Langevin model. The preferential frequency of the Brownian force is found. (iii) The size effect relating to the Brownian forces is not monotonically varying with $1/Kn=R/\unicode[STIX]{x1D706}$ and is also found. It first increases and then decreases after it reaches the maximum value at $1/Kn\approx 250$. The variation process for $1/Kn<250$ observed in the present work has not been reported in previous research to date.

Direct numerical simulations of particle-laden spatially developing turbulent boundary layers over a flat plate have been performed to investigate the effect of inertial particles on turbulence modulation, using the Eulerian–Lagrangian point-particle approach with two-way coupling. The particles are smaller than the Kolmogorov length scale of the dilute flow, and inter-particle collisions are not considered. The simulation results show that the addition of small solid particles increases the mean streamwise fluid velocity, which in turn leads to a reduction in the boundary layer integral parameters and an increase in the skin-friction drag. These effects become more pronounced as the particle Stokes number and mass loading increase. The streamwise turbulence intensity is slightly enhanced in the close vicinity of the wall but damped in the outer layer. In contrast, the Reynolds stress and the turbulence intensities in the wall-normal and spanwise directions are substantially attenuated across the entire boundary layer, and the levels of attenuation increase monotonically with both particle Stokes number and mass loading. The exchange of kinetic energy between particles and fluid indicates that particle–fluid interactions cause extra energy dissipation, which plays a crucial role in turbulence modulation.