Table of Contents
Fetching ...

Trace Regularity PINNs: Enforcing $\mathrm{H}^{\frac{1}{2}}(\partial Ω)$ for Boundary Data

Doyoon Kim, Junbin Song

TL;DR

This work introduces TRPINN, a physics-informed neural network that enforces boundary data in the Sobolev–Slobodeckij space $\mathrm{H}^{1/2}(\partial\Omega)$ for divergence-form elliptic PDEs. By discretizing only the essential portion of the $\mathrm{H}^{1/2}$ semi-norm and avoiding problematic denominators, TRPINN achieves boundary regularity consistent with theory while maintaining comparable computational cost to standard PINNs. Theoretical results show that zero loss implies convergence in $\mathrm{H}^1(\Omega)$, and NTK analysis indicates faster convergence rates relative to conventional PINNs, especially for high-frequency boundary data. Numerical experiments on Laplace problems with oscillatory and sharply peaked boundary data demonstrate substantial accuracy gains and robustness to weight choices, validating the practical impact of enforcing $\mathrm{H}^{1/2}(\partial\Omega)$ boundary regularity in PINN frameworks.

Abstract

We propose an enhanced physics-informed neural network (PINN), the Trace Regularity Physics-Informed Neural Network (TRPINN), which enforces the boundary loss in the Sobolev-Slobodeckij norm $H^{1/2}(\partial Ω)$, the correct trace space associated with $H^1(Ω)$. We reduce computational cost by computing only the theoretically essential portion of the semi-norm and enhance convergence stability by avoiding denominator evaluations in the discretization. By incorporating the exact $H^{1/2}(\partial Ω)$ norm, we show that the approximation converges to the true solution in the $H^{1}(Ω)$ sense, and, through Neural Tangent Kernel (NTK) analysis, we demonstrate that TRPINN can converge faster than standard PINNs. Numerical experiments on the Laplace equation with highly oscillatory Dirichlet boundary conditions exhibit cases where TRPINN succeeds even when standard PINNs fail, and show performance improvements of one to three decimal digits.

Trace Regularity PINNs: Enforcing $\mathrm{H}^{\frac{1}{2}}(\partial Ω)$ for Boundary Data

TL;DR

This work introduces TRPINN, a physics-informed neural network that enforces boundary data in the Sobolev–Slobodeckij space for divergence-form elliptic PDEs. By discretizing only the essential portion of the semi-norm and avoiding problematic denominators, TRPINN achieves boundary regularity consistent with theory while maintaining comparable computational cost to standard PINNs. Theoretical results show that zero loss implies convergence in , and NTK analysis indicates faster convergence rates relative to conventional PINNs, especially for high-frequency boundary data. Numerical experiments on Laplace problems with oscillatory and sharply peaked boundary data demonstrate substantial accuracy gains and robustness to weight choices, validating the practical impact of enforcing boundary regularity in PINN frameworks.

Abstract

We propose an enhanced physics-informed neural network (PINN), the Trace Regularity Physics-Informed Neural Network (TRPINN), which enforces the boundary loss in the Sobolev-Slobodeckij norm , the correct trace space associated with . We reduce computational cost by computing only the theoretically essential portion of the semi-norm and enhance convergence stability by avoiding denominator evaluations in the discretization. By incorporating the exact norm, we show that the approximation converges to the true solution in the sense, and, through Neural Tangent Kernel (NTK) analysis, we demonstrate that TRPINN can converge faster than standard PINNs. Numerical experiments on the Laplace equation with highly oscillatory Dirichlet boundary conditions exhibit cases where TRPINN succeeds even when standard PINNs fail, and show performance improvements of one to three decimal digits.
Paper Structure (18 sections, 6 theorems, 52 equations, 8 figures, 5 tables)

This paper contains 18 sections, 6 theorems, 52 equations, 8 figures, 5 tables.

Key Result

Lemma 2.1

Let $\Omega \subset \mathbb{R}^2$. For any fixed $\delta >0$, for all $g \in \mathrm{H}^{\frac{1}{2}}(\partial \Omega)$, where the semi-norm equivalence depends on $\delta$.

Figures (8)

  • Figure 1: The integration regions for (a)–(c) correspond to the two-dimensional case ($d = 2$). (a) The integration region corresponding to Eq. \ref{['eq01']}. (b) The integration region used for the numerical integration of the boundary loss in the DRM with AQR method proposed in liu2023deep, where the squares composed of 9 small cells correspond to $\mathcal{E}_D$, and the small cells correspond to $\mathcal{E}_D^k$ (see liu2023deep). (c) The integration region for the numerical integration of the boundary loss in the proposed TRPINN method.
  • Figure 2: The left column shows the magnitudes of the eigenvalues for PINN (blue line) and TRPINN (orange line), computed by extracting the largest one million. The right column shows the absolute value of their differences. From top to bottom, the boundary sampling methods tested are linspace, randomized, and uniform random.
  • Figure 3: Harmonic extension of $\sin(20\theta)$. This is used as the exact solution to compute the relative error.
  • Figure 4: Top-left panel shows the boundary prediction of the vanilla PINN with weights $[1,1]$; all remaining panels show TRPINN. The weights follow the same order as the panels—from the top-left to the bottom-right—as listed in \ref{['V=20_rate1']}.
  • Figure 5: Relative $\mathrm{H}^1(\Omega$) error. Top row: 3 layers; middle row: 5 layers; bottom row: 7 layers. Left column: 30 units; center column: 50 units; right column: 100 units. The x-axis denotes iterations: 0--49,999 correspond to Adam training, while iterations $\ge 50,000$ correspond to L-BFGS training.
  • ...and 3 more figures

Theorems & Definitions (13)

  • Lemma 2.1
  • Theorem 5.1
  • Remark 5.2
  • Example 5.3
  • Definition 6.1
  • Remark 6.2
  • Remark 6.3
  • Lemma 6.4
  • Remark 6.5
  • Remark 6.6
  • ...and 3 more