Microgravity droplet combustion: effect of tethering fiber on burning rate and flame structure

Droplets tethering on fibers has become a well established technique for conducting droplet combustion experiments in microgravity conditions. The effects of these supporting fibers are frequently assumed to be negligible and are not considered in the experimental analysis or in numerical simulations. In this work, the effect of supporting fibers on the characteristics of microgravity droplet combustion has been investigated numerically; a priori predictions have then been compared with published experimental data. The simulations were conducted using a transient one-dimensional spherosymmetric droplet combustion model, where the effect of the supporting fiber was implicitly taken into account. The model applied staggered convective flux finite volume method combined with high-order implicit time integration. Thermal radiation was evaluated using a statistical narrow band radiation model. Chemical kinetics and thermophysical properties were represented in rigorous detail. Tether fiber diameter, droplet diameter, ambient pressure and oxygen concentration were varied over a range for n-decane droplets in the simulations. The results of the simulations were compared to previously published experiments conducted in the Japan Microgravity Center (JAMIC) 10 second drop tower and the NASA Glenn Research Center (GRC) 5.2 second drop tower. The model reproduces closely nearly all aspects of tethered n-decane droplet burning phenomena, which included droplet burning history, transient and average burning rate, and flame standoff ratio. The predictions show that the presence of the tethering fiber significantly influences the observed burning rate, standoff ratio, and extinction.     

Asymptotic analysis of quasi-steady n-heptane droplet combustion supported by cool-flame chemistry

A skeletal chemical-kinetic mechanism for n-heptane cool flames is simplified to the maximum extent possible by introduction of steady-state approximations for intermediaries, following procedures employed previously in addressing two-stage ignition. A pair of ordinary differential equations in mixture-fraction space is thereby obtained, describing the quasi-steady structures of the temperature and heptylketohydroperoxide fields. Application of activation-energy asymptotics for the partial-burning regime to this pair of equations is shown to provide convenient expressions for flame structures and the extinction condition associated with maximally reduced chemistry. With the mixture-fraction co-ordinate related to radius, these results are used to address droplet-combustion experiments that have been performed in the International Space Station. Droplet diameters at extinction are predicted as functions of the oxygen concentration in the atmosphere and are compared with experiment. While the results are encouraging concerning qualitative predictions of dependences of extinction diameters on atmospheric conditions, there are noticeable quantitative differences that point to deficiencies in the analysis, likely resulting from a number of oversimplifications. Further investigation is therefore recommended.     

Identification of droplet burning modes in lean, partially prevaporized swirl-stabilized spray flames

Experimental and numerical investigations of single droplet burning modes in a lean, partially prevaporized swirl-stabilized spray flame are reported. In the experiment single droplet flames have been visualized by CH-PLIF and simultaneous recording of the Mie signal. Two single droplet burning modes were identified: the envelope flame is a spherical diffusion flame burning at near-stoichiometric conditions. The wake flame is a potentially lean, partially premixed flame located downstream of the droplet. The droplet burning mode is of practical relevance, since it has significant impact on NO formation due to incomplete prevaporization.The droplet burning mode is determined by the ratio of chemical and convective time scales. The convective time scale is related to the droplet slip velocity. The impact of turbulent gas phase velocity fluctuations on droplet mechanics and droplet burning is discussed, based on a previous numerical investigation. In the present study the droplet slip velocity was measured with the 3D Phase Doppler (3D-PD) technique. For the measured slip velocities and ambient conditions in the hot gas region of the spray flame, simulations of single droplet burning were performed utilizing detailed models for chemical reaction, diffusive transport and vaporization. An agreement between the droplet burning modes predicted by the simulation and the droplet burning modes observed in the experiments was found.     

A hybrid droplet vaporization-chemical surrogate approach for emulating vaporization,physical properties,and chemical combustion behavior of multicomponent fuels

The complex nature of multicomponent aviation fuels presents a daunting task for accurately simulating combustion behavior without incurring impractical computational costs. To reduce computation time, chemical fuel surrogates comprised of only a few species are used to emulate the combustion of complex pre-vaporized fuels. These surrogates are often unable to match the vaporization behavior and physical properties of the real fuel and fail to capture the effect of preferential vaporization on combustion behavior. Therefore, a computationally efficient, hybrid droplet vaporization-chemical surrogate approach has been developed which emulates both the physical and chemical properties of a multicomponent kerosene fuel. The droplet vaporization/physical portion of the hybrid uses the Coupled Algebraic–Direct Quadrature Method of Moments with delumping to accurately solve for the evolution of every discrete species in a vaporizing fuel droplet with the computational efficiency of a continuous thermodynamic model. The chemical surrogate portion of the hybrid is linked to the vaporization model using a functional group matching method, which creates an instantaneous surrogate composition to match the distribution of chemical functional groups (CH₂, (CH₂)_n, CH₃ and Benzyl-type) in the vaporization flux of the full fuel. The result is a hybrid method which can accurately and efficiently predict time-dependent, distillation-resolved combustion property targets of the vaporizing fuel and can be used to investigate the effects of preferential vaporization on combustion behavior.     

Ignition and combustion of a falling, single sodium droplet

Ignition and combustion of a single sodium droplet has been studied experimentally, by use of a falling droplet. It is found that the ignition delay time increases first gradually and then rapidly, with decreasing initial temperature and/or oxygen concentration, and reaches the limit of ignitability, because of the suppression of surface reaction in the ignition stage. It is also found that with decreasing droplet diameter, the ignition delay time first decreases gradually, because of the decrease in the droplet mass to be heated, and then increases steeply, because of the enhancement of heat loss from the droplet surface. As for the effect of the relative speed, it is found that the ignition delay time increases with increasing relative speed, because of the enhanced heat loss. Experimental comparisons with the analytical results have also been conducted to elucidate dominant parameters, and it is confirmed that a set of comprehensive parameters in the literature can be useful in correlating dominant parameters that influence the ignition delay and/or the limit of ignitability. Furthermore, the analysis has been extended to determine the critical size for the ignition and that for the minimum ignition delay time. Combustion behavior after the ignition has also been examined, and it is found that d²-law can hold for the sodium droplet combustion. In addition, it is found that the burning rate-constant without forced convection has nearly the same value as those for usual hydrocarbon droplets, while the sodium combustion in air is quite similar to that of the usual hydrocarbon fuel in an oxidizer-rich environment.     

Observation of droplet motion during flame spread on three-fuel-droplet array with a pendulum suspender

Flame spread on a fuel droplet array has been studied as a simple model of spray combustion. A three-fuel-droplet array with a pendulum suspender was employed to investigate interactions between flame spread and droplet motion in the axial direction. Initial droplet diameter was 0.8 mm, and fuel was n-heptane. A silicon carbide pendulum suspender of 15 μm in diameter and 30 mm in length was used for the third droplet. The first fixed droplet was ignited by electric spark. Behavior of the flame and the third droplet was observed using a high-speed video camera with an image intensifier. Particle tracking velocimetry (PTV) measurements were performed to explain the behavior of the third movable droplet. The dimensionless droplet span, which is the average of droplet-to-droplet distances divided by the average initial diameter of the three droplets, was varied from 2.5 to 8 for observing flame spread, and fixed at 5.5 for PTV measurements. It was observed that the third droplet moved away from the second droplet before the flame spread to the third droplet. The displacement of the third droplet is remarkable when the dimensionless droplet span is close to the limit of flame spread. This implies that the movement of the droplet decreases the dimensionless span of the flame spread limit and the flame spread speed near the flame spread limit. Results of PTV measurements suggest that the heat expansion wave, caused by ignition of the premixture which was accumulated around the second droplet, and the burned gas flow from the second droplet pushed away the third droplet; then natural convection, induced by the flames of the first and second droplets, drew the third droplet to the second droplet. The heat expansion wave and the burned gas flow of the second droplet reached nearly 12 in dimensionless span.     

A semi-detailed kinetic model of char combustion with consideration of thermal annealing

The present paper presents a semi-detailed kinetic model of coal char combustion which embodies consideration of thermal annealing as a mechanism leading to the loss of char combustion reactivity along burn off. The distinctive feature of this model is that deactivation induced by thermal annealing is followed along with combustion. Thermodeactivation is modelled according to the power-law equation proposed by Senneca and Salatino [1]. A semi-detailed combustion mechanism was taken after Hurt and Calo [2] and includes three steps: formation of carbon–oxygen complexes (chemisorption), switch-over of surface oxides and desorption of oxygen complexes to yield combustion products. Computation results allow to discuss the impact of thermal annealing on char combustion under conditions of practical interest.     

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.     

Cobalt sorption properties of MgO prepared by solution combustion

Water effect on the combustion preparation of MgO is presented. The obtained materials are characterized through their specific surface area, morphology, particle shape, fractal dimension and Co²⁺ sorption. The surface fractal dimension of the combustion prepared sample was 2.0 but in those where water was included it decreased to 1.8. In the sample prepared by calcination it was 2.3. A linear correlation between the fractal dimension and Co²⁺ sorption was found.     

Synthesis of bismuth titanate with urea as fuel by solution combustion route and its dielectric and ferroelectric properties

Preparation of ferroelectric bismuth titanate (Bi₄Ti₃O₁₂) is carried out by solution combustion route with urea as fuel at much lower calcinations temperatures. The single phase bismuth titanate was obtained after calcinations at 800 °C. SEM micrographs of the calcined powders show agglomerated, flaky and foamy morphology, which is typical of combustion synthesis and that of sintered ceramics shows the grain formation. Behavior of dielectric constant and dielectric loss as a function of temperature of as-prepared sample are reported in this communication. Ferroelectric to paraelectric phase transition occurs at the temperature T_c ∼ 660 °C. Its remnant polarization (2P_r) is very less of the order of 0.012 μC/cm².     

超疏水表面液滴的振动特性及其与液滴体积的关系

超疏水表面液滴的振动特性与接触线的移动、液滴体积、基底振幅等因素密切相关.本文在基底振幅较小且恒定的条件下,研究了超疏水表面液滴的共振振幅、模式区间、共振频率等振动特性及其与液滴体积(20—500μL)的关系.此外,将基于一般性疏水表面建立的Noblin共振频率计算模型应用于超疏水表面,并提出虚驻点的概念,借此对模...     

Dielectric properties of Bismuth Titanate (Bi4Ti3O12) synthesized using solution combustion route

Ferroelectric Bismuth Titanate (Bi₄Ti₃O₁₂) was prepared by solution combustion route with glycine as fuel. The single phase Bismuth Titanate was obtained after calcination at 800 °C, which was confirmed with the help of X-ray diffraction studies and EDS analysis. SEM micrographs of the calcined powders show agglomerated particles, which is typical of combustion synthesis. Behavior of dielectric constant and dielectric loss as a function of temperature of as prepared sample are reported here. Ferroelectric to paraelectric phase transition occurs at the temperature T_c∼650 °C. Impedance studies were made in the frequency range from 1 KHz to 1 MHz. The semicircles observed in the complex impedance diagrams indicate deviation from the Debye behavior. Activation energy of the sample around T_c is found to be ∼0.35 eV and below T_c is ∼0.13 eV, which was calculated using the Arrhenius plots.     

Phase transitions of perfluorocarbon nanoemulsion induced with ultrasound: A mathematical model

While ultrasound has been used in many medical and industrial applications, only recently has research been done on phase transformations induced by ultrasound. This paper presents a numerical model and the predicted results of the phase transformation of a spherical nanosized droplet of perfluorocarbon in water. Such a model has applications in acoustic droplet vaporization, the generation of gas bubbles for medical imaging, therapeutic delivery and other biomedical applications.The formation of a gas phase and the subsequent bubble dynamics were studied as a function of acoustic parameters, such as frequency and amplitude, and of the physical aspects of the perfluorocarbon nanodroplets, such as chemical species, temperature, droplet size and interfacial energy. The model involves simultaneous applications of mass, energy and momentum balances to describe bubble formation and collapse, and was developed and solved numerically. It was found that, all other parameters being constant, the maximum bubble size and collapse velocity increases with increasing ultrasound amplitude, droplet size, vapor pressure and temperature. The bubble size and collapse velocity decreased with increasing surface tension and frequency. These results correlate with experimental observations of acoustic droplet vaporization.     

Large eddy simulation of piloted pulverised coal combustion using extended flamelet/progress variable model

An extended flamelet/progress variable (EFPV) model for simulating pulverised coal combustion (PCC) in the context of large eddy simulation (LES) is proposed, in which devolatilisation, char surface reaction and radiation are all taken into account. The pulverised coal particles are tracked in the Lagrangian framework with various sub-models and the sub-grid scale (SGS) effects of turbulent velocity and scalar fluctuations on the coal particles are modelled by the velocity-scalar joint filtered density function (VSJFDF) model. The presented model is then evaluated by LES of an experimental piloted coal jet flame and comparing the numerical results with the experimental data and the results from the eddy break up (EBU) model. Detailed quantitative comparisons are carried out. It is found that the proposed model performs much better than the EBU model on radial velocity and species concentrations predictions. Comparing against the adiabatic counterpart, we find that the predicted temperature is evidently lowered and agrees well with the experimental data if the conditional sampling method is adopted.     

Oscillatory thermal-diffusive instability of combustion waves in a model with chain-branching reaction and heat loss

In this paper we investigate the properties and the linear stability of premixed combustion waves in a non-adiabatic thermal-diffusive model with a two-step chain-branching reaction mechanism. Here we focus only on the emergence of the pulsating instabilities, and the stability analysis is carried out for Lewis numbers for fuel greater than one, and various values of Lewis number for radicals. We consider the problem in two spatial dimensions to allow perturbations of a multidimensional nature. It is demonstrated that the flame speed as a function of the parameters is a double-valued C-shaped function, i.e. for a given set of parameter values there are either two solutions, fast and slow solution branches, propagating with different speed, or the combustion wave does not exist. The extinction of combustion waves occurs at finite values of the parameters and non-zero flame speed. The flame structure demonstrates a slow recombination regime behaviour with negligible fuel leakage for the fast solution branch away from the extinction condition. For parameter values close to the extinction condition and on the slow solution branch, the fuel leakage is significant and a fast recombination regime is observed. It is demonstrated that two types of instabilities emerge in the model: the uniform planar and the travelling instability. The slow solution branch is always unstable due to the uniform perturbations. The fast solution branch is either stable or loses stability due to the travelling or uniform perturbations. The switching between the onset of various regimes of instability is due to the bifurcation of co-dimension two. In the adiabatic limit this bifurcation is found for Lewis number for fuel equal to one, whereas in the non-adiabatic case it moves towards values above unity. The properties of the travelling instability are studied in detail.     

Optical properties of InGaN/GaN double quantum wells with varying well thickness

The optical properties and recombination kinetics of the InGaN/GaN double quantum well (DQW) structures with different well thickness (L_w) have been studied by means of photoluminescence (PL), time-resolved PL, and cathodoluminescence (CL) measurements. With increasing quantum well thickness up to 4 nm, the PL emission energy decreases and the blueshift of the PL emission energy increases with increasing excitation density. On the other hand, the PL emission energy of the DQWs with L_w=16 nm is higher than that of the DQWs with L_w=4 nm, and is independent of the excitation density. With increasing L_w from 1 to 4 nm, the PL decay times increase. In contrast, the decay times of 16 nm DQWs are faster than those of 4 nm DQWs. These different results for 16 nm DQWs such as the blueshift of the emission energy, the decrease of the excitation density dependence, and the increase of recombination rate can be ascribed to the relaxation of the piezoelectric field. We also observed the inhomegeneity in the CL spectra of the DQWs with L_w=1 nm on 1 μm scale.     

Prediction of soot formation characteristics in a pulverized-coal combustion field by large eddy simulations with the TDP model

In this study, the soot formation characteristics in a pulverized-coal combustion field formed by a 4 kW Central Research Institute of Electric Power Industry (CRIEPI) jet burner were predicted by large eddy simulation (LES) employing a tabulated-devolatilization-process model (TDP model) [N. Hashimoto et al., Combust. Flame 159 (2012) 353–366]. This model enables to take into account the effect of coal particle heating rate on coal pyrolysis. The coal-derived soot formation model proposed by Brown and Fletcher [A. L. Brown and T. H. Fletcher, Energy Fuels 12 (1998) 745–757] was employed in the LES. A comparison between the data predicted by LES and the soot volume fraction distribution data measured by laser induced incandescence confirmed that the soot formation characteristics in the coal combustion field of the CRIEPI burner can be accurately predicted by LES. A detailed analysis of the data predicted by LES showed that the soot particle distribution in this burner is narrow because the net soot formation rate is negative on both sides of the base of the soot volume fraction. At these positions, soot particles diffused from the peak position of soot volume fraction are oxidized due to a relatively high oxygen concentration. Finally, the effect of soot radiation on the predicted gas temperature distribution was examined by comparing the simulation results obtained with and without soot radiation. This comparison showed that the maximum gas temperature predicted by the simulation performed with soot radiation was over 100 K lower than that predicted by the simulation performed without soot radiation. From result strongly suggests the importance of considering a soot formation model for performing numerical simulations of a pulverized-coal combustion filed.     

Effect of Fe-doping on the structural,morphological and optical properties of ZnO nanoparticles synthesized by solution combustion process

The effect of Fe-doping on the structural, morphological and optical properties of ZnO nanoparticles synthesized by simple solution combustion process are reported. The powder XRD pattern indicates that the Fe-doped ZnO samples exhibit primary and secondary phases. The primary phase indicates the hexagonal wurtzite structure with the average crystalline size of around 25–50 nm and the secondary phase is associated with the face centered cubic structure of magnetite iron oxide. The elemental composition of pure and Fe-doped samples are evaluvated by EDX. The results of FE-SEM and HR-TEM cleary show that particles morphology have changed with respect to the incorporation of doping agent and particles are in aggregating nature. The vibrational properties of the synthesized ZnO nanoparticles are investigated by Raman scattering technique and it exhibits that the influence of Fe-doping significantly modify the lattice vibrational characteristics in ZnO sites. The optical properties of the Fe-doped ZnO nanoparticles are carried out by UV–vis absorption and PL spectra. The results of PL spectra show the near-band edge related emission as well as strong blue emissions in the Fe-doped ZnO nanoparticles.     

Combustion synthesis and luminescent properties of a blue-green emitting phosphor: (Ba1.95, Eu0.05)ZnSi2O7:B3+

A blue-green emitting phosphor (Ba_1.95, Eu_0.05)ZnSi₂O₇: B _x ³⁺ was prepared by combustion synthesis and an efficient blue-green emission under near-ultraviolet was observed. The luminescence, crystallinity and particle sizes were investigated by using luminescence spectrometry, X-ray diffractometry (XRD) and scanning electron microscopy (SEM). The emission spectrum shows a single band centered at 503 nm, which corresponds to the 4f ⁶5d ¹ →4f ⁷ transition of Eu²⁺. The excitation spectrum is a broad band extending from 260 to 465 nm, which matches the emission of ultraviolet light-emitting diodes. The optical absorption spectra of the (Ba_1.95, Eu_0.05)ZnSi₂O₇: B _0.06 ³⁺ exhibited band-gap energies of 3.9 eV. The results showed that boric acid was effective in improving the luminescence intensity of (Ba_1.95, Eu_0.05)ZnSi₂O₇, and the optimum molar ratio of boric acid to zinc nitrate was about 0.06. The phosphor (Ba_1.95, Eu_0.05)ZnSi₂O₇: B_0.06³⁺ synthesized by combustion method showed 1.5 times improved emission intensity compared with that of the Ba_1.95ZnSi₂O₇: Eu_0.05²⁺ phosphor under λ _ex = 353 nm.      

A spray flamelet/progress variable approach combined with a transported joint PDF model for turbulent spray flames

A spray flamelet/progress variable approach is developed for use in spray combustion with partly pre-vaporised liquid fuel, where a laminar spray flamelet library accounts for evaporation within the laminar flame structures. For this purpose, the standard spray flamelet formulation for pure evaporating liquid fuel and oxidiser is extended by a chemical reaction progress variable in both the turbulent spray flame model and the laminar spray flame structures, in order to account for the effect of pre-vaporised liquid fuel for instance through use of a pilot flame. This new approach is combined with a transported joint probability density function (PDF) method for the simulation of a turbulent piloted ethanol/air spray flame, and the extension requires the formulation of a joint three-variate PDF depending on the gas phase mixture fraction, the chemical reaction progress variable, and gas enthalpy. The molecular mixing is modelled with the extended interaction-by-exchange-with-the-mean (IEM) model, where source terms account for spray evaporation and heat exchange due to evaporation as well as the chemical reaction rate for the chemical reaction progress variable. This is the first formulation using a spray flamelet model considering both evaporation and partly pre-vaporised liquid fuel within the laminar spray flamelets. Results with this new formulation show good agreement with the experimental data provided by A.R. Masri, Sydney, Australia. The analysis of the Lagrangian statistics of the gas temperature and the OH mass fraction indicates that partially premixed combustion prevails near the nozzle exit of the spray, whereas further downstream, the non-premixed flame is promoted towards the inner rich-side of the spray jet since the pilot flame heats up the premixed inner spray zone. In summary, the simulation with the new formulation considering the reaction progress variable shows good performance, greatly improving the standard formulation, and it provides new insight into the local structure of this complex spray flame.