An apparatus-independent extinction strain rate in counterflow flames

Resistance to extinction by stretch is a key property of any flame, and recent work has shown that this property controls the overall structure of several important types of turbulent flames. Multiple definitions of the critical strain rate at extinction (ESR) have been presented in the literature. However, even if the same definition is used, different experiments report different extinction strain rates for flames burning the same fuel-air mixture at very similar temperatures using similarly constructed opposed-flow instruments. Here we show that at extinction, all these flames are essentially identical, so one would expect that each would be assigned the same value of a parameter representing its intrinsic resistance-to-stretch-induced-extinction, regardless of the specifics of the experimental apparatus. A similar situation arises in laminar flame speed measurements since different apparatuses could result in different strain rate distributions. In that instance, the community has agreed to report the unstretched laminar flame speed, and methods have been developed to translate the experimental (stretched) flame speed into a universal unstretched laminar flame speed. We propose an analogous method for translating experimental measurements for stretch-induced extinction into an unambiguous and apparatus-independent quantity (ESR_∞) by extrapolating to infinite opposing burner separation distance. The uniqueness of the flame at extinction is shown numerically and supported experimentally for twin premixed, single premixed, and diffusion flames at Lewis numbers greater than and less than one. A method for deriving ESR_∞ from finite-boundary experimental studies is proposed and demonstrated for methane and propane experimental diffusion and premixed single flame data. The two values agree within the range of ESR differences typically observed between experimental measurements and simulation results for the traditional ESR definition.     

Effect of Soret diffusion on lean hydrogen/air flames at normal and elevated pressure and temperature

The influence of Soret diffusion on lean premixed flames propagating in hydrogen/air mixtures is numerically investigated with a detailed chemical and transport models at normal and elevated pressure and temperature. The Soret diffusion influence on the one-dimensional (1D) flame mass burning rate and two-dimensional (2D) flame propagating characteristics is analysed, revealing a strong dependency on flame stretch rate, pressure and temperature. For 1D flames, at normal pressure and temperature, with an increase of Karlovitz number from 0 to 0.4, the mass burning rate is first reduced and then enhanced by Soret diffusion of H₂ while it is reduced by Soret diffusion of H. The influence of Soret diffusion of H₂ is enhanced by pressure and reduced by temperature. On the contrary, the influence of Soret diffusion of H is reduced by pressure and enhanced by temperature. For 2D flames, at normal pressure and temperature, during the early phase of flame evolution, flames with Soret diffusion display more curved flame cells. Pressure enhances this effect, while temperature reduces it. The influence of Soret diffusion of H₂ on the global consumption speed is enhanced at elevated pressure. The influence of Soret diffusion of H on the global consumption speed is enhanced at elevated temperature. The flame evolution is more affected by Soret diffusion in the early phase of propagation than in the long run due to the local enrichment of H₂ caused by flame curvature effects. The present study provides new insights into the Soret diffusion effect on the characteristics of lean hydrogen/air flames at conditions that are relevant to practical applications, e.g. gas engines and turbines.     

On disintegration of near-limit cellular flames

A strongly non-linear geometrically-invariant model for the dynamics of near-limit cellular flame is proposed, where the flame evolution is governed by a system of equations for the flame interface and its temperature. The model generalizes its earlier weakly non-linear version pertinent to a mildly perturbed planar flame. Numerical simulations of the new model show that at sufficiently high levels of heat losses the cellular flame resulting from the diffusive instability exhibits a tendency toward self-fragmentation, quite in line with direct numerical simulations of the associated reaction-diffusion system.     

Uncertainty assessment of species measurements in acetone counterflow diffusion flames

To quantitatively understand the uncertainty of intrusive species sampling measurements using a microprobe, velocity and speciation profiles of acetone counterflow diffusion flames have been experimentally investigated with cross validations using non-intrusive particle image velocimetry (PIV) and laser induced fluorescence (LIF) measurements. It is shown that the separation distance between the fuel and oxidizer nozzles needs to be sufficiently large to achieve uniform radial velocity profiles at the nozzle exit and accurate measurements of fuel concentration distributions in flames. The impacts of the diffusion flame location relative to the stagnation plane and the diffusion flame thickness on quantitative species sampling are investigated by varying the fuel to oxygen ratio as well as nitrogen and helium as fuel diluents. The results show that the diffusion flame needs to be located on the fuel side far from the stagnation plane in order to obtain reliable speciation measurements of fuel oxidation-related species. For helium dilution in the fuel side, a significant deviation from the model prediction is found due to the excessively fast diffusion velocity of helium. The impact of the intrusive probe on the flow field and the structure of the counterflow diffusion flame are identified by acetone and OH LIF measurements. The uncertainty in the speciation measurement associated with flow perturbations by the probe is quantified and found to be comparable to the outer diameter of the probe, ±0.3 mm. A simple Reynolds number analysis shows that the flow near the probe is just on the outskirts of the Stokes regime. Finally, the structure of the acetone diffusion flame is measured quantitatively with species measurements of ethane, ethylene, and acetylene. The comparison between predictions and measurements indicate that the current C₂ kinetic mechanism needs to be improved for quantitative prediction of the acetone flame structures.     

Effects of gas and soot radiation on soot formation in counterflow ethylene diffusion flames

Numerical study of soot formation in counterflow ethylene diffusion flames at atmospheric pressure was conducted using detailed chemistry and complex thermal and transport properties. Soot kinetics was modelled using a semi-empirical two-equation model. Radiation heat transfer was calculated using the discrete-ordinates method coupled with an accurate band model. The calculated soot volume fractions are in reasonably good agreement with the experimental results in the literature. The individual effects of gas and soot radiation on soot formation were also investigated.     

Pressure and temperature dependence of soot in highly controlled counterflow ethylene diffusion flames

Soot volume fraction and dispersion index were measured by pyrometry in a series of highly controlled counterflow diffusion flames, with peak temperatures, T_max, spanning a few hundred degrees and pressure covering the 0.1–0.8 MPa range. An unprecedented level of control was implemented by selecting flames with a self-similar structure to ensure that the normalized temperature-time history experienced by the reactants was the same, regardless of pressure. The self-similarity was verified by suitably rescaling the transverse coordinate with respect to a characteristic diffusion length. At constant T_max, the soot volume fraction increases approximately by two orders of magnitude as the pressure is raised from 1 atm to 4 atm, and by one to two additional orders of magnitude with an additional doubling of the pressure to 8 atm. At constant pressure, the soot load spans two to three orders of magnitude and soot formation exhibits increased sensitivity to temperature as the pressure is raised. Soot inception occurs near the flame, with an increase in soot concentration that becomes steeper at higher T_max. The increase is accompanied by a decrease in the dispersion exponent that is suggestive of dehydrogenation and aging of the particles and is sharper at higher T_max. Soot experiences continuous growth in a monotonically decreasing temperature field until it is convected away radially at the stagnation plane, with essentially no opportunity for oxidation. Evidence of two distinct mechanisms for soot formation was found: the classic high temperature, high activation energy process affecting soot formed in the vicinity of the flame and followed by dehydrogenation; and a relatively low-temperature, zero activation energy process, associated with the increase in volume fraction at low-temperatures in proximity of the stagnation plane. The latter is tentatively attributed to dimerization of aromatics, as revealed by the concurrent increase in the dispersion index corresponding to an increase in the particle hydrogen content.     

Fuel-dilution effect on differential molecular diffusion in laminar hydrogen diffusion flames

Laminar flame calculations have been made for a Tsuji counterflow geometry to investigate salient features caused by the differential diffusion effect in nitrogen-diluted hydrogen diffusion flames. A strong dependence of the differential diffusion parameter z_H on fuel dilution is found, where z_H is the difference of the mixture fractions based on H and O elements. The strain rate, however, appears to have a relatively minor impact on z_H. A simplified transport equation for the z_H parameter has been derived to explain qualitatively the behaviours exhibited in the numerical solutions. Two source terms of z_H are identified in the transport equation; one is due to mixing among species of different diffusion coefficients and the other one is associated with chemical reactions of H₂. More importantly, the second source term is found to be dominant in reacting flows, and it increases with inert gas dilution. This feature causes the differential diffusion parameter to increase with the amount of fuel dilution. The z_H values at the stoichiometric position are shown to correlate well with the ratio, Y_H2O|_max/(Z_H,1?Z_H,2), which may be useful for quantifying the influence of chemical reactions on the differential diffusion effect. For flames at low strain rates, the scalar dissipation rate exhibits a local minimum near the stoichiometric position. This peculiar feature is found to be caused by the differential diffusion effect modulated by chemical reactions. The local minimum in the scalar dissipation rate disappears at high strain rates when the convective transport overwhelms the molecular diffusion.     

Soot and NO formation in counterflow ethylene/oxygen/nitrogen diffusion flames

Formation of soot and NO in counterflow ethylene/oxygen/nitrogen diffusion flames was numerically investigated. Detailed chemistry and complex thermal and transport properties were used. A simplified two-equation soot model was adopted. The results indicate that NO emission has negligible influence on soot formation. However, soot formation affects the emission of NO through the radiation induced thermal effect and the reaction induced chemical effect. When the oxygen index of the oxidant stream is lower, the relative influence of chemical reaction caused by soot on NO emission is more important, while the relative influence of the radiation induced thermal effect becomes more important for the flame with a higher oxygen index in the oxidant stream.     

Carbon nanotube synthesis on catalytic metal alloys in methane/air counterflow diffusion flames

Various morphologies of multi-walled carbon nanotubes (MWNTs) are grown catalytically on metal-alloy probes in counterflow diffusion flames using methane as fuel. Carbon nanotube (CNT) properties and morphologies are investigated as functions of local gas-phase temperatures, C-related species concentrations (e.g. C₂H₂, CO), sampling positions, C₂H₂ adding to the fuel, and metal-alloy compositions (i.e., Fe, Fe/Cr, Ni/Cu, Ni/Ti, Ni/Cr, Ni/Cr/Fe). MWNTs grow optimally in non-sooty regions of the flames. C₂H₂ addition is found to promote direct synthesis of vertically well-aligned MWNTs with uniform diameters from Ni/Cr/Fe and Ni/Ti alloys.     

The influence of water on extinction and ignition of hydrogen and methane flames

The influence of water vapor on critical conditions of extinction and autoignition of premixed and nonpremixed flames is investigated. The fuels tested are hydrogen (H₂) and methane (CH₄). Studies on premixed systems are carried out by injecting a premixed reactant stream made up of fuel, oxygen (O₂), and nitrogen (N₂) from one duct, and an inert-gas stream of N₂ from the other duct. Critical conditions of extinction are measured for various amounts of water vapor added to the premixed reactant stream. The ratio of fuel to oxygen is maintained at a constant value, and the amounts of water vapor and nitrogen are so chosen that the adiabatic temperature remains the same. This ensures that the physical influence of water is the same for all cases. Therefore, changes in values for the critical conditions of extinction are attributed to the chemical influence of water vapor. Studies on nonpremixed systems are carried out by injecting a fuel stream made up of fuel and N₂ from one duct ,and an oxidizer stream made up of O₂ and N₂ from the other duct. Critical conditions of extinction are measured with water vapor added to the oxidizer stream. The concentrations of reactants are so chosen that the adiabatic temperature and the flame position stay the same for all cases. Critical conditions of autoignition are measured by preheating the oxidizer stream of the nonpremixed system. Water vapor is added to the oxidizer stream. Numerical calculations are performed using a detailed chemical-kinetic mechanism and compared with measurements. Experimental and numerical studies show that addition of water makes the premixed and nonpremixed flames easier to extinguish and harder to ignite. The chemical influence of water is attributed to its enhanced chaperon efficiency in three body reactions.     

Investigation of droplet combustion in strained counterflow diffusion flames using planar laser-induced fluorescence

Experimental investigation of an isolated droplet burning in a convective flow is reported. Acetone droplets were injected in a steady laminar diffusion counterflow flame operating with methane. Planar laser-induced fluorescence measurements applied to OH radical and acetone was used to measure the spatial distribution of fuel vapour and the structure of the flame front around the droplet. High-magnification optics was used in order to image flow areas with a ratio of 1:1.2. The different combustion regimes of an isolated droplet could be observed from the configuration of the envelope flame to that of the boundary-layer flame, and occurrence of these regimes was found to depend on the droplet Reynolds number. Experimental results were compared with 1D numerical simulations using detailed chemistry for the configuration of the envelope flame. Good agreement was obtained for the radial profile of both OH radical and fuel vapour. Influence of droplet dynamics on the counterflow flame front was also investigated. Results show that the flame front could be strongly distorted by the droplet crossing. In particular, droplets with high velocity led to local extinction of the flame front whereas droplets with low velocity could ignite within the flame front and burn on the oxidiser side. PACS 33.50.-j; 42.62.-b; 47.55.D-; 47.70.Pq; 47.80.Jk     

Synthesis of carbon nanotubes on Ni-alloy and Si-substrates using counterflow methane–air diffusion flames

We have conducted an experimental study to investigate the synthesis of multi-walled carbon nanotubes (CNTs) in counterflow methane–air diffusion flames, with emphasis on effects of catalyst, temperature, and the air-side strain rate of the flow on CNTs growth. The counterflow flame was formed by fuel (CH₄ or CH₄ + N₂) and air streams impinging on each other. Two types of substrates were used to deposit CNTs. Ni-alloy (60% Ni + 26% Cr + 14% Fe) wire substrates synthesized curved and entangled CNTs, which have both straight and bamboo-like structures; Si-substrates with porous anodic aluminum oxide (AAO) nanotemplates synthesized well-aligned, self-assembled CNTs. These CNTs grown inside nanopores had a uniform geometry with controllable length and diameter. The axial temperature profiles of the flow were measured by a 125 μm diameter Pt/10% Rh–Pt thermocouple with a 0.3 mm bead junction. It was found that temperature could affect not only the success of CNTs synthesis, but also the morphology of synthesized CNTs. It was also found, against previous general belief, that there was a common temperature region (1023–1073 K) in chemical vapor deposition (CVD) and counterflow diffusion flames where CNTs could be produced. CNTs synthesized in counterflow flames were significantly affected by air-side strain rate not through the residence time, but through carbon sources available for CNTs growth. Off-symmetric counterflow flames could synthesize high-quality CNTs because with this configuration carbon sources at the fuel side could easily diffuse across the stagnation surface to support CNTs growth. These results show the feasibility of using counterflow flames to synthesize CNTs for particular applications such as fabricating nanoscale electronic devices.     

Experimental and numerical investigation of microscale hydrogen diffusion flames

Characteristics of microscale hydrogen diffusion flames produced from sub-millimeter diameter (d = 0.2 and 0.48 mm) tubes are investigated using non-intrusive UV Raman scattering coupled with LIPF technique. Simultaneous, temporally and spatially resolved point measurements of temperature, major species concentrations (O₂, N₂, H₂O, and H₂), and absolute hydroxyl radical concentration (OH) are made in the microflames for the first time. The probe volume is 0.02 × 0.04 × 0.04 mm³. In addition, photographs and 2-D OH imaging techniques are employed to illustrate the flame shapes and reaction zones. Several important features are identified from the detailed measurements of microflames. Qualitative 2-D OH imaging indicates that a spherical flame is formed with a radius of about 1 mm as the tube diameter is reduced to 0.2 mm. Raman/LIPF measurements show that the coupled effect of ambient air leakage and pre-heating enhanced thermal diffusion of H₂ leads to lean-burn conditions for the flame. The calculated characteristic features and properties indicate that the buoyancy effect is minor while the flames are in the convection–diffusion controlled regime because of low Peclet number. Also, the effect of Peclet number on the flame shape is minor as the flame is in the convection–diffusion controlled regime. Comparisons between the predicted and measured data indicate that the trends of temperature, major species, and OH distributions are properly modeled. However, the code does not properly predict the air entrainment and pre-heating enhanced thermal-diffusive effects. Therefore, thermal diffusion for light species and different combustion models might need to be considered in the simulation of microflame structure.     

Chemical structure of atmospheric pressure premixed laminar formic acid/hydrogen flames

The work presents an experimental and kinetic modeling study of laminar premixed formic acid [HC(O)OH]/H₂/O₂/Ar flames at different equivalence ratios (φ=0.85, 1.1 and 1.3) stabilized on a flat burner at atmospheric pressure, as well as laminar flame speed of HC(O)OH/O₂/Ar flames (φ=0.5–1.5) at 1 atm. Flame structure as well as laminar flame speed were simulated using three different detailed chemical kinetic mechanisms proposed for formic acid oxidation. The components in the fuel blends show different consumption profiles, namely, hydrogen is consumed slower than formic acid. According to kinetic analysis, the reason of the observed phenomenon is that the studied flames have hydrogen as a fuel but also as an intermediate product formed from HC(O)OH decomposition. Comparison of the measured and simulated flame structure shows that all the mechanisms satisfactorily predict the mole fraction profiles of the reactants, main products, and intermediates. It is noteworthy that the mechanisms proposed by Glarborg et al., Konnov et al. and the updated AramcoMech2.0 adequately predict the spatial variations in the mole fractions of free radicals, such as H, OH O and HO₂. However, some drawbacks of the mechanisms used were identified; in particular, they predict different concentrations of CH₂O. As for laminar flame speed simulations, the Konnov et al. mechanism predicts around two times higher values than in experiment, while the Glarborg et al. and updated AramcoMech2.0 show good agreement with the experimental data.     

Modeled quenching limits of spherical hydrogen diffusion flames

Hydrogen–air diffusion flames were modeled with an emphasis on kinetic extinction. The flames were one-dimensional spherical laminar diffusion flames supported by adiabatic porous burners of various diameters. Behavior of normal (H₂ flowing into quiescent air) and inverse (air flowing into quiescent H₂) configurations were considered using detailed H₂/O₂ chemistry and transport properties with updated light component diffusivities. For the same heat release rate, inverse flames were found to be smaller and 290 K hotter than normal flames. The weakest normal flame that could be achieved before quenching has an overall heat release rate of 0.25 W, compared to 1.4 W for the weakest inverse flame. There is extensive leakage of the ambient reactant for both normal and inverse flames near extinction, which results in a premixed flame regime for diffusion flames except for the smallest burners with radii on the order of 1 μm. At high flow rates H + OH(+M) → H₂O(+M) contributes nearly 50% of the net heat release. However at flow rates approaching quenching limits, H + O₂(+M) → HO₂(+M) is the elementary reaction with the largest heat release rate.     

Effect of hydrogen enrichment of laminar ethylene diffusion flames on thermal structure and soot yields at pressures up to 10 bar

A series of high-pressure experiments were conducted to assess the influence of hydrogen enrichment of laminar diffusion flames of nitrogen-diluted ethylene on the thermal flame structure and soot yields at pressures above atmospheric. In parallel experiments, added hydrogen is replaced by helium, either in equal mole fractions or in mass fractions, to evaluate the thermal, dilution, and direct chemical interaction effects of hydrogen in soot formation. Experiments covered pressures from atmospheric to 10 bar. In the first set of experiments, conducted at 3, 6, and 10 bar pressure, base fuel was an ethylene-nitrogen mixture with 33.3% ethylene and 66.7% nitrogen (by mole as well as by mass). This base fuel was doped with either hydrogen or helium such that hydrogen and helium mass fractions and mole fractions in the fuel stream are matched in two cases. In the second set of experiments, which were conducted at 1.2 bar pressure with ethylene as the base fuel, hydrogen or helium is added such that additive mole fraction in the fuel stream was 44%. Temperature measurements in the first set of experiments indicate that, when hydrogen is added to nitrogen-diluted ethylene, the changes in the temperature field of the co-flow diffusion flames are negligible, except at lower in the flame where hydrogen added flames display slightly higher temperatures. When helium is added instead of hydrogen, however, the temperatures were measurably lower than those of the base fuel. Results show that, once the dilution effects are accounted for, the hydrogen addition to ethylene does not suppress soot formation by direct chemical interaction at elevated pressures. These findings, which are not in agreement with the previous experimental results obtained at atmospheric pressure, are discussed in terms of the higher molecular diffusivity of hydrogen and shorter residence times of high-pressure flames.     

Temperature measurements in sooting counterflow diffusion flames using laser-induced fluorescence of flame-produced nitric oxide

Laser-based diagnostic methods are often used for non-intrusive studies of delicate processes of soot formation. When soot particles are heated by the laser pulse, their size distribution can be estimated from the cooling rate, provided that the local gas temperature is known. However, strong light absorption, scattering and fluorescence in sooting environment hinder non-intrusive laser-based temperature measurements. Methods based on fitting of laser-induced fluorescence spectra work well in stationary flames but usually require temperature tracer seeded into the flame. We have shown that in counterflow diffusion flames, often used for soot-formation studies, enough nitric oxide is produced for two-dimensional temperature imaging. Measured temperature profiles agree very well with chemical kinetic calculations for a variety of fuels if laser intensity is reduced to keep NO excitation in the linear regime. Gas composition affects line shapes at temperatures below 600 K and should be taken into account for accurate measurements.     

Dynamics and burning limits of near-limit hot,warm, and cool diffusion flames of dimethyl ether at elevated pressures

The near-limit diffusion flame regimes and extinction limits of dimethyl ether at elevated pressures and temperatures are examined numerically in the counterflow geometry with and without radiation at different oxygen concentrations. It is found that there are three different flame regimes—hot flame, warm flame, and cool flame—which exist, respectively, at high, intermediate, and low temperatures. Furthermore, they are governed by three distinct chain-branching reaction pathways. The results demonstrate that the warm flame has a double reaction zone structure and plays a critical role in the transition between cool and hot flames. It is also shown that the cool flame can be formed in several different ways: by either radiative extinction or stretch extinction of a hot flame or by stretch extinction of a warm flame. A warm flame can also be formed by radiative extinction of a hot flame or ignition of a cool flame. A general €-shaped flammability diagram showing the burning limits of all three flame regimes at different oxygen mole fractions is obtained. The results show that thermal radiation, reactant concentration, temperature, and pressure all have significant impacts on the flammable regions of the three flame regimes. Increases in oxidizer temperature, oxygen concentration, and pressure shift the cool flame regime to higher stretch rates and cause the warm flame to have two extinction limits. At elevated temperatures, it is found that there is a direct transition between the hot flame and warm flame at low stretch rates. The results also show that, unlike the hot flame, the cool flame structure cannot be scaled by using pressure-weighted stretch rates due to the its significant reactant leakage and strong dependence of reactivity on pressure. The present results advance the understanding of near-limit flame dynamics and provide guidance for experimental observation of different flame regimes.     

The two-dimensional structure of low strain rate counterflow nonpremixed-methane flames in normal and microgravity

The structure and extinction of low strain rate nonpremixed methane–air flames was studied numerically and experimentally. A time-dependent axisymmetric two-dimensional (2D) model considering buoyancy effects and radiative heat transfer was developed to capture the structure and extinction limits of normal gravity (1-g) and zero gravity (0-g) flames. For comparison with the 2D modelling results, a one-dimensional (1D) flamelet computation using a previously developed numerical code was exercised to provide information on the 0-g flames. A 3-step global reaction mechanism was used in both the 1D and 2D computations to predict the measured extinction limit and flame temperature. Photographic images of flames undergoing the process of extinction were compared with model calculations. The axisymmetric numerical model was validated by comparing flame shapes, temperature profiles, and extinction limits with experiments and with the 1D computational results. The 2D computations yielded insight into the extinction mode and flame structure. A specific maximum heat release rate was introduced to quantify the local flame strength and to elucidate the extinction mechanism. The contribution by each term in the energy equation to the heat release rate was evaluated to investigate the multi-dimensional structure and radiative extinction of the 1-g flames. Two combustion regimes depending on the extinction mode were identified. Lateral heat loss effects and multi-dimensional flame and flow structure were also found. At low strain rates in 1-g flames (‘regime A’), the flame is extinguished from the weak outer edge of the flame, which is attributed to a multi-dimensional flame structure and flow field. At high strain rates, (‘regime B’), the flame extinction initiates near the flame centreline owing to an increased diluent concentration in the reaction zone, similar to the extinction mode of 1D flames. These two extinction modes can be clearly explained by consideration of the specific maximum heat release rate.     

Analysis of steady state polydisperse counterflow spray diffusion flames in the large Stokes number limit

A new mathematical analysis is presented of certain aspects of the behavior of opposed flow polydisperse spray diffusion flames within the framework of a model in which large slip is permitted between the droplets and their host surroundings. The sectional approach is used to model the polydisperse spray. Operating conditions are identified under which the inverses of sectional Stokes numbers are small spray-related parameters to be used in a perturbation analysis of the liquid phase governing equations. The steady state equations and their solutions are similar in form to the equivalent equations considered in previous work of the authors in which dynamical equilibrium of the droplets with the carrier phase was assumed. However, here there is much more mathematical complexity involved in the spray equations solution. A hybrid Eulerian–Lagrangian approach is also suggested to get an insight into the phenomenon of reversal in the motion of the droplets that has been reported in independent experimental and computational research. Computed results based on the analytical solutions up to the 1st order of approximation reveal the influence of large droplet slip on the droplets velocity field and on the spray diffusion flame’s thermal field, for which appreciable heterogeneous combustion can occur under the operating conditions considered.