Dynamics of variations in size and composition of supercritical droplet at binary condensation

A general analytical solution is found in quadratures for the radius and concentration of a solution droplet, which isothermally grows or evaporates in a diffusion or free-molecular regime in a binary mixture of vapors. The obtained solution describes the dynamics of variations in the size and composition of a super-critical droplet during the binary condensation in mixed vapors at an arbitrary initial droplet composition. It is shown that, at small (linear) deviations of the growth regime and droplet composition from the stationarity, these quadratures lead to the results that were recently obtained for the composition relaxation in a growing droplet. Moreover, it is demonstrated that, in terms of the nonlinear theory, when the deviation of solution concentration in a droplet from its stationary value is not small, it is invalid to use the law of stationary variations in the size of a droplet with time to describe the relaxation process for its chemical composition.     

Stationary concentration establishment in the growing or evaporating droplet of ideal binary solution

The establishment of stationary solution concentration in a growing or evaporating droplet of an ideal binary solution (binary droplet) placed in a vapor mixture of constituting substances and passive gas is described analytically. Relations defining time dependences of solution concentration in a droplet, the number of molecules of each constituting component, and droplet radius are derived at known parameters of the vapor-gas mixture and the initial composition of a binary droplet. The results of calculations of time dependences of aforementioned values are reported for several variants of the initial composition of a droplet and a vapor mixture.     

Thermal effects accompanying stationary binary condensation of vapors into overcritical droplet

A set of equations is derived to calculate the stationary temperature and concentration of a solution in a overcritical droplet with regard to the heat release accompanying the condensation of a binary mixture of vapors in a diffusion or free-molecular regime. In the approximation of an ideal solution, relations are found for the stationary temperature of droplet growing under the conditions of strong and weak thermal effects. For the general case and the cases of strong and weak thermal effects, the temperature and concentration of the droplet and the coefficient of the thermal deceleration of the droplet growth are calculated as functions of the density of a passive gas. The influence of the condensation heat values of the first and second components of the mixture on the stationary temperature and concentration of the solution in the growing droplet is investigated separately.     

Regularities of binary condensation of vapors when one of them is undersaturated

Simple analytical expressions are derived for the stationary concentration of a binary solution in a markedly supercritical droplet growing exothermically in diffusion or free-molecular regimes in mixed vapors when one of condensing vapors is supersaturated and present in a small amount and another vapor is slightly undersaturated and present in a large amount. The condensation of sulfuric acid and water vapors on a droplet under the conditions of Earth atmosphere is considered as an example of practical importance. Under isothermic conditions, analytical expressions are obtained for the time of establishing a power law for variations in a droplet radius with time under the diffusion and free-molecular regimes of the droplet growth. The power laws are derived in an explicit form under these regimes, which describe the rapid establishment of a stationary concentration of a solution in a growing droplet.     

Evaporation dynamics of a binary sessile droplet: Theory and comparison with experimental data on a droplet of a sulfuric-acid solution

Diffusion evaporation of a sessile binary droplet in an atmosphere of a noncondensable carrier gas has been considered. For a droplet consisting of two infinitely miscible liquids, a relation between the current values of solution concentration and volume of the droplet has been derived in an explicit form under the ideal solution approximation. It has been shown that the volume of a sessile binary droplet may, as well as the volume of a free binary droplet, vary nonmonotonically with time. The evaporation of a droplet of an aqueous sulfuric-acid solution has been considered in detail taking into account the nonideality of the solution. Time variations in the volume, base area, and contact angle have been experimentally measured for the sessile droplet of an aqueous sulfuric-acid solution on a hydrophobized substrate. The experimental data obtained at different initial humidities of water-vapor and droplet-solution concentrations have been analyzed within the theory of the stationary isothermal diffusion evaporation of a sessile binary droplet.     

Nonstationary vapor concentration fields near the growing droplet of binary solution

Nonstationary vapor concentration fields near the droplet of binary solution growing in the vapor—gas mixture are revealed using the concepts of similarity. The revealed fields are determined with the exact account of the motion of droplet surface and refer to the times at which the droplet reaches sizes that provide for the diffusion regime of droplet growth. To obtain the self-similar solution of the problem of binary condensation, it is necessary to ensure a constant (in time) concentration of binary solution in the growing droplet. The velocities of an increase in the number of molecules and the radius of two-component droplet with time are found with allowance for the equation ensuring this solution. The conditions for the transformation of the self-similar solution of the problem of the condensation of two-component mixture into the solution, which was derived previously for the condensation of one component, are elucidated.     

Dynamics of variations in size and composition of a binary droplet in a mixture of two condensing vapors and a passive gas under arbitrary initial conditions

The dynamics of variations in the size and composition of a droplet in a mixture of two vapors and a noncondensable carrier gas has been studied at an arbitrary initial droplet size and two limiting initial concentrations of the solution in the droplet corresponding to a pure first or pure second component. The conditions for nonmonotonic variations in the droplet radius with time have been analyzed. The physical situation has been investigated for the course of binary condensation, in which, at an initial stage, before the droplet begins to grow and the stationary concentration of the solution in it is established, the droplet size markedly decreases. The opposite situation is also considered, in which the droplet grows at the initial stage and then passes to the regime of monotonic evaporation.     

Diffusion of Vapor in the Presence of a Growing Droplet

As is shown, the solution to the diffusion equation for the concentration of vapor in the presence of a droplet growing in it, derived for the usual initial condition and equilibrium boundary conditions at the droplet surface, fails to ensure an equality between the numbers of molecules that have left the vapor due to diffusion by the current moment and those that have been included in the growing droplet. The difference between the total numbers of vapor molecules at the initial moment (when the vapor had a given uniform concentration) and at the current moment (when the size of the growing droplet is much larger than its initial size) differs from the total number of molecules in the droplet by a factor of 3/2. By substituting the usual boundary condition at the droplet surface by a time-dependent boundary condition at the surface of a constant-radius sphere with the center in the center of the growing droplet, a solution to the diffusion problem for the vapor concentration is derived. This solution describes the evolution of the vapor concentration field, which agrees with the rate of the vapor absorption by the growing droplet and with the law of the conservation of matter.     

Self-similar theory of nonisothermal vapor condensation on a droplet growing in a vapor-gas medium

Rigorous self-similar solutions to the joint problems of vapor diffusion toward a droplet growing in a vapor-gas medium and the removal of heat released during vapor condensation are found. An equation for the temperature of a droplet ensuring the existence of a self-similar solution is derived. This equation sets the constancy of the temperature of a droplet throughout the time of its growth and unambiguously determines this temperature. In the case of the strong heat effects, when the rate of droplet growth decreases substantially, the analytical solution to this equation is obtained. This temperature coincides precisely with the temperature, which is established in the droplet at the diffusion regime of its growth. At the found droplet temperature, interconnected fields of vapor concentration and temperature of vapor-gas medium around the droplet are expressed through the initial (prior to the droplet nucleation) parameters of a vapor-gas medium. These parameters are used to express the dependence of the radius of a droplet on the time at the diffusion regime of its growth and the time required to establish the diffusion regime of droplet growth. The case of weak heat effects is also studied.     

Self-propelled motion of a fluid droplet under chemical reaction

We study self-propelled dynamics of a droplet due to a Marangoni effect and chemical reactions in a binary fluid with a dilute third component of chemical product which affects the interfacial energy of a droplet. The equation for the migration velocity of the center of mass of a droplet is derived in the limit of an infinitesimally thin interface. We found that there is a bifurcation from a motionless state to a propagating state of droplet by changing the strength of the Marangoni effect.     

Power-law stage of slow relaxation in solutions with spherical micelles

General (independent of models selected for surfactant molecular aggregates) analytical relations are derived to describe the initial stage of slow relaxation in micellar solutions with spherical micelles. This stage precedes the final stage of the relaxation occurring via an exponential decay of disturbances with time. The relations obtained are applicable throughout the interval of micellar solution concentrations from the first to the second critical micellization concentration. It is shown that the initial stage is characterized by power laws of variations in the concentrations of monomers and micelles with time, these laws being different for the relaxation processes proceeding from above and below toward equilibrium values of micellar solution parameters. Relations are derived for the duration of this stage, and the effect of initial conditions is studied. Characteristic times of the power-law stage are determined and compared with the characteristic time of the final exponent-law relaxation stage. The behavior of these times is investigated at surfactant solution concentrations in the vicinity of, and noticeably above, the first critical micellization concentration. On the basis of the droplet and quasi-droplet thermodynamic models of surfactant molecular aggregates, numerical solutions are found for nonlinearized equations of slow relaxation for the time dependence of surfactant monomer concentrations at all stages of the slow relaxation. Numerical results obtained from the models are compared with the results of a general analytical study.     

Kinetics of the Establishment of Quasi-Steady-State Regime of Overcoming Activation Barrier of Nucleation on the Macroscopic Wettable Nuclei

The solution of the kinetic equation of nucleation on macroscopic wettable condensation nuclei was constructed for the initial (incubation) stage. The solution thus constructed determines the times of relaxation to quasi-steady-state distribution of droplets generating on droplet nuclei in the vicinity of maximum of the work of droplet heterogeneous formation as well as the relaxation to quasi-equilibrium droplet distribution throughout the entire region located to the left of this vicinity at the droplet size axis. The dependence of relaxation times on the height of activation barrier of nucleation, size of nuclei, their nature, and characteristics of matter comprising condensate was elucidated. It was shown when the non-steady-state rate of nucleation becomes actually equal to the quasi-steady-state rate of nucleation.     

On the evolution of a multicomponent droplet during nonisothermal diffusion growth or evaporation

A set of equations has been derived for the size, composition, and temperature of a multicomponent droplet of a nonideal solution during its diffusion nonisothermal condensation growth or evaporation in a multicomponent mixture of vapors with an incondensable carrier gas. In addition to complete equations for material and heat transfer in the vapor-gas medium surrounding the droplet, the derived set, in the general case, describes the nonstationary growth or evaporation of the droplet under arbitrary initial conditions (initial size and temperature of the droplet and the concentrations of the nonideal multicomponent solution in it) and the establishment of the stationary values of the composition, temperature, and the rate of variations in the size of the droplet with allowance for heat effects and diffusion and thermodiffusion material transfer, Stefan flux, motion of the droplet surface, and the nonideality of the solution in the droplet. A simplified set of equations obtained without taking into account the contributions from the flow, cross effects, and thermal expansion in the equations of the material and heat transfer in the vapor-gas medium has been considered. Equations describing growth/evaporation in the stationary regime have been analyzed for droplets of ideal multicomponent solutions.     

Analysis of the effects of Marangoni stresses on the microflow in an evaporating sessile droplet

We study the effects of Marangoni stresses on the flow in an evaporating sessile droplet, by extending a lubrication analysis and a finite element solution of the flow field in a drying droplet, developed earlier. The temperature distribution within the droplet is obtained from a solution of Laplace's equation, where quasi-steadiness and neglect of convection terms in the heat equation can be justified for small, slowly evaporating droplets. The evaporation flux and temperature profiles along the droplet surface are approximated by simple analytical forms and used as boundary conditions to obtain an axisymmetric analytical flow field from the lubrication theory for relatively flat droplets. A finite element algorithm is also developed to solve simultaneously the vapor concentration, and the thermal and flow fields in the droplet, which shows that the lubrication solution with the Marangoni stress is accurate for contact angles as high as 40 degrees. From our analysis, we find that surfactant contamination, at a surface concentration as small as 300 molecules/microm(2), can almost entirely suppress the Marangoni flow in the evaporating droplet.     

Micellization kinetics with allowance for fussion and fission of spherical and cylindrical micelles: 1. Set of nonlinear equations describing slow relaxation

Based on the general kinetic equation that describes the aggregation and fragmentation of surfactant molecular aggregates, a closed set of nonlinear equations is derived for the slow relaxation of surfactant monomer concentration and the total concentrations of coexisting spherical and cylindrical micelles to the equilibrium state of a micellar solution. Both the transitions accompanied by the emission and capture of surfactant monomers by micelles and the transitions resulting from the fussion and fission of micelles, are taken into account. The derived set of equations describes all stages of the slow relaxation from the initial perturbance to the final equilibrium state of a micellar solution.     

Mass dissolution of polydisperse set of particles in a flow system

A mathematical model of continuous dissolution of salts based on a dynamic equation for crystal size distribution function was considered. A stationary solution was obtained by the moments method. For a nonstationary state, an evolution equation for the undersaturation of solution was derived. An approximate analytical method of solving this equation, which is in good agreement with the numerical solution, was proposed. The notion of the coefficient of continuous dissolution efficiency was introduced and analytical expressions for this value were obtained. The stability of the stationary solution was analyzed.Translated from Kolloidnyi Zhurnal, Vol. 66, No. 6, 2004, pp. 793–801.Original Russian Text Copyright © 2004 by Moshinskii.     

Boltzmann distributions and slow relaxation in systems with spherical and cylindrical micelles

Equilibrium and nonequilibrium distributions of molecular aggregates in a solution of a nonionic surfactant are investigated at the total surfactant concentration above the second critical micelle concentration (CMC2). The investigation is not limited by the choice of a specific micellar model. Expressions for the direct and reverse fluxes of molecular aggregates over the potential humps of the aggregation work are derived. These aggregation work humps set up activation barriers for the formation of spherical and cylindrical micelles. With the aid of the expressions for molecular aggregate fluxes, a set of two kinetic equations of micellization is derived. This set, along with the material balance equation, describes the molecular mechanism of the slow relaxation of micellar solution above the CMC2. A realistic situation has been analyzed when the CMC2 exceeds the first critical micelle concentration, CMC1, by an order of magnitude, and the total surfactant concentration varies within the range lying markedly above the CMC2 but not by more than 2 orders of magnitude. For such conditions, an equation relating the parameters of the aggregation work of a cylindrical micelle to the observable ratio of the total surfactant concentration and the monomer concentration is found for an equilibrium solution. For the same conditions, but in the nonequilibrium state of materially isolated surfactant solution, a closed set of linearized relaxation equations for total concentrations of spherical and cylindrical micelles is derived. These equations determine the time development of two modes of slow relaxation in micellar solutions markedly above the CMC2. Solving the set of equations yields two rates and two times of slow relaxation.     

Self-similar solution to the problem of vapor diffusion toward the droplet nucleated and growing in a vapor-gas medium

The problem of vapor diffusion toward a droplet nucleated and growing in the diffusion regime is exactly solved using the similarity theory. The surface motion of droplets is taken into account in the solution. The constructed nonstationary concentration field of vapor satisfies the diffusion equation, the boundary condition of equilibrium on the surface of growing droplet, and the initial homogeneous condition. According to the found solution, the radius of a droplet is proportional to the square root of the time of its growth. Far from the critical point, at a low ratio between the densities of excess vapor and a liquid droplet, the proportionality coefficient coincides with that resulting from an approximate solution. The balance between the numbers of molecules removed from vapor and those composing a growing droplet exactly corresponds to the obtained solution.     

Postnucleation droplet growth in supersaturated gas with arbitrary vapor concentration

This work concerns the reexamination and extension of the current theory of phase transition dynamics for liquid droplets growing on soluble aerosols from a supersaturated gas mixture for the general case of arbitrary value of vapor concentration. We found that the inconsistency in the common treatment of the vapor diffusion, due to an implicit assumption of the constancy of gas density in the vicinity of a droplet by neglecting its dependency on temperature and vapor concentration, leads to the obvious discrepancy in the Maxwell expression for the growth rate regarding droplets of near critical size. Restoring the correct treatment of the vapor diffusion in terms of the mass concentration of water vapor and taking into the consideration variations of gas density in the vicinity of a droplet in compliance with the equation of state of moist air, we have obtained a new expression for the droplet growth rate valid for an arbitrary value of vapor concentration. The limitations imposed by the molecular kinetic fluxes to postnucleation diffusional growth of small droplets with a large Knudsen number are also reevaluated to include previously neglected physical effects. In particular, the essential contribution of the vapor molecular energy flux into the total kinetic molecular heat flux as well as the temperature variations of mean thermal velocities of air and vapor molecules in the vicinity of the droplet interface have been taken into consideration. Surprisingly significant differences have been found in new expressions derived for the droplet growth rate and droplet temperature, even in the limit of small vapor concentration, if comparing with commonly used results. These findings could help with better interpretation of experimental measurements to infer more reliable data for the mass and thermal accommodations coefficients.     

Nonlinear chromatography Recent theoretical and experimental results

The theory of nonlinear chromatography has been advanced by the incorporation of recent results obtained by the theory of partial differential equations. The system of equations of the ideal model has been solved analytically in the case of a single component for which the equilibrium isotherm between the mobile and the stationary phases is given by a Langmuir equation. A series of computer programs has been written which permits the calculation of numerical solutions of the semi-ideal model. The properties of the solutions obtained are described and discussed for a one-component system (profile of high concentration bands of a pure compound eluted by a pure solvent), several two-component systems (elution of a pure compound band by a binary mobile phase, separation of a binary mixture eluted by a pure mobile phase), and three-component systems (separation of a binary mixture eluted by a binary solvent, displacement and separation of a binary mixture). Experimental results are reported which validate the conclusions derived from the numerical integration of the model. The conclusions of the work apply to all high-performance chromatographic procedures, i.e., to those where the kinetics of mass transfer are fast enough for the mobile and stationary phases always to be near equilibrium. More specifically, the contribution from the kinetics of the retention mechanism to the mass transfer resistance must itself be negligible. This clearly excludes affinity chromatography.