Global random walk for one-dimensional one-phase Stefan-type moving-boundary problems: Simulation results

Nicolae Suciu,Surendra Nepal,Yosief Wondmagegne,Magnus Ögren,Adrian Muntean

2024-10-16

Abstract:This work presents global random walk approximations of solutions to one-dimensional Stefan-type moving-boundary problems. We are particularly interested in the case when the moving boundary is driven by an explicit representation of its speed. This situation is usually referred to in the literature as moving-boundary problem with kinetic condition. As a direct application, we propose a numerical scheme to forecast the penetration of small diffusants into a rubber-based material. To check the quality of our results, we compare the numerical results obtained by global random walks either using the analytical solution to selected benchmark cases or relying on finite element approximations with a priori known convergence rates. It turns out that the global random walk concept can be used to produce good quality approximations of the weak solutions to the target class of problems.

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is how to use the Global Random Walk (GRW) method to approximately solve the weak solution of the one - phase Stefan - type moving boundary problem with kinetic conditions in the one - dimensional Stefan - type moving boundary problem. Specifically, the paper focuses on the situation when the velocity of the moving boundary is driven by the explicit representation of its velocity, which is usually referred to as the moving boundary problem with kinetic conditions in the literature. The paper proposes a numerical scheme for predicting the penetration of small diffusive substances into rubber - based materials, and checks the quality of the numerical results by comparing with the analytical solutions of selected benchmark cases or the finite - element approximation results depending on the a priori known convergence rate. The research shows that the global random walk concept can be used to generate high - quality approximations of the weak solutions of the target problem class.

In short, the core problem of the paper is to explore to what extent the global random walk method can approximately solve the weak solution of the one - phase Stefan - type moving boundary problem with kinetic conditions. This involves developing a new numerical method to deal with such problems more effectively, especially when it comes to specific scenarios in engineering applications, such as the penetration of small diffusive substances into hyperelastic materials, and this need is particularly important.