SUBJ‘IARMGRLC GENEPAHON “‘3 AN ACUUSTIC FixfiiiiY ‘ PERDT ENKRFERO MEYER Thesis far the Degree of M. S. WWGM! STAY}: UNIVERSEW mm A. BAMBI-1R8 1%? f w LIBRARY Michigan State University TH ESYS ABSTRACT SUBHARMONIC GENERATION IN AN ACOUSTIC FABRY-PEROT INTERFEROMETER by Jack A. Bamberg Above a certain energy density threshold in an intense ultra- sonic standing wave, subharmonics of the driving frequency are sdmetimes present. To study this effect, an acoustic Fabry-Perot interferometer was constructed, consisting of two air-backed quartz transducers submerged in water. One of the transducers is driven at frequencies around its 3 MHz resonance. Electronic analysis of the waveforms present in the cavity is accomplished with a spectrum analyser connected to the other quartz transducer. Optical methods are used as a supplementary means of analysis. It was observed that subharmonics exist only at frequencies for which the cavity is in resonance and always occur in pairs, such that the sum of their frequencies is equal to the driving frequency. These discrete frequencies appear only in regions near integral sub- multiples of the driving frequency. The pressure in the cavity at the onset of subharmonic oscillation may be as low as 1.5 atmosphere. In this study, subharmonics were obtained only when the cavity reflector had a frequency dependent reflection coefficient. SUBHARMONIC GENERATION IN AN ACOUSTIC FABRY-PEROT INTERFEROMETER by Jack A. Bamberg A THESIS Submitted to Michigan State University in partial fulfillment of the requirements for the degree of MASTER OF SCIENCE Department of Physics 1967 ACKNOWLEDGMENT The author wishes to eXpress his gratitude to Professor E. A. Hiedemann and Dr. B. D. Cook for the assistance received in the formulation and solution of this problem. The discussions with Dr. W. R. Klein, Dr. F. Ingenito, and the other members of the ultrasonics group at Michigan State University also have been very helpful. The financial support of the Office of Naval Research is gratefully acknowledged. ii TABLE OF CONTENTS ACKNOWLEDGEMENT . . Chapter I. INTRODUCTION . II. APPARATUS . . III. PROCEDURE . . IV. RESULTS . V. DISCUSSION . . APPENDIX: Ultrasonic Light BIBLIOGRAPHY Diffraction 111 Theory Page ii 10 13 29 33 hl LIST OF TABLES Table Page 1.1 Terminology for various frequencies Of the system C O O O O O O O O O O O O O 3 h.1 Typical subharmonic data sample . . . . . . 1h iv Figure 2.1 h.8 h.9 h.10 h.ll LIST OF FIGURES Experimental arrangement of electronic equipment Experimental arrangement of optical equipment Electronic pulse alignment arrangement A typical subharmonic threshold curve Graph of observed subharmonics as a function of frequency, with a 2.0 MHz reflector . . . . . . . . . . . . Graph of subharmonic thresholds as a function of frequency, with a 2.0 MHz reflector . . . . . . . . . . . . . Graph of observed subharmonics; 1.0 MHz reflector . . . . . . . . . Graph of subharmonic thresholds; 1.0 MHz reflector . . . . . . . . . . . . . . Spectrum analyser tracing of subharmonics around f/2 for odd m . . . . . . . . Bands of observed subharmonics in terms of the driving frequency f Brass reflector . . . . . . Typical cavity resonance, showing regions for which subharmonics occur Response curve for cavity, showing relative sound pressure . . . . . . . . Diffraction pattern photographs obtained during subharmonic generation . . . . Page 11 15 16 18 19 20 21 21 22 23 25 27 Figure A.l Graph of average light intensities in the first few diffraction orders, for a stationary wave . . . . . . . vi Page 36 CHAPTER I INTRODUCTION It is sometimes possible to solve problems in acoustics by the elimination of nonlinearities in the equations describing the system; this approximation often leads to valuable information about the system. However, there are also many systems for which linear treatments will not suffice. For example, the deviation from a linear Hooke's Law relationship becomes significant when large amplitude vibrations occur in an elastic medium. A more accurate description may be obtained in such instances, by including a few nonlinear terms as perturbations on the linear situa- tion. In other cases, new phenomena occur in a nonlinear system which can not occur in a linear system. An example of this type is the occur- rence of oscillations at a lower frequency than the resonant frequency of the system, i.e. subharmonics. The above approximate techniques are not useful here, since the effect depends upon the very terms we need to neglect. It is difficult to treat such problems analytically. The existence of subharmonics has long been known; the sub- ject was discussed as early as 1877 by Lord Rayleigh [1]. Since then, numerous examples have appeared in the literature. Most of these exhibit a mechanical nonlinearity leading to an equation of motion of the form m x + c x + f(x) F Out) (1.1) where x is the vibration amplitude, F is the external driving force, m and c are constants of the system, and Solutions of Eq. 1.1 predict subharmonic reSponse under certain conditions. Stoker [2], Hayashi [3], and others have eXplored these solutions in some detail. Recently, nonlinear wave interactions have been shown to produce subharmonic waves at ultrasonic frequencies in liquids. In 1965, Kuljis [h] observed subharmonic generation in water at around 300 kHz in a cavity, using two Spherical bowl transducers to focus the sound. The cavitation threshold was very low, so that subharmonics were present only during the nonlinear cavitation process. Cavitation and focusing, which immensely complicate an analytical analysis, are not required for subharmonic generation. Korpel and Adler [5], using an acoustic Fabry-Perot interferometer Operating around 5 MHz in water, observed subharmonics at intensities well below the cavitation threshold. Analysis of the waveforms present in the cavity was accomplished by observation of the optical diffraction pattern produced by passing monochromatic light through the standing wave sound field. From Eq. A.h in the Appendix, it can be seen that the spacing of diffraction orders is directly proportional to the sound frequency. Hence, the appearance of additional orders between the orders due to the fundamental wave was attributed to the presence of subharmonic waves in the cavity. If the cavity was driven at a multiple of the TABLE 1.1 Terminology for various frequencies of the system. SYMBOL TERM f fundamental (or pump) frequency f fundamental cavity resonant frequency fS subharmonic frequency ftr transducer resonant frequency EXPLANATION The frequency used to excite the cavity. The lowest frequency at which a standing wave can exist in the cavity, i.e. the effective cavity length is equal to x /2, where Kb is the wavelength 8f a wave of frequency fo in the fluid. Frequencies in the cavity lower than f, due to nonlinear coupling to the pump frequency. Resonant frequency of transducer which is driven at the pump frequency. fundamental cavity resonant frequency (see Table 1.1), they observed subharmonics above a certain energy density threshold in the water. The subharmonics occurred in pairs such that the sum of the two fre- quencies was always equal to the pump frequency. Also, the subharmonics were observed to be multiples of f0, and consequently are at resonance in the cavity. The subharmonic frequencies were all in the vicinity of f/Z. In 1966, McCluney [6] and Breazeale [7], utilizing the optical diffraction pattern and electronic means to determine the subharmonic frequencies, concluded that subharmonics occur at integral sub-multiples of the pump frequency. Since a given subharmonic could be obtained at several successive cavity resonances, the subharmonic is very seldom at resonance in the cavity, a result incompatible with the results of Korpel and Adler. Several questions arise from these two studies. Are the subharmonics at a cavity resonance or not? What are the acoustic pressures of the fundamental and subharmonics at threshold? What factors influence the thresholds for the various subharmonics? How does the type of reflector and its thickness effect subharmonic genera- tion? In this study, initiated late in 1965, subharmonic generation in an acoustic Fabry-Perot cavity is systematically investigated, in an attempt to answer these and other questions. CHAPTER II APPARATUS The acoustic Fabry-Perot interferometer is composed of two air-backed, x-cut quartz transducers, submerged in water, and accurately aligned such that the front surfaces are parallel. Each transducer assembly is mounted on an optical bench, which allows the separation between them to be varied from 0.5 cm. to 30 cm. Each transducer I may be easily replaced with other transducers of varying thickness. rs TWO Spring loaded adjustment screws, pivot the reflecting surfaces about two perpendicular axes in the plane of the transducer faces to within 0.05 degree. A stable, variable frequency transmitter with a maximum output power of 100 watts excites one transducer, whose resonant frequency in water is 2975 kHz. The other transducer serves as a reflector and as.a receiver for electronic analysis of waveforms present in the cavity. Both optical and electronic techniques are used to study the waveforms present in the cavity. Electronic analysis is accomp- lished with a Singer Metrics SB 12 spectrum analyser connected to the quartz reflector, as in Fig. 2.1. The analyser displays amplitude as a function of frequency, causing a sinusoidal signal to appear as a vertical "pip" of height proportional to its amplitude. In order to compensate for the non-uniform response of the quartz to various fre- quencies, the reaponse curves as a function of frequency of each quartz transducer to be used as a reflector, was ascertained. A frequency counter measures the transmitter frequency directly; the frequencies of .ucmamasvm ofiaouuomam mo ucoamwamuum Housmawumaxm HzmszDhQ< Hmemm :mHm: MOHummno mo :mmuw . m.: shaman O m H m: m u a H a com 8H 8H oi ofi 62 ow 8 _ _ _ _ _ . . . . u . .I 0mm ..... ”m: M\N MK." m... "Hm m.m l .....n. ...w... ...m ....... .....u an 1 SN im .u. an“. “......“ .... . ......... ”a...” +1 H Am...“ .mmm 1 6mm ........ .... O O.m .I- on" omm _H “\W H a H I I I I 8.38% 926%? Mafia I I 6mm 3sz wUZMDommm M\m m? M\H MW. 8 ... .. 1.1.1.. .. 2. .....n ...“n 1 com ....— H S m.m I Eu O.m "mszmq MHH>v 5:11:11 aaaezmeom Nemaso .1 06m _ .1 co: wozmnomme sz<26mmm NH <60 .1 06m so 01m "mauzma >9H>ummno mo samuu Om ate": a ens omH 66H . 4.: snowflm am 1 N.m TI o.m Assay wozmaommm w.m m.m m\m so o.m macadamon um: H :3. M\ H up N\H wozmpommm azazOmmm m\m ”a "mHUzmA MHH>v _ ueHazmeom Nemflumaou wcflsonm Azufl>mo com pmoswmcmnu Nuuwsv Mom mo>uso omaoamom . OH.: ouswwm Ana: szpomme lam m.m o.m we m.m 1- -L ._ n .7 A _. MHH>HH = 31% p(x,t) . (1.2) where L is the width of the sound beam, and K is the piezo- optic coefficient of the medium. For a single stationary wave in the cavity, the acoustical pressure may be written as * * . * * p(x,t) = p0 [sinfiu t - k=x) + sin.fin t + k x)] , (A.3) where wf is the angular frequency of the sound, and k* is the wave constant of the sound. Solving Eq. A.l, leads to a discrete light distribution with light occurring only at angles n A sin 9n = * : (A'h) A where n is an integer. The time averaged light intensities are given by +oo - 2 2 I..- E Jr (v0) JM (v0) . (15> r=-oo where Jr is the rth order Bessel Function, and v0 = §£k%_§ po . 35 Equation A.h is more general than the above indicates, and holds for any periodic sound field provided the light impinges normally. Figure A.l shows how the first few diffraction orders vary in intensity with pressure [11]. Note that if I; is measured experi- mentally, and the piezo-optic coefficient is known, the acoustic pressure may be determined. When subharmonics are present, Eq. A.3 must be generalized. For the common case of only two subharmonics with f (1) + f (2) = f s s we may write * * * p(x,t) = E pj[sin(mj + - kJ x) + sinQfi: t + kj x)], (A.6) j=O where the subscripts 0, 1, and 2 refer to the fundamental, sub- harmonic fa(1), and subharmonic f (2), respectively. Combining a s Eqs. A.6 and A.2, results in 2 ‘x7_1 v(x,t) = E vj [ sin aj + sin Bj] (A.7) J=O' where 2n LrK J ' A j ’ * * (1 = <0. t — k x , J J J and * * .o>o3 unwcowuwum a you Amuovuo aofluomuMMHo Bow umuww osu OH mofiufimcmucw unwafi owouo>s mo sumac . H.< mesmHh 36 'H 0.0 H.0 N.0 m.o :6 m.o H II—I 37 Using the identity Exp [iha sin C] = E Jr (a) exp[i r O], (A.8) r=-°o where Jr is the rth order Bessel Function, and inserting Eq. A.7 into Eq. A.1, we obtain +00 +00 A(e) = C e100t Z Z Jm(vo) Jfi(vo) Jq(v1) Jr(v1) J8(v2) Ju(v2) m: n: q r: 3:“ =-m =-oo D/2 X 61:1) [i ( 2:11 -D/2 sin 9 + m as + n 80 + qai + r81 + 80% + u82)] dx . (A-9) Evaluating the integral in Eq. A.9, we obtain +00 +00 11(9) = c em}: ZJm(vo) Jn(vo) ngl) Jr(v1) JSIV2) Ju(v2) m) n) q r: 3) u =-00 =-m 2 sin{[§£ sin e + k: (n-m) + k: (r-q) + k; (11-3)] D} X »2fl . 'X- * * ‘ A sin 6 + k0 (n-m) + k1 (r-q) + k2 (u-s). X exp {it [(1): (m+n) + a): (q+r) +5; (s+u)]} . (A.10) 38 1 : V2 = O, and A(e) = 8(9). Using this to evaluate the normalization constant, we obtain With no sound present, v0 = v C = —— . (A.11) If the quantities j = m~- n, k = q - r, .2 = s - u, and 'X' sin[(gflsine-ij-kki-Xk2)D] 3k! — 2n * * * [ A sin 9 j ko k k1 -.2 k2 ] D are defined, Eq. A.10 reduces to +m Me) = E Jm(vo> J.” (vonq (v1) Jq_k (v1) 13 (v2) 1“ (v2) J9m)k:Q: 1 9 s='m ijk‘ exp {ithi—w: (2m- j) +039]? (Zq - k) +u): (28 ‘2) I} . '- (A.12) For the limiting case as D approaches infinity, ij1 1G 0 only if 2n . * * * A sin 6 - Jko - kk1 -‘R k2 — 0 Rewriting Eq. A.13 as sin 9 = -A:- j + -A; k + —;;.2 , "o "1 >‘2 and using the fact that 1 l — —_¥’ + * A A ' l-‘ .nl 39 finally gives the relation sin 6"k' = -L; j' + 5* k' , (A.1h) J KO >"1 where j' = j +32 and k' = k.-.2 . The amplitudes of the discrete orders are Ajkl := 2E:: Jm(vo) Jm-j(vo) Jq(vl) Jq-k(vl) Js(v2) J8-1(v2) m) (I) s =-m X exp {1th. + 1): (2m - 1') + w: (2q-k) + w: (2s -.¢ )1] . (11.15) The time averaged light intensities are .fj-kl = EJ:(VO) “Ti-j (v0) 2 J:(v1) Ji-k(vl) m=-°o qzooo X E 1:072) 1:162) . (A. 16) Equations A.1h and A.16 describe the Fraunhofer diffraction pattern which results from the interference of light with a funda- mental and two subharmonic standing waves which obey Eq. h.1. Equation A.16 is true only under a phase grating assumption. This assumption is not entirely valid for the acoustical pressures observed in the cavity. Therefore, intensities calculated with Eq. A.16 are only approximate. For a detailed discussion of the validity of the phase grating assumption, refer to Klein and Cook [12]. [1] [2] [3] [It] [6] I7] [9] [10] [11] [12) BIBLIOGRAPHY Lord Rayleigh, The Theory of Sound, Volume l, London, Macmillan and CO., 1877, pp.78-85. J. J. Stoker, Nonlinear Vibrations, Interscience Publishers, Inc., N. Y., 1950. C. Hayashi, Nonlinear Oscillations 12_Physical Systems, McGraw- Hill, N. Y., 196A. M. Kuljis, "Experimental Study of Acoustic Parametric Open Cavity Array-Amplifier", Electrotehnika, l, 1965, pp.3-12. A. Korpel and R. Adler, "Parametric Phenomena Observed on Ultrasonic Waves in Water", Appl. Phys. Lett., 1, 15 August 1965, pp. 106-1080 W. R. McCluney, "An Investigation of Subharmonic Generation in an Ultrasonic Resonant Cavity", M.S. Thesis, University of Tennessee, 1966. (Also published as Technical Report No. 2, AD 6hl 60h, Office of Naval Research.) M. A. Breazeale and W. R. McCluney, "Subharmonic Generation in an Ultrasonic Resonant Cavity", J. Acoust. Soc. Am. £0, November, 1966, p. 1262. C. V. Raman and N. S. N. Nath, "The Diffraction of Light by High Frequency Sound Waves: Part II", Proc. Indian Acad. Sci., 2, 1935, h13-h20. W. G. Mayer, "Reflection and Refraction of Mechanical Waves at Solid-Liquid Boundaries. II", J. Appl. Phys. 33, November, 1963, pp. 3286-3288. W. Keck and R. T. Beyer, "Frequency Spectrum of Finite Amplitude Ultrasonic Waves in Liquids", Phys. Fluids, 3, May-June, 1960, pp. 3h6-3h7. B. D. Cook and E. A. Hiedemann, "Diffraction of Light by Ultrasonic Waves of Various Standing Wave Ratios", J. Acoust. Soc. Am., 33, July, 1961, pp. 9h5-9h8. W. R. Klein and B. D. Cook, to be published in IEEE Trans. on Sonics and Ultrasonics. hl IIIIIIIIIIIIIIIIIIIIIIIIIIIIIIII IllIIII][IIIIIIIIIIIIIIIILIIIIIIIIII 93 0308