Table of Contents
Fetching ...

Spatio-Temporal Uncertainty-Modulated Physics-Informed Neural Networks for Solving Hyperbolic Conservation Laws with Strong Shocks

Darui Zhao, Ze Tao, Fujun Liu

Abstract

Physics-Informed Neural Networks (PINNs) frequently encounter difficulties in accurately resolving shock waves within high-speed compressible flows, a failure largely attributed to the "gradient pathology" arising from extreme stiffness at discontinuities. To overcome this limitation, we propose the Spatio-Temporal Uncertainty-Modulated PINN (UM-PINN), a probabilistic framework that reinterprets the training process as a multi-task learning problem governed by homoscedastic aleatoric uncertainty. By integrating a gradient-based spatial mask with learnable variance parameters, our method dynamically balances the conflicting contributions of Partial Differential Equation (PDE) residuals and initial conditions across the spatiotemporal domain, further stabilized by Quasi-Monte Carlo Sobol sampling. We validate the framework against challenging benchmarks, including the one-dimensional (1D) Sod shock tube, the high-frequency Shu-Osher problem, and the complex two-dimensional (2D) Riemann interaction, where standard gradient-based weighting schemes typically fail. Experimental results demonstrate that UM-PINN achieves orders of magnitude improvement in accuracy and shock resolution compared to baseline methods, establishing a robust new paradigm for mesh-free Computational Fluid Dynamics in hyperbolic systems.

Spatio-Temporal Uncertainty-Modulated Physics-Informed Neural Networks for Solving Hyperbolic Conservation Laws with Strong Shocks

Abstract

Physics-Informed Neural Networks (PINNs) frequently encounter difficulties in accurately resolving shock waves within high-speed compressible flows, a failure largely attributed to the "gradient pathology" arising from extreme stiffness at discontinuities. To overcome this limitation, we propose the Spatio-Temporal Uncertainty-Modulated PINN (UM-PINN), a probabilistic framework that reinterprets the training process as a multi-task learning problem governed by homoscedastic aleatoric uncertainty. By integrating a gradient-based spatial mask with learnable variance parameters, our method dynamically balances the conflicting contributions of Partial Differential Equation (PDE) residuals and initial conditions across the spatiotemporal domain, further stabilized by Quasi-Monte Carlo Sobol sampling. We validate the framework against challenging benchmarks, including the one-dimensional (1D) Sod shock tube, the high-frequency Shu-Osher problem, and the complex two-dimensional (2D) Riemann interaction, where standard gradient-based weighting schemes typically fail. Experimental results demonstrate that UM-PINN achieves orders of magnitude improvement in accuracy and shock resolution compared to baseline methods, establishing a robust new paradigm for mesh-free Computational Fluid Dynamics in hyperbolic systems.

Paper Structure

This paper contains 28 sections, 46 equations, 18 figures, 2 tables.

Figures (18)

  • Figure 1: Schematic architecture of the Uncertainty-Modulated Physics-Informed Neural Network (UM-PINN). The network predicts the primitive variables $(\rho, u, p)$ from spatio-temporal coordinates $(t, x)$. The core innovation is the Uncertainty-Modulated (UM) Module
  • Figure 2: Comparison of predicted profiles for density, velocity, and pressure against analytical solutions for the 1D Sod shock tube at $t=0.5$
  • Figure 3: Quantitative error metrics (RMSE, Relative $L_2$ error, and $L_\infty$ error) comparing the Baseline PINN and the proposed UM-PINN for the Sod shock tube problem.
  • Figure 4: Comparative analysis of the Shu-Osher problem at $t=1.80$. (a) The proposed UM-PINN (left) successfully resolves the intricate oscillations and matches the reference solution with high accuracy. (b)The Standard Baseline PINN (right) fails to capture high-frequency entropy waves due to spectral bias.
  • Figure 5: Qualitative comparison of the steady-state Riemann problem results. (a) Density ($\rho$) fields; (b) Pressure ($p$) fields. In each panel, the top row displays the sharp interfaces captured by the proposed UM-PINN, while the bottom row shows the results from the Baseline method, which suffer from numerical diffusion.
  • ...and 13 more figures

Theorems & Definitions (2)

  • Remark 1
  • Remark 2