{"title":"Monte Carlo physics-informed neural networks for multiscale heat conduction via phonon Boltzmann transport equation","authors":"Qingyi Lin , Chuang Zhang , Xuhui Meng , Zhaoli Guo","doi":"10.1016/j.jcp.2025.114364","DOIUrl":null,"url":null,"abstract":"<div><div>The phonon Boltzmann transport equation (BTE) is widely used for the description of multiscale heat conduction (from nm to <span><math><mi>μ</mi></math></span>m or mm) in solid materials. Developing numerical approaches to solve this equation is challenging since it is a 7-dimensional integral-differential equation. In this work, we propose the Monte Carlo physics-informed neural networks (MC-PINNs), which provide an effective way to combat the <em>“curse of dimensionality”</em> in solving the phonon Boltzmann transport equation for modeling multiscale heat conduction in solid materials. In MC-PINNs, we utilize a deep neural network to approximate the solution to the BTE and encode the BTE as well as the corresponding boundary/initial conditions using automatic differentiation. In addition, we propose a novel two-step sampling approach to address the issues of inefficiency and inaccuracy in the widely used sampling methods in PINNs. In particular, we first randomly sample a certain number of points in the temporal-spatial space (Step I) and then draw another number of points randomly in the solid angular space (Step II). The training points at each step are constructed based on the data drawn from the above two steps using the tensor product. The two-step sampling strategy enables the MC-PINNs (1) to model the heat conduction from ballistic to diffusive regimes, and (2) to be more memory efficient compared to the conventional numerical solvers or existing PINNs for BTE. A series of numerical examples, including quasi-one-dimensional (quasi-1D) steady/unsteady heat conduction in a film, and the heat conduction in quasi-two-dimensional (quasi-2D) and three-dimensional (3D) domains, are conducted to justify the effectiveness of the MC-PINNs for heat conduction spanning diffusive and ballistic regimes. Finally, we perform a comparison on the computational time and the memory usage between the MC-PINNs and one of the state-of-the-art numerical methods to demonstrate the potential of MC-PINNs for large-scale problems in real-world applications.</div></div>","PeriodicalId":352,"journal":{"name":"Journal of Computational Physics","volume":"542 ","pages":"Article 114364"},"PeriodicalIF":3.8000,"publicationDate":"2025-09-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Journal of Computational Physics","FirstCategoryId":"101","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0021999125006461","RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q2","JCRName":"COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS","Score":null,"Total":0}
引用次数: 0
Abstract
The phonon Boltzmann transport equation (BTE) is widely used for the description of multiscale heat conduction (from nm to m or mm) in solid materials. Developing numerical approaches to solve this equation is challenging since it is a 7-dimensional integral-differential equation. In this work, we propose the Monte Carlo physics-informed neural networks (MC-PINNs), which provide an effective way to combat the “curse of dimensionality” in solving the phonon Boltzmann transport equation for modeling multiscale heat conduction in solid materials. In MC-PINNs, we utilize a deep neural network to approximate the solution to the BTE and encode the BTE as well as the corresponding boundary/initial conditions using automatic differentiation. In addition, we propose a novel two-step sampling approach to address the issues of inefficiency and inaccuracy in the widely used sampling methods in PINNs. In particular, we first randomly sample a certain number of points in the temporal-spatial space (Step I) and then draw another number of points randomly in the solid angular space (Step II). The training points at each step are constructed based on the data drawn from the above two steps using the tensor product. The two-step sampling strategy enables the MC-PINNs (1) to model the heat conduction from ballistic to diffusive regimes, and (2) to be more memory efficient compared to the conventional numerical solvers or existing PINNs for BTE. A series of numerical examples, including quasi-one-dimensional (quasi-1D) steady/unsteady heat conduction in a film, and the heat conduction in quasi-two-dimensional (quasi-2D) and three-dimensional (3D) domains, are conducted to justify the effectiveness of the MC-PINNs for heat conduction spanning diffusive and ballistic regimes. Finally, we perform a comparison on the computational time and the memory usage between the MC-PINNs and one of the state-of-the-art numerical methods to demonstrate the potential of MC-PINNs for large-scale problems in real-world applications.
期刊介绍:
Journal of Computational Physics thoroughly treats the computational aspects of physical problems, presenting techniques for the numerical solution of mathematical equations arising in all areas of physics. The journal seeks to emphasize methods that cross disciplinary boundaries.
The Journal of Computational Physics also publishes short notes of 4 pages or less (including figures, tables, and references but excluding title pages). Letters to the Editor commenting on articles already published in this Journal will also be considered. Neither notes nor letters should have an abstract.