DeepAI AI Chat
Log In Sign Up

Solving Inverse PDE Problems using Grid-Free Monte Carlo Estimators

by   Ekrem Fatih Yılmazer, et al.

Modeling physical phenomena like heat transport and diffusion is crucially dependent on the numerical solution of partial differential equations (PDEs). A PDE solver finds the solution given coefficients and a boundary condition, whereas an inverse PDE solver goes the opposite way and reconstructs these inputs from an existing solution. In this article, we investigate techniques for solving inverse PDE problems using a gradient-based methodology. Conventional PDE solvers based on the finite element method require a domain meshing step that can be fragile and costly. Grid-free Monte Carlo methods instead stochastically sample paths using variations of the walk on spheres algorithm to construct an unbiased estimator of the solution. The uncanny similarity of these methods to physically-based rendering algorithms has been observed by several recent works. In the area of rendering, recent progress has led to the development of efficient unbiased derivative estimators. They solve an adjoint form of the problem and exploit arithmetic invertibility to compute gradients using a constant amount of memory and linear time complexity. Could these two lines of work be combined to compute cheap parametric derivatives of a grid-free PDE solver? We investigate this question and present preliminary results.


page 5

page 6

page 7

page 8

page 9


Grid-Free Monte Carlo for PDEs with Spatially Varying Coefficients

Partial differential equations (PDEs) with spatially-varying coefficient...

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

We introduce the walk-on-boundary (WoB) method for solving boundary valu...

Reversible random number generation for adjoint Monte Carlo simulation of the heat equation

In PDE-constrained optimization, one aims to find design parameters that...

A nonsmooth primal-dual method with simultaneous adaptive PDE constraint solver

We introduce an efficient first-order primal-dual method for the solutio...

FDTD: solving 1+1D delay PDE

We present a proof of concept for adapting the finite-difference time-do...

Monte Carlo Gradient in Optimization Constrained by Radiative Transport Equation

Can Monte Carlo (MC) solvers be directly used in gradient-based methods ...

FDTD: solving 1+1D delay PDE in parallel

We present a proof of concept for solving a 1+1D complex-valued, delay p...