Distribution of Colloidal Particles during Solvent Evaporation

Requires a Wolfram Notebook System

Interact on desktop, mobile and cloud with the free Wolfram CDF Player or other Wolfram Language products.

Requires a Wolfram Notebook System

Edit on desktop, mobile and cloud with any Wolfram Language product.

This Demonstration examines the underlying processes that take place when a colloidal dispersion of uniform thickness consolidates as solvent is removed. The overlying gas phase is assumed to be unsaturated with vapor, thereby providing the driving force for solvent evaporation to take place. The evaporation rate can be expressed as , where is an appropriate mass transfer coefficient and and are the vapor pressures at the interface and in the bulk air, respectively. As the solvent evaporates (), the film thins and particles begin to consolidate at the gas–liquid interface . The distribution of particles at any instant in time is denoted by .


The vertical distribution of particles in the uniformly consolidating film is controlled by the Péclet number , where is the Stokes–Einstein diffusion coefficient for the colloidal particles. For , diffusion of particles relative to convection of particles is strong, and an almost uniform particle distribution during drying is therefore expected (see Snapshot 1). On the other hand, for , we expect to see a skin-like formation with a random close packing of particles having a volume fraction in the neighborhood of the advancing interface (see Snapshot 2).

To solve this problem, it is convenient to transform the moving domain to a fixed domain by introducing the dimensionless variables and .

The species balance for volume fraction becomes:


subject to the following boundary conditions:

at ,

at .

We solve the partial differential equation using the spectral collocation method with the number of collocation points .

Use the sliders to vary the initial volume fraction , the Péclet number and the evaporation time .

The green curve shows the evolution of the interface volume fraction, in the drying film for . At , this value is equal to (shown by the dotted red horizontal line). At , this curve intersects the blue curve (i.e. the curve representing the particle density profile expressed in terms of the volume fraction or versus ) at a point marked in cyan. During drying (i.e. increased values of ), the volume fraction increases as the solvent evaporates and decreases. For and near the interface, approaches its limiting value, , as is evident from Snapshot 2.


Contributed by: Brian G. Higgins and Housam Binous (September 2016)
Open content licensed under CC BY-NC-SA




[1] A. F. Routh and W. B. Zimmerman, "Distribution of Particles during Solvent Evaporation from Films," Chemical Engineering Science, 59(14), 2004 pp. 2961–2968. doi:10.1016/j.ces.2004.04.027.

Feedback (field required)
Email (field required) Name
Occupation Organization
Note: Your message & contact information may be shared with the author of any specific Demonstration for which you give feedback.