D 2 Law For Liquid Droplet Vaporization

1. D2 Law For Liquid Droplet Vaporization • References: • Combustion and Mass Transfer, by D.B. Spalding (1979, Pergamon Press). • “Recent advances in droplet vaporization and combustion”, C.K. Law, Progress in Energy and Combustion Science, Vol. 8, pp. 171-201, 1982. • Fluid Dynamics of Droplets and Sprays, by W.A. Sirignano (1999, Cambridge University Press).

2. Heat Fuel Vapour Droplet With Internal Circulation Gas-Phase Streamlines Spherical Symmetry Model Buoyancy and droplet size destroy spherical symmetry evaporation Droplet With Internal Circulation Gas-Phase Streamlines

3. Without/micro gravity With gravity

4. Mass Transfer I DEFINITIONS IN USE: • density – mass of mixture per unit volume ρ [kg/m3] • species - chemically distinct substances, H2O, H2, H, O2, etc. • partial density of A – mass of chemical compound (species) A per unit • volume ρA [kg/m3] • mass fraction of A – ρA/ρ = mA • note: • ρA + ρB + ρC + … = ρ • mA + mB + mC + … = 1

5. DEFINITIONS IN USE: • total mass velocity of mixture in the specified direction (mass flux density) • – mass of mixture crossing unit area normal to this direction in unit time • GTOT [kg/m2s], GTOT = ru (density x velocity) • total mass velocity of A in the specified direction = GTOT,A[kg/m2s] • note: GTOT,A + GTOT,B + GTOT,C …= GTOT • convective mass velocity of A in the specified direction • mAGTOT = (ρA/ρ) GTOT= GCONV,A • note: GCONV,A + GCONV,B + GCONV,C …= GTOT • but generally, GCONV,A ≠ GTOT,A • diffusive mass velocity of A in the specified direction • GTOT,A – GCONV,A = GDIFF,A • note: GDIFF,A + GDIFF,B + GDIFF,C + … = 0

6. DEFINITIONS IN USE: • velocity of mixture in the specified direction = GTOT/r [m/s] • concentration – a word used for partial density or for mass fraction • (or for mole fraction, partial pressure, etc.) • composition of mixture – set of mass fractions

8. mass flux in - mass flux out mass accumulated =

9. The d2 Law - assumptions • Spherical symmetry: forced and natural convection are neglected. This reduces the analysis to one-dimension. • No spray effect: the droplet is an isolated one immersed in an infinite environment. • Diffusion being rate controlling. The liquid does not move relative to the droplet center. Rather, the surface regresses into the liquid as vaporization occurs. Therefore heat and mass transfer in the liquid occur only because of diffusion with a moving boundary (droplet surface) but without convection. • Isobaric processes – constant pressure. • Constant gas-phase transport properties. This causes the major uncertainty in estimation the evaporation rate (can vary by a factor of two to three by using different, but reasonable, averaged property value – specific heats, thermal conductivity, diffusion coefficient, vapour density, etc). • Gas-phase quasi-steadiness. Because of the significant density disparity between liquid and gas. Liquid properties at the droplet surface (regression rate, temperature, species concentration) changes at rates much slower than those of gas phase transport processes. This assumption breaks down far away from the droplet surface where the characteristic diffusion time is of the same order as the surface regression time.

10. Gas-phase QUASI-steadiness – characteristic times analysis. In standard environment the gas-phase heat and mass diffusivities, ag and dg are of the same order of 100 cm2s-1, whereas the droplet surface regression rate, K = -d(D02)/dt is of the order of 10-3cm2s-1 for conventional hydrocarbon droplet vaporizing in standard atmosphere. Thus, there ratio is of the same order as the ratio of the liquid-to-gas densities, . It means that gas mass and heat diffusion occurs much faster than droplet surface regression time. If we further assume that properties of the environment also change very slowly, then during the characteristic gas-phase diffusion time the boundary locations and conditions can be considered to be constant. Thus the gas-phase processes can be treated as steady (time independent because the surface almost “freezes”), with the boundary variations occurring at longer time scales. When (at which value of D∞) this assumption breaks down, i.e. when the diffusion time is equal to the surface regression time? When surface regression characteristic time becomes equal to gas mass/heat diffusion time, i.e. when D∞2/ dg ≈ D02/K? Remembering that must still be valid (it doesn’t depend on the distance from the droplet), we can conclude that the steady assumption breaks down when . For standard atmospheric conditions it breaks down at It means that our model will be valid for the distances less than this one.

11. The d2 Law – assumptions (vii) Single fuel species. Thus it is unnecessary to analyze liquid-phase mass transport (no diffusion term). (viii) Constant and uniform droplet temperature. This implies that there is no droplet heating. Where all the heat goes? Combined with (vii), we see that liquid phase heat and mass transport processes are completely neglected. Therefore the d2 Law is essentially a gas-phase model. (ix) Saturation vapour pressure at droplet surface. This is based on the assumption that the phase-change process between liquid and vapour occurs at a rate much faster than those for gas-phase transport. Thus, evaporation at the surface is at thermodynamic equilibrium, producing fuel vapour which is at its saturation pressure corresponding to the droplet surface temperature. (x) No Soret (mass flow because of the temperature gradient), Dufour (heat flow because of the concentration gradient) and radiation effects (how this effects the validity of the model?).

12. Rate of accumulation of mass of component j Rate of generation of mass of component j from reaction Rate of depletion of mass of component j from reaction Mass flow rate of component jinto the system Mass flow rate of component jout of system

14. Droplet evaporation I (no energy concerns) • The phenomenon considered: • A small sphere of liquid in an infinite gaseous atmosphere vaporizes and • finally disappears. • What is to be predicted? • Time of vaporization as a function of the properties of liquid, vapor and environment. • Assumptions: • spherical symmetry (non-radial motion is neglected) • (quasi-) steady state in gas • ΓVAP independent of radius • large distance between droplets • no chemical reaction

15. G = GTOT,VAP Go ro r Vapor concentration distribution mVAP in the gas.

17. 1 mVAP,0 Why the curve doesn’t start at mVAP=1? mVAP mVAP,∞ Droplet surface r = r0 r

20. Example

21. Limitations • mVAP,0 has a strong influence, but is not usually known, it depends on temperature. • relative motion of droplet and air augments the evaporation rate (inner circulation of the liquid) by causing departures from spherical symmetry. • the vapour field of neighboring droplets interact • mVAP,0 and mVAP,∞may both vary with time. • ΓVAP usually depends on temperature and composition.

22. The Energy Flux DEFINITIONS IN USE:

23. S E +dE E Dx

25. 0

26. 0, for the case of Stefan flow

27. Droplet evaporation II Go G = GTOT,VAP ro E r Qo heat flow to gas phase close to liquid surface

29. y = -x 1 0 So, a positive G0 reduces the rate of heat transfer at the liquid surface. It means that if the heat is transferred to some let us say solid surface, that we want to prevent from heating up, we should eject the liquid to the thermal boundary layer (possibly through little holes). This liquid jets will accommodate a great part of the heat on vaporization of the liquid. Thus, we’ll prevent the surface from heating – transpiration cooling. The smaller the holes the smaller a part of heat towards the liquid interior and, subsequently towards the solid surface.

31. Clausius-Clayperon equation for

33. Linkage of equations

34. Equilibrium vaporization – droplet is at such a temperature that the heat transfer to its surface from the gas is exactly equals the evaporation rate times the latent heat of vaporization: This implies: What if –Q0 ≠ G0L

36. Unity Le number

37. ! Temperature and concentration profiles look exactly the same.

39. Cases of interest: • When T∞ is much greater than the boiling-point temperature TBOILING, mVAP,0 is close to 1 and T0 is close to TBOILING. Then the vaporization rate is best calculated from: • When T∞ is low, and mVAP,∞ is close to zero, T0 is close to T∞. This implies T0≈T∞. Thus, mVAP,0 is approximately equal to the value given by setting T0=T∞ in and the vaporization rate can be calculated by: As in example with water droplet evaporating at 100C

40. The choice depends on whether T0 or mVAP,0 is easier to estimate Evaporation rate [m2/s]

41. Qualitative results for D2-Law

42. Droplet heat up effect on temperature and lifetime We can divide the droplet evaporation process into two stages. At first, while the droplet is cold (evaporation is slow), all the heat from the hot environment will be used for droplet interior heat up. As the droplet temperature will approach its steady state value, droplet heat up will slow down, while evaporation will accelerate. r0=r0(t)

43. Diffusion limit Distillation limit D2 Law Droplet heat up effect on temperature and lifetime Fastest limit Slowest limit (r0/r0,INITIAL)2 T 1 Diffusion limit Distillation limit 380 D2 Law Surface Temperature Center Temperature 0 300 (aLIQ/r0,INITIAL2)t (aLIQ/r0,INITIAL2)t 0.1 0.2 0.1