A Diffuse Domain Approximation with Transmission-Type Boundary Conditions II: Gamma--Convergence
Toai Luong, Tadele Mengesha, Steven M. Wise, Ming Hei Wong
TL;DR
The paper analyzes a diffuse-domain method for PDEs on domains with interfaces by introducing a phase-field $\phi_\varepsilon$ and studying the energy functionals $F_\varepsilon$ that encode the diffuse problem. It proves that, for the Neumann case ($\kappa=0$), $F_\varepsilon$ $\Gamma$-converges to the sharp-interface functional $F_0$ as $\varepsilon\to0$ in any dimension, and that minimizers $u_\varepsilon$ converge strongly to the sharp-interface minimizer $u_0$ in $H^1(\Omega)$ up to subsequences. In 1D, the results extend to general Robin-type transmission ($\kappa\ge0$), with the interfacial term converging to the boundary contribution and strong $H^1(\Omega)$ convergence still holding. The analysis relies on weighted diffuse-interface lemmas, compactness, and explicit liminf/limsup constructions, providing a rigorous justification of the diffuse-domain approach for two-sided transmission problems. The findings support reliable numerical and analytical use of DDMs in complex geometries and highlight dimensional nuances, especially in 1D where the Robin case is tractable.
Abstract
Diffuse domain methods (DDMs) have gained significant attention for solving partial differential equations (PDEs) on complex geometries. These methods approximate the domain by replacing sharp boundaries with a diffuse layer of thickness $\varepsilon$, which scales with the minimum grid size. This reformulation extends the problem to a regular domain, incorporating boundary conditions via singular source terms. In this work, we analyze the convergence of a DDM approximation problem with transmission-type Neumann boundary conditions. We prove that the energy functional of the diffuse domain problem $Γ$--converges to the energy functional of the original problem as $\varepsilon \to 0$. Additionally, we show that the solution of the diffuse domain problem strongly converges in $H^1(Ω)$, up to a subsequence, to the solution of the original problem, as $\varepsilon \to 0$.
