arxiv:physics/ v1 [physics.flu-dyn] 17 Sep 2003

Transcription

1 A New Damping Mechanism in Non-linear Bubble Dynamics Ahma Moshaii 1,3, asool Saighi-Bonabi 1,2, an Mohamm Taeibi-ahni 4 1 Department of Physics, Sharif University of Technology, P.O. Box: , Tehran, I.. Iran 2 Bonab esearch enter, P.O. Box: , Bonab, Azarbayejan Province, I.. Iran 3 Institute for Stuies in Theoretical Physics an Mathematics, P.O. Box: , Tehran, I.. Iran an 4 Department of Aerospace Engineering, Sharif University of Technology, P.O. Box: , Tehran, I.. Iran arxiv:physics/ v1 [physics.flu-yn] 17 Sep 2003 Non-linear equations of raial motion of a gas bubble in a compressible viscous liqui have been moifie consiering effects of viscosity an compressibility more complete than all previous works. A new set of equations has been erive incluing new terms resulte from consieration of the viscosity an compressibility not only at the bubble interface, but also in the bulk of liqui. The new equations are two non-linear couple equations, which can not be merge into one equation unlike all previously erive equations. Numerical calculations have been performe consiering effects of heat an mass transfer at the bubble interface. The results inicate that the new terms exhibit an important amping role at the collapse, so that their consieration ramatically weakens the bubble rebouns after the collapse. Depenence of this new amping mechanism to amplitue an frequency of the eriving pressure has been investigate. PAS numbers: Bx, Yw, y, Mq I. INTODUTION When a small isolate gas bubble, immerse in a liqui, experiences a high amplitue spherical soun fiel, it grows an contracts non-linearly. Description of the ynamics of such non-linear motion is an ol challenging problem. The complexities of the problem arise from the effects of heat conuction, mass iffusion, compressibility, viscosity, an energy losses involve in ampe oscillations of the bubble. So far, many people have investigate this problem, each concentrating on ifferent aspects of its complexities. However, a rather complete escription has not been presente yet. Lor ayleigh [1] was the first who analytically solve the collapse of an empty bubble in an incompressible liqui. Plesset [2] subsequently stuie the ynamics of a vapor-fille bubble in a time-epenent pressure fiel. Noltingk an Neppiras [3] were the first people erive a ifferential equation for the gas bubble motion uner the influence of such pressure fiel. The extension of this equation to the bubble motion in a compressible liqui has been stuie by many authors; Herring [4], Trilling [5], Gilmore [6], Keller an Koloner [7], Hickling an Plesset [8], Jahsman [9], Flynn [10], Lastman an Wentzell [11], Keller an Miksis [12]. On the other han, heat conuction effects were presente in the works of Hickling [13], Fujikawa an Akumatsu [14], an Yasui [15]. In aition, the works of Löfstet et al. [16], Nigmatulin et al. [17], an Brujan [18] are of the subsequent papers aressing this problem. Between all previous works, the most complete escription of the bubble ynamics equations was presente by Prosperetti an Lezzi [19]. They use a singularperturbation metho of the bubble-wall Mach number Electronic aress: moshaii@mehr.sharif.eu an erive the following one-parameter family of equations escribing the bubble motion in the first orer approximation of compressibility [19]: ( ( 1 (η + 1Ṙ ( t (P l P (3η + 1Ṙ (1 ηṙ Ṙ 2 = (Pl P, (1 where,,, P, an are the bubble raius, liqui soun spee, liqui pressure at infinity, an liqui ensity, respectively. Also, η is an arbitrary parameter. Equation (1 must be supplemente by a bounary conition equation at the bubble interface to relate the liqui pressure, P l, to the gas pressure insie the bubble. Like all previous authors, Prosperetti an Lezzi [19] use the following incompressible equation for this purpose: P l = P g 4µṘ 2σ, (2 where, P g, µ, an σ are the gas pressure at the bubble interface, liqui viscosity, an surface tension, respectively. Most of the previously erive equations belong to this single parameter family of equations, corresponing to ifferent values of η. Moreover, η = 0 yiels results in closest agreement with the numerical simulation of full partial ifferential equations [19]. Two specific approximations have been use in the erivation of Eq ns. (1 an (2. The first approximation is the negligence of the viscosity effects for the liqui motion aroun the bubble, which has been use in the erivation of Eq n. (1. In fact, Eq n. (1 has been erive from the Euler equation, in which the viscosity is eliminate. Note that the viscous term of Eq n. (2 has been resulte from the liqui viscosity at the bubble interface, but not from the bulk of liqui.

2 2 The secon approximation is the incompressibility assumption of the liqui an the gas at the bubble interface, which has been use in the erivation of Eq n. (2. All of the effects of the liqui compressibility in the work of Prosperetti an Lezzi, as well as in all other previous works, have been resulte from the liqui motion aroun the bubble, but not from the bubble bounary conition equation. In fact, all previous authors, on one han took into account the compressibility of the liqui motion aroun the bubble, but on the other han neglecte its consieration at the bubble interface. Although, the two mentione approximations have been use in the erivations of the existing bubble ynamics equations [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19], but the applicability of these approximations for all times of the bubble motion nees to be clarifie. Especially, at the en of the collapse, when the bubble motion is significantly compressible, these approximations may be inapplicable. In this paper, by eliminating the above mentione approximations, we have originally moifie the bubble ynamics equations. A new set of equations have been erive incluing all effects of the viscosity an compressibility of both the liqui an the gas. These equations contain new terms resulte from the effects of two coefficients of viscosity of both the liqui an the gas. The influence of the ae new terms arisen from the liqui viscosity has been numerically investigate. The results clearly inicate that the aition of the new terms consierably affects the bubble motion at the collapse time an uring the bubble rebouns. II. DEIVATION OF THE BUBBLE DYNAMIS EQUATIONS To erive the bubble ynamics equations, we assume that the motion of the bubble interface an the infinitely extene surrouning liqui are always spherically symmetric. Uner this circumstance, the continuity an the momentum equations are as follows [20]: [ ] u t + u u [ 1 t + u ] = u 2u r = p + µ [ + (λ + µ [ 1 r 2 ( r 2 u 2u r 2 =, (3 ( 1 ( r 2 r 2 u ] ], (4 where,, u, p, an are ensity, velocity, pressure, an ivergence of the velocity, respectively. Also, µ an λ are first an secon coefficients of viscosity. The two brackets of the right han sie of Eq n. (4 have the same ifferential forms. Hence: ( u t + u u = p + (λ + 2µ. (5 Diviing Eq n. (5 by an integrating it with respect to r from the bubble interface to infinity at a fixe instant of time an assuming that the two coefficients of viscosity are constant, we get: u t Ṙ2 r 2 = 1 1 p + (λ + 2µ, (6 where, the liqui velocity at the far fiel is assume to be sufficiently small. When the liqui ensity oes not strongly change, Eq n. (6 can be written as: u Ṙ2 r = t 2 P P l ( l. (7 where, l an are the ivergence of the liqui velocity at the bubble interface an at infinity, respectively. The liqui pressure at the far fiel, P, is: P (t = P 0 + P a (t, (8 where, P 0 an P a (t are ambient an riving pressures, respectively. From Eq n. (3 the ivergence of the velocity can be written as: = 1 t = 1 p c 2 t, (9 where the soun spee, c, which is assume to be constant in the liqui, is efine as: c 2 = p/. Thus, Eq n. (7 becomes: u Ṙ2 r = t 2 P P l 2 2 t (P l P. (10 The viscous term in Eq n. (10 has been resulte from the simultaneous effects of the liqui viscosity an compressibility. This term has been neglecte in the erive equations of the previous authors. To convert Eq n. (10 to a ifferential equation, the integral term must be approximate. Because of the irrotationality assumption of the bubble motion the velocity potential, ϕ(r, t, can be introuce into this equation as: ϕ t (, t = Ṙ2 2 + P l P 2 2 t (P l P, (11 where subscript t enotes the temporal erivation. In the simplest approximation of ϕ t (, t, the liqui motion is assume to be incompressible. So, the liqui velocity at any istance r will be: u(r = 2 Ṙ r 2. (12

3 3 Accoring to this approximation, ϕ t (, t can be written as: [ ( ] 2 Ṙ ϕ t (, t = = t r + 2Ṙ2. (13 r= Inserting Eq n. (13 into Eq n. (11 an eliminating the liqui viscous term ue to incompressibility approximation yiels the well-known ayleigh-plesset equation: + 3 2Ṙ2 = P l P. (14 To introuce the compressibility effects into approximation of ϕ t (, t, we assume that the liqui potential fulfills the following acoustic equation for the spherical waves: t (rϕ + (rϕ = 0. (15 As we have shown in the following, this approximation is equivalent to the first orer compressibility consieration of ef. [19]. Double raial ifferentiation of Eq n. (15 results: t (r + (r = 0. (16 Through the efinition of total erivative of the velocity along with Eq n. (3 it can be obtaine that: u ( u ϕ t (, t= t r = 2u2 u + r. (17 t r By a partial integration: ϕ t (, t = + 2Ṙ2 Ṙ l + Applying Eq n. (16 in Eq n. (18 yiels: r r. (18 t ϕ t (, t = + 2Ṙ2 Ṙ l + l. (19 Substituting Eq n. (18 into Eq n. (11 results: ( + 2Ṙ2 3 = 1 Ṙ P l t + P l P 2 2 t (P l P. (20 This equation is the moifie form of Flynn equation [10] along with viscosity consieration. The compressibility effects can be introuce in a ifferent manner. Differentiating Eq n. (15 with respect to t an substituting ϕ t (, t in it from Eq n. (11 yiels: Ṙ ( u + ( p t t ( u Ṙ + P P l ( p λ + 2µ 2 2 P t + + 2Ṙ2 1 = t (P l P, (21 in which, the spatial an the temporal partial erivatives of the viscous term have been neglecte ue to smallness as well as the spatial erivative of the liqui pressure at infinity. Inserting the relation = p/c 2 into Eq n. (3 an using this equation an the momentum equation along with the total erivatives of the pressure an the velocity at the bubble interface, the four partial erivatives in Eq n. (21 can be obtaine as: ( u = t + Ṙ P l 2 t + 2Ṙ2, (22 ( p = Ṙ t + P l t, (23 ( u = 1 P l 2 t 2Ṙ, (24 ( p =. (25 Note, in the erivation of the Eq ns. (22-25, the spatial an the temporal partial erivatives of the two viscous terms in the Navier-Stokes equation, Eq n. (4, have been neglecte. Inserting Eq ns. (22-25 into Eq n. (21 an retaining only the terms up to the orer of Ṙ/ results the moifie form of Herring-Trilling equation [4, 5] along with viscosity effects: ( ( 1 2Ṙ + 1 4Ṙ 3 + P l P 3 2Ṙ2 = 2 2 t (P l P t (P l P. (26 It is straight forwar to obtain the moifie form of Eq n. (1 by the aition of Eq n. (26 an the prouct of (1 ηṙ/ with Eq n. (14 as: ( ( 1 (η + 1Ṙ + 1 (3η + 1 Ṙ 3 2Ṙ2 = 3 ( t (P P l P l P (1 ηṙ 2 2 t (P l P. (27 The new Eq n. (27 provies a suitable escription of the simultaneous effects of liqui compressibility an liqui viscosity in the bubble motion. The ae viscous new term of this equation arises from the liqui compressibility.

4 4 To complete the argument, it is necessary to moify the bounary conition Eq n. (2 with the compressibility effects. The raial component of the stress tensor is: ( u T rr = p + λ u + 2µ. (28 Inserting the velocity ivergence from Eq n. (3, into this equation yiels: ( u T rr = p + (λ + 2µ + 2u 4 µu r r = p + (λ + 2µ 4 µu r. (29 The bounary continuity requirement at the bubble interface is: T rr (liqui = T rr (gas +2 σ. (30 Substituting Eq ns. (29 an (9 into Eq n. (30 leas to: P l + 4 µṙ ( λ + 2µ + P l 2 = P g + 4 µ gṙ t ( λg + 2µ g g + 2 σ t, (31 g where, g, µ g, an λ g are the gas ensity, the first an the secon coefficients of viscosity of the gas, respectively. The new Eq n (31 provies the most complete bounary conition equation at the bubble interface, which contains all effects of compressibility an viscosity of both the liqui an the gas. Here, we concentrate on the effects of liqui viscous new term of Eq n (31 an eliminate the gas viscosity terms as previous works [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19]. Therefore, Eq n. (31 becomes: ( λ + 2µ P l P l + 2 = P g 4 µṙ t 2 σ. (32 omparison of Eq ns. (2 an (32 inicates the existence of a new viscous term in Eq n. (32 ue to the liqui compressibility. The set of Eq ns. (27 an (32 present the most complete form of the bubble ynamics equations containing effects of the liqui compressibility an viscosity not only at the bubble interface, but also in the bulk of liqui. In fact, these equations account for the viscosity of a compressible liqui. While, all previous equations accounte for the viscosity of an incompressible liqui an compressibility, separately. To generalize the argument, Eq ns. (27 an (32 are expresse in imensionless forms. The imensionless variables of this problem are efine as: = 0, = Ṙ, P g = P g 2, t = t 0, P l = P l 2, P = P 2, (33 where, 0 is ambient raius of the bubble. Substituting the imensionless variables into Eq ns. (27 an (32, the imensionless equations are obtaine as: ( 1 (η + 1Ṙ + 3 ( (3η + 1Ṙ Ṙ 2 ( = 1 + (1 ηṙ (Pl Pa P0 + t (P l P a, (34 Pl = Pg Ṙ 4µ 2 σ (λ + 2µ P l t. (35 The quantities σ, µ, an λ are imensionless surface tension an imensionless liqui viscosity coefficients, which are efine as: σ = σ/ 0 2, µ = µ/ 0, an λ = λ/ 0. These imensionless numbers, which are basically inverse of Weber Number an inverse of eynols Number, characterize significance of the surface tension an the liqui viscosity in the bubble ynamics. III. BUBBLE INTEIO EVOLUTION To quantify effects of the new viscous terms on the bubble ynamics, evolution of the gas pressure at the bubble interface, P g, must be specifie. It can be etermine from simultaneous solution of the conservation equations for the bubble interior an the bubble raius equations [23, 24, 25, 26, 27, 28]. Also, heat conuction an mass exchange between the bubble an the surrouning liqui affect the bubble evolution. These complexities were consiere in a complete gas ynamics simulation by Storey an Szeri [29]. On the other han, strong spatial inhomogeneities insie the bubble are not remarkably reveale, unless at the en of an intense collapse [26, 27]. Therefore, the uniformity assumption for the bubble interior seems to be useful an provies many features of the bubble motion [21, 22]. Using this assumption Toegel et al. presente an ODE moel [30], in which effects of heat an mass transfer at the bubble interface have been consiere. This moel accurately escribes various experimental phase iagrams [31] an provies a goo agreement with the complete irect numerical simulation of Storey an Szeri [29]. In this paper, we use the Toegel et al. s moel [30] for specifying the bubble interior evolution. We escribe an argon bubble in water uner the conitions of Single Bubble Sonoluminescence [21, 22]. The gas evolution moel can be summarize as follows. The gas pressure is moele by the van er Waals equation of state: P g = N totkt V N tot B (36

5 5 where, N tot = N Ar + N H2O is total number of particles insie the bubble. The covolume constant B = m 3 is assume to be equal for both water vapor an argon [31]. The value of N tot changes with time because of evaporation an conensation of the water vapor molecules at the bubble interface. The rate of change can be moele as [30]: Ṅ tot =4π 2 D n 0 H2 O n H2 O l, l = min ( D Ṙ,, π (37 where, n an n H2 O 0 H2 O are the instantaneous an equilibrium concentration of the water vapor molecules, respectively. The iffusion coefficient, D, is given by [32]: D = D 0 (n 0 /n tot, where D 0 = m 2 /s an n 0 = m 3. The quantity l is thickness of iffusive bounary layer (see ef. [30] for more etails. The equilibrium concentration of the water vapor molecules is given by the number ensity of the saturate vapor pressure at the ambient liqui temperature T 0 ; n 0H2 O = P ν (T 0 /kt 0. In our calculations: n 0H2 O = m 3. Similar to Eq n. (37, the heat exchange at the bubble interface can be approximate by [30]: ( Q = 4π 2 κ T 0 T χ, l th = min l th Ṙ,, (38 π where, T an T 0 are the gas temperature an the ambient liqui temperature, respectively. Also, κ = W/mK [32] is thermal conuctivity coefficient of the gas content, l th is thickness of thermal bounary layer an χ is thermal iffusivity coefficient, which is given by: χ = κ/c p, with c p = ( 5 2 n Ar n H 2Ok the constant pressure heat capacity per unit volume of the gas. Applying the basic energy equation for the bubble content evolution results the rate of change of the gas temperature [30, 31]: T = Q P g V [ + 4T 0 3T T ( ] θ i /T k Ṅ tot ν ν e θi/t 1 ν (39 ν = 3 2 N Ark + [3 + ( (θ i /T 2 e θi/t ] kn H2O, (40 (e θi/t 1 2 where, three ifferent values for θ i correspon to characteristic vibration temperatures of H 2 O: θ 1 = 2295K, θ 2 = 5255K, an θ 3 = 5400K [31]. In this paper for simplicity we neglecte the effects of chemical reactions [31], which are important only at an extremely short time interval when the bubble temperature is high enough (more than 5000K. Equations (36-40 along with the bubble ynamics equations are the set of equations, which totally escribe the evolution of the bubble characteristics. Uner FIG. 1: Time variations of the bubble raius, accoring to the new (soli an the ol (ashe sets of equations. Graphs (a an (b shows the bubble evolution in a complete perio an uring the bubble rebouns, respectively. The equilibrium raius is 0 = 4.0 µm an the eriving pressure amplitue is P a = 1.4 atm. these circumstances, time variations of the bubble properties have been numerically calculate for both the new (Eq ns. 34,35 an the ol (Eq ns. 1,2 bubble ynamics equations (for η = 0. IV. NUMEIAL ANALYSIS The calculations were carrie out for a perioic riving pressure; P a (t = P a sin (ωt, with ω = 2π 26.5 khz. The constants an parameters of the bubble ynamics equations were set for the water at room temperature, T 0 = K, an atmospheric ambient pressure, P 0 = 1.0 atm; = kg/m 3, = m/s, µ = kg/ms, σ = kgs 2 [34]. The secon coefficient of viscosity of water at room temperature was set to be λ = kg/ms [35]. Figures (1-3 illustrate the results of our calculations for P a = 1.4 atm an 0 = 4.0 µm. Similar values for these parameters have been reporte in recent experi-

6 6 FIG. 2: Time variations of the number of H 2O particles insie the bubble accoring to the new (soli an the ol (ashe sets of equations. Graphs (a an (b show the evolution in a complete perio an uring the bubble rebouns, respectively. The parameters an constants are the same as Fig. (1. mental works of Ketterling an Apfel [36], Simon et al. [37], an Vazquez et al. [38]. Figure (1 shows the variations of the bubble raius for the new an the ol bubble ynamics equations. It is observe that the aition of the new viscous terms consierably changes the bubble evolution after the collapse. The bubble motion is remarkably compressible uring the collapse. Therefore, the new viscous terms, which have been arisen from the liqui compressibility, are important in this time interval. These terms exhibit a amping role an their consieration reuces the amplitue of the bubble rebouns. Also, the perio of the rebouns ecreases with the aition of the new terms. Details of our calculations at the en of the collapse show that the minimum raius for the new case is about 4% greater than that of the ol one. The ifference between the two cases for the evolution of the number of H 2 O particles has been illustrate in Fig. (2. As this figure inicates, a large amount of water vapor evaporates into the bubble uring the expansion. Inee, at maximum raius, most of the bubble content is water vapor. During the collapse, the water FIG. 3: The gas temperature evolution uring the bubble rebouns (a an at the en of the collapse (b for the the new (soli an the ol (ashe cases for the same parameters an constants as Fig. (1. vapor molecules rapily conense to the liqui. It is seen that the remarkable ifference between the two cases also appears on the variations of the number of H 2 O particles after the collapse. Note that, the ifference graually isappears as the bubble rebouns weaken. In Fig. (3, the gas temperature evolution uring the bubble rebouns an at the en of the collapse have been emonstrate. Damping feature of the new term is clearly observe by consierable ecrease of the peak temperature at the collapse time as well as remarkable reuce of the seconary peak temperatures. Also, the pulse with of the main peak temperature ecreases with the aition of the new terms. Figures (4 an (5 illustrate epenence of the bubble properties at the en of the collapse to the variation of the amplitue an the frequency of riving pressure for the two cases. The ifferent values of P a, corresponing to a specific values of 0, can be experimentally obtaine by varying the concentration of the issolve gas in the liqui [21, 36, 37, 38]. Also, the frequency range in the calculations of Fig. (5 is in the range of the experimental results [21]. Figure (4 shows the variation of the peak tempera-

7 7 FIG. 4: The bubble characteristics at the time of collapse as a function of the riving pressure amplitue for the new (soli an the ol (ashe bubble ynamics equations; peak temperature (a, peak pressure (b, an minimum raius (c. The equilibrium raius was fixe ( 0 = 4.5 µm. Other constants are the same as Figs. (1-3. FIG. 5: The epenence of the bubble characteristics at the en of the collapse to the frequency of riving pressure for the new (soli an ol (ashe bubble ynamics equations; peak temperature (a, peak pressure (b, an minimum raius (c. The equilibrium raius an the pressure amplitue were fixe ( 0 = 4.5 µm an P a = 1.35 atm. Other constants are the same as Figs. (1-3. ture, the peak pressure, an the minimum raius as a function of the riving pressure amplitue. The ambient raius was fixe ( 0 = 4.5 µm. The amping feature of the new viscous terms is clearly observe in this figure. For the peak temperature in Fig. 4(a, the maximum relative ifference appears aroun P a = 1.25 atm (about 35%. The ifference between the two cases ecreases for higher riving pressures. However, for the peak pressure in Fig. 4(b, the ifference increases with the increase of the amplitue. The values of min for the two cases are similar in high riving pressures because of the effects of the exclue volume, which prevent the bubble compression [33].

8 8 Figure (5 represents the epenence of the bubble characteristics at the en of the collapse to the variations of the riving pressure frequency. The ambient raius an the amplitue were fixe in these calculations ( 0 = 4.5 µm an P a = 1.35 atm. As Fig. (4, the amping feature of the new terms is also seen here. The ifference between the two cases for the peak temperature an the minimum raius increases when the frequency is increase. While, for the peak pressure the ifference reuces with the increase of the frequency. A major eficiency of the ol bubble ynamics equations is that for strongly riven bubbles, such as sonoluminescence bubbles, large amplitue rebouns are prouce after the collapse, so that they often last until the next acoustic cycle of the perioic riving pressure. This is in contrast with the experimental results, which show rapily ampe rebouns [21, 39]. By introucing a amping term arisen from the gas compressibility, Moss et. al [39] provie a typical solution for this problem. The effects of the suggeste term by Moss et. al is very similar to the amping effects of the new terms in this paper, (compare Fig. 1(b with Figs. (3 an (4 of ef. [39]. It seems the amping feature of the the new terms is a better way for solving the mentione problem. The reason is that Eq ns. (34 an (35 have been erive irectly from the basic equations of flui mechanics, on the contrary to Eq n. (3.2 of ef. [39], which was erive by an approximate metho. V. ONLUSIONS The equations of the bubble motion in a compressible viscous liqui were newly erive from the full Navier- Stokes equations. These equations contain two similar new terms resulte from the simultaneous effects of the liqui viscosity an compressibility. These new terms have a consierable amping role at the collapse, when the bubble motion is significantly compressible. This new amping mechanism ramatically changes the bubble properties at the en of the collapse an uring the bubble rebouns. The results of this work inicate that, the neglect of the new terms in the previous works is not reasonable for the collapse time an the new effects shoul be consiere for the preiction of the quantities relate to the collapse, e.g., the value of light emission by a single sonoluminescing bubble as well as the bubble stability limits. AKNOWLEDGEMENTS This work was supporte by Sharif University of Technology an Bonab esearch enter. Partial support of this work by Institute for Stuies in Theoretical Physics an Mathematics is appreciate. The authors thank Anrea Prosperetti for his helpful comments. [1] L. ayleigh, On the pressure evelope in a liqui uring the collapse of a spherical cavity. Philos. Mag., 34, 94 (1917. [2] M. S. Plesset, The ynamics of cavitaion bubbles. J. Appl. Mech. 16, 277 (1949. [3] B. E. Noltingk an E. A. Neppiras, avitation prouce by ultrasonics. Proc. Phys. Soc. Lonon B63, 674 (1950. [4]. Herring, Theory of the pulsations of the gas bubble prouce by an unerwater explosion. OSD ep. No. 236 (ND 4-sr (1941. [5] L. Trilling, The collapse an reboun of a gas bubble. J. Appl. Phys. 23, 14 (1952. [6] F.. Gilmore, The growth or collapse of a spherical bubble in a viscous compressible liqui. ep. No. 26-4, Hyroyn. Lab., alif. Inst. Tech. (1952. [7] J. B. Keller an I. I. Koloner, Damping of unerwater explosion bubble oscillations. J. Appl. Phys. 27, 1152 (1956. [8]. Hickling an M. S. Plesset, ollapse an reboun of a spherical bubble in water. Phys. Fluis 7, 7 (1964. [9] W. E. Jahsman, ollapse of a gas-fille spherical cavity. Trans. ASME E: J. Appl. Mech. 35, 579 (1968. [10] H. G. Flynn, avitation ynamics. I. A mathematical formulation. J. Acoust. Soc. Am. 57, 1379 (1975. [11] G. J. Lastman an. A. Wentzell, avitation of a bubble in an invisci compressible liqui. Phys. Fluis 22, 2259 (1979. [12] J. B. Keller an M. Miksis, Bubble oscillations of large amplitue. J. Acoust. Soc. Am. 68, 628 (1980. [13]. Hickling, Effects of thermal conuction in sonoluminescence. J. Acoust. Soc. Am. 35, 967 (1963. [14] S. Fujikawa an T. Akamatsu, Effects of the nonequlibrium conensation of vapour on the pressure wave prouce by the collapse of a bubble in a liqui. J. Flui Mech. 97, 481 (1980. [15] K. Yasui, Alternative moel of single-bubble sonoluminescence. Phys. ev. E 56, 6750 (1997. [16]. Löfstet, B. P. Barber, an S. J. Putterman, Towar a hyroynamic theory of sonoluminescence. Phys. Flui A 5, 2911 (1993. [17]. I. Nigmatulin, I. SH. Akhatov, N. K. Vakhitova, an. T. Lahey, On the force oscillations of a small gas bubble in a spherical liqui-fille flask. J. Flui Mech. 414, 47 (2000. [18] E. A. Brujan, The equation of bubble ynamics in a compressible linear viscoelastic liqui. Flui Dynamics esearch 29, 287 (2001. [19] A. Prosperetti an A. Lezzi, Bubble ynamics in a compressible liqui. 1. First orer theory. J. Flui Mech. 168, 457 (1986; A Lezzi an A. Prosperetti, Bubble ynamics in a compressible liqui. 2. Secon orer theory. J. Flui Mech. 185, 289 (1987. [20]. Aris, Vectors, tensors, an the basic equations of flui mechanics, Prentice-Hall, IN., Englewoo liffs, (1962. [21] B. P. Barber,. A. Hiller,. Löfstet, S. J. Putterman, an K.. Weninger, Defining the unknowns of sonoluminescence. Phys. ep. 281, 65 (1997.

9 9 [22] M. P. Brenner, S. Hilgenfelt, an D. Lohse, Singlebubble sonoluminescence. ev. Mo. Phys.74, 425 (2002. [23].. Wu an P. H. oberts, Shock-wave propagation in a sonoluminescing gas bubble. Phys. ev. Lett. 70, 3424 (1993. [24] W.. Moss, D. B. larke, J. W. White, an D. A. Young, Hyroynamics simulation of bubble collapse an picosecon sonoluminescence. Phys. Fluis 6, 2979 (1994; W.. Moss, D. B. lark, an D. A. Young, alculate pulse withs an spectra of a single sonoluminescing bubble. Scince, 276, 1398 (1997. [25] L. Konic, J. I. Gersten, an. Yuan, Theoretical stuies of sonoluminescence raiations: raiative transfer an parametric epenence. Phys. ev. E 52, 4976 (1995. [26] V. Q. Voung an A. J. Szeri, Sonoluminescence an iffusive transport. Phys. Fluis 8, 2354 (1996. [27] L. Yuan, H. Y. heng, M.-. hu, an P. T. Leung, Physical parameters affecting sonoluminescence: A selfconsistent hyroynamics stuy. Phys. ev. E 57, 4265 (1998. [28] N. Xu,. Apfel, A. Khong, X. Hu, an L. Wang, Water vapor iffusion effects on gas ynamics in a sonoluminescing bubble. Phys. ev. E 68, (2003. [29] B. D. Storey an A. J. Szeri, Water vapour, sonoluminescence an sonochemistry. Proc. oy. Soc. Lonon, Ser. A 456, 1685 (2000; [30]. Toegel, B. Gompf,. Pecha, an D. Lohse, Does water vapor prevent upscaling sonoluminescence. Phys. ev. Lett. 85, 3165 (2000. [31]. Toegel, D. Lohse, Phase iagrams for sonoluminescing bubbles: A comparison between theory an experiment. J. hem. Phys. 118, 1863 (2003. [32] X. Lu, A. Prosperetti,. Toegel, an D. Lohse, Harmonic enhancement of single-bubble sonoluminescence. Phys. ev. E. 67, (2003. [33]. Toegel, S. Hilgenfelt, an D. Lohse, Suppressing issociation in sonoluminescing bubbles: The effect of exclue volume. Phys. ev. Lett. 88, (2002. [34] Hanbook of hemistry an Physics, eite by D. Lie, Press, Boca aton, (1995. [35] S. M. Karim, L. osenhea, The secon oefficient of viscosity of liqui an gases. ev. Mo. Phys.24, 108 (1952. [36] J. A. Ketterling an. E. Apfel, Extensive experimental mapping of sonoluminescence parameter space. Phys. ev. E 61, 3832 (2000. [37] G. Simon, I. sabai, A. Horvath, an F. Szalai, Sonoluminescence an phase iagrams of single bubbles at low isolve air concentration. Phys. ev. E 63, (2001. [38] G. Vazquez,. amara, S. J. Putterman, an K. Weninger, Blackboy spectra for sonoluminescing hyrogen bubbles. Phys. ev. Lett. 88, (2002. [39] W.. Moss, J. L. levatin, an A. J. Szeri, A new amping mechanism in strongly collapsing bubbles. Proc. oy. Soc. Lonon, Ser. A 456, 2983 (2000.