Solving Inverse PDE Problems using Grid-Free Monte Carlo Estimators

08/03/2022
by   Ekrem Fatih Yılmazer, et al.
0

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.

READ FULL TEXT

page 5

page 6

page 7

page 8

page 9

research
01/31/2022

Grid-Free Monte Carlo for PDEs with Spatially Varying Coefficients

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

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

We introduce the walk-on-boundary (WoB) method for solving boundary valu...
research
02/06/2023

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

In PDE-constrained optimization, one aims to find design parameters that...
research
09/16/2022

Solving Stochastic PDEs Using FEniCS and UQtk

The intrusive (sample-free) spectral stochastic finite element method (S...
research
09/25/2022

Monte Carlo Gradient in Optimization Constrained by Radiative Transport Equation

Can Monte Carlo (MC) solvers be directly used in gradient-based methods ...
research
07/19/2017

FDTD: solving 1+1D delay PDE

We present a proof of concept for adapting the finite-difference time-do...
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...

Please sign up or login with your details

Forgot password? Click here to reset