Logarithmic Resonance at Mixed Boundary Junctions in Gradient-Dependent Semilinear Equationsv

Abstract

The regularity of solutions to elliptic partial differential equations degrades severely at mixed Dirichlet-Neumann boundary junctions, characterized classically by an $\mathcal{O}(r^{1/2})$ leading singular function. While this linear behavior is well documented, the introduction of gradient-dependent semilinear perturbations alters the local asymptotic profile. This article proves that gradient penalties scaling linearly with $|\nabla u|$ induce a highly localized $\mathcal{O}(r^{-1/2})$ source term that resonates with the half-integer spectrum of the principal homogeneous differential operator. This non-orthogonal resonance causes standard separable polynomial assumptions to fail at order $\mathcal{O}(r^{3/2})$. We establish a generalized resonance theorem that forces the emergence of a logarithmic anomaly, providing the exact analytical formulation of the resulting $r^{3/2}\ln (r)$ profile alongside rigorous local Sobolev regularity bounds for the remainder. By calculating the exact logarithmic coefficients and angular offsets for ${\ell }_1$, ${\ell }_2$, ${\ell }_{\infty }$, and arbitrary ${\ell }_q$-norm penalties, we demonstrate the universality of this obstruction. Finally, we formalize the corresponding enriched continuous Galerkin space (XFEM), establish its quasi-optimality, and present finite element experiments that confirm the predicted localized pollution, recover the predicted logarithmic coefficients across the ${\ell }_1$, ${\ell }_2$, and ${\ell }_{\infty }$ penalties on a curvature-free flat junction, and show that resolving the junction recovers optimal degree-of-freedom efficiency in standard finite element solvers; we close by outlining the targeted software architectures required for a fully enriched implementation.

Work in progress