# Nonlinear Resonance of Cavities Filled with Bubbly Liquids: A Numerical Study with Application to the Enhancement of the Frequency Mixing Effect

Authors: María Teresa Tejedor Sastre; Christian Vanhille
Journal: Shock and Vibration
 (2018)

---

## Abstract

This paper studies the nonlinear resonance of a cavity filled with a nonlinear biphasic medium made of a liquid and gas bubbles at a frequency generated by nonlinear frequency mixing. The analysis is performed through numerical simulations by mixing two source signals of frequencies well below the bubble resonance. The finite-volume and finite-difference based model developed in the time domain simulates the nonlinear interaction of ultrasound and bubble dynamics via the resolution of a differential system formed by the wave and Rayleigh–Plesset equations. Some numerical results, consistent with the literature, validate our procedure. Other results reveal the existence of a frequency shift of the cavity resonance at the difference-frequency component, which rises with pressure amplitude and evidences the global changes undergone by the bubbly medium under finite amplitudes. Finally, this work shows the enhancement of the amplitude of the difference-frequency component generated by parametric excitation using the nonlinear resonance shift, which is more pronounced when the second primary frequency is constant, the first one is varied to match the nonlinear resonance, and both have the same amplitude.


---

## Body

## 1. Introduction

Adding bubbles to a liquid modifies its acoustic properties [1–5]. The nonlinear parameter increases by several orders of magnitude. The sound speed, attenuation coefficient, compressibility, and nonlinear parameter acquire dispersive dependence on bubble resonance. The nonlinear interaction of ultrasound and bubble oscillations must be understood to take advantage of these properties in different applied frameworks such as sonochemistry [6], medicine [7], and others [8, 9]. Lauterborn, in [10], studies the nonlinear behavior of a single bubble in an acoustic field to analyze the effect of the pressure amplitude on the bubble resonance and concludes that a shift of the bubble resonance exists and is dependent on pressure amplitude.The nonlinearity of the medium is responsible for the generation of harmonics from the fundamental frequency and generates combinations of frequencies by nonlinear frequency mixing (sum frequency and difference frequency) when several ultrasonic signals travel through the medium [11]. These effects have multiple applications. Medical imaging can be generated from higher harmonic components [12]. Underwater exploration or transmission and nondestructive testing are fields where the difference-frequency signal has a huge interest because of its low attenuation, good directivity, and high penetration [13, 14]. Characterization and detection of bubbles are also attractive applications of the frequency mixing phenomenon [15–18].Several studies based on linear models have been performed to understand the behavior of ultrasonic waves in bubbly liquids inside a cavity [19–21]. Omta studied the behavior of a bubbly liquid cloud in [22] showing that the nonlinear response emitted from the cloud, much lower than the bubble resonance, is determined mainly its total gas content. Other studies that analyze the behavior of standing ultrasonic waves are based on nonlinear models [23, 24]. In those papers, both the sound speed and the resonance frequencies are calculated without taking into account the amplitude of the waves [2, 3]. In this paper, we aim at showing that the pressure amplitude of the signal changes the resonance of the cavity (and the sound speed).The dependence of the resonance frequency on drive amplitude has been observed in solids, for which the nonlinear features of ultrasound are used in areas as damage diagnostics in materials [25], granular media and dynamic earthquake triggering [26], and fluids in closed tubes of variable cross section [27]. Omta also analyzed in [22] the signal emitted from a bubbly liquid cloud as a function of the amplitude of the acoustic perturbation, concluding that the frequency of this signal undergoes a variation that is amplitude dependent. Up to our knowledge, that paper, and more specifically its Figures 4–6, was the very first demonstration of the shift of the resonance of a bubbly cloud with pressure amplitude. Matsumoto and Yoshizawa, in [28], also detected the shift with pressure amplitude of the resonance of a cluster containing a bubbly liquid. This effect has also been studied in bubbly liquids for a resonance frequency associated to the multiple scattering of bubbles that changes as a function of the amplitude of an incident Gaussian pulse [29].The objective of this work is to study the variation with pressure amplitude of the resonance of a one-dimensional resonator filled with a fluid made of a liquid and gas bubbles when working at nonlinear regime by mixing two finite-amplitude continuous excitation signals. Frequencies well below the bubble resonance are used to take advantage of the nonlinearity of the dispersive medium with a relative low attenuation.In Section2, we present the physical problem and the corresponding mathematical model used in this work. Several numerical experiments performed by varying the amplitude at the source are shown in Section 3. They allow us to observe the nonlinear resonance phenomenon of the cavity at the difference-frequency component generated by nonlinear frequency mixing. This resonance frequency shift is used to maximize its amplitude. Similarities with classic results are also commented. Section 4 gives the conclusions of this work.


## 2. Materials and Methods

We consider a one-dimensional cavity of lengthL filled with a mixture of water and air bubbles. Under the Rayleigh–Plesset approximation, we suppose that, among others, the bubbles are spherical and have the same size. We also assume that they are evenly distributed in the liquid. The model assumes that bubbles are the only source of attenuation, dispersion, and nonlinearity. The buoyancy and Bjerknes and viscous drag forces are not considered in this work. The interaction between the acoustic pressure px,t and the volume variation of the bubbles vx,t=Vx,t−v0g is modeled by the wave equation, Equation (1), and a Rayleigh–Plesset equation, Equation (2) [3, 30], where x is the one-dimensional space coordinate, t is the time, V is the current volume of the bubble, and v0g=4/3πR0g3 is the initial bubble volume, with R0g as the initial radius.(1)∂2p∂x2−1c0l2∂2p∂t2=−ρ0lNg∂2v∂t2,0<x<L,0<t<Tt,(2)∂2v∂t2+δω0g∂v∂t+ω0g2v+ηp=av2+b2v∂2v∂t2+∂v∂t2,0≤x≤L,0<t<Tt.In Equation (1), c0l and ρ0l are the sound speed and the density at the equilibrium state of the liquid. Ng is the density of bubbles, i.e., the bubble number per m3. In Equation (2), δ=4νl/ω0gR0g2 is the viscous damping coefficient of the bubbly fluid, in which νl is the cinematic viscosity of the liquid and ω0g=3γgp0g/ρ0lR0g2 is the resonance frequency of the bubbles, where γg is the specific heats ratio of the gas, p0g=ρ0gc0g2/γg is its atmospheric pressure, and ρ0g and c0g are the density and sound speed at the equilibrium state of the gas. The parameter η=4πR0g/ρ0l and the nonlinear coefficients a=γg+1ω0g2/2v0g and b=1/6v0g are constant. The numerical experiments last a total time Tt. In the following studies, Section 3, the value of this parameter Tt is high enough to ensure that the steady state of the waves is reached. The system is closed by supposing that the liquid and the bubbles are unperturbed at the onset of the studies:(3)px≠0,0=0,vx,0=0,∂p∂tx≠0,0=0,∂v∂tx,0=0,0≤x≤L,and the resonator is excited by a time-dependent pressure source st placed at x=0:(4)px=0,t=st,0≤t≤Tt,and a free-wall condition is imposed at the reflector:(5)pL,t=0,0≤t≤Tt.To solve this differential system, we use a numerical model developed in [24] based on a finite volume method in the space dimension and a finite difference method in the time domain. The frequency components of the time-dependent solution used in the next section are obtained by applying a fast Fourier transform.


## 3. Results

The objective of this section is, using the phenomenon known as nonlinear frequency shift [22, 28, 29], to show the enhancement of the difference-frequency component generated in a cavity that contains a bubbly liquid by nonlinearly mixing two signals of different frequencies.The following data for the bubbly liquid are set into the model:c0l=1500m/s, ρ0l=1000kg/m3, and νl=1.43×10−6m2/s for the liquid (water) and c0g=340m/s, ρ0g=1.29kg/m3, and γg=1.4 for the gas (air). We use bubbles of radius R0g=2.5μm, and the bubble density is Ng=5×1011m−3.Although very few studies exist in the literature, the dependence of the resonance of a bubble cloud on pressure amplitude is a phenomenon, known as nonlinear frequency shift, that has been observed previously in seminal papers by Omta [22], Matsumoto and Yoshizawa [28], and Doc et al. [29]. In Appendix, we show some results obtained with the model described in the above section that corroborate the conclusions of these papers: (i) the increase of pressure amplitude induces the nonlinear frequency shift of the cavity resonance; (ii) this effect relies on the softening of the bubbly liquid, that is due to the variation of the average volume of bubbles; and (iii) this nonlinear effect is more pronounced at higher void fraction in the cavity.We focus now on the discussion on the application of the softening behavior of the bubbly liquid by taking advantage of the nonlinear frequency shift to strengthen the amplitude of the difference-frequency component generated in the context of the nonlinear frequency mixing of two signals of different frequencies [1, 31]. To this purpose, the analysis is performed by means of a comparison of several numerical experiments for which we search the highest response by parametric emission, i.e., the maximal response of the system at the difference frequency fd=f2−f1: (1) by setting the first primary source frequency f1 at a constant value and moving the second primary source frequency f2, with the same constant source amplitude for both the first and the second primary signals p01=p02; (2) by setting f1 at a constant value and moving f2 and considering two subcases, (a) a constant p02 and a varied p01 and (b) a constant p01 and a varied p02; (3) by setting f2 at a constant value and moving f1, with the same constant value p01=p02.To this end, in this section the pressure source we use isst=p01sinω1t+p02sinω2t, where ω1=2πf1 and ω2=2πf2. The cavity length is set to fit the linear resonance at the difference frequency fdL=200kHz, L=λdL/2=0.0031m, where λdL=cdL/fdL is the wavelength and cdL=1222.8m/s is the sound speed in this biphasic and dispersive medium at this frequency [3]. We study the nonlinear resonance shift of the difference-frequency component of the signal pressure in the cavity. For each numerical experiment 1 to 3, simulations are performed varying the source amplitude. For each amplitude, we apply a frequency sweep, in such a way that the difference frequency fd is around the linear resonance fdL, by increments δf=10Hz to evaluate the highest difference-frequency pressure amplitude reached in the cavity pdm at each frequency, and the maximal value pdmax over the frequency range is then localized. Note that we work at primary frequencies chosen to be close to half the resonance frequency of the bubbles since the nonlinearity is high at this frequency in the dispersive medium [3].Case 1. The first primary frequency is constant, f1=700kHz, whereas the second primary source frequency f2 is moved from 896kHz up to 902kHz. The source amplitude is varied from p01=p02=1kPa up to p01=p02=6.5kPa. Figure 1 shows the result, i.e., pdm as a function of frequency (around fdL) over the amplitude range. It is seen here that the behavior of the difference-frequency component is assimilated to others previously observed through other frequencies, i.e., it shows the same main properties when amplitude is raised as the ones described in the literature and in Appendix [22, 28, 29], which means that the amplitude-dependent behavior of the medium can also be characterized by the behavior of fd. The resonance of the cavity at the difference frequency clearly undergoes a dependence on pressure amplitude, i.e., a nonlinear frequency shift exists. This means that the softening of the bubbly liquid in the cavity with pressure amplitudes also affects the difference-frequency component. In this case, at p01=p02=6.5kPa, the resonance is at fd=197.6kHz, denoted by fdNL, the frequency shift is ΔfdNL=2.4kHz, and the highest value is pdNL=19.794kPa, which is 304.5% of the source amplitude. Since L is constant, this frequency shift means that the sound speed in the medium is cdNL=2LfdNL=1208.12m/s. Also, the symmetry of the curves around the linear resonance observed at the lowest amplitudes is lost when the latest rise, the nonlinear attenuation reduces the ratio of pdmax to source amplitude.Figure 1
Maximum difference-frequency pressure amplitude in the cavitypdm vs. frequency (around fdL) for different source amplitudes in Case 1.Figure2 presents the frequency shift Δf obtained as a function of pressure amplitude at the source p01=p02 (a), as a function of average bubble volume increase Δv (b), and the maximal value of difference-frequency pressure amplitude reached in the cavity pdmax over the frequency range as a function of pressure amplitude at the source p01=p02 (c). The fitting curves are also displayed (green color). A 4th degree polynomial fit is obtained for Figure 2(a). This means that the frequency shift increases hugely as the pressure amplitude rises. The linear dependence of the frequency shift observed in Figure 2(b) proves that the softening of the medium is due to the increase of the mean bubble volume that raises the compressibility, and thus, the nonlinearity of the bubbly medium. Figure 2(c) is a consequence of the two other diagrams, and a 3rd degree polynomial fit is obtained, which means that the nonlinearity of the medium due to the increase of the pressure amplitude generates a huge difference-frequency amplitude.Figure 2
Fitting curves in Case 1 of frequency shiftΔf vs. pressure amplitude at the source p01=p02, denoted by p0 in (a) and vs. average volume increase Δv (b) and maximal value of difference-frequency pressure amplitude pdmax over the frequency range vs. pressure amplitude at the source p01=p02, denoted by p0 in (c).
(a)
(b)
(c)The source amplitude for this study isp01=p02=6.5kPa, and we analyze the difference-frequency generation by comparing the results when we do take into account the nonlinear resonance frequency shift and when we do not. The study is thus performed at two difference frequencies, fdL=200kHz, which is the linear resonance, and fdNL=197.6kHz, which is the nonlinear resonance found above at this source amplitude (Figure 1), i.e., the frequency that produces the highest response by parametric emission in this case. The primary frequencies at the source are set at f1L=700kHz and f2L=f1L+fdL=900kHz for the linear resonance case and at f1NL=700kHz and f2NL=f1NL+fdNL=897.6kHz for the nonlinear resonance case. Figure 3 shows the pressure amplitude distribution along the cavity of the primary frequencies f1L, f2L, the difference frequency fdL (continuous lines), the primary frequencies f1NL, f2NL, and the difference frequency fdNL (dashed lines). Different amplitudes are observed for fdL and fdNL. The response at the difference frequency is much better by taking into account the nonlinear frequency shift. Whereas the maximum pressure for fdL is pdL=9.217kPa (141.8% of the source amplitude), the corresponding value for fdNL is pdNL=19.794kPa (304.5% of the source amplitude, which is a very high value for parametric emission). The benefit drawn in terms of difference-frequency amplitude is 162.7%. It is also interesting to note that by moving one of the primary source frequencies, f2 from f2L to f2NL, its amplitude decreases, pf2NL<pf2L, whereas the amplitude of the other primary signal increases, pf1NL>pf1L. This results clearly shows that, by taking into account the resonance frequency shift, the second source frequency is the one that undergoes a strong loss of energy to feed the difference-frequency component that acquires intensity and becomes much stronger, pdNL>pdL. This behavior is most likely due to the fact that f2 is closer to the bubble resonance (see Case 3).Figure 3
Pressure amplitude distribution of frequency components along the resonator in Case 1 atp01=p02=6.5kPa: f1L=700kHz, f2L=900kHz, fdL=200kHz (blue, green, and red continuous lines, respectively) and f1NL=700kHz, f2NL=897.6kHz, fdNL=197.6kHz (blue, green, and red dashed lines, respectively).Case 2. The first primary frequency is constant, f1=700kHz, whereas the second primary source frequency f2 is moved from 896kHz up to 902kHz. Two configurations are considered here. For the first one, the source amplitude of the second primary component is constant, p02=6.5kPa, whereas the source amplitude of the first primary component is varied from p01=6.5kPa up to p01=8.125kPa. Figure 4(a) shows pdm as a function of frequency (around fdL) for three amplitude values. At p01=8.125kPa, the resonance is at fdNL=196.54kHz, the frequency shift is ΔfdNL=3.46kHz, and the highest value is pdNL=22.276kPa, which is 274.2% of the source amplitude. Since L is constant, this frequency shift means that the sound speed in the medium is cdNL=2Lf=1201.65m/s. For the second configuration, the source amplitude of the first primary component is constant, p01=6.5kPa, whereas the source amplitude of the second primary component is varied from p02=6.5kPa up to p02=8.125kPa. Figure 4(b) shows pdm as a function of frequency (around fdL) for three amplitude values. At p02=8.125kPa, the resonance is at fdNL=196.24kHz, the frequency shift is ΔfdNL=3.76kHz, and the highest value is pdNL=23.674kPa, which is 291.4% of the source amplitude. Since L is constant, this frequency shift means that the sound speed in the medium is cdNL=2LfdNL=1199.81m/s.Figure 4
Maximum difference-frequency pressure amplitude in the cavitypdm vs. frequency (around fdL) in Case 2 (constant f1, varied f2) for (a) p02=6.5kPa and three values of p01 and (b) p01=6.5kPa and three values of p02.
(a)
(b)It must be noted that the nonlinear curves are not symmetric. Also, the maximal values obtained for Case 2 are lower, relatively to the source amplitude (even at more source amplitude), than the one obtained for Case 1, for which the same amplitude at the source is applied to both primary frequencies.Case 3. The second primary frequency is constant, f2=900kHz, whereas the first primary source frequency f1 is moved from 698kHz up to 704kHz. The source amplitude is set at p01=p02=6.5kPa. Figure 5 shows the comparison of pdm as a function of frequency (around fdL) over the amplitude range in this case and in Case 1. The general behavior observed in Case 1 becomes apparent here as well, especially a clear nonlinear frequency shift of resonance at the difference frequency compared to the linear resonance. However, the efficiency of the mixing-frequency process is higher in Case 3, giving a higher amplitude of the difference-frequency component, pdNL=20.825kPa (320.4% of the source amplitude, at fdNL=197.38kHz with ΔfdNL=2.62kHz) instead of pdNL=19.794kPa (304.5% of the source amplitude, at fdNL=197.6kHz with ΔfdNL=2.4kHz) in Case 1. In Case 3, the sound speed in the medium is cdNL=2LfdNL=1206.83m/s instead of cdNL=2LfdNL=1208.12m/s in Case 1.Figure 5
Maximum difference-frequency pressure amplitude in the cavitypdm vs. frequency (around fdL) for p01=p02=6.5kPa in Case 1 (constant f1, varied f2, green line) and Case 3 (varied f1, constant f2, red line).This effect is the most likely due to the following. At finite amplitude, the bubble resonance also undergoes a variation from its linear value,f0g=ω0g/2π=1.34MHz, toward lower frequencies, and this frequency shift increases with amplitudes [10]. Thus, when f2 is constant, at high amplitudes the difference between the nonlinear bubble resonance and the primary frequency f2 is significantly reduced, whereas when f1 is constant that difference depends on f2 that is moving (toward lower values when fitting the nonlinear resonance). Since the primary frequency f2 is the one that gives most energy to the difference-frequency component [24], the f2 component of the source then tends to excite the bubbles with more intensity in such a way that the bubbles oscillate around a mean volume that is higher in Case 3 than in Case 1. The softening of the bubbly medium that takes place in the cavity is thus more pronounced in Case 3 than in Case 1 (see Appendix and Refs. [22, 28, 29]), and this gives way to a frequency shift of the difference frequency in the cavity that is higher and produces more intensity in Case 3 than in Case 1.The results obtained in this section suggest the following new points that are of great interest in the framework of nonlinear ultrasound in bubbly liquids: the nonlinear frequency shift of the cavity resonance (decrease of sound speed, softening of the medium by increase of the effective bubble volume with pressure amplitude) can be applied to the mixing of two signals of different frequencies to strengthen the nonlinear generation of the difference frequency component; the comparison of Cases 1 and 2 and Cases 1 and 3 performed above, with primary frequencies well below the bubble resonance, evidences that the enhancement of the difference-frequency component at the nonlinear resonance is very effective (in relation to the source amplitude) when the second primary frequency is constant, whereas the first one is varied to match the nonlinear resonance, and both primary component amplitudes are set at the same value; an unbalanced contribution, in terms of amplitudes, of the primary signals limits the necessary equilibrium to maximize the difference frequency but promotes the generation of harmonics of the strongest primary signal; and a variation of the second primary frequency instead of the first one limits the efficiency of the nonlinear mixing frequency.


## 4. Conclusions

This work shows that a frequency shift, which grows with pressure amplitudes (nonlinear resonance effect), of a system composed by a bubbly liquid in a cavity exists at the difference-frequency component generated by nonlinear frequency mixing of two primary signals at frequencies well below the bubble resonance. This numerical study also analyzes different ways to enhance the intensity of the difference-frequency signal using this nonlinear resonance effect and suggests the use at the source of a constant second primary frequency combined with a varied first primary frequency to adjust the difference frequency at the nonlinear resonance, both at the same amplitude.



---

*Source: 1570508-2018-12-06.xml*