Simultaneous Measurement of Condensation and Thermal Accommodation Coefficients for Cloud Droplet Growth in Due Consideration of a New Moving Surface-Boundary Effect

Norihiko Fukuta Department of Meteorology, University of Utah, Salt Lake City, Utah

Search for other papers by Norihiko Fukuta in
Current site
Google Scholar
PubMed
Close
and
Marcus N. Myers Department of Meteorology, University of Utah, Salt Lake City, Utah

Search for other papers by Marcus N. Myers in
Current site
Google Scholar
PubMed
Close
Full access

Abstract

A droplet growth theory that describes a new effect of vapor and temperature field shift due to the growth-based movement of droplet surface boundary (moving boundary effect) was derived and found to enhance the growth rate as a function of supersaturation (S − 1) and droplet radius (a) in second order. The theory plays the role of a missing link and resolves the gap that exists among the measured data; the high (S − 1) used was found to overestimate thermal accommodation and condensation coefficients, α and β. The theory provided the basis for bettering the measurement and correction for data analyzed without the moving boundary effect, and suggested the broadening effect in the droplet size spectrum by accelerating the growth on the larger end together with the growth slowdown effect of α, β on the smaller.

Measurements employing a horizontal-flow thermal diffusion chamber and the Mie scattering method using a He–Ne laser for droplet size estimation at temperature T = 277 K, (S − 1) = 0.32%, a ≈ 2 μm under two pressures of modulation, p = 100 and 860 hPa, with correction of the moving boundary effect for nuclei of NaCl and (NH4)2SO4, resulted in the simultaneous determination of α = 0.81 ± 0.07 and β = 0.043 ± 0.016 on average. For paper smoke, α = 0.68 and β = 0.022. Comparison with other data of the low moving boundary effect suggests a slight β increase with temperature lowering. Based on the result of the present measurement and other theoretical and experimental works, the mechanism of β being much smaller than α was suggested mostly as a kinetic process controlled by the removal rate of the released latent heat at the moment of molecular impact on the surface and vice versa for evaporation.

Corresponding author address: Norihiko Fukuta, Department of Meteorology, University of Utah, 135 S 1460 E, Rm. 819, Salt Lake City, UT 84112-0110. Email: nfukuta@met.utah.edu

Abstract

A droplet growth theory that describes a new effect of vapor and temperature field shift due to the growth-based movement of droplet surface boundary (moving boundary effect) was derived and found to enhance the growth rate as a function of supersaturation (S − 1) and droplet radius (a) in second order. The theory plays the role of a missing link and resolves the gap that exists among the measured data; the high (S − 1) used was found to overestimate thermal accommodation and condensation coefficients, α and β. The theory provided the basis for bettering the measurement and correction for data analyzed without the moving boundary effect, and suggested the broadening effect in the droplet size spectrum by accelerating the growth on the larger end together with the growth slowdown effect of α, β on the smaller.

Measurements employing a horizontal-flow thermal diffusion chamber and the Mie scattering method using a He–Ne laser for droplet size estimation at temperature T = 277 K, (S − 1) = 0.32%, a ≈ 2 μm under two pressures of modulation, p = 100 and 860 hPa, with correction of the moving boundary effect for nuclei of NaCl and (NH4)2SO4, resulted in the simultaneous determination of α = 0.81 ± 0.07 and β = 0.043 ± 0.016 on average. For paper smoke, α = 0.68 and β = 0.022. Comparison with other data of the low moving boundary effect suggests a slight β increase with temperature lowering. Based on the result of the present measurement and other theoretical and experimental works, the mechanism of β being much smaller than α was suggested mostly as a kinetic process controlled by the removal rate of the released latent heat at the moment of molecular impact on the surface and vice versa for evaporation.

Corresponding author address: Norihiko Fukuta, Department of Meteorology, University of Utah, 135 S 1460 E, Rm. 819, Salt Lake City, UT 84112-0110. Email: nfukuta@met.utah.edu

1. Introduction

It is now well accepted that the development of cloud properties and subsequent precipitation processes becomes a primary influence during the early stage of formation. When moist air lifts past the cloud base to continually expand, cool, and generate supersaturation, the cloud condensation nuclei (CCN) activate from the low supersaturation end of the activity spectrum. The nucleated and growing droplet begins to remove water vapor from and to return the generated heat of condensation to the air environment and reduce the rate of supersaturation increase. During increasing supersaturation, more and more CCN nucleate, and the formed droplets grow. The reduction rate of supersaturation eventually catches up to the generation rate to a balance, and at this supersaturation maximum, the CCN activation ceases. There, the initial number density of droplets is decided. The process is thus a nucleation–growth interaction in the thermodynamically changing air environment. A similar interaction takes place in the process of aerosol generation.

In the nucleation–growth interaction above the cloud base, the contribution of the CCN activity spectrum is well recognized (Twomey 1977). In the droplet growth that follows the nucleation, incoming vapor flow toward and outgoing flow of the generated condensation heat from the surface establish a balance. The two balanced counterflows connect the respective values of saturation at the surface. Nevertheless, the processes became known to depend on the efficiency of exchanging both properties at the droplet surface, for heat the thermal accommodation coefficient, α, and for the mass of vapor the condensation coefficient, β (sometimes called the sticking or mass accommodation coefficient). In the early days, such inefficiencies were not considered, and the droplet growth process is now referred to as the Maxwellian (Maxwell 1890), which assumes continuities of the two properties at the droplet surface. The inefficiency appears when the free molecular gas kinetic flows just above the surface come into play as a part of continuum diffusional flows outside. The growth process that involves the two serial flows is called the “diffusion-kinetic” (Smirnov 1971). Early diffusion-kinetic theories involved one of the two flows only; that is, heat or vapor mass flow. When the first analytic solution including diffusion-kinetic flows of both vapor and heat was obtained (Fukuta and Walter 1970), the significant contributions in cloud microphysical processes were realized, such as the nucleation in gaseous-phase homogeneous condensation in particular (Fukuta 2004). The droplet size spectrum spreads by slowing down growth more on the smaller side, and the number density increases to influence cloud lifetime, albedo, and precipitation (Zou and Fukuta 1999).

In the application of the diffusion-kinetic theory for measurement and analysis, the thermal accommodation coefficient used to be assumed as unity. Measurements were difficult and values of the condensation coefficient reported in the last half a century scatter over two orders of magnitude (Fukuta and Walter 1970; Pruppacher and Klett 1978, 1997; Mozurkewich 1986). The condensation coefficient seems to group about two values around unity and 0.04. The problem thus boils down to accurately measure both coefficients in a properly designed experiment.

We have found during the theoretical development that contributions of atmospheric pressure to the condensation and thermal accommodation coefficient terms are different. Two sets of conditions can be used to solve for the two independent accommodation coefficients simultaneously. Measurements were designed and carried out under conditions close to atmospheric cloud formation, selecting the droplet size range and the method of measurement sensitive to distinguish the diffusion-kinetic effect. In the analysis of measured data and comparison with the work of others, it was necessary to clarify the previously neglected effect arising from the shift of steady-state temperature and vapor fields, caused by the movement of growing droplet surface boundary. The attempt for the theoretical derivation was fruitful. It appeared to be the missing link between the two groups of reported coefficients dissolving the gap. The correction was obtained for original data of the present measurement. It even showed a new effect for broadening of droplet size spectrum. It is the purpose of this work to report the result.

2. Theory

For experimental determination of the condensation and thermal accommodation coefficients, the use of a growing or evaporating droplet system is advantageous since all the fluxes of heat and vapor are in the steady-state (or quasi steady state) balance in between the outer boundary of air environment and the inner boundary at the droplet surface.

a. Behavior of the diffusion-kinetic droplet growth equation

In the steady growth of droplets in air to determine α and β, the analytic solution exposes the relationships among all the variables and constants. Direct application of the primitive transport equations without obtaining the solution and adopting them in the curve fitting method tends to miss some hidden relationships. While linearized suitably for natural cloud conditions, the solution of Fukuta and Walter (1970, hereafter FW70) with a minor improvement for treatment of both accommodation coefficients (Fukuta and Xu 1996, hereafter FX96) gives the basis for analysis of droplet growth under supersaturated vapor of the environment.

For a droplet with a radius, a, and sufficiently dilute solution with a negligible Kelvin effect (for notations, see appendix A), FW70 and FX96 give
i1520-0469-64-3-955-e1
Equation (1) is now found to be valid in the entire range of a and integrable under the assumed condition. The integral form, with a = 0 at t = 0, takes the form
i1520-0469-64-3-955-e2
or
i1520-0469-64-3-955-e3
where C1 and C2 are constants for the air environment in question; α and β are thermal accommodation and condensation coefficients defined by Knudsen (1934) and included in continuum heat and vapor fluxes at the droplet surface (Chapman and Cowling 1970; FX96); α′ and β′ are representative values of α and β that appear in the free molecular flux density of Hertz–Knudsen equation as well as FW70’s.

In Eqs. (1)(3), it may be seen that only when a → ∞ (and fictitious condition of α′ = β′ = ∞ or α = β = 2), fα = fβ = 1, and Eq. (1) approaches to the limit; that is, the Maxwellian with C2 = 0 in Eq. (3). In Eq. (3), C2 ≠ 0 describes the diffusion-kinetic. This is to say that both equations take the same parabolic Maxwellian-functional style and they would be indistinguishable if an arbitrary parallel shift were applied in a curve fit analysis. The parallel shift for the diffusion-kinetic is restricted with t = 0, a = 0, which is accompanied by another unreal condition, t = 0, a < 0 (see Fig. 1).

The slowdown effect of the diffusion-kinetic droplet growth rate from the Maxwellian is a function of the droplet size, a, and it is desirable to plan the measurement in the a range where the ratio between Maxwellian and diffusion-kinetic is sufficiently large. Figure 2 shows r and dm/dt ratios between Maxwellian and diffusion-kinetic growths for a case of α′ = 1, β′ = 0.03, at 0°C, 1 atm, and (S − 1) = 0.003. Under these assumed conditions, the ratios are expected to be larger than 2 when aMW < 4 μm.

The fall velocity of droplet thus growing may be expressed by the Stokes–Cunningham equation
i1520-0469-64-3-955-e4
where C3 is a constant near unity. The equation provides a basis to assess the range of fall distance of the growing droplet to stay within the width of the laser beam of the measurement.

It should be pointed out that the diffusion-kinetic equation used above does not include the Kelvin (or the curvature) and solute effects of the nucleus particle, which produce a single vapor pressure maximum, reduce the vapor pressure difference or the effective (S − 1) between the environment and lower the droplet growth rate. The maximum leads to saturation as the droplet grows. The maximum becomes smaller as the dry mass of the nucleus increases. Through selection of the nucleus mass to keep the maximum sufficiently low and use of (S − 1) in the growth environment sufficiently high, one can find a combination to make the effects essentially negligible on the condition used in the above diffusion-kinetic growth equation. The details of the combination actually applied are given in section 2d.

b. Droplet growth environment

In clouds, a droplet grows in supersaturated air caused by the adiabatic expansion of lifting. In the laboratory, the adiabatic expansion method is apt to cause a damping oscillation in pressure. In addition, supersaturation deviations from the aimed value due to wall disturbances may occur. The thermal diffusion chamber method offers a steady profile of supersaturation, both in static and flow types, if and when both warmer top and cooler bottom plates are securely saturated with respect to water. The temperature field in the chamber establishes linearly with respect to the height. The saturation ratio, S, and supersaturation, (S − 1), are a function of the vapor pressure ratio between the real and saturation values or S = e/es. The es may be easily determined from the temperature since es = f (T). The e profile varies depending on the property assumed in the linear steady-state profile, such as e (Katz and Mirabel 1975) or ρ (Chodes et al. 1974). It is, however, most reasonable to take the linear vapor flux density profile with height or a constant rate of molecular transport in the existing temperature field (Fukuta and Gramada 2003). Between the parallel top and bottom plates of the chamber under the steady state,
i1520-0469-64-3-955-e5
and w2T and ne/T in the gas kinetics lead to
i1520-0469-64-3-955-e6
or
i1520-0469-64-3-955-e7

The e value at the midpoint of the chamber estimated based on relationship (7) between the top and bottom plates gives (S − 1) = 0.0032 at 4°C with a temperature difference of 2.5°C. Under the conditions, assumption of linear profiles with height for e and ρ leads to an overestimation of (S − 1) by 6% and an underestimation by 8%, respectively.

c. Method for simultaneous determination of thermal accommodation and condensation coefficients

The thermal accommodation and condensation coefficients in Eqs. (1) and (2) appear combined with terms of different pressure dependence. Measurements of droplet growth under two different pressures, therefore, yield two sets of growth equations to be simultaneously solved for the coefficients. By exposing the pressure dependence, the solutions may be obtained as
i1520-0469-64-3-955-eq1
and
i1520-0469-64-3-955-e8
where
i1520-0469-64-3-955-e9
α′ and β′ are representative values used in FW70 and others and are related to the true values, α and β, as shown in Eq. (8). When they are small, both values coincide. The method thus requires experimental determination of both da/dt and a or alternatively t values under two pressures, p1 and p2, at constant temperature.

d. The initial method for droplet growth rate estimation

There are three methods commonly used for droplet size determination: the fall velocity measurement in the Stokes regime with analysis by Eq. (4), the electrodynamic levitation method, and the Mie peak identification method of monochromatic light scattering. The fall velocity measurement offers a simultaneous determination of r in Eq. (4), but when droplets are small, the measurement will receive the updraft/downdraft effect from the heat of condensation/evaporation as well as the diffusio- and thermophoretic effect if used in a horizontal thermal diffusion chamber. The electrodynamic levitation method compensates the fall of a droplet to a standstill, but for the large droplet size that is usable in levitation, the diffusion-kinetic effect diminishes. We selected the Mie-peak detection method of light scattering utilizing a He–Ne laser (λ = 0.6328 μm) in a flow-type horizontal thermal diffusion chamber.

Figure 3 shows the computed intensity of scattered light for a He–Ne laser at 90° plotted as a function of droplet radius. Upon an inspection of the figure, it is clear that the interval between two adjacent scattering peaks, Δa(= 0.162 μm), is constant so that by measuring the time elapsed between the two adjacent flashes of light or peaks, Δt, the radial growth rate,
i1520-0469-64-3-955-e10
may thus be easily measured. The second value needed in Eqs. (8) and (9), a, may be determined by measuring t, the time elapsed until the moment of (da/dt) measurement, and applying it to Eq. (1).

The conditions used in the present measurement are typically (S − 1) = 3.2 × 10−3 and T = 277 K with two naturally occurring CCN compounds, NaCl and (NH4)2SO4, prepared by drying nebulized a 0.1% aqueous solution with an average particle mass at about 10−14 g. Equations used here are without Kelvin effect for the increase of the droplet vapor pressure and the solute effect for depression of the pressure, which tend to overestimate the effective supersaturation to the droplet. With the NaCl nucleus of 10−14 g, the maximum supersaturation of the solution droplet exists at about (S − 1) = 3.9 × 10−4 and a = 2 μm (Mason 1957). Taking an average value for the entire growth process at 75% of the maximum or effective (S − 1) = 2.9 × 10−4 with Δa = 0.162 μm, the growing droplet experiences the environmental value less this amount or 3.2 × 10−3 − 2.9 × 10−4 = 2.91 × 10−3, to alter β by a factor of 1/0.87. It was found that an increase of Δa to 0.17 μm from 0.162 μm by holding the effective (S − 1) = 2.91 × 10−3 would change β by a factor of 0.88. This is to say that [(S − 1), Δa (μm)] = (2.91 × 10−3, 0.162), the realistic condition, leads to about the same β with (3.2 × 10−3, 0.17). The data analysis was carried out based on the latter. This correction, however, was smaller than the error of measurement as will be shown later.

e. Effect of transitional steady state of a moving boundary in droplet growth kinetics and the correction

During the data analysis of measurements based on the above method and comparison with others, it was realized that an important effect was missing in the contemporary theoretical structure of small droplet growth. The effect is necessary for the steady-state vapor and temperature fields in air surrounding the droplet to adjust to new levels owing to the growth-based movement of the boundary at the droplet surface by transferring the field difference existing between them before and after the growth. Although a complete formulation of the theory is beyond the scope of this paper and will be done elsewhere, the summary derivation, and a corresponding correction factor γ for both α′ and β′, due to the transitional steady state of the moving boundary were obtained (see appendix B) as
i1520-0469-64-3-955-eq2
where
i1520-0469-64-3-955-eq3
and
i1520-0469-64-3-955-eq4
The factor, fMB, becomes very large as a reduces to and below about the 1-μm level. When the data were obtained by the method described in section 2c, the vapor flux or da/dt becomes a known constant. Since fMB affects the outer continuum vapor field only, the factor and β in the free molecular process at the droplet surface, both affect the measured vapor flux or da/dt. Based on this effect, the data of the initial measurement were reinterpreted for β and similarly for α.

3. Experimental apparatus and measurement

A double-decked, low-height, flow-type thermal diffusion chamber was constructed for the measurement (see Fig. 4). The measuring device consisted of a saturation cell followed by a diffusion cell. The walls of the cells were constructed of 0.67-cm (1/4″) copper plates for strength under low pressure operation and uniform temperature distribution. The saturation cell was 0.5 cm high, 7.6 cm wide, and 40.6 cm long. The bottom plate of the cell was covered with fiberglass filter paper, which was adhered by Dow Corning RTV Sealant silicone glue and treated with heat for degassing. It was wetted with water during the tests. The glass fiber filter paper and silicone glue combination was chosen because of the inertness of glass and to secure adhesion between them and copper plates under steam cleaning, heat treatment for degassing, and prolonged wet conditions. The saturation cell was held at the same temperature as the bottom plate of the diffusion cell. The air entered into the diffusion cell and established the steady-state profile in about 1 s. The airflow rate was typically 300 mL s−1 with the airspeed 1.2 cm s−1 in the middle of the cell.

The diffusion cell was the same dimensions as the saturation cell except that the effective height was 0.8 cm. Both the top and bottom plates were covered with glass fiber filter paper that held the water. While the rest of the system was constructed of copper, the sidewalls of the diffusion cell were fabricated from polycarbonate to provide visibility in the cell and insulation between the plates comprising the bottom and top of the cell. A 0.3-cm-diameter port covered with optical glass was placed in the end of the cell to allow a laser beam (Oriel Helium-Neon Laser 79245) to traverse the length of the cell. A second port 3.8 cm long and 0.3 cm high with a glass cover was cut in the middle of the sidewall centered at 30.5 cm from the end of the cell to obtain droplet observation in the cell. Silicone rubber “O” ring seals were used to make the cells accessible. Ports were also provided to allow the introduction of water to the top and bottom of the diffusion cell and the bottom of the saturation cell, with outlet ports to remove excess water. The system was tilted at an angle of approximately 1° to ensure water would flow down the fiberglass filter paper, keeping it wet. The water supply and removal ports were sealed during the tests at reduced pressure. The top plate temperature was controlled by a Neslab Endocal Model ULT-80 low temperature bath circulator. The bottom plate and the saturation cell temperature control utilized a Forma Scientific Model 2160 bath and circulator. The temperature of the cell was measured with copper-constantan thermocouples in ports drilled in the middle of upper and lower plates. These were located 25 cm from the end of the cell. The thermocouple readings were taken using a Yokogawa Model LR421OE recorder. Prior to use, the entire system was steam cleaned for approximately 1 h, followed by heating to ∼80°C on a hotplate over a weekend with a slow flow of air passing through the system.

The introduction of nucleus samples was made through an 18-gauge (0.125 cm O.D., approximately 0.076 cm I.D.) hypodermic needle that had a hole drilled in the side and the end plugged. The needle port was located at 20.3 cm from the inlet of the cell. The needle was positioned so that the hole was in the center of the cell and the laser beam. With extreme care, sample introduction provided a thin stream of nuclei in low concentration, and droplets were not influencing the growth of adjacent ones to allow the identification of individual droplets. The introduction was accomplished in two ways: Nuclei were produced using a small solution nebulizer blowing into a 250 mL Erlenmyer flask that held clean, dry air. For the room pressure tests, samples of the flask contents were taken with a 10-mL syringe and very slowly introduced into the system through the needle. An alternative of this method was to let the nucleus sample be sucked into the chamber air under slight negative pressure. The sample air was estimated to reach the supersaturation of the surrounding air in about 2 × 10−2 s. For the low pressure tests, the flask containing the suspended nuclei was connected to the introduction needle by a small silicon rubber tube that had a variable restrictor in the line as well as an open or closed valve. The pressure in the flask was reduced somewhat and the restrictor adjusted to allow a very small leakage into the needle. The use of the valve enabled a very small amount of sample to be injected by the difference in pressure between the system and the sample flask by quickly opening and closing the valve, usually much less than one second during this process. The flashing points of red laser light by the droplet growth were recorded using a Sony Model No. DC-TRV-18 digital camera recorder and analyzed using Pinnacle Studio DV, Version 7 software.

The air supply to the system was prepared by simply passing the room air through an electrostatic precipitator to remove the existing particles after treatment of the carrier room air by a liquid nitrogen trap, bubbling through an oxidizer bath followed by electrostatic particle removal or use of nitrogen gas from the liquid reservoir. No recognizable difference in condensation coefficient was observed from that measured in air treated through the electrostatic particle precipitator alone. In a separate estimation, we found that there did not exist a sufficient amount of adsorbable organic molecules to form a monolayer on the surface of a formed droplet for growth impediment (Fuentes et al. 2000). Airflow was obtained using a dry vane vacuum pump and controlled by needle valves placed after the system for the atmospheric pressure runs and before the system for the reduced pressure tests.

Deionized water was used in the tests for keeping the filter paper wet and for making solutions. The water production system had a charcoal trap followed by a reverse osmosis system, which, in turn, was followed by mixed ion exchange resin beds. All tubing used in the system was Teflon or silicon, which gave no volatiles and contained no plasticizers.

Analysis of the data was done by observing the position, distance, and time elapsed between two adjacent Mie scattering peaks or flashes for a given particle. The time between the two flashes gave Δa = 0.162 μm. The airflow rate was given by the distance to the middle of two flashes and the time between them. The distance to the middle of the two flash positions gave the average time from nuclei injection [cf. Eq. (2)] using the flow speed to convert distance to time. The average droplet radius when measured was ∼2 μm. The temperature at the center of the cell was approximately 277 K, with a temperature difference across the cell of ∼2.5°C. Particles introduced at the top of the beam would fall a maximum of ∼2 mm during their traverse from the injection point to the end of the observation. They were kept within the fall distance limit by controlling the flow speed.

4. Results of measurement and discussion

In the course of the preliminary test, under some conditions that we never pursued nor identified, the scattered beams from the He–Ne laser were observed colored other than the original red. The reason for this was not known to us. It appeared that possibly harmonics of the laser light occurred in the scattering. The original red laser beam did scatter at the position of the measurement.

a. Measured data and the mechanism of α and β appearance

The measurements were first carried out based on the pressure modulation method of the diffusion-kinetic theory (FW70; FX96). As the theoretical structure applied would play a crucial role, the necessity of solving the unexplored moving boundary effect became clear because it was closely associated with the main, longest lasting steady-state process and was as durable as the latter. The theoretical derivation was successfully carried out, as briefly shown in appendix B, and a correction factor [Eq. (B8)] was obtained and applied to the result of the initial analysis of the diffusion-kinetic theory alone without the moving boundary effect. The moving boundary effect increases the droplet growth rate and gives an overestimate for α and β. It is supersaturation dependent. Under the low supersaturation of the present measurement, a small correction was given; ∼−16% for β and ∼−10% in α.

Results of 190 measurements carried out under a combination of conditions at around T = 277 K, (S − 1) = 0.32%, p1 = 100 hPa, p2 = 860 hPa, a ≃ 2 μm with nuclei of NaCl, (NH4)2 SO4, and paper smoke (Kimwipe) with and without correction of the moving boundary effect are listed in Table 1 showing some ranges of observational errors at 99th percentile of total runs. With corrections of Kelvin, solute, and moving boundary effects, α = 0.81 and β = 0.043, averaged between two typical CCN, NaCl, and (NH4)2SO4, should represent the pure water surface under the conditions of measurement: β for paper smoke does give a value slightly smaller than the above. However, in general, impurities such as organic molecules do not necessarily adsorb at the droplet surface to lower β. Furthermore, existing concentrations of organic molecules (Fuentes et al. 2000) are insufficient to form a monolayer at the surface by adsorption.

The α of the droplet growth system represents the efficiency of exchanging sensible heat between noncondensable air molecules and the liquid water surface. On the other hand, β is the efficiency of condensation for condensable water molecules at the liquid surface with the release and the removal of the latent heat of condensation. Without the removal of heat, the β process becomes the same as that of α.

Orientation of the surface layer or the presence of adsorbed molecules may intervene in the condensation process, as evidenced with smaller β in paper smoke. Nevertheless, such an intervention is also expected in the α process. The nearly 1.5 orders of magnitude difference between α and β in the table, with the former near unity, clearly suggests a significant difference in the mechanism, that is, the dissipation of the released latent heat in the β process through conduction to the surroundings. The amount of the latent heat released upon the molecular capture is about 20 times more than the kinetic energy involved in a specular reflection. Furthermore, this heat removal is a rate process and takes some finite amount of time. If the “residence” time of impinging water molecules is much shorter than that of the heat removal, the molecules will reflect back, reducing the capture rate or β. The reverse involving the “evaporation” coefficient is expected to hold for molecules that are about to leave the water surface by gaining the necessary kinetic energy from the surroundings by heat conduction to give the physical basis for equality of β and the evaporation coefficient.

For prenucleation clusters with small numbers of molecules, the β values are expected to fall even further. In the new homogeneous condensation nucleation theory with a population of molecular clusters, it was indeed the necessary condition to describe the stability of the supersaturated phase below the nucleation threshold (Fukuta 2004).

Under the influence of prior works, Pruppacher and Klett (1978) classified the measured data into two groups: one was based on evaporation from a quasi-quiescent water surface that reported lower average values of β and the other from a rapidly renewing water surface with higher average values. Physical basis for this classification does not exist; regardless of the surface renewing or not, the water surface, even under equilibrium, is always rough due to molecular bombardment.

b. The moving boundary effect of droplet growth and data comparison

The theoretical development of the moving boundary effect of droplet growth happened to be, as mentioned earlier, posterior to the present measurement. It appears to close the last large remaining gap in the growth kinetics of droplet and to provide the missing link among the existing results of measurements. There are two typical experimental conditions: One has the moving boundary effect attached to the main diffusion-kinetic process expressed in Eq. (B7). This condition represents the growth process with no expansion or continuous expansion like the growing cloud droplets in the updraft. The continuous expansion process holds steady-state growth by continuous generation of supersaturation within the cell, the space claimed by a droplet, until the supersaturation generation rate becomes overcome by the increased removal rate with larger droplets and the steady-state field collapses. The other is the growth under instantaneous expansion in which the established boundary of the cell permits neither mass nor heat flux to cross, and the process also does not allow the continuous generation of supersaturation. Failing to establish diffusion-kinetic growth in the instantaneous expansion process, (1 + fMB) of the no-expansion process in Eq. (B7) goes to fMB.

Figure 5 summarizes the moving boundary effect in two typical cases for no expansion and instantaneous expansion in the ratio
i1520-0469-64-3-955-e11
In the figure, the moving boundary effect increases RMW for (S − 1) from about a = 10−5 to 10−2 cm. For both no expansion and instantaneous expansion cases, when (S − 1) ≥ 0.3, the ratio approaches and exceeds the Maxwellian growth rate ratio, RMW = 1 before falling back. This effect is clearly visible in the unsuspected expansion experiment of Vietti and Schuster (1973) in comparison with FW70 under very high supersaturations. It is important to point out that the moving boundary effect on the diffusion-kinetic process, due to the growth rate increase, thus tends to give a deceptive impression as if the growth behavior were truly Maxwellian. For (S − 1) = 0.0032, the value taken in the present measurement, the RMW increase from the diffusion-kinetic process is sufficiently small. However, for the upper limit of natural cloud processes, say (S − 1) ≈ 0.032, the deviation becomes quite recognizable. This is likely to enhance the broadening effect of the droplet size spectrum by increasing the growth rate on the larger side, in addition to another broadening effect due to the growth rate reduction of the diffusion-kinetic theory on the smaller end, the latter being the original contention of FW70 as mentioned in Hagen et al. (1989). It may enhance early formation of large droplets to let the collision–coalescence process proceed toward raindrop generation or graupel and hail processes, especially in the convective clouds of high updraft.

From the above viewpoint of the moving boundary effect riding on the diffusion-kinetic process, Fig. 5 also suggests the experimental condition most suitable for α and β determination; the moving boundary effect is sufficiently suppressed while a well-established diffusion-kinetic process prevails. This is to say that measurement in the no-expansion condition like a thermal diffusion chamber with small (S − 1), the moving boundary effect and droplet size for enhanced diffusion-kinetic effect, still being outside the Kelvin and solute effects, should be advantageous over the others to avoid a large correction.

In comparison of the present result with others, the common demographic approach in the values of α and β does not solve the problem and is even misleading. Reported data are influenced by the method of measurement as well as the theory applied in the analysis. Their choice is crucial to the comparison. As far as the measurement methods are concerned, the one that employs droplet growth has a distinct advantage with measurable macroscopic environmental properties, which are continuously connected to the droplet surface, the place where α and β effects appear. For this reason and the process being a direct reproduction of what happens in the clouds, here data of droplet growth-based measurements are selected for comparison. For the theory of data analysis, almost all of the past measurements applied the diffusion-kinetic theory of FW70 and FX96, obviously without consideration of the moving boundary effect. The moving boundary effect is hence the key point for the comparison, and the latter in turn may test the effect. The selected data of measurements are categorized into three groups based on the type of chambers or the method used as discussed below.

1) Expansion chamber

The chamber was used in two different modes—instantaneous expansion and continuous expansion. In the instantaneous expansion mode of operation applying the Mie scattering method of the He–Ne laser for droplet size determination under S = 1.49–3.45, Vietti and Schuster (1973) obtained β′ = 0.0065 with the α′ = 1 assumption.2 A very large positive deviation of droplet radius from that of the diffusion-kinetic theory was observed, which later switched into a negative. The reported β′ value is therefore questionable but is consistent and explainable with the behavior of the moving boundary effect under high (S − 1) in Fig. 5.

The instantaneous expansion mode was also used in the recent work of Winkler et al. (2004), under the condition of T = 250–290 K, S = 1.3–1.5 with droplet size determination of the Mie scattering of laser and data analysis based on a curve-fitting method for simultaneous determination of α′ and β′ under pressure modulation. The Maxwellian behavior of α′ = β′ = 1 was reported. Whereas, the measured data were showing da/dt ≠ 0 at t = 0, the diffusion-kinetic behavior, as shown in Fig. 1. The high (S − 1) of the measurement most likely placed the condition in the zone close to RMW = 1 for moving boundary effect in the instantaneous expansion mode (see Fig. 5).

The results of extensive expansion chamber measurement in the continuous expansion mode were reported by Hagen et al. (1989) using the Mie scattering method with an argon continuous wave (CW) laser for the droplet size detection with the starting T = 16°C, cooling rate of 10°C min−1 corresponding to 17 m s−1 cloud-base updraft, which would reach nominal (S − 1) = 0.3 in 30 s in our estimate. The reported maximum supersaturation was around 5% or 6%, and droplets were measured to a = 12–15 μm. Droplet growth was found to be more rapid during the first test, possibly indicating impurity accumulation in the following test. The main result is that, under the α′ = 1 assumption, β′ was found to decrease in the neighborhood of 0.01 from a value near unity at the outset of the test. Nevertheless, with a careful examination of the reported data, one can find a clear increase of β′ with a size in the range below a = 2–5 μm. Furthermore, for the size range, data points may be seen crammed at the β′ = 1 level, an indication of growth exceeding the Maxwellian behavior. These behaviors are unexplainable with the purity or lack of accumulation of contamination in the early stage of the growth. They are, however, compatible with the moving boundary effect of the diffusion-kinetic theory.

In Fig. 5, the upward deviation of droplet growth by the moving boundary effect from the diffusion-kinetic line, with α′ and β′ of our present measurement, is shown by iso-(S − 1) lines and other diffusion-kinetic lines with larger α′ and β′ (not shown) lying above. The maxima of the moving boundary effect for no expansion that correspond to a = 2–5 μm may be found with the (S − 1) = 32% line and can be seen in a portion of the line exceeding RMW = 1 or the Maxwellian. The zone of (S − 1) = 5%–6% should exist close to this. Under continuous expansion and increasing (S − 1) before the maximum, a droplet grows with the moving boundary effect and passes through the iso-(S − 1) lines to approach the Maxwellian and give a β′ increase if interpreted by the diffusion-kinetic theory alone and vice versa as it grows past the maximum. This tendency always exists with growth along (S − 1) = const.

We conclude that, while some contamination could have existed in this work to reduce the β′ value in the later stage of measurement, β′ near unity at the beginning of the growth and the steady decrease thereafter toward 0.01 are the deceptive results caused by data reduction with diffusion-kinetic theory alone without the moving boundary effect. The decrease may also be attributable to the loss of the supersaturated vapor to the wall during the rather prolonged measurement.

2) Thermal diffusion chamber

Three operational modes exist with the thermal diffusion chamber; the static mode and the vertical, as well as the horizontal, flow mode. Chodes et al. (1974) made use of the thermal diffusion chamber in the static mode at T ≈ 20°C, (S − 1) = 0.0040–0.0075, with the droplet size directly measured by the fall velocity in the chopped light beam, to report β′ = 0.033 on average. The released latent heat of condensation during growth of the droplet could have slowed down the fall velocity to give a slight smaller β′. Phoretic effects were not considered but are expected to be small. The moving boundary effect to increase the growth rate and therefore β′ is also expected to be small due to small (S − 1). The reported data should thus be regarded close to the true value.

Sinnarwalla et al. (1975), on the other hand, applied the chamber vertically with flow at (S − 1) = 0.005–0.01, T = 20°C, and photographic determination of the fall velocity and the size, to obtain β′ = 0.026 under the α′ = 1 assumption. Constraints here are basically the same as those of Chodes et al. (1974), and the reported β′ can be considered sufficiently close to reality. The present work is an example of α′ and β′ measurements that make use of a horizontal flow thermal diffusion chamber.

3) Electrodynamic levitation

The electrodynamic levitation method uses a charged particle, few tens of micrometers in diameter, which happened to be in the zone of low sensitivity for the diffusion-kinetic effect (cf. Fig. 2) and also low moving boundary effect (cf. Fig. 5), suspended in the quadrupole, and through the fall behavior with balancing voltage, to measure growth or evaporation. Sageev et al. (1986) measured the growth rate of an aqueous (NH4)2SO4 solution drop of diameter between 10 and 20 μm to report α′ ≈ 1. At T = −34.9°C, Shaw and Lamb (1999) could determine the possible range for α′ = 0.1–1 and β′ = 0.04–0.1 applying the pressure modulation method for evaporation of a = 20 μm drop. While the accuracy of these measurements is low, the results are generally in line with those of the thermal diffusion chamber method.

5. Conclusions

A new droplet growth theory of moving boundary effect with second-order (S − 1) dependence was derived and found to play an important role as a “missing link” between the two controversial groups of data for α, β measurements. The theory explained and closed the gap. In turn, the data confirmed validity for the moving boundary effect of the theory. The high (S − 1) of expansion chamber methods, not the rapid renewal of water surface, was found to enhance the moving boundary effect for faster growth and lead to deceptive interpretation toward larger α, β. The theory provided the basis for bettering the measurement, correction for data analyzed without the moving boundary effect. It suggests a new broadening effect in the droplet size spectrum by accelerating the growth on the larger end together with the diffusion-kinetic effect of growth slowdown on the smaller, toward onset of the collision–coalescence process and precipitation development.

Measurements employing a flow-type horizontal thermal diffusion chamber and the Mie scattering method using a He–Ne laser at T = 277 K, (S − 1) = 0.32%, a ≈ 2 μm under two pressures of modulation, p = 100 and 860 hPa, with correction of the moving boundary effect for nuclei of NaCl and (NH4)2S04, resulted in the simultaneous determination of α = 0.81 and β = 0.043 on average. For paper smoke, α = 0.68 and β = 0.022.

Comparison with other data of low moving boundary effect suggests a slight β increase with temperature lowering. Based on the result of the present measurement and other theoretical and experimental works, the mechanism of β being much smaller than α is suggested mostly as a kinetic process controlled by the removal rate of the released latent heat at the moment of molecular impact at the surface for β and without the process for α. An identical situation applies to evaporation.

Acknowledgments

This work was partially supported by NSF under Grant ATM-0094909.

REFERENCES

  • Carslaw, H. S., and J. C. Jaeger, 1959: Conduction of Heat in Solids. Oxford University Press, 510 pp.

  • Chapman, S., and T. G. Cowling, 1970: The Mathematical Theory of Non-Uniform Gases. Cambridge University Press, 422 pp.

  • Chodes, N. J., J. Warner, and A. Gagin, 1974: A determination of the condensation coefficient of water from the growth rate of small cloud droplets. J. Atmos. Sci., 31 , 13511357.

    • Search Google Scholar
    • Export Citation
  • Fuentes, J. D., and Coauthors, 2000: Biogenic hydrocarbons in the atmospheric boundary layer: A review. Bull. Amer. Meteor. Soc., 81 , 15371575.

    • Search Google Scholar
    • Export Citation
  • Fukuta, N., 1992: Theories of competitive cloud droplet growth and their application to cloud physics studies. J. Atmos. Sci., 49 , 11071114.

    • Search Google Scholar
    • Export Citation
  • Fukuta, N., 2004: Restructuring nucleation kinetics based on new size-dependent surface free energy. Nucleation and Atmospheric Aerosols 2004, M. Kasahara and M. Kulmala, Eds., Kyoto University Press, 178–181.

    • Search Google Scholar
    • Export Citation
  • Fukuta, N., and L. A. Walter, 1970: Kinetics of hydrometeor growth from a vapor-spherical model. J. Atmos. Sci., 27 , 11601172.

  • Fukuta, N., and N. Xu, 1996: Nucleation-droplet growth interactions and microphysical property development in convective clouds. Atmos. Res., 41 , 122.

    • Search Google Scholar
    • Export Citation
  • Fukuta, N., and C. M. Gramada, 2003: Vapor pressure measurement of supercooled water. J. Atmos. Sci., 60 , 18711875.

  • Grinin, A. P., A. K. Shchekin, F. M. Kuni, E. A. Grinina, and H. Reiss, 2004: Non-steady profile of vapor concentration around the growing droplet: Balance of the condensing vapor and the moving boundary of the droplet. Nucleation and Atmospheric Aerosols 2004, M. Kasahara and M. Kulmala, Eds., Kyoto University Press, 375–379.

    • Search Google Scholar
    • Export Citation
  • Hagen, D. E., J. Schmitt, M. Trueblood, J. Carstens, D. R. White, and D. J. Alofs, 1989: Condensation coefficient measurement for water in the UMR cloud simulation chamber. J. Atmos. Sci., 46 , 803816.

    • Search Google Scholar
    • Export Citation
  • Katz, J. L., and P. Mirabel, 1975: Calculation of supersaturation profile in thermal diffusion cloud chambers. J. Atmos. Sci., 32 , 646648.

    • Search Google Scholar
    • Export Citation
  • Knudsen, M., 1934: The Kinetic Theory of Gases: Some Modern Aspects. John Wiley & Sons, 64 pp.

  • Mason, B. J., 1957: The Physics of Clouds. Oxford University Press, 481 pp.

  • Maxwell, J. C., 1890: Theory of the wet bulb thermometer. The Scientific Papers of James Clerk Maxwell, Vol. 2, W. D. Niven, Ed., Dover Publications, 636–640.

    • Search Google Scholar
    • Export Citation
  • Mozurkewich, M., 1986: Aerosol growth and the condensation coefficient for water: A review. Aerosol Sci. Technol., 5 , 223236.

  • Pruppacher, H. R., and J. D. Klett, 1978: Microphysics of Clouds and Precipitation. D. Reidel, 714 pp.

  • Pruppacher, H. R., and J. D. Klett, 1997: Microphysics of Clouds and Precipitation. Kluwer Academic, 954 pp.

  • Sageev, G., R. C. Flagan, J. H. Seinfeld, and S. Arnold, 1986: Condensation rate of water on aqueous droplets in the transition regime. J. Colloid Interface Sci., 113 , 421429.

    • Search Google Scholar
    • Export Citation
  • Shaw, R. A., and D. Lamb, 1999: Experimental determination of the thermal accommodation and condensation coefficients of water. J. Chem. Phys., 111 , 1065910663.

    • Search Google Scholar
    • Export Citation
  • Sinnarwalla, A. M., D. J. Alofs, and J. C. Carstens, 1975: Measurements of growth rate to determine condensation coefficients for water drops growth on natural cloud nuclei. J. Atmos. Sci., 32 , 592599.

    • Search Google Scholar
    • Export Citation
  • Smirnov, V. I., 1971: The rate of quasi-steady growth and evaporation of small drops in a gaseous medium. Pure Appl. Geophys., 86 , 187194.

    • Search Google Scholar
    • Export Citation
  • Twomey, S., 1977: Atmospheric Aerosols. Elsevier, 302 pp.

  • Vietti, M. A., and B. G. Schuster, 1973: Laser scattering measurement of droplet growth in binary mixtures: Part I. H2O and air. J. Chem. Phys., 58 , 434441.

    • Search Google Scholar
    • Export Citation
  • Winkler, P. M., A. Vrtala, P. E. Wagner, M. Kulmala, K. E. J. Lehtanen, and T. Versala, 2004: An experimental study on thermal and mass accommodation coefficients for the condensation of water vapour. Nucleation and Atmospheric Aerosols 2004, M. Kasahara and M. Kulmala, Eds., Kyoto University Press, 144–146.

    • Search Google Scholar
    • Export Citation
  • Zou, Y-S., and N. Fukuta, 1999: The effect of diffusion kinetics on the supersaturation in clouds. Atmos. Res., 52 , 115141.

APPENDIX A

List of Symbols and Abbreviations

  • a  Droplet radius

  • a′  Term defined by Eq. (1)

  • A  Term defined by Eq. (9)

  • b′  Term defined by Eq. (1)

  • B  Term defined by Eq. (9)

  • C  Constant

  • CCN  Cloud condensation nuclei

  • D  Diffusivity of water vapor in air

  • D0  D at standard temperature and pressure

  • [DK]  Square bracket term in Eq. (1)

  • e  Vapor pressure of water

  • E  Term defined by Eq. (9)

  • fc  Correction factor for cell boundary controlled kinetics (Fukuta 1992)

  • fα  Term defined by Eq. (1)

  • fβ  Term defined by Eq. (1)

  • fMB  The moving boundary factor defined by Eq. (B6)

  • FW70  Fukuta and Walter (1970)

  • FX96  Fukuta and Xu (1996)

  • g  Gravity acceleration

  • K  Thermal conductivity of air

  • l  Mean free path of air

  • lα  Term defined by Eq. (1)

  • lβ  Term defined by Eq. (1)

  • L  Specific latent heat of condensation

  • m  Mass of water droplet

  • n  Number density of water molecules

  • p  Air pressure

  • p0  Standard air pressure

  • r  Radial distance

  • R  Radial distance to the cell boundary

  • Ra  Specific gas constant of air

  • Rυ  Specific gas constant of water vapor

  • RMW  Ratio defined by Eq. (11)

  • S  Saturation ratio

  • (S − 1)  Supersaturation

  • t  Time

  • T  Temperature

  • T0  Standard temperature

  • Ta  Temperature of air

  • w  Vertical velocity of gas molecules

  • wd  Fall velocity of droplet

  • z  Vertical distance

  • α  Thermal accommodation coefficient

  • α′  Representative thermal accommodation coefficient defined by Eq. (1)

  • β  Condensation coefficient

  • β′  Representative condensation coefficient defined by Eq. (1)

  • γ  Correction factor to α′ and β′ of diffusion-kinetic estimation

  • ρ  Density of water vapor

  • ρs  Saturation density of water vapor

  • ρp  Density of particle

  • ρΔ  Density of water vapor just above the droplet

  • ρ  Density of water vapor in the environment

  • η  Dynamic viscosity of air

  • λ  Wavelength of light

APPENDIX B

Transitional Steady-State Kinetics of Droplet Growth with the Moving Surface Boundary and the Correction Factor

The growth of a cloud droplet involves diffusion of water vapor to and conduction of the released latent heat away from it. These transport phenomena occur generally in two processes, the changing or non–steady state and the true steady state (Carslaw and Jaeger 1959). The former processes are short lived, while the latter remain as long as the necessary condition holds. The contemporary theories of droplet growth, such as FW70, are formulated based on the assumption of the long-lasting steady state. Nevertheless, when the droplet growth takes place, the surface or the inner boundary of the vapor or temperature field moves outward (the problem of moving boundary), violating the condition necessary to maintain true steady state. The outward movement of growing droplet surface causes shifts in both the vapor and temperature fields in the larger surrounding air space, and the respective differences appearing in the processes must be transferred to and from the droplet by altering their gradients there. The process is essentially in nonsteady state but is closely attached to the ongoing steady-state process. So, it is termed here as the transitional steady-state process of moving boundary. This problem has been known for sometime. However, the attempt to solve it has never been successful, for the starting process of heat conduction in solids for modification deals with no boundary movement and, therefore, no space exists for treatment at least for the transitional steady state–moving boundary aspect (Grinin et al. 2004).

The diffusion-kinetic theory of droplet growth in steady state consists of the following two parts: the continuous process of vapor or heat transfer in the air space and the free molecular processes at the droplet surface. The aforementioned transitional steady-state process of a moving boundary directly affects the former, and in the solution, the effect is conveyed simultaneously to the later.

Figure B1 illustrates the droplet growth system of the transitional steady state–moving boundary effect. For the reason stated above, the excess vapor density to be removed from the steady-state solution for da growth in dt time is (cf. FW70)
i1520-0469-64-3-955-eb1
The inward mass transfer across the surface at a distance r from the droplet for dt time is then expressed by integrating Eq. (B1) from r to R,
i1520-0469-64-3-955-eb2
Whereas, equating Fick’s law of diffusion to the drop volume increase rate with ρL or the mass growth rate in steady state at the droplet surface gives
i1520-0469-64-3-955-eb3
where, the steady-state solution of FW70,
i1520-0469-64-3-955-eb4
is used, and [DK] is the square bracket term in Eq. (1).
The steady-state form of Eq. (B2) for r = a is
i1520-0469-64-3-955-eb5
Then, combination of Eqs. (B2) and (B5) with Eqs. (B3) and (B4) yields the mass flux rate ratio between the moving boundary and steady state or the moving boundary correction factor,
i1520-0469-64-3-955-eb6
The chief contributors to fMB are thus (S − 1) and a.
The correction factor, fMB combined with the preexisting steady-state flux increases the steady-state continuous flow of vapor and heat in a manner identical to fc of the cell boundary controlled kinetics of droplet growth (Fukuta 1992) which also involves steady state and transitional steady-state flows. Then, for both steady state and transitional steady state of moving boundary processes occurring simultaneously, Eq. (1) changes to
i1520-0469-64-3-955-eb7
where 1 + fMB = fc since fMB is the fraction to the diffusion-kinetic process.
In the initial analysis of measurement, Eq. (1) was used to estimate α′ and β′, while Eq. (B7) takes account of an additional transitional steady state–moving boundary process and, therefore, should be more accurate and preferable for the application. Equations (1) and (B7), when compared in the same mass or radial growth-rate form, differ in the square bracket of the denominator. However, in the analysis both are equated to the observed data. The square bracket term of Eq. (1) is constant under a given condition of the initial observation. Equation (B7), which is the correct form, has reduced the first term in the bracket of Eq. (1) by (1 + fMB) > 1. To make both bracket values equal, an increase of the second term in the same bracket is necessary. Since α′ and β′ estimated based on Eq. (1) are to be corrected and appear both in the denominators of lα and l β, a common and direct correction factor, γ, is introduced to them. Then,
i1520-0469-64-3-955-eqb1
or
i1520-0469-64-3-955-eb8
From the product of γ and α′ and β′ of Eq. (1), correct α and β may be obtained.

Fig. 1.
Fig. 1.

Comparison of droplet radius, a, between Maxwellian and diffusion-kinetic theories as a function of time, t.

Citation: Journal of the Atmospheric Sciences 64, 3; 10.1175/JAS3834.1

Fig. 2.
Fig. 2.

Ratios of droplet radius, a, and mass growth rate, dm/dt, between the Maxwellian (MW) and diffusion-kinetic (DK) growths plotted as a function of the Maxwellian radius, aMW.

Citation: Journal of the Atmospheric Sciences 64, 3; 10.1175/JAS3834.1

Fig. 3.
Fig. 3.

Mie scattering peaks at 90° for the He–Ne laser (λ = 0.6328 μm) plotted as a function of droplet radius (from Oregon Medical Laser Center, http://omlc.ogi.edu/calc/mie_calc.html).

Citation: Journal of the Atmospheric Sciences 64, 3; 10.1175/JAS3834.1

Fig. 4.
Fig. 4.

The double-decked, flow-type thermal diffusion chamber used in the present measurement.

Citation: Journal of the Atmospheric Sciences 64, 3; 10.1175/JAS3834.1

Fig. 5.
Fig. 5.

The mass growth rate ratio between the diffusion-kinetic with the moving boundary effect and the Maxwellian, RMW, plotted as a function of droplet radius, a, for different supersaturation, (S − 1): NE stands for no expansion and instantaneous expansion for instantaneous expansion, respectively.

Citation: Journal of the Atmospheric Sciences 64, 3; 10.1175/JAS3834.1

i1520-0469-64-3-955-fb01

Fig. B1. The system of a growing droplet in transitional steady state under the influence of moving surface boundary.

Citation: Journal of the Atmospheric Sciences 64, 3; 10.1175/JAS3834.1

Table 1.

List of the condensation coefficient, β, and the thermal accommodation coefficient, α, measured for three different CCN compounds at temperature T ≈ 277 K and supersaturation (S − 1) = 0.0032 with and without the moving boundary effect corrections.

Table 1.

1

The factor of 2 increase in the thermal (and vapor) flux gradient in appendix B of Hagen et al. (1989) and in Mozurkewich (1986) is incorrect because the gradient at the droplet surface is based on l distance, while in the ordinary gas kinetic formulation, it is on 2l distance. The Stefan flow correction cited and ignored in Hagen et al. (1989) is a misconception. The flow was thought to be hydrodynamic but is, in effect, diffusional. Hydrodynamic flow is given by a shared average velocity component between the vapor and air in the direction of the flow but, in diffusion across the droplet surface, no shared component exists since an average air velocity is zero.

2

Constants α and β used in most works are α′ and β′ in the present formulation, and changes are made here accordingly.

Save
  • Carslaw, H. S., and J. C. Jaeger, 1959: Conduction of Heat in Solids. Oxford University Press, 510 pp.

  • Chapman, S., and T. G. Cowling, 1970: The Mathematical Theory of Non-Uniform Gases. Cambridge University Press, 422 pp.

  • Chodes, N. J., J. Warner, and A. Gagin, 1974: A determination of the condensation coefficient of water from the growth rate of small cloud droplets. J. Atmos. Sci., 31 , 13511357.

    • Search Google Scholar
    • Export Citation
  • Fuentes, J. D., and Coauthors, 2000: Biogenic hydrocarbons in the atmospheric boundary layer: A review. Bull. Amer. Meteor. Soc., 81 , 15371575.

    • Search Google Scholar
    • Export Citation
  • Fukuta, N., 1992: Theories of competitive cloud droplet growth and their application to cloud physics studies. J. Atmos. Sci., 49 , 11071114.

    • Search Google Scholar
    • Export Citation
  • Fukuta, N., 2004: Restructuring nucleation kinetics based on new size-dependent surface free energy. Nucleation and Atmospheric Aerosols 2004, M. Kasahara and M. Kulmala, Eds., Kyoto University Press, 178–181.

    • Search Google Scholar
    • Export Citation
  • Fukuta, N., and L. A. Walter, 1970: Kinetics of hydrometeor growth from a vapor-spherical model. J. Atmos. Sci., 27 , 11601172.

  • Fukuta, N., and N. Xu, 1996: Nucleation-droplet growth interactions and microphysical property development in convective clouds. Atmos. Res., 41 , 122.

    • Search Google Scholar
    • Export Citation
  • Fukuta, N., and C. M. Gramada, 2003: Vapor pressure measurement of supercooled water. J. Atmos. Sci., 60 , 18711875.

  • Grinin, A. P., A. K. Shchekin, F. M. Kuni, E. A. Grinina, and H. Reiss, 2004: Non-steady profile of vapor concentration around the growing droplet: Balance of the condensing vapor and the moving boundary of the droplet. Nucleation and Atmospheric Aerosols 2004, M. Kasahara and M. Kulmala, Eds., Kyoto University Press, 375–379.

    • Search Google Scholar
    • Export Citation
  • Hagen, D. E., J. Schmitt, M. Trueblood, J. Carstens, D. R. White, and D. J. Alofs, 1989: Condensation coefficient measurement for water in the UMR cloud simulation chamber. J. Atmos. Sci., 46 , 803816.

    • Search Google Scholar
    • Export Citation
  • Katz, J. L., and P. Mirabel, 1975: Calculation of supersaturation profile in thermal diffusion cloud chambers. J. Atmos. Sci., 32 , 646648.

    • Search Google Scholar
    • Export Citation
  • Knudsen, M., 1934: The Kinetic Theory of Gases: Some Modern Aspects. John Wiley & Sons, 64 pp.

  • Mason, B. J., 1957: The Physics of Clouds. Oxford University Press, 481 pp.

  • Maxwell, J. C., 1890: Theory of the wet bulb thermometer. The Scientific Papers of James Clerk Maxwell, Vol. 2, W. D. Niven, Ed., Dover Publications, 636–640.

    • Search Google Scholar
    • Export Citation
  • Mozurkewich, M., 1986: Aerosol growth and the condensation coefficient for water: A review. Aerosol Sci. Technol., 5 , 223236.

  • Pruppacher, H. R., and J. D. Klett, 1978: Microphysics of Clouds and Precipitation. D. Reidel, 714 pp.

  • Pruppacher, H. R., and J. D. Klett, 1997: Microphysics of Clouds and Precipitation. Kluwer Academic, 954 pp.

  • Sageev, G., R. C. Flagan, J. H. Seinfeld, and S. Arnold, 1986: Condensation rate of water on aqueous droplets in the transition regime. J. Colloid Interface Sci., 113 , 421429.

    • Search Google Scholar
    • Export Citation
  • Shaw, R. A., and D. Lamb, 1999: Experimental determination of the thermal accommodation and condensation coefficients of water. J. Chem. Phys., 111 , 1065910663.

    • Search Google Scholar
    • Export Citation
  • Sinnarwalla, A. M., D. J. Alofs, and J. C. Carstens, 1975: Measurements of growth rate to determine condensation coefficients for water drops growth on natural cloud nuclei. J. Atmos. Sci., 32 , 592599.

    • Search Google Scholar
    • Export Citation
  • Smirnov, V. I., 1971: The rate of quasi-steady growth and evaporation of small drops in a gaseous medium. Pure Appl. Geophys., 86 , 187194.

    • Search Google Scholar
    • Export Citation
  • Twomey, S., 1977: Atmospheric Aerosols. Elsevier, 302 pp.

  • Vietti, M. A., and B. G. Schuster, 1973: Laser scattering measurement of droplet growth in binary mixtures: Part I. H2O and air. J. Chem. Phys., 58 , 434441.

    • Search Google Scholar
    • Export Citation
  • Winkler, P. M., A. Vrtala, P. E. Wagner, M. Kulmala, K. E. J. Lehtanen, and T. Versala, 2004: An experimental study on thermal and mass accommodation coefficients for the condensation of water vapour. Nucleation and Atmospheric Aerosols 2004, M. Kasahara and M. Kulmala, Eds., Kyoto University Press, 144–146.

    • Search Google Scholar
    • Export Citation
  • Zou, Y-S., and N. Fukuta, 1999: The effect of diffusion kinetics on the supersaturation in clouds. Atmos. Res., 52 , 115141.

  • Fig. 1.

    Comparison of droplet radius, a, between Maxwellian and diffusion-kinetic theories as a function of time, t.

  • Fig. 2.

    Ratios of droplet radius, a, and mass growth rate, dm/dt, between the Maxwellian (MW) and diffusion-kinetic (DK) growths plotted as a function of the Maxwellian radius, aMW.

  • Fig. 3.

    Mie scattering peaks at 90° for the He–Ne laser (λ = 0.6328 μm) plotted as a function of droplet radius (from Oregon Medical Laser Center, http://omlc.ogi.edu/calc/mie_calc.html).

  • Fig. 4.

    The double-decked, flow-type thermal diffusion chamber used in the present measurement.

  • Fig. 5.

    The mass growth rate ratio between the diffusion-kinetic with the moving boundary effect and the Maxwellian, RMW, plotted as a function of droplet radius, a, for different supersaturation, (S − 1): NE stands for no expansion and instantaneous expansion for instantaneous expansion, respectively.

  • Fig. B1. The system of a growing droplet in transitional steady state under the influence of moving surface boundary.

All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 572 400 40
PDF Downloads 159 44 2