Boundary Element Method for Electrostatic Fields
Boundary Element Method for Electrostatic Fields: Theoretical Foundations
Electrostatic Field BEM
Professor, can BEM be used for electrostatic fields too?
BEM is optimal for open-space problems. FEM requires meshing the surrounding space, but BEM only needs the surfaces of conductors/dielectrics.
Green's function $G = 1/(4\pi r)$ (3D). Solve for surface potential and normal electric field as unknowns.
Infinite domains can be handled with just surface meshes!
Correct. It demonstrates its power in analyzing electric fields around power lines, shielding effects of lightning surges, and open-space analysis for EMC problems.
Comparison of FEM and BEM
Summary
- Open space can be solved with only surface meshes — The greatest advantage of BEM
- Complementary with FEM — FEM for closed domains, BEM for open space
- FEM-BEM Coupling — Hybrid method with internal FEM + external BEM
The Origin of Boundary Integral Equations—The Lonely Discovery of 19th-Century Mathematician Green
George Green (1793–1841), who published "Green's theorem," the theoretical foundation of electrostatic BEM, was the son of a Nottingham baker and self-taught in mathematics. The paper he self-published in 1828 was largely ignored during his lifetime, but was rediscovered posthumously by Lord Kelvin and became a cornerstone of electromagnetism. It is an astonishing story that his Green's function continues to live at the heart of BEM nearly 200 years later.
Computational Methods for Boundary Element Method for Electrostatic Fields
BEM Discretization
Discretize the surface with triangular/quadrilateral elements. Construct a system of equations with potential $\phi$ and normal electric flux density $D_n$ on each element as unknowns:
Influence matrices $[H], [G]$ are surface integrals of the Green's function. Dense matrices, so $O(N^2)$ memory.
FMM Acceleration
For large-scale problems, accelerate to $O(N\log N)$ using FMM (Fast Multipole Method). Implemented in FastCap and COMSOL's BEM module.
Summary
- $[H]\{\phi\} = [G]\{D_n\}$ — Fundamental BEM equation
- Dense matrix → Accelerated by FMM
BEM's "Green's Function" is Not Straightforward
The core of the Boundary Element Method is the Green's function—a mathematical expression for "the influence a point charge placed at a certain point has on the surrounding potential." For electrostatic fields, it takes the relatively simple form $1/r$ (inverse distance), but it becomes extremely complex for anisotropic dielectrics or multilayer media. For multilayer substrates in semiconductor packages, different Green's functions are needed for each dielectric layer, and numerical techniques to accurately perform these integrals are at the forefront of research.