Analysis of Dielectric Materials
Theory and Physics
What is a Dielectric?
Professor, are dielectrics the same as insulators?
They are the same in the sense that they do not conduct electricity, but the term 'dielectric' focuses on the fact that they polarize in an electric field.
$\mathbf{P}$: Polarization vector, $\varepsilon_r$: Relative permittivity. Polarization weakens the internal electric field.
Major Dielectric Materials
| Material | $\varepsilon_r$ | Applications |
|---|---|---|
| Vacuum | 1.0 | Reference |
| Air | 1.0006 | ≈Vacuum |
| Teflon (PTFE) | 2.1 | High-frequency substrates |
| Polyimide | 3.4 | Flexible substrates |
| Epoxy (FR-4) | 4.5 | PCB substrates |
| SiO₂ | 3.9 | Semiconductor gate oxide films |
| BaTiO₃ | 1000 to 10000 | Ceramic capacitors |
| Water | 80 | Biological/Chemical |
The permittivity of BaTiO₃ is orders of magnitude different.
Ferroelectrics possess spontaneous polarization, resulting in extremely high $\varepsilon_r$. This contributes to the miniaturization of MLCCs. However, $\varepsilon_r$ varies with temperature, electric field, and frequency.
Summary
- $\mathbf{D} = \varepsilon_0 \varepsilon_r \mathbf{E}$ — Constitutive relation for linear dielectrics
- $\varepsilon_r$ is material-specific — Ranges from 1 (vacuum) to 10000 (ferroelectrics)
- Temperature, frequency, and electric field dependence — Note nonlinear effects
The Four Mechanisms of Dielectric Polarization—Which One Dominates Depends on Frequency
There are four mechanisms by which dielectrics polarize in response to an electric field: ① Electronic polarization (GHz~THz), ② Ionic polarization (IR band), ③ Orientational polarization (MHz~GHz), and ④ Interfacial polarization (low frequency). The dominant mechanism varies by material, appearing as characteristic "steps" in the frequency spectrum of the permittivity. Looking at a material's permittivity graph allows one to infer its molecular structure—permittivity is like a fingerprint of the material.
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) produces voltage in a nearby coil by rotating a magnet—a direct application of this law stating 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's law. 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$ becomes non-negligible, describing electromagnetic wave radiation.
- Gauss's law $\nabla \cdot \mathbf{D} = \rho_v$: States that electric charge is the source of divergence of electric flux. 【Everyday Example】Rubbing a plastic sheet on hair makes hair stand up due to static electricity—electric field lines radiate 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 following Gauss's law.
- Magnetic flux conservation $\nabla \cdot \mathbf{B} = 0$: Indicates the absence 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 a pair. This means magnetic field lines form "closed loops with no start or end points." In numerical analysis, the formulation using vector potential $\mathbf{B} = \nabla \times \mathbf{A}$ is used to satisfy this condition, 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 edge effects can be ignored
- Isotropic assumption: Direction-specific property definitions needed for anisotropic materials (e.g., rolling direction of silicon steel)
- Non-applicable cases: Additional constitutive relations 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 to 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
Dielectrics in FEM
For dielectrics, simply assign $\varepsilon_r$ to the FEM elements. At interfaces between different dielectrics, the continuity condition for $D_n$ is automatically satisfied.
Anisotropic Dielectrics
For crystals or laminated substrates, permittivity is a tensor:
Tensor input is possible in COMSOL or Maxwell.
Nonlinear Dielectrics
The $\varepsilon(E)$ dependence of ferroelectrics is handled by Newton-Raphson iteration. Same technique as magnetic material analysis with B-H curves.
Summary
- Linear: Assign $\varepsilon_r$ to elements
- Anisotropic: Tensor input
- Nonlinear: Newton-Raphson iteration
The "Frequency Dependence" of Permittivity—Debye Relaxation and Kramers-Kronig
The relative permittivity of a dielectric changes with frequency (dispersion). At low frequencies, dipole orientation keeps up with electric field changes, resulting in high permittivity; at high frequencies, it cannot keep up, and permittivity decreases. This relaxation phenomenon is described by the Debye model (1913). Furthermore, there is a constraint called the Kramers-Kronig relation between the real and imaginary parts of permittivity, allowing calculation of one from measurement data of the other. In CAE modeling of dielectric materials, obtaining permittivity data for the operating frequency band is the starting point.
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 analysis.
FEM vs BEM (Boundary Element Method)
FEM: Handles nonlinear materials and inhomogeneous media. BEM: Naturally handles infinite domains (open region problems). Hybrid FEM-BEM is also effective.
Nonlinear Convergence (Related Topics
なった
詳しく
報告