Xiaojun Lei , Tongxiang Gu , Xiaowen Xu , Hengbin An , Yanzhong Yao
{"title":"2D 3T PINNs for solving 2-D 3-T heat conduction equations based on physics-informed neural networks","authors":"Xiaojun Lei , Tongxiang Gu , Xiaowen Xu , Hengbin An , Yanzhong Yao","doi":"10.1016/j.cpc.2025.109572","DOIUrl":null,"url":null,"abstract":"<div><div>The 2-D 3-T heat conduction equations can approximate the propagation of energy in the material, as well as the energy exchange process of electrons, ions and photons. In many cases, the computational time to solve these equations accounts for a large proportion (more than 80%) of that of the entire simulation of radiation hydrodynamics. PINNs is a promising way to solve partial differential equations (PDEs). Although numerical methods have been successful in solving 2-D 3-T heat conduction equations, the PINNs method also has some advantages, such as mesh-free, suitable to high dimension and complex domain problems. But the original PINNs cannot solve the 2-D 3-T heat conduction equations to a reasonable precision. This work aims to explore which techniques need to be added to PINNs and to what extent it can address the challenges posed by strong nonlinearity and multi-scale phenomena present in 2-D 3-T heat conduction equations. Hence, we adopt guaranteed positive constraint to the outputs so that the network can be trained, give a relatively large weight to the initial loss and Dirichlet boundary loss, take the logarithm of the initial loss, use transfer learning and Fourier feature embedding to improve accuracy. We call our improved approach 2D 3T PINNs. Numerical experiments show that the relative <span><math><msub><mrow><mi>L</mi></mrow><mrow><mn>2</mn></mrow></msub></math></span> and the absolute <span><math><msub><mrow><mi>L</mi></mrow><mrow><mo>∞</mo></mrow></msub></math></span> error between the 2D 3T PINNs prediction and the reference solution is of the order of <span><math><msup><mrow><mn>10</mn></mrow><mrow><mo>−</mo><mn>2</mn></mrow></msup><mo>−</mo><msup><mrow><mn>10</mn></mrow><mrow><mo>−</mo><mn>3</mn></mrow></msup></math></span>.</div></div>","PeriodicalId":285,"journal":{"name":"Computer Physics Communications","volume":"312 ","pages":"Article 109572"},"PeriodicalIF":7.2000,"publicationDate":"2025-03-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Computer Physics Communications","FirstCategoryId":"101","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S001046552500075X","RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q1","JCRName":"COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS","Score":null,"Total":0}
引用次数: 0
Abstract
The 2-D 3-T heat conduction equations can approximate the propagation of energy in the material, as well as the energy exchange process of electrons, ions and photons. In many cases, the computational time to solve these equations accounts for a large proportion (more than 80%) of that of the entire simulation of radiation hydrodynamics. PINNs is a promising way to solve partial differential equations (PDEs). Although numerical methods have been successful in solving 2-D 3-T heat conduction equations, the PINNs method also has some advantages, such as mesh-free, suitable to high dimension and complex domain problems. But the original PINNs cannot solve the 2-D 3-T heat conduction equations to a reasonable precision. This work aims to explore which techniques need to be added to PINNs and to what extent it can address the challenges posed by strong nonlinearity and multi-scale phenomena present in 2-D 3-T heat conduction equations. Hence, we adopt guaranteed positive constraint to the outputs so that the network can be trained, give a relatively large weight to the initial loss and Dirichlet boundary loss, take the logarithm of the initial loss, use transfer learning and Fourier feature embedding to improve accuracy. We call our improved approach 2D 3T PINNs. Numerical experiments show that the relative and the absolute error between the 2D 3T PINNs prediction and the reference solution is of the order of .
期刊介绍:
The focus of CPC is on contemporary computational methods and techniques and their implementation, the effectiveness of which will normally be evidenced by the author(s) within the context of a substantive problem in physics. Within this setting CPC publishes two types of paper.
Computer Programs in Physics (CPiP)
These papers describe significant computer programs to be archived in the CPC Program Library which is held in the Mendeley Data repository. The submitted software must be covered by an approved open source licence. Papers and associated computer programs that address a problem of contemporary interest in physics that cannot be solved by current software are particularly encouraged.
Computational Physics Papers (CP)
These are research papers in, but are not limited to, the following themes across computational physics and related disciplines.
mathematical and numerical methods and algorithms;
computational models including those associated with the design, control and analysis of experiments; and
algebraic computation.
Each will normally include software implementation and performance details. The software implementation should, ideally, be available via GitHub, Zenodo or an institutional repository.In addition, research papers on the impact of advanced computer architecture and special purpose computers on computing in the physical sciences and software topics related to, and of importance in, the physical sciences may be considered.