Random walks and moving boundaries: Estimating the penetration of diffusants into dense rubbers

Abstract

For certain materials science scenarios arising in rubber technology, one-dimensional moving boundary problems (MBPs) with kinetic boundary conditions are capable of unveiling the large-time behavior of the diffusants penetration front, giving a direct estimate on the service life of the material. In this paper, we propose a random walk algorithm able to lead to good numerical approximations of both the concentration profile and the location of the sharp front. Essentially, the proposed scheme decouples the target evolution system in two steps: (i) the ordinary differential equation corresponding to the evaluation of the speed of the moving boundary is solved via an explicit Euler method, and (ii) the associated diffusion problem is solved by a random walk method. To verify the correctness of our random walk algorithm we compare the resulting approximations to results based on a finite element approach with a controlled convergence rate. Our numerical experiments recover well penetration depth measurements of an experimental setup targeting dense rubbers.

arXiv comment15 pages, 10 figures, 2 tables

References in corpus (8)

• Realistic boundary conditions for stochastic simulations of reaction-diffusion processes
• A second order scheme for a Robin boundary condition in random walk algorithms
• Global random walk solvers for fully coupled flow and transport in saturated/unsaturated porous media (extended version)
• Local boundary conditions for NMR-relaxation in digitized porous media
• A moving boundary approach of capturing diffusants penetration into rubber: FEM approximation and comparison with laboratory measurements
• Analysis of a fully discrete approximation to a moving-boundary problem describing rubber exposed to diffusants
• Stochastic solutions of Stefan problems
• Numerical simulations of NMR relaxation in chalk using local Robin boundary conditions