EFFECTS OF ACOUSTIC SCATTERING AT ROUGH SURFACES ON THE

Similar documents
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

SURFACE ROUGHNESS EFFECTS IN POROSITY ASSESSMENT BY ULTRASONIC

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

MEASUREMENT OF REFLECTANCE FUNCTION FOR LAYERED STRUCTURES USING FOCUSED ACOUSTIC WAVES INTRODUCTION

IMPROVEMENT OF TIME REVERSAL PROCESSING IN TITANIUM INSPECTIONS

A COMPARISON OF DIFFERENT METHODS FOR THE DETECTION OF A WEAK ADHESIVE/ADHEREND INTERFACE IN BONDED JOINTS

APPLICATION-DIRECTED MODELING OF RADIATION AND PROPAGATION OF ELASTIC WAVES IN ANISOTROPIC MEDIA: GPSS AND OPOSSM

ULTRASONIC INSPECTION OF TITANIUM ALLOYS WITH A TIME REVERSAL

(2) RC<m) is [1]: Z3 V Z3 Z2 V MEASUREMENT OF CHANGES IN SURFACE ROUGHNESS USING UL TRASONIC REFLECTION COEFFICIENT

A SINGLE TRANSDUCER BROADBAND TECHNIQUE FOR LEAKY LAMB WAVE DETECTION. P. B. Nagy*, W. R. Rose**, and L. Adler

THEORY-ACOUSTIC PRESSURE FIELD AND BEAM DIRECTIVITY

LECTURE NO. 4-5 INTRODUCTION ULTRASONIC * PULSE VELOCITY METHODS

ATTENUATION AND POROSITY ESTIMATION USING THE FREQUENCY-INDEPENDENT PARAMETER Q

Lamb Wave Behavior in Bridge Girder Geometries

Frequency Dependence of Ultrasonic Wave Scattering from Cracks

Shear waves in solid-state materials

CHAPTER 4 BASICS OF ULTRASONIC MEASUREMENT AND ANFIS MODELLING

ULTRASONIC WAVE PROPAGATION IN DISSIMILAR METAL WELDS APPLICATION OF A RAY-BASED MODEL AND COMPARISON WITH EXPERIMENTAL RESULTS

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

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

ULTRASONIC TESTING OF RAILS INCLUDING VERTICAL CRACKS-

DETERMINATION OF ELASTIC CONSTANTS OF ANISOTROPIC MATERIALS FROM OBLIQUE

Optimization of Curved Broadband Arrays for Pipe Inspection

IMPROVED METHODOLOGY FOR INSPECTION RELIABILITY ASSESSMENT FOR DETECTING SYNTHETIC HARD ALPHA INCLUSIONS IN TITANIUM

SIMULATION OF THE INSPECTION OF PLANAR NON MAGNETIC MATERIALS WITH ELECTRO MAGNETIC ACOUSTIC TRANSDUCERS

Angular Spectrum Decomposition Analysis of Second Harmonic Ultrasound Propagation and its Relation to Tissue Harmonic Imaging

SIMULATION OF THE INSPECTION OF PLANAR NON MAGNETIC MATERIALS WITH ELECTRO MAGNETIC ACOUSTIC TRANSDUCERS

v. K. Kinra, Y. Wang and C. Zhu Center for Mechanics of Composites Department of Aerospace Engineering Texas A&M University College Station, TX 77843

J. K. Wang+ and C. S. Tsai

A GEOMETRICAL INTERPRETATION OF THE ECHO FORMATION OF SHORT PULSES ON SIMPLE SHAPED TARGETS. D. De Vadder, R. Raillon and A.

THE ULTRASONIC FIELD OF FOCUSED TRANSDUCERS THROUGH A LIQUID- M. EI Amrani Intercontrole 13 rue du Capricome, Rungis Cedex, France

A model for the ultrasonic field radiated by an Electro-Magnetic Acoustic Transducer in a ferromagnetic solid

Let us consider a typical Michelson interferometer, where a broadband source is used for illumination (Fig. 1a).

Study of the accuracy of ultrasonic flowmeters for liquid

SPECULAR REFLECTION BY CONTACTING CRACK FACES. J.D. Achenbach and A.N. Norris. The Technological Institute Northwestern University Evanston, IL 60201

1532 J. Acoust. Soc. Am. 102 (3), September /97/102(3)/1532/8/$ Acoustical Society of America 1532

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

Lamb Waves in Plate Girder Geometries

ELASTIC MODULI OF SILICON CARBIDE PARTICULATE REINFORCED ALUMINUM METAL MATRIX COMPOSITES

ULTRASONIC ATTENUATION RESULTS OF THERMOPLASTIC RESIN COMPOSITES UNDERGOING THERMAL AND FATIGUE LOADING

Phased Array Inspection at Elevated Temperatures

6th NDT in Progress Evaluation of metal weld-zone by using Poisson s ratio distribution

ULTRASONIC A TTENUA TION RESULTS OF THERMOPLASTIC RESIN COMPOSITES UNDERGOING THERMAL AND FATIGUE LOADING

Numerical and Experimental analysis of long range guided waves for NonDestructiveTesting of pipes.

ULTRASONIC SPECTROSCOPY OF TWO PARALLEL IMPERFECT INTERFACES

Proceedings of Meetings on Acoustics

Physical Characterization of Solid-Liquid Slurries at High Weight Fractions Using Optical and Ultrasonic Methods Grant # /28/02

5. LIGHT MICROSCOPY Abbe s theory of imaging

MECHANICAL RECIPROCITY PRINCIPLES AND ULTRASONIC MEASUREMENT MODELS

Model-Assisted Probability of Detection for Ultrasonic Structural Health Monitoring

EFFICIENT MODELING OF FINITE ACOUSTIC BEAM EXCITATION AND DETECTION OF INTERFACE AND BULK WAVES ON PLANAR AND CYLINDRICAL

Enhanced backscattering due to total internal reflection at a dielectric-air interface

c. P. Hsieh, C-H. Chou, and B. T. Khuri-Yakub

ULTRASONIC REFLECTION BY A PLANAR DISTRIBUTION OF SURFACE BREAKING CRACKS

PROPERTY STUDY ON EMATS WITH VISUALIZATION OF ULTRASONIC PROPAGATION

Pulse-echo method can't measure wave attenuation accurately

Principle and application of ultrasonic wave

Effects of material microstructure and surface geometry on ultrasonic scattering and flaw detection

IR SPECTROSCOPY FOR BONDING SURFACE CONTAMINATION CHARACTERIZATION

Application of Ultra-Sonic Pulse Velocity Techniques for Concrete Quality Measurements

Smart Pigs. Adrian Belanger March 7,

Simulation of Contrast Agent Enhanced Ultrasound Imaging based on Field II

SUBSURFACE WAVES IN SOLIDS WITH CURVED SURFACE AND ACOUSTICAL IMPEDANCE ON IT

Non-contact evaluation of thickness reduction of plates and pipes using EMAT-generated guided wave

a~ a~ a~ a~ DEVELOPMENTS IN THE FINITE ELEMENT MODELING OF ULTRASONIC NDT PHENOMENA R. Ludwig and w. Lord Department of Electrical Engineering

TESTING FOR NONGAUSSIAN FLUCTIJA TIONS IN GRAIN NOISE

ULTRASONIC INVESTIGATION OF THE STIFFNESS OF GRAPHITE-

MODELING THROUGH TRANSMISSION ULTRASONICS FOR CONTROLLING THE FUSION

ACOUSTIC TRANSMISSION WITH MODE CONVERSION PHENOMENON

POD(a) = Pr (Y(a) > '1').

6th NDT in Progress Lamb waves in an anisotropic plate of a single crystal silicon wafer

WAVE PROPAGATION AND SCATTERING IN RANDOM MEDIA

BROADBAND MIMO SONAR SYSTEM: A THEORETICAL AND EXPERIMENTAL APPROACH

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

Negative refraction of ultrasonic waves in 2D phononic crystals.

George A. Alers and Alfonso Manzanares

Electrical Engineering Department Colorado State University Fort Collins, CO 80523

Mathematical modeling of the radiated acoustic field of ultrasonic transducers

10. OPTICAL COHERENCE TOMOGRAPHY

PHYSICS PRACTICAL (CBSE) - X

STUDIES OF ELECTROMAGNETIC SOUND GENERATION FOR NDE* Thomas J. Moran Department of Physics Wayne State University Detroit, Michigan

B 2 P 2, which implies that g B should be

On-stream Inspection for High Temperature Hydrogen Attack

Optimisation using measured Green s function for improving spatial coherence in acoustic measurements

SEAFLOOR MAPPING MODELLING UNDERWATER PROPAGATION RAY ACOUSTICS

ULTRASONIC NDE RHEALOGICAL MEASUREMENT TOOLS FOR INDUSTRIAL PROCESS CONTROL

ANALYSIS OF FREQUENCY CHARACTERISTICS ON NON-INVASIVE ULTRASONIC-DOPPLER FLOW MEASUREMENT FOR METAL PIPES

UL masonic VELOCITY MEASUREMENTS APPLIED TO PROCESS CONTROL

Finite element simulation of the critically refracted longitudinal wave in a solid medium

Measurement of Elastic Constants Using Ultrasound

Electromagnetic Acoustic Transducers for In and Out of plane Ultrasonic Wave Detection

INVERSION ASSUMING WEAK SCATTERING

Today s menu. Last lecture. Ultrasonic measurement systems. What is Ultrasound (cont d...)? What is ultrasound?

Phased Array Rotation Scanner Probe System for Ultrasonic Testing of Sleeve Shafts

Caution! Pay close attention to the special operation and safety instructions in the manual of the ultrasonic echoscope.

A suggested method to be used to measure scattering coefficients of full scale samples.

Single Emitter Detection with Fluorescence and Extinction Spectroscopy

EFIT SIMULATIONS FOR ULTRASONIC NDE

Practical Results of Ultrasonic Imaging by Inverse Wave Field Extrapolation

Transcription:

EFFECTS OF ACOUSTIC SCATTERING AT ROUGH SURFACES ON THE SENSITIVITY OF ULTRASONIC INSPECTION Peter B. Nagy and Laszlo Adler Department of Welding Engineering The Ohio State University Columbus, Ohio 4321 James H. Rose Center for NDE Iowa State University Ames, Iowa 511 INTRODUCTION Ultrasonic inspection of ordinary samples with more or less rough surfaces is an everyday problem in industrial NDE. Contact techniques require flat or other regular (e. g., cylindrical) surfaces of negligible roughness with respect to the acoustic wavelength. Immersion techniques are less susceptible to surface topography, but they still require that the surface radius be larger than the beam diameter and the surface roughness be comparable or less than the wavelength in the immersion fluid. This difference is due to the fact that in immersion inspection surface irregularities do not significantly reduce the energy transmission into the specimen but rather randomize the field through incoherent scattering. Figure 1 shows the schematic diagram of ultrasonic inspection of a rough specimen by the immersion method. The probability of detection of a given flaw is ultimately limited by the signal-to-noise ratio produced at the receiver. The flaw signal results from coherent reflection from a single, relatively large and strong scatterer. In comparison, the noise is incoherent scattering from a large number of randomly distributed, relatively small and weak scatterers such as material inhomogeneities or geometrical irregularities. Surface roughness can substantially reduce the signal-to-noise ratio with respect to an otherwise similar smooth sample. First, surface roughness attenuates the coherent flaw signal much more than the incoherent material noise (1,2]. Second, surface roughness increases the overall noise level by adding another incoherent component to the material noise. This paper discusses the adverse effect of the excess surface noise on ultrasonic flaw detection in rough samples. Review of Progress in Quantitative Nondestructive Evaluation, Vol. 12 Edited by D.O. Thompson and D.E. Chimenti. Plenum Press. New York. 1993 1775

Transducer Rough Surface Flaw Fig. 1 Schematic diagram of ultrasonic inspection of a rough specimen by the immersion method. Under ordinary conditions, material noise is essential unaffected by surface roughness and the degradation of the signal-to-noise ratio is mainly due to the loss of the coherent signal. For all three coherent components, the scattering loss can be written in the following form: 2 2 A refl, shear, long (Binc) h CO K refl, shear, long(binc), (1) where h is the r. m. s. roughness of the surface, co is the angular frequency. K denotes the so-called loss coefficient which can be easily calculated, at least in the phase-screen approximation, for the reflected and shear or longitudinal transmitted waves from the known sound velocities and the angle of incidence, Binc [1]. In some cases, when the material noise is relatively weak and the flaw to be detected is in close proximity of the rough surface, the background noise is dominated by direct backscattering from the rough surface. This effect is demonstrated in Figure 2 showing the r. m. s. noise level in a lowcarbon steel specimen as a function of propagation time. Both the grain size and the surface roughness were approximately 2 11m. The measurement was made by averaging the square of the received rf signal over a -by-2" area in a 3 MHz wide frequency range centered around 1 MHz. There is only weak electrical noise before the arrival of the front reflection. From the smooth side of the specimen, the material noise slightly decreases between the front- and back-wall echoes, which is partly due to the spread of the acoustic beam and partly to scattering induced attenuation in the sample. From the rough side, the overall noise level is much higher. Close to the surface, the additional surface noise is significantly stronger than the inherent material noise. Due to its faster decay, surface noise becomes negligible with respect to the material noise at large depths. Acoustic scattering at a randomly rough liquid-solid interface has been studied for a considerable time. A comprehensive review of the relevant literature was recently published by Ogilvy [3]. The scattering process includes first-, second-, and higher-order effects 1776

1,...,.8 :>....6 -Cl..4... Cl) ril Z.2 front-wall echo back-wall echo 5 1 15 2 Fig. 2 Time [] Example of increased near-surface noise in a low-carbon steel specimen. caused by single-, double-, or higher-order interactions between the acoustic wave and the irregular surface. In the weak-scattering limit, i. e., when the r. m. s. roughness is much smaller than the acoustic wavelength (h «A), the first-order scattering process is much stronger than any of its higher-order counterparts. Still, because of their slower decay, higher-order components become dominant at longer times after the initial interaction. In the following, we shall analyze separately the contributions of first- and higher-order scattering to the additive surface noise. FIRST -ORDER BACKSCATTERING Figure 3 shows the schematic diagrams of first- second-, and third-order scattering from a randomly rough liquid-solid interface. The length of the first-order scattering process, t 1> is limited by the usually rather small variation of the distance between any two points within the aperture of the transmitter on one side and the insonified spot on the surface on the other side: d2 2d tg(9inc) tl --+, 2zv v (2) where d is the diameter of the ultrasonic beam, z is the distance between the transducer and the surface, and v is the sound velocity in the fluid. At normal incidence, the strong but very short first-order scattering is usually overshadowed by the coherent reflection from the surface. At oblique incidence, the first-order scattering is much longer and, at least in the customary pitch-catch arrangement, not affected by the coherent reflection. The problem of first-order acoustic scattering from a randomly rough liquid-solid interface was studied in detail by debilly et ai. by using the potential method [4,5]. Essentially the same results can be obtained by using the much simpler phase-screen approximation originally introduced by Eckhart [6]. The first-order backscattered power WI can be expressed as follows 1777

(a) (b) (c) I \ I \ \ I I \ I \ - [1 [2 [1 [2 [3 [1 Fig. 3 Schematic diagram of (a) first-, (b) second-, and (c) third-order scattering from a randomly rough liquid-solid interface. 2 2 3 4k h cos (Sine) - C [2k sin (Sine )], W;ne At (3) where Wine denotes the incident power, k is the wave number in the immersion fluid, At is the cross-sectional area of the acoustic beam, and C[kr] is the two-dimensional Fourier transform of the correlation function, c(r). In order to assess the additive surface noise, we have to find the best model for the correlation function and determine the r. m. s. roughness and the correlation length of the surface. The roughness is most easily assessed from the surface roughness induced loss of the coherent reflection. As an example, Figure 4 shows the frequency-dependent attenuation of the coherent reflection from a sand-blasted glass specimen at normal incidence. From the best-fitting f2 curve (solid-line), the r. m. s. roughness was determined as h 1 Ilm. The correlation length can be assessed from the angular-dependence of the backscattered signal at a given frequency. Figure 5 shows the normalized backscattering as a function of the angle of incidence for the same sand-blasted glass sample at 12 MHz and d = 1 mm. Apparently, the experimental data can be fit much better by assuming an exponential correlation function rather than the often-used Gaussian one. From the best fitting theoretical curve, the correlation length can be determined as L 1 Ilm. By determining h, L, and the type of the correlation function for the roughness, we have solved the inverse problem of surface characterization. In order to assess the excess surface noise, we have to substitute these data into Eq. 3. As an example, Figure 6 shows the comparison between the measured and predicted backscattered noise from the same glass specimen at Sine = 22. In good agreement with our expectation, at high frequencies, the exponential correlation function fits our measured data much better than the Gaussian one. Considering that the theoretical curve is a very simple approximate prediction for the mean level of the backscattered surface noise without any adjustable parameter, the agreement with the measured data is fairly good. Beside the absolute strength of the surface backscattering, another important factor in the degradation of the signal-to-noise ratio during ultrasonic inspection is the time-depen- 1778

3 o experimental data 25 - best fitting f2 (h=1 1JIll) 2...!::: :g 15 Q) tl < 1 5 5 1 15 2 25 Frequency [MHz] Fig. 4 Scattering loss of the reflected signal at normal incidence from sand-blasted glass...., -3 c:o ' ';i; -35!::: c Q) u til...i<: -4 u -45 CIl c:o ' Q) -5.!::l c; o experimental data -55 -- Gaussian Z - Exponential -6 \ 5 1 15 225 335445 Angle of Incidence [deg] Fig. 5 Normalized backscattering versus angle of incidence at 12 MHz. dence of the scattered signal. For weak scattering, the time-dependence is fully determined by the beam profile, u(f), produced by the ultrasonic transducer. For first-order scattering, (4) where the integration is carried out over the insonified spot on the surface. The timedependence can be numerically calculated from Eq. 4 by substituting r = vt (where the -3 " D experimental - Exponential '";:' -35 -- Gaussian -4 u til...i<: u CIl -45 " N -5-55 Z -6 5 1 15 2 Frequency [MHz] Fig. 6 Normalized backscatter versus frequency from the same sand-blasted glass specimen at Sine = 22. 1779

coordinate vector originates from the center of the transducer's aperture). Figure 7 shows the normalized backscatter as a function of time at six different frequencies along with the calculated result for an ideal rectangular beam profile (constant within the beam diameter and zero eveiywhere else) for d == 1 mm and z == 5 mm. In agreement with our expectation, the length of the backscattered pulse does not depend on frequency and can be easily estimated from the beam diameter by using Eq. 2. Figure 8 shows the shape of the backscattered pulse at different angles of incidence at 8 MHz. Below Bine == 4, the firstorder back scattering is very short and overshadowed by the coherent reflection. At normal incidence, the range of interest for ultrasonic inspection is covered by the lagging "tail" of the backscattered signal, which is weaker but much longer than the first-order initial part. Even for weak scattering (h «A), this slowly decaying part is dominated by higher-order backscattering from the surface. HIGHER-ORDER BACKSCATTERING There seems to be no simple theoretical result available in the literature for the strength of the second- or higher-order backscattering from a randomly rough liquid-solid interface. On the other hand, for weak scattering, the time-dependence of the backscattered signal can be still readily predicted for a given order. For example, the second- and thirdorder components can be expressed as follows: (8) and (9) 1 co 5 1)....5-5 G) () til -1..>:: () co '" -15 IU ' > -2 '" cg -25-3 2 25 3 I:::. /:). 6 MHz 8 MHz 1 MHz 12 MHz 14 MHz 16MHz 35 4-3 co -4... G) () -5 til..>:: () co '" "C IU N -6-7 E -8 Z -9 4,8,12,16,2,24 5 1 15 2 Time [llsl Time [llsl Fig. 7 Relative backscatter versus time at different frequencies. Fig. 8 Normalized back scattering versus time at different angles of incidence. 178

The length of the second-order backscattering, t2 dlv, is mainly limited by the lateral dimension of the beam, while the third-order component is infinitely long even for a finitediameter beam. Figure 9 shows the time-dependence of the second-order back-scattering calculated from Eq. 8 for d = 1 mm and z = 5 mm. In order to demonstrate the effect of edge. ditlfaction on the backscattered signal, results for two other profiles are also shown beside the previously used rectangular one. These other models for the beam profile were derived from the rectangular one by smoothening the distribution by 1 and 2 % of the radius. The principal slope of the decay is 2.6 dbiils and the ditlfaction correction has little effect on this slope although the smoothened profiles predict a slightly longer "exponential" part. Figure 1 shows the measured time-dependence of the backscattered signal from the same sand-blasted glass specimen at normal incidence at seven different frequencies. The backscattered signal greatly increases with frequency, but the its slope remains essentially the same up to approximately 8 MHz. Figure 11 shows that the measured slope agrees very well with the calculated value of 2.6 dbiils at low frequencies, where the weak scattering approximation is acceptable. At higher frequencies, the decay becomes faster as the scattering induced attenuation of the backscattering becomes more significant. Our simple theory predicts that the slope of the second-order backscattering increases with decreasing beam diameter. In a limited range, this behavior was experimentally verified. However, for small beam diameters, the initial "exponential" part became very short and the slope decreased with time indicating a slower than exponential decay. For focused beams, the lateral dimension of the insonified spot on the surface is so small that the second-order scattering becomes very short, too. Figure 12 shows the measured time-dependence of the backscattered signal for a sharply focused ultrasonic beam ofd =.25" and F = 1.25". The backscattering appears to be inversely proportional to the third-power of time. One possible explanation is that third-order scattering dominates the process. For the very small focal spot, r12 r23 vt/2, where t is the time delay after the initial interaction with the surface. In this case, the backscattered intensity aw 3 I at t-3, i. e. Eq. 9 predicts that the backscattered signal decreases proportionally to the third power of time. This prediction is in good agreement with our experimental data shown in Fig. 12. Time Delay [J.IS] Fig. 9 Calculated time-dependence of the second-order backscattering for d = 1 mm. Time Delay [J.IS] Fig. 1 Measured time-dependence of the backscattering at normal incidence. 1781

5-r-----------, o experimental data 4 - theory n "\.J LJ o o o 2 4 6 8 1 12 Frequency [MHz] -3 iii' -4... <I) u -5...: '" u <II -6 CO " -7 <ii -8 Z -9 45 5 55 Time [IJS] 6 9 MHz o 8 o 7 o 6 o 5 6 65 Fig. 11 Slope of the backscattered signal from sand-blasted glass at normal incidence. Fig. 12 Backscattering versus time for a focused beam (d ==.25", F == 1.25"). CONCLUSIONS Direct ultrasonic backscattering from a rough liquid-solid interface produces excess surface noise in addition to the inherent material noise. In weakly scattering materials and near the surface, this excess noise may significantly reduce the signal-to-noise ratio and thereby the probability of detection for a given flaw. At oblique incidence, the surface noise is primarily due to strong first-order scattering. At normal incidence, the first-order scattering becomes very short and the excess noise is primarily due to higher -order scattering. We showed that the crucial time-dependence of different orders can be determined from purely geometrical considerations. ACKNOWLEDGEMENTS This work was sponsored by the N. S. F. under Grant No. ECO-98272. REFERENCES 1. 2. 3. 4. 5. 6. P. B. Nagy, L. Adler, and 1. H. Rose, in Review of Progress in Quantitative Nondestructive Evaluation, (Plenum, New York, 1992), Vol. lib, pp. 1,71-1,78. 1. H. Rose, M. Bilgen, P. B. Nagy, and L. Adler, in Review of Progress in Quantitative Nondestructive Evaluation, (Plenum, New York, 1992), Vol. lib, pp. 1,693-1,7. 1. A. Ogilvy, Theory of Wave Scattering from Random Rough Surfaces (Adam Hilger, Bristol, 1991). M. debilly and G. Quentin, J. Acoust. Soc. Am. 72, 591 (1982). M. deb illy, F. Cohen-Tenoudji, A. Jungman, and G. Quentin, IEEE Trans. Sonics Uitrason. SU-23, 356 (1976). C. Eckhart, 1. Acoust. Soc. Am. 25, 556 (1953). 1782