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

Surendra Nepal,Magnus Ogren,Yosief Wondmagegne,Adrian Muntean

DOI: https://doi.org/10.1016/j.probengmech.2023.103546

2023-05-15

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.

Numerical Analysis,Probability

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the kinetic process of how diffusing substances (such as chlorides) penetrate into dense rubber in rubber materials. Specifically, the author focuses on the fact that in some materials science scenarios, the one - dimensional moving boundary problems (MBPs) combined with kinetic boundary conditions can reveal the large - time behavior of the penetration front of diffusing substances and directly estimate the service life of the material. The paper proposes a random walk algorithm that can provide a good numerical approximation of the concentration distribution and the position of the sharp front. This algorithm decouples the target evolution system into two steps: (i) Solve the ordinary differential equations related to the evaluation of the moving boundary velocity using the explicit Euler method; (ii) Solve the related diffusion problem using the random walk method. The correctness of the proposed random walk algorithm is verified by comparing with the results of the finite element method. The numerical experimental results are in good agreement with the penetration depth measured by the experimental device for dense rubber.

The core of the paper lies in solving a one - dimensional moving boundary problem describing the penetration of diffusing substances into dense rubber through the random walk method, especially dealing with problems under kinetic boundary conditions. This method can not only simulate the concentration distribution of diffusing substances, but also accurately track the dynamic changes of the moving boundary, thus providing a new tool for understanding the penetration characteristics of rubber materials and predicting their service life.