Grid-Free Monte Carlo for PDEs with Spatially Varying Coefficients

by   Rohan Sawhney, et al.

Partial differential equations (PDEs) with spatially-varying coefficients arise throughout science and engineering, modeling rich heterogeneous material behavior. Yet conventional PDE solvers struggle with the immense complexity found in nature, since they must first discretize the problem – leading to spatial aliasing, and global meshing/sampling that is costly and error-prone. We describe a method that approximates neither the domain geometry, the problem data, nor the solution space, providing the exact solution (in expectation) even for problems with extremely detailed geometry and intricate coefficients. Our main contribution is to extend the walk on spheres (WoS) algorithm from constant- to variable-coefficient problems, by drawing on techniques from volumetric rendering. In particular, an approach inspired by null-scattering yields unbiased Monte Carlo estimators for a large class of 2nd-order elliptic PDEs, which share many attractive features with Monte Carlo rendering: no meshing, trivial parallelism, and the ability to evaluate the solution at any point without solving a global system of equations.



page 1

page 2

page 11

page 12

page 13

page 14

page 15

page 16


Multilevel Monte Carlo estimators for elliptic PDEs with Lévy-type diffusion coefficient

General elliptic equations with spatially discontinuous diffusion coeffi...

KO-PDE: Kernel Optimized Discovery of Partial Differential Equations with Varying Coefficients

Partial differential equations (PDEs) fitting scientific data can repres...

Toeplitz Monte Carlo

Motivated mainly by applications to partial differential equations with ...

Strong L^p-error analysis of nonlinear Monte Carlo approximations for high-dimensional semilinear partial differential equations

Full-history recursive multilevel Picard (MLP) approximation schemes hav...

A note on a PDE approach to option pricing under xVA

In this paper we study partial differential equations (PDEs) that can be...

Robust Optimization of PDEs with Random Coefficients Using a Multilevel Monte Carlo Method

This paper addresses optimization problems constrained by partial differ...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.