Trapping of Planar Brownian Motion: Full First Passage Time Distributions by Kinetic Monte-Carlo, Asymptotic and Boundary Integral Methods

12/13/2021
by   Jake Cherry, et al.
0

We consider the problem of determining the arrival statistics of unbiased planar random walkers to complex target configurations. In contrast to problems posed in finite domains, simple moments of the distribution, such as the mean (MFPT) and variance, are not defined and it is necessary to obtain the full arrival statistics. We describe several methods to obtain these distributions and other associated quantities such as splitting probabilities. One approach combines a Laplace transform of the underlying parabolic equation with matched asymptotic analysis followed by numerical transform inversion. The second approach is similar, but uses a boundary integral equation method to solve for the Laplace transformed variable. To validate the results of this theory, and to obtain the arrival time statistics in very general configurations of absorbers, we introduce an efficient Kinetic Monte Carlo (KMC) method that describes trajectories as a combination of large but exactly solvable projection steps. The effectiveness of these methodologies is demonstrated on a variety of challenging examples highlighting the applicability of these methods to a variety of practical scenarios, such as source inference. A particularly useful finding arising from these results is that homogenization theories, in which complex configurations are replaced by equivalent simple ones, are remarkably effective at describing arrival time statistics.

READ FULL TEXT

page 1

page 11

page 12

page 19

page 20

page 21

page 23

research
10/10/2019

Time-Domain Boundary Integral Methods in Linear Thermoelasticity

This paper is concerned with the application of time domain boundary int...
research
11/20/2019

A Scrambled Method of Moments

Quasi-Monte Carlo (qMC) methods are a powerful alternative to classical ...
research
05/22/2022

On the inversion of the Laplace transform (In Memory of Dimitris Gatzouras)

The Laplace transform is a useful and powerful analytic tool with applic...
research
11/12/2021

Nyström methods for high-order CQ solutions of the wave equation in two dimensions

We investigate high-order Convolution Quadratures methods for the soluti...
research
07/29/2019

A Path Integral Monte Carlo Method based on Feynman-Kac Formula for Electrical Impedance Tomography

A path integral Monte Carlo method (PIMC) based on Feynman-Kac formula f...
research
07/22/2021

On Integral Theorems: Monte Carlo Estimators and Optimal Functions

We introduce a class of integral theorems based on cyclic functions and ...
research
07/09/2018

Computer Assisted Localization of a Heart Arrhythmia

We consider the problem of locating a point-source heart arrhythmia usin...

Please sign up or login with your details

Forgot password? Click here to reset