Computational and experimental study of JP-8, a surrogate, and its components in counterflow diffusion flames

Non-sooting counterflow diffusion flames have been studied both computationally and experimentally, using either JP-8, or a six-component JP-8 surrogate mixture, or its individual components. The computational study employs a counterflow diffusion flame model, the solution of which is coupled with arc length continuation to examine a wide variety of inlet conditions and to calculate extinction limits. The surrogate model includes a semi-detailed kinetic mechanism composed of 221 gaseous species participating in 5032 reactions. Experimentally, counterflow diffusion flames are established, in which multicomponent fuel vaporization is achieved through the use of an ultrasonic nebulizer that introduces small fuel droplets into a heated nitrogen stream, fostering complete vaporization without fractional distillation. Temperature profiles and extinction limits are measured in all flames and compared with predictions using the semi-detailed mechanism. These measurements show good agreement with predictions in single-component n-dodecane, methylcyclohexane, and iso-octane flames. Good agreement also exists between predicted and measured variables in flames of the surrogate, and the agreement is even better between the experimental JP-8 flames and the surrogate predictions.     

Numerical investigation of the effects of polydisperse water sprays on extinction conditions of counterflow methane non-premixed flames

Water, sprayed in the form of tiny droplets, has emerged as a potential fire suppressant after the halon compounds such as trifluorobromomethane (CF₃Br, Halon 1301) were banned by the Montreal protocol. The size distribution of the water droplet plays a crucial role in the effectiveness of the water spray in fire suppression. A numerical investigation of the influence of size distribution of a polydisperse water spray on extinction of counterflow diffusion flames is presented in this paper. This study uses laminar finite rate model with reduced CHEMKIN chemistry for numerical simulations. The discrete phase, namely the water spray, is simulated using Lagrangian Discrete Phase Modelling approach. In this work, the polydispersity of water spray is taken into account in the numerical simulation by a suitable Rosin–Rammler distribution. Results obtained from numerical simulation are validated with the experimental results reported in the literature. This study demonstrates that the representation of the polydisperse spray by a monodisperse spray (with droplet diameter same as the SMD of the polydisperse spray) in numerical simulations is not always justified and it leads to deviation from the experimental results. The effects of number mean diameter and spread parameter on the efficacy of flame suppression are investigated for polydisperse sprays. A comprehensive comparison is done between the effectiveness of monodisperse and polydisperse water sprays. An optimum droplet diameter is obtained for monodisperse sprays for which the effectiveness of the spray is maximum. The effects of evaporation Damköhler number and Stokes number of water droplets on flame suppression have also been explained.     

A three mixture fraction flamelet model for multi-stream laminar pulverized coal combustion

A three mixture fraction flamelet model is proposed for multi-stream laminar pulverized coal combustion. The technique of coordinate transformation is utilized to map the flamelet solutions from a unit pyramid space into a unit cubic space to improve the stability of the simulation. The validity of the three mixture fraction flamelet model was assessed on different configurations, including a laminar counterflow pulverized coal/methane flame and a laminar piloted pulverized coal jet flame. The flamelet predictions were compared to the reference results of the detailed chemistry solutions. For the counterflow flame, it was found that the flame temperature and major species mass fractions are correctly predicted by the three mixture fraction flamelet model. However, discrepancies are observed for combustion-mode-sensitive species such as CO and H₂ in the premixed combustion region. The thermo-chemical quantities in the char surface reaction zone cannot be correctly predicted if the mixing between the char off-gas stream and other streams is neglected. For the piloted jet flame, it was shown that the stable thermo-chemical variables can be correctly predicted at the upper and middle stream locations. However, at the downstream location, discrepancies can be observed in certain regions. Overall, the validity of the three mixture fraction flamelet model for multi-stream pulverized coal combustion is confirmed and its performance in turbulent pulverized coal combustion will be tested in future work.     

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.     

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.     

Transient interactions between a premixed double flame and a vortex

The interaction between a laminar flame and a vortex is an important study for understanding the fundamentals of turbulent combustion. In the past, however, flame-vortex interactions have been investigated only for high-temperature flames. In this study, the impact of a vortex on a premixed double flame, which consists of a coupled cool flame and a hot flame, is examined experimentally and computationally using dimethyl ether/oxygen/ozone mixtures. The double flame is first shown to occur near the extinction limit of the hot flame. The differences between steady-state cool flames, double flames, and hot flames are explored in a one-dimensional counterflow configuration. The transient interactions between double flames and impinging vortices are then investigated experimentally using a micro-jet and numerically in two-dimensional transient modeling. It is seen that the vortex can extinguish the near-limit hot flame locally, resulting in a lone cool flame. At higher vortex intensities, the cool flame may also be extinguished after the extinction of the hot flame. It is found that there can be three different transient flame structures coexisting at the same time: an extinguished flame hole, a cool flame, and a double flame. Moreover, flame curvature is shown to play an important role in determining whether the vortex weakens or strengthens the cool flame and double flame.     

Forced Convection with Laminar Pulsating Counterflow in a Saturated Porous Circular Tube

An analytical solution is obtained for forced convection in a circular tube occupied by a core–sheath-layered saturated porous medium with counterflow produced by pulsating pressure gradients. The case of the constant heat-flux boundary conditions is considered, and the Brinkman model is employed for the porous medium. A perturbation approach is used to obtain analytical expressions for the velocity, temperature distribution, and transient Nusselt number for convection produced by an applied pressure gradient that fluctuates with small amplitude harmonically in time about a non-zero mean. It is shown that the fluctuating part of the Nusselt number alters in magnitude and phase as the dimensionless frequency increases. The magnitude increases from zero, goes through a peak, and then decreases to zero. The height of the peak depends on the values of various parameters. The phase (relative to that of the steady component) decreases from π/2 to − π/2 as the frequency increases.     

Isotachophoresis with counterflow in an open capillary: Computer simulation and experimental validation

The purpose of applying a countercurrent flow to isotachophoretic migration is to increase the effective separation channel length during ITP. However, severe dispersion induced by applying a counterflow can be detrimental to ITP. This paper uses numerical simulations in a 2D axisymmetric domain to investigate the dispersion caused by a parabolic counterflow in open‐capillary ITP. Counterflow in these simulations was generated by applying a back pressure to stop the isotachophoretic stack, i.e., forming stationary ITP zones. It is found that dispersion is strongly related to analyte molecular diffusivity: R‐phycoerythrin, due to its small diffusivity, showed ～20‐fold increase in zone width in stationary counterflow ITP, compared to ITP in the absence of counterflow, while fluorescein only had ～10% increase in zone width under similar operating conditions. Applying the Taylor–Aris dispersion formula in counterflow ITP simulations provided only a rough estimate of the dispersion, e.g., overestimation of analyte zone widths. Experiments on counterflow ITP were conducted in a silica capillary that was covalently and dynamically coated to exclude electroosmosis effect. The counterflow was generated by adjusting the relative height of the fluids in the two reservoirs at the capillary ends. Good qualitative agreement between simulations and experiments was found.     

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.     

Bilayer counterflow transport at filling factor 1 in the strong interacting regime

We review magneto-transport properties of interacting GaAs bilayer hole systems, with very small inter-layer tunneling, in a geometry where equal currents are passed in opposite directions in the two, independently contacted layers (counterflow). In the quantum Hall state at total bilayer filling ν=1 both the longitudinal and Hall counterflow resistances tend to vanish in the limit of zero temperature, suggesting the existence of a superfluid transport mode in the counterflow geometry. As the density of the two layers is reduced, making the bilayer more interacting, the counterflow Hall resistivity (ρ_xy) decreases at a given temperature while the counterflow longitudinal resistivity (ρ_xx), which is much larger than ρ_xy, hardly depends on density. Our data suggest that the counterflow dissipation present at any finite temperature is a result of mobile vortices in the superfluid created by the ubiquitous disorder in this system.