For dense particle aggregates, pore network models are used to describe the void space by a network of geometrically similar pores, and discrete transport rules are set up to describe drying as a sequence of pore emptying. Thus the evolution of liquid distributions can be simulated. In a one-way coupling scheme (liquid - solid), capillary forces are computed over time from the filling state of pores and applied as loads on each particle in discrete element method. If bond strength between particles is exceeded, individual contacts will break up and particles may reorganize, leading to more or less sever structural damage. The influence of liquid phase distributions during drying and material properties on mechanical response (shrinkage and crack) is presented.
In case of highly porous particle aggregates, instead of approximating the void space by a pore network, a discretization of the full domain of interest - consisting of solid, liquid, and gas - is employed: The phase distributions are described by time-dependent cell volume fractions on a stationary cubic mesh. The solid phase volume fractions are computed from an arbitrary collection of spherical primary particles. The volume of fluid method is used to track the liquid-gas interface over time. Local evaporation rates are computed from a finite difference solution of a vapor diffusion problem in the gas phase, and the liquid-gas interface dynamics is described by volume-conserving mean curvature flow, with an additional equilibrium contact angle condition along the three-phase contact lines. Numerical results illustrate the capillary effects commonly observed in experiments. The evolutions of the liquid distribution over time for different wetting properties of the solid surface as well as binary liquid bridges between solid particles are presented.