BEM, (engl. boundary element method), numerische Behandlung der Randintegralmethode zur Lösung partieller Differentialgleichungen.

Sie entsteht unter Verwendung eines Finite-Elemente-Ansatzes über dem Integrationsbereich. Dabei wird die gesuchte Funktion ersetzt durch eine Approximation in einem endlichdimensionalen Unterraum der Basisfunktionen {(ϕ_i)}. Man versucht dann entweder, die entstehenden Gleichungen in allen Knotenpunkten der finiten Elemente zur erfüllen (Kollokation), oder durch zusätzliche Integration mit den Basisfunktionen über dem Rand durch die Galerkin-Methode zu lösen.

Bei einer Elementgröße h entstehen bei Problemen im ℝ^d mit der Randelementmethode lineare Gleichungssysteme der Ordnung O(h^1−d), während sie bei einem Finite-Elemente-Ansatz der ursprünglichen Differentialgleichung O(h^−d) ist. Allerdings sind die Matrizen bei der BEM i. allg. voll besetzt.