A technique that reformulates boundary value problems for linear partial differential equations as integral equations posed only on the boundary, using layer potentials and fundamental solutions to reduce dimensionality and incorporate radiation conditions.

Principle

Represent the solution in the domain as single and/or double layer potentials built from a Green's function; apply trace and jump relations on the boundary to obtain Fredholm integral equations whose solvability yields the interior or exterior solution.

Demonstration

For the exterior Helmholtz Dirichlet problem, represent the field as a single‑layer potential with unknown density on the obstacle boundary; the resulting boundary integral equation is solved for the density and the potential reconstructs the scattered field satisfying Sommerfeld radiation.

Misapplication

Applying the boundary integral formulation at frequencies that are interior eigenfrequencies without using combined field formulations, or discretizing singular kernels with naive quadrature, which leads to ill‑conditioned systems and spurious solutions.

Consequence

Provides dimension reduction, exact enforcement of radiation conditions for unbounded domains, spectral insight for fast solvers (fast multipole, H‑matrix) and high accuracy for smooth boundaries with appropriate quadrature and regularization.

Reversal

Domain methods such as finite element or finite difference discretizations that solve volume equations are the reversal: they handle arbitrary materials and nonlinearities more directly but require artificial truncation or absorbing conditions for unbounded domains.

Boundary

Requires knowledge of an explicit fundamental solution and sufficient boundary regularity; best suited to linear constant‑coefficient operators or problems where layer potentials are available, and less directly applicable to strongly nonlinear PDEs or heterogeneous media without adapted Green's functions.

Semantic Tension

Tension arises between boundary integral approaches (dimension reduction, global kernels) and local domain discretizations (FEM): tradeoffs include handling of inhomogeneities, adaptivity, conditioning at resonant frequencies, and implementation of fast algorithms.

Synthesis

The boundary integral method converts boundary value problems into boundary integral equations through layer potentials and jump relations, enabling efficient, accurate treatment of exterior and transmission problems when fundamental solutions and careful discretization are available.