Thayer Watkins Silicon Valley & Tornado Alley USA |
---|
Due to Diffusion |
In cloud physics there is an equation, usually called the growth
equation, which describes how a cloud droplet grows. Rogers and
Yau in their A Short Course in Cloud Physics provide a derivation
of this equation. Although overall Rogers and Yau's text is excellent,
in the matter of the derivation of the growth equation their presentation
is incomplete and points confusing because it misleads the reader. What
follows is an attempt to provide a more comprehensive derivation. The
presentation of the derivation differs from that of Rogers and Yau's in
notation. The generic radial distance variable is denoted by the lower
case symbol r and upper case R is used to denote the particular radius of
the droplet. This is the opposite of Rogers and Yau. The environmental
levels of water vapor density and temperature are identified by a
subscript e, in contrast to Rogers and Yau who denote the environmental
levels by the absence of a subscript. The values of water vapor
density and temperature are denoted here by a subscript R.
The Analysis
Cloud droplets may grow by condensation or diminish by evaporation. Both condensation and evaporation are diffusion processes described by the Fick diffusion equation; i.e., the flux of water vapor passing through a plane M is given by:
where D is the diffusion coefficient for water vapor and ρ is the density of water vapor in the space surrounding the droplet. The negative sign in the above equation reflects the fact that the mass flow is in the opposite direction from the gradient of water vapor density.
When water vapor diffuses into the droplet it brings heat energy which raises the temperature of the droplet and that heat energy may diffuse out of the droplet into the surrounding moist air. Such diffusion of heat energy obeys the same form of equation as does the diffusion of the water vapor. The heat energy flux Q passing through a plane is given by
When the net flow of water vapor or heat into or out of an infinitesimal volume is considered they must generate a change in the density and the temperature within that infinitesimal volume; i.e.,
Under steady-state conditions (∂ρ/∂t = 0 and ∂T/∂t = 0) the vapor density and temperature fields in the space surrounding the droplet must satisfy the equations:
In spherical coordinates under conditions of spherical symmetry the Laplacian ∇2 is:
where r is the distance from the center of the droplet.
For a droplet of radius R the boundary conditions are:
where ρe and Te are the environmental vapor density and temperature, respectively.
The solution to the steady-state diffusion equation reduces to integrating
once with respect to r to get
where C1 is a constant.
Dividing by r2 gives the equation
Integating this with respect to r gives
The integration constants must be chosen to satisfy the boundary conditions.
Subtracting the second equation from the first gives
The solution to the boundary value problem for the steady-state diffusion equation is therefore
This equation reflects the fact that if ρe is greater than ρR then ρ is decreasing smoothly as r decreases toward R.
Temperature obeys the same equations as water vapor density so
The latter form reflects the fact that if TR is greater than Te then T increases smoothly as r decreases to R.
At the surface of the droplet
Therefore the net amounts of water vapor mass and heat energy diffusing into and out of the droplet surface are:
The heat energy supplied to the droplet from the condensation of water vapor into the droplet is:
For equilibrium this must match the heat energy Q diffused out of the droplet; i.e.,
This above equation may be solved for the unknown ρR in terms of the unknown TR as
The water vapor density at the surface of the droplet, ρR, also obeys the ideal gas law so
where RV is the gas constant for water vapor and e's(TR) is the saturated vapor pressure at the temperature TR taking into account the Kelvin effect due to the curvature of the droplet.
From the Kelvin Equation and the Raoult effect e's is given by
where es is the vapor pressure of water over an infinite plane.
Combining the two equations for ρR gives
Dividing by (es/RVTR) gives:
The RHS of the above equation is a function of the known parameters and TR; the LHS is a cubic equation in (1/R). Thus for given ρe and Te the solution for R for any value of TR can be found. From the value of TR the value of ρR can be found. It is, of course, the inverse relationships of TR and ρR as functions of R that are desired but that is a mere mathematical transformation.
Let ρVS and eS be the density and vapor pressure of water vapor under saturated conditions. The Clasius-Clapyeron equation is usually expressed in the form:
where L is the latent heat of fusion and the α's refer to the specific densities of water vapor and liquid water under conditions of saturation. The specific volumes are just the reciprocals of the densities. Therefore the Clasius-Clapyeron equation can be expressed as:
Since ρVS in negligible compared to ρL the above equation reduces to:
But water vapor satisfies the ideal gas law
Equating the two expressions for deS/dT and dividing both sides by ρVSRVT gives:
The integration of this latter differential equation from Te to TR yields
At this point Mason makes use of an approximation:
Thus if (y/x) is approximately 1
Applying this approximation to the previous equation yields
Since (TR-Te) = L(dm/dt)/(4πKR) the previous equation can be expressed as
It was previously established that
The two equations above can be added together. The addition of the LHS gives:
Under the assumption that ρR=ρRS the above reduces to
where S is called the saturation ratio and (S-1) the supersaturation ratio.
The RHS of the two equations which were added together is
Solving for dm/dt gives
Since
It thus follows that
The density ρeS can be replaced by eeS/RVTe to obtain the growth equation:
the final modification is to assume TR to be equal to Te. Thus the final version of the growth equation is:
Sources: