Date: Friday, October 04, 2019
Location: 1084 East Hall (3:00 PM to 4:00 PM)
Title: Multiscale numerical modeling of sorption kinetics
Abstract: The trapping of diffusing particles by either a single or a distribution of moving traps is an interesting topic that has been employed to model a variety of different real problems in chemistry, physics and biology. Here we study the dynamics of diffusing particles in a domain with an oscillating bubble. Laboratory experiments provide evidence of a non monotone behavior in time of the concentration of particles by a detector located behind the bubble, under suitable experimental condition. A comprehensive explanation of the phenomenon is not yet fully available. The particles are attracted and trapped near the surface of the bubble. The basic mathematical model is a driftdiffusion model, where the particles diffuse and feel the potential of the bubble when they are near its surface. A tentative explanation of the mechanism is based on twocarrier dynamics.
The numerical simulation of the system presents two multiscale challenges. One is spatial: the range of the bubble potential is confined within a few microns at the bubble surface, while the bubble radius is of the order of a millimeter, so a fully resolved solution would be too expensive. The second challenge is on the time scale: the bubble oscillates with a frequency of the order of 100 Hz, while the diffusion time scale is of the order of 1000 seconds, this requiring at least one million time steps to fully resolve the problem in time.
A reduced model is derived to solve the multiscale problem it space for the single carrier dynamics: the interaction with the bubble is modeled as a very thin layer, with a particle surface density proportional to the local density in the bulk, near the bubble. In the rest of the domain the particle density satisfies just a diffusion equation, with suitable boundary conditions on the bubble, deduced from conservation properties.
The model is carefully tested on problems in 1D, 2D planar and 3D axissymmetric geometry. The equation is discretized on a regular Cartesian mesh, using a ghostpoint approach, and solved by CrankNicolson scheme. The implicit step is efficiently solved by a suitably adapted multigrid method. The amplitude of the bubble oscillations is small compared to the bubble radius. We take advantage of this fact by replacing the time dependent position by a suitable time dependent velocity at the bubble surface. Because of the low Reynolds number, the velocity distribution is computed by Stokes approximation. The multiscale challenge in time, as well as the multiscale model for multicarrier dynamics are still under investigation.
Files:
Speaker: Giovanni Russo
Institution: University of Catania and University of Michigan
Event Organizer: Robert Krasny krasny@umich.edu
