Journal of Combustion

Journal of Combustion / 2012 / Article

Research Article | Open Access

Volume 2012 |Article ID 138619 |

Etim S. Udoetok, "Modified Quasi-Steady Fuel Droplet Combustion Model", Journal of Combustion, vol. 2012, Article ID 138619, 8 pages, 2012.

Modified Quasi-Steady Fuel Droplet Combustion Model

Academic Editor: Evangelos G. Giakoumis
Received19 Dec 2011
Accepted10 Mar 2012
Published17 May 2012


The quasi-steady model of the combustion of a fuel droplet has been modified. The approach involved the modification of the quasi-steady model to reflect the difference in constant properties across the flame front. New methods for accurately estimating gas constants and for estimating Lewis number are presented. The proposed theoretical model provides results that correlate favorably with published experimental results. The proposed theoretical model also eliminates the need for unguided adjustment of thermal constants or the complex analysis of the variation of thermal properties with temperature and can serve as a basis for analysis of other combustion conditions like droplets cloud and convective and high-pressure conditions.

1. Introduction

Fuel droplet models are used to describe the influence of droplet size and ambient conditions on fuel combustion in devices such as diesel engines, rocket engines, gas turbines, oil fired boilers, and furnaces [1, 2]. The simple quasi-steady model, which is the focus of this paper, has its origin in the 1950s [3, 4] and is widely accepted as the theoretical model of fuel droplet combustion [5]. The important results of this theory are as follows [1, 4ā€“7]:Ģ‡š‘š=4šœ‹š‘˜š‘”š‘Ÿš‘ š‘š‘š‘”ln(1+šµ),(1) where šµ=Ī”ā„Žš‘/šœ+š‘š‘š‘”ī€·š‘‡āˆžāˆ’š‘‡š‘ ī€øš‘žš‘–āˆ’š‘™+ā„Žš‘“š‘”,š‘‡š‘“=ī€·š‘žš‘–āˆ’š‘™+ā„Žš‘“š‘”ī€øš‘š‘š‘”([]1+šœ)šœšµāˆ’1+š‘‡š‘ ,š‘Ÿš‘“š‘Ÿš‘ =ln(1+šµ),š·ln(1+1/šœ)2=š·02š‘”āˆ’š¾š‘”,š‘‘=š·02š¾,(2) where š¾=8š‘˜š‘”šœŒš‘™š‘š‘š‘”ln(1+šµ)=2Ģ‡š‘ššœ‹šœŒš‘™š‘Ÿš‘ .(3) Note that though the simple quasi-steady model theory has its roots in the 1950s, it has been developed upon through the 1950s to the 1970s in order to arrive at the results presented above [1ā€“7]. However, simple quasi-steady model provides unsatisfactory results in comparison with several experimental observations [1, 5, 8ā€“12]. The simple quasi-steady model is best at predicting fuel mass flow rate Ģ‡š‘š and the š‘‘2 law but does not accurately predict the flame to droplet radius ratio š‘Ÿš‘“/š‘Ÿš‘  and flame temperature š‘‡š‘“ [5]. The simple quasi-steady model typically predicts š‘Ÿš‘“/š‘Ÿš‘  higher than observed values, and š‘‡š‘“ is usually lower than that experimentally observed. Experiments have also shown that the values of š‘Ÿš‘“/š‘Ÿš‘  and š‘‡š‘“ may not be constant. Law et al. [10] showed that fuel vapor accumulation causes transient effects in the values of š‘Ÿš‘“/š‘Ÿš‘  and š‘‡š‘“, but Raghunandan and Mukunda [5] later showed that the condensed-phase unsteadiness lasts for 20ā€“25% of burning time and concluded that the discrepancies between experiments and simple quasi-steady model cannot be attributed to condensed-phase unsteadiness. A study by He et al. [8] revealed that the flame front motion has the effect of causing unsteadiness and variation of š‘Ÿš‘“/š‘Ÿš‘  and š‘‡š‘“ during combustion. The analytical model of Raghunandan and Mukunda [5] for quasi-steady droplet combustion with variable thermodynamic and transport properties, and nonunity Lewis number gave accurate prediction of š¾ and significant improvement in š‘Ÿš‘“/š‘Ÿš‘  and š‘‡š‘“ compared to the simple quasi-steady model. Puri and Libby [13] used a detailed expression for the heat transfer and transport properties and came up with a complex model of the quasi-steady fuel droplet combustion. Their model is best solved numerically. Filho [14] solved the quasi-steady fuel droplet combustion problem in a way similar to Puri and Libbyā€™s solution. Filhoā€™s solution was less complicated and involved the removal of nonlinearity in the heat transfer and transport properties coefficients. Imaoka and Sirignano [15ā€“17] solved the fuel droplet combustion problem for the case of a droplet cloud using unity Lewis number assumption. They acknowledged that the use of unity Lewis number significantly overestimates š‘Ÿš‘“/š‘Ÿš‘ . Imaoka and Sirignano [15ā€“17] focused on the variation of constants from one droplet to another in the droplet cloud. The vaporization rate was found for each droplet in the droplet cloud because they assumed that the solutions for each droplet are not equal.

In this paper, the simple quasi-steady fuel droplet combustion model is modified for higher accuracy by assuming discontinuity in the heat transfer and transport coefficients across the flame sheet and nonunity Lewis number for the inner and outer region. A method for estimating property constants for the two regions is recommended. Note that while the discontinuity in Imaoka and Sirignano [15ā€“17] solution is from droplet to droplet, the discontinuity assumed in this paper is from the inner region to the outer region of a burning droplet.

2. Method

In the derivation of the classical droplet combustion model, the following assumptions are made [1, 4ā€“7].(i) Burning droplet is spherical and surrounded by a spherically symmetric flame in a quiescent infinite medium.(ii)Burning process is quasi-steady.(iii)Fuel is a single component and pressure is uniform and constant.(iv)Gaseous species are of 3 types: fuel vapor, oxidizer, and combustion products.(v)Stoichiometric proportions of fuel-oxidizer are at flame.(vi) Unity Lewis number is assumed.(vii)Radiation heat transfer is negligible.(viii)No soot or liquid water is present.(ix)Uniform species thermal constants: š‘š‘ and š‘˜.

These assumptions are good, but the following assumptions changes will be made in order to improve the accuracy of the model.(i)Unity Lewis is assumed only at the source of diffusing species, and nonunity Lewis number is assumed in the outer and inner regions. This assumption is made because, at the sources of diffusing specie, the generation of the diffusing specie causes the thermal diffusivity to balance the mass diffusivity, while, away from the source of diffusing species, the thermal diffusivity and mass diffusivity have different values depending on the species concentration, species properties, and temperature profile. (ii)The property of the inner region is different from the property of the outer region. This assumption is made because the average temperature in the outer region is different from the average temperature in the inner region and the species composition in the outer region is different from the species composition in the inner region.

The new assumption that the property of the inner region is different from the property of the outer region is shown in Figure 1. This new assumption makes it necessary to have two average temperatures since there are two different sets of temperature extremes in the two regions.

In order to relate the outer constants to the inner constants, let š‘+=š‘›š‘āˆ’,(4) where š‘› is a constant and (4) implies that š‘›=šæš‘’+š‘+š‘š‘”š‘˜āˆ’š‘”šæš‘’āˆ’š‘āˆ’š‘š‘”š‘˜+š‘”.(5) The mass flow rate, Ģ‡š‘š, is treated as a constant and independent of radius, š‘Ÿ, since quasi-steady burning is assumed. In the inner region, Fickā€™s law can be presented in the form [1] Ģ‡š‘šāˆ£fuel=āˆ’4šœ‹š‘Ÿ2šœŒD1āˆ’š‘Œš¹š‘‘š‘Œš¹š‘‘š‘Ÿ(6) with boundary conditions (BCs) š‘Œš¹ī€·š‘Ÿš‘ ī€ø=š‘Œš¹,š‘ ī€·š‘‡š‘ ī€ø,š‘Œš¹ī€·š‘Ÿš‘“ī€ø=0.(7) Integration of (6) and application of BCs (7) gives š‘Œš¹,š‘ ī€·=1āˆ’expāˆ’š‘āˆ’Ģ‡š‘š/š‘Ÿš‘ ī€øī€·expāˆ’š‘āˆ’Ģ‡š‘š/š‘Ÿš‘“ī€ø,(8) whereš‘=1/4šœ‹šœŒD=šæš‘’ā‹…š‘š‘š‘”/4šœ‹š‘˜š‘”.

In the outer region, Fickā€™s law in terms of the constant fuel mass flow rate can be presented in the form [1] Ģ‡š‘šāˆ£fuel=āˆ’4šœ‹š‘Ÿ2šœŒDšœāˆ’š‘Œš‘‚š‘„š‘‘š‘Œš‘‚š‘„š‘‘š‘Ÿ(9) with BCs š‘Œš‘‚š‘„ī€·š‘Ÿš‘“ī€øš‘Œ=0,š‘‚š‘„(š‘ŸāŸ¶āˆž)=1.(10) Integration of (9) and application of BCs (10) gives the relation between Ģ‡š‘š and š‘Ÿš‘“ as ī‚µš‘exp+Ģ‡š‘šš‘Ÿš‘“ī‚¶=šœ+1šœ.(11) Equation (4) in (11) gives ī‚µexpš‘›š‘āˆ’Ģ‡š‘šš‘Ÿš‘“ī‚¶=šœ+1šœ.(12) In order to find the temperature profiles in the inner and outer region, the Shvab-Zeldovich form of the energy equation [1] is used, that is, š‘‘ī€·š‘Ÿ2ī€ø(š‘‘š‘‡/š‘‘š‘Ÿ)š‘‘š‘Ÿ=š‘Ģ‡š‘šš‘‘š‘‡š‘‘š‘Ÿ(13) with two sets of BCs for the inner and outer regions BCāˆ’ī‚»š‘‡ī€·š‘Ÿš‘ ī€ø=š‘‡š‘ š‘‡ī€·š‘Ÿš‘“ī€ø=š‘‡š‘“,BC+ī‚»š‘‡ī€·š‘Ÿš‘“ī€ø=š‘‡š‘“,š‘‡(š‘ŸāŸ¶āˆž)=š‘‡āˆž.(14) Integration of (13) and application of inner region BCs (14) gives š‘‡āˆ’=ī€·š‘‡(š‘Ÿ)š‘ āˆ’š‘‡š‘“ī€øexp(āˆ’š‘āˆ’Ģ‡š‘š/š‘Ÿ)+š‘‡š‘“ī€·expāˆ’š‘āˆ’Ģ‡š‘š/š‘Ÿš‘ ī€øāˆ’š‘‡š‘ ī€·expāˆ’š‘āˆ’Ģ‡š‘š/š‘Ÿš‘“ī€øī€·expāˆ’š‘āˆ’Ģ‡š‘š/š‘Ÿš‘ ī€øī€·āˆ’expāˆ’š‘āˆ’Ģ‡š‘š/š‘Ÿš‘“ī€ø(15) for the inner region and š‘‡+ī€·š‘‡(š‘Ÿ)=š‘“āˆ’š‘‡āˆžī€øī€·expāˆ’š‘+ī€øĢ‡š‘š/š‘Ÿ+š‘‡āˆžī€·expāˆ’š‘+Ģ‡š‘š/š‘Ÿš‘“ī€øāˆ’š‘‡š‘“ī€·expāˆ’š‘+Ģ‡š‘š/š‘Ÿš‘“ī€øāˆ’1(16) for the outer region, and (4) into (16) gives š‘‡+=ī€·š‘‡(š‘Ÿ)š‘“āˆ’š‘‡āˆžī€øexp(āˆ’š‘›š‘āˆ’Ģ‡š‘š/š‘Ÿ)+š‘‡āˆžī€·expāˆ’š‘›š‘āˆ’Ģ‡š‘š/š‘Ÿš‘“ī€øāˆ’š‘‡š‘“ī€·expāˆ’š‘›š‘āˆ’Ģ‡š‘š/š‘Ÿš‘“ī€ø.āˆ’1(17) At the surface of the droplet, the heat conducted to it balances the heat used to vaporize and heat up the droplet. Hence, the energy balance at the droplet surface [1] can be written in the form š‘˜āˆ’š‘”4šœ‹š‘Ÿš‘ 2š‘‘š‘‡āˆ’|||š‘‘š‘Ÿš‘Ÿš‘ ī€·ā„Ž=Ģ‡š‘šš‘“š‘”+š‘žš‘–āˆ’š‘™ī€ø.(18) Differentiating (15) and substituting into (18) give the energy balance at the droplet as 4šœ‹š‘˜āˆ’š‘”š‘āˆ’ī€·š‘‡š‘“āˆ’š‘‡š‘ ī€øī€·š‘žš‘–āˆ’š‘™+ā„Žš‘“š‘”ī€øā‹…ī€·expāˆ’š‘āˆ’Ģ‡š‘š/š‘Ÿš‘ ī€øī€ŗī€·expāˆ’š‘āˆ’Ģ‡š‘š/š‘Ÿš‘ ī€øī€·āˆ’expāˆ’š‘āˆ’Ģ‡š‘š/š‘Ÿš‘“ī€øī€»+1=0.(19) At the flame sheet, the heat of combustion is conducted away by both the inner region and outer region gases. Therefore, the energy balance at the flame sheet can be written in the form [1]Ģ‡š‘šĪ”ā„Žš‘=š‘˜āˆ’š‘”4šœ‹š‘Ÿš‘“2š‘‘š‘‡āˆ’|||š‘‘š‘Ÿš‘Ÿš‘“āˆ’š‘˜+š‘”4šœ‹š‘Ÿš‘“2š‘‘š‘‡+|||š‘‘š‘Ÿš‘Ÿš‘“.(20) Differentiating (15) and (17) and substituting into (20) give the simplified energy balance at the flame sheet as1=4šœ‹Ī”ā„Žš‘ā‹…š‘˜āˆ’š‘”š‘āˆ’ī€·š‘‡š‘“āˆ’š‘‡š‘ ī€øī€·expāˆ’š‘āˆ’Ģ‡š‘š/š‘Ÿš‘“ī€øī€ŗī€·expāˆ’š‘āˆ’Ģ‡š‘š/š‘Ÿš‘ ī€øī€·āˆ’expāˆ’š‘āˆ’Ģ‡š‘š/š‘Ÿš‘“āˆ’ī€øī€»4šœ‹Ī”ā„Žš‘ā‹…š‘˜+š‘”š‘›š‘āˆ’ī€·š‘‡āˆžāˆ’š‘‡š‘“ī€øī€·expāˆ’š‘›š‘āˆ’Ģ‡š‘š/š‘Ÿš‘“ī€øī€ŗī€·1āˆ’expāˆ’š‘›š‘āˆ’Ģ‡š‘š/š‘Ÿš‘“.ī€øī€»(21) Solving (12), (19), and (21) for Ģ‡š‘š, š‘Ÿš‘“, and š‘‡š‘“ gives Ģ‡š‘š=4šœ‹š‘˜āˆ’š‘”š‘Ÿš‘ šæš‘’āˆ’š‘āˆ’š‘š‘”Ć—īƒÆīƒ¬1lnī€·š‘˜1+āˆ’š‘”/š‘›šœš‘˜+š‘”ī€ø+ī€·š‘‡āˆžāˆ’š‘‡š‘ ī€øšæš‘’āˆ’š‘āˆ’š‘š‘”/ī€·ī€·š‘˜1+āˆ’š‘”/š‘›šœš‘˜+š‘”ī€øī€ø+Ī”ā„Žš‘/ī€·ī€·š‘›šœš‘˜+š‘”/š‘˜āˆ’š‘”ī€øī€ø+1ī€·š‘žš‘–āˆ’š‘™+ā„Žš‘“š‘”ī€øīƒ­+1š‘›ī‚ƒlnšœ+1šœī‚„īƒ°.(22) Let 1š“=ī€·š‘˜1+āˆ’š‘”/š‘›šœš‘˜+š‘”ī€øā‹…ī‚€šœ+1šœī‚1/š‘›,ī€·š‘‡šµ=āˆžāˆ’š‘‡š‘ ī€øšæš‘’āˆ’š‘āˆ’š‘š‘”/ī€·ī€·š‘˜1+āˆ’š‘”/š‘›šœš‘˜+š‘”ī€øī€ø+Ī”ā„Žš‘/ī€·ī€·š‘›šœš‘˜+š‘”/š‘˜āˆ’š‘”ī€øī€ø+1ī€·š‘žš‘–āˆ’š‘™+ā„Žš‘“š‘”ī€øā‹…ī‚€šœ+1šœī‚1/š‘›,(23) then Ģ‡š‘š=4šœ‹š‘˜āˆ’š‘”š‘Ÿš‘ šæš‘’āˆ’š‘āˆ’š‘š‘”š‘‡ln(š“+šµ),(24)š‘“=ī€·š‘žš‘–āˆ’š‘™+ā„Žš‘“š‘”ī€øšæš‘’āˆ’š‘āˆ’š‘š‘”ī‚øī‚€šœ(š“+šµ)ā‹…ī‚šœ+11/š‘›ī‚¹āˆ’1+š‘‡š‘ ,š‘Ÿš‘“š‘Ÿš‘ =š‘›ln(š“+šµ)[].ln(šœ+1)/šœ(25) The new mass flow rate from (24) was used to obtain burning constant and droplet life time as š¾=8š‘˜āˆ’š‘”šœŒš‘™š‘āˆ’š‘š‘”šæš‘’āˆ’š‘”ln(š“+šµ),š‘‘=š·02š¾.(26)

3. Estimation of Thermal Property Constants

The species properties for the inner region are estimated as follows:š‘‡āˆ’ī€·š‘‡=0.5š‘ +š‘‡š‘“ī€ø,š‘˜(27)āˆ’š‘”=0.4š‘˜š¹ī‚€š‘‡āˆ’ī‚+0.6š‘˜š‘‚š‘„ī‚€š‘‡āˆ’ī‚š‘,(28)āˆ’š‘š‘”=š‘Œāˆ’š¹š‘š‘š¹ī‚€š‘‡āˆ’ī‚+š‘Œāˆ’š‘‚š‘„š‘š‘š‘‚š‘„ī‚€š‘‡āˆ’ī‚.(29) Equation (27) is average temperature in the inner region. Equation (28) has been directly adapted from law and Williamā€™s suggestion [9, 18] since it was experimentally derived, and it replaces the complex estimation of thermal conductivity for a mixture of gas. Equation (29) is the specific heat constant of the gaseous mixture in the inner region. The mass fractions in (29) are estimated by assuming a linear fuel mole fraction from the fuel surface to the region close to the flame where there is a stoichiometric mixture of fuel and oxidizer. Therefore, at the fuel droplet surface, šœ’š¹,š‘ šœ’ā‰ˆ1,š‘‚š‘„,š‘ ā‰ˆ0.(30) At the region close to the flame and for the case of hydrocarbon fuel droplet combustion in air, the reaction equation is typically Cš‘„Hš‘¦+ī‚€š‘¦š‘„+4ī‚ī€·O2+3.76N2ī€øāŸ¶š‘„CO2+(š‘¦/2)H2ī‚€š‘¦O+3.76š‘„+4ī‚N2.(31) The mole fractions can be estimated from the reactants as šœ’š¹,š‘Ÿš‘“āˆ’=1,šœ’1+š‘„+š‘¦/4š‘‚š‘„,š‘Ÿš‘“āˆ’=š‘„+š‘¦/4.1+š‘„+š‘¦/4(32) Therefore, the average mole fractions for the inner region arešœ’āˆ’š¹ī‚€šœ’=0.5š¹,š‘ +šœ’š¹,š‘Ÿš‘“āˆ’ī‚=0.5+0.5,šœ’1+š‘„+š‘¦/4āˆ’š‘‚š‘„ī‚€šœ’=0.5š‘‚š‘„,š‘ +šœ’š‘‚š‘„,š‘Ÿš‘“āˆ’ī‚=0.5š‘„+š‘¦/8.1+š‘„+š‘¦/4(33) The mole fractions can then be used to estimate the mass fractionsš‘Œāˆ’š¹=šœ’āˆ’š¹š‘€š¹šœ’āˆ’š¹š‘€š¹+šœ’āˆ’š‘‚š‘„š‘€š‘‚š‘„,š‘Œāˆ’š‘‚š‘„=šœ’āˆ’š‘‚š‘„š‘€š‘‚š‘„šœ’āˆ’š¹š‘€š¹+šœ’āˆ’š‘‚š‘„š‘€š‘‚š‘„.(34) And the specie constants for the outer region are estimated as follows: š‘‡+ī€·š‘‡=0.5āˆž+š‘‡š‘“ī€ø,š‘+š‘š‘”=š‘š‘š‘‚š‘„ī‚€š‘‡+ī‚,š‘˜+š‘”=š‘˜š‘‚š‘„ī‚€š‘‡+ī‚.(35) The oxidizer (which is usually air) mainly dominates the outer region, so the outer region constants are evaluated directly from the oxidizer properties. Additionally, it is reasonable to assume that š‘žš‘–āˆ’š‘™=0 and š‘‡š‘ =š‘‡boil, since the droplet is burning vigorously after an initial transient heat up.

4. Estimation of Lewis Number

For the estimation of Lewis number in the two regions, unity Lewis number is assumed at the source of diffusing specie. Therefore, in the inner region where fuel diffuses from the droplet surface, unity Lewis number is assumed at the droplet surface šæš‘’āˆ’ī€·š‘‡š‘ ī€ø=1.(36) By definition, šæš‘’āˆ’=š›¼āˆ’Dāˆ’.(37) Assuming ideal-gas behavior, the pressure and temperature dependence of diffusion coefficient [19] is given as š‘‡Dāˆ3/2š‘ƒ(38) which implies that š‘‡D=š¶3/2š‘ƒ=š¶š‘‡3/2,(39) where is š¶ is a constant to be found and š‘ƒ has been absorbed into the constant because the combustion takes place at constant pressure. By applying the boundary condition šæš‘’āˆ’(š‘‡š‘ )=1, š¶ is found as š›¼š¶=āˆ’ī€·š‘‡š‘ ī€øš‘‡š‘ 3/2.(40) Combining (40), (39), and (37) gives šæš‘’āˆ’=š›¼āˆ’ī‚€š‘‡āˆ’ī‚š›¼āˆ’ī€·š‘‡š‘ ī€øā‹…ī‚µš‘‡š‘ š‘‡āˆ’ī‚¶3/2,(41) where š›¼āˆ’ī€·š‘‡š‘ ī€ø=š‘˜š¹ī€·š‘‡š‘ ī€øšœŒš¹āˆ’vaporī€·š‘‡š‘ ī€øā‹…š‘š‘š¹ī€·š‘‡š‘ ī€ø.(42) Fuel dominates at the vapor-surface interface, so density is estimated as šœŒš¹āˆ’vaporī€·š‘‡š‘ ī€ø=š‘ƒš‘‡š‘ š‘…š¹=š‘ƒš‘‡š‘ ā‹…ī€·š‘…š‘¢/š‘€š¹ī€ø,š›¼āˆ’ī‚€š‘‡āˆ’ī‚=š‘˜āˆ’š‘”šœŒāˆ’ā‹…š‘āˆ’š‘š‘”,(43) where š‘˜āˆ’š‘” and š‘āˆ’š‘š‘” are given by (28) and (29), and šœŒāˆ’ is estimated as šœŒāˆ’=š‘ƒš‘‡āˆ’š‘…āˆ’=š‘ƒš‘‡āˆ’ā‹…ī€·š‘…š‘¢/š‘€āˆ’ī€ø,(44) with š‘€āˆ’=š‘Œāˆ’š¹š‘€š¹+š‘Œāˆ’š‘‚š‘„š‘€š‘‚š‘„.(45) Secondly, in the outer region where the combustion products diffuse from the flame sheet, so unity Lewis number is assumed at the flame sheet šæš‘’+ī€·š‘‡š‘“ī€ø=1.(46) Using estimation method similar to that done for the inner region gives šæš‘’+=š›¼+ī‚€š‘‡+ī‚š›¼+ī€·š‘‡š‘“ī€øā‹…ī‚µš‘‡š‘“š‘‡+ī‚¶3/2.(47) At the flame sheet, the mixture fractions can be estimated from the reaction equation ((31) for the combustion of CxHy droplet in air). Ī±s in (47) are estimated as š›¼+ī€·š‘‡š‘“ī€ø=š‘˜Prš‘œductī€·š‘‡š‘“ī€øšœŒPrš‘œductī€·š‘‡š‘“ī€øā‹…š‘š‘Prš‘œductī€·š‘‡š‘“ī€ø,š›¼+ī‚€š‘‡+ī‚=š‘˜+š‘”šœŒ+ā‹…š‘+š‘š‘”.(48) For the case of combustion in air, N2 dominates the product and outer region, and the property of air can be used to estimate Ī±sšæš‘’+ā‰ˆš›¼Airī‚€š‘‡+ī‚š›¼Airī€·š‘‡š‘“ī€øā‹…ī‚µš‘‡š‘“š‘‡+ī‚¶3/2.(49) Equation (49) provides a valid approximation for Lewis number of the outer region and is recommended, since tabulated values of š›¼Air are readily available [1, 20ā€“23].

5. Results and Discussions

As an example, calculation of combustion variables for the case of š‘›-heptane (C7H16) droplet combustion in air was done. Both the simple quasi-steady model and the proposed new model (modified quasi-steady model) were used. Ambient conditions were used, that is, š‘ƒ=1atm and š‘‡āˆž=298K. It was assumed that š‘‡š‘ =š‘‡boil, and droplet heating is negligible, that is, š‘žš‘–āˆ’š‘™=0. š‘‡š‘ =š‘‡boil is assumed because it has been experimentally observed that the droplet boils vigorously during the combustion after an initial and brief heatup, and the heat used to heat the droplet from its initial temperature, š‘žš‘–āˆ’š‘™, is usually negligible and has negligible effect on the model result [1, 2]. Initial guess used for š‘‡š‘“ is 2100ā€‰K since tabulated values of adiabatic flame temperature of common hydrocarbon fuels are approximately 2000ā€‰K.

The calculation results using the simple quasi-steady model are given in Table 1. Law and Williamsā€™s suggestion was used to evaluate species properties for use in the simple quasi-steady model [1, 9, 18], and the results for each iteration step are presented. Iteration was repeated till the solution converged to a difference of 2ā€‰K or less between the guess flame temperature and the calculated flame temperature. Similarly, results obtained by using the proposed new model are presented in Table 2. Additionally, calculation results for the case of combustion of hexane droplet in ambient conditions using the proposed new model are shown in Table 3. The sample calculations showed that proposed new model predicts realistic flame temperature, evaporation constant, and flame to droplet radius ratio compared to the simple quasi-steady model. In the next section, the results of the proposed new model will be compared to experimental results published.

Iteration š‘‡ š‘“ šŗ [K] š‘‡ [K] š‘˜ š‘‚ š‘„ [W/m ā‹… K] š‘˜ š¹ [W/m ā‹… K] š‘˜ g [W/m ā‹… K] š‘ š‘ š¹ [kJ/kg ā‹… K] šµ š‘‡ š‘“ [K] š‘Ÿ š‘“ / š‘Ÿ š‘  K [103ā€‰cm2/s] | š‘‡ š‘“ š‘” āˆ’ š‘‡ š‘“ | [K]


Fixed variables: š‘‡ š‘  = š‘‡ b o i l = 3 7 1 . 5 K , š‘‡ āˆž = 2 9 8 K , šœŒ š‘™ = 6 8 4 ā€‰kg/ m 3 , Ī” ā„Ž š‘ = 4 4 9 2 6 ā€‰kJ/kg, ā„Ž š‘“ š‘” = 3 1 6 ā€‰kJ/kg, and šœ = 1 5 . 0 8 .

Iteration š‘‡ š‘“ šŗ [K] š‘‡ āˆ’ [K] š‘˜ āˆ’ g [ W / m ā‹… K ] š‘ āˆ’ š‘ g [ k J / k g ā‹… K ] š‘‡ + [K] š‘˜ + g [ W / m ā‹… K ] š‘ + š‘ š‘” [ k J / k g ā‹… K ] šæ š‘’ āˆ’ šæ š‘’ + š‘› š“ šµ š‘‡ š‘“ [K] š‘Ÿ š‘“ / š‘Ÿ š‘  K [ 1 0 3 c m 2 / s ] | š‘‡ š‘“ šŗ āˆ’ š‘‡ š‘“ | [K]


Fixed variables: š‘‡ š‘  = š‘‡ b o i l = 3 7 1 . 4 K , š‘‡ āˆž = 2 9 8 K , šœŒ š‘™ = 6 8 4 ā€‰kg/ m 3 , Ī” ā„Ž š‘ = 4 4 9 0 7 ā€‰kJ/kg, ā„Ž š‘“ š‘” = 3 1 6 ā€‰kJ/kg, šœ = 1 5 . 0 8 , šœ’ āˆ’ š¹ = 0 . 5 4 , šœ’ āˆ’ š‘‚ š‘„ = 0 . 4 6 , š‘Œ āˆ’ š¹ = 0 . 8 1 , š‘Œ āˆ’ š‘‚ š‘„ = 0 . 1 9 .

Iteration š‘‡ š‘“ šŗ [K] š‘‡ āˆ’ [K] š‘˜ āˆ’ g [W/m ā‹… K] š‘ āˆ’ š‘ g [kJ/kg ā‹… K] š‘‡ + [K] š‘˜ + g [W/m ā‹… K] š‘ + š‘ š‘” [kJ/kg ā‹… K] šæ š‘’ āˆ’ šæ š‘’ + š‘› š“ šµ š‘‡ š‘“ [K] š‘Ÿ š‘“ / š‘Ÿ š‘  K [ 1 0 3 c m 2 / s ] | š‘‡ š‘“ šŗ āˆ’ š‘‡ š‘“ | [K]


Fixed variables: š‘‡ š‘  = š‘‡ b o i l = 3 4 2 K , š‘‡ āˆž = 2 9 8 K , šœŒ š‘™ = 6 5 9 ā€‰kg/ m 3 , Ī” ā„Ž š‘ = 4 5 0 9 1 ā€‰kJ/kg, ā„Ž š‘“ š‘” = 3 5 5 kJ/kg, šœ = 1 5 . 2 , šœ’ āˆ’ š¹ = 0 . 5 5 , šœ’ āˆ’ š‘‚ š‘„ = 0 . 4 5 , š‘Œ āˆ’ š¹ = 0 . 7 9 , š‘Œ āˆ’ š‘‚ š‘„ = 0 . 2 1 .

6. Comparison with Published Experimental Results

The proposed new model calculation results are compared favorably with published experimental results on the combustion of š‘›-heptane. These experimental results and other models calculation results for the combustion of an š‘›-heptane droplet in air are summarized in Table 4.

Experimental results [2, 3, 9, 24ā€“27] š‘ƒ = 1 atm and š‘‡ āˆž = 293ā€“300ā€‰KProposed model š‘ƒ = 1 a t m and š‘‡ āˆž = 2 9 8 K 1999 model [14] š‘ƒ = 1 a t m and š‘‡ āˆž = 2 9 8 K 1991 model [13] š‘ƒ = 1 a t m and š‘‡ āˆž = 2 9 8 K 1950sā€“1970s model [1ā€“7] š‘ƒ = 1 a t m and š‘‡ āˆž = 2 9 8 K

š‘‡ š‘“ (K)23052559.2263126311142.5
K (103ā€‰cm2/s)7.2ā€“11.47.9314.47.94.74
š‘Ÿ š‘“ / š‘Ÿ š‘  3ā€“1010.0933.2935.2

The proximity between calculated and measured values has been greatly improved by the new model proposed in this paper. The proposed new model accurately predicts š¾ and š‘Ÿš‘“/š‘Ÿš‘ . It predicts š‘‡š‘“ slightly higher and yet closer to the experimentally observed range compared to previous models. The proposed new model also predicted flame temperature closer to the values predicted by the 1999 and 1991 models compared to the original simple quasi-steady model of the 1950s. The š¾ and š‘Ÿš‘“/š‘Ÿš‘  values obtained for the combustion of š‘›-hexane and š‘›-heptane were approximately equal, while š‘‡š‘“ for š‘›-heptane was slightly higher than that for hexane.

The estimated flame temperature seems to have the greatest error, and the most probable source of error is in the estimation of š‘˜š‘“.š‘˜š‘“ at š‘‡āˆ’ was estimated by extrapolation, which is not good for estimating thermal conductivity of vapor, because the available data used ranged up to 1000ā€‰K only. However, since the estimation of flame temperatures by the newest models is approximately 2600ā€‰K, the error may come from experimental error since it may be more difficult to capture the spiked temperature of the flame sheet.

Most of the available tabulated fuel vapor thermal conductivities ranges up to 500ā€‰K, and this points out the need for having thermal conductivities tables or curve fits that ranges up to 1500ā€‰K or higher in order to use and achieve results with less error.

7. Conclusions

The simple quasi-steady model of a fuel droplet was modified to reflect the difference in constant properties across the flame sheet. Two average temperatures were used: one for the inner region and the other for the outer region. The two average temperatures were used to evaluate the assumed constant specific heat and thermal conductivities for the two regions. Nonunity Lewis number was assumed for the two regions while unity Lewis number was assumed at the source of diffusing species, which implies that unity Lewis number was assumed at the flame sheet for the outer region and at the liquid-vapor interface for the inner region. The Lewis numbers obtain in the sample calculation falls within the range that has been observed experimentally [5]. Sample calculations and comparison with experimental results showed that the new model accurately modeled the droplet combustion than the simple quasi-steady model. The new model performance shows that the quasi-steady model of fuel droplet combustion when appropriately applied is a good approximation of the combustion results. The new model also eliminates the need for unguided adjustment of thermal constants and eliminates the need for complex analysis of specific heat and thermal conductivity variation with temperature. The proposed new model is slightly more complex than the original simple quasi-steady model; however, it does not require complex numerical computation for its solutions. The result of the theoretical models of the droplet combustion is best estimated by iteration as shown in the sample calculations. The new model was derive following the process used to derived the old model; hence, it can be noted that if š‘›=1, š‘+š‘š‘”=š‘āˆ’š‘š‘”=š‘š‘š‘”, š‘˜+š‘”=š‘˜āˆ’š‘”=š‘˜š‘”, and šæš‘’+=šæš‘’āˆ’=1 are substituted into the new model equations, the old quasi-steady model will be obtained. The new model can serve as a basis for analysis of other droplet combustion conditions like droplet cloud and convective and high pressure conditions.

Symbols Used

š“:Smaller transfer number (~1)
šµ:Transfer or Spalding number
š‘š‘š‘”:Specific heat constant of gas [J/kgā‹…K]
š·:Droplet diameter [m]
D:Mass diffusivity [m2/s]
ā„Žš‘“š‘”:Latent heat of vaporization [J/kg]
š¾:Evaporation rate constant [m2/s]
š‘˜š‘”:Thermal conductivity of gas [W/mā‹…K]
š‘€:Molecular weight [kg/kmol]
Ģ‡š‘šĢ‡:Fuel mass flow rate [kg/s]
šæš‘’:Lewis number
š‘›:Constant: Z ratio
š‘ƒ:Pressure [atm]
š‘žš‘–āˆ’š‘™:Interface to liquid heat transfer per unit mass (droplet heating) [J/kg]
š‘…:Gas constant [J/Kgā‹…K]
š‘Ÿ:Radius [m]
š‘‡:Temperature [atm]
š‘”:Time [s]
š‘”š‘‘:Droplet life time [s]
š‘„:Number of carbon atoms in fuel molecule
š‘Œ:Mass fraction [kg/kg]
š‘¦:Number of hydrogen atoms in fuel molecule
Ī”ā„Žš‘:Enthalpy of combustion [J/kg]
Ī±:Thermal diffusivity [m2/s]
šœ:Oxidizer-to-fuel stoichiometric mass ratio [kg/kg]
šœŒ:Density [kg/m3]
šœ’:Mole fraction [kmol/kmol].
+:Outer region
āˆ’:Inner region.
0:Initial condition
āˆž:Free streamā€”far from surface
boil:Boiling point
š‘ :Droplet surface


  1. S. R. Turns, An Introduction to Combustion: Concepts and Applications, McGraw Hill, 2nd edition, 2000.
  2. M. Goldsmith and S. S. Penner, ā€œOn the burning of single drops of fuel in oxidizing atmosphere,ā€ Jet Propulsion, vol. 24, pp. 245ā€“251, 1954. View at: Google Scholar
  3. G. A. E. Godsave, ā€œStudies of the combustion of drops in a fuel spray-the burning of single drops of fuel,ā€ Symposium (International) on Combustion, vol. 4, no. 1, pp. 818ā€“830, 1953. View at: Publisher Site | Google Scholar
  4. D. B. Spalding, ā€œThe combustion of liquid fuels,ā€ Symposium (International) on Combustion, vol. 4, no. 1, pp. 847ā€“864, 1953. View at: Publisher Site | Google Scholar
  5. B. N. Raghunandan and H. S. Mukunda, ā€œThe problem of liquid droplet combustion-A reexamination,ā€ Combustion and Flame, vol. 30, no. C, pp. 71ā€“84, 1977. View at: Google Scholar
  6. J. W. Aldred, J. C. Patel, and A. Williams, ā€œThe mechanism of combustion of droplets and spheres of liquid n-heptane,ā€ Combustion and Flame, vol. 17, no. 2, pp. 139ā€“148, 1971. View at: Google Scholar
  7. H. Wise, J. Lorell, and B. J. Wood, ā€œThe effects of chemical and physical parameters on the burning rate of a liquid droplet,ā€ Symposium (International) on Combustion, vol. 5, no. 1, pp. 132ā€“141, 1955. View at: Publisher Site | Google Scholar
  8. L. He, S. D. Tse, and C. K. Law, ā€œRole of flamefront motion and criterion for global quasi-steadiness in droplet burning,ā€ Symposium (International) on Combustion, vol. 2, pp. 1943ā€“1950, 1998. View at: Google Scholar
  9. S. Kumagai, T. Sakai, and S. Okajimah, ā€œCombustion of free fuel droplets in a freelly falling chamber,ā€ Symposium (International) on Combustion, vol. 30, pp. 778ā€“785, 1971. View at: Google Scholar
  10. C. K. Law, S. H. Chung, and N. Srinivasan, ā€œGas-phase quasi-steadiness and fuel vapor accumulation effects in droplet burning,ā€ Combustion and Flame, vol. 38, no. C, pp. 173ā€“198, 1980. View at: Google Scholar
  11. C. K. Law, ā€œRecent advances in droplet vaporization and combustion,ā€ Progress in Energy and Combustion Science, vol. 8, no. 3, pp. 169ā€“199, 1982. View at: Google Scholar
  12. F. A. Williams, Combustion Theory, Benjamin/Cummings, Menlo Park, Calif, USA, 2nd edition, 1985.
  13. I. K. Puri and P. A. Libby, ā€œThe influence of transport properties on droplet burning,ā€ Combustion Science and Technology, vol. 76, pp. 67ā€“80, 1991. View at: Google Scholar
  14. F. F. Filho, ā€œAn analytical solution for the quasi-steady droplet combustion,ā€ Combustion and Flame, vol. 116, no. 1-2, pp. 302ā€“306, 1998. View at: Publisher Site | Google Scholar
  15. R. T. Imaoka and W. A. Sirignano, ā€œTransient vaporization and burning in dense droplet arrays,ā€ International Journal of Heat and Mass Transfer, vol. 48, no. 21-22, pp. 4354ā€“4366, 2005. View at: Publisher Site | Google Scholar
  16. R. T. Imaoka and W. A. Sirignano, ā€œA generalized analysis for liquid-fuel vaporization and burning,ā€ International Journal of Heat and Mass Transfer, vol. 48, no. 21-22, pp. 4342ā€“4353, 2005. View at: Publisher Site | Google Scholar
  17. R. T. Imaoka and W. A. Sirignano, ā€œVaporization and combustion in three-dimensional droplet arrays,ā€ in Proceedings of the 30th International Symposium on Combustion, vol. 30, pp. 1981ā€“1989, July 2004. View at: Publisher Site | Google Scholar
  18. C. K. Law and F. A. Williams, ā€œKinetics and convection in the combustion of alkane droplets,ā€ Combustion and Flame, vol. 19, no. 3, pp. 393ā€“405, 1972. View at: Google Scholar
  19. R. H. Perry, D. W. Green, and J. O. Maloney, Perry’s Chemical Engineers Handbook, McGraw-Hill, New York, NY, USA, 6th edition, 1984.
  20. T. L. Bergman, A. S. Lavine, F. P. Incropera, and D. P. Dewitt, Fundamentals of Heat and Mass Transfer, John Wiley and Sons, 7th edition, 2011.
  21. American Society of Heating, Refrigeration and Air-Conditioning Engineers, ASHRAE Handbook: Fundamentals, ASHRAE, 1981.
  22. N. B. Vargaftik, Tables of Thermophysical Properties of Liquids and Gases, Hemisphere Publishing, New York, NY, USA, 2nd edition, 1975.
  23. A. Kamei, S. W. Beyerlein, and R. T. Jacobsen, ā€œApplication of nonlinear regression in the development of a wide range formulation for HCFC-22,ā€ International Journal of Thermophysics, vol. 16, no. 5, pp. 1155ā€“1164, 1995. View at: Publisher Site | Google Scholar
  24. S. Kumagai and H. Isoda, ā€œCombustion of fuel droplets in a falling chamber,ā€ Symposium (International) on Combustion, vol. 6, no. 1, pp. 726ā€“731, 1957. View at: Publisher Site | Google Scholar
  25. M. T. Monaghan, R. G. Siddall, and M. W. Thring, ā€œThe influence of initial diameter on the combustion of single drops of liquid fuel,ā€ Combustion and Flame, vol. 12, no. 1, pp. 45ā€“53, 1968. View at: Google Scholar
  26. J. W. Aldred and A. Williams, ā€œThe burning rates of drops of n-alkanes,ā€ Combustion and Flame, vol. 10, no. 4, pp. 396ā€“398, 1966. View at: Google Scholar
  27. M. Goldsmith, ā€œExperiments on the burning of single drops of fuel,ā€ Jet Propulsion, vol. 26, pp. 172ā€“178, 1956. View at: Google Scholar

Copyright © 2012 Etim S. Udoetok. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly as possible. We will be providing unlimited waivers of publication charges for accepted research articles as well as case reports and case series related to COVID-19. Review articles are excluded from this waiver policy. Sign up here as a reviewer to help fast-track new submissions.