Experimental study of the effect of helium/nitrogen concentration and initial droplet diameter on nonane droplet combustion with minimal convection

In this paper, we study the influence of inert concentration and initial droplet diameter on nonane (C₉H₂₀) droplet combustion in an environment that promotes spherical droplet flames. The oxygen concentration is fixed while the inert is varied between nitrogen and helium. A range of initial droplet diameters (D_o) are examined in each ambient gas: 0.4 mm < D_o < 0.8 mm; and an oxidizing ambiance consisting of 30% oxygen (fixed) and 70% inert (fixed), with the inert in turn composed of mixtures of nitrogen and helium in concentrations of 0, 25, 50, 75, and 100% N₂. The experiments are carried out at normal atmospheric pressure in a cold ambiance (room temperature) under low gravity to minimize the influence of convection and promote spherical droplet flames. For burning within a helium inert (0% N₂), the droplet flames are entirely blue and there is no influence of initial droplet diameter on the local burning rate (K). With increasing dilution by nitrogen, droplet flames show significant yellow luminosity indicating the presence of soot and the individual burning histories show K reducing with increasing D_o. The evolution of droplet diameter D(t) is nonlinear for a given D_o in the presence of either helium or nitrogen inerts indicating that soot formation has little to do with nonlinear burning. A correlation is presented of the data in the form where the effective burning rate, K′, and ε are concentration-dependent. Correlations for these parameters are presented in the paper.     

Nonane droplet combustion with and without buoyant convection: Flame structure, burning rate and extinction in air and helium

The burning and extinction characteristics of isolated small nonane droplets are examined in a buoyant convective environment and in an environment with no external axial convection (as created by doing experiments at low gravity) to promote spherical droplet flames. The ambience is air and a mixture of 30%O₂/70%He to assess the influence of soot formation. The initial droplet diameter (D_o) ranges from 0.4 to 0.95 mm. Measurements are reported of the extinction diameter and time to extinction, and of the evolution of droplet diameter, flame diameter, soot shell diameter, burning rate, and broadband radiative emissions.In a buoyancy-free environment for air larger droplets burn slower than smaller droplets for the range of D_o examined, which is attributed to the influence of soot. In the presence of a buoyant flow in air, no influence of D_o is observed on the burning rate while the buoyant flames are still heavily sooting. The effect of D_o is believed to be due to a combination of dominance of the nonluminous, nonsooting, portion of the buoyant flame around the forward half of the droplet on heat transport and the secondary role of the luminous wake portion of the flame. In a non-sooting helium inert at low gravity, no effect of D_o is found on the evolution of droplet diameter.Flame extinction is observed only in the 30%O₂/70%He ambience. For all of the observations, extinction appears to occur before the disappearance of the droplet which is then followed by a period of evaporation. The extinction diameter and time to extinction increases with D_o and an empirical correlation is presented for these two variables.     

Effects of fine fuel droplets on a laminar flame stabilized in a partially prevaporized spray stream

A partially prevaporized spray burner was developed to investigate the interaction between fuel droplets and a flame. Monodispersed partially prevaporized ethanol sprays with narrow diameter distribution were generated by the condensation method using rapid pressure reduction of a saturated ethanol vapor–air mixture. A tilted flat flame was stabilized at the nozzle exit using a hot wire. Particle tracking velocimetry (PTV) was applied to measurements of the droplet velocity; the laminar burning velocity was obtained from gas velocity derived from the droplet velocity. Observations were made of flames in partially prevaporized spray streams with mean droplet diameters of 7 μm and the liquid equivalence ratios of 0.2; the total equivalence ratio was varied. In all cases, a sharp vaporization plane was observed in front of the blue flame. Flame oscillation was observed on the fuel-rich side. At strain rates under 50 s⁻¹, the change in the burning velocity with the strain rate is small in fuel-lean spray streams. In spray streams of 0.7 and 0.8 in the total equivalence ratio, burning velocity increases with strain rates of greater than 50 s⁻¹. However, in spray streams with 0.9 and 1.0 in the total equivalence ratio, burning velocity decreases as the strain rate increases. At strain rates greater than 80 s⁻¹, burning velocity decreases with an increased gas equivalence ratio. The effect of mean droplet diameter, and the entry length of droplets into a flame on the laminar burning velocity, were also investigated to interpret the effect of the strain rate on the laminar burning velocity of partially prevaporized sprays.     

Burning velocity of triple flames with gentle concentration gradient

We investigated the local flame speed of a two-dimensional, methane-air triple flame in a rectangular burner. The velocity fields and the concentration profiles were measured with particle image velocimetry and the Rayleigh scattering method, respectively. There was a requisite combination of initial velocity and initial concentration gradient for consistency of the local concentration gradient at the leading edge of the flame. In these cases, the flame curvatures were also consistent. Accordingly, the burning velocity, defined as local flow velocity at the triple point, was determined by the flame curvature. The burning velocity increased with increasing flame curvature, when the curvature was near zero. After that, the burning velocity decreased with increasing curvature. The peak value thus exceeded the adiabatic one-dimensional laminar burning velocity. Comparing the effects of the measured flame stretch rate on the flow strain κ_s and flame curvature κ_c, κ_s is larger and increases more rapidly than κ_c for flame curvatures satisfying 1/R_f < 250 m⁻¹ and then becomes constant while κ_c still increases for 250 m⁻¹ < 1/R_f, so that κ_c becomes much larger than κ_s. There is also a peak in burning velocity at roughly the transition in flame curvature specified above. Therefore, the burning velocity for a low concentration gradient correlates with the flame stretch rate.     

Spatiotemporal measurements of flame stretch and propagation rates for lean and rich CH4/air premixed flames interacting with a turbulent-wake

A 1.5 m long turbulent-wake combustion vessel with a 0.15 m × 0.15 m cross-sectional area is proposed for spatiotemporal measurements of curvature, strain, dilatation and burning rates along a freely downward-propagating premixed flame interacting with a parallel row of staggered vortex pairs having both compression (negative) and extension (positive) strains simultaneously. The wanted wake is generated by rapidly withdrawing an electrically-controlled, horizontally-oriented sliding plate of 5 mm thickness for flame–wake interactions. Both rich and lean CH₄/air flames at the equivalence ratios  = 1.4 and  = 0.7 with nearly the same laminar burning velocity are studied, where flame–wake interactions and their time-dependent velocity fields are obtained by high-speed, high-resolution DPIV and laser-tomography. Correlations among curvature, strain, stretch, and dilatation rates along wrinkled flame fronts at different times are measured and thus their influences on front propagation rates can be analyzed. It is found that strain-related effects have significant influence on front propagation rates of rich CH₄/air (diffusionally stable) flames even when the curvature weights more in the total stretch than the strain rate does. The local propagation rates along the wrinkled flame front are more intense at negative strain rates corresponding to positive peak dilatation rates but the global propagation rate averaged along the rich flame front remains constant during all period of flame–wake interaction. For lean CH₄/air (diffusionally unstable) flames, the curvature becomes a dominant parameter influencing the structure and propagation of the wrinkled flame front, where both local and global propagation rates increase significantly with time, showing unsteady flame propagation. These experimental results suggest that the theory of laminar flame stretch can be applicable to a more complex flame–wake interaction involving unsteadiness and multitudinous interactions between vortices.     

Droplet combustion in standing sound waves

Interaction between droplet combustion and acoustic oscillation is clarified. As the simplest model, an isolated fuel droplet is combusted in a standing sound wave. Apart from the conventional idea that oscillatory component of flow influences heat and mass transfer and promotes combustion, a new model that a secondary flow dominates combustion promotion is examined. The secondary flow, found by the authors in the previous work, is driven by acoustic radiation force due to Reynolds normal stress, and named as thermo-acoustic streaming. Since the force is described by the same equation as buoyancy, i.e., F = ΔρVg, the nature of the streaming is thought to be the same as natural convection. The flow patterns of the streaming are analyzed and its influence on burning rate of a droplet is predicted. Experimental investigation was mainly done with burning droplets located in the middle of node and anti-node of standing sound waves. This location realizes the strongest streaming. By varying sound pressure level, ambient pressure, and acoustic frequency, the strength of the streaming was controlled. Flame configuration including soot and burning rate were examined. Microgravity conditions were employed to clarify the influence of acoustic field through the streaming, since it is similar to and must be distinguished from natural convection. Experiments using microgravity conditions confirmed the new combustion promotion model and the way to quantify it. By introducing a new non-dimensional number Gr_a, that is the ratio of acoustic radiation force to viscosity, burning rate constants for various ambient and sound conditions are rearranged. As a result, it was found that the excess burning rate (k/k₀ − 1) is proportional to or , for weak sound and for strong sound, respectively.     

Nonlinear kinematic response of premixed flames to harmonic velocity disturbances

This paper analyzes the nonlinear dynamics of premixed flames responding to harmonic velocity disturbances. These nonlinear dynamics were studied by solving a constant flame speed front tracking equation for the flame’s response to harmonically oscillating velocity disturbances. The solution to these equations is used to quantify the transfer function relating the ratio of the normalized flame area to velocity fluctuations, G = (A′/A_o)/(u′/u_o), upon the amplitude of velocity oscillations, ε = u′/u_o. Due to nonlinearities, the amplitude of this transfer function relative to its linear value decreases with increasing amplitude of velocity oscillation, u′/u_o. In contrast, the transfer function phase exhibits almost no amplitude dependence. The velocity amplitude where transfer function nonlinearities become significant depends strongly upon three parameters: a Strouhal number, St = ωL_f/u_o (where L_f is the flame length), the ratio of the flame length to width, β = L_f/R, and the flame shape in the absence of perturbations (i.e., conical, inverted wedge, etc.). In the linear case, the transfer function, G, depends only upon an algebraic combination of the first two parameters, given by St₂ = St (1 + β²)/β². In general, however, G exhibits a distinct dependence upon both parameters St and β. In particular, we show that the nonlinear response of G is an intrinsically dynamic phenomenon; i.e., its quasi-steady response (St 1) is purely linear. As such, nonlinearity is enhanced with increasing Strouhal numbers. In contrast, nonlinearity is suppressed at large β values; as such, the response of a long flame remains quite similar to its linear value, even at large ε values where the flame front exhibits substantial corrugation and cusping. Finally, we show that the response of conical flames remains much more linear at comparable disturbance amplitudes than for “V” or wedge-shaped flames. These predictions are shown to be consistent with available experimental data.     

Oscillatory cool flame combustion behavior of submillimeter sized n-alkane droplet under near limit conditions

This paper reports simulation results of oscillatory cool flame burning of an isolated, submillimeter sized n-heptane (n-C₇H₁₆) droplet in a selectively ozone (O₃) seeded nitrogen-oxygen (N₂-O₂) environments at atmospheric pressure. An evolutionary one-dimensional droplet combustion code encompassing relevant physics and detailed chemistry was employed to explore the roles of low-temperature chemistry, O₃ seeding, and dynamic flame structure on burning behaviors. For X_O2= 21% and a range of selective ozone seeding, near-quasi-steady cool flame burning is achieved directly (without requiring hot flame initiation and radiative extinction). Under low oxygen index conditions, but with significant O₃ seeding (X_O3 = 5%), a nearly quasi-steady cool flame is initially established that then transitions to a dynamically oscillating cool flame burning mode which continues until the droplet is completely consumed. It is found that the oscillation occurs as result of a initial depletion of fuel vapor-oxidizer layer evolving near the droplet surface and its dynamic re-establishment through liquid vaporization and vapor/oxidizer transport. A kinetic analysis indicates that the dynamic competition between the reaction classes- (a) degenerate chain branching and (b) chain termination/propagation - along with continuous fuel and oxygen leakage through the flame location contributes to an oscillatory burning phenomena of ever-increasing amplitude. Analysis based on single full-cycle of oscillatory burning shows that the reaction progression matrices (evolution of heat and species) for QOOH➔chain propagation/termination reactions (here, Q = C₇H₁₄-) directly scales with the gas phase temperature field. On the contrary, the QOOH➔degenerate branching reactions undergoes three distinct stages within the same oscillatory cycle. The coupled flame dynamics and kinetics suggest that in the oscillatory burning mode, kinetic processes dynamically cross through conditions characterizing the negative temperature coefficient (NTC) turnover temperature, separating low temperature and NTC kinetic regimes. In addition, a parametric study is conducted to determine the role of O₃ seeding level on the observed oscillation phenomena.     

Hyperacceleration effects on turbulent combustion in premixed step-stabilized flames

Experimental results are presented from an investigation of the effects of large transverse accelerations on flame propagation and blowout limits in premixed step-stabilized flames. The accelerations, which exceed ±10,000 g in the present study, induce large body forces on the high-density reactants and low-density products. These body forces can substantially alter the flame propagation mechanisms and dramatically increase the flame blowout limits. Sustained centripetal accelerations a_c ≡ U²/R are created by flowing a premixed propane–air reactant stream with equivalence ratios 0.7  Φ  1.9 at various speeds U through a semicircular channel with radius R. A backward-facing step of height h on the radially outer (a_c > 0) or inner (a_c < 0) wall stabilizes the flame. For a_c > 0 the acceleration acts to force high-density reactants into the recirculation zone and low-density products into the reactant stream, while a_c < 0 forces hot products into the recirculation zone and impedes cold reactants from entering this zone. An otherwise identical straight channel provides corresponding baseline (a_c = 0) results for comparison. The flow speed U, equivalence ratio Φ, and step height h are systematically varied for a_c = 0, a_c > 0, and a_c < 0. Shadowgraph and chemiluminescence imaging show that as a_c→ +∞ the propagation of the flame across the channel becomes independent of the flame burning velocity and instead is primarily due to large-scale “centrifugal pumping” driven by the induced body forces. For a_c → −∞ the body forces effectively segregate reactants and products to produce a nearly flat flame. In both cases, for large |a_c| values the resulting blowout limits can be substantially higher than those at a_c = 0.     

Karhunen–Loève decomposition of peripheral blood flow signal

The Karhunen–Loève expansion is applied to scalar signals and the effect of window length (t_w), time lag (τ) and embedding dimension (d) is analysed for periodic signals and for signals modeled by the Lorenz equations. For τ≠k/2f_i (f_i are characteristic frequencies of the signal, k is positive integer), we obtain 2m modes from an m-periodic signal. For a large set of parameters a finite number of modes was not obtained from the Lorenz system. It is further shown that, on the time scale of a minute, the peripheral blood flow signal contains oscillatory modes that occur in pairs thereby confirming that the blood flow through the cardiovascular system is oscillatory. Some of the difficulties of applying Karhunen–Loève expansion to scalar signals are pointed out.     

Combustion characteristics and detailed simulations of surrogates for a Tier II gasoline certification fuel

An experimental and numerical study of combustion of a gasoline certification fuel (‘indolene’), and four (S4) and five (S5) component surrogates for it, is reported for the configurations of an isolated droplet burning with near spherical symmetry in the standard atmosphere, and a single cylinder engine designed for advanced compression ignition of pre-vaporized fuel. The intent was to compare performance of the surrogate for these different combustion configurations and to assess the broader applicability of the kinetic mechanism and property database for the simulations. A kinetic mechanism comprised of 297 species and 16,797 reactions was used in the simulations that included soot formation and evolution, and accounted for unsteady transport, liquid diffusion inside the droplet, radiative heat transfer, and variable properties. The droplet data showed a clear preference for the S5 surrogate in terms of burning rate. The simulations showed generally very good agreement with measured droplet, flame, and soot shell diameters. Measurements of combustion timing, in-cylinder pressure, and mass-averaged gas temperature were also well predicted with a slight preference for the S5 surrogate. Preferential vaporization was not evidenced from the evolution of droplet diameter but was clearly revealed in simulations of the evolution of mixture fractions inside the droplets. The influence of initial droplet diameter (D_o) on droplet burning was strong, with S5 burning rates decreasing with increasing D_o due to increasing radiation losses from the flame. Flame extinction was predicted for D_o =3.0 mm as a radiative loss mechanism but not predicted for smaller D_o for the conditions of the simulations.     

Extinguishment of propane/air co-flowing diffusion flames by fine water droplets

In the present study, extinguishment of propane/air co-flowing diffusion flame by fine water droplets was investigated experimentally. Water droplets are generated by piezoelectric atomizers with the maximum droplets flow rate of 1500 ml/h. When the fuel injection velocity U_f is low, an attached laminar diffusion flame with a premixed flame at the base is stabilized. At some distance from the burner rim, a transition from laminar to turbulent diffusion flame occurs, and a turbulent diffusion flame is formed in the downstream region. When the fuel injector rim is thin (δ = 0.5 mm), the flame stability deteriorates with increase of the co-flowing air stream velocity U_a and the water droplets flow rate Q_m. The stability mechanism can be explained by the balance of the gas velocity and the burning velocity of premixed flame formed at the base. However, when the injector rim is thick (δ = 5 mm), a recirculation zone is produced downstream of the injector rim. The dependence of the quenching distance H_q on U_f and Q_m is relatively weak, and the stability diagram shows curious features. It was shown that U_a is crucially important since it determines flow residence time; if U_a < 0.4 m/s, water droplets can evaporate when they go by the recirculation zone, and the water vapor can diffuse into the recirculation zone. However, if U_a > 0.4 m/s, the water droplets should pass by the recirculation zone without sufficiently evaporated and are not so effective to extinguish the flame. The supply velocity of droplet-laden air should be low enough so that water droplets can evaporate and water vapor can diffuse into the premixed region at the base to obtain sufficient effectiveness of water droplets for fire suppression.     

Enhancement of the Debye frequency due to a charged-boson exchange model of high temperature superconductors

A simple model of 1-2-3 superconductors in which electrons (holes) in CuO₂ planes interact via exchange with two kinds of bosons is considered. Namely, via one-phonon exchange (weak coupling-Cooper pairing), and via paired holes on oxygen O⁰ from Cu-O chains. The mechanisms of paired holes exchange (“charged bosons”-“O⁰” exchange) considered here in strong coupling leads to the enhancement of the Fröhlich constant g_f (g²_F→Kg²F), and as a consequence to the enhancement of the Debye frequency ω_D^K=f^Kω_D, f^K 1. In the proposed model the exact expression for the constant K is derived.     

Fundamental burning velocities of meso-scale propagating spherical flames with H2, CH4 and C3H8 mixtures

This study is performed to experimentally examine the fundamental burning velocity characteristics of meso-scale outwardly propagating spherical laminar flames in the range of flame radius r_f approximately from 1 to 5 mm for hydrogen, methane and propane mixtures, in order to make clear a method for improving combustion of micro–meso scale flames. Macro-scale laminar flames with r_f > 7 mm are also examined for comparison. The mixtures have nearly the same laminar burning velocity (S_L0 = 25 cm/s) for unstretched flames and different equivalence ratios ?. The radius r_f and the burning velocity S_Ll of meso-scale flames are estimated by using sequential schlieren images recorded under appropriate ignition conditions. It is found that S_Ll of hydrogen and methane premixed meso-scale flames at the same r_f or the Karlovitz number Ka shows a tendency to increase with decreasing ?, whereas S_Ll of propane flames increases with ?. However, S_Ll tends to decrease with the Lewis number Le and the Markstein number Ma, irrespective of the type of fuel and ?. It also becomes clear that the optimum flame size and Ka to improve the burning velocity exist for some mixtures depending on Le and fuel types.     

Concentration-dependent luminescence of Tb ions in high calcium aluminosilicate glasses

This study deals with the results on the concentration-dependent fluorescence properties of Tb³⁺-doped calcium aluminosilicate (CAS) glasses of composition (100−x)(58SiO₂–23CaO–5Al₂O₃–4MgO–10NaF in mol%)-x Tb₂O₃ (x=0, 0.25, 0.5, 1, 2, 4, 8, 16, 24, 32, 40 in wt%). The FTIR reflectance spectra suggested the role of dopant ions as network modifiers in the glass network. The fluorescence spectra of low Tb³⁺-doped glasses have revealed prominent blue and green emissions from ⁵D₃ and ⁵D₄ excited levels to ⁷F_j ground state multiplet, respectively. The glass with 2 wt% of Tb₂O₃ has exhibited maximum intensity of blue emission from ⁵D₃ level, while green emission from ⁵D₄ level has increased linearly up to 24 wt% and showed reduction in the rate of increase for higher Tb₂O₃ concentrations. The concentration quenching of blue emission (⁵D₃→⁷F_j) is attributed mainly to the resonant energy transfer (RET) assisted cross-relaxation (CR) among the excited and nearest neighbour unexcited Tb³⁺ ions in the glass matrix. The decline in rate of increase of green emission (⁵D₄→⁷F_j) at higher concentrations has been explained due to a possible occurrence of cooperative energy transfers leading to 4f⁸→4f⁷5d transition interactions. The blue and green emission decay kinetics have been recorded to compute the excited level (⁵D₃ and ⁵D₄) lifetimes, which confirmed the Tb³⁺ concentration quenching of the blue emission in these glasses.     

Burning velocity correlation of methane/air turbulent premixed flames at high pressure and high temperature

Turbulent burning velocities for methane/air mixtures at pressures ranging from atmospheric pressure up to 1.0 MPa and mixture temperatures of 300 and 573 K were measured, which covers the typical operating conditions of premixed-type gas-turbine combustors. A bunsen-type flame stabilized in a high-pressure chamber was used, and OH-PLIF visualization was performed with the pressure and mixture temperature being kept constant. In addition to a burner with an outlet diameter of 20 mm for the high-pressure experiments, a large-scale burner with an outlet diameter of 60 mm was used at atmospheric pressure to extend the turbulence Reynolds number based on the Taylor microscale, R_λ, as a common parameter to compare the pressure and temperature effects. It was confirmed that R_λ over 100 could be attained and that u′/S_L could be extended even at atmospheric pressure. Based on the contours of the mean progress variable c = 0.1 determined using OH-PLIF images, turbulent burning velocity was measured. S_T/S_L was also found to be greatly affected by pressure for preheated mixtures at 573 K. The bending tendency of the S_T/S_L curves with u′/S_L was seen regardless of pressure and mixture temperature and the R_λ region where the bending occurs corresponded well to the region where the smallest scale of flame wrinkling measured as a fractal inner-cutoff approaches the characteristic flame instability scale and becomes almost constant. A power law of S_T/S_L with (P/P₀)(u′/S_L) was clearly seen when S_T was determined using c = 0.1 contours, and the exponent was close to 0.4, indicating agreement with the previous results using the mean flame cone method and the significant pressure effects on turbulent burning velocity.     

Premixed turbulent flame front structure investigation by Rayleigh scattering in the thin reaction zone regime

Premixed turbulent flames of methane–air and propane–air stabilized on a bunsen type burner were studied using planar Rayleigh scattering and particle image velocimetry. The fuel–air equivalence ratio range was from lean 0.6 to stoichiometric for methane flames, and from 0.7 to stoichiometric for propane flames. The non-dimensional turbulence rms velocity, u′/S_L, covered a range from 3 to 24, corresponding to conditions of corrugated flamelets and thin reaction zones regimes. Flame front thickness increased slightly with increasing non-dimensional turbulence rms velocity in both methane and propane flames, although the flame thickening was more prominent in propane flames. The probability density function of curvature showed a Gaussian-like distribution at all turbulence intensities in both methane and propane flames, at all sections of the flame.The value of the term Dκ, the product of molecular diffusivity evaluated at reaction zone conditions and the flame front curvature, has been shown to be smaller than the magnitude of the laminar burning velocity. This finding questions the validity of extending the level set formulation, developed for corrugated flames region, into the thin reaction zone regime by increasing the local flame propagation by adding the term Dκ to laminar burning velocity.     

Stabilization mechanism of turbulent premixed flames in strongly swirled flows

The stabilization of turbulent premixed flames in strongly swirled flows undergoing vortex breakdown is studied in the case of the ALSTOM En-Vironmental (EV) double cone burner using a simple one-dimensional boundary layer type model and computational fluid dynamics, mainly at the level of large-eddy simulation. The analysis shows that, due to flame curvature effects, the flame speed on the combustor axis is 2 D _t/R _F lower than the turbulent burning rate, where D _t is a characteristic turbulent diffusion coefficient and R _F the flame radius of curvature. Flame propagation with negative speed observed in the experiments, i.e. the flame completely embedded in the central recirculation zone on the symmetry axis, is explained with the one-dimensional model as caused by the factor 2 D _t/R _F being larger than the characteristic turbulent burning rate. A peculiar sudden displacement of the flame anchoring location deep into the burner, which takes place experimentally at a critical value of the equivalence ratio, cannot however be explained with the present one-dimensional approach due to the modelling assumptions. The mathematical analysis is supported in this case with large-eddy simulation which can accurately reproduce the flame behaviour across the full operating range. It is finally shown that steady RANS methods cannot cope with the problem due to their inability to correctly predict the velocity flowfield in this burner.     

Regression fronts in random sphere packs: Application to composite solid propellant burning rate

The burning rate of a composite solid propellant may be estimated by global modeling, such as the widely used BDP model. The backbone of such models is the “mixture law” that links the propellant burning rate r_p with the burning rate of its own components, i.e., oxidizer r_ox and binder r_b. However, different laws are available in literature which all read: 1/r_p = q(ξ)/r_ox + (1 − q(ξ))/r_b, with q(ξ) a function of oxidizer volume fraction ξ. This work attempts in analyzing numerically the validity of those empirical formulations by surface regression computation. Composite propellants are modeled by a random packing of monomodal spheres and the evolution of the regression front is computed via the resolution of Hamilton–Jacobi equations. It is shown that the popular choice q(ξ) = ξ is fairly valid but only provided that burn rate ratio Z = r_ox/r_b is about 1. When Z > 1, combustion surface is no longer plane and global burning rate deviates from postulated laws. A special regime is also noticed for high rate ratio Z (typically Z > 5) because combustion then preferentially takes place through adjacent oxidizer particles. Computed results occur to be correctly modeled by percolation theory. This hints that percolation is a common feature of propellant combustion and a critical percolation threshold on volume fraction is numerically found to be about ξ_c  0.2. First validations show encouraging correlations with experimental data.     

The response of turbulent premixed flames to normal acoustic excitation

To model the thermo-acoustic excitation of flames in practical combustion systems, it is necessary to know how a turbulent flame front responds to an incident acoustic wave. This will depend partly on the way in which the burning velocity responds to the wave. In this investigation, the response of CH₄/air and CH₄/H₂/air mixtures has been observed in a novel flame stabilisation configuration, in which the premixture of fuel and air is made to decelerate under controlled conditions in a wide-angle diffuser. Control is provided by an annular wall-jet of air and by turbulence generators at the inlet. Ignition from the outlet of the diffuser allows an approximately flat flame to propagate downwards and stabilise at a height that depends on the turbulent burning velocity. When the flow is excited acoustically, the ensemble-averaged height oscillates. The fluctuations in flow velocity and flame height are monitored by phase-locked particle image velocimetry and OH-planar laser induced fluorescence, respectively. The flame stabilised against a lower incident velocity as the acoustic amplitude increased. In addition, at the lowest frequency of 52 Hz, the fluctuations in turbulent burning velocity (as represented by the displacement speed) were out-of-phase with the acoustic velocity. Thus, the rate of displacement of the flame front relative to the flow slowed as the flow accelerated, and so the flame movement was bigger than it would have been if the burning velocity had not responded to the acoustic fluctuation. With an increase in frequency to 119 Hz, the relative flame movement became even larger, although the phase-difference was reduced, so the effect on burning velocity was less dramatic. The addition of hydrogen to the methane, so as to maintain the laminar burning velocity at a lower equivalence ratio, suppressed the response at low amplitude, but at a higher amplitude, the effect was reversed.