Solving partial differential equations in participating media

Bailey Miller	Rohan Sawhney	Keenan Crane	Ioannis Gkioulekas

ACM Trans. Graph. (2025)

We consider the problem of solving partial differential equations (PDEs) in domains with complex microparticle geometry that is impractical, or intractable, to model explicitly. Drawing inspiration from volume rendering, we propose tackling this problem by treating the domain as a participating medium that models microparticle geometry stochastically, through aggregate statistical properties (e.g., particle density). We first introduce the problem setting of PDE simulation in participating media. We then specialize to exponential media and describe the properties that make them an attractive model of microparticle geometry for PDE simulation problems. We use these properties to develop two new algorithms, volumetric walk on spheres and volumetric walk on stars, that generalize previous Monte Carlo algorithms to enable efficient and discretization-free simulation of linear elliptic PDEs (e.g., Laplace) in participating media. We demonstrate experimentally that our algorithms can solve Laplace boundary value problems with complex microparticle geometry more accurately and more efficiently than previous approaches, such as ensemble averaging and homogenization.

Bailey Miller, Rohan Sawhney, Keenan Crane, Ioannis Gkioulekas (2025). Solving partial differential equations in participating media. ACM Trans. Graph., 44(4).

@article{Miller:VWoS:2025,
author = {Bailey Miller and Rohan Sawhney and Keenan Crane and Ioannis Gkioulekas},
title = {Solving partial differential equations in participating media},
year = {2025},
issue_date = {July 2025},
publisher = {Association for Computing Machinery},
address = {New York, NY, USA},
volume = {44},
number = {4},
issn = {0730-0301},
url = {https://doi.org/10.1145/3731152},
doi = {10.1145/3731152},
journal = {ACM Trans. Graph.},
month = {jul},
numpages = {20},
}