Monte Carlo PDE Solvers for Nonlinear Radiative Boundary Conditions

2026-04-22

AI Summary

This work addresses the challenge of efficiently handling nonlinear radiative boundary conditions in Monte Carlo solvers for partial differential equations. The authors propose a novel framework that integrates Picard-type fixed-point iteration with heteroscedastic regression-based denoising. Specifically, they introduce fixed-point iteration into Monte Carlo PDE solvers for the first time to manage nonlinear radiation boundaries and develop a boundary-aware heteroscedastic regression technique to reduce estimation variance, thereby filling a critical gap in existing variance reduction methods for this setting. Experimental results demonstrate that the proposed approach achieves strong stability and empirical convergence on both synthetic benchmarks and realistic thermal radiation simulations involving complex geometries, significantly outperforming conventional linearization strategies in accuracy.

Technology Category

• Machine Learning: Calibration & Uncertainty Quantification
• Reasoning under Uncertainty: Stochastic Optimization
• Search and Optimization: Sampling/Simulation-based Search

Application Category

• Graph Algorithms and Modeling for the Web: Algorithms and analysis for incomplete, noisy, or partially observed Web-related graphs
• Search and Retrieval-Augmented AI: Web learning to rank, online learning, and counterfactual learning for ranking
• User Modeling, Personalization and Recommendation: User modeling and simulation for interactive and conversational systems

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.

Problem

• Monte Carlo PDE solvers
• nonlinear boundary conditions
• thermal radiation
• radiative boundary conditions

Innovation

• Monte Carlo PDE solvers
• nonlinear radiative boundary conditions
• Picard iteration
• heteroscedastic regression denoising
• variance reduction