Boundary Element Method for Electrostatic Fields
Theory and Physics
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.
Physical Meaning of Each Term
- Electric Field Term $\nabla \times \mathbf{E} = -\partial \mathbf{B}/\partial t$: Faraday's law of electromagnetic induction. A time-varying magnetic flux density generates an electromotive force. 【Everyday Example】A bicycle dynamo (generator) generates voltage in a nearby coil by rotating a magnet—a direct application of this law that a changing magnetic field induces an electric field. Induction cooking (IH) heaters also use the same principle, where high-frequency magnetic field changes induce eddy currents in the pot bottom, heating it via Joule heating.
- Magnetic Field Term $\nabla \times \mathbf{H} = \mathbf{J} + \partial \mathbf{D}/\partial t$: Ampère-Maxwell's law. Electric current and displacement current generate a magnetic field. 【Everyday Example】When current flows through a wire, a magnetic field is created around it—this is Ampère's law. Electromagnets operate on this principle, passing current through a coil to create a strong magnetic field. Smartphone speakers also apply this law: current → magnetic field → force on the diaphragm. At high frequencies (e.g., GHz-band antennas), the displacement current $\partial D/\partial t$ cannot be ignored, describing electromagnetic wave radiation.
- Gauss's Law $\nabla \cdot \mathbf{D} = \rho_v$: Indicates that electric charge is the divergence source of electric flux. 【Everyday Example】Rubbing a plastic sheet on hair makes hair stand up due to static electricity—electric field lines radiate outward from the charged sheet (charge), exerting force on the light hair. In capacitor design, the electric field distribution between electrodes is calculated using this law. ESD (electrostatic discharge) countermeasures are also based on electric field analysis derived from Gauss's law.
- Magnetic Flux Conservation $\nabla \cdot \mathbf{B} = 0$: Expresses the non-existence of magnetic monopoles. 【Everyday Example】Even if you cut a bar magnet in half, you cannot create a magnet with only an N pole or only an S pole—they always exist as an N-S pair. This means magnetic field lines form "closed loops with no start or end points." In numerical analysis, to satisfy this condition, the formulation using vector potential $\mathbf{B} = \nabla \times \mathbf{A}$ is used, automatically guaranteeing magnetic flux conservation.
Assumptions and Applicability Limits
- Linear Material Assumption: Permeability and permittivity do not depend on magnetic/electric field strength (nonlinear B-H curve needed in saturation region)
- Quasi-Static Approximation (Low Frequency): Displacement current term can be ignored ($\omega \varepsilon \ll \sigma$). Common in eddy current analysis
- 2D Assumption (Cross-Sectional Analysis): Effective when current direction is uniform and end effects can be ignored
- Isotropic Assumption: Direction-specific property definitions needed for anisotropic materials (e.g., rolling direction of silicon steel sheets)
- Non-Applicable Cases: Additional constitutive laws needed for plasma (ionized gas), superconductors, nonlinear optical materials
Dimensional Analysis and Unit Systems
| Variable | SI Unit | Notes / Conversion Memo |
|---|---|---|
| Magnetic Flux Density $B$ | T (Tesla) | 1T = 1 Wb/m². Permanent Magnets: 0.2–1.4T |
| Magnetic Field Strength $H$ | A/m | Horizontal axis of B-H curve. Conversion from CGS Oe (Oersted): 1 Oe = 79.577 A/m |
| Current Density $J$ | A/m² | Calculated from conductor cross-section and total current. Note non-uniform distribution due to skin effect |
| Permeability $\mu$ | H/m | $\mu = \mu_0 \mu_r$. In vacuum $\mu_0 = 4\pi \times 10^{-7}$ H/m |
| Electrical Conductivity $\sigma$ | S/m | Copper: approx. 5.96×10⁷ S/m. Decreases with temperature rise |
Numerical Methods and Implementation
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.
Edge Elements (Nedelec Elements)
Elements specialized for electromagnetic field analysis. Automatically guarantee continuity of tangential components and eliminate spurious modes. Standard for 3D high-frequency analysis.
Nodal Elements
Used for scalar potential formulations. Effective for scalar potential methods in magnetostatics and electrostatic field analysis.
FEM vs BEM (Boundary Element Method)
FEM: Handles nonlinear materials and inhomogeneous media. Related Topics
なった
詳しく
報告