Residual-Based Adaptive Refinement in Hybrid Variational-Collocation PINNs for Singular Elliptic PDEs with Mixed Dirichlet-Neumann Boundary Conditions
Abstract
Solutions of second-order elliptic PDE problems on nonconvex or cracked domains with mixed Dirichlet-Neumann conditions typically lie outside H2(Ω). Near re-entrant corners, crack tips and points where the boundary condition changes type, the gradient grows like rλ−1 with λ<1. We study residual-based adaptive refinement in hybrid variational-collocation PINNs. A multiscale Fourier-feature network is trained on a composite loss. The loss combines a variational residual tested against Legendre bubble functions, a strong-form collocation term and penalties for the Dirichlet and Neumann data. The variational term needs only first derivatives of the network, whereas the collocation term requires second derivatives and provides a residual that can be evaluated at any point. After each training stage, candidate points selected by this residual are added to the collocation set. The method is tested on five two dimensional benchmarks, namely a Mode~III crack, a cross-shaped domain with four re-entrant corners with unequal amplitude, a reaction-diffusion problem with four localized peaks, a monotone semilinear problem with a boundary layer and the Motz problem. Compared with uniform point insertion under the same collocation and iteration budget, residual-based refinement reduced the relative L2 error by up to a factor of $32$, the relative L∞ error by up to a factor of $38$ and the relative H1 semi-norm error by up to a factor of $15$.
Related
- Residual-based Adaptive Refinement of hp-Variational PINNs for elliptic PDEs with geometric singularities and localized features
- Uniqueness results for the exterior Dirichlet and Neumann problem in two dimensions for the Euler-Lagrange equations of variational integrals with linear growth
- Residual Calibration via Local Feature-Space Refinement
- A volume-conserving ALE continuous Galerkin method with residual-based Dirichlet--Neumann fluxes for nonlinear free-surface waves
- Variational Inference for Functional Data Clustering via Dirichlet Process Mixtures with Correlated Errors
BibTeX
@misc{giri2026residualbased,
title = {{Residual-Based Adaptive Refinement in Hybrid Variational-Collocation PINNs for Singular Elliptic PDEs with Mixed Dirichlet-Neumann Boundary Conditions}},
author = {Somnath Giri and Muslim Malik and Rishikesh Yadav},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.26458},
eprint = {2610.26458},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.26458},
note = {CC BY-NC-SA 4.0},
}
RIS
TY - UNPB
TI - Residual-Based Adaptive Refinement in Hybrid Variational-Collocation PINNs for Singular Elliptic PDEs with Mixed Dirichlet-Neumann Boundary Conditions
AU - Somnath Giri
AU - Muslim Malik
AU - Rishikesh Yadav
PY - 2026
DA - 2026/10/23
PB - arXiv
DO - 10.48550/arxiv.2610.26458
UR - https://arxiv.org/abs/2610.26458
AB - Solutions of second-order elliptic PDE problems on nonconvex or cracked domains with mixed Dirichlet-Neumann conditions typically lie outside H^2(Ω). Near re-entrant corners, crack tips and points where the boundary condition changes type, the gradient grows like r^λ−1 with λ<1. We study residual-based adaptive refinement in hybrid variational-collocation PINNs. A multiscale Fourier-feature network is trained on a composite loss. The loss combines a variational residual tested against Legendre bubble functions, a strong-form collocation term and penalties for the Dirichlet and Neumann data. The variational term needs only first derivatives of the network, whereas the collocation term requires second derivatives and provides a residual that can be evaluated at any point. After each training stage, candidate points selected by this residual are added to the collocation set. The method is tested on five two dimensional benchmarks, namely a Mode~III crack, a cross-shaped domain with four re-entrant corners with unequal amplitude, a reaction-diffusion problem with four localized peaks, a monotone semilinear problem with a boundary layer and the Motz problem. Compared with uniform point insertion under the same collocation and iteration budget, residual-based refinement reduced the relative L^2 error by up to a factor of $32$, the relative L^∞ error by up to a factor of $38$ and the relative H^1 semi-norm error by up to a factor of $15$.
ER -