International Journal of Heat and Mass Transfer xxx (2010) xxx–xxx
Contents lists available at ScienceDirect
International Journal of Heat and Mass Transfer
journal homepage: www.elsevier.com/locate/ijhmt
A simplified model for bi-component droplet heating and evaporation
S.S. Sazhin a,*, A. Elwardany a, P.A. Krutitskii b, G. Castanet c, F. Lemoine c, E.M. Sazhina a, M.R. Heikal a
a
Sir Harry Ricardo Laboratories, Centre for Automotive Engineering, School of Environment and Technology, Faculty of Science and Engineering,
University of Brighton, Brighton BN2 4GJ, UK
b
Keldysh Institute for Applied Mathematics, Department 4, Miusskaya Sq. 4, Moscow 125047, Russia
c
LEMTA, Nancy-Université, CNRS UMR 7563, 2, Avenue de la Forêt de Haye, BP 160, 54504 Vanduvre-lès-Nancy, France
a r t i c l e
i n f o
Article history:
Received 15 December 2009
Accepted 22 June 2010
Available online xxxx
Keywords:
Droplets
Multi-component fuel
Heating
Evaporation
Diffusion equation
Acetone
Ethanol
a b s t r a c t
A simplified model for bi-component droplet heating and evaporation is developed and applied for the
analysis of the observed average droplet temperatures in a monodisperse spray. The model takes into
account all key processes, which take place during this heating and evaporation, including the distribution of temperature and diffusion of liquid species inside the droplet and the effects of the non-unity
activity coefficient (ideal and non-ideal models). The effects of recirculation in the moving droplets on
heat and mass diffusion within them are taken into account using the effective thermal conductivity
and the effective diffusivity models. The previously obtained analytical solution of the transient heat conduction equation inside droplets is incorporated in the numerical code alongside the original analytical
solution of the species diffusion equation inside droplets. The predicted time evolution of the average
temperatures is shown to be reasonably close to the measured one, especially in the case of pure acetone
and acetone-rich mixture droplets. It is shown that the temperatures predicted by the simplified model
and the earlier reported vortex model are reasonably close. Also, the temperatures predicted by the ideal
and non-ideal models differ by not more than several degrees. This can justify the application of the simplified model with the activity coefficient equal to 1 for the interpretation of the time evolution of temperatures measured with errors more than several degrees.
Crown Copyright Ó 2010 Published by Elsevier Ltd. All rights reserved.
1. Introduction
In a number of papers a new model for heating and evaporation
of monocomponent droplets, based on the analytical solution of
the heat conduction equation inside droplets, has been described
[1–3]. The effects described by this model were shown to be particularly important for modelling fuel spray autoignition processes
[4,5]. It was shown to be more CPU efficient compared with the
model based on the numerical solution of the heat conduction
equation inside droplets [6], and was validated against the results
of available experimental measurements of droplet temperatures
[7,8].
The practical application of the above-mentioned model, however, is limited by the fact that most real-life fuels are multi-component. The models of multi-component droplet heating and
evaporation could be subdivided into two main groups: those
based on the analysis of individual components [9–15], applicable
in the case when a small number of components needs to be taken
into account, and those based on the probabilistic analysis of a
large number of components (e.g. continuous thermodynamics ap* Corresponding author. Tel.: +44 (0)1273642677; fax: +44 (0)1273642309.
E-mail address: S.Sazhin@brighton.ac.uk (S.S. Sazhin).
proach [16–18] and the Distillation Curve Model [19]). In the second family of models a number of additional simplifying
assumptions were used, including the assumption that species inside droplets mix infinitely quickly. A model containing features of
both these groups of models has been suggested in [20].
The focus of this paper is on the extension of the model developed in [2,3] for monocomponent droplets to the case of multicomponent droplets. Only the first group of multi-component
droplet heating and evaporation models is considered, which allows us to perform the deterministic analysis of individual species.
Moreover, at this stage only the simplest case of bi-component
droplets is considered, as in [12,13,15], which allows us to develop
a better understanding of the underlying physics of the processes
involved. The new model is much simpler than the previously suggested models (e.g. [15]) which makes it potentially attractive for
implementation into computational fluid dynamics (CFD) codes.
The predictions of the model are compared with the measured
time evolution of droplet temperatures in monodisperse bi-component (ethanol/acetone) droplet streams.
The basic equations and approximations of the new model are
described in Section 2. The experimental set-up for measurements
of droplet temperatures in monodisperse bi-component (ethanol/
acetone) droplet streams is briefly described in Section 3. In
0017-9310/$ - see front matter Crown Copyright Ó 2010 Published by Elsevier Ltd. All rights reserved.
doi:10.1016/j.ijheatmasstransfer.2010.06.044
Please cite this article in press as: S.S. Sazhin et al., A simplified model for bi-component droplet heating and evaporation, Int. J. Heat Mass Transfer (2010),
doi:10.1016/j.ijheatmasstransfer.2010.06.044
2
S.S. Sazhin et al. / International Journal of Heat and Mass Transfer xxx (2010) xxx–xxx
Nomenclature
BM
BT
c
C
D
F
GE
h
h0
k
kB
L
m
M
n
Nu
p
Pe
Pr
qn
QL
Qn
R
Ru
Rd
Re
Sc
Sh
t
T
Te 0 ðRÞ
vA
X
Y
Spalding mass transfer number
Spalding heat transfer number
specific heat capacity
distance parameter
diffusion coefficient
function introduced in Eq. (8) or (22)
Gibbs free energy per unit mole
convection heat transfer coefficient
h
R 1 or Dal Rd 1
kl d
thermal conductivity
Boltzmann constant
latent heat of evaporation
mass
molar mass
molar concentration
Nusselt number
pressure
Peclet number
Prandtl number
parameter introduced in Eq. (3)
power spent on droplet heating
parameter introduced by Eq. (14)
distance from the droplet centre
universal gas constant
droplet radius
Reynolds number
Schmidt number
Sherwood number
time
temperature
RTd0(R)
molar volume
molar fraction
mass fraction
b
d
c
ev,a
g
j
jR
k, kn
l
l0(t)
m
q
r
u
v
vY
Greek symbols
parameter defined by Eq. (11) or introduced in Eq. (44)
parameter introduced in Eq. (44)
parameter introduced in Eq. (44)
activity coefficient
evaporation rate
Lennard–Jones energy
function defined by Eq. (28)
thermal diffusivity
parameter introduced in Eq. (3)
eigenvalues
dynamic viscosity
variable introduced in Eq. (3)
kinematic viscosity
density
Lennard–Jones length
parameter defined by Eq. (24) or associated factor of solvent
parameter defined by Eq. (7)
parameter defined by Eq. (18)
Subscripts
ac
acetone
amb
ambient
b
boiling
d
droplet
eff
effective
eth
ethanol
g
gas
i
species
iso
isolated
l
liquid
p
constant pressure
s
surface
v
vapour
0
initial
1
ambient conditions
a
Section 4, some results of the numerical solution of the basic equations for the values of parameters relevant to the experimental setup described in Section 3 are presented. The main results of the paper are summarised in Section 5.
2. Basic equations and approximations
The model developed in the paper is based on the equations
describing, liquid phase heating, species diffusion in the liquid
phase, droplet evaporation, and species mass fractions at the surface of the droplets. These equations and their approximations
and analytical solutions, where appropriate, are presented and discussed below.
2.1. Droplet heating
hðT g T s Þ ¼ kl
@T
;
@RR¼Rd 0
ð2Þ
and the initial condition T(t = 0) = Td0(R), where Ts = Ts(t) is the droplet’s surface temperature, Tg = Tg(t) is the ambient gas temperature,
h is the convection heat transfer coefficient, linked with the Nusselt
number Nu via the following equation:
Nu ¼ 2Rd h=kg ;
As in the case of monocomponent droplets [1–3], the process of
heating (or cooling) for stationary spherically-symmetric multicomponent droplets is described by the following transient heat
conduction equation for the temperature T T(t, R) in the liquid
phase [21,22]:
!
@T
@ 2 T 2 @T
;
¼j
þ
@t
@R2 R @R
where j = kl/(clql) is the liquid thermal diffusivity, kl, cl, and ql are
the liquid thermal conductivity, specific heat capacity, and density,
respectively, R is the distance from the centre of the spherical droplet, t is time.
Eq. (1) inside the droplet is to be solved for t > 0 and 0 6 R < Rd,
where Rd is the droplet radius, with the following boundary
condition:
ð1Þ
where kg is the gas thermal conductivity.
Remembering the physical background to the problem, we are
interested only in a solution which is continuously differentiable
twice in the whole domain. This implies that T should be bounded
for 0 6 R < Rd.
Assuming that h = const., the solution of Eq. (1), subject to the
above boundary and initial conditions, can be presented as [2]:
Please cite this article in press as: S.S. Sazhin et al., A simplified model for bi-component droplet heating and evaporation, Int. J. Heat Mass Transfer (2010),
doi:10.1016/j.ijheatmasstransfer.2010.06.044
S.S. Sazhin et al. / International Journal of Heat and Mass Transfer xxx (2010) xxx–xxx
TðR; tÞ ¼
(
!
1
R2 sin kn
1X
l0
qn exp jR k2n t d 2 2 l0 ð0Þ 0 2
R n¼1
jR kn
kv n k kn
)
2
0
R
sin
k
l
R
n
þ T g ðtÞ; ð3Þ
exp jR k2n t d 2 4 0 sin kn
Rd
kv n k kn jR
where kn are solutions of the equation:
k cos k þ h0 sin k ¼ 0;
!
Rd
sin 2kn
Rd
h0
2
¼
kv n k ¼
;
1
1þ 2
2
2kn
2
h0 þ k2n
Z Rd
1
R
Te 0 ðRÞ sin kn
dR;
qn ¼
2
Rd
kv n k 0
jR ¼
kl
2
l Rd
cl q
l0 ðtÞ ¼
;
T eff ¼ T g þ
h
;
BT ¼
ð5Þ
cpv ðT g T s Þ
Leff
FðBT Þ ¼ ð1 þ BT Þ0:7
QL
Leff ¼ L þ m
_ ¼
d
P
lnð1 þ BT Þ
;
BT
i i Li
þ PQ Lm_ ; Q L is the power spent on droplet heati
i
ing, cpv is the specific heat capacity of fuel vapour, i = i(t) are the
P
_ d ¼ im
_ d ðm
_ i Þ.
_ i ¼ i m
evaporation rates of species i; m
This model is equally applicable to one-component and multicomponent droplets. The liquid thermal conductivity of multicomponent droplets is calculated as described in Appendix A.
2.2. Species diffusion in the liquid phase
Assuming that the processes inside droplets are spherically
symmetric (no convection), equations for mass fractions of liquid
species Yli Yli(t, R) there can be presented in the following form:
!
@Y li
@ 2 Y li 2 @Y li
þ
;
¼ Dl
@t
R @R
@R2
ð9Þ
where i = 1, 2, Dl is the liquid mass diffusivity.
Eq. (9) is solved with the following boundary condition [26,5]:
aði Y lis Þ ¼ Dl
@Y li
;
@R R¼Rd 0
ð10Þ
and the initial condition Yli(t = 0) = Yli0(R), where Ylis = Ylis(t) are liquid components’ mass fractions at the droplet’s surface,
RdðnewÞ ¼ RdðoldÞ þ R_ d Dt;
where the value of R_ d is controlled by the droplet evaporation model (see Eq. (27) in Section 2.3). Note that updating of liquid density
can lead to Rd(new) > Rd(old) (thermal swelling, see Eq. (27)).
The average surface temperature in a moving droplet can still
be correctly predicted by Eq. (1), with appropriate boundary conditions, and its solution (3) if the liquid thermal conductivity kl is replaced by the so-called effective thermal conductivity keff [23]:
ð6Þ
where the coefficient v varies from 1 (at droplet Peclet number Ped(l) = Red(l)Prl < 10) to 2.72 (at Ped(l) > 500) and can be approximated
as:
v ¼ 1:86 þ 0:86 tanh½2:225log10 PedðlÞ =30 :
ð8Þ
is the Spalding heat transfer number,
where L is the latent heat of evaporation, the value of R_ d can be taken from the previous time step and estimated based on Eq. (27). Rd
is assumed constant in the analytical solutions, but is updated at the
end of the time step Dt:
keff ¼ vkl ;
lnð1 þ BT Þ
BT
h
i
0
1
1
ð1 þ Red Prd Þ1=3 max 1; Re0:077
d
A;
@1 þ
2FðBT Þ
Nuiso ¼ 2
where
e 0 ðRÞ ¼ RT d0 ðRÞ. The solution of Eq. (4) gives a set
h0 ¼ ðhRd =kl Þ 1; T
of positive eigenvalues kn numbered in ascending order (n = 1, 2, . . .).
The trivial solution k = 0 is not considered.
Note that the norm kvnk2 differs from the norm chosen in [2] by
the factor Rd (see their Eq. (A4)). This does not affect the final
solution.
Solution (3) is valid for h0 > 1, which is satisfied, remembering
the physical background of the problem (h > 0). The condition
h = const. is valid for sufficiently small time steps.
When deriving Solution (3) it was assume that it is applied to
individual short time steps. In this case we can ignore the time
dl0 ðtÞ
dependence
D
Eof dt during the time steps and assume that
dl0 ðtÞ
dl0 ðtÞ
0
l
0 . Note that l0 depends on time via Tg only.
dt
dt
To take into account the effect of droplet evaporation in analytical solution (3), gas temperature should be replaced by the socalled effective temperature defined as:
ql LR_ d
tion of Solution (3) to the case of moving droplets by replacing kl
with keff, assuming that h = const.
The value of Nu for an isolated moving droplet is estimated
based on the following equation [23]:
ð4Þ
hT g ðtÞRd
;
kl
3
ð7Þ
Liquid fuel transport properties were used for calculating Ped(l). The
model based on Eqs. (6) and (7) is known as the effective thermal
conductivity model. This model was developed mainly for the estimation of the average surface temperature of droplets, which controls droplet evaporation (cf. [24,25]). It cannot predict adequately
the details of the distribution of temperature inside droplets, which
include vortex structures for non-zero droplet velocities, but these
are not required in most practical engineering applications. Hence,
the applicability of this model can be justified. It allows the applica-
a¼
_ dj
jm
4pql R2d
;
ð11Þ
_ d is the droplet evaporation rate, the calculation of which is dism
cussed in the next section (see Eq. (20)).
Assuming that species concentrations in the ambient gas are
equal to zero, the values of i can be found from the following relation [9–11]:
Y
i ¼ P v is
i Y v is
;
ð12Þ
where the subscript v indicates the vapour phase. We assume that i
is constant and is determined by the values of Yvjs (j refers to all species) at the beginning of the time step. The conditions i = const. and
a = const. can always be guaranteed for sufficiently small time
steps.
As in the case of Eq. (1), we are interested only in a solution
which is continuously differentiable twice in the whole domain.
This implies that Yli should be bounded for 0 6 R < Rd. Moreover,
the physical meaning of Yli, as the mass fraction, implies that
0 6 Yli 6 1.
Eq. (9) with boundary conditions (10) have essentially the same
structure as Eq. (1) with boundary condition (2). The former equations can be obtained from the latter if we replace T with Yli, j with
Please cite this article in press as: S.S. Sazhin et al., A simplified model for bi-component droplet heating and evaporation, Int. J. Heat Mass Transfer (2010),
doi:10.1016/j.ijheatmasstransfer.2010.06.044
4
S.S. Sazhin et al. / International Journal of Heat and Mass Transfer xxx (2010) xxx–xxx
Dl, Tg with i, Ts with Ylis and kl/h with Dl/a. However, although Eq.
(9) look rather similar to Eq. (1), the solution of Eq. (1) (see Expression (3)) cannot be used for (9). The reason for this is that solution
(3) is valid only for h0 > 1. At the same time the boundary condition for Eq. (9) at R = Rd (see Eq. (10)) leads to the situation in which
h0 < 1. The solution of Eq. (9) in this case can be presented as
("
" #
2
1
k0
R
Y li ¼ i þ
t ½qi0 Q 0 i sinh k0
exp Dl
Rd
R
Rd
"
#
"
#)
2
1
X
kn
R
exp Dl
t ½qin Q n i sin kn
;
þ
Rd
Rd
n¼1
ð13Þ
where k0 and kn (n P 1) are solutions of equations
tanh k ¼
k
h0
and
tan k ¼
k
:
h0
2.3. Droplet evaporation
In the case of isolated moving monocomponent droplets, their
evaporation rate is given by the following equation [5]:
_ d ¼ 2pRd Dv qtotal BM Shiso ;
m
where Dv is the binary diffusion coefficient of vapour in air and BM is
the Spalding mass transfer number defined as:
q qv 1 Y v s Y v 1
BM ¼ v s
¼
;
qgs
1 Y vs
Shiso
h
i
0
1
ð1 þ Red Scd Þ1=3 max 1; Re0:077
1
d
lnð1 þ BM Þ @
A;
¼2
1þ
BM
2FðBM Þ
ð22Þ
l
Qn ¼
qin ¼
>
:
2
kv 0 k
1
kv n k2
1
kv n k2
Z
2
Rd
k0
1
Rd
kn
2
ð1 þ h0 Þ sinh k0
ð1 þ h0 Þ sin kn
when n ¼ 0;
when n P 1;
ð14Þ
Rd
0
RY li0 ðRÞv n ðRÞdR;
ð15Þ
v 0 ðRÞ ¼ sinh
k0
R
;
Rd
R
v n ðRÞ ¼ sin kn ; n P 1:
Rd
ð16Þ
The average mass fraction of species in a moving droplet can
still be correctly predicted by Eq. (9), with appropriate boundary
conditions, if the liquid diffusivity Dl is replaced with the so-called
effective diffusivity Deff [27,26]:
Deff ¼ vY Dl ;
ð17Þ
where the coefficient vY varies from 1 to 2.72 and can be approximated as:
vY ¼ 1:86 þ 0:86 tanh½2:225log10 ðRel Scl =30Þ;
ð18Þ
Scl is the liquid Schmidt number defined as:
Scl ¼
ml
Dl
;
mair
Dv
is the Schmidt number,
FðBM Þ ¼ ð1 þ BM Þ0:7
lnð1 þ BM Þ
:
BM
BT and BM are linked by the following equation [23]:
where
If Eq. (13) is applied to a short time step it can be simplified to:
("
2 #
1
k0
R
Y li ¼ i þ
t ½qi0 Q 0 i sinh k0
1 þ Dl
Rd
R
Rd
"
#
)
2
1
X
kn
R
1 Dl
t ½qin Q n i sin kn
:
þ
Rd
Rd
n¼1
Scd ¼
BT ¼ ð1 þ BM Þu 1;
n P 0, where
ð21Þ
Yv is the vapour mass fraction and Shiso is the Sherwood number
approximated for the isolated droplets by the following expression:
respectively, h0 ¼ 1 þ aDRd .
8
>
<
ð20Þ
ð19Þ
ml is the liquid kinematic viscosity.
Eq. (17) allows the application of Solution (13) to the case of
moving droplets by replacing Dl with Deff, assuming that a = const.
Following [26] this model is called the effective diffusivity model.
As in the case of the effective thermal conductivity model, this
model cannot describe the details of species mass fractions inside
droplets, including vortex structures in the moving droplets, but
this information is not required in most engineering applications.
A more complex approach based on the analysis of Hill vortices
is discussed in a number of publications including [26,15].
The liquid kinematic viscosity and diffusivity in multi-component droplets are calculated as described in Appendix A.
u¼
cpv
cp;
air
1
;
Le
ð23Þ
ð24Þ
Le = kg/(cp,airqtotalDv) is the Lewis number.
In the case of multi-component droplets the problem of modelling droplet evaporation is complicated by the fact that different
species diffuse at different rates, and the evaporation rate of one
of the species is affected by the evaporation rate of other species.
In [15] the analysis of evaporation of multi-component droplets
_ d:
led to the following expression for m
_ d ¼ 2pRd Div qtotal BMi ShisoðiÞ ;
m
ð25Þ
where BMi is the species Spalding mass transfer number defined as:
BMi ¼
Y iv s Y iv 1
;
i Y iv s
ð26Þ
Div is the diffusion coefficient of species i in air, Shiso(i) is defined by
(22) replacing BM by BMi and Dv by Div. Remembering (12) one can
see that BM = BMi. Hence, for stationary droplets this leads to the
_ d is predicted by Eq. (25) for differparadox that the same value of m
ent Div. This paradox is resolved by the fact that although Eq. (25) is
correct, the value of Shiso(i) cannot be approximated by the analogue
of Eq. (22) which is implicitly based on the assumption that the
evaporating species do not affect each other.
_ d is based on Eq. (20), assuming that the mixOur analysis of m
ture of vapour species can be treated as a separate gas, similar to
treating the mixture of nitrogen, oxygen and carbon dioxide as
P
air ðY v s ¼ i Y v is Þ. The value of Dv is estimated as described in
Appendix A.
When calculating the value of R_ d we took into account both
droplet evaporation during the time step, and the change in their
density during this time step (see Eq. (20) in [3]):
2
3
!1=3
R
ðT
Þ
q
ðT
Þ
0 4
0
d
R_ d ¼
þ
15;
Dt
qðT 1 Þ
4pR2d ql
_d
m
ð27Þ
where T 0 and T 1 are average droplet temperatures at the beginning
t = t0 and the end of the time step t = t1, Dt = t1 t0.
Please cite this article in press as: S.S. Sazhin et al., A simplified model for bi-component droplet heating and evaporation, Int. J. Heat Mass Transfer (2010),
doi:10.1016/j.ijheatmasstransfer.2010.06.044
S.S. Sazhin et al. / International Journal of Heat and Mass Transfer xxx (2010) xxx–xxx
The corrections to Nuiso and Shiso due to the finite distance between droplets are calculated as in [7,35]:
Sh
Nu
¼
Shiso
Nuiso
1 0:57 exp½0:13ðC 6Þ
;
¼ 1 0:57 1
1 þ 0:57 exp½0:13ðC 6Þ
gðCÞ ¼
ð28Þ
where C is the ratio of the distance between droplets and their
diameters (distance parameter).
2.4. Species mass fractions at the surface of the droplets
_ i and the values
To calculate the species mass evaporation rate m
of the evaporation rate of species i based on Eq. (12) we need to
calculate first the values of Yvis. The latter depends on the partial
pressure of species i in the vapour state in the immediate vicinity
of the droplet surface [28]:
pv i ¼ X lis ci pv i ;
ð29Þ
where Xlis is the molar fraction of the ith species in the liquid near
the droplet surface, pv i is the partial vapour pressure of the ith species in the case when Xlis = 1, ci is the activity coefficient. The subscript s is omitted in Eq. (29) and all the following equations in
this subsection.
In the limit when ci = 1, Eq. (29) describes the Raoult law. A
more accurate approximation for ci for the ethanol–acetone mixture is described in Appendix B, following [29,30,15]. This approximation will be used in our analysis.
Remembering the Clausius–Clapeyron equation [28], Eq. (29)
can be re-written as:
pv i ¼ X lis ci pamb exp
Li M i 1
1
;
Ru T bi T s
ð30Þ
where Mi is the molar mass and Tbi is the boiling temperature of the
ith species, pamb is the ambient pressure. When deriving Eq. (30) it
was taken into account that pv i is equal to the ambient pressure
when Ts = Tbi, Li is the latent heat of evaporation of species i. Note
that there are typos in Eqs. (135) and (136) of [5].
The numerical values of LiMi at T = Tbi, values of Tbi for ethanol
and acetone at atmospheric pressure, and Mi are presented in Table 1 [31].
The values of Li for ethanol and acetone (in J/kg) at arbitrary T
are calculated based on the following equations [30,34]:
Leth ¼ 120:91 103 ð516:2 T Þ0:38 ;
0:38
508:1 T
;
Lac ¼ 489 103
508:1 T b;ac
ð31Þ
ð32Þ
where T is in K, Tb,ac is the boiling temperature of acetone given in
Table 1.
Eq. (30) can be re-written for the molar fractions of species i in
the vapour phase:
X v i ¼ X lis ci exp
Li M i 1
1
:
Ru T bi T s
ð33Þ
3. Experimental set-up and input parameters
5
generated by a Rayleigh-type disintegrating liquid jet. The fuel
temperature was regulated in the injector body and was measured
near the outlet by a thermocouple. Downstream distance from the
injector was converted into time with the help of the space evolution of the droplet velocity measured using the Laser Doppler Anemometry (LDA) technique. Droplets were injected into a quiescent
atmosphere at room temperature. Droplet temperatures were
measured using the technique described in [7,15,30]. The input
parameters for the model were the initial droplet temperature
(assuming that this temperature is initially homogeneous), ambient gas temperature (assuming that this temperature remains constant during the experiment), the distance parameter (ratio of the
distances between droplets and their diameters) and the droplet
velocities.
The measured time evolution of the droplet velocities for pure
acetone is shown in Fig. 1. In the same figure, the linear approximation of the experimental results is shown. Similar plots were obtained for pure ethanol and various mixtures of acetone and
ethanol (not shown). In all cases the linear approximations of the
experimental results, similar to the one shown in Fig. 1, were used.
These are summarised in Table 2.
The two-colour laser-induced fluorescence (2cLIF) technique,
previously used for the measurement of droplet temperatures,
was further developed to include the third colour band. This additional band is required for the analysis of bi-component droplets
[32]. In this new technique, droplets were seeded with a low concentration of a fluorescent dye and illuminated by laser beams having a wavelength tuned on the absorption spectrum of the tracer.
The fluorescence signal was detected over specific spectral bands
of emissions. The ratio of the intensity of these bands depends
on the temperature and, to a lesser extent, on the droplet composition, while the fluorescence dependencies on tracer concentration, probe volume dimensions, laser intensity and optical layout
were eliminated [32]. The probe volume (intersection between
the laser beams and the detection field of view) was about
150 150 1200 lm3. It was larger than the droplet in order to
provide a global excitation of the whole droplet volume. The signal
was averaged over the total time of droplet transit in the probe
volume.
In the presence of thermal gradients inside droplets, the temperature obtained based on this signal could differ from the volume-averaged droplet temperature. This is related to the fact
that the laser intensity is not distributed uniformly inside the
droplet, mainly due to light refraction at the droplet surface. Also,
the receiving optics collect the fluorescent signal from the droplet
interior with a spatial distortion. These effects could be taken into
account with the help of an efficient ray tracing code, based on the
geometric optics approximation [33]. This analysis, however, is beyond the scope of this paper. We would just mention that the 2cLIF
measurements allow us to estimate the volume-averaged droplet
temperature with systematically more weighting in the zones near
the centre of the droplet rather than at its edge. We cannot quantify this effect at this stage and we assume that the measured temperatures are close to the droplet average temperatures.
4. Results
The experimental set-up used for validation of the model is the
same as described in [7,15,30]. A monodisperse droplet stream was
Table 1
The values Li M i at T ¼ T bi ; T bi and M i for ethanol and acetone.
Substance
LiMi (kJ/mol)
Tbi (K)
Mi (kg/kmol)
Ethanol (C2H6O)
38.56
351.80
46.069
Acetone (C3H6O)
29.10
329.22
58.080
The plots of the time dependence of the temperatures for pure
acetone and ethanol are shown in Figs. 2 and 3, respectively. The
experimentally observed average droplet temperatures are shown
by filled triangles in these figures. The results of calculations are
shown by solid, dashed and dashed dotted curves referring to the
central, average and surface temperatures. The calculations started
with the first observed droplets approximately 1 ms after the start
of injection. At the earlier times, the liquid fuel formed an unstable
Please cite this article in press as: S.S. Sazhin et al., A simplified model for bi-component droplet heating and evaporation, Int. J. Heat Mass Transfer (2010),
doi:10.1016/j.ijheatmasstransfer.2010.06.044
6
S.S. Sazhin et al. / International Journal of Heat and Mass Transfer xxx (2010) xxx–xxx
13
Experiment
Approximation
Ud (m/s)
12
11
100% acetone
10
0
1
2
3
4
5
6
7
Time (ms)
Fig. 1. Experimentally observed velocities for pure acetone droplets (squares). These were approximated by the straight line Ud = 12.81 0.316 t, where t is measured from
the moment of injection.
evolution of droplet radii were not available. Note that the effective
thermal conductivity and effective diffusivity models, on which our
analysis is based, are primarily designed to predict correctly the
average surface temperature and species mass fractions of droplets, but not their average temperature and species mass fractions.
In the case of ethanol, shown in Fig. 3, the experimentally observed temperatures lie between average and surface temperatures predicted by the model. Note that the temperature scale in
Fig. 3 is much finer than in Fig. 2. Hence the overall agreement between experimental data and predictions of the model looks reasonably good. This level of agreement between the model and
experimental data is similar to the one reported in [15] for the case
of the vortex model (see their Fig. 7). Hence, the application of our
simplified model can be justified in this case. Note that despite the
overall agreement between the results of modelling and experimental data, the observed rate of temperature decrease at the later
times is lower than predicted by the model.
Plots similar to those shown in Figs. 2 and 3, but for the mixtures of ethanol and acetone, are presented in Figs. 4–6. The calcu-
Table 2
The approximations of U drop .
Substance
Approximation of Udrop in m/s (t is in ms)
100% acetone
100% ethanol
25% ethanol + 75% acetone
50% ethanol + 50% acetone
75% ethanol + 25% acetone
12.81–0.316t
12.30–0.344t
12.75–0.370t
12.71–0.448t
12.28–0.306t
jet the temperature evolution of which cannot be interpreted using
our model. Also, we assumed that there is no temperature gradient
inside droplets at the initial moment of time used in calculations.
As follows from these figures, all three temperatures are well
separated for both acetone and ethanol. Hence, the difference between them needs to be taken into account in the analysis of
experimental data (cf. [36,37]). In the case of acetone there seems
to be a reasonable agreement between the values of average temperature of the droplets and experimental data. Data on the time
40
Experiment
Ts
Ts
35
Tav
Tav
Tc
Tc
T (C)
30
25
20
100% acetone
15
10
0
1
2
3
4
5
6
7
Time (ms)
Fig. 2. The time evolution of droplet surface, average and centre temperatures (Ts, Tav, and Tc) and experimentally observed temperatures for pure acetone droplets with
initial diameter 143.4 lm and homogeneous temperature 35.1 °C in an ambient gas at temperature equal to 21.5 °C.
Please cite this article in press as: S.S. Sazhin et al., A simplified model for bi-component droplet heating and evaporation, Int. J. Heat Mass Transfer (2010),
doi:10.1016/j.ijheatmasstransfer.2010.06.044
7
S.S. Sazhin et al. / International Journal of Heat and Mass Transfer xxx (2010) xxx–xxx
40
Experiment
TsTs
Tav
Tav
T (C)
35
TcTc
30
25
100% ethanol
20
0
1
2
3
4
Time (ms)
5
6
7
Fig. 3. The same as Fig. 2 but for pure ethanol droplets with initial diameter 140.8 lm, homogeneous temperature 38.0 °C and ambient gas temperature equal to 22.0 °C.
lations were performed for the cases of the ideal mixture (ci = 1 in
Eqs. (29) and (30)) and the non-ideal mixture (ci in Eqs. (29) and
(30) was calculated based on Eq. (45)). As can be seen from these
figures, in all cases the predictions of the temperatures by the ideal
and non-ideal models are noticeably different (by up to several degrees), especially at later times. However, both these models predict about the same trend in the evolution of temperature with
time. The ideal model can be used if the prediction errors of several
degrees can be tolerated. This seems to be our case where the random errors of the estimates of droplet temperatures seem to be
about 2–3°.
In the case of an acetone dominated mixture (25% ethanol–75%
acetone: see Fig. 4), the agreement between the observed and predicted average droplet temperatures, for both ideal and non-ideal
models, is reasonably good, although the scatter of experimental
data in this case is more noticeable than in the case of pure acetone
(see Fig. 2). In the case of the 50% ethanol–50% acetone mixture
(see Fig. 5), the experimentally observed temperatures lie between
the average and surface temperatures predicted by both ideal and
non-ideal models, as in the case of pure ethanol shown in Fig. 3.
These temperatures are closer to the average temperatures predicted by the non-ideal model than to those predicted by the ideal
model. In the case of the 75% ethanol–25% acetone mixture (see
Fig. 6), the experimentally observed temperatures lie well below
the average temperatures predicted by both ideal and non-ideal
models, although they are closer to the average temperatures predicted by the non-ideal model than to those predicted by the ideal
model, as in the case of the 50% ethanol–50% acetone mixture. The
reason for this deviation between the measured and predicted
temperatures, which could reach up to 5 K, is not clear to us.
In Fig. 7 we have compared the predictions of our simplified
model and the predictions of the vortex model, reported in [15].
As one can see from this figure, the experimental results agree better with the predictions of the average temperature by the vortex
model than by our simplified model. It seems, however, that this
was achieved by the choice of the lower initial droplet temperature
in the vortex model (which could be used as a fitting parameter). If
the values of this temperature were taken to be the same, one
40
25 % ethanol -75 % acetone
35
T (C)
30
25
Experiment
20
15
non-ideal-Ts
Ts (non-ideal)
ideal-Ts
Ts (ideal)
non-ideal-Tav
Tav (non-ideal)
Tav (ideal)
ideal-Tav
Tc (non-ideal)
non-ideal-Tc
Tc (ideal)
ideal-Tc
10
0
1
2
3
4
Time (ms)
5
6
7
Fig. 4. The same as Figs. 2 and 3 but for the 25% ethanol–75% acetone mixture droplets with initial diameter 133.8 lm, homogeneous temperature 32.5 °C and ambient gas
temperature equal to 21.1 °C. The results of calculations based on the ideal (ci = 1) and non-ideal models are presented.
Please cite this article in press as: S.S. Sazhin et al., A simplified model for bi-component droplet heating and evaporation, Int. J. Heat Mass Transfer (2010),
doi:10.1016/j.ijheatmasstransfer.2010.06.044
8
S.S. Sazhin et al. / International Journal of Heat and Mass Transfer xxx (2010) xxx–xxx
40
50% ethanol + 50% acetone
35
T (C)
30
non-ideal-Ts
Ts (non-ideal)
non-ideal-Tav
Tav (non-ideal)
25
non-ideal-Tc
Tc (non-ideal)
Ts (ideal)
ideal-Ts
20
Experiment
Tav (ideal)
ideal-Tav
Tc (ideal)
ideal-Tc
15
0
1
2
3
4
Time (ms)
5
6
7
Fig. 5. The same as Fig. 4 but for the 50% ethanol–50% acetone mixture droplets with initial diameter 142.7 lm, homogeneous temperature 37.5 °C and ambient gas
temperature equal to 20.8 °C.
40
75% ethanol + 25% acetone
T (C)
35
30
non-ideal-Ts
Ts (non-ideal)
non-ideal-Tav
Tav (non-ideal)
non-ideal-Tc
Tc (non-ideal)
25
ideal-Ts
Ts (ideal)
Experiment
Tav (ideal)
ideal-Tav
Tc (ideal)
ideal-Tc
20
0
1
2
3
4
Time (ms)
5
6
7
Fig. 6. The same as Figs. 4 and 5 but for the 75% ethanol–25% acetone mixture droplets with initial diameter 137.1 lm, homogeneous temperatures 38.6 °C and ambient gas
temperature equal to 21.6 °C.
would expect that the predictions of the simplified and vortex
models would be very close.
The plots of time evolution of droplet radius for the same mixture as in the case of Figs. 6 and 7, inferred from our simplified
model and the numerical results reported in [15], are shown in
Fig. 8. As follows from this figure, the vortex model predicts a
slightly higher evaporation rate than the simplified model (the difference in predicted radii is less than 0.3 lm). In both cases, the
non-ideal model predicts a slightly higher evaporation rate than
the ideal one.
5. Conclusions
A simplified model for bi-component droplet heating and evaporation is suggested. This model takes into account droplet heating
by convection from the ambient gas, the distribution of tempera-
ture inside the droplet, diffusion of liquid species inside the droplet, droplet swelling or contraction due to changing average
temperature, effects of the non-unity activity coefficient (ideal
and non-ideal models) and the effects of interaction between moving droplets due to the finite distance parameter (ratio of the average distance between droplets to their diameter). The effects of
recirculation in the moving droplets on heat and species mass diffusion within them are taken into account using the effective thermal conductivity and effective diffusivity models. The previously
obtained analytical solution of the transient heat conduction equation has been incorporated in the numerical code alongside the original analytical solution of the species diffusion equation inside
droplets.
The predicted time evolutions of surface, average and central
droplet temperatures have been compared with the results of direct measurements of droplet average temperatures for the case
Please cite this article in press as: S.S. Sazhin et al., A simplified model for bi-component droplet heating and evaporation, Int. J. Heat Mass Transfer (2010),
doi:10.1016/j.ijheatmasstransfer.2010.06.044
9
S.S. Sazhin et al. / International Journal of Heat and Mass Transfer xxx (2010) xxx–xxx
40
75% ethanol + 25% acetone
T (C)
35
30
Experiment
non-ideal-Ts
Ts (non-ideal)
25
Tav (non-ideal)
non-ideal-Tav
Numerical results of [15] - Ts (non-ideal)
Numerical results of [15] - Tav (non-ideal)
20
0
1
2
3
4
Time (ms)
5
6
7
Fig. 7. The same as Fig. 6 but including additional results of calculations of Ts and Tav (based on the vortex model, as reported in [15]) instead of the results referring to Tc.
68.8
75% ethanol + 25% acetone
68.6
68.4
Rd (µm)
68.2
68
67.8
New model (non-ideal)
67.6
New model (ideal)
67.4
Numerical results of [15] (non-ideal)
67.2
Numerical results of [15] (ideal)
67
0
1
2
3
4
Time (ms)
5
6
7
Fig. 8. The plots of Rd versus time predicted by the ideal (ci = 1) and non-ideal, simplified and vortex models for the same droplets as in Figs. 6 and 7.
of various mixtures of ethanol and acetone. There is a general
agreement between the predicted and observed average temperatures in the case of pure acetone and acetone-rich mixtures. In the
case of ethanol, 50% ethanol–50% acetone and 75% ethanol–25%
acetone mixture droplets the predicted average droplet temperature was several degrees (up to 5 K) higher compared with the observed one. It was shown that the temperatures predicted by the
simplified model and the earlier reported vortex model were reasonably close. Also, the temperatures predicted by the ideal and
non-ideal models differ by not more that several degrees. This
can justify the application of the simplified model with the activity
coefficient equal to 1 for the interpretation of the time evolution of
temperatures measured with similar errors.
Acknowledgements
The authors are grateful to the European Regional Development
Fund Franco-British INTERREG IVA (Project C5, Reference 4005)
and the Royal Society (UK) (Project JP71790) for the financial support of this project.
Appendix A. Transport coefficients for mixtures
A.1. Liquid thermal conductivity for a mixture
A number of approximate formulae for thermal conductivity in
multi-component liquids are reviewed in [31]. We use the most
simple Filippov equation, valid for binary mixtures:
kl ¼ Y 1 kl1 þ Y 2 kl2 0:72Y 1 Y 2 jkl1 kl2 j;
ð34Þ
where Y1 and Y2 are mass fractions for species 1 and 2, kl1 and kl2 are
thermal conductivities of species 1 and 2. The thermal conductivities of ethanol (kl,eth) and acetone (kl,ac) (in W/(m K)) were calculated based on the following approximations [38,39,30]:
Please cite this article in press as: S.S. Sazhin et al., A simplified model for bi-component droplet heating and evaporation, Int. J. Heat Mass Transfer (2010),
doi:10.1016/j.ijheatmasstransfer.2010.06.044
10
S.S. Sazhin et al. / International Journal of Heat and Mass Transfer xxx (2010) xxx–xxx
kl;eth ¼ 0:61572 0:24127 102 T þ 0:31333 105 T 2 ;
to
where T is in K everywhere.
D12 ¼ X 1 D012 þ X 2 D021 :
3
kl;ac ¼ 0:3133614225 0:8163 10
5
T þ 0:1 10
are close to 1, Formula (38) can be simplified
and a
Since both m
2
T ;
ð39Þ
D012
Among various approximations for
and
we have chosen
the Wilke–Chang approximation given by the following formula:
A.2. Liquid viscosity for a mixture
A number of approximate formulae for dynamic viscosity of
multi-component liquids (ll) are reviewed in [31]. It was shown
that the approach suggested by Grunberg and Nissan is the most
accurate and convenient in practical applications. In the case of a
bi-component liquid, this approach is based on the following
equation:
ln ll ¼ X 1 ln ll1 þ X 2 ln ll2 G12 X 1 X 2 ;
ð35Þ
where X1,2 are molar fractions of species 1 and 2. The values of G12
depend on the type of species involved [31]. In our case, following
[30], we use a simplified version of Eq. (35) assuming that G12 = 0.
In this case Eq. (35) is simplified to:
ll ¼ exp X 1 ln ll1 þ X 2 ln ll2 :
ð36Þ
The kinematic viscosity is obtained from Eq. (36) by dividing ll by
ql ¼
Y1
ql1
þ
Y2
ql2
1
ð37Þ
;
where q1,2 are densities of the components. The specific heat capacity of the mixture cl was estimated as
cl ¼ Y 1 cl1 þ Y 2 cl2 :
The values of ll (in Pa s), ql (in kg/m3) and cl (in J/(kg K)) for ethanol and acetone were estimated based on the following approximations [38,39,30]:
ll;eth ¼ 10ð
686:645:282
T
pffiffiffiffiffiffiffiffiffiffiffi
uM B T
lB V 0:6
A
ð40Þ
;
where D0AB is the mutual diffusion coefficient of solute A at very low
concentrations in solvent B, m2/s, MB is the molar mass of solvent B,
kg/kmol, T temperature in K, lB dynamic viscosity of solvent B, cP
(1 cP = 103 kg m1 s1), VA is the molar volume of solute A at its
normal boiling temperature, cm3/mol (this can be recalculated from
density if required, taking into account the molar mass of the substance), u is the associated factor of solvent B (u = 1.5 if B is ethanol, u = 1 if B is acetone).
VA (in cm3/mol) is calculated from the following equation (cf.
Eq. 11-4.2 of [31]):
V A ¼ ðrA =1:18Þ3 ;
ð41Þ
where the Lennard–Jones length rA in A is taken from Table A1.
When calculating D12 based on Eq. (38) the values of X1, X2 and
T, averaged over the whole volume of droplets, were used. These
parameters were updated at each timestep.
A.4. Gas diffusivity for a mixture
The value of Dv is estimated from the Wilke and Lee formula
[31]:
h
3:03 0:98=M 1=2
va
i
2
pM 1=2
v a rv a XD
ð107 ÞT 3=2
ð42Þ
;
Mv a ¼ 2½ð1=M v Þ þ ð1=M air Þ1
þ 0:223333 107 T 2 ;
ql;eth ¼ 1053:6 0:925 T;
ql;ac ¼ 986:53035880:6014966034T 0:2754046133103 T 2 ;
cl;eth ¼ 15039 130:53 T þ 0:4143 T 2 0:39583 T 3 ;
cl;ac ¼ 2165:234225 2:963 T þ 0:01 T 2 :
A.3. Liquid diffusivity for a mixture
As in the case of liquid thermal conductivity and liquid viscosity, various approximations for liquid diffusivity for a mixture were
discussed in [31]. One of the simplest approximations was given by
the Sanchez and Clifton formula:
ð38Þ
is to be found from one mixture datum
where the parameter m
is the thermodynamic factor defined as:
point, a
@a1
@a2
¼
;
@X 1 T;p @X 2 T;p
¼
7:4 1012
where Dv is in m2/s, T is temperature in K,
ll;ac ¼ 0:3183313525 102 0:1629735179 104 T
¼
a
D0AB
Dv ¼
Þ;
h
i
Þ;
þm
a
D12 ¼ X 1 D012 þ X 2 D021 ð1 m
D021
a1,2 are the activities of species 1 and 2, D012 and D021 are diffusivities
of dilute solute 1 in solvent 2, and dilute solute 2 in solvent 1.
Mv and Mair are molar masses of vapour and air, respectively, p is
pressure in bar, rva = (rv + ra)/2, rv and ra are characteristic Lennard–Jones lengths for vapour and air, respectively, measured in
Angstrom (A), XD is the function of the normalised temperature
pffiffiffiffiffiffiffiffiffiffiffi
T* = kBT/eva given by Eq. (B6) in [4], ev a ¼ ev ea , ev and ea are characteristic Lennard–Jones energies for vapour and air, respectively, kB
is the Boltzmann constant. The values of rv,a and ev,a/kB for air, ethanol and acetone are given in Table A1 [31]:
In the case of the mixture of vapour components, all input
parameters (molar masses and Lennard–Jones parameters) are calculated as molar averaged, taking into account their relative molar
concentrations Xeth/(Xeth + Xacet) and Xacet/(Xeth + Xacet).
A.5. Gas thermal conductivity and viscosity
The thermal conductivity (kv) and dynamic viscosity (lv) for the
mixture of air, ethanol vapour and acetone vapour was calculated
following the procedure described in [4] (see their Appendix B).
Table A1
The Lennard–Jones lengths and energies for air, ethanol and acetone.
ev,a/kB (K)
Substance
rv,a (A)
Air
3.711
78.6
Ethanol
4.530
362.6
Acetone
4.600
560.2
Please cite this article in press as: S.S. Sazhin et al., A simplified model for bi-component droplet heating and evaporation, Int. J. Heat Mass Transfer (2010),
doi:10.1016/j.ijheatmasstransfer.2010.06.044
S.S. Sazhin et al. / International Journal of Heat and Mass Transfer xxx (2010) xxx–xxx
The following approximations for kv (in W/(m K)) and lv (in Pa s)
for ethanol and acetone vapour were used [38,39,30]:
kv ;eth ¼ 1:8037 103 þ 7:419 106 T þ 1:1536 107 T 2 ;
20:3428433650106 T 2
T
kv ;ac ¼ 0:01143468
;
273:15
lv ;eth ¼ 0:19757 106 þ 0:29211 107 T;
lv ;ac ¼ 0:1579884447 105 þ 0:3090958433 107 T
0:449329 1011 T 2 :
Appendix B. Activity coefficient for the ethanol–acetone
mixtures
The activity coefficient (ci) for any liquid mixture is related to
the excess of Gibbs free energy per unit mole (GE) by the following
formula [28]:
ln ci ¼
E
1 @ nl;total G
;
Ru T
@nli
ð43Þ
where Ru is the universal gas constant, T is temperature in K,
P
nl;total ¼ i nli ; nli is the molar concentration of the ith component
in the liquid phase.
The following approximation for the ethanol (i = 1)–acetone
(i = 2) mixture, based on fitting experimental data, was suggested
[29]:
"
#
nl1 nl2
anl2
bnl1
dnl1 nl2
nl;total G ¼ Ru T
þ
;
nl1 þ nl2 nl1 þ nl2 nl1 þ nl2 ðnl1 þ nl2 Þ2
E
ð44Þ
where
546:3
0:9897;
T
543:3
0:9483;
b¼
T
15:63
d¼
þ 0:0759:
T
a¼
Having substituted (44) into (43) the following expressions for
c1 = ceth and c2 = cac are obtained:
h
h
ii 9
=
c1 ¼ exp X 2l2 a þ 2ðb a dÞX l1 þ 3dX 2l1 >
h
h
ii
:
;
c2 ¼ exp X 2l1 b þ 2ða b dÞX l2 þ 3dX 2l2 >
ð45Þ
Expressions (45) were used in Eq. (29).
References
[1] S.S. Sazhin, P.A. Krutitskii, A conduction model for transient heating of fuel
droplets, in: H.G.W. Begehre, R.P. Gilbert, M.W. Wong (Eds.), Progress in
analysis, Proceedings of the 3rd International ISAAC (International Society for
Analysis, Applications and Computations) Congress (August 20–25, 2001,
Berlin), vol. II, World Scientific, Singapore, 2003, pp. 1231–1239.
[2] S.S. Sazhin, P.A. Krutitskii, W.A. Abdelghaffar, S.V. Mikhalovsky, S.T. Meikle,
M.R. Heikal, Transient heating of diesel fuel droplets, Int. J. Heat Mass Transfer
47 (2004) 3327–3340.
[3] S.S. Sazhin, W.A. Abdelghaffar, P.A. Krutitskii, E.M. Sazhina, M.R. Heikal, New
approaches to numerical modelling of droplet transient heating and
evaporation, Int. J. Heat Mass Transfer 48 (2005) 4215–4228.
[4] S.S. Sazhin, T. Kristyadi, W.A. Abdelghaffar, M.R. Heikal, Models for fuel droplet
heating and evaporation: comparative analysis, Fuel 85 (2006) 1613–1630.
[5] S.S. Sazhin, Advanced models of fuel droplet heating and evaporation, Prog.
Energy Combust. Sci. 32 (2006) 162–214.
[6] C. Bertoli, M. Migliaccio, A finite conductivity model for diesel spray
evaporation computations, Int. J. Heat Fluid Flow 20 (1999) 552–561.
11
[7] C. Maqua, G. Castanet, F. Grisch, F. Lemoine, T. Kristyadi, S.S. Sazhin,
Monodisperse droplet heating and evaporation: experimental study and
modelling, Int. J. Heat Mass Transfer 51 (2008) 3932–3945.
[8] T. Kristyadi, V. Deprédurand, G. Castanet, F. Lemoine, S.S. Sazhin, A. Elwardany,
E.M. Sazhina, M.R. Heikal, Monodisperse monocomponent fuel droplet heating
and evaporation, Fuel 89 (2010), in press.
[9] G.M. Faeth, Evaporation and combustion of sprays, Prog. Energy Combust. Sci.
9 (1983) 1–76.
[10] A.Y. Tong, W.A. Sirignano, Multicomponent transient droplet vaporization
with internal circulation: integral equation formulation, Numer. Heat Transfer
10 (1986) 253–278.
[11] G. Continillo, W.A. Sirignano, Unsteady, spherically-symmetric flame
propagation through multicomponent fuel spray clouds, in: G. Angelino, L.
De Luca, W.A. Sirignano (Eds.), Modern Research Topics in Aerospace
Propulsion, Springer-Verlag, 1991, pp. 173–198.
[12] M. Klingsporn, U. Renz, Vaporization of a binary unsteady spray at high
temperature and high pressure, Int. J. Heat Mass Transfer 37 (Suppl. 1) (1994)
265–272.
[13] P.L.C. Lage, C.M. Hackenberg, R.H. Rangel, Nonideal vaporization of dilating
binary droplets with radiation absorption, Combust. Flame 101 (1995) 36–
44.
[14] D.J. Torres, P.J. O’Rourke, A.A. Amsden, Efficient multi-component fuel
algorithm, Combust. Theory Model. 7 (2003) 67–86.
[15] C. Maqua, G. Castanet, F. Lemoine, Bi-component droplets evaporation:
temperature measurements and modelling, Fuel 87 (2008) 2932–2942.
[16] J. Tamim, W.L.H. Hallett, Continuous thermodynamics model for
multicomponent vaporization, Chem. Eng. Sci. 50 (1995) 2933–2942.
[17] A.M. Lippert, R.D. Reitz, Modelling of Multicomponent Fuels using Continuous
Distributions with Application to Droplet Evaporation and Sprays, SAE
Technical Paper 972882, 1997.
[18] G.-S. Zhu, R.D. Reitz, A model for high-pressure vaporization of droplets of
complex liquid mixture using continuous thermodynamics, Int. J. Heat Mass
Transfer 45 (2002) 495–507.
[19] M. Burger, R. Schmehl, K. Prommersberger, O. Schäfer, R. Koch, S. Wittig,
Droplet evaporation modelling by the distillation curve model: accounting for
kerosene fuel and elevated pressures, Int. J. Heat Mass Transfer 46 (2003)
4403–4412.
[20] C. Laurent, Dévelopment et validation de modèles d’évaporation multicomposant, Thèse de l’Institut Supériour de l’Aéronautique et de l’Espace de
Toulouse, 2008.
[21] H.S. Carslaw, J.C. Jaeger, Conduction of Heat in Solids, Clarendon Press, Oxford,
1986.
[22] E.M. Kartashov, Analytical Methods in Heat Transfer Theory in Solids,
Vysshaya Shkola, Moscow, 2001 (in Russian).
[23] B. Abramzon, W.A. Sirignano, Droplet vaporization model for spray
combustion calculations, Int. J. Heat Mass Transfer 32 (1989) 1605–1618.
[24] B. Abramzon, S. Sazhin, Droplet vaporization model in the presence of thermal
radiation, Int. J. Heat Mass Transfer 48 (2005) 1868–1873.
[25] B. Abramzon, S. Sazhin, Convective vaporization of fuel droplets with thermal
radiation absorption, Fuel 85 (2006) 32–46.
[26] W.A. Sirignano, Fluid Dynamics and Transport of Droplets and Sprays,
Cambridge University Press, Cambridge, 1999.
[27] J.-P. Delplanque, R.H. Rangel, W.A. Sirignano, Liquid-waste incineration in a
parallel-stream configuration: effect of auxiliary fuel, Progress in Aeronautics
Astronautics 132 (1991) 164–184.
[28] P. Atkins, J. de Paula, Atkins’ Physical Chemistry, seventh ed., Oxford University
Press, 2002.
[29] S.W. Campbell, R.A. Willsak, G. Thodos, Vapor–liquid equilibrium
measurements for the ethanol–acetone system at 372.7, 397.7 and 422.6 K,
J. Chem. Eng. Data 32 (1987) 357–362.
[30] C. Maqua, Contribution à la compréhension de l’évaporation de gouttes de
combustible bi-composant à l’aide de méthodes optique, Thèse, NancyUniversité, 2007.
[31] B.E. Polling, J.M. Prausnitz, J.P. O’Connell, The Properties of Gases and Liquids,
McGraw-Hill, 2000.
[32] C. Maqua, G. Castanet, N. Doué, G. Lavergne, F. Lemoine, Temperature
measurements of binary droplets using three color laser-induced
fluorescence, Exp. Fluids 40 (2006) 786797.
[33] G. Castanet, A. Delconte, F. Lemoine, L. Mees, G. Grehan, Evaluation of
temperature gradients within combusting droplets in linear stream using two
colors laser-induced fluorescence, Exp. Fluids 39 (2005) 431–440.
[34] K.M. Watson, Thermodynamics of the liquid state, Ind. Eng. Chem. 35 (1943)
398–406.
[35] G. Castanet, C. Maqua, M. Orain, F. Grisch, F. Lemoine, Investigation of heat and
mass transfer between the two phases of an evaporating droplet stream using
laser-induced fluorescence techniques: comparison with modelling, Int. J. Heat
Mass Transfer 50 (2007) 3670–3683.
[36] L.A. Dombrovsky, S.S. Sazhin, A parabolic temperature profile model for
heating of droplets, ASME J. Heat Transfer 125 (2003) 535–537.
[37] L.A. Dombrovsky, S.S. Sazhin, A simplified non-isothermal model for droplet
heating and evaporation, Int. Commun. Heat Mass Transfer 30 (2003) 787–
796.
[38] R.C. Reid, J.M. Prausnitz, T.K. Sherwood, The Properties of Gases and Liquids,
third ed., Mc Graw Hill, 1977.
[39] R.C. Reid, J.M. Prausnitz, B.E. Poling, The Properties of Gases and Liquids, fourth
ed., Mc Graw Hill, 1987.
Please cite this article in press as: S.S. Sazhin et al., A simplified model for bi-component droplet heating and evaporation, Int. J. Heat Mass Transfer (2010),
doi:10.1016/j.ijheatmasstransfer.2010.06.044