A statistical formulation is developed describing the composition in an evaporating multicomponent-fuel liquid drop and in the gas phase surrounding it. When a complementary discrete-component model is used, it is shown that, when drops are immersed in a carrier gas containing fuel vapor, condensation of species onto the drop results in the development of a minor peak in the liquid composition probability distribution function (PDF). This peak leads to a PDF shape that can be viewed as a combination of two gamma PDFs, which is determined by ve parameters. A model is developed for calculating the parameters of the two combined gamma PDFs. Extensive tests of the model for both diesel and gasoline show that the PDF results replicate accurately the discrete model predictions. Most important, the mean and variance of the composition at the drop surface are in excellent agreement with the discrete model. Results from the model show that although the second peak is minor for the liquid PDF, its corresponding peak for the vapor distribution at the drop surface has a comparable magnitude to and sometimes exceeds that corresponding to the rst peak. Four-parameter models are also exercised, and it is shown that they are unable to capture the physics of the problem. Superscripts.s/ = drop surface conditions .1/ = far-eld conditions

A representation of multicomponent-liquid (MC-liquid) composition as a linear combination of two single-Gamma probability distribution functions (PDFs) is used to describe a large number of MC-liquid drops evaporating in a gas flow. The PDF, called the double-Gamma PDF, depends on the molar mass. The gas-phase conservation equations are written in an Eulerian frame and the drops are described in a Lagrangian frame. Gas conservation equations for mass, momentum, species and energy are combined with differential conservation equations for the first four moments of the gascomposition PDF and coupled to the perfect gas equation of state. Source terms in all conservation equations account for the gas/drop interaction. The drop governing equations encompass differential conservation statements for position, mass, momentum, energy and four moments of the liquid-composition PDF. Simulations are performed for a three-dimensional mixing layer whose lower stream is initially laden with drops colder than the surrounding gas. Initial perturbations excite the layer to promote the double pairing of its four initial spanwise vortices to an ultimate vortex. During the layer evolution, the drops heat and evaporate. The results address the layer evolution, and the state of the gas and drops when layers reach a momentum-thickness maximum past the double vortex pairing. Of interest is the influence of the liquid composition on the development of the vortical features of the flow, on the vortical state reached after the second pairing, and on the gas temperature and composition. The MC-liquid simulations are initiated with a single-Gamma PDF composition so as to explore the development of the double-Gamma PDF. Examination of equivalent simulations with n-decane, diesel and three kerosenes as the liquid, permits assessment of the single-species versus the MC-liquid aspect, and of mixture composition specific effects. Global layer growth and global rotational characteristics are unaffected by liquid specificity; however, the global mixing is highly liquid-specific. Also liquidspecific is the evolution of the ensemble-averaged drop characteristics and of the volumetric averages representing the gas composition. Visualized rotational characteristics show that the small-scale vortical activity increases with increased fuel volatility, which is confirmed by analysis of the vorticity budgets. Homogeneous-planeaverage budgets of the vorticity and vorticity-magnitude equations indicate that the stretching and tilting, and momentum-source terms are responsible for the difference among simulations. For all MC liquids, the gas displays a high level of composition heterogeneity, which can directly be traced to the original PDF representing the MC-liquid composition. Under most conditions, the single-Gamma PDF develops into a double-Gamma PDF; however, the extent of this transformation, indicative of vapour condensation onto drops, is not readily parametrized by the liquid volatility, initial carrier-gas temperature or trace vapour in the initial ...

A model of a temporal three-dimensional mixing layer laden with fuel drops of a liquid containing a large number of species is derived. The fuel model is based on continuous thermodynamics, whereby the composition is statistically described through a distribution function parametrized on the species molar weight. The drop temperature is initially lower than that of the carrier gas, leading to drop heat up and evaporation. The model describing the changes in the multicomponent ͑MC͒ fuel drop composition and in the gas phase composition due to evaporation encompasses only two more conservation equations when compared with the equivalent single-component ͑SC͒ fuel formulation. Single drop results of a MC fuel having a sharply peaked distribution are shown to compare favorably with a validated SC-fuel drop simulation. Then, single drop comparisons are performed between results from MC fuel and a representative SC fuel used as a surrogate of the MC fuel. Further, two mixing layer simulations are conducted with a MC fuel and they are compared to representative SC-fuel simulations conducted elsewhere. Examination of the results shows that although the global layer characteristics are generally similar in the SC and MC situations, the MC layers display a higher momentum-thickness-based Reynolds number at transition. Vorticity analysis shows that the SC layers exhibit larger vortical activity than their MC counterpart. An examination of the drop organization at transition shows more structure and an increased drop-number density for MC simulations in regions of moderate and high strain. These results are primarily attributed to the slower evaporation of MC-fuel drops than of their SC counterpart. This slower evaporation is due to the lower volatility of the higher molar weight species, and also to condensation of already-evaporated species on drops that are transported in regions of different gas composition. The more volatile species released in the gas phase earlier during the drop lifetime reside in the lower stream while intermediary molar weight species, which egress after the drops are entrained in the mixing layer, reside in the mixing layer and form there a very heterogeneous mixture; the heavier species that evaporate later during the drop lifetime tend to reside in regions of high drop number density. This leads to a segregation of species in the gas phase based on the relative evaporation time from the drops. The ensemble-average drop temperature becomes eventually larger/smaller than the initial drop temperature in MC/SC simulations. Neither this species segregation nor the drop temperature variation with respect to the initial temperature or as a function of the mass loading can be captured by the SC-fuel simulations.

No abstract

An extension of the Turbulent Flame Speed Closure model rendering the model applicable to multiphase flow and ignition is presented. As formerly no coupling between reaction progress variable and enthalpy was existent, except through the temperature dependency of the laminar flame speed, an adaptation is proposed which offers an interface to initiate the combustion process. The modification to incorporate multiphase conditions is achieved by substituting the mixture fraction variable as representation of the composition in the original implementation of the Turbulent Flame Speed Closure model with independent species. Source terms to correlate the species progress to the reaction progress variable are derived in this work. The additional transport equations serve a higher generality of the model and enable the proper treatment of vaporizing fuel droplets. It is demonstrated that limitations which arise in the standard formulation of the model, stemming from differences in the transport equation for the reaction progress variable and the mixture fraction, are addressed and resolved by the new approach. Regarding the initiation of the flame, an additional source term for the reaction progress variable is introduced, which relates the reaction progress to the auto-ignition time. This allows the development of the flame without imposing artificial boundary conditions. The correct model behavior is established by means of a series of widely used test cases. The results of these simulations show that the model's potential to predict flame growth and more generally the flame evolution as a function of time and space is preserved. At the same time more sophisticated test case boundary conditions involving multiphase conditions and variable inflows in terms of composition can be incorporated. As a thorough assessment of the extended model capabilities, a multiphase lab scale setup , which provides a comprehensive data set, is presented. The good agreement of the obtained results underline the range of applicability of the extended model and its accuracy, albeit its simplicity, for multiphase conditions.