Modeling of heat and mass transfer processes in phase transformation cycle of sprayed water into gas: 3. Energy and thermal states analysis of slipping droplet in a humid air flow.

Miliauskas, G. ; Maziukiene, M. ; Ramanauskas, V. 等

Nomenclature [B.sub.T]--Spalding transfer number; D--mass diffusivity, [m.sup.2]/s;

[F.sub.o]--Fourier number; g--evaporation velocity, kg/s;

[k.sup.-.sub.c]--effective thermal conductivity parameter; L--latent heat of evaporation, J/kg; m--vapour mass flux, kg/([m.sup.2]s);

Nu--Nusselt number; p--pressure, Pa; P--symbol of free parameter in heat-mass transfer; [bar.P]--dimensionless parameter in heat-mass transfer; [??]--average parameter in heat-mass transfer; Pe--Peclet number; Pr--Prandtl number; q--heat flux, W/[m.sup.2]; r--radial coordinate, m; Re Reynolds number; S--area, [m.sup.2]; T--temperature, K; [eta]--non-dimensional radial coordinate; [lambda]--thermal conductivity, W/(m K); [mu]- molecular mass, kg/kmol; [rho]- density, kg/[m.sup.3]; [tau]- time, s;

subscripts

c--convection; "c"--convective heating; "c + r"--convective-radiative heating; C--droplet centre; co--condensation; e--equilibrium evaporation; f--phase change; g--gas; i--time index in a digital scheme; it--number of iteration; IT--index of closing iteration; I--index of control time; j--index of radial coordinate; J--index of droplet surface; k--conduction; " k"--conductive heating; l--liquid; m--mass average; r--radiation; rt--dew point; R--droplet surface; v--vapor; vg--vapor-gas mixture; 0--initial state; [SIGMA]--total; [infinity]--far from a droplet;

superscripts

+--external side of a droplet surface; ---internal side of a droplet surface.

[B.sub.T]--Spalding transfer number; D--mass diffusivity, [m.sup.2]/s;

[F.sub.o]--Fourier number; g--evaporation velocity, kg/s;

[k.sup.-.sub.c]--effective thermal conductivity parameter; L--latent heat of evaporation, J/kg; m--vapour mass flux, kg/([m.sup.2]s);

Nu--Nusselt number; p--pressure, Pa; P--symbol of free parameter in heat-mass transfer; [bar.P]--dimensionless parameter in heat-mass transfer; [??]--average parameter in heat-mass transfer; Pe--Peclet number; Pr--Prandtl number; q--heat flux, W/[m.sup.2]; r--radial coordinate, m; Re Reynolds number; S--area, [m.sup.2]; T--temperature, K; [eta]--non-dimensional radial coordinate; [lambda]--thermal conductivity, W/(m K); [mu]- molecular mass, kg/kmol; [rho]- density, kg/[m.sup.3]; [tau]- time, s;

subscripts

c--convection; " c"--convective heating; "c + r"--convective-radiative heating; C--droplet centre; co--condensation; e--equilibrium evaporation; f--phase change; g--gas; i--time index in a digital scheme; it--number of iteration; IT--index of closing iteration; I--index of control time; j--index of radial coordinate; J--index of droplet surface; k--conduction; " k"--conductive heating; l--liquid; m--mass average; r--radiation; rt--dew point; R--droplet surface; v--vapor; vg--vapor-gas mixture; 0--initial state; [SIGMA]--total; [infinity]--far from a droplet;

superscripts

+--external side of a droplet surface;- --internal side of a droplet surface.

1. Introduction

Many thermal technologies are based on by droplets heat transfer and phase transformations. Therefore, researchers' attention for water heat transfer processes and for heat exchange and phase transformations of hydrocarbon and other liquids dispersed into droplets do not decrease [1-6]. A liquid droplets phase transformation cycle combines condensing, unsteady and equilibrium evaporation modes [tau] [equivalent to] 0 / [[tau].sub.co] / [[tau].sub.nf] / [[tau].sub.f] An intensive interaction of heat and mass transfer processes sets when a heat supplying conditions for droplets are rapidly changing. Droplet environment provides thermal energy by processes of heat exchange and liquid vapour condensation on the surface of the droplet. In condensation mode a supplying heat warmth the liquid in a droplet. At unsteady evaporation mode heat is being provided by heat transfer, where a part of heat vaporizes a liquid. At equilibrium evaporation mode all heat that is provided for the droplet participates in liquid evaporation. Peculiarities of interaction between phases determine the droplet thermal state, which has a reversible impact for a droplet energy state variation in phase transformation cycle.

Droplets phase transformation cycle is determined by sprayed liquid and carrier gas flow parameters. Liquid spray dispersity, a two-phase flow velocity and a temperature can be considered the essential parameters of sprayed liquid. Parameter of liquid vapour condensation in gas mixture is very important for heat and mass transfer intensity. A phase transformation cycle starts with condensing phase transformation mode, when sprayed pure liquid temperature is lower than a dew point temperature [[bar.T].sub.0] = [T.sub.rt]/[T.sub.0] > 1. In condensing phase mode a droplet is heated up and therefore grows rapidly. In condensing mode droplet heats up to dew point temperature and in unsteady evaporation mode droplets dispersity is defined. A droplet slipping in gas flow has a bright effect for condensing phase transformation mode that causes a longer term of liquid vapour condensation on a droplet surface by enhanced convective heating conditions [5]. A longer duration of condensing phase transformation means extended droplet surface heating to dew point temperature. A bit more liquid vapour condenses on the droplet surface therefore a droplet is larger at the end of condensing mode. Peculiarity of slipping droplet in condensing phase transformation mode can be based by intensive inner layers warming, compared with the case of a stable liquid, which is caused by fluid circulation. This requires assessment of droplet energy state in phase transformation cycle. It necessitate analysis of heat and mass flows on the droplet surface by conditions of a complex transfer processes. A problem of temperature gradient in a droplet rises when defining the intensity of the processes. This gradient is highly dependent from a droplet definition methodology. When a liquid is stable in a droplet, then temperature gradient in a droplet can be defined by integral type model that is convenient for numerical research [7]. Integral type model combines heat spread in droplet by conduction and radiation.

However, liquid can circulate in a droplet. These presumptions are made by non-isothermal in a droplet and by friction forces that arises on the droplet surface [1, 8, 9]. Archimedes forces rises in non-isothermal droplet. Their effect for liquid stability depends from droplet dispersity. In common occurring technologies such as liquid fuel burning, fireplaces or air conditioning systems droplets are quite small. Rising Archimedes forces in droplets are week and do not cause a fluid circulation [7]. A forced fluid circulation is arising in a slipping droplet. This circulation is described by Navier-Stokes equation system. It is impossible to solve this system analytically. Its numerical solution was analyzed in work [10]. This solution is receptive for machine counting and is difficult to apply for iterative scheme.

This article presents an energy state balance assessment methodology of sliding droplet in gas. In humid air slipping water droplet case is analyzed as well as energy and thermal states need for complex analysis is justified.

2. Research method

Energy state change is defined by connection of droplet internal and external heat transfer and by phase transformation processes that is ongoing on its surface. In case of condition [T.sub.g] > [T.sub.i] a droplet external heat transfer ensures continuous heat supply for a droplet in all phase transformation cycle. Heating intensity is described by density function of total heat flow on a droplet surface external side:

[q.sup.+.sub.[SIGMA],"c+r"]([tau]) = [q.sup.+.sub.r]([tau]) + [q.sup.+.sub.c]([tau]) (1)

In non-radiant environment for warming droplet case [q.sub.r] = 0, so [q.sup.+.sub.[SIGMA]] [equivalent to] [q.sup.+.sub.c]. In droplet phase transformation cycle a convective heating heat flow is described by expression:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (2)

where phase transformation impact for convective heating is taken into account by the Spalding transfer parameter [B.sub.T] [5]. A droplet slip velocity [DELTA][w.sub.l] = [w.sub.l] - [w.sub.g] defines a convection heating intensity, that determines Reynolds number Re = 2R [absolute value of [DELTA]w] [[rho].sub.g] /[[mu].sub.vg], where a gas density is chosen according to temperature [T.sub.d], and liquid vapor and gas mixture dynamic viscosity coefficient is selected by temperature according to "1/3" rule: [T.sub.vg] [equivalent to] [T.sub.R] + ([T.sub.g] - [T.sub.R])/ 3.

A combination of empirical model [11] for effective thermal conductivity and integral model of heat spread in droplet by conductivity-radiation [7] are adapted to define the heat exchange in a droplet. A heat flow that is leaded to the droplet is described by modified Fourier law of heat spread:

[q.sup.-.sub.c] ([tau]) = -[k.sup.-.sub.c] ([tau]) [[lambda].sub.1]([tau]) [partial derivative][T.sub.1](3)

A forced liquid circulation influence for heat spread in a droplet is evaluated by effective thermal conductivity parameter, which is defined by methodology [11] as function [k.sup.-.sub.c] = f ([Pe.sub.l]) of Peclet number. The temperature gradient in equation (3) is defined by methodology [7].

Heat flow which participates in phase transformations is calculated by water vapor flow density on the surface of a droplet:

[q.sup.-.sub.f]([tau]) = [m.sup._.sub.v]([tau])L([tau]) (4)

On a droplet surface, expressions (2)-(4) is applied to concretize an energy balance condition [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] for all modeled droplet heating cases. A balance condition requires matching of heat flows that flows in and flows from the surface of the droplet. In outspread form a balance condition becomes a transcendental system of algebraic and integral equations, because its solutions are unambiguous only then, when the droplet surface temperature function of time [T.sub.R]([tau]) is defined. This function is defined by iterative technique according to method [5]. However, the temperature distribution inside the droplet T(r < R, [tau]), can only be defined by case assumptions of non-circulating fluid. At numerical experiment attention is taken into account to phase transformation in technologies of heat recovery from humid combustion products. Water injection takes an important place for gas mixture cleaning and for additional irrigation [7]. Peculiarities of convective heating was highlight by modeling a water droplets of 278 K temperature at unsteady phase transformation mode in 500 K temperature humid gas flow ([[bar.p].sub.v,[infinity]] = 0.3, when p = 0.1 MPa). Droplets primary slipping is defined by Reynolds number values 0, 5, 10, 20, 40, and 80 for 2[R.sub.0] = 150 x [10.sup.-6] m droplet diameter, when a gas flow velocity is [w.sub.g] [equivalent to] 10 m/s. These [Re.sub.0] number values were ensured by primaries velocities [w.sub.l,0] m/s: 10, 11.12, 12.25, 14.49, 18.98 and 27.97, respectively. Fourier number application in time scale was carried out in numerical experiment. The main attention is laid on for unsteady phase transformation mode 0 / [Fo.sub.nf]. In numeric control scheme for condensing mode 21 control time moments are provided for time grid formation. Condensing mode duration [FO.sub.co] is defined [F.sub.O] [equivalent to] 0 / [FO.sub.co] (Fig. 1).

[FIGURE 1 OMITTED]

Unchanging individual time step [DELTA][Fo.sub.i] = [Fo.sub.co] /([I.sub.co] - 1), when [I.sub.co] = 21 was kept in each modeled cycle. A droplet surface temperature function [T.sub.R]([F.sub.o]) is defined according to (10) - (18) schemes in work [5]. Other parameters functions P(Fo) of heat and mass transfer are defined in parallel. At numerical experiment process a condition matching of heat flows that flows in and flows from a droplet surface is satisfied (Fig. 2.) by a strict requirement:

[absolute value of 1 - [q.sup.-.sub.c](Fo) + [q.sup.+.sub.f](Fo)/ [q.sup.+.sub.c] (Fo)] x 100% < [[delta].sub.R] [equivalent to] 0.01% (5)

[FIGURE 2 OMITTED]

In phase transformation cycle a gas surrounding gives the heat for a droplet by heat transfer (together with condensing heat or without a heat that is necessary for evaporation process):

[Q.sup.+.sub.J,co] + [Q.sup.+.sub.J,ng] + [Q.sup.+.sub.J,eg] = [Q.sup.+.sub.J] (6)

According to expression (6) described heat quantity can be defined according to heat flow density average values [[??].sup.-.sub.c,i] in a droplet that is distinguish at time intervals

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (7)

Because of physical nature change in phase transformation modes at heat flow in a droplet [7], peculiarities of phase transformation cycle modes should be taken into account when defining [[??].sup.-.sub.c,i]. At condensing phase transformation mode a gas surrounding gives thermal energy in heat transfer process together with vapour phase transformation heat that condensates on the surface of a droplet. At unsteady evaporation mode a convective heat warmth and evaporates a droplet. In equilibrium evaporation mode a convectional heat can only evaporates the droplet and enthalpy of cooling down droplet can participate in evaporation:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (9)

[FIGURE 3 OMITTED]

[FIGURE 4 OMITTED]

A droplet slipping in gas flow is decreasing because of the impact of resistance forces. The fastest droplet velocity change is in condensing phase transformation mode, while in evaporation mode a droplet slipping becomes negligible in term [approximately equal to] 5 x [[bar.F]o.sub.oco] (Fig. 3). Then a case of heating by conductivity approach in practice. Therefore at a droplet energy state analysis the main attention is given for condensing and unsteady evaporation modes.

3. Energy evaluation of unsteady phase transformation mode

A droplet energy state graphical interpretation, in freely chosen time Azi intervals at condensing and unsteady evaporation modes is reflected in Fig 4. In condensation phase transformation mode a droplet diameter increases from 2[R.sub.0] to 2[R.sub.co], therefore [R.sub.i] > [R.sub.i-1]. At evaporation mode a droplet diameter decreases from 2[R.sub.co] to zero, therefore [R.sub.i] < [R.sub.i-1].

[FIGURE 5 OMITTED]

[FIGURE 6 OMITTED]

[FIGURE 7 OMITTED]

[FIGURE 8 OMITTED]

At condensation phase transformation mode droplet convective heating intensity exceeds analogous solid particle heating case, when phase transformation does not expose on its surface. In unsteady evaporation mode a droplet heating intensity is weaker than solid particle, while in regime change moment, matches it black points (Fig. 5. Droplet slipping in gas flow has a bright impact for convective heat flow that warmth a droplet (Fig. 6), and also has an impact for phase transformation heat flow (Fig. 7). In unsteady phase transformation mode at a droplet liquid is heated up intensively (8 Fig. a), its surface temperature increases, and temperature difference [DELTA]T = [T.sub.g] - [T.sub.R] that reflects external heat exchange driving force, decreases (Fig. 9). Therefore convective heating intensity suffocates in unsteady condensing and evaporation phase transformation modes (Fig. 6).

[FIGURE 9 OMITTED]

In condensation mode a growing heat exchange area has a bright impact for convective heat flow density that is provided for a droplet (Fig. 10.). A warming liquid expansion and liquid vapour condensation on the surface causes a droplet surface area growth in condensation mode. Surface water evaporation compensates liquid expansion in unsteady evaporation mode. At balance moment of these factors extreme point forms in droplet surface variation curve therefore a droplet surface area start to decrease.

[FIGURE 10 OMITTED]

[FIGURE 11 OMITTED]

At equilibrium evaporation mode a droplet surface decreases linearly and satisfies a well known law [D.sup.2] [12]. For a more extensive evaluation of droplet phase transformation peculiarities needs to make a broader numerical research. Gas humidity and temperature, as well temperature of sprayed water and heat transfer of a droplet would be angular aspects. This is one of the future research topics. In unsteady phase transformation mode a heat that is provided for a droplet by convection from a gas surrounding is reflected in Fig. 11, while heat change that participates in phase transformation process is presented in Fig. 12.

[FIGURE 12 OMITTED]

Heat is being provided for a droplet by convection in whole phase transformation cycle, therefore a provided heat amount is consistently growing. In phase transformation balance assessment the evolved heat considered to be positive in condensation process, while a heat that is required for water evaporation considered to be negative. Therefore at the end of condensation mode in curve [Q.sup.+.sub.J,f] (Fo) of phase transformation heat dynamics an extreme point observes. Liquid evaporation starts from this point. Evaporation process continues due to a part of heat that is provided by convection. In condensation process, a released heat of phase transformation gets a bright effect from droplet slipping in gas: [Q.sup.+.sub.J,f] ([Re.sub.0], Fo [equivalent to] [Fo.sub.co]) x [10.sup.4], J is 3.195, 3.433, 4.159, 5.033, 5.649, 5.893, when Re is 0, 5, 10, 20, 40 and 80, respectively. The same heat amount is used at initial state of unsteady evaporation, that is defined by time moment [Q.sup.+.sub.J,f]([Re.sub.0], Fo) [equivalent to] 0 (Fig. 12 points x).

Therefore in condensation phase transformation mode a heat amount that is provided for a droplet is [Q.sup.+.sub.J,f] ([Re.sub.0] > [Fo.sub.co])j [Q.sup.+.sub.J,f] ([Re.sub.0] > [Fo.sub.co]). It is leaded to droplet by internal heat convection [Q.sup.-.sub.J,f] ([Re.sub.0], [Fo.sub.co]) x [10.sup.4], J: 4.985, 6.063, 7.546, 9.403, 10.933, and 11.962 when Re is 0, 5, 10, 20, 40 and 80, respectively. In a droplet water is heated from [T.sub.l,0] = 278 to [T.sub.l,m] (Fo [equivalent to] [Fo.sub.co]). Heat amount that is leaded to a droplet by internal heat convection (Fig. 13) is defined by condition [Q.sup.-.sub.J,f] (Fo) [equivalent to] [Q.sup.-.sub.J,f](fo [equivalent to] [Fo.sub.nf]), which enables to define unsteady phase transformation duration [Fo.sub.nf]([Re.sub.0]) from scientific assessment: 2.68, 1.739, 1.849, 1.999, 2.128, 2.186, when Re is 0, 5, 10, 20, 40 and 80, respectively (Fig. 13, b). In equilibrium evaporation mode a liquid in a droplet is not heated anymore. For droplet that is heated up by conductively this is reflected by condition [[bar.Q].sup.-.sub.J,f] (Fo >[Fo.sub.nf])= 1.

Internal heat convection in slipping droplet can still going on, but convection leads out heat from cooling down droplet to its surface and stimulates an evaporation process. This is reflected by the condition [[bar.Q].sup.-.sub.J,f] (Fo > [Fo.sub.nf])< 1.

For integrated droplet energy and thermal state change analysis it is required to define a droplet warming dynamics in all phase transformation cycle. A thermal state analysis of circulating water in a droplet is quite complicated and would extend this work. It would be one of the future research topics.

[FIGURE 13 OMITTED]

4. Conclusions

Water droplets slipping have a maximum impact for droplets heat transfer in humid gas at condensation phase transformation mode. A set analysis of energy state balance for different speeds moving droplets shows, that essential role falls far liquid circulation that is caused by friction forces that arises on the surface of the slipping droplet. For this reason, heat is leaded out to central layers of the droplet more intensively. Therefore droplet heats up gradually; a temperature gradient in surface layers is lower than at conductivity heating case and droplet surface heats up till a dew point longer. On the surface of a droplet slipping in humid gas liquid vapor condenses more intensively therefore a droplet diameter grows in condensation mode.

References

[1.] Sirignano, William A. 2014. Advances in droplet array combustion theory and modeling, Progress in Energy and Combustion Science 42: 54-86. http://dx.doi.org/10.1016/j.pecs.2014.01.002.

[2.] Feng, Y.; Kleinstreuer, C.; Rostami, A. 2015. Evapo-ration and condensation of multicomponent electronic cigarette droplets and conventional cigarette smoke particles in an idealized G3-G6 triple bifurcating unit, Journal of Aerosol Science 80: 58-74. http://dx.doi.org/10.1016/j.jaerosci.2014.11.002.

[3.] Irfan Javed; Seung Wook Baek; Khalid Waheed. 2014. Effects of dense concentrations of aluminum nanoparticles on the evaporation behavior of kerosene droplet at elevated temperatures: The phenomenon of microexplosion, Experimental Thermal and Fluid Science 56: 4433-4446. http://dx.doi.org/10.1016/j.expthermflusci.2013.11.006.

[4.] Wei An; Tong Zhu; NaiPing Gao. 2015. Accelerative iteration for coupled conductive-radiative heat transfer computation in semitransparent media, Int. J. Heat Mass Transfer 82: 503-509. http://dx. doi. org/ 10.1016/j.ij heatmasstransfer.2014.11.035.

[5.] Miliauskas, G; Maziukiene, M. 2015. Modeling of heat and mass transfer processes in phase transformation cycle of sprayed water into gas. 2. Phase transformation peculiarities of a droplet heated up by convection in humid gas flow, Mechanika 21(2): 123-128. http://dx.doi.org/10.5755/j01.mech.2L2.9774.

[6.] Chungen Yin. 2015. Modelling of heating and Evaporation of n-Heptane droplets: Towards a generic model for fuel droplet/particle conversion, Fuel 141: 64-73. http://dx.doi.org/10.1016/j.fuel.2014.10.031.

[7.] Miliauskas, G. 2003. Interaction of the transfer processes in semitransparent liquid droplets, Int. J. Heat Mass Transfer 46: 4119-41138. http://dx.doi.org/10.1016/S0017-9310(03)00231-X.

[8.] Wegener, M.; Paul, N.; Kraume, M. 2014. Fluid dynamics and mass transfer at single droplets in liquid/liquid systems, Int. J. Heat Mass Transfer 71: 475495. http://dx.doi.org/10.1016/j.ijheatmasstransfer.2013.12.024.

[9.] Deepak Kumar Mandal; Shamit Baksh. 2012. Internal circulation in a single droplet evaporating in a closed chamber, International Journal of Multiphase Flow 42: 42-51. http://dx.doi.org/10.1016/j.ijmultiphaseflow.2012.01.008.

[10.] Roland F. Engberg; Mirco Wegener; Eugeny Y. Kenig. 2014. The impact of Marangoni convection on fluid dynamics and mass transfer at deformable single rising droplets--A numerical study, Chemical Engineering Science 116: 208-222. http://dx.doi.org/10.1016/j.ces.2014.04.023.

[11.] Abramzon, B.; Sirignano, W. A. 1989. Droplet vaporization model for spray combustion calculations, Int. J. of Heat and Mass Transfer 32: 1605-1618. Available from internet: ftp://ftp.demec.ufpr.br/CFD/bibliografia/propulsao/abramzon_et_al_1989.pdf.

[12.] Faeth, GM. 1983. Evaporation and combustion of sprays, Prog. Energy combust. 9: 1-76. http://dx.doi.org/10.1016/0360-1285(83)90005-9.

Received April 29, 2015

Accepted June 23, 2015

G. Miliauskas *, M. Maziukiene *, V. Ramanauskas **

* Kaunas university of technology, K.Donelaicio 20, Kaunas, LT-44239, Lithuania, E-mail: [email protected]

** UAB "Enerstena", Raktazoliq 21, Kaunas, LT-52181, Lithuania, E-mail: [email protected]

http://dx.doi.org/10.5755/j01.mech.21.5.12167