Modeling of heat and mass transfer processes in phase transformation cycle of sprayed water into gas: 1. The calculation peculiarities of droplet phase transformation parameters/Silumokaitos ir mases pernasos procesu modeliavimas dujose ispurksto vandens faziniu virsmu cikle. 1. Laselio faziniu virsmu parametru apskaiciavimo savitumai.

Miliauskas, G. ; Maziukiene, M.

Nomenclature

a--thermal diffusivity, [m.sup.2]/s; B--Spalding transfer number; [c.sub.p]--mass specific heat, J/(kg K); D--mass diffusivity, [m.sup.2]/s; Fo--Fourier number; G--mass flow rate, kg/s; L--latent heat of evaporation, J/kg; m--vapour mass flux, kg/([m.sup.2]s); p--pressure, Pa; [bar.P]--symbol of free parameter in heat-mass transfer; [bar.P]-droplet dimensionless parameter in heat-mass transfer; q--heat flux, W/[m.sup.2]; r-radial coordinate, m; [R.sub.[mu]]--universal gas constant J/(kmol K); 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.

[Subscripts.bar]: C--droplet centre; c--convective; e equilibrium evaporation; f--phase change; g--gas; i--time index in a digital scheme; it--number of iteration; I--index of control time; j--index of radial coordinate; J--index of droplet surface; k--conduction; l--liquid; m--mass average; r--radiation; R--droplet surface; sor--source; v-vapor; 0--initial state; [infinity]--far from a droplet; [SIGMA]-total.

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

1. Introduction

Water is widely used in industry and energy sector. Water spraying is the first step towards a modern micro-systems, by which technologies of transfer processes are based on. In these technologies contact area is strongly developed between one-piece carrier (often gas) and discretionary (often liquid droplets or solid particles) medium. Fast change of heat and mass transfer between phases can be organized. The latter is defined by intensity of transfer processes. Therefore, a management of water sprayed technologies requires a deep-knowledge in droplets and in gas two-phase flow transformation processes. Heat and mass transfer processes are closely linked and interact with each other. This interaction is operated by conditions of heat and mass transfer [1]. Primary parameters of gas and water are important to define these conditions, because parameters concretizes droplets phase transformation cycle. Droplet phase transformation cycle of pure liquid includes condensation, unsteady evaporation and equilibrium evaporation modes [2]. For a condensing phase transformation mode proceeding, the most important water vapor component must be ([[bar.p].sub.v,[infinity]] = [p.sub.v,[infinity]] / p > 0) in carrier gas mixture, and a droplet surface temperature must be lower than a dew point temperature ([[bar.T].sub.R] [equivalent to] [T.sub.R] / [T.sub.rt] < 1). At condensing phase transformation mode, droplet is heated up by the warmth from the gas in heat exchange process (case [T.sub.d] > [T.sub.R]) and by the warmth of phase transitions, which is released at water vapor condensation time. Therefore, droplet surface layers heats up rapidly. When they heat up till temperature of dew point, phase transformation regime on a surface of the droplet changes into unsteady evaporation. This happens at time moment, when parameter is [[bar.T].sub.R] = 1. At unsteady evaporation mode a part of warmth, which is provided for a droplet by a heat transfer, is used for heating. The other part of warmth participates in process of water evaporation. A droplet which evaporates unsteady heats up till temperature [T.sub.e], which describes equilibrium evaporation conditions. Equilibrium evaporation is comprehensible phase transformation case, when all warmth, which is given to a droplet at heat transfer process, evaporates water. Temperature of equilibrium evaporation mode is defined by gas parameters [T.sub.d] and [[bar.p].sub.v,[infinity]] as well as droplet heating conditions [3]. In condensating and unsteady evaporation modes a droplet is non-isothermal, and its thermal state at equilibrium evaporation mode is influenced by peculiarities of heat transfer between gas flux and a droplet [4].

A liquid vapor flux intensity factor, is very important in sprayed liquid technologies. This factor describes efficiency of phase transformation. At liquid fuel case vapor flux that spreads from droplets defines a combustion process. Peculiarities of water sprayed technology defines a vapor flux that is required for a water case. Cases of air conditioning, sudden gas cooling and heat recovering of phase transformation from removed smoke can be distinguish. Droplet evaporation process defines efficiency of irrigation and cooling, while a vapor condensation process in smoke determines a heat recovery. Phase transformations that occur on droplets surface raise a structural change of two-phase flow. In order to define them, dynamics of droplets phase transformation schould be taken into account, linking warming and weight change of droplet.

This article presents a defining specifics of phase transformation in droplet lifetime cycle, when droplet warms in different gas conditions.

2. Research method

Droplet mass dynamics is defined by vapor flux function on droplet surface:

dM([tau])/d[tau] = -[g.sup.+.sub.v]([tau]) (i)

In Eq. (1) vapor flux is positive at liquid evaporation mode, when droplet mass decreases. In condensing phase transformation mode a vapor flux, that flows condensates on the surface of the droplet, considered to be negative, then mass of droplet grows. Liquid vapor flux density can be described by analytical model [5]:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (2)

Then vapor flux function can be defined spherically:

[g.sup.+.sub.v]([tau]) = 4[pi][R.sup.2] ([tau])[m.sup.+.sub.v]([tau]). (3)

Vapor flux function density [m.sup.+.sub.v]([tau]) calculation according Eq. (2) is directly related with functions definition and describes the change of droplet surface temperature [T.sub.R]([tau]) and its dimension R(t) in a droplet phase transformation cycle. Droplet warming and ongoing phase transformation processes on its surface are closely related and influence one to another. This impact is indirect and revels through function [T.sub.R]([tau]), R([tau]) as well [m.sup.+.sub.v]([tau]) influence for heat fluxes dynamics on the surface of a droplet.

In droplet phase transformation process, heat flux intensity is defined by a vapor flux density on the droplet surface: [q.sup.+.sub.f]([tau]) = [m.sup.+.sub.v]([tau])L([tau]). Surrounded energy impact for a droplet is defined by compound heat transfer of a total heat flux density [q.sup.+.sub.[SIGMA]]([tau]) = [q.sup.+.sub.c]([tau]) [q.sup.+.sub.r]([tau]) in thermal technology at a common case, when [T.sub.d] > [T.sub.R]. Total heat flux of convective and radiative components impact for a droplet thermal state have a different mechanisms: convection warmth is given for a surface of a droplet and in radiation heat flux is absorbed in semitransparent droplet [6-8]. The radiation heat flux directly warmth droplet inner layers, but practically has no direct effect for a droplet surface temperature function [T.sub.R]([tau]). A convective heat flux impact for a droplet surface temperature is direct, but for internal thermal state it is only based on by a partial spreading in a droplet by inner heat transfer. The rest of external convection heat flux participates directly in liquid evaporation process on the surface of a droplet. Heat radiation, that is absorbed by a droplet, can only participate in evaporation process when is leaded out to the surface on the droplet by internal heat transfer. This is possible when a negative temperature gradient field forms [7].

Heat fluxes, that interact on the surface of a droplet, defines function [T.sub.R]([tau]). Its description is a key for heat and mass transfer task, which is known as a "droplet" problem. For that droplet internal and external heat transfer and heat fluxes [q.sup.-.sub.[SIGMA]]([tau]), [q.sup.+.sub.[SIGMA]]([tau]) as well [q.sup.+.sub.f]([tau]), that describes intensity of phase transformation are combined by energy flow balance on the droplet surface according to request of transfer processes in quasi-steady state:

[[??].sup.+.sub.[SIGMA]]([tau]) + [[??].sup.+.sub.[SIGMA]]([tau]) + [[??].sup.+.sub.f] = 0. (4)

Eq. (4) requires matching of heat fluxes that flows in and flows out from a droplet surface. But Eq. (4) is a formal expression, therefore it is necessary to specify expression in order to make a "droplet" task numerical solution algorithm. This requires analysis of four droplet heat and mass transfer tasks: 1) case of compound heat exchange in semitransparent droplet; 2) spectral radiation absorption in semitransparent droplet; 3) phase transformation that are set on the surface of a droplet; 4) droplet heating by external convection in phase transitions movement conditions. The first problem solution are functions of unsteady temperature field T(r, [tau]) and total heat flux [q.sup.-.sub.[SIGMA]](r), respectively; for the second--a local radiation flux function [q.sub.r](r) in a droplet at defined time moments [tau]; the third--vapor flux on the surface of a droplet, vapor flux density and heat flux of phase transformations functions [g.sup.+.sub.v]([tau]), [m.sup.+.sub.v]([tau]), [q.sup.+.sub.f]([tau]); the fourth--a heat flux density function of external convection [q.sup.+.sub.c](t). Each solution of discussed problem requires for already defined answers of remaining problems. This is taken into account during unambiguousity boundary conditions formulating. This produces only conditionally independent solutions of separate tasks. At case of combined heat transfer by conduction-radiation in a droplet (4) expression integral mathematical model is developed in [7] work. It is solved numerically. Initial conditions are defined by parameters [T.sub.i,0], [R.sub.0], [T.sub.d], p, [p.sub.v,[infinity]], heating time [[tau].sub.I] of droplet is provided, and in time change interval 0 - [[tau].sub.I] index of control time I is selected that meets condition:

[I.summation over (i = 2)]([[tau].sub.i] - [T.sub.i-1]) [equivalent to] [[tau].sub.I]. (5)

For each time moment Ti in iterative cycle it [equivalent to] 1 - IT a droplet surface temperature [T.sub.R,i] [equivalent to] [T.sub.R,i,it [equivalent to] IT] is selected. Whole number IT defines the end of iterative cycle. This number provides (4) condition requiring accuracy of the balance [[delta].sub.q] Iterative cycle is exercised in fastest descent method [2-4, 7, 9]. Rising a high requiring accuracy [[delta].sub.q] < 0.1 %, a droplet surface temperature changes of hundredth row must be taken into account [9]. Condition (4) sensitivity for a droplet surface temperature is defined by Eq. (2), which is described as vapor flux logarithmic dependence from generated saturated water vapor pressure on the surface of droplet (in case of condensation from saturated vapor pressure). The latter is a droplet surface function [p.sub.v,R] [equivalent to] [p.sub.s]([T.sub.R]). Therefore, water droplet phase transformation cycle is modeled by applying a strict Gerry empirical correlation between saturated water vapor and temperature [10].

Defending heat fluxes in Eq. (4) at selected time moments [T.sub.i], it is necessary to dispose with dimension of droplet size. At spherical assumptions case, this is droplet diameter 2[R.sub.i] defined by droplet rays [R.sub.i]. An inconvenience of numerical scheme formation is that droplet volume, in rigorous assessment, can only be defined when iterative cycle ends by modified expression (6):

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (6)

A problem of droplet dimension definition in common iteration is solved by making assumptions. It is necessary to avoid potential iterative cycle instability problem. The latter is related with hardly programmed by (2) expression vapor flux values, when random temperatures [T.sub.R,i,it] it are selected at each iterative cycle initial state. When [R.sub.i=1] = [R.sub.0], and in i > 1 cases [R.sub.i-1] are defined in previous iterative cycles, it is popular to apply a droplet dimension stability assumption [R.sub.i,it] = [R.sub.i-1]. Then a change of droplet diameter is calculated just before next iterative cycle When droplet surface temperature [T.sub.R,i,it] and radius [R.sub.i,it] parameter are selected for ordinary iteration, then heat fluxes in Eq. (4) are described by defined functions [3, 7, 11]:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (7)

Phase transformation heat fluxes and external convection intensity of heat fluxes are inversely proportional to radius of spherical droplet, when translucent droplet decreases a radiation absorption suffocate in them [3, 7, 12]. Therefore, in "droplet" studies for assessment of assumption of droplet dimension stability requires for comprehensive numerical research of phase transformation cycle and systemic evaluation of results. Droplet heat transfer and phase transformation parameters is appropriate to combine in groups of thermal [P.sub.T], energetic [P.sub.q], dynamic [P.sub.d] and phase transitions [P.sub.f] [13]. Different impact of assumption [R.sub.i,it] [equivalent to] [R.sub.i-1] is expected for each group parameters. In order to highlight this impact a consistent way of evaluation is suppositional, when [R.sub.i,it] [equivalent to] [R.sub.i,it-i] assumption is applied in numeric scheme calculating for selected group parameters, and in other groups a condition [R.sub.i,it] [equivalent to] [R.sub.i-1] remains the same. Phase transformation parameters are often important in thermal technologies of water spraying. A liquid vapor flux on the surface of the droplet [g.sup.+.sub.v] and its density [m.sup.+.sub.v] are assigned for them, as well they defining a droplet mass M, volume V and radius R dynamics in droplet phase transformation cycle. A ray dynamics numerical scheme is formed in Eq. (6) basis:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (8)

where:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (9)

here [R.sub.i=1] = [R.sub.0] when parameter [rho] is non-isothermal droplet mass average density, which defines a droplet mass average temperature function [rho] [equivalent to] [[rho].sub.l]([T.sub.m]), when [R.sub.i=1] = [R.sub.0]:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (10)

Vapor flux density on the surface of a droplet is calculated according to scheme:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (11)

The droplet mass is specified at each iterative cycle ends it = 1 + IT according to scheme:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (12)

where [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII].

Evaluation of assumptions [R.sub.i,it] [equivalent to] [R.sub.i-1] and

[R.sub.i,it] [equivalent to] [R.sub.i,it-1] operation for calculated parameters [P.sub.f] in unsteady phase transformation mode is given below.

3. Results and discussion of numerical research of droplets unsteady phase transformation

Condensing and unsteady evaporation modes were modelled according [2] methods to highlight functions, [g.sup.+.sub.v]([tau])-water vapour, [m.sup.+.sub.v]([tau])--density, M ([tau])--droplet mass, V([tau]) and R([tau])--dimension, that describes calculation peculiarities of phase transformation of sprayed water. Air temperature was [T.sub.g] = 500 K, pressure p = 0.1 MPa, air humidity [p.sub.g,[infinity]] / p = 0.3, diameter of sprayed droplet was 2[R.sub.0] = 150 x [10.sup.-6] m, and temperature [T.sub.0] = 278 K. It is considered, that [G.sub.l,0]/[G.sub.d,0] [right arrow] 0, therefore heat transfer and droplet phase transformation do not change gas flow parameters. Assuming, that air is provided for droplets by conductivity ("k" heat transfer case), (4) expression is concretized for conditions Nu [equivalent to] 2 and [q.sub.r] [equivalent to] 0. Droplet that are heating by conduction--radiation (combined case "k+r") radiation flux density [q.sub.r](r), unsteady temperature field T(r,[tau]) and its gradient gradT(r [equivalent to] R, [tau]) are described by integral expressions according [7] methodology. Iterative numerical scheme grid is formed by applying universal droplet radial coordinate [eta] = r/R, which in case [eta] = 0 indicates a spherical droplet centre, and in case [eta] = 1 defines universal unitary droplet radius. By selected number of control cross-sections J = 61, a droplet single dimension liner splitting, satisfies the condition:

[J.summation over (j = 2)]([[eta].sub.i] - [[eta].sub.i-1]) [equivalent to] 1. (13)

Because of valid condition [T.sub.0]/[T.sub.rt] < 1, therefore water droplet phase transformation cycle begins at condensing mode, which will be consistently replaced by unsteady and equilibrium evaporation modes. The time scale is formed for phase transformation modes in numeric schemes cycle [tau] [equivalent to] 0 - [T.sub.k0] - [T.sub.nf] - [T.sub.f]. They are taken into account individually, but equivalent modes importance in cycle treatment principle is applied. Time scale is expressed by Fourier number [F.sub.0] = ([a.sub.0] / [R.sup.2.sub.0]) x z and unsteady phase transformation cycle is [F.sub.0] [equivalent to] 0 - [F.sub.0ko] - [Fo.sub.nf] - [Fo.sub.f] modelled. Heat transfer cycle in "k" case is considered to be supporting for affect evaluation of droplets heat transfer conditions

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (14)

It is important, because sprayed water temperature [T.sub.0] and gas parameters [T.sub.d] as well [[bar.p].sub.v,[infinity]] in defined cases are independent from a droplet dispersity. Here impact of Knudsen layer for heat transfer and phase transformation can be denied [9]. [Fo.sub.fv,"k"] criteria determinates time of phase transformation mode cycle and serves for creation of unit time phase transformation cycle [[bar.F].sub.O"k"] = 0 - 1 - 2 - 3, where [[bar.F].sub.O"k"], = f ([Fo.sub.], [Fo.sub.fv,"k"]) [2]. Modeled phase transition mode values of [[bar.F].sub.Oi,"k"], = (i - l)/(l - 1) are defined by whole number I of control time index, where [[bar.F].sub.Oi=1,"k"] = 0 describes the beginning of phase transformation mode and [[bar.F].sub.Oi=1,"k"] = 1 refers its end. Condensing and unsteady evaporation modes are modeled for case I = 41. Individual time grid for each mode are defined according to scheme [[bar.F].sub.Oi,"k"] [right arrow] [Fo.sub.i,"k"], [right arrow] [[tau].sub.i]. Necessary values of numbers [Fo.sub.fv,"k"] were guessed. They can be selected according to [14, 15] recommendations.

[FIGURE 1 OMITTED]

Assumption of a droplet dimension stability influence in iterative cycle is primary evaluated in condensing phase transformation mode at "k" heat transfer case. At discussed initial conditions case a dew point temperature is [T.sub.rt] = 342.275 K. Reference condensing phase transition regime duration [Fo.sub.ko,"k"] [approximately equal to] 0.72 is selected according to the parameter [T.sub.0]/[T.sub.rt] = 0.812. Two cases were modeled, for droplet phase transformation parameters Pf in numerical schemes when conditions [R.sub.i,it] [equivalent to] [R.sub.i-1] or [R.sub.i,it] [equivalent to] [R.sub.i,it-1] are selected for temperature [T.sub.R,i,"k"] definition. For other parameters condition [R.sub.i,it] [equivalent to] [R.sub.i-1] was kept the same. Results of numerical research are graphically summarized in Figs. 1-7. In both cases, the fastest descent method provides a high Eq. (4) condition requirement (Fig. 1).

[FIGURE 2 OMITTED]

[FIGURE 3 OMITTED]

In iterative cycles for ensuring Eq. (4) conditions requirement by accuracy of [[delta].sub.i,IT] < [+ or -] 0.04% (Fig. 1, 2), a droplet surface temperature change at hundredth degree series must be taken into account (Fig. 3). In primary stage of condensing phase mode [F.sub.0] [equivalent to] 0 - 4 a calculated intensity of droplet warming heat flux [q.sup.-.sub.k] = [q.sup.+.sub.k] + [q.sup.+.sub.f] is higher in case [R.sub.i,it] = [R.sub.i,it-1], while in final stage this intensity is higher at [R.sub.i,it] [equivalent to] [R.sub.i-1] case (Fig. 4, a). This initial phase stage for energy assessment is more significant in condensing mode, therefore in case [R.sub.i,it] [equivalent to] [R.sub.i,it-1] a droplet surface warmth first up to dew point temperature (Fig. 4, b).

Durations of condensing phase transformation mode were qualified by additional numerical research: [Fo.sub.ko,"k"] = 0.722, when [R.sub.i,it] [equivalent to] [R.sub.i,it-1] and, [Fo.sub.ko,"k"] = 0.771 when [R.sub.i,it] [equivalent to] [R.sub.i-1].

[FIGURE 4 OMITTED]

In both cases a droplet central and surface layers warming dynamics is different: a droplet surface layers warming rate rapidly suffocates, while central layers warming rate grows at the beginning of condensing mode and at the moment [Fo.sub.] [approximately equal to] 0.2 it reaches maximum (Fig. 5). Then gradually suffocates, remaining greater than surface warming rate. Therefore a bright non-isothermal observes at condensation mode in a droplet (Fig. 4, b). Droplet heating (Fig. 5) and phase transformation (Fig. 6) speed dependence from droplet dimension selection method in iterative cycle causes a bright calculated droplet diameter changes for assumptions [R.sub.i,it] [equivalent to] [R.sub.i,it-1] and [R.sub.i,it] [equivalent to] [R.sub.i-l] cases. (Fig. 7).

Unsteady phase transformation mode 0--1--2 is modeled in its numeric scheme for condensing mode maintaining the previous time grid, and for nonstationary evaporation time grid is formed in whole number I = 30 basis. Value [[tau].sub.i=41] represents the end of condensation mode, which coincides with the beginning of unsteady evaporation regime. By equilibrium evaporation conditions approaching asymptotically.

[FIGURE 5 OMITTED]

[FIGURE 6 OMITTED]

Therefore a beginning of equilibrium evaporation definition requires for agreement of remaining condition. When this condition is being satisfied evaporation is considered to be equilibrium. This condition can be defined according calculated droplet surface temperature deviation from theoretical temperature, that ensures equilibrium evaporation [T.sub.R,e] - [T.sub.i] < [DELTA][T.sup.leist.sub.R,e] or by permissible nonisothermally [T.sub.R,i] - [T.sub.C,i] < [DELTA][T.sup.leist.sub.R,C] [15]. At discussed primary conditions droplets that are heated by conductivity in unsteady phase transition mode warmth to [T.sub.m,e] [equivalent to] [T.sub.R,e] = 348.565 K. Condition [DELTA][T.sup.leist.sub.R,e] = 0.01 K. is provided for equilibrium evaporation. For guessed [Fo.sub.nf,"k"] = 3 the duration of unsteady phase transformation [Fo.sub.nf,"k"] is defeined graphically (Fig. 8): 2.27 and 2.07 for conditions [R.sub.i,it] [equivalent to] [R.sub.i-1] and [R.sub.i,it] [equivalent to] [R.sub.i,it-1], respectively. Then unsteady evaporation duration [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] is defined: 1.5 and 1.35, respectively.

[FIGURE 7 OMITTED]

[FIGURE 8 OMITTED]

Droplet phase transformation cycle 0 - 1 - 2 - 3 is modeled in numeric scheme for condensing and unsteady evaporation mode maintaining the previous time grid. For equilibrium evaporation mode a time grid is formed in base of whole number I = 50. Here [[tau].sub.i-71] represents the beginning of equilibrium evaporation mode, which matches with the end of unsteady phase transformation mode. The duration of equilibrium evaporation [Fo.sub.eg,"k"] is 23.53 for [R.sub.i,it] [equivalent to] [R.sub.i-1] and 34.63 for [R.sub.i,it] [equivalent to] [R.sub.i,it-1], respectively. These durations were selected in numerical research in order that in i = 121 iterative cycle must ensure a droplet evaporation.

A phase transformation cycle is defined according to "k" heat transfer modeling results: 0 - 0.771 - 2.27 - 25.8 when [R.sub.i,it] [equivalent to] [R.sub.i-1] and 0--0.722--2.07-36.7, when [R.sub.i,it] [equivalent to] [R.sub.i,it-1]. Droplet dimension selection method influence for comparative evaluation of the main phase transformation parameters for [R.sub.i,tt] [equivalent to] [R.sub.i-1] and [R.sub.i,it] [equivalent to] [R.sub.i,1] cases are given in (Table). Nonvalidity of assumption [R.sub.i,it] [equivalent to] [R.sub.i-1] application for droplet phase transformation cycle in "k" heat transfer case modeling is highlight.

It was state that in combined droplet heating case by conductivity-radiation, absolutely black body radiation source radiates air temperature. Due to peculiarities of radiant flux absorption in droplets, at combined heat transfer case a droplet phase transformation cycle depends from droplets dispersity. A quantitative indicators of the cycle that are listed below valid only for modeled 2[R.sub.0] [equivalent to] 150 x [10.sup.-6] m droplet diameter. A comparative analysis of modeled "k" and "k+r" cases gives a qualitative assessment. For combined heat transfer case a time grid of individual phase transformation in numerical schemes formulation methodology was kept the same as for conduction heating.

Due to additional heat flux from radiation source a condensing mode becomes shorter in "k+r" case:[Fo.sub.ko,"k+r"] = 0.764 when [R.sub.i,it] [equivalent to] [R.sub.i-1] and [Fo.sub.ko,"k+r"] = 0.713 when [R.sub.i,it] [equivalent to] [R.sub.i,it-1]. Because of intense surface layers warming, bright non-isothermally is proper for droplet in condensation mode. Peculiarities of droplet surface and central layers warming rate (Fig. 5) causes a formation of non-isothermal peak at the beginning of unsteady evaporation (Fig. 9). It observes in all heat transfer modes and is the same at qualitative assessment. At compound heat transformation case a radiation heat flux that is absorbed in a droplet warmth more droplet central layers (Fig. 10). This forms assumptions for the second non-isothermal peak (Fig. 10). In case of "k+r" droplet begins to evaporate equilibrium at the second non-isothermal peak formation time, when [DELTA][T.sub."k+r"]([tau] [equivalent to] [[tau].sub.e]) [equivalent to] [DELTA][T.sub.max,2]. According to that Fourier number [Fo.sub.nf,"k+r"] is defined: 3.16, when [R.sub.i,it] [equivalent to] [R.sub.i-1] and 2.91, when [R.sub.i,tt] [equivalent to] [R.sub.i,it-1]. Then time of unsteady evaporation mode [Fo.sub.ng,"k+r"] = [Fo.sub.nf,"k+r"] - [Fo.sub.ko,"k+r"] is 2.396 and 2.197, respectively.

In the end of unsteady evaporation a temperature gradient inside a droplet provides full heat that is absorbed by radiation output to droplet surface by conductivity; therefore all external compound heat flux, starts participate in evaporation process whose intensivity is reflected by its density [q.sup.+.sub.[SIGMA]] = [q.sup.+.sub.k] + [q.sub.+.sub.r] (Fig. 11).

[FIGURE 9 OMITTED]

[FIGURE 10 OMITTED]

Phase tansformation cycle duration [Fo.sub.f,"k+s"] is defined by numerical research: 23.53, when [R.sub.i,it] [equivalent to] [R.sub.i-1] and 34.63, when [R.sub.i,it] [equivalent to] [R.sub.i,it-1]. Then equilibrium evaporation mode [Fo.sub.eg,"k+r"] = [Fo.sub.f,"k+r"] - [Fo.sub.f,"k+r"] lasts 22.86 and 34, respectively. At equilibrium evaporation mode droplets decreases rapidly (Fig. 12), therefore the influence of radiation source fades away and non-isotermality suffocates at a droplet.

According to generalized "k+r" heat transfer modeling results, as well in case "k", a weakness of assumption [R.sub.i,it] [equivalent to] [R.sub.i-1] application for phase cycle modeling, was highlight.

The results were confirmed by different calculating methods of droplet mass dynamics in Fig. 12, curves 4 and 5: curves 4 are formed by numerical scheme that is based on by a model of droplet mass change Eq. (13) and curves 5 are based on by linking droplet mass and volume of non-isothermal droplet mass medium density [M.sub.l] [equivalent to] [V.sub.l] x [[rho].sub.l,m].

[FIGURE 11 OMITTED]

[FIGURE 12 OMITTED]

Water density is selected by expression (11), which describes temperature [[rho].sub.l,m] = [[rho].sub.l]([T.sub.m]). Functions [M.sub.l]([Fo.sub.]) graphs matches on assumptions [R.sub.i,it] = [R.sub.i,it-1] case and for [R.sub.i,it] = [R.sub.i-1] case are markedly different (Fig. 12).

4. Conclusions

A droplet diameter selection method in iterative cycle importance at humid air sprayed water phase transformation cycle impact for modeling results of droplet heating by conduction and combined heating by conduction-radiation, was highlighted.

Assumption [R.sub.i,it] [equivalent to] [R.sub.i,it-1] that was applied in numeric schemes for phase transformation parameters [P.sub.f] can be extended in thermal [P.sub.T], energy [P.sub.q] and dynamic [P.sub.d] droplet parameters at numerical calculation algorithms.

Receved August 25, 2014

Accepted December 15, 2014

References

[1.] Sirignano, W.A. 1999. Fluid dynamics and transport of droplets and sprays. Cambridge: Cambridge University Press. 311 p. http://dx.doi.org/10.1017/CBO9780511529566.

[2.] Miliauskas, G 2014. Universal cycle of the phase transformation regimes for pure liquid droplets 1. Cycle formation method and numerical modeling framework for droplet transfer processes, Power Engineering 60-2: 77-95. http://dx.doi.org/10.6001/energetika.v60i2.2933.

[3.] Miliauskas, G; Sabanas, V. 2006. Interaction of transfer processes during unsteady evaporation of water droplets, Int. J. Heat Mass Transfer 49: 1790-1803. http://dx.doi.org/10.1016/j.ijheatmasstransfer.2005.11.0 15.

[4.] Miliauskas, G; Sinkunas, S. 2008. Influence of droplets heating manner to the sprayed liquid thermal state, CHT'08: proceedings of the International Symposium of the Advances in Computational Heat Transfer (Moroco): 1-17. http://dx.doi.org/10.1016/j.ijheatmasstransfer.2005.11.0 15.

[5.] Zuev, V.E., Kopitin, U.D., Kuzikovskii, A.V. 1982. Nonlinear optical effects in aerosols, Novosibirsk, Nauka. 237 p.

[6.] Tuntomo, A.; Tien, C.L.; Park, S.H. 1991. Internal distribution of radiant absorption in a spherical particle, J. Heat Transfer 113: 407-412. http://dx.doi.org/10.1115/L2910576.

[7.] Miliauskas, G 2001. Regularities of unsteady radiative-conductive heat transfer in evaporating semitransparent liquid droplets, Int. J. Heat Mass Transfer 44: 785-798. http://dx.doi.org/10.1016/S0017-9310(00)00127-7.

[8.] Tseng, C.C.; Viskanta, R. 2006. Enhancement of water droplet evaporation by radiation absorption, Fire Safety J. 41: 236-247. http://dx.doi.org/10.1016/j.firesaf.2006.01.001.

[9.] Miliauskas, G; Sabanas, V.; Bankauskas, R.; Miliauskas, G; Sankauskaite, V. 2008. The

peculiarities of sprayed liquid's thermal state change, as droplets are heated by conduction, Int. J. of Heat and Mass Transfer 51: 4145-4160. http://dx.doi.org/10.1016/j.ijheatmasstransfer.2008.01.0 16.

[10.] Vukalovitch, M.P. Thermodynamics properties of water and steam, VEB Verlag Technik, Berlin, 1958. 168 p.

[11.] Abramzon, B., Sirignano, W.A. 1989. Droplet vaporization model for spray combustion calculations, Int. J. of Heat and Mass Transfer 32: 1605-1618. http://dx.doi.org/10.1016/0017-9310(89)90043-4.

[12.] Harpole, GM. 1980. Radiative absorption by evaporating droplets, Int. J. Heat Mass Transfer 22: 1726. http://dx.doi.org/10.1016/0017-9310(80)90134-9

[13.] Miliauskas, G., Norvaisien?, K. 2013. A systematic evaluation of the unsteady transfer process interaction in evaporating droplets, Power Engineering 59-1: 2641. http://dx.doi.org/10.6001/energetika.v59i1.2617.

[14.] Miliauskas, G., Sink?nas, S., Miliauskas, G. 2010. Eva-poration and condensing augmentation of water droplets in flue gas, Int. J. of Heat and Mass Transfer 53: 1220-1230. http://dx.doi.org/10.1016/j.ijheatmasstransfer.2009.10.0 05.

[15.] Miliauskas, G., Norvaisien?, K. 2013. The lifetime cycle of sprayed water and the universal duration diagram of droplets unsteady evaporation in dry air, Power Engineering 59-2: 57-68. http://dx.doi.org/10.6001/energetika.v59i2.2661.

G. Miliauskas, M. Maziukiene

Kaunas University of Technology, K. Donelaicio 20, Kaunas, LT-44239, Lithuania, E-mail: gimil@ktu.lt

http://dx.doi.Org/ 10.5755/j01.mech.20.6.8748

Table

Comparative evaluation of the main phase
transformation parameters

                                            [R.sub.i,it]
                                            [euqivalent to]
                                            [R.sub.i-1]

[R.sub.ko] x [10.sup.-6], m                 83.31
[M.sub.ko,"k"] x [10.sup.10], kg            19.078
[Fo.sub.ng,"k"]                             1.5
[R.sub.nf,"k"] [10.sup.-6], m               82.3
[m.sup.+.sub.v,nf,"k"], kg/([m.sup.2]s):    0.02388
[g.sup.+.sub.v,nf,"k"] x [10.sup.7], kg/s   0.0203
[M.sub.nf,"k"] x [10.sup.10], kg            18.008
[FO.sub.eg,"k"]                             23.53

                                            [R.sub.i,it]
                                            [euqivalent to]
                                            [R.sub.i,it-1]

[R.sub.ko] x [10.sup.-6], m                 77.35
[M.sub.ko,"k"] x [10.sup.10], kg            18.979
[Fo.sub.ng,"k"]                             1.35
[R.sub.nf,"k"] [10.sup.-6], m               76.19
[m.sup.+.sub.v,nf,"k"], kg/([m.sup.2]s):    0.0258
[g.sup.+.sub.v,nf,"k"] x [10.sup.7], kg/s   0.0188
[M.sub.nf,"k"] x [10.sup.10], kg            18.083
[FO.sub.eg,"k"]                             34.63