Walk on Heat Stars provides a boundary-integral Monte Carlo solver for parabolic PDEs with Neumann conditions via exact heat-ball sampling that yields unbiased estimators.
Monte Carlo PDE Solvers for Nonlinear Radiative Boundary Conditions
1 Pith paper cite this work. Polarity classification is still indexing.
abstract
Monte Carlo PDE solvers have become increasingly popular for solving heat-related partial differential equations in geometry processing and computer graphics due to their robustness in handling complex geometries. While existing methods can handle Dirichlet, Neumann, and linear Robin boundary conditions, nonlinear boundary conditions arising from thermal radiation remain largely unexplored. In this paper, we introduce a Picard-style fixed-point iteration framework that enables Monte Carlo PDE solvers to handle nonlinear radiative boundary conditions. While strict theoretical convergence is not generally guaranteed, our method remains stable and empirically convergent with a properly chosen relaxation coefficient. Even with imprecise initial boundary estimates, it progressively approaches the correct solution. Compared to standard linearization strategies, the proposed approach achieves significantly higher accuracy. To further address the high variance inherent in Monte Carlo estimators, we propose a heteroscedastic regression-based denoising technique specifically designed for on-boundary solution estimates, filling a gap left by prior variance reduction methods that focus solely on interior points. We validate our approach through extensive evaluations on synthetic benchmarks and demonstrate its effectiveness on practical heat radiation simulations with complex geometries.
fields
cs.GR 1years
2026 1verdicts
UNVERDICTED 1representative citing papers
citing papers explorer
-
Walking on Heat Stars for Parabolic Heat Equations with Neumann Boundary Conditions
Walk on Heat Stars provides a boundary-integral Monte Carlo solver for parabolic PDEs with Neumann conditions via exact heat-ball sampling that yields unbiased estimators.