Wave equations in elastic media, stress-strain tensor relationships, Christoffel equations for acoustic anisotropy, diffraction phenomena, and advanced attenuation modeling that underpin Level III technique design and procedure development.
Elasticity Theory and the Wave Equation
The Wave Equation in Elastic Solids
At Level III, you must understand wave propagation not merely as a velocity-through-material concept but as a solution to the fundamental equations of motion in elastic media. This understanding is essential for developing procedures in anisotropic materials, evaluating technique limitations, and providing technical direction when standard approaches fail.
Stress, Strain, and Hooke's Law
In an elastic solid, the stress tensor σᵢⱼ relates to the strain tensor εₖₗ through the stiffness tensor Cᵢⱼₖₗ:
σᵢⱼ = Cᵢⱼₖₗ × εₖₗ
For an isotropic material, the 81 components of Cᵢⱼₖₗ reduce to just two independent elastic constants: the Lamé parameters λ and μ (where μ is the shear modulus G).
The longitudinal velocity derives from:
V_L = √((λ + 2μ) / ρ)
The shear velocity from:
V_S = √(μ / ρ)
The ratio V_L/V_S depends only on Poisson's ratio ν:
V_L/V_S = √((2 - 2ν) / (1 - 2ν))
For steel (ν ≈ 0.29): V_L/V_S ≈ 1.83
For aluminum (ν ≈ 0.33): V_L/V_S ≈ 1.96
The Navier Equation
The equation of motion for an isotropic elastic solid combines Newton's second law with the stress-strain relationship:
ρ × ∂²u/∂t² = (λ + μ) × ∇(∇·u) + μ × ∇²u
Where u is the displacement vector field. This equation has two types of solutions:
1. Longitudinal waves: ∇ × u = 0 (irrotational, curl-free)
2. Shear waves: ∇ · u = 0 (solenoidal, divergence-free)
These two solution types propagate independently in unbounded isotropic media. At boundaries and interfaces, they couple through mode conversion - a Level III must understand this coupling to predict beam behavior in complex geometries.
Elastic Constants Relationships - Reference
| Relationship | Formula |
|---|---|
| Young's modulus from Lamé | E = μ(3λ + 2μ)/(λ + μ) |
| Bulk modulus | K = λ + 2μ/3 |
| Poisson's ratio | ν = λ / 2(λ + μ) |
| V_L from E | V_L = √(E(1-ν) / ρ(1+ν)(1-2ν)) |
| V_S from E | V_S = √(E / 2ρ(1+ν)) |
| V_L/V_S ratio | √((2-2ν)/(1-2ν)) |
Velocity Ratio V_L/V_S for Common Materials:
| Material | ν | V_L/V_S | V_L (m/s) | V_S (m/s) |
|---|---|---|---|---|
| Carbon steel | 0.29 | 1.83 | 5,900 | 3,230 |
| Stainless steel | 0.28 | 1.80 | 5,740 | 3,130 |
| Aluminum | 0.33 | 1.96 | 6,320 | 3,130 |
| Titanium | 0.32 | 1.93 | 6,070 | 3,120 |
| Copper | 0.34 | 2.00 | 4,700 | 2,260 |
| Inconel 625 | 0.31 | 1.90 | 5,820 | 3,020 |
| Glass | 0.22 | 1.65 | 5,640 | 3,280 |
| PMMA (Plexiglas) | 0.40 | 2.45 | 2,730 | 1,430 |
Christoffel Equation for Anisotropic Media:
For a plane wave propagating in direction n̂ in an anisotropic solid:
Γᵢₖ × pₖ = ρV² × pᵢ
Where Γᵢₖ = Cᵢⱼₖₗ × nⱼ × nₗ (Christoffel matrix)
The eigenvalues give three possible wave velocities; eigenvectors give polarization directions. In general anisotropic media, all three waves are quasi-longitudinal or quasi-shear - pure modes exist only along symmetry axes.
When Elasticity Theory Matters in Practice
As a Level III, you will encounter situations where isotropic assumptions fail:
Austenitic Weld Metal: Columnar dendritic grain structure creates transverse isotropy. The velocity along the dendrite axis differs from perpendicular directions by 5-15%. Standard beam path calculations using isotropic velocity give incorrect depth readings and beam steering predictions. You must specify technique qualifications that account for this anisotropy.
Rolled Plate Products: Rolling creates crystallographic texture. Carbon steel plates typically show <1% velocity anisotropy between rolling and transverse directions. Titanium and zirconium can show 2-4% anisotropy. For precision thickness measurements on textured materials, specify which direction to measure and which velocity to use.
Single Crystal Components: Gas turbine blades and other single-crystal or directionally solidified components have extreme anisotropy. Velocity can vary by 30%+ depending on propagation direction relative to the crystal axes. Standard UT procedures are inadequate - specialized techniques with crystallographic orientation data are required.
Composite Materials: Fiber-reinforced polymers have velocity ratios (parallel vs perpendicular to fibers) of 2:1 or more. Through-transmission is often the only practical technique because pulse-echo beam behavior is unpredictable in strongly anisotropic layups.
Level III Wave Physics Errors
1. Applying isotropic beam spread formulas to anisotropic weld metal - The standard beam spread formula sin(θ) = 1.22λ/D assumes isotropic velocity. In anisotropic media, the beam spreads asymmetrically and the energy direction (group velocity) differs from the phase velocity direction. Procedures for DMW and austenitic weld examination must account for this.
2. Confusing phase velocity with group velocity - In dispersive media (Lamb waves, guided waves), the velocity at which the wave pattern moves (phase velocity) differs from the velocity at which energy propagates (group velocity). Time-of-flight measurements give group velocity, but refraction angles follow phase velocity. This distinction is critical when developing guided wave procedures.
3. Ignoring frequency-dependent velocity in dispersive modes - Lamb wave velocity depends on the frequency-thickness product. A procedure qualified at one frequency may not work at another because the mode velocity changes, altering beam angles and detection sensitivity.
4. Assuming linear superposition in high-amplitude fields - Near the transducer face and at focal points, acoustic intensity can be high enough that nonlinear effects become significant. This affects harmonic generation and can create artifacts in techniques that rely on frequency analysis.
Evaluating Technique Applicability for Anisotropic Materials
When asked to approve or develop a UT procedure for an anisotropic material, the Level III must systematically evaluate whether the proposed technique can work:
Step 1: Characterize the Anisotropy
- What is the material? (wrought, cast, welded, composite)
- What is the grain structure? (equiaxed, columnar, textured, single crystal)
- What is the degree of anisotropy? (velocity variation with direction)
- Is the anisotropy uniform or position-dependent?
Step 2: Assess Beam Behavior
- Will beam steering occur? (If velocity varies >3% with direction, steering is significant)
- Will beam skewing occur? (Energy deflected out of the examination plane)
- What is the effective beam spread in the anisotropic medium?
- Are there directions where the beam focuses or defocuses?
Step 3: Evaluate Detection Capability
- Signal-to-noise ratio: will grain scattering prevent detection at the required sensitivity?
- Is the beam reaching the intended examination volume?
- Are there shadow zones where the beam cannot penetrate due to steering?
- Can the technique achieve the required sizing accuracy?
Step 4: Select Appropriate Technique
- Mild anisotropy (<3% velocity variation): Standard techniques with frequency optimization
- Moderate anisotropy (3-10%): Refracted longitudinal wave, low frequency, PAUT with beam steering compensation
- Severe anisotropy (>10%): Specialized techniques, possibly TOFD at low frequency, or alternative methods (RT, ET)
Step 5: Qualification Requirements
- Technique must be demonstrated on representative test specimens
- Specimens must replicate the actual anisotropic microstructure
- Detection, sizing, and characterization capability must be proven
- Document all limitations and the valid range of the technique
Diffraction, Scattering, and Advanced Attenuation
Diffraction Phenomena in Ultrasonic Testing
Diffraction is the bending of waves around obstacles and through apertures. In UT, diffraction is fundamental to both flaw detection (tip diffraction signals used in TOFD) and limitations (beam spreading past the geometric shadow boundary).
Huygens' Principle
Every point on a wavefront acts as a source of secondary spherical wavelets. The new wavefront is the envelope of these wavelets. This principle explains:
- Why the beam spreads beyond the near field (wavelets at the edge of the transducer radiate into the geometric shadow)
- Why diffraction signals are generated at flaw tips (the tip acts as a point source of cylindrical waves)
- Why small flaws scatter energy in all directions (the flaw dimension is comparable to the wavelength)
Tip Diffraction - The Physical Basis for TOFD
When a wave encounters a crack tip, the tip acts as a secondary source radiating a diffracted wave in all directions. The diffracted wave amplitude is typically 20-30 dB below the specular reflection from the crack face, but it propagates in all directions including back toward the transmitting transducer.
The TOFD technique exploits tip diffraction by using separate transmit and receive transducers positioned symmetrically about the weld. The time difference between diffracted signals from the upper and lower crack tips gives the through-wall extent:
Δt = (2/c) × √(d² + (s/2)²)
Where d is the depth to the tip, c is the longitudinal velocity, and s is the probe separation.
Kirchhoff Diffraction Theory
The Kirchhoff approximation treats the flaw surface as a combination of reflecting and diffracting regions. For a flat-bottomed hole (FBH):
- The face reflects specularly (strong signal when beam is perpendicular)
- The rim diffracts (weaker signals detectable from multiple angles)
For a crack:
- The crack face reflects specularly (strong signal when beam angle matches crack orientation)
- The crack tips diffract (weaker signals from all angles - used by TOFD)
The ratio of diffracted to reflected amplitude depends on the ka product (k = 2π/λ, a = flaw dimension). When ka >> 1 (flaw much larger than wavelength), diffraction effects are minor. When ka ≈ 1, diffraction is significant and amplitude-based sizing becomes unreliable.
Scattering Regimes and Their Impact on Procedure Design
| Regime | Condition | Scattering Dependence | Impact on UT |
|---|---|---|---|
| Rayleigh | D << λ (D/λ < 0.1) | α_s ∝ D³f⁴ | Minimal; standard procedures work |
| Stochastic | D ≈ λ (0.1 < D/λ < 10) | α_s ∝ Df² | Significant noise; frequency selection critical |
| Diffusion | D >> λ (D/λ > 10) | α_s ∝ 1/D | Beam loses coherence; UT may not be viable |
Where D = grain diameter, λ = wavelength, f = frequency.
Practical Frequency Selection by Material:
| Material | Typical Grain Size | Recommended Frequency | Scattering Regime |
|---|---|---|---|
| Fine-grain carbon steel | 20-50 μm | 2.25-5 MHz | Rayleigh |
| Coarse-grain carbon steel | 100-200 μm | 1.0-2.25 MHz | Rayleigh-Stochastic |
| Austenitic weld metal | 1-5 mm | 0.5-1.0 MHz | Stochastic |
| Cast stainless steel (CCSS) | 5-30 mm | 0.25-1.0 MHz | Stochastic-Diffusion |
| Inconel weld overlay | 0.5-2 mm | 1.0-2.25 MHz | Stochastic |
| Aluminum forgings | 50-500 μm | 2.25-10 MHz | Rayleigh-Stochastic |
| Titanium forgings | 100-500 μm | 2.25-5 MHz | Rayleigh-Stochastic |
Attenuation Coefficient Measurement Methods:
| Method | Applicable To | Accuracy | Notes |
|---|---|---|---|
| Back wall echo decay | Parallel surfaces | ±0.5 dB/inch | Most common field method |
| Through-transmission | Access to both sides | ±0.2 dB/inch | Requires two transducers |
| Spectral analysis | Lab conditions | ±0.1 dB/inch | Frequency-dependent data |
| Buffer rod technique | Lab conditions | High | Eliminates coupling variability |
Procedure: Determining Maximum Examinable Thickness for a Given Material/Frequency Combination
Purpose: Establish the maximum material thickness at which reliable flaw detection can be achieved for a specific material condition and transducer frequency.
Step 1: Measure Material Attenuation
- Select at least 3 representative samples of the material in the condition to be examined (same heat treatment, grain structure)
- Using the intended examination frequency, measure the attenuation coefficient using back wall echo decay or through-transmission
- Record α (dB/mm or dB/inch) at the examination frequency
Step 2: Determine System Dynamic Range
- Measure the maximum usable gain of the instrument (gain at which electronic noise reaches 5% FSH)
- Measure the gain needed to bring the reference reflector to evaluation level at the calibration distance
- System dynamic range = maximum usable gain - reference gain
Step 3: Calculate Maximum Examinable Distance
- Subtract beam spread losses (use DAC curve extrapolation or beam spread formula)
- Remaining dB budget = dynamic range - beam spread loss at maximum distance
- Maximum one-way path: d_max = remaining budget / (2α) (factor of 2 because sound travels there and back)
- For angle beam: convert beam path to thickness using cos(refracted angle)
Step 4: Verify with SNR Measurement
- At the calculated maximum distance, verify that SNR ≥ 6 dB
- If SNR < 6 dB, reduce the maximum examinable thickness until SNR ≥ 6 dB
Step 5: Document
- Record: material type, grain size (if known), frequency, attenuation coefficient, system dynamic range, calculated maximum thickness, verified SNR at maximum thickness
- This becomes part of the procedure qualification record
Case Study: Procedure Qualification Failure - DMW Beam Angle Inadequacy
A nuclear facility developed a UT procedure for examining dissimilar metal welds (DMW) connecting carbon steel nozzles to stainless steel safe ends. The procedure used conventional 45° and 60° shear wave angles at 2.25 MHz - the same technique qualified for carbon steel welds.
Qualification Test: During blind demonstration on a DMW mock-up containing five implanted fatigue cracks, the procedure detected only 2 of 5 flaws.
Root Cause Analysis by the Level III:
1. The DMW weld metal had columnar dendritic grains 8-15mm long oriented perpendicular to the fusion boundary, creating severe acoustic anisotropy (velocity variation of 12% between grain-parallel and grain-perpendicular directions).
2. Shear wave beam steering: The Christoffel equation analysis showed that 45° shear waves entering the columnar structure were steered 10-15° from their geometric path. The beam never reached the predicted flaw positions within the weld volume.
3. Beam skewing: Energy was deflected out of the examination plane, further reducing the signal from in-plane reflectors.
4. Scattering loss: At 2.25 MHz in the columnar structure (D/λ ≈ 5-10), stochastic scattering reduced SNR below 6 dB within the weld volume.
Resolution:
- Replaced shear wave technique with refracted longitudinal wave at 1.0 MHz (longer wavelength moves scattering toward Rayleigh regime)
- Refracted L-waves experience less beam steering in columnar microstructure than shear waves
- Added encoded scanning with beam modeling to compensate for remaining steering
- Revised procedure successfully detected all 5 implanted cracks
- Level III documented the metallurgical basis for technique selection in the procedure qualification record
Level III Lesson: Standard shear wave techniques cannot be blindly applied to acoustically anisotropic materials. The Level III must evaluate the material microstructure, predict beam behavior using anisotropic wave theory, and select technique parameters that maintain detection capability. Procedure qualification on representative specimens is not optional - it is the only way to verify that the technique works in the actual material.
Standards References - Advanced Wave Physics
ASME V Article 14 - Examination System Qualification: Requires demonstration that examination techniques can detect, characterize, and size flaws in the material and geometry of the actual component. Provides framework for performance demonstration testing.
ASME Section XI Appendix VIII - Performance Demonstration: Specifies qualification requirements for UT systems, including personnel, equipment, and procedures. Requires blind testing on specimens containing representative flaws. Directly addresses the need for technique qualification in challenging materials.
ASTM E2192 - Standard Guide for Planar Flaw Height Sizing by Ultrasonics: Provides guidance on sizing techniques including tip diffraction methods. References the theoretical basis for diffraction-based sizing.
ASTM E2375 - Standard Practice for UT of Wrought Products: Addresses material-specific examination requirements including attenuation characterization and frequency selection for different grain structures.
ISO 16810 - NDT - Ultrasonic Testing - General Principles: International framework aligning with the wave physics concepts. References mode conversion, beam behavior, and material property effects.
EPRI Performance Demonstration Initiative (PDI): Industry-specific qualification program for nuclear ISI that directly addresses the challenges of examining dissimilar metal welds and cast stainless steel components.