2023/07/23 by Zheyuan Hu, Khemraj Shukla, George Em Karniadakis +1 · 1 voice · 160 citations
Engineering · Physics and Astronomy · #Artificial intelligence #Artificial neural network #Computer science #Curse #Curse of dimensionality #Machine learning #Model Reduction and Neural Networks #Nuclear Engineering Thermal-Hydraulics #Nuclear reactor physics and engineering #Physics #Sociology #Statistical physics
paper · pdf · open access · doi:10.1016/j.neunet.2024.106369
published in Neural Networks 176, 106369 (Elsevier BV)
openalex publication_date 2024/05/07 · openalex created_date 2025/10/10 · openalex updated_date 2026/08/03
The curse-of-dimensionality taxes computational resources heavily with exponentially increasing computational cost as the dimension increases. This poses great challenges in solving high-dimensional partial differential equations (PDEs), as Richard E. Bellman first pointed out over 60 years ago. While there has been some recent success in solving numerical PDEs in high dimensions, such computations are prohibitively expensive, and true scaling of general nonlinear PDEs to high dimensions has never been achieved. We develop a new method of scaling up physics-informed neural networks (PINNs) to solve arbitrary high-dimensional PDEs. The new method, called Stochastic Dimension Gradient Descent (SDGD), decomposes a gradient of PDEs' and PINNs' residual into pieces corresponding to different dimensions and randomly samples a subset of these dimensional pieces in each iteration of training PINNs. We prove theoretically the convergence and other desired properties of the proposed method. We demonstrate in various diverse tests that the proposed method can solve many notoriously hard high-dimensional PDEs, including the Hamilton-Jacobi-Bellman (HJB) and the Schrödinger equations in tens of thousands of dimensions very fast on a single GPU using the PINNs mesh-free approach. Notably, we solve nonlinear PDEs with nontrivial, anisotropic, and inseparable solutions in less than one hour for 1000 dimensions and in 12 h for 100,000 dimensions on a single GPU using SDGD with PINNs. Since SDGD is a general training methodology of PINNs, it can be applied to any current and future variants of PINNs to scale them up for arbitrary high-dimensional PDEs.