Eddy Current Testing
Theory and Physics
Principles of Eddy Current Testing (ECT)
Professor, how does eddy current testing work?
An alternating current is passed through a probe coil, inducing eddy currents in the test object. If a defect (crack, corrosion) is present, the flow of eddy currents is disturbed, causing a change in the probe's impedance.
The location and size of the defect are estimated from the impedance change $\Delta Z$.
So the inspection depth changes with frequency, right?
Correct. The standard penetration depth $\delta = \sqrt{2/(\omega\mu\sigma)}$ serves as a guideline for inspection depth. Low frequencies inspect deeper regions, high frequencies inspect the surface. In multi-frequency ECT, multiple frequencies are used simultaneously to improve depth resolution.
Summary
- Impedance change $\Delta Z$ — Principle of defect detection
- Inspection depth $\approx \delta$ — Controlled by frequency
- Non-contact & high-speed — Applicable to inline inspection on production lines
Crack Detection in Jet Engine Blades—How Eddy Current Testing Finds "Invisible Flaws"
Have you ever seen a scene where a maintenance technician runs a small probe over the surface of an aircraft engine blade during inspection? That's the eddy current testing (ECT) field in action. When an alternating magnetic field is applied, eddy currents flow on the metal surface beneath the coil. However, if a crack or corrosion is present, the eddy current flow is disturbed, changing the coil's impedance. By picking up this minute impedance change, invisible cracks smaller than 0.1mm can be detected non-contact and non-destructively. Unlike X-rays or UT (ultrasonic testing), ECT has the strengths of being "specialized for conductor surface/subsurface flaws" and "probes are small and light, easily mounted on robots," making it an indispensable inspection method in the aerospace, nuclear, and automotive industries.
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 a 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. An IH (induction heating) cooking heater also uses the same principle, where a high-frequency magnetic field change induces 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. 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. An electromagnet operates on this principle, passing current through a coil to create a strong magnetic field. A smartphone speaker also applies 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$: 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—the charged sheet (electric charge) radiates electric field lines outward, 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 grounded in 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—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, 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
- Isotropy assumption: For anisotropic materials (e.g., rolling direction of silicon steel sheets), direction-specific property definitions are needed
- Non-applicable cases: Additional constitutive laws are required 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 unit 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
ECT Simulation Using FEM
What is the significance of simulating ECT using FEM?
It is used for probe design optimization, defect signal prediction, and optimization of inspection conditions (frequency, lift-off). 3D eddy current analysis calculates the interaction between the probe and the defect.
The probe is scanned, and the impedance change $\Delta Z$ at each position is plotted (Lissajous figure).
Isn't the computational cost high for 3D?
Combine the A-V method with edge elements and use fine meshes only near the defect. Also, reduce computational scale using symmetry or Fourier decomposition. COMSOL's Parametric Sweep can automatically scan probe positions.
Summary
- 3D Eddy Current Analysis — Probe-defect interaction
- Impedance Change Calculation — Defect determination via Lissajous figure
- Parametric Sweep — Automatic sweep of probe position
Impedance Plane Diagram—The Technique of Deciphering the "Fingerprint" of Eddy Current Testing
Eddy current testing signals are often displayed as an "impedance plane diagram (Lissajous figure)." Plotting the resistance component on the X-axis and the reactance component on the Y-axis, cracks, conductivity changes, and lift-off (probe lift) each draw vectors in different directions. Experienced inspectors read this "direction and length of the trajectory" to determine the depth and type of the flaw. Recently, research using machine learning for automatic determination has become active, with cases where convolutional neural networks show judgment accuracy equal to or better than veteran inspectors. Signal data from CAE simulations for "flaw-free cases" and "cases with various defects" are used as training data.
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. BEM: Naturally handles infinite domains (open region problems). Hybrid FEM-BEM is also effective.
Nonlinear Convergence (Magnetic Saturation)
Nonlinearity of the B-H curve is handled by the Newton-Raphson method. Residual criterion: $||R||/||R_0|| < 10^{-4}$ is typical.
Frequency Domain Analysis
Reduced to a steady-state problem by assuming time-harmonic conditions. Requires complex number operations, but broadband characteristics are obtained via time-domain analysis.
Time Domain Time Step
A time step smaller than 1/20 of the highest frequency component is required. Implicit time integration allows larger steps but requires attention to accuracy.
Choosing Between Frequency Domain and Time Domain
Frequency domain analysis is like "tuning a radio to a specific frequency"—it can efficiently calculate the response at a single frequency. Time domain analysis is like "recording all channels simultaneously"—it can reproduce transient phenomena containing all frequency components, but computational cost is high.
Practical Guide
Application in Practice
Related Topics
なった
詳しく
報告