Flaw Scattering Models

Similar documents
Center for NDE and the Department of Engineering Science and Mechanics and the Engineering Research Institute Iowa State University Ames, Iowa 50011

Signal Loss. A1 A L[Neper] = ln or L[dB] = 20log 1. Proportional loss of signal amplitude with increasing propagation distance: = α d

Gaussian Beams and Transducer Modeling

Physics 451/551 Theoretical Mechanics. G. A. Krafft Old Dominion University Jefferson Lab Lecture 18

Optical Imaging Chapter 5 Light Scattering

Basic concepts to start Mechanics of Materials

Lecture notes 5: Diffraction

Earthscope Imaging Science & CIG Seismology Workshop

Linear Electrical and Acoustic Systems Some Basic Concepts

Surface Waves and Free Oscillations. Surface Waves and Free Oscillations

Attenuation and dispersion

PEAT SEISMOLOGY Lecture 9: Anisotropy, attenuation and anelasticity

ULTRASONIC INSPECTION, MATERIAL NOISE AND. Mehmet Bilgen and James H. Center for NDE Iowa State University Ames, IA 50011

ULTRASONIC SIGNALS FROM "WORST-CASE" HARD-ALPHA INCLUSIONS

MODELLING AND MEASUREMENT OF BACKSCATTERING FROM PARTIALLY WATER-FILLED CYLINDRICAL SHELLS

PLANAR KINETIC EQUATIONS OF MOTION (Section 17.2)

MEASUREMENTS OF ULTRASONIC SCATTERING FROM NEAR-SURFACE FLAWS. D. K. Hsu*, T. A. Gray and R. B. Thompson

Solutions Ph 236b Week 1

SHAPE OF CONVEX VOIDS. James H. Rose* Ames Laboratory Iowa State University Ames, IA Jon L. Opsal**

Wave Propagation in Heterogeneous Media: Born and Rytov Approximations. Chris Sherman

Hydrogeophysics - Seismics

Waveform Parameters. Correlation of Count Rate and Material Properties. Emission Hits and Count. Acoustic Emission Event Energy

Chapter 12. Recall that when a spring is stretched a distance x, it will pull back with a force given by: F = -kx

ME equations. Cylindrical symmetry. Bessel functions 1 kind Bessel functions 2 kind Modifies Bessel functions 1 kind Modifies Bessel functions 2 kind

Handout 6: Rotational motion and moment of inertia. Angular velocity and angular acceleration

EFIT SIMULATIONS FOR ULTRASONIC NDE

6. LIGHT SCATTERING 6.1 The first Born approximation

Solutions: Homework 5

Plane waves and spatial frequency. A plane wave

Steady and unsteady diffusion

On the study of elastic wave scattering and Rayleigh wave velocity measurement of concrete with steel bar

MATH2000 Flux integrals and Gauss divergence theorem (solutions)

E&M. 1 Capacitors. January 2009

BEAM DEFLECTION THE ELASTIC CURVE

Summary of various integrals

Doppler echocardiography & Magnetic Resonance Imaging. Doppler echocardiography. History: - Langevin developed sonar.

AE301 Aerodynamics I UNIT B: Theory of Aerodynamics

Structural Acoustics Applications of the BEM and the FEM

. (70.1) r r. / r. Substituting, we have the following equation for f:

1. Propagation Mechanisms

VI. Local Properties of Radiation

MATH 2083 FINAL EXAM REVIEW The final exam will be on Wednesday, May 4 from 10:00am-12:00pm.

Lamb Wave Behavior in Bridge Girder Geometries

Modelling I. The Need for New Formulas Calculating Near Field, Lateral Resolution and Depth of Field D. Braconnier, E. Carcreff, KJTD, Japan

Lecture 6 Surface Diffusion Driven by Surface Energy

Gauss s Law. Lecture 3. Chapter Course website:

A few principles of classical and quantum mechanics

Anton Nagl, Y. J. Stoyanov*,J. V. Subrahmanyam*, and H. llberall. Department of Physics, Catholic University, Washington, DC 20064

Physics 451/551 Theoretical Mechanics. G. A. Krafft Old Dominion University Jefferson Lab Lecture 18

Lecture 5 Scattering theory, Born Approximation. SS2011: Introduction to Nuclear and Particle Physics, Part 2 2

(a) The points (3, 1, 2) and ( 1, 3, 4) are the endpoints of a diameter of a sphere.

3/22/2016. Chapter 27 Gauss s Law. Chapter 27 Preview. Chapter 27 Preview. Chapter Goal: To understand and apply Gauss s law. Slide 27-2.

MATH 332: Vector Analysis Summer 2005 Homework

Plane waves and spatial frequency. A plane wave

Electricity & Magnetism Qualifier

Advances in Acoustics Metamaterials and Acoustic Cloaks

Topic 7. Electric flux Gauss s Law Divergence of E Application of Gauss Law Curl of E

Part 7. Nonlinearity

Estimating Permeability from Acoustic Velocity and Formation Resistivity Factor

Gauss s Law. Lecture 4. Chapter 27. Channel 61 (clicker) Physics II

Exam 4 Solutions. a. 1,2,and 3 b. 1 and 2, not 3 c. 1 and 3, not 2 d. 2 and 3, not 1 e. only 2

Engineering Mechanics Prof. U. S. Dixit Department of Mechanical Engineering Indian Institute of Technology, Guwahati Kinematics

MATH 280 Multivariate Calculus Fall Integrating a vector field over a curve

Chapter 2 Basic Optics

PHYSICS. Chapter 24 Lecture FOR SCIENTISTS AND ENGINEERS A STRATEGIC APPROACH 4/E RANDALL D. KNIGHT

(1) Recap of Differential Calculus and Integral Calculus (2) Preview of Calculus in three dimensional space (3) Tools for Calculus 3

CBE 6333, R. Levicky 1. Orthogonal Curvilinear Coordinates

Quantum Condensed Matter Physics Lecture 5

PEAT SEISMOLOGY Lecture 12: Earthquake source mechanisms and radiation patterns II

SOLAR MHD Lecture 2 Plan

1.1 Single Variable Calculus versus Multivariable Calculus Rectangular Coordinate Systems... 4

Fourier Optics - Exam #1 Review

If the wavelength is larger than the aperture, the wave will spread out at a large angle. [Picture P445] . Distance l S

Physical principles of Harmonic Imaging Min Joo Choi, PhD

Tracking of Spread Spectrum Signals

SEAFLOOR MAPPING MODELLING UNDERWATER PROPAGATION RAY ACOUSTICS

Lecture 6 Scattering theory Partial Wave Analysis. SS2011: Introduction to Nuclear and Particle Physics, Part 2 2

Oscillatory Motion. Solutions of Selected Problems

Chapter 6. Secondary source theories

Let b be the distance of closest approach between the trajectory of the center of the moving ball and the center of the stationary one.

Some negative results on the use of Helmholtz integral equations for rough-surface scattering

Longitudinal Waves. waves in which the particle or oscillator motion is in the same direction as the wave propagation

Moments of Inertia (7 pages; 23/3/18)

SURFACE WAVE DISPERSION PRACTICAL (Keith Priestley)

Surface Cracks Modelling of Ultrasonic Non-Destructive Testing

Single Particle Motion

EITN90 Radar and Remote Sensing Lecture 5: Target Reflectivity

Development of PC-Based Leak Detection System Using Acoustic Emission Technique

Sound Propagation in the Nocturnal Boundary Layer. Roger Waxler Carrick Talmadge Xiao Di Kenneth Gilbert

Scattering. 1 Classical scattering of a charged particle (Rutherford Scattering)

Cepstral Deconvolution Method for Measurement of Absorption and Scattering Coefficients of Materials

Solution Set of Homework # 2. Friday, September 09, 2017

ENGI Gradient, Divergence, Curl Page 5.01

Topic 4: Waves 4.3 Wave characteristics

14. Rotational Kinematics and Moment of Inertia

Amplitude preserving AMO from true amplitude DMO and inverse DMO

2. Modeling of shrinkage during first drying period

(Refer Slide Time: 1:58 min)

MULTIVARIABLE INTEGRATION

PRACTICE PROBLEMS. Please let me know if you find any mistakes in the text so that i can fix them. 1. Mixed partial derivatives.

Transcription:

Flaw Scattering Models

Learning Objectives Far-field scattering amplitude Kirchhoff approximation Born approximation Separation of Variables Examples of scattering of simple shapes (spherical pore, flat crack, side-drilled hole)

Flaw Scattering Incident plane wave p e i Fluid Model e s y p pressure amplitude flaw x s r s n At many wavelengths away from the flaw the scattered waves are spherical waves p scatt ( y, ω ) = p A( e ; e ) i s exp ( ikr ) r s s scattered pressure A is called the plane wave far-field scattering amplitude

Flaw Scattering Incident plane wave p e i Fluid Model e s y p pressure amplitude flaw x s r s n plane wave far field scattering amplitude of the flaw 1 p% A ik p ik ds 4π n S f p = p ( x s,ω ) ( e ; e ) = + ( e n) % exp( x e ) ( x ) i s s s s s p

Flaw Scattering 1 p% A ik p ik ds 4π n ( e ; e ) = + ( e n) % exp( x e ) ( x ) i s s s s s S f ( ) p = p x s,ω p To obtain the far field scattering amplitude, we need to know the pressure and velocity on the surface of the flaw due to the incident and scattered waves These quantities can be found by solving a flaw scattering boundary value problem.

u scatt ( y ) scattered displacement from a flaw Flaw Scattering ( β P) Elastic Solid ( ) ( β S) A ei ; es A ei ; es, ω = U exp ik r + U exp ik r ( ) p s s s rs rs P-wave S-wave U β ei U displacement amplitude for wave of type β (β = P,S) in problem (1) flaw x s n r s α es Scattered wave of type α (α = P,S)

Flaw Scattering ( β α ; ) i s Ae e Vector scattering amplitude for a scattered wave of type α due to a plane wave of unit displacement amplitude and type β ( β ) ( ; β ; = ) A e e f e e P P P P i s s s ( β ) ; β ( ; β ; = ) A e e f f e e S S S S S i s s s ( α ) ( ) αβ ; αβ ; lkpj l p α α l l k α sk p j s s s 4πρc α x S j ( β = P, S) elastic constants C e u% f = f e = n + ik e n u exp ik ds 2 % x e x displacements and displacement gradients ( α = P,S) ( β = P, S)

Flaw Scattering 3-D scattering amplitude u scatt ( y ) ( β P) ( ) ( β S) A ei ; es A ei ; es, ω = U exp ik r + U exp ik r ( ) p s s s rs rs The received voltage in some measurement models is related to a specific component of the vector scattering amplitude given by ( ω ) = ( e ; ) ( ; )( ) i es = A ei es d A A β α β α α polarization vector of an incident wave from receiving transducer (to be discussed shortly)

Flaw Scattering A(ω) can be obtained from: Analytical/Numerical methods Separation of variables Boundary Elements (BEM) Finite Elements (FEM) Method of Optimal Truncation (MOOT) T-matrix Kirchhoff Approximation Born Approximation or A(ω) can be found experimentally through deconvolution: V R ( ω) = E( ω)a ω ( ) A( ω) = VR E * ( ω) E ( ω) 2 2 ( ω) + n measure model, measure

Flaw Scattering Kirchhoff Approximation for a Volumetric Flaw incident wave e i tangent plane P n e r D x s reflected wave x O S lit on the "lit" surface: fields are those from plane waves interacting with a plane surface on the rest of the surface the fields are assumed to be zero

Flaw Scattering It has been recently shown that the pulse-echo scattering amplitude response, A(ω), for a stress-free flaw (void or crack) in an elastic solid is identical to the scalar (fluid) scattering amplitude: ik A = A = ikβ ds 2π ( ) ( ) β ω e β ; e β ( e β n) exp( 2 x e β ) ( x ) i i i s i s S β ei lit e β s = e β i

Flaw Scattering-Spherical Void For the pulse-echo response of a spherical void of radius a a A( ω) A( ei; ei, ω) = exp( ika) exp( ika) 2 sin ( ka) ka.7.6 2a A(ω) a = 1 mm c = 59m/s.5.4.3.2.1 5 1 15 2 25 3 frequency, MHz

Flaw Scattering-Spherical Void Comparison with Kirchhoff solution with exact separation of variables solution for the P-P scattering amplitude of a void in an elastic solid (pulse-echo) leading edge response 1.4 1.2 exact solution (solid line) 1 "creeping" wave 2 A /a.8.6.4 Kirchhoff solution (dotted line).2 5 1 15 2 25 3 frequency, MHz

Flaw Scattering-Spherical Void Inverting the scattering amplitude (fluid model) into the time domain a 2a c 2a A() t A( ei; ei, t) = δ t U,; t 2 + c 2a c (, ; t) U t t 1 2 1 = t < t < t 1 2 otherwise "lit" surface response A(t) c 4 leading edge response 2a c t

Flaw Scattering-Spherical Void Exact solution for a void in an elastic solid (time domain) vs the Kirchhoff approximation solution 2.5 2 1.5 1.5 creeping wave -.5-1 -1.5-2 2 4 6 8 1 time, μsec

Flaw Scattering-Spherical Void Although the Kirchhoff approximation is formally a high frequency approximation ( ka >>1), numerical tests have shown it capable of producing good agreement (<1 db difference) for the peak-to-peak time domain response of the spherical void down to ka =1 provided the bandwidth is sufficient white: differences <1 db gray: >1 db but < 1.5 db black: > 1.5 db ka 5 4.5 4 3.5 3 2.5 2 1.5 1.5 1 2 3 4 5 6 7 8 9 bandwidth, %

Flaw Scattering-Inclusion Leading Edge Response - General Convex Inclusion (scattered mode same as incident mode) e s tangent plane n α θ i e i x s stat α e s d S lit High frequency, stationary phase approximation θi = α RR 1 2 A( ei; es, ω ) = R12 exp ik i s d 2 e e Plane wave reflection coefficient Gaussian curvature of flaw surface at stationary phase point

Flaw Scattering-Inclusion Leading Edge Response - General Convex Inclusion (scattered mode same as incident mode) e s tangent plane n α θ i x s stat α e s S lit e i In the time domain d + A () t 1 = ( ) exp( ) 2 A ω i ω t d ω π RR A i s t R12 t i s d c 2 ( e ; e, ) = 1 2 δ ( + e e / )

Flaw Scattering-Crack Kirchhoff Approximation Flat Elliptical Cracks e i x 3 e s n e a q r e a 2 1 u x 2 u 2 1 plane normal to e q e q = e i e s e i e s x 1 ( ) A( e i ;e s )= ia 1a 2 e i n e i e s r e J ( 1 k e i e s r ) e r e = a 2 ( 1 e q u ) 2 1 + a 2 ( 2 e q u ) 2 2

Flaw Scattering-Crack General pitch-catch crack response-circular crack 1.8 1.6 1.4 1.2 1.8.6.4.2 5 1 15 2 25 3 frequency, MHz The magnitude of the P-wave pulse-echo far-field scattering amplitude versus frequency for a 1 mm radius circular crack in steel with an angle of incidence of from the crack normal.

Flaw Scattering-Crack Comparison of the Kirchoff approximation and MOOT for a.381 mm radius circular crack at an incident angle of 45 degrees (pulse-echo).25 inc = 45 KIR MOOT.2 amplitude.15.1.5 5 1 15 2 25 frequency, MHz

When e q is parallel to the crack normal, n : A( e i ;e s )= ik ( e i n)a 1 a 2 2 Special case - pulse echo: A( e i ; e i )= ika 1 a 2 2 16 14 12 1 8 6 4 2 Flaw Scattering-Crack 5 1 15 2 25 3 frequency, MHz The magnitude of the P-wave pulse-echo far-field scattering amplitude versus frequency for a 1 mm radius circular crack in steel at normal incidence. e i e i e q n e s n e s

Flaw Scattering-Crack Comparison of the Kirchoff approximation and MOOT for a.381 mm radius crack at normal incidence (pulse-echo) 2 1.8 KIR MOOT amplitude 1.6 1.4 1.2 1 inc= numerical errors in MOOT.8.6.4.2 5 1 15 2 25 frequency, MHz

Inverting these results into the time domain e i n () ( ) aac ( e n) 1 2 i t t < e 2 2 2 i 2 i; s, π ei es re i s e / ( ) At Ae e t = e e r c t otherwise flash points Flaw Scattering-Crack e s re c Example: pulse-echo t

Flaw Scattering-Crack At normal incidence, for pulse-echo A e ; e, t = ( ) i i aa 2c 1 2 dδ dt t ( ) A(e i ; e s ) t

Flaw Scattering-Crack A planar crack is a very specular reflector θ β ei 3 2.5 2 1.5 1.5 1 2 3 4 5 6 7 8 9 angle, degrees

Flaw Scattering-Crack Although the Kirchhoff approximation is formally a high frequency approximation ( ka >>1), numerical tests have shown it capable of producing good agreement (<1 db difference) for the peak-to-peak time domain response of a circular crack at normal incidence down to ka =1.5 provided the bandwidth is sufficient white: differences <1 db gray: >1 db but < 1.5 db black: > 1.5 db ka 5 4.5 4 3.5 3 2.5 2 1.5 1.5 1 2 3 4 5 6 7 8 9 bandwidth, %

Flaw Scattering-Crack Comparison of synthesized waveforms scattered from a.381 mm radius crack using the Kirchhoff approximation and MOOT using frequencies from -25 MHz approximately (pulse echo) 1.9.8.7 amplitude.6.5.4.3.2.1 5 1 15 2 25 frequency, MHz

Flaw Scattering-Crack -15 degrees from crack normal 2 5 KIR MOOT 15 1 1 15 5-5 -1-15 -2

Flaw Scattering-Crack 2-35 degrees from crack normal 4 2 25 KIR MOOT 3 2 3 35 1-1 -2-3 -4-5

Flaw Scattering-Crack 4-55 degrees from crack normal 1.5 4 45 KIR MOOT 1 5 55.5 -.5-1 -1.5

Flaw Scattering-Crack 6-7 degrees from crack normal 1.8 6 65 KIR MOOT.6 7.4.2 -.2 -.4 -.6 -.8-1

Flaw Scattering-Crack 75-85 degrees from crack normal.8.6 75 KIR MOOT 8 85.4.2 -.2 -.4 -.6 -.8

Flaw Scattering-Crack Comparison of peak-to-peak values of MOOT and Kirchhoff versus angle of incidence for the.381mm radius crack (pulse-echo) The arrow shows where agreement is within 1 db 25 MOOT KIR 2 15 1 5 1 2 3 4 5 6 7 8 9

Now consider the scattering amplitude with narrow band Gaussian window The central frequency is 1 MHz, and the bandwidth is 1 MHz Radius of the crack a =.381mm 1.9.8.7 Flaw Scattering-Crack amplitude.6.5.4.3.2.1 5 1 15 2 25 frequency, MHz

Flaw Scattering-Crack the range of angles where the agreement is good is significantly reduced 3.5 KIR MOOT 3 2.5 2 1.5 1.5 1 2 3 4 5 6 7 8 9

Now consider the scattering amplitude with wider band Gaussian window The central frequency is 1 MHz, and the bandwidth is 6 MHz Radius of the crack a =.381mm 1.9.8.7 Flaw Scattering-Crack amplitude.6.5.4.3.2.1 5 1 15 2 25 frequency, MHz

Flaw Scattering-Crack The range of excellent agreement now is back to angles as great as 45 degrees or more 18 16 KIR MOOT 14 12 Amplitude 1 8 6 4 2 1 2 3 4 5 6 7 8 9 Angle, degree

Flaw Scattering-Crack Maximum incident angle where the peak-to-peak time domain agreement between the Kirchhoff approximation and the exact solution for a circular crack is less than 1 db for ka =5. (other ka values shown same trend). max. incident angle, degrees 6 55 5 45 4 35 3 25 2 15 1 1 2 3 4 5 6 7 8 9 bandwidth, %

Flaw Scattering-SDH Kirchhoff approximation for pulse-echo scattering of a side-drilled hole (incident direction in plane perpendicular to the hole axis) β ei L 2b ( β ) kbl β β A( ei ; ei ) = J1( 2kβb) is1( 2kβb) + 2 ( β ) ikbl π Bessel function Struve function

Flaw Scattering-SDH Comparison with exact 2-D separation of variables solution for P-waves (pulse-echo).9.8.7.6.5.4.3.2.1 2 4 6 8 1 The three-dimensional normalized pulse-echo P-wave scattering amplitude versus normalized wave number for a side drilled hole in the Kirchhoff approximation (solid line) and from the exact two-dimensional separation of variables solution (dashed line).

Flaw Scattering-SDH Comparison with exact 2-D separation of variables solution for S-waves (pulse-echo_) 1.1 1.8.6.4.2 2 4 6 8 1 The three-dimensional normalized pulse-echo SV-wave scattering amplitude versus normalized wave number for a side drilled hole in the Kirchhoff approximation (solid line) and from the exact two-dimensional separation of variables solution (dashed line).

Flaw Scattering Kirchhoff approximation - Summary For volumetric flaws -the Kirchhoff approximation properly models the leading edge signal as long as ka >1 approximately but does not model other waves (creeping waves, etc.) For cracks the Kirchhoff approximation models the flash point signals properly in pulse-echo as long as ka > 1 approximately and the incident angle is less than about 5 degrees for wide band responses. For narrow band responses this angle is considerably reduced to as little as 15-2 degrees.

Flaw Scattering -Inclusion The Born Approximation The Born approximation assumes that the material of an inclusion differs little from the host material so that to first order the incident wave passes through the inclusion unchanged (weak scattering, low frequency approximation) e in c incident wave

Flaw Scattering -Inclusion The Born approximation generally is developed from a volume integral expression for the scattering amplitude density difference difference in elastic constants α d u% A u ik e C ik dv 4πρcα x V f j β α q ( ) 2 α m α ei ; es = exp 2 Δ ρω % q + α skδ kqm j ( αx es ) ( x) fields in the flaw replaced by incident fields u% = u% q incident q u% m u% = x x j incident m j

Flaw Scattering -Inclusion For a spherical inclusion in pulse-echo A ( β β e ) i ei ; = 4k b F ( k b) j 2 3 1 2 β β 2kb β spherical Bessel function where F 1 Δρ = + 2 ρ Δc c β relative density difference relative wave speed difference

Flaw Scattering -Inclusion Corresponding time domain impulse response for a spherical inclusion (pulse-echo) front surface response back surface response 2b c α 2b c α t

Comparison with "exact" separation of variables solution for the pulse-echo response of a spherical inclusion by synthesizing a time-domain response from frequencies ranging from -2 MHz approximately for a 1 mm radius inclusion (1 % differences in properties) 1.5 Flaw Scattering -Inclusion 1 amplitude.5 -.5-1 -.5.5 1 time, μsec The time domain pulse-echo P-wave response of a 1 mm radius spherical inclusion in steel where the density and compressional wave speed are both ten percent higher than the host steel. Solid line: Born approximation, dashed line: separation of variables solution.

Flaw Scattering -Inclusion Comparison with "exact" separation of variables solution for the pulse-echo response of a spherical inclusion by synthesizing a time-domain response from frequencies ranging from -2 MHz approximately for a 1 mm radius inclusion (5 % differences in properties) 8 amplitude errors amplitude 6 4 2 time of arrival errors -2-4 -1 -.5.5 1 time, μsec The time domain pulse-echo P-wave response of a 1 mm radius spherical inclusion in steel where the density and compressional wave speed are both fifty percent higher than the host steel.solid line: Born approximation, dashed line: separation of variables solution.

A Doubly Distorted Born Approximation β β ( ) replace wave speeds of host material by that of the flaw in the Born approximation for the pulse-echo response of a spherical inclusion β β ( ) ( kb β ) j 2 3 1 2 β 2 β e ; e = 4k b F 1 Δρ kb F = + 2 ρ Born i i Flaw Scattering -Inclusion ( f β ) j k b 2 3 1 2 i ; i = 4 DDBA fβ 2kf β A e e k b F% b Δc c β where F% 1 Δρ = + ρ Δc 2 f c f β

Front surface amplitude response is improved, and relative time of arrival of back surface response is correct, but there is an error in absolute time of arrivals (5% differences in properties) 8 6 Flaw Scattering -Inclusion amplitude 4 2-2 -4-1 -.5.5 1 time, μsec The time domain pulse-echo P-wave response of a 1 mm radius spherical inclusion in steel where the density and compressional wave speed are both fifty percent higher than the host steel. Solid line: Doubly Distorted Born approximation, dashed line: separation of variables solution.

Flaw Scattering -Inclusion Why is the front surface response amplitude improved by the Doubly Distorted Born Approximation? 2 1.8 1.6 1.4 1.2 F 1.8.6.4.2 R 12 F % 1 1.2 1.4 1.6 1.8 2 2.2 2.4 2.6 2.8 3 Δρ Δc = ρ c ; Answer: because F % R ββ (plane wave reflection coefficient) 12

Flaw Scattering -Inclusion This suggests we apply a phase correction to the doubly distorted Born approximation and replace F % ; by R β β 12 resulting in the modified Born approximation (MBA) for the pulse-echo response of a spherical inclusion β β 2 3 β; β ( ei ; ei ) = 4 2 f 12 exp β 2 f β ( 1 fβ / β ) A k b R ik b c c MD 1 ( 2 f β ) j k b 2k f β b phase correction

Flaw Scattering -Inclusion The MBA (5 % differences in properties) 6 5 4 amplitude 3 2 1-1 -2-1 -.5.5 1 time, μsec The time domain pulse-echo P-wave response of a 1 mm radius spherical inclusion in steel where the density and compressional wave speed are both fifty percent higher than the host steel.solid line: MBA approximation, dashed line: separation of variables solution.

Flaw Scattering -Inclusion The MBA (1 % differences in properties) 2 15 amplitude 1 5-5 -1-1.5-1 -.5.5 1 1.5 2 time, μsec The time domain pulse-echo P-wave response of a 1 mm radius spherical inclusion in steel where the density and compressional wave speed are both one hundred percent higher than the host steel. Solid line: MBA approximation, dashed line: separation of variables solution.

Flaw Scattering The Method of Separation of Variables The sphere and the cylinder are the only two geometries where we can obtain exact separation of variables solutions for elastic wave scattering problems. These are commonly used as "exact" solutions to test more approximate theories and numerical methods.

Flaw Scattering-Void Example 1: pulse-echo P-wave scattering of a spherical void A 1 ; 1 n p p ( e e ) = ( ) A i i n ik p n= A n = E E E E E E E E 3 42 4 32 31 42 41 32 { } 2 2 2 ( 2 1 ) ( /2) s n( pb) 2 p n+ ( p ) ( 2 1) {( 1) n( pb) p n+ ( p )} 2 2 2 () 1 () 1 ( /2) s n ( p ) 2 p n+ ( p ) ( ) () 1 () 1 1 n ( p ) p n+ ( p ) E = n+ n n k b j k + k bj kb 3 1 E = n+ n j k k bj k b 4 1 E = n n k b h k b + k bh k b 31 1 E = n h k b k bh k b 41 1 ( 1) ( 1) () ( ) () ( ) ( ) () ( ) 1 1 E32 = n n+ n hn ksb ksbhn+ 1 ksb 2 2 2 1 () 1 1 s /2 n s s n+ ( kb s ) E = n k b h k b k bh 42 1

Flaw Scattering-Void The normalized P-wave pulse-echo scattering amplitude 2A/b for a spherical void of radius b. 1.4 1.2 2 A / b 1.8.6.4.2 5 1 15 2 25 3 frequency, MHz

Flaw Scattering-Void Using this separation of variables solutions at many frequencies to synthesize a P-wave impulse time domain solution (pulseecho) 6 amplitude -4-8 -12-1 -.5.5 1 1.5 time, μsec The time-domain pulse-echo P-wave response of a.5 mm radius spherical void in steel ( c p = 59 m/s, c s = 32 m/sec) obtained by applying a low-pass cosine-squared windowing filter between 1 and 2 MHz to the separation of variables solution and then inverting the result into the time domain with the inverse Fourier transform.

Flaw Scattering-Void Example 2: pulse-echo SV-wave scattering of a spherical void A 1 = ik s s ( ei; ei ) B n s n= 1 n ( 1) ( 2n+ 1) 2 H J H J J = H H H H H 13 42 43 12 41 13 42 43 12 41 ( 1) ( 1) ( ) ( ) J12 = n n+ n jn ksb ksb jn+ 1 ksb ( 1) ( 1) () ( ) B n () ( ) 1 1 H12 = n n+ n hn ksb ksbhn+ 1 ksb () 1 () 1 ( /2) s n ( p ) 2 p n+ ( p ) 2 2 2 H13 = n n kb h kb + kbh 1 kb J = n j k b k b j k b ( 1) n( s ) s n+ ( s ) ( ) () 1 () 1 1 ( ) ( ) 41 1 H = n h k b k bh k b 41 n s s n+ 1 s ( ) s n( s ) s n+ ( s ) 2 2 2 J42 = n 1 k b /2 j k b + k b j 1 k b 2 2 2 1 1 H 42 n 1 ks b /2 hn ksb + ksbhn+ 1 ksb () = ( ) ( ) ( ) () 1 () 1 1 n ( p ) p n+ ( p ) H = n h k b k bh k b 43 1 () ( )

Flaw Scattering-Void The SV-wave pulse-echo scattering amplitude.7.6.5 amplitude.4.3.2.1 2 4 6 8 1 12 14 16 18 2 frequency, MHz The magnitude of the pulse-echo SV-wave response,, versus frequency for a.5 mm radius spherical void in steel ( c p = 59 m/s, c s = 32 m/sec) as calculated by the method of separation of variables

Flaw Scattering-Void Using this separation of variables solutions at many frequencies to synthesize an SV-wave impulse time domain solution (pulseecho) 6 amplitude -4-8 -12-1 -.5.5 1 1.5 time, μsec The time domain pulse-echo SV-wave response for the same void considered in the P-wave case by applying a low-pass cosinesquared windowing filter between 1 and 2 MHz to the separation of variables solution and then inverting the result into the time domain with the inverse Fourier transform.

Flaw Scattering- SDH Example 3: pulse-echo P-wave scattering of a cylindrical void A p p ( e ; e ) 3D i i L i = n 2π n= ( 2 δ )( 1) n F n A 2D ( ω) 2iπ = kα 2 1/2 A 3D ( ω) L δ n 1 n = = otherwise F n = 1+ ( ) ( ) ( ) ( ) () ( ) () ( ) ( ) ( ) ( ) ( ) () ( ) () ( ) C k b C k b D k b D k b 2 1 2 1 n p n s n p n s C kb C kb D kb D kb 1 1 1 1 n p n s n p n s 2 () ( ) ( ) () 2 ( ) ( ) () ( ) ( ) () ( ) () () ( ( ) + 1 ( )) i i i i C x = n + n k b 2 H x 2nH x xh x n s n n n () () ( ( ) + 1 ( )) i i i i D x = n n+ 1 H x n 2nH x xh x n n n n

Flaw Scattering-SDH Recall, the scattering amplitude for the pulse-echo P-wave case was: A3 D / L.9.8.7.6.5.4.3.2.1 2 4 6 8 1 kb p

Flaw Scattering- SDH Using this separation of variables solutions at many frequencies to synthesize a P-wave impulse time domain solution (pulseecho) 3 amplitude 2 1-1 -2 creeping wave -3-1 -.5.5 1 1.5 time, μsec The time-domain pulse-echo P-wave response of a.5 mm radius cylindrical void in steel (c p = 59 m/s, c s = 32 m/sec) obtained by applying a low-pass cosine-squared windowing filter between 1 and 2 MHz to the separation of variables solution and then inverting the result into the time domain with the inverse Fourier transform.

Flaw Scattering - SDH Example 4: pulse-echo SV-wave scattering of a cylindrical void A sv sv ( e ; e ) 3D i i L δ n i = n 2π n= 1 n = = otherwise ( 2 δ )( 1) n G n G n = 1+ ( ) ( ) 2 ( 1 ) ( ) ( 2 ) ( ) ( 1 ) ( ) () 1 ( ) () 1 ( ) () 1 ( ) () 1 ( ) C k b C k b D k b D k b n s n p n s n p C kb C kb D kb D kb n p n s n p n s

Flaw Scattering -SDH Again, the scattering amplitude for the pulse-echo SVwave case was: 1.1 1.8 A3 D / L.6.4.2 2 4 6 8 1 kb s

Flaw Scattering - SDH Using this separation of variables solutions at many frequencies to synthesize an SV-wave impulse time domain solution (pulse-echo) 3 2 creeping wave amplitude 1-1 -2-3 -4-1 -.5.5 1 1.5 time, μsec The time-domain pulse-echo SV-wave response of a.5 mm radius cylindrical void in steel (c p =59 m/s, c s = 32 m/sec) obtained by applying a low-pass cosine-squared windowing filter between 1 and 2 MHz to the separation of variables solution and then inverting the result into the time domain with the inverse Fourier transform.

Flaw Scattering - SDH Experimentally determined scattering amplitude by deconvolution (side-drilled hole) V R ( ) G( ) ( ω) A ω = ω L G ω = s ω E ω ( ) ( ) ( ) () ( ) ˆ 1 ( ( ) ˆ 2 ) 4πρ2c α 2 E ω = V z, ω V ( z, ω) dz Ta ; ik L α 2Zr

Flaw Scattering -SDH.9.8.7.6 A/L.5.4.3.2.1 2 4 6 8 1 12 14 16 18 2 frequency, MHz