A Practical Walk-on-Boundary Method for Boundary Value Problems

05/08/2023
by   Ryusuke Sugimoto, et al.
0

We introduce the walk-on-boundary (WoB) method for solving boundary value problems to computer graphics. WoB is a grid-free Monte Carlo solver for certain classes of second order partial differential equations. A similar Monte Carlo solver, the walk-on-spheres (WoS) method, has been recently popularized in computer graphics due to its advantages over traditional spatial discretization-based alternatives. We show that WoB's intrinsic properties yield further advantages beyond those of WoS. Unlike WoS, WoB naturally supports various boundary conditions (Dirichlet, Neumann, Robin, and mixed) for both interior and exterior domains. WoB builds upon boundary integral formulations, and it is mathematically more similar to light transport simulation in rendering than the random walk formulation of WoS. This similarity between WoB and rendering allows us to implement WoB on top of Monte Carlo ray tracing, and to incorporate advanced rendering techniques (e.g., bidirectional estimators with multiple importance sampling, the virtual point lights method, and Markov chain Monte Carlo) into WoB. WoB does not suffer from the intrinsic bias of WoS near the boundary and can estimate solutions precisely on the boundary. Our numerical results highlight the advantages of WoB over WoS as an attractive alternative to solve boundary value problems based on Monte Carlo.

READ FULL TEXT

page 2

page 8

page 9

page 10

page 11

page 12

page 13

research
02/23/2023

Boundary Value Caching for Walk on Spheres

Grid-free Monte Carlo methods such as walk on spheres can be used to sol...
research
02/23/2023

Walk on Stars: A Grid-Free Monte Carlo Method for PDEs with Neumann Boundary Conditions

Grid-free Monte Carlo methods based on the walk on spheres (WoS) algorit...
research
08/03/2022

Solving Inverse PDE Problems using Grid-Free Monte Carlo Estimators

Modeling physical phenomena like heat transport and diffusion is crucial...
research
06/28/2020

Simplest random walk for approximating Robin boundary value problems and ergodic limits of reflected diffusions

A simple-to-implement weak-sense numerical method to approximate reflect...
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
01/31/2022

Grid-Free Monte Carlo for PDEs with Spatially Varying Coefficients

Partial differential equations (PDEs) with spatially-varying coefficient...
research
05/26/2023

Inferring the Future by Imagining the Past

A single panel of a comic book can say a lot: it shows not only where ch...

Please sign up or login with your details

Forgot password? Click here to reset