Yilin Ye,Adrien Chaigneau,Denis S. Grebenkov
Abstract:We propose an efficient numerical approach to simulate the boundary local time, as well as the time and position of the associated reaction event on a smooth boundary of a Euclidean domain. This approach combines the standard walk-on-spheres algorithm in the bulk with the approximate solution of the escape problem in the boundary layer. In this way, the most time-consuming simulation of reflected Brownian motion near the boundary is replaced by an equivalent escape event. We validate the proposed escape-from-a-layer approach by comparing simulated statistics of the boundary local time with exact results known for simple domains (a disk, a circular annulus, a sphere, a spherical shell) and with the numerical results obtained by a finite-element method in more sophisticated domains. This approach offers a powerful tool for simulating diffusive processes in confinements and for solving the related partial differential equations. Its applications in the context of diffusion-controlled reactions in chemical physics are discussed.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is how to efficiently simulate the boundary local time in the Euclidean domain, as well as the time and location of the related reaction events. Specifically, the author proposes a new numerical method - the "escape - from - a - layer" (EFL) method, to replace the detailed and time - consuming simulation of the reflected Brownian motion near the boundary in traditional methods. This method combines the standard walk - on - spheres (WOS) algorithm and the approximate solution to the escape problem within the boundary layer, thereby significantly improving the computational efficiency.
### Background and Problem Description
In the study of many physical, chemical, and biological processes, the boundary local time is an important concept, especially in cases involving diffusion - controlled reactions. The boundary local time \(\ell_t\) describes the residence time of particles near the boundary and is crucial for understanding reaction kinetics and particle behavior. However, for domains with complex geometric structures, directly simulating the boundary local time and its related statistics is very difficult because a large number of reflection events need to be processed, which leads to extremely high computational costs.
### Main Contributions of the Paper
1. **Proposing the EFL Method**: This method greatly reduces the computation time by introducing an escape event within the boundary layer to replace the detailed trajectory simulation. Specifically, when a particle enters the boundary layer, instead of simulating its detailed trajectory within the layer, the escape time, escape location, and the boundary local time obtained during this process are directly generated.
2. **Theoretical Basis**: The author uses the encounter - based approach to derive the spectral expansion form of the joint probability density of the escape time, escape location, and boundary local time, and provides the exact expressions of these quantities in the Laplace domain.
3. **Numerical Verification**: The effectiveness and accuracy of the EFL method are verified by comparing with the exact results in simple domains (such as discs, annular regions, spheres, spherical shells) and the results of the finite element method in complex domains.
4. **Application Prospects**: This method is not only applicable to simulating diffusion processes but can also be used to solve related partial differential equations, and is especially suitable for the study of diffusion - controlled reactions in chemical physics.
### Method Overview
- **Walk - on - Spheres Algorithm (WOS)**: Use the WOS algorithm in the area far from the boundary to quickly simulate large - scale displacements.
- **Layer Escape Problem**: When a particle enters the boundary layer, replace the detailed simulation by generating an escape event. The escape event includes the escape time \(\tau\), the escape location \(X_\tau\), and the boundary local time \(\ell_\tau\).
- **Approximation Methods**:
- **Flat Layer Approximation (FLA)**: Assume that the boundary layer is flat to simplify the calculation.
- **Curved Layer Approximation (CLA)**: Consider the curvature of the boundary layer to improve the accuracy.
### Conclusion
The EFL method proposed in this paper provides a new tool for efficiently simulating the boundary local time, especially when dealing with domains with complex geometric structures. This method not only improves the computational efficiency but also provides a new perspective for the study of diffusion - controlled reactions.