Analysis of Mixing Models for use in Simulations of Turbulent Spray Combustion

The present study concerns the investigation of different mixing models for use in the transported probability density function (PDF) modeling of turbulent (reacting) spray flows. The modeling of the turbulent mixing and other characteristic scalar variables such as gas enthalpy using transported (joint) PDFs has become an important method to describe turbulent (reacting) spray flows since the evaporation process causes the PDF of the mixture fraction to deviate from the widely used β function, which is typically used in models for turbulent gas flows. In the PDF transport equation, the molecular mixing does not appear in closed form so that modeling strategies are required. For gas combustion, the interaction-by-exchange-with-the-mean (IEM) model, the modified Curl (MC) model, and the Euclidean minimum spanning tree (EMST) models are used. More recently, a new mixing model, the PSP model, which is based on parameterized scalar profiles has been developed. The present study focuses on the use and analysis of the IEM, MC and PSP models for turbulent spray flames. For this purpose, the models are reconsidered with respect to the evaporation process that must be included and evaluated if spray combustion is considered. For model evaluation, turbulent ethanol/air spray flames are simulated, and the results are compared to experimental data by A. Masri, University of Sydney, Australia.     

Effects of Turbulence, Evaporation and Heat Release on the Dispersion of Droplets in Dilute Spray Jets and Flames

The dispersion characteristics of a selection of non-evaporating non-reacting, evaporating non-reacting, and reacting dilute spray jets issuing in ambient air (Gounder et al, Combust Sci Technol 182:702–715, 2010; Masri and Gounder, Combust Flame 159:3372–3397, 2010) and in a hot coflow (Oloughlin and Masri, Flow Turbul Combust 89:13–35, 2012) are analysed. Other than the cases found in those contributions, two additional sprays of kerosene have been investigated in order to systematically study the effects of evaporation. The burners are well designed such that boundary conditions may be accurately measured for use in numerical simulations. The dynamics and dispersion characteristics are analysed by conditioning results on the droplet Stokes numbers and by systematically investigating changes in dispersion and dynamics as a function of carrier air velocity, liquid loading, ignition method, and location within the flame or spray jet. The tendency for droplet dispersion defined by the ratio of radial rms velocity to axial mean velocity varies significantly between reacting and non-reacting flows. However, dispersion is found to be largely unaffected by evaporation. The total particle concentration, or number density of droplets within the spray has also been used as a direct measure of spray dispersion with the effect of evaporation on a turbulent polydisperse spray being isolated by investigating acetone and kerosene sprays with similar boundary conditions. The rate of change of droplet size with radial position is almost identical for the kerosene and acetone cases. The dispersion characteristics, closely related to the ‘fan spreading’ phenomenon are dependant on the carrier air velocity and axial location within the spray.     

Modelling and simulation of sprays in laminar flames

In the present paper we model and numerically simulate steady, laminar, premixed spray flames. The gasphase is described in Eulerian form by the equations governing the conservation of overall mass, momentum, energy and species mass. The liquid phase is described in Lagrangian form by the overall continuity equation, which reduces to an equation for the droplet size, the equations of motion, the energy equation and a droplet density function transport equation. The latter is the so-called spray equation, which, at any position in the chemically reacting flowfield, describes the joint distribution of droplet size, droplet velocity and droplet temperature. Herein the spray equation is solved using a Monte Carlo method. Detailed models of the exchange of mass, momentum and energy between the gaseous and the liquid phase are taken into account. The results presented in this paper are for an octane-air flame, where small amounts of liquid octane in form of a liquid spray are added to a fresh, unburnt gaseous octane-air mixture.Presented at Euromech Colloquium 324: The Combustion of Drops, Sprays and Aerosols, 25th–27th July 1994, Marseilles, France.     

A Comprehensive Study of Effects of Mixing and Chemical Kinetic Models on Predictions of n-heptane Jet Ignitions with the PDF Method

The composition Probability Density Function (PDF) model is coupled with a Reynolds-averaged k???ε turbulence model and three computationally efficient, yet widely used chemical mechanisms to simulate transient n-heptane spray injection and ignition in a high temperature and high density ambient fluid. Molecular diffusion is modelled by three mixing models, namely the interaction by exchange with the mean (IEM), modified Curl (MC) and Euclidean minimum spanning trees (EMST) models. The liquid phase is modelled by a discrete phase model (DPM). This represents among the first applications of the PDF method in practical diesel engine conditions. A non-reacting case is first considered, with the focus on the ability of the model to capture the spray structure, e.g., vapour penetration and liquid length, fuel mixture fraction and its variance. Reacting cases are then investigated to compare and evaluate the three different chemical mechanisms and the three mixing models. It is concluded that the EMST mixing model in conjunction with a reduced chemical kinetic model (Lu et al., Combust Flame 156(8):1542–1551, 2009) performs the best among the options considered. The sensitivity of the results to the choice of the mixing constant is also studied to understand its effect on the flame ignition and stabilisation. Finally, the PDF model is compared to a well-mixed model that assumes turbulent fluctuations are negligible, which has been widely used in the diesel spray combustion community. Significant structural differences in the modelled flame are revealed comparing the PDF method with the well-mixed model. Quantitatively, the PDF model exhibits excellent agreement with the measurements and shows much better results than the well-mixed model in all ambient O₂ and temperature conditions.     

Simulation of Dilute Acetone Spray Flames with LES-CMC Using Two Conditional Moments

Large-eddy simulations (LES) have been coupled with a conditional moment closure (CMC) method for the computation of a series of turbulent spray flames. An earlier study by Ukai et al. (Proc. Combust. Inst. 34(1),1643–1650, 2013) gave reasonable results for the prediction of temperature and velocity profiles, but some limitations of the method became apparent. These limitations are primarily related to the upper limit in mixture fraction space. In order to enhance the applicability of the LES-CMC model, this paper proposes a two-conditional moment approach to account for the existence of pre-evaporated fuel by introducing two sets of conditional moments based on different mixture fractions. The two-conditional moment approach is first tested for a non-reacting test case. The results indicate that the spray evaporation induces relatively large conditional fluctuations within a CMC cell, and one set of conditional moments might not be sufficient. The upper limit of the mixture fraction space is dynamically selected for the solution of the second set of conditional moments, and the corresponding CMC solution in a CFD cell is estimated by interpolation between the two conditional moments weighted by the amount of vapour emitted within the domain. The cell-filtered value is given by integration of the conditional moment across mixture fraction space using a bounded β-FDF for the distribution of the scalar. As a result, the fuel concentration profiles given by LES and the two-conditional moment approach are shown to agree well. Then, the two-conditional moment approach is applied to four different flame configurations. The comparison of LES cell quantities and conditionally averaged moments indicates that the two sets of conditional moments are necessary for accurate predictions in zones where gas phase mixture fraction is significantly increased by droplet evaporation within the computational domain. The unconditional temperature profiles clearly show that the new approach improves the predictions of mean temperature especially along the centerline. Also, the better predictions of the temperature field improve the accuracy of the predicted mean axial droplet velocities. Overall, good agreement with the experimental results is found for all four cases, and the methodology is shown to be applicable to flames with a relatively wide range of fuel vapour concentrations.     

Modeling of Turbulent Natural Gas and Biogas Flames of the Delft Jet-in-Hot-Coflow Burner: Effects of Coflow Temperature,Fuel Temperature and Fuel Composition on the Flame Lift-Off Height

The structure of a turbulent non-premixed flame of a biogas fuel in a hot and diluted coflow mimicking moderate and intense low dilution (MILD) combustion is studied numerically. Biogas fuel is obtained by dilution of Dutch natural gas (DNG) with CO₂. The results of biogas combustion are compared with those of DNG combustion in the Delft Jet-in-Hot-Coflow (DJHC) burner. New experimental measurements of lift-off height and of velocity and temperature statistics have been made to provide a database for evaluating the capability of numerical methods in predicting the flame structure. Compared to the lift-off height of the DNG flame, addition of 30 % carbon dioxide to the fuel increases the lift-off height by less than 15 %. Numerical simulations are conducted by solving the RANS equations using Reynolds stress model (RSM) as turbulence model in combination with EDC (Eddy Dissipation Concept) and transported probability density function (PDF) as turbulence-chemistry interaction models. The DRM19 reduced mechanism is used as chemical kinetics with the EDC model. A tabulated chemistry model based on the Flamelet Generated Manifold (FGM) is adopted in the PDF method. The table describes a non-adiabatic three stream mixing problem between fuel, coflow and ambient air based on igniting counterflow diffusion flamelets. The results show that the EDC/DRM19 and PDF/FGM models predict the experimentally observed decreasing trend of lift-off height with increase of the coflow temperature. Although more detailed chemistry is used with EDC, the temperature fluctuations at the coflow inlet (approximately 100K) cannot be included resulting in a significant overprediction of the flame temperature. Only the PDF modeling results with temperature fluctuations predict the correct mean temperature profiles of the biogas case and compare well with the experimental temperature distributions.     

Statistical Analysis of Turbulent Flame-Droplet Interaction: A Direct Numerical Simulation Study

Turbulent combustion of mono-disperse droplet-mist has been analysed based on three-dimensional Direct Numerical Simulations (DNS) in canonical configuration under decaying turbulence for a range of different values of droplet equivalence ratio (?_d), droplet diameter (a_d) and root-mean-square value of turbulent velocity (u^′). The fuel is supplied in liquid phase and the evaporation of droplets gives rise to gaseous fuel for the flame propagation into the droplet-mist. It has been found that initial droplet diameter, turbulence intensity and droplet equivalence ratio can have significant influences on the volume-integrated burning rate, flame surface area and burning rate per unit area. The droplets are found to evaporate predominantly in the preheat zone, but some droplets penetrate the flame front, reaching the burned gas side where they evaporate and some of the resulting fuel vapour diffuses back towards the flame front. The combustion process in gaseous phase takes place predominantly in fuel-lean mode even for ?_d > 1. The probability of finding fuel-lean mixture increases with increasing initial droplet diameter because of slower evaporation of larger droplets and this predominantly fuel-lean mode of combustion exhibits the attributes of low Damköhler number combustion and gives rise to thickening of flame with increasing droplet diameter. The chemical reaction is found to take place under both premixed and non-premixed modes of combustion and the relative contribution of non-premixed combustion to overall heat release increases with increasing droplet size. The statistical behaviours of the flame propagation and mode of combustion have been analysed in detail and detailed physical explanations have been provided for the observed behaviour.     

Acetone Droplet Behavior in Reacting and Non Reacting Turbulent Flow

Acetone droplet characteristics in reacting and non-reacting turbulent flow are predicted and compared to experimental data. Investigations are conducted to study the effects of surrounding environment properties on the velocities, dispersion, and evaporation of a relatively volatile spray fuel that featured a wide range of Stokes numbers. The simulations are performed in the framework of Reynolds Averaged Navier Stokes equations along with the Eulerian-Lagrangian approach in which 12 different classes of the dispersed phase. The phase transition is modeled by the Langmuir-Knudsen law that accounts for non equilibrium effects based on a consistent determination of the molar mass fraction on the droplet surfaces. For the droplet dispersion, the Markov sequence model is improved by adding a correction drift term to the fluid fluctuation velocity at the parcel position along the droplet trajectory. This correction term aimed at accounting for the non-homogeneity effects in the turbulent flow. The combustion is captured using the Bray-Moss-Libby model that is extended to account for the partially premixed spray combustion. The chemistry is described with the flamelet model using a recent detailed reaction mechanism that involves 84 species and 409 reactions for which the Lewis number is not set to the unity. Mean droplet velocities for reacting and non-reacting test cases are compared with experimental data. Good agreement is observed. The spray is interacting with the nozzle edge developing new classes and relatively dense region. Hence the RMS-velocities close to the nozzle exit plan demonstrate discrepancies. The droplets group combustion effect is found to be important in the modeling of the burning velocity which influences the flame propagation. Reasonable agreements between the numerical and the experimental results are also observed in the spray flux and temperature profiles.     

Droplet motion in an electrohydrodynamic fine spray

The results of a combined experimental and numerical study on droplet behavior within an electrohydrodynamic fine spray are presented. The fine spray exists in the transition region between the multiple cone-jet and rim emission spray modes. Experiments were conducted specifically to characterize the motion of droplets within the spray. Light-sheet visualizations and measurements of droplet speed and velocity using laser-based, single-particle counters were obtained. Additionally, a numerical simulation of the droplet motion within the spray was made and compared to the experimental results. The electrohydrodynamic fine spray of ethanol droplets ( 1 to 40 m diameter) was generated using a typical capillary-plate configuration, with a capillary tip electric field intensity of 10⁶ V/m and a spray charge density of 70 C/m³. Acquired images of the spray revealed a zone of rapid expansion near the capillary followed by a more gradual expansion farther from the capillary. In situ laser-diagnostic measurements confirmed these observations. Measured droplet speeds decreased rapidly with increasing axial distance from the capillary, but then increased beyond the spray's axial mid-plane as a result of a change in the sign of the axial internal electric field. Droplet axial velocity components behaved similarly. The radial velocity components exhibited a maximum value off of the spray's centerline in the near-capillary region. Farther away from the capillary, they increased monotonically with increasing radial position. These trends identified the significant role that the radial internal electric field plays in spray expansion. The numerical simulation of the normal spray verified the inferred change in sign of the axial internal field and underscored the dominant contribution of the external electric field in the near-capillary region and of the internal electric field farther away.     

A Probability Density Function Method for Modelling Liquid Fuel Sprays

A stochastic model is proposed for modelling the turbulent dispersion of liquid fuel sprays. The approach adopted is based on the evolution equation for the joint probability density function (PDF) of the droplet properties. Turbulent dispersion is described by Langevin's equation in which a Weiner process is used to represent the stochastic force term. The method leads to plausible results when applied to a kerosene spray flame and provides a rational framework for the incorporation of liquid film break up and droplet formation processes. This revised version was published online in July 2006 with corrections to the Cover Date.     

Effects of Spray and Turbulence Modelling on the Mixing and Combustion Characteristics of an n-heptane Spray Flame Simulated with Dynamic Adaptive Chemistry

Accurate modelling of spray combustion process is essential for efficiency improvement and emissions reduction in practical combustion engines. In this work, both unsteady Reynolds-averaged Navier-Stokes (URANS) simulations and large eddy simulations (LES) are performed to investigate the effects of spray and turbulence modelling on the mixing and combustion characteristics of an n-heptane spray flame in a constant volume chamber at realistic conditions. The non-reacting spray process is first simulated with URANS to investigate the effects of entrainment gas-jet model on the penetration characteristics and fuel vapor distributions. It is found that the droplet motion near the nozzle has significant influence on the fuel vapor distribution, while the liquid penetration length is controlled by the evaporation process and insensitive to gas-jet model. For the case considered, both URANS with the gas-jet model and large eddy simulations can properly predict the vapor penetration. For the combustion characteristics, it is found that LES yields better predictions in the global combustion characteristics. The URANS with gas jet model yields a comparable flame length and lift-off-length (LOL) to LES, but results in a larger ignition delay time compared to the experimental data. Another focus of this work is to qualify the convergence characteristics of the dynamic adaptive chemistry (DAC) method in these transient combustion simulations, where DAC is applied to reduce the mechanism locally and on-the-fly to accelerate chemistry calculations. The instantaneous flame structures and global combustion characteristics such as ignition delay time, flame lift-off length and emissions are compared between simulations with and without DAC. For URANS, good agreements are observed both on instantaneous flame structures and global characteristics. For LES, it is shown that the errors incurred by DAC are small for scatter distributions in composition space and global combustion characteristics, while they may significantly affect instantaneous flame structures in physical space. The study reveals that for DAC application in transient simulations, global or statistic information should be used to assess the accuracy, such as manifolds in composition space, conditional quantities and global combustion characteristics. For the cases investigated, a speed-up factor of more than two is achieved by DAC with a 92-species skeletal mechanism with less than 0.2 % and 3.0 % discrepancy in ignition delay and LOL, respectively.     

Large Eddy Simulation of a Polydisperse Ethanol Spray Flame

Ethanol is identified as an interesting alternative fuel. In this regards, the predictive capability of combustion Large Eddy Simulation approach coupled to Lagrangian droplet dynamic model to retrieve the turbulent droplet dispersion, droplet size distribution, spray evolution and combustion properties is investigated in this paper for an ethanol spray flame. Following the Eulerian-Lagrangian approach with a fully two way coupling, the Favre-filtered low Mach number Navier-Stokes equations are solved on structured grids with dynamic sub-grid scale models to describe the turbulent carrier gas phase. Droplets are injected in polydisperse manner and generated in time dependent boundary conditions. They evaporate to form an air-fuel mixture that yields spray flame. Part of the ethanol droplets evaporates within the prevaporization area before reaching the combustion zone, making the flame to burn in a partially premixed regime. The chemistry is described by a tabulated detailed chemistry based on the flamelet generated manifold approach. The fuel, ethanol, is modeled by a detailed reaction mechanism consisting of 56 species and 351 reversible reactions. The simulation results including excess gas temperature, droplet velocities and corresponding fluctuations, droplet mean diameters and spray volume flux at different distances from the exit plane show good agreement with experimental data. Analysis of combustion spray features allows gaining a deep insight into the two-phase flow process ongoing.     

Joint Scalar versus Joint Velocity-Scalar PDF Simulations of Bluff-Body Stabilized Flames with REDIM

Two transported PDF strategies, joint velocity-scalar PDF (JVSPDF) and joint scalar PDF (JSPDF), are investigated for bluff-body stabilized jet-type turbulent diffusion flames with a variable degree of turbulence–chemistry interaction. Chemistry is modeled by means of the novel reaction-diffusion manifold (REDIM) technique. A detailed chemistry mechanism is reduced, including diffusion effects, with N ₂ and CO ₂ mass fractions as reduced coordinates. The second-moment closure RANS turbulence model and the modified Curl’s micro-mixing model are not varied. Radiative heat loss effects are ignored. The results for mean velocity and velocity fluctuations in physical space are very similar for both PDF methods. They agree well with experimental data up to the neck zone. Each of the two PDF approaches implies a different closure for the velocity-scalar correlation. This leads to differences in the radial profiles in physical space of mean scalars and mixture fraction variance, due to different scalar flux modeling. Differences are visible in mean mixture fraction and mean temperature, as well as in mixture fraction variance. In principle, the JVSPDF simulations can be closer to physical reality, as a differential model is implied for the scalar fluxes, whereas the gradient diffusion hypothesis is implied in JSPDF simulations. Yet, in JSPDF simulations, turbulent diffusion can be tuned by means of the turbulent Schmidt number. In the neck zone, where the turbulent flow field results deteriorate, the joint scalar PDF results are in somewhat better agreement with experimental data, for the test cases considered. In composition space, where results are reported as scatter plots, differences between the two PDF strategies are small in the calculations at hand, with a little more local extinction in the joint scalar PDF results.     

The influence of gas phase velocity fluctuations on primary atomization and droplet deformation

The effects of grid-generated velocity fluctuations on the primary atomization and subsequent droplet deformation of a range of laminar liquid jets are examined using microscopic high-speed backlit imaging of the break-up zone and laser Doppler anemometry of the gas phase separately. This is done for fixed gas mean flow conditions in a miniature wind tunnel experiment utilizing a selection of fuels, turbulence-generating grids and two syringe sizes. The constant mean flow allows for an isolated study of velocity fluctuation effects on primary atomization in a close approximation to homogeneous decaying turbulence. The qualitative morphology of the primary break-up region is examined over a range of turbulence intensities, and spectral analysis is performed in order to ascertain the break-up frequency which, for a case of no grid, compares well with the existing literature. The addition of velocity fluctuations tends to randomize the break-up process. Slightly downstream of the break-up region, image processing is conducted in order to extract a number of metrics, which do not depend on droplet sphericity, and these include droplet aspect ratio and orientation, the latter quantity being somewhat unconventional in spray characterization. A turbulent Weber number $We^{\prime}$ which takes into account gas phase fluctuations is utilized to characterize the resulting droplet shapes, in addition to a mean Weber number <We _d>. Above a $We^{\prime}>0.05$ a clear positive relationship exists between the mean aspect ratio of droplets and the turbulent Weber number where $We^{\prime}$ is varied by altering all relevant variables including the velocity root mean square, the initial droplet diameter, the surface tension and the density.     

Joint PDF Closure of Turbulent Premixed Flames

In this paper, a novel model for turbulent premixed combustion in the corrugated flamelet regime is presented, which is based on transporting a joint probability density function (PDF) of velocity, turbulence frequency and a scalar vector. Due to the high dimensionality of the corresponding sample space, the PDF equation is solved with a Monte-Carlo method, where individual fluid elements are represented by computational particles. Unlike in most other PDF methods, the source term not only describes reaction rates, but accounts for “ignition” of reactive unburnt fluid elements due to propagating embedded quasi laminar flames within a turbulent flame brush. Unperturbed embedded flame structures and a constant laminar flame speed (as expected in the corrugated flamelet regime) are assumed. The probability for an individual particle to “ignite” during a time step is calculated based on an estimate of the mean flame surface density (FSD), latter gets transported by the PDF method. Whereas this model concept has recently been published [21], here, a new model to account for local production and dissipation of the FSD is proposed. The following particle properties are introduced: a flag indicating whether a particle represents the unburnt mixture; a flame residence time, which allows to resolve the embedded quasi laminar flame structure; and a flag indicating whether the flame residence time lies within a specified range. Latter is used to transport the FSD, but to account for flame stretching, curvature effects, collapse and cusp formation, a mixing model for the residence time is employed. The same mixing model also accounts for molecular mixing of the products with a co-flow. To validate the proposed PDF model, simulation results of three piloted methane-air Bunsen flames are compared with experimental data and very good agreement is observed.     

Quadrature method of moments for modeling multi-component spray vaporization

This article puts forward the quadrature method of moments (QMoM) for modeling droplet composition during the spray vaporization process. This method is implemented for solving the Continuous Thermodynamic Model (CTM) of multi-component droplet vaporization, an advantageous alternative to the classical Discrete Component Model (DCM) when the droplet is formed of a great number of components. The CTM approach consists in modeling the droplet’s composition using a probability density function (PDF). This method was first tried out for vaporizing droplets by Hallett, who assumed a Gamma-function for the PDF. However, Harstadt et al. underlined some problems in the case of vapor condensation on the droplet surface, since the Gamma-PDF model presumes the PDF’s mathematical form. The QMoM which does not require this hypothesis is studied in this article, according to Lage’s research dealing with QMoM application to phase equilibria. The numerical features of QMoM are investigated in detail, and then the method is implemented for the difficult test case of vapor condensation. The results are analyzed to illustrate the application of QMoM to multi-component droplet vaporization modeling and to provide a better understanding of the QMoM main advantages and limitations.     

Experimental and numerical analysis of spray-atomization characteristics of biodiesel fuel in various fuel and ambient temperatures conditions

The purpose of this work is to reveal the effects of fuel temperatures and ambient gas conditions on the spray-atomization behavior of soybean oil methyl ester (SME) fuel. The spray-atomization behavior was analyzed through spray parameters such as the axial distance from the nozzle tip, local and overall Sauter mean diameter (SMD). These parameters were obtained from a spray visualization system and a droplet measuring system. In addition, the experimental results were compared with the numerical results calculated by the KIVA-3V code. It was revealed that the increase of the fuel temperature (from 300 K to 360 K) little affects the spray liquid tip penetration. The increase of the ambient gas temperature (from 300 K to 450 K) caused a increase in the spray liquid tip penetration. Also, biodiesel fuel evaporation actively occurred due to the increase in the fuel temperature and the ambient gas temperature. Of special significance was that the highest vapor fuel mass concentration was observed at the center region of the spray axis. In the results of the microscopic characteristics, the detected local droplet size at the axial direction and overall droplet size at the axial and radial direction in a control volume increased when the fuel temperature increased. This is believed to be due to an increase in the number of small droplets that quickly evaporated. In addition, the increased fuel temperature caused the decrease of the number of droplets and the increase of the vapor fuel mass. The mean axial velocity of droplets decreased with increasing fuel temperature.     

Accounting for Polydispersion in the Eulerian Large Eddy Simulation of the Two-Phase Flow in an Aeronautical-type Burner

A major issue for the simulation of two-phase flows in engines concerns the modeling of the liquid disperse phase, either in the Lagrangian or the Eulerian approach. In the perspective of massively parallel computing, the Eulerian approach seems to be more suitable, as it uses the same algorithms as the gaseous phase solver. However taking into account the whole physics of a turbulent spray, especially in terms of polydispersity, requires an additional modeling effort. The Mesoscopic Eulerian Formalism (MEF) [13] accounts for the effect of turbulence on the disperse phase, and was extended to the Large Eddy Simulation framework [41], but is limited to monodisperse flows. In [38], the influence of polydispersity on resolved and unresolved turbulent motions of the disperse phase was highlighted, and a first model was proposed, based on size-conditioned statistics. Starting from this idea, a coupling between the MEF and the Multifluid Approach (MA) [30] is proposed. The MA decomposes the Eulerian phase into several fluid classes called sections, and corresponding to size intervals. Each section uses then size-conditioned closures. The original idea of this work is to use the MEF closures independently in each section, taking into account the mean droplet size of this section. This new approach, called Multifluid Mesoscopic Eulerian Formalism (MMEF), is then able to capture polydispersion with associated size-conditioned turbulent dynamics. First, the importance of polydispersity and the ability of MMEF to capture it are highlighted with a 0D evaporation case and a 2D vortex case, showing its impact on dynamics in both size and physical spaces. Then, the MMEF is applied to the MERCATO configuration of ONERA [18]. Results are compared to monodisperse Eulerian, Lagrangian and experimental results.     

Transported PDF Modelling of a High Velocity Bluff-Body Stabilised Flame (HM2) Using Detailed Chemistry

A transported probability density function (PDF) approach closed at the joint scalar level was used to model a bluff body stabilised turbulent diffusion flame (HM2) investigated experimentally by Masri and co-workers. The current effort extends a previous study of HM1 (Re?=?15,800) to a flame with a higher degree of local extinction (Re?=?23,900). The impact of an algebraic model that accounts for local Damköhler number effects on the time-scale ratio of scalar to mechanical turbulence is also evaluated along with the impact of improved thermochemistry. The computations have been performed using a hybrid Monte Carlo/finite volume algorithm and a systematically reduced H/C/N/O mechanism featuring 300 reactions, 20 solved and 28 steady-state species. The joint scalar PDF equations were solved using moving particles in a Lagrangian framework and the velocity field was closed at the second moment level. The redistribution terms were modelled using the Generalized Langevin model of Haworth and Pope. Results show that scalar fields are reproduced with encouraging accuracy and that the revised time scale model improves agreement with experimental data. A high sensitivity to the NO chemistry was observed and encouraging agreement was obtained for the first two moments following adoption of updated reaction rates proposed in an earlier study.     

Transient analysis of sub-critical evaporation of fuel droplet in non-isothermal stagnant gaseous mixtures: effects of radiation and thermal expansion

Transient sub-critical droplet evaporation in non-isothermal stagnant gaseous mixtures taking into account the effects of radiation, liquid volumetric expansion and droplet heating is investigated numerically. We obtained equations for Stefan velocity and the rate of change of the droplet radius taking into account liquid volumetric expansion, and derived the boundary conditions taking into account the effect of liquid thermal expansion. It is shown that in the case of sub-critical evaporation neglecting the liquid volumetric expansion causes underestimation of the evaporation rate at the initial stage and overestimation of the evaporation rate at the final stage of droplet evaporation.