Modeling desorption of fluids from disordered mesoporous materials

The desorption mechanism of fluids in disordered mesoporous glasses is studied by Monte Carlo simulations of a coarse-grained lattice model with realistic matrix configurations representative of Vycor. Two methods of simulation are considered: grand canonical ensemble Monte Carlo simulations and dynamic Monte Carlo simulations which mimic the diffusion of the fluid in and out of the material using Kawasaki dynamics. In the grand canonical simulations, cavitation via nucleation of bubbles inside the pores plays the dominant role in determining the fluid configurations along the desorption isotherm. The Kawasaki dynamics simulations indicate that such configurations are achieved dynamically via the gradual advancement of macroscopic front interfaces toward the interior. This is made possible by the bubble nucleation mechanism operating on a length scale that is determined by both the typical pore size and the strength of the solid-fluid interaction.