Effect of Moisture Content on High Strain Rate Compressibility and Particle Breakage in Loose Sand

Soil-filled gabion structures are widely used to protect against the effects of blast and fragmentation. It is known that moisture content significantly affects the capability of such structures, but the behaviour of partially-saturated soils is not well characterised at the strain rates and stresses experienced in these events. In particular, little data is available for loose soils, whose compaction behaviour can have a substantial impact on structural stability and ballistic performance. This paper describes the use of split Hopkinson pressure bar experiments to characterise the pre- and post-saturation compressibility of a loose quartz sand at moisture contents of up to 15.0%. In contrast to dense soils, increases in moisture content between 0.0% and 7.5% led to a decrease in the stiffness of the sand. Above 7.5% moisture content, specimens reached full saturation during the experiment: the additional water had no further effect on the pre-saturation stiffness, but post-saturation behaviour was dominated by the stiffness of the pore water. Full saturation occurred at lower dry densities as moisture content increased, leading to a decrease in particle breakage.     

Direct Numerical Simulation of Head-On Quenching of Statistically Planar Turbulent Premixed Methane-Air Flames Using a Detailed Chemical Mechanism

A three-dimensional compressible Direct Numerical Simulation (DNS) analysis has been carried out for head-on quenching of a statistically planar stoichiometric methane-air flame by an isothermal inert wall. A multi-step chemical mechanism for methane-air combustion is used for the purpose of detailed chemistry DNS. For head-on quenching of stoichiometric methane-air flames, the mass fractions of major reactant species such as methane and oxygen tend to vanish at the wall during flame quenching. The absence of \(\text {OH}\) at the wall gives rise to accumulation of carbon monoxide during flame quenching because \(\text {CO}\) cannot be oxidised anymore. Furthermore, it has been found that low-temperature reactions give rise to accumulation of \(\text {HO}_{2}\) and \(\mathrm {H}_{2}\mathrm {O}_{2}\) at the wall during flame quenching. Moreover, these low temperature reactions are responsible for non-zero heat release rate at the wall during flame-wall interaction. In order to perform an in-depth comparison between simple and detailed chemistry DNS results, a corresponding simulation has been carried out for the same turbulence parameters for a representative single-step Arrhenius type irreversible chemical mechanism. In the corresponding simple chemistry simulation, heat release rate vanishes once the flame reaches a threshold distance from the wall. The distributions of reaction progress variable c and non-dimensional temperature T are found to be identical to each other away from the wall for the simple chemistry simulation but this equality does not hold during head-on quenching. The inequality between c (defined based on \(\text {CH}_{4}\) mass fraction) and T holds both away from and close to the wall for the detailed chemistry simulation but it becomes particularly prominent in the near-wall region. The temporal evolutions of wall heat flux and wall Peclet number (i.e. normalised wall-normal distance of \(T = 0.9\) isosurface) for both simple and detailed chemistry laminar and turbulent cases have been found to be qualitatively similar. However, small differences have been observed in the numerical values of the maximum normalised wall heat flux magnitude \(\left ({\Phi }_{\max } \right )_{\mathrm {L}}\) and the minimum Peclet number \((Pe_{\min })_{\mathrm {L}}\) obtained from simple and detailed chemistry based laminar head-on quenching calculations. Detailed explanations have been provided for the observed differences in behaviours of \(\left ({\Phi }_{\max }\right )_{\mathrm {L}}\) and \((Pe_{\min })_{\mathrm {L}}\). The usual Flame Surface Density (FSD) and scalar dissipation rate (SDR) based reaction rate closures do not adequately predict the mean reaction rate of reaction progress variable in the near-wall region for both simple and detailed chemistry simulations. It has been found that recently proposed FSD and SDR based reaction rate closures based on a-priori DNS analysis of simple chemistry data perform satisfactorily also for the detailed chemistry case both away from and close to the wall without any adjustment to the model parameters.     

Effects of Lewis Number on Head on Quenching of Turbulent Premixed Flames: A Direct Numerical Simulation Analysis

The head on quenching of statistically planar turbulent premixed flames by an isothermal inert wall has been analysed using three-dimensional Direct Numerical Simulation (DNS) data for different values of global Lewis number Le(0.8, 1.0 and 1.2) and turbulent Reynolds number Re_t. The statistics of head on quenching have been analysed in terms of the wall Peclet number Pe (i.e. distance of the flame from the wall normalised by the Zel’dovich flame thickness) and the normalised wall heat flux Φ. It has been found that the maximum (minimum) value of Φ(Pe) for the turbulent Le=0.8 cases are greater (smaller) than the corresponding laminar value, whereas both Pe and Φ in turbulent cases remain comparable to the corresponding laminar values for Le=1.0 and 1.2. Detailed physical explanations are provided for the observed Le dependences of Pe and Φ. The existing closure of mean reaction rate \(\overline {\dot {\omega }}\) using the scalar dissipation rate (SDR) in the near wall region has been assessed based on a-priori analysis of DNS data and modifications to the existing closures of mean reaction rate and SDR have been suggested to account for the wall effects in such a manner that the modified closures perform well both near to and away from the wall.     

Does Density Ratio Significantly Affect Turbulent Flame Speed?

In order to experimentally study whether or not the density ratio σ substantially affects flame displacement speed at low and moderate turbulent intensities, two stoichiometric methane/oxygen/nitrogen mixtures characterized by the same laminar flame speed S _L = 0.36 m/s, but substantially different σ were designed using (i) preheating from T _u = 298 to 423 K in order to increase S _L , but to decrease σ, and (ii) dilution with nitrogen in order to further decrease σ and to reduce S _L back to the initial value. As a result, the density ratio was reduced from 7.52 to 4.95. In both reference and preheated/diluted cases, direct images of statistically spherical laminar and turbulent flames that expanded after spark ignition in the center of a large 3D cruciform burner were recorded and processed in order to evaluate the mean flame radius \(\bar {R}_{f}\left (t \right )\) and flame displacement speed \(S_{t}=\sigma ^{-1}{d\bar {R}_{f}} \left / \right . {dt}\) with respect to unburned gas. The use of two counter-rotating fans and perforated plates for near-isotropic turbulence generation allowed us to vary the rms turbulent velocity \(u^{\prime }\) by changing the fan frequency. In this study, \(u^{\prime }\) was varied from 0.14 to 1.39 m/s. For each set of initial conditions (two different mixture compositions, two different temperatures T _u , and six different \(u^{\prime })\), five (respectively, three) statistically equivalent runs were performed in turbulent (respectively, laminar) environment. The obtained experimental data do not show any significant effect of the density ratio on S _t . Moreover, the flame displacement speeds measured at u′/S _L = 0.4 are close to the laminar flame speeds in all investigated cases. These results imply, in particular, a minor effect of the density ratio on flame displacement speed in spark ignition engines and support simulations of the engine combustion using models that (i) do not allow for effects of the density ratio on S _t and (ii) have been validated against experimental data obtained under the room conditions, i.e. at higher σ.     

Analysis of the Combined Modelling of Sub-grid Transport and Filtered Flame Propagation for Premixed Turbulent Combustion

Flame surface density (FSD) based reaction rate closure is an important methodology of turbulent premixed flame modelling in the context of Large Eddy Simulations (LES). The transport equation for the Favre-filtered reaction progress variable needs closure of the filtered reaction diffusion imbalance (FRDI) term (i.e. filtered value of combined reaction rate and molecular diffusion rate) and the sub-grid scalar flux (SGSF). A-priori analysis of the FRDI and SGSF terms has in the past revealed advantages and disadvantages of the specific modelling attempts. However, it is important to understand the interaction of the FRDI and SGSF closures for a successful implementation of the FSD based closure. Furthermore, it is not known a-priori if the combination of the best SGSF model with the best FRDI model results in the most suitable overall modelling strategy. In order to address this question, a variety of SGSF models is analysed in this work together with one well established and one recent FRDI closure based on a-priori analysis. It is found that the success of the combined FRDI and SGSF closures depends on subtle details like the co-variances of the FRDI and SGSF terms. It is demonstrated that the gradient hypothesis model is not very successful in representing the SGSF term. However the gradient hypothesis provides satisfactory performance in combination of a recently proposed FRDI closure, whereas unsatisfactory results are obtained when used in combination with another existing closure, which was shown to predict the FRDI term satisfactorily in several previous analyses.     

Measurement of local forcing on a turbulent boundary layer using PIV

An experimental study was carried out to investigate the effect of periodic blowing and suction on a turbulent boundary layer. Particle image velocimetry (PIV) was used to probe the characteristics of the flow. Local forcing was introduced to the boundary layer via a sinusoidally-oscillating jet issuing from a thin spanwise slot. Three forcing frequencies (f⁺=0.44, 0.66 and 0.88) with a fixed forcing amplitude (A⁺=0.6) were employed at Re _θ =690. The effect of three different forcing angles (α=60°, 90° and l20°) was investigated under a fixed forcing frequency (f⁺=0.088). The PIV results showed that the wall-region velocity decreases on imposition of the local forcing. Inspection of the phase-averaged velocity profiles revealed that spanwise large-scale vortices are generated downstream of the slot and persist farther downstream. The highest reduction in skin friction was achieved at the highest forcing frequency (f⁺=0.088) and a forcing angle of α=120°. The spatial fraction of the vortices was examined to analyze the skin friction reduction.     

Numerical Study of Turbulent Jet Ignition in a Lean Premixed Configuration

Direct numerical simulations (DNS) of a hot combustion product jet interacting with a lean premixed hydrogen-air coflow are conducted to fundamentally investigate turbulent jet ignition (TJI) in a three-dimensional configuration. TJI is an efficient method for initiating and controlling combustion in ultra-lean combustion systems. Fully compressible gas dynamics and species equations are solved with high order finite difference methods. The hydrogen-air reaction is simulated with a reliable detailed chemical kinetics mechanism. The physical processes involved in the TJI-assisted combustion are investigated by considering the flame heat release, temperature, species concentrations, vorticity, and Baroclinc torque. The complex turbulent flame and flow structures are delineated in three main: i) hot product jet, ii) burned-mixed, and iii) flame zones. In the TJI-assisted combustion, the flow structures and the flame features such as flame speed, temperature, and species distribution are found to be quite different than those in “standard” turbulent premixed combustion due to the existence of a high energy turbulent hot product jet. The flow structures and statistics are also found to be different than those normally seen in non-isothermal non-reacting jets.     

Heat Release Imaging in Turbulent Premixed Ethylene-Air Flames Near Blow-off

The focus of this work is to visualise the regions of CH₂O and heat release (HR) of an unconfined turbulent premixed bluff body stabilised ethylene-air flame at conditions approaching lean blow-off using simultaneous imaging of OH- and CH₂O-PLIF. The HR regions are estimated from the product of the OH and CH₂O profiles. At conditions near blow-off, wide regions of CH₂O are observed inside the recirculation zone (RZ). The presence of CH₂O and HR inside the RZ is observed to follow fragmentation of the downstream flame parts near the top of the RZ. The presence of wide regions void of both OH and CH₂O inside the RZ at conditions very close to blow-off indicates the possible entrainment of un-reacted gases into the RZ. The behaviour of the lean ethylene-air flame with Lewis number (Le) greater than 1 is compared to that of a lean methane-air flame with Le of approximately 1. For both fuels, qualitatively similar observations of flame fragmentation downstream followed by build-up of CH₂O and HR inside the RZ are observed at conditions near lean blow-off. Also, a similar trend of flame front curvature conditioned on HR was observed for both the ethylene-air and methane-air flames, where the magnitude of HR was observed to increase with the absolute value of curvature.     

Pre-Chamber Ignition Mechanism: Simulations of Transient Autoignition in a Mixing Layer Between Reactants and Partially-Burnt Products

The structure of autoignition in a mixing layer between fully-burnt or partially-burnt combustion products from a methane-air flame at ? = 0.85 and a methane-air mixture of a leaner equivalence ratio has been studied with transient diffusion flamelet calculations. This configuration is relevant to scavenged pre-chamber natural-gas engines, where the turbulent jet ejected from the pre-chamber may be quenched or may be composed of fully-burnt products. The degree of reaction in the jet fluid is described by a progress variable c (c = taking values 0.5, 0.8, and 1.0) and the mixing by a mixture fraction ξ (ξ = 1 in the jet fluid and 0 in the CH₄-air mixture to be ignited). At high scalar dissipation rates, N₀, ignition does not occur and a chemically-frozen steady-state condition emerges at long times. At scalar dissipation rates below a critical value, ignition occurs at a time that increases with N₀. The flame reaches the ξ = 0 boundary at a finite time that decreases with N₀. The results help identify overall timescales of the jet-ignition problem and suggest a methodology by which estimates of ignition times in real engines may be made.     

Experimental and Theoretical Sensitivity Analysis of Turbulent Flow Past a Square Cylinder

The Siemens SGT-800 3^rd generation DLE burner fitted to an atmospheric combustion rig has been numerically investigated. Pure methane and methane enriched by 80 vol% hydrogen flames have been considered. A URANS (Unsteady Reynolds Averaged Navier-Stokes) approach was used in this study along with the k ? ω SST and the k ? ω SST-SAS models for the turbulence transport. The chemistry is coupled to the turbulent flow simulations by the use of a laminar flamelet library combined with a presumed PDF. The effect of the mesh density in the mixing and the flame region and the effect of the turbulence model and reaction rate model constant are first investigated for the methane/air flame case. The results from the k ? ω SST-SAS along with flamelet libraries are shown to be in excellent agreement with experimental data, whereas the k ? ω SST model is too dissipative and cannot capture the unsteady motion of the flame. The k ? ω SST-SAS model is used for simulation of the 80 vol% hydrogen enriched flame case without further adjusting the model constants. The global features of the hydrogen enrichment are very well captured in the simulations using the SST-SAS model. With the hydrogen enrichment the time averaged flame front location moves upstream towards the burner exit nozzle. The results are consistent with the experimental observations. The model captures the three dominant low frequency unsteady motion observed in the experiments, indicating that the URANS/LES hybrid model indeed is capable of capturing complex, time dependent, features such as an interaction between a PVC and the flame front.     

Damköhler-Shelkin Paradox in the Theory of Turbulent Flame Propagation,and a Concept of the Premixed Flame at the Intermediate Asymptotic Stage

We formulated a paradox in the theory of turbulent premixed flame in the flamelet regime: discrepancy between the Damköhler (1940) and Shelkin (1943) estimate of the turbulence flame speed \(U_{t} \sim {u}^{\prime }\) in the case of strong turbulence (\({u}^{\prime }>>S_{L} \)) and numerous experiments that show a strong dependence of U_t on the speed of the instantaneous flame S_L. We name this discrepancy the Damköhler-Shelkin paradox. The first aim of the research is to validate and clarify this estimate, which is based on intuitive considerations, as the paradox must be a statement that seems contradictory to observations but is actually true. We analysed the turbulent flame in the context of the original hyperbolic combustion equation that directly describes the leading edge of the flame, which is a locus of the Zel’dovich “leading points” controlling the speed of the turbulent flame. Analysis of the corresponding characteristic equations results in the expression for speed on the steady-state turbulent flame \(U_{t} ={u}^{\prime }\sqrt {1+(S_{L} /{u}^{\prime })^{2}} \), which is the case when \({u}^{\prime }>>S_{L} \) becomes \(U_{t} \cong {u}^{\prime }\). This result confirms and improves the Damköhler-Shelkin estimate \(U_{t} \sim {u}^{\prime }\). The second aim is to resolve the Damköhler-Shelkin paradox. We explain the discrepancy with observations by the fact that turbulent flames are transient due to insufficient residence time in the real burners to reach statistical equilibrium of wrinkle structures of the random flame surface. We consider the transient flame in the intermediate asymptotic stage when the small-scales wrinkles are in statistical equilibrium, while at the same time the large-scale wrinkles are far from equilibrium. The expressions for the flame speed and width, which we deduce using the dimensional analysis and general properties of the ransom surface, \(U_{t} \sim ({u}^{\prime }S_{L})^{1/2}\) and \(\delta _{t} \sim ({u}^{\prime }Lt)^{1/2}\), show that this transient flame is in fact a turbulent mixing layer travelling with constant speed U_t depending on S_L, the intermediate steady propagation (ISP) flame. Qualitative estimations of the times required for the small-scale and large-scale wrinkles to reach statistical equilibrium show that the turbulent Bunsen- and V-flames correspond to the intermediated asymptotic stage, and the turbulent flames with a complete equilibrium structure of the wrinkled flamelet surface are not attainable under laboratory conditions. We present the results of numerical simulations of the impingent flames, which count in favour of the belief that these flames are also transient.     

Assessment of an Evaporation Model in CFD Simulations of a Free Liquid Pool Fire Using the Mass Transfer Number Approach

In the present paper an evaporation model is implemented and assessed in a Computational Fluid Dynamics (CFD) code named ISIS. First, the influence of the cell size and time step on the temperature field is studied via simulations with a prescribed fuel Mass Loss Rate (MLR). Then, the evaporation model is assessed using predictive simulations. The experimental scenario is a 30 cm-diameter heptane pool fire. The average fuel Mass Loss Rate Per Unit Area (MLRPUA) is predicted within 5.5% deviation from the experimental value. In addition, an analysis of the temperature and heat fluxes at the surface of the liquid, the mass transfer coefficient and the temperature inside the liquid is performed.     

Experimental Control of Turbulent Boundary Layers with In-plane Travelling Waves

The experimental control of turbulent boundary layers using streamwise travelling waves of spanwise wall velocity, produced using a novel active surface, is outlined in this paper. The innovative surface comprises a pneumatically actuated compliant structure based on the kagome lattice geometry, supporting a pre-tensioned membrane skin. Careful design of the structure enables waves of variable length and speed to be produced in the flat surface in a robust and repeatable way, at frequencies and amplitudes known to have a favourable influence on the boundary layer. Two surfaces were developed, a preliminary module extending 152 mm in the streamwise direction, and a longer one with a fetch of 2.9 m so that the boundary layer can adjust to the new surface condition imposed by the forcing. With a shorter, 1.5 m portion of the surface actuated, generating an upstream-travelling wave, a drag reduction of 21.5% was recorded in the boundary layer with Re _τ =?1125. At the same flow conditions, a downstream-travelling produced a much smaller drag reduction of 2.6%, agreeing with the observed trends in current simulations. The drag reduction was determined with constant temperature hot-wire measurements of the mean velocity gradient in the viscous sublayer, while simultaneous laser Doppler vibrometer measurements of the surface recorded the wall motion. Despite the mechanics of the dynamic surface resulting in some out-of-plane motion (which is small in comparison to the in-plane streamwise movement), the positive drag reduction results are encouraging for future investigations at higher Reynolds numbers.     

Effect of Multilateral Jet Mixing on Stability and Structure of Turbulent Partially-Premixed Flames

An experimental study was conducted to investigate the effects of multilateral jet mixing, using both three and four side-jets, on the structure and stability of turbulent partially-premixed flames. Particle Image Velocimetry and OH*-chemiluminescence were used to study the effects of geometry and operating conditions on the resulting flow-field and reaction zone structures, respectively. These effects were compared under varying ratios of side-jet to primary flow momentum, whilst keeping the bulk flow constant. It was found that the mixing regimes upstream of the nozzle exit affect the flame characteristics, i.e. an impinging regime is likely to generate a lifted flame whilst a backflow regime is likely to generate an attached flame. Unlike the 4 side-jets cases, the OH* images and v _{r m s} profiles for the 3 side-jets cases show distinct asymmetry, with intense OH* and low velocity fluctuations on the opposite sides of the fuel injection. It was also found that the flow and scalar fields become independent of the upstream conditions, for both 3 and 4 side-jets, after one diameter downstream of the nozzle exit.     

Miscible Thermo-Viscous Fingering Instability in Porous Media. Part 1: Linear Stability Analysis

The development of the thermo-viscous fingering instability of miscible displacements in homogeneous porous media is examined. In this first part of the study dealing with stability analysis, the basic equations and the parameters governing the problem in a rectilinear geometry are developed. An exponential dependence of viscosity on temperature and concentration is represented by two parameters, thermal mobility ratio β _T and a solutal mobility ratio β _C , respectively. Other parameters involved are the Lewis number Le and a thermal-lag coefficient λ. The governing equations are linearized and solved to obtain instability characteristics using either a quasi-steady-state approximation (QSSA) or initial value calculations (IVC). Exact analytical solutions are also obtained for very weakly diffusing systems. Using the QSSA approach, it was found that an increase in thermal mobility ratio β _T is seen to enhance the instability for fixed β _C , Le and λ. For fixed β _C and β _T , a decrease in the thermal-lag coefficient and/or an increase in the Lewis number always decrease the instability. Moreover, strong thermal diffusion at large Le as well as enhanced redistribution of heat between the solid and fluid phases at small λ is seen to alleviate the destabilizing effects of positive β _T . Consequently, the instability gets strictly dominated by the solutal front. The linear stability analysis using IVC approach leads to conclusions similar to the QSSA approach except for the case of large Le and unity λ flow where the instability is seen to get even less pronounced than in the case of a reference isothermal flow of the same β _C , but β _T  = 0. At practically, small value of λ, however, the instability ultimately approaches that due to β _C only.     

Effects of Fuel Lewis Number on Localised Forced Ignition of Turbulent Mixing Layers

The influences of fuel Lewis number Le _F on localised forced ignition of inhomogeneous mixtures are analysed using three-dimensional compressible Direct Numerical Simulations (DNS) of turbulent mixing layers for Le _F  = 0.8, 1.0 and 1.2 and a range of different root-mean-square turbulent velocity fluctuation u′ values. For all Le _F cases a tribrachial flame has been observed in case of successful ignition. However, the lean premixed branch tends to merge with the diffusion flame on the stoichiometric mixture fraction isosurface at later stages of the flame evolution. It has been observed that the maximum values of temperature and reaction rate increase with decreasing Le _F during the period of external energy addition. Moreover, Le _F is found to have a significant effect on the behaviours of mean temperature and fuel reaction rate magnitude conditional on mixture fraction values. It is also found that reaction rate and mixture fraction gradient magnitude \(\vert \nabla \xi \vert \) are negatively correlated at the most reactive region for all values of Le _F explored. The probability of finding high values of \(\vert \nabla \xi \vert \) increases with increasing Le _F . For a given value of u′, the extent of burning decreases with increasing Le _F . A moderate increase in u′ gives rise to an increase in the extent of burning for Le _F  = 0.8 and 1.0, which starts to decrease with further increases in u′. For Le _F  = 1.2, the extent of burning decreases monotonically with increasing u′. The extent of edge flame propagation on the stoichiometric mixture fraction ξ = ξ _st isosurface is characterised by the probability of finding burned gas on this isosurface, which decreases with increasing u′ and Le _F . It has been found that it is easier to obtain self-sustained combustion following localised forced ignition in case of inhomogeneous mixtures than that in the case of homogeneous mixtures with the same energy input, energy deposition duration when the ignition centre is placed at the stoichiometric mixture. The difficultly to sustain combustion unaided by external energy addition in homogeneous mixture is particularly prevalent in the case of Le _F  = 1.2.     

A RANS simulation toward the effect of turbulence and cavitation on spray propagation and combustion characteristics

A multidimensional computational fluid dynamic code was developed and integrated with probability density function combustion model to give the detailed account of multiphase fluid flow. The vapor phase within injector domain is treated with Reynolds-averaged Navier–Stokes technique. A new parameter is proposed which is an index of plane-cut spray propagation and takes into account two parameters of spray penetration length and cone angle at the same time. It was found that spray propagation factor (SPI) tends to increase at lower r/d ratios, although the spray penetration tends to decrease. The results of SPI obtained by empirical correlation of Hay and Jones were compared with the simulation computation as a function of respective r/d ratio. Based on the results of this study, the spray distribution on plane area has proportional correlation with heat release amount, NO _x emission mass fraction, and soot concentration reduction. Higher cavitation is attributed to the sharp edge of nozzle entrance, yielding better liquid jet disintegration and smaller spray droplet that reduces soot mass fraction of late combustion process. In order to have better insight of cavitation phenomenon, turbulence magnitude in nozzle and combustion chamber was acquired and depicted along with spray velocity.     

Spray-Flame Dynamics in a Rich Droplet Array

We numerically study spray-flame dynamics. The initial state of the spray is schematized by alkane droplets located at the nodes of a face-centered 2D-lattice. The droplets are surrounded by a gaseous mixture of alkane and air. The lattice spacing s reduced by the combustion length scale is large enough to consider that the chemical reaction occurs in a heterogeneous medium. The overall spray equivalence ratio is denoted by ?_T, with ?_T = ?_L + ?_G, where ?_G corresponds to the equivalence ratio of the gaseous surrounding mixture at the initial saturated partial pressure, while ?_L is the so-called liquid loading. To model such a heterogenous combustion, the retained chemical scheme is a global irreversible one-step reaction governed by an Arrhenius law, with a modified heat of reaction depending on the local equivalence ratio. ?_T is chosen in the range 0.9 ≤ ?_T ≤ 2. Three geometries (s = 3, s = 6, s = 12) and four liquid loadings, ?_L = 0.3, ?_L = 0.5, ?_L = 0.7, ?_L = 0.85 are studied. In the rich sprays, our model qualitatively retrieves the recent experimental measurements: the rich spray-flames can propagate faster than the single-phase flames with the same overall equivalence ratio. To analyse the conditions for this enhancement, we introduce the concept of “spray Peclet number”, which compares the droplet vaporization time with the combustion propagation time of the single-phase flame spreading in the fresh surrounding mixture.     

Linearized Stability for Semiflows Generated by a Class of Neutral Equations,with Applications to State-Dependent Delays

We prove a principle of linearized stability for semiflows generated by neutral functional differential equations of the form x′(t) = g(? x _t , x _t ). The state space is a closed subset in a manifold of C ²-functions. Applications include equations with state-dependent delay, as for example x′(t) = a x′(t + d(x(t))) + f (x(t + r(x(t)))) with \({a\in\mathbb{R}, d:\mathbb{R}\to(-h,0), f:\mathbb{R}\to\mathbb{R}, r:\mathbb{R}\to[-h,0]}\).