Solving partial differential equations in participating media

2025-06-09

AI Summary

Conventional PDE solvers struggle with domains containing complex microstructural geometries—difficult to model explicitly—facing dual bottlenecks of inaccurate explicit modeling and insufficient homogenization accuracy. Method: We propose a mesh-free stochastic walk framework: treating the medium as a participating medium, its microstructure is implicitly encoded via statistical attributes (e.g., particle density); we introduce voxelized spherical and voxelized star-shaped random walk algorithms for efficient, high-accuracy solution of linear elliptic PDEs (e.g., Laplace boundary-value problems). The method integrates volume-rendering-inspired medium modeling, exponential-medium theoretical analysis, and Monte Carlo sampling—eliminating mesh discretization. Contribution/Results: Experiments demonstrate significant improvements in both accuracy and computational efficiency over ensemble averaging and homogenization methods, while exhibiting strong robustness to variations in particle geometry.

Technology Category

• Reasoning under Uncertainty: Stochastic Optimization
• Search and Optimization: Mixed Discrete/Continuous Search
• Machine Learning: Learning with Manifolds

Application Category

• Graph Algorithms and Modeling for the Web: Efficient manipulation of static and dynamic Web-related graphs
• Web Mining and Content Analysis: Robustness and generalizability of Web mining methods
• Systems and Infrastructure for Web, Mobile and WoT: Experiences and lessons learnt from Web-based algorithms and system deployments

Abstract

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.

Problem

• Solving PDEs in complex microparticle geometry domains
• Modeling domains as stochastic participating media
• Developing efficient Monte Carlo algorithms for PDEs

Innovation

• Treats domain as stochastic participating media
• Introduces volumetric walk on spheres
• Generalizes Monte Carlo for elliptic PDEs