Gauss's Law (Electrostatics)
Theory and Physics
Overview
Professor, Gauss's Law is one of Maxwell's equations, right? How is it used in electrostatic field analysis?
Gauss's Law states that "the total electric flux through a closed surface is equal to the total charge enclosed within that surface." It corresponds to the first of Maxwell's equations.
In differential form, it is
Here, $\mathbf{D} = \varepsilon \mathbf{E}$ is the electric flux density, and $\rho$ is the volume charge density.
How does this relate to FEM electrostatic field analysis?
Substituting $\mathbf{E} = -\nabla\phi$ yields the Poisson equation.
The FEM solver discretizes and solves this. In other words, Gauss's Law is the fundamental equation underlying FEM electrostatic field analysis.
Analytical Solutions Using Symmetry
There are cases where the electric field can be found directly from Gauss's Law, right?
It can only be found when there is high symmetry.
| Symmetry | Gaussian Surface | Electric Field |
|---|---|---|
| Spherical (Point Charge) | Concentric Spherical Surface | $E = Q/(4\pi\varepsilon r^2)$ |
| Cylindrical (Line Charge) | Coaxial Cylindrical Surface | $E = \lambda/(2\pi\varepsilon r)$ |
| Planar (Surface Charge) | Rectangular Box | $E = \sigma/(2\varepsilon)$ |
These analytical solutions are essential for validating FEM results. Whether in COMSOL or Ansys Maxwell, the golden rule is to first confirm agreement with theoretical values for symmetric problems before moving on to complex ones.
"The Shape of the Closed Surface Can Be Anything" – The Exquisite Freedom of Gauss's Law
The interesting aspect of Gauss's Law is that "regardless of the shape of the closed surface, the integral value is determined solely by the enclosed charge." Sphere, cube, weird shape – all are OK. Utilizing this, calculations become dramatically simpler for problems with symmetry (like spherical charges or infinite cylindrical charges). For example, in the cross-section of a high-voltage cable – a core wire at the center and a shield on the outside – a structure with almost perfect cylindrical symmetry, the electric field distribution can be found in one shot by simply taking a Gaussian surface as a coaxial cylinder. Of course, FEM solves it numerically, but this intuitive solving method is invaluable in the field when you want to "check if it matches the analytical solution."
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 heat.
- Magnetic Field Term $\nabla \times \mathbf{H} = \mathbf{J} + \partial \mathbf{D}/\partial t$: Ampère-Maxwell law. Current and displacement current generate a magnetic field. 【Everyday Example】When current flows through a wire, a magnetic field is generated 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 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 rooted in Gauss's Law.
- Magnetic Flux Conservation $\nabla \cdot \mathbf{B} = 0$: Expresses the non-existence of magnetic monopoles. 【Everyday Example】Cutting a bar magnet in half does not create a magnet with only an N pole or only an S pole – N and S poles always exist as a 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 are independent of 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-Section 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-sectional area 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
Details of Numerical Methods
How is the Poisson equation derived from Gauss's Law formulated in FEM?
The Galerkin method for weak formulation is fundamental. Multiplying by a test function $w$ and integrating by parts yields
Discretization gives $[K]\{\phi\} = \{f\}$. This is exactly the same form as the stiffness equation in structural FEM.
So permittivity corresponds to Young's modulus, right?
Exactly. That's why if you have experience with structural FEM, electrostatic FEM is easy to understand.
Charge Calculation Using Gauss's Law
How do you calculate the charge on a conductor from FEM results in post-processing?
Use Gauss's Law to integrate the electric flux density on the conductor surface.
In COMSOL, this can be calculated directly using the surface integral post-processing feature. From this value, the capacitance $C = Q/V$ is calculated.
Extraction of Capacitance Matrix
How is the capacitance matrix for a multi-conductor system obtained?
Set conductor $j$ to 1V and all others to 0V, then calculate the induced charge on each conductor. Repeating this for all conductors yields the capacitance matrix $C_{ij}$. In Ansys Q3D, this operation is fully automated and widely used for parasitic capacitance extraction in PCB wiring.
The Choice of Gaussian Surface Affects FEM Accuracy
When using Gauss's Law to verify charge quantities in numerical analysis, "where to place the Gaussian surface" significantly changes the result accuracy. Setting a Gaussian surface that cuts through a region with coarse mesh leads to large integration errors and charge imbalance. In practice, the technique of "placing the Gaussian surface after ensuring a region with sufficiently fine mesh between the electrode surface and the analysis space" is used. Making the shape closer to a sphere or rectangular box can sometimes simplify integration, and mastering the use of Gauss's Law for post-analysis verification allows checking "whether the FEM has correctly calculated the charge."
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
Related Topics
なった
詳しく
報告