No Arabic abstract
A series of Direct Numerical Simulations (DNS) of lean methane/air flames was conducted in order to investigate the enhancement of the turbulent flame speed and modifications to the reaction layer structure associated with the systematic increase of the integral scale of turbulence $l$ while the Karlovitz number and the Kolmogorov scale are kept constant. Four turbulent slot jet flames are simulated at increasing Reynolds number and up to $Re approx 22000$, defined based on the bulk velocity, slot width, and the reactants properties. The turbulent flame speed $S_T$ is evaluated locally at select streamwise locations and it is observed to increase both in the streamwise direction for each flame and across flames for increasing Reynolds number, in line with a corresponding increase of the turbulent integral scale. In particular, the turbulent flame speed $S_T$ increases exponentially with the integral scale for $l$ up to about 6 laminar flame thicknesses, while the scaling becomes a power-law for larger values of $l$. These trends cannot be ascribed completely to the increase in the flame surface, since the turbulent flame speed looses its proportionality to the flame area as the integral scale increases; in particular, it is found that the ratio of turbulent flame speed to area attains a power-law scaling $l^{0.2}$. This is caused by an overall broadening of the reaction layer for increasing integral scale, which is not associated with a corresponding decrease of the reaction rate, causing a net enhancement of the overall burning rate. This observation is significant since it suggests that a continuous increase in the size of the largest scales of turbulence might be responsible for progressively stronger modifications of the flames inner layers even if the smallest scales, i.e., the Karlovitz number, are kept constant.
Diffusive transport of mass occurs at small scales in turbulent premixed flames. As a result, multicomponent mass transport, which is often neglected in direct numerical simulations (DNS) of premixed combustion, has the potential to impact both turbulence and flame characteristics at small scales. In this study, we evaluate these impacts by examining enstrophy dynamics and the internal structure of the flame for lean premixed hydrogen-air combustion, neglecting secondary Soret and Dufour effects. We performed three-dimensional DNS of these flames by implementing the Stefan-Maxwell equations in the code NGA to represent multicomponent mass transport, and we simulated statistically planar lean premixed hydrogen-air flames using both mixture-averaged and multicomponent models. The mixture-averaged model underpredicts the peak enstrophy by up to 13% in the flame front. Comparing the enstrophy budgets of these flames, the multicomponent simulation yields larger peak magnitudes compared to the mixture-averaged simulation in the reaction zone, showing differences of 17% and 14% in the normalized stretching and viscous effects terms. In the super-adiabatic regions of the flame, the mixture-averaged model overpredicts the viscous effects by up to 13%. To assess the effect of these differences on flame structure, we reconstructed the average local internal structure of the turbulent flame through statistical analysis of the scalar gradient field. Based on this analysis, we show that large differences in viscous effects contribute to significant differences in the average local flame structure between the two models.
Implementing multicomponent diffusion models in numerical combustion studies is computationally expensive; to reduce cost, numerical simulations commonly use mixture-averaged diffusion treatments or simpler models. However, the accuracy and appropriateness of mixture-averaged diffusion has not been verified for three-dimensional, turbulent, premixed flames. In this study we evaluated the role of multicomponent mass diffusion in premixed, three-dimensional high Karlovitz-number hydrogen, n-heptane, and toluene flames, representing a range of fuel Lewis numbers. We also studied a premixed, unstable two-dimensional hydrogen flame due to the importance of diffusion effects in such cases. Our comparison of diffusion flux vectors revealed differences of 10-20% on average between the mixture-averaged and multicomponent diffusion models, and greater than 40% in regions of high flame curvature. Overall, however, the mixture-averaged model produces small differences in diffusion flux compared with global turbulent flame statistics. To evaluate the impact of these differences between the two models, we compared normalized turbulent flame speeds and conditional means of species mass fraction and source term. We found differences of 5-20% in the mean normalized turbulent flame speeds, which seem to correspond to differences of 5-10% in the peak fuel source terms. Our results motivate further study into whether the mixture-averaged diffusion model is always appropriate for DNS of premixed turbulent flames.
The excitation and further sustenance of large-scale magnetic fields in rotating astrophysical systems, including planets, stars and galaxies, is generally thought to involve a fluid magnetic dynamo effect driven by helical magnetohydrodynamic turbulence. While this scenario is appealing on general grounds, it however currently remains largely unconstrained, notably because a fundamental understanding of the nonlinear asymptotic behaviour of large-scale fluid magnetism in the astrophysically-relevant but treacherous regime of large magnetic Reynolds number $Rm$ is still lacking. We explore this problem using local high-resolution simulations of turbulent magnetohydrodynamics driven by an inhomogeneous helical forcing generating a sinusoidal profile of kinetic helicity, mimicking the hemispheric distribution of kinetic helicity in rotating turbulent fluid bodies. We identify a transition at large $Rm$ to an asymptotic nonlinear state, followed up to $Rmsimeq 3times 10^3$, characterized by an asymptotically small resistive dissipation of magnetic helicity, by its efficient spatial redistribution across the equator through turbulent fluxes driven by the hemispheric distribution of kinetic helicity, and by the presence in the tangled dynamical magnetic field of plasmoids typical of reconnection at large $Rm$.
In this video, we present the dynamics of an array of falling particles at intermediate Reynolds numbers. The film shows the vorticity plots of 3, 4, 7, 16 falling particles at $Re = 200$. We highlight the effect of parity on the falling configuration of the array. In steady state, an initially uniformly spaced array forms a convex shape when $n=3$, i.e the middle particle leads, but forms a concave shape when $n = 4$. For larger odd numbers of particles, the final state consists of a mixture of concave and convex shapes. For larger even numbers of particles, the steady state remains a concave shape. Below a threshold of initial particle spacing, particles cluster in groups of 2 to 3.
Propagation of weakly stretched spherical flames in partially pre-vaporized fuel sprays is theoretically investigated in this work. A general theory is developed to describe flame propagation speed, flame temperature, droplet evaporation onset and completion locations. The influences of liquid fuel and gas mixture properties on spherical spray flame propagation are studied. The results indicate that the spray flame propagation speed is enhanced with increased droplet mass loading and/or evaporation heat exchange coefficient (or evaporation rate). Opposite trends are found when the latent heat is high, due to strong evaporation heat absorption. Fuel vapor and temperature gradients are observed in the post-flame evaporation zone of heterogeneous flames. Evaporation completion front location considerably changes with flame radius, but the evaporation onset location varies little relative to the flame front when the flame propagates. For larger droplet loading and smaller evaporation rate, the fuel droplet tends to complete evaporation behind the flame front. Flame bifurcation occurs with high droplet mass loading under large latent heat, leading to multiplicity of flame propagation speed, droplet evaporation onset and completion fronts. The flame enhancement or weakening effects by the fuel droplet sprays are revealed by enhanced or suppressed heat and mass diffusion process in the pre-flame zone. Besides, for heterogeneous flames, heat and mass diffusion in the post-flame zone also exists. The mass diffusion for both homogeneous and heterogeneous flames is enhanced with decreased Lewis number. The magnitude of Markstein length is considerably reduced with increased droplet loading. Moreover, post-flame droplet burning behind heterogeneous flame influences the flame propagation speed and Markstein length when the liquid fuel loading is relatively low.