# Acoustic Characterization of Fluorinert FC-43 Liquid with Helium Gas Bubbles: Numerical Experiments

Authors: Christian Vanhille; Cristian Pantea; Dipen N. Sinha
Journal: Shock and Vibration
 (2017)

---

## Abstract

In this work, we define the acoustic characteristics of a biphasic fluid consisting of static helium gas bubbles in liquid Fluorinert FC-43 and study the propagation of ultrasound of finite amplitudes in this medium. Very low sound speed and high sound attenuation are found, in addition to a particularly high acoustic nonlinear parameter. This result suggests the possibility of using this medium as a nonlinear enhancer in various applications. In particular, parametric generation of low ultrasonic frequencies is studied in a resonator cavity as a function of driving pressure showing high conversion efficiency. This work suggests that this medium could be used for applications such as parametric arrays, nondestructive testing, diagnostic medicine, sonochemistry, underwater acoustics, and ultrasonic imaging and to boost the shock formation in fluids.


---

## Body

## 1. Introduction

This paper deals with the propagation of ultrasound in a biphasic bubbly fluid comprised of Fluorinert FC-43 liquid and helium gas bubbles: (FC-43)-He. The primary objective of this paper is to conduct a theoretical study of this interesting medium to determine whether it is a good candidate for generating low ultrasonic frequency signal from the frequency mixing of two high frequency ultrasonic waves. Therefore, we first evaluate the acoustic parameters, such as sound speed, sound attenuation coefficient, compressibility coefficient, and acoustic nonlinear parameter of this bubbly fluid through numerical simulations. These parameters then allow us to determine the frequency mixing efficiency of this unique medium in creating low frequency ultrasonic waves using nonlinear parametric theory and numerical experiments.Bubbles created by acoustic cavitation are employed in sonochemistry and other industrial processes, such as cleaning, mixing, and water treatment [1–6]. Sonoluminescence has been observed under some particular circumstances [7]. Bubbles are also used in ultrasonic diagnosis and therapy [8, 9].Fluorinert FC-43 (3M™) was chosen for this study due to its low sound speed (646 m/s) and relatively large (for a liquid) acoustic nonlinearity parameterβ (6.8–7.6) [10–13]. Fluorinert is an electrically insulating, stable hydrocarbon-based fluid used in various cooling applications of electronic components. It has several very useful properties, such as a very low viscosity, excellent chemical and thermal stability below its boiling point (~165°C), and compatibility with sensitive materials, and the liquid is practically nontoxic. The dielectric constant of Fluorinert is approximately 10 times than that of air, deeming it safe for nonlinear acoustics studies, which typically require high excitation power. Fluorinert FC-43 was used in the past to generate a collimated, narrow beam of low frequency (15–120 kHz), using primary frequencies around 1 MHz [11, 14].A mixture of Fluorinert and gas bubbles was chosen for the present study because of the well-known acoustic nonlinearity enhancement found in other bubbly fluids. For example, water’s acoustic nonlinearity increases drastically when air bubbles are added to the liquid [15–19].Helium (He-4) was chosen because of the lower acoustic attenuation coefficient of the resulting bubbly fluid as compared to the fluid mixture made of air bubbles (see Table1). Helium gas also allows the generation of much smaller bubbles with ease than with air. The difference frequency signal that can be generated using this novel bubbly fluid through frequency mixing because of its enhanced acoustic nonlinearity is expected to be much stronger than that is possible with air bubbles. The combination of enhanced nonlinearity, low sound speed, and low acoustic attenuation can facilitate the generation of highly collimated low frequency ultrasonic beams [11, 14]. It is worth stressing the potential for practical applications of this unique fluid medium that include parametric arrays, nondestructive testing, and diagnostic medicine [1, 2, 6, 9]. The low sound speed of FC-43 allows the development of smaller length devices as the wavelength is considerably shorter than in usual liquids used. For example, the wavelength in FC-43 is approximately 3 times shorter than in water.Table 1
Values of sound speed, attenuation coefficient, and compressibility coefficient of FC-43 liquid with helium gas bubbles (bold font values), of FC-43 liquid with air bubbles, and of water with air bubbles. Driving frequency isfr=20 kHz, bubble radius is R0r=7 μm, and bubble density in the liquid is Nr=1×1011/m3.
AtNr - R0r - fr c B f (m/s) α B f (/m) κ B f (/Pa) FC-43 & helium 501.12 0.23 2.14 × 10−9 FC-43 & air 488.3 0.29 2.26 × 10−9 Water & air 842.28 0.082 1.41 × 10−9For the reasons mentioned above, it is important to first obtain the acoustic characteristics of this novel bubbly fluid (FC-43)-He as no such data currently exist in the literature. In this paper, we present a theoretical approach to determine such characteristics of this bubbly fluid. Section2 presents the mathematical model used to perform the theoretical analysis. In Section 3.1, the sound speed, attenuation coefficient, and compressibility coefficient of the bubbly liquid are determined. The nonlinear propagation of finite-amplitude ultrasound in the medium is presented in Section 3.2, and the evaluation of the acoustic nonlinear parameter of the medium that turns out to be very large is discussed in Section 3.3. This larger nonlinearity is employed to analyze the generation of the difference frequency component from frequency mixing in Section 3.4. Section 4 presents the conclusions of this work.


## 2. Materials and Methods

We consider the following differential system formed by the Rayleigh-Plesset and the wave equations, valid for the modeling of nonlinearly oscillating gas bubbles interacting in a liquid with a nonlinear ultrasonic field:(1)∂2v∂t2+δω0∂v∂t+ω02v=av2+b2v∂2v∂t2+∂v∂t2-ηp,0≤x≤L,0<t<T,∂2p∂x2-1cl2∂2p∂t2=-ρlN∂2v∂t2,0<x<L,0<t<T,p=∂p∂t=0,t=0,0<x≤L,v=∂v∂t=0,t=0,0≤x≤L,px=0=p0sin⁡ωt,0<t≤T,px=L=gLt,0<t≤T,where v and p are the coupled bubble volume variation and acoustic pressure: both are functions of the one-dimensional space coordinate x and time t. L is the spatial dimension of the study performed until the last instant T.In the wave equation,cl is the small amplitude sound speed of liquid, ρl is the equilibrium density of liquid, and N is the bubble density in the liquid.In the Rayleigh-Plesset equation,ω0=3γgpg/ρlR02 is the bubble resonance, in which γg is the specific heat ratio of gas, R0 is the initial radius of bubbles, and pg=ρgcg2/γg is the atmospheric pressure of gas, where ρg and cg are the equilibrium density and small amplitude sound speed of gas, δ=4μl/ω0R02 is the viscous damping coefficient of the bubbly liquid, for which μl is the kinematic viscosity of the liquid, a=γg+1ω02/2v0, b=1/6v0, where v0 is the initial volume of bubbles, and η=4πR0/ρl.The source of frequencyf, pulsation ω=2πf, and amplitude p0 excites the fluid from x=0. The time-dependent function gLt describes the behavior of ultrasound at x=L, depending on the kind of problem we solve (standing or propagating waves).The system takes into account the nonlinearity, dispersion, and dissipation due to the gas bubbles in the liquid. The model is described in detail in [15, 16, 18]. Acoustic streaming, buoyancy, Bjerknes forces, and bubble collapse are not considered [19].Two theoretical approaches are used here. The first approach corresponds to the application of an analytic perturbation method to the second order that allows us to calculate the sound speed, attenuation coefficient, and compressibility coefficient in the bubbly liquid and this is presented in Section3.1. This approach is also described in detail in [17, 18]. The second approach corresponds to the application of the Snow-Bl code [15, 20] that allows us to describe the propagation of nonlinear ultrasound in the bubbly liquid, by tracking variables p and v in time, and is discussed in Section 3.2. This is to evaluate the nonlinear parameter of the medium, presented in Section 3.3, and also to study the medium’s capability to create a difference frequency component from two finite-amplitude signals by nonlinear mixing, discussed in Section 3.4. The numerical model is described in [15, 20].


## 3. Results and Discussion

This section presents the results obtained from applying the two approaches described in Section2 to (FC-43)-He. These are presented in Section 3.1, and in Sections 3.2–3.4, respectively.


### 3.1. Acoustic Characteristics of the Bubbly Fluid

It is well-known that the addition of even a relatively low concentration of gas bubbles with homogeneous density distribution into a liquid drastically modifies the acoustic properties of the liquid but hardly changes its density. In particular, the bubbly liquid becomes dispersive and its compressibility increases to a much higher value than for the homogeneous liquid in certain frequency ranges. In this section, we evaluate the sound speedcBf, attenuation coefficient αBf, and compressibility coefficient κBf of (FC-43)-He mixture, for several bubble sizes and over a large frequency range by conserving the same bubble density in the liquid (Nr=1×1011/m3) (Figure 1(a)). In particular, we show the values for a specific bubble radius (R0r=7 μm) (Figure 1(b)) and for a specific frequency (fr=20 kHz) (Table 1, bold font values). These specific values will be used as reference values in the rest of the paper. The evaluation of the parameters follows the first approach described in Section 2 [17, 18]. Table 1 also lists the corresponding results for Fluorinert FC-43 liquid with air bubbles and for water with air bubbles. The parameters used in this paper are the following: sound speed, density, and kinematic viscosity of Fluorinert FC-43 liquid as cl=646 m/s, ρl=1860 kg/m3, and μl=2.5×10-6 m2/s and sound speed, density, and specific heat ratio of gas helium as cg=1016 m/s, ρg=0.1637 kg/m3, and γg=5/3, respectively.Figure 1
Sound speed, attenuation coefficient, and compressibility coefficient of the bubbly liquid (a) versus bubble size and driving frequency and (b) versus driving frequency for the specific bubble radiusR0r=7 μm. Bubble density is Nr=1×1011/m3.
(a) (b)It can be seen in Figure1 that cBf, αBf, and κBf are frequency dependent. This dispersive nature affects the nonlinear ultrasonic field since each of its harmonic components travels at its own sound speed with its associated attenuation. It can be seen that at a frequency fr=20 kHz, situated well below the resonance frequency of the bubbles present, the frequency range of largest attenuation in the system under consideration can be avoided. This frequency corresponds to a low sound speed and high compressibility of the medium and this combination enhances its nonlinearity [10, 17–19].


### 3.2. Ultrasonic Propagation in the Bubbly Fluid

The compressibility change obtained by adding helium gas bubbles to Fluorinert FC-43 liquid (Figure1 and Table 1) implies a higher nonlinearity of the medium at certain frequency ranges [17–19]. This modification affects the behavior of ultrasound propagation through the bubbly liquid. The simulation of the nonlinear propagation of ultrasound in (FC-43)-He is performed here by the second approach described in Section 2 [15, 20]. In this section, we calculate the response obtained in the medium for bubbles of radius R0r that are evenly distributed in the liquid at density Nr for a continuous ultrasonic signal of amplitude 20 kPa driven at fr. Figure 2 shows the pressure waveforms obtained at several distances x from the source with the corresponding frequency content.Figure 2
(a) Pressure waveform and (b) frequency content at several distances from the source. Bubble radius isR0r=7 μm, bubble density is Nr=1×1011/m3, driving frequency is fr=20 kHz, and pressure amplitude at the source is p0=20 kPa.
(a) (b)As the wave travels through the bubbly liquid, the strong attenuation due to the bubbles reduces the amplitude and the signal undergoes a strong nonlinear distortion that is made evident by the asymmetry between compression and rarefaction amplitudes and by the creation of higher harmonic components.


### 3.3. Determination of the Nonlinear Parameter of the Bubbly Fluid

The numerical results presented in the previous section revealed the strong nonlinear character of ultrasound at quite high amplitude in (FC-43)-He mixture. In this section, we evaluate the nonlinear parameterβ of this bubbly liquid by applying the second approach described in Section 2 [15, 20]. The acoustic nonlinear parameter, β=1+B/2A, measures the capacity of a medium to distort nonlinearly an acoustic wave, where B/A is a measure of the nonlinearity in the equation of state of a material [19, 21–24]. To this purpose, we use the finite-amplitude technique described in [25]. The amplitude values of the fundamental and second harmonic are calculated in the configuration shown in Section 3.2 (Figure 2) but for a lower amplitude at the source (1800 Pa). This is to avoid the formation of higher harmonic components (from third to higher) that can adversely affect the application of this technique. Figure 3 shows the second harmonic, p2, as a function of the square of the fundamental, p1, up to a distance of x=30.4 cm from the source, as well as a linear fit of the data. The slope of the linear fit is S=3.71×10-4/mPa. The nonlinear parameter β is defined by β=SρBfcBf3/πf [25–27], where ρBf is the density of the bubbly liquid and cBf is the sound speed in the bubbly liquid at frequency f. The β value is thus 1380. This nonlinear parameter value is several orders of magnitude higher than the β value measured in homogeneous Fluorinert FC-43 liquid, which is 7.6 [11].Figure 3
Second harmonic of pressure as a function of the square of the fundamental versus distance from the source. Bubble radius isR0r=7 μm, bubble density is Nr=1×1011/m3, driving frequency is fr=20 kHz, and pressure amplitude at the source is p0=1800 Pa.


### 3.4. Difference Frequency Generation in the Bubbly Fluid

The acoustic nonlinearity of the bubbly medium (FC-43)-He reflected by the very high nonlinear parameterβ obtained in Section 3.3 leads to the prospect of using this medium for generating a low difference frequency signal from two high frequency signals by nonlinear mixing [11, 14, 28]. A cavity filled with the bubbly medium (FC-43)-He is considered. Here, the second approach described in Section 2 is used [29]. The difference frequency fd is chosen to be fr by setting the primary frequencies at the source to f1=800 kHz and f2=820 kHz. The length of the cavity is set to fit (3/4)th of the wavelength at fd in the bubbly liquid. The pressure amplitude at the source p0 is set at 20 kPa. Figure 4(a) displays the ratio of the difference frequency amplitude pd to p0 expressed in percent versus distance x from the source in the resonator. Figure 4(b) represents the maximal value of pd found in the resonator when p0, expressed in kPa, is varied from very small values up to 50 kPa.Figure 4
(a) Ratio ofpd to p0 (%) versus distance from the source in the resonator and (b) maximal value of pd in the resonator versus p0 (in kPa). Bubble radius is R0r=7 μm, bubble density is Nr=1×1011/m3, primary frequencies are f1=800 kHz and f2=820 kHz, and difference frequency is fd=20 kHz. Pressure amplitude at the source is p0=20 kPa in (a) and is varied up to 50 kPa in (b).
(a) (b)Figure4(a) shows a nonhomogeneous distribution in the cavity in terms of amplitude of difference frequency and shows that the point of maximal amplitude is located at 13 mm from the source. Figure 4(b) indicates that at low amplitude no difference frequency is present in the cavity. The difference frequency component is due to the nonlinearity of the medium that is decisive at finite-amplitude pressures. This means that the difference frequency amplitude increases with pressure amplitude at the source. The fitting curve of this component indicates that the response of the system follows a fourth-degree polynomial curve.


## 3.1. Acoustic Characteristics of the Bubbly Fluid

It is well-known that the addition of even a relatively low concentration of gas bubbles with homogeneous density distribution into a liquid drastically modifies the acoustic properties of the liquid but hardly changes its density. In particular, the bubbly liquid becomes dispersive and its compressibility increases to a much higher value than for the homogeneous liquid in certain frequency ranges. In this section, we evaluate the sound speedcBf, attenuation coefficient αBf, and compressibility coefficient κBf of (FC-43)-He mixture, for several bubble sizes and over a large frequency range by conserving the same bubble density in the liquid (Nr=1×1011/m3) (Figure 1(a)). In particular, we show the values for a specific bubble radius (R0r=7 μm) (Figure 1(b)) and for a specific frequency (fr=20 kHz) (Table 1, bold font values). These specific values will be used as reference values in the rest of the paper. The evaluation of the parameters follows the first approach described in Section 2 [17, 18]. Table 1 also lists the corresponding results for Fluorinert FC-43 liquid with air bubbles and for water with air bubbles. The parameters used in this paper are the following: sound speed, density, and kinematic viscosity of Fluorinert FC-43 liquid as cl=646 m/s, ρl=1860 kg/m3, and μl=2.5×10-6 m2/s and sound speed, density, and specific heat ratio of gas helium as cg=1016 m/s, ρg=0.1637 kg/m3, and γg=5/3, respectively.Figure 1
Sound speed, attenuation coefficient, and compressibility coefficient of the bubbly liquid (a) versus bubble size and driving frequency and (b) versus driving frequency for the specific bubble radiusR0r=7 μm. Bubble density is Nr=1×1011/m3.
(a) (b)It can be seen in Figure1 that cBf, αBf, and κBf are frequency dependent. This dispersive nature affects the nonlinear ultrasonic field since each of its harmonic components travels at its own sound speed with its associated attenuation. It can be seen that at a frequency fr=20 kHz, situated well below the resonance frequency of the bubbles present, the frequency range of largest attenuation in the system under consideration can be avoided. This frequency corresponds to a low sound speed and high compressibility of the medium and this combination enhances its nonlinearity [10, 17–19].


## 3.2. Ultrasonic Propagation in the Bubbly Fluid

The compressibility change obtained by adding helium gas bubbles to Fluorinert FC-43 liquid (Figure1 and Table 1) implies a higher nonlinearity of the medium at certain frequency ranges [17–19]. This modification affects the behavior of ultrasound propagation through the bubbly liquid. The simulation of the nonlinear propagation of ultrasound in (FC-43)-He is performed here by the second approach described in Section 2 [15, 20]. In this section, we calculate the response obtained in the medium for bubbles of radius R0r that are evenly distributed in the liquid at density Nr for a continuous ultrasonic signal of amplitude 20 kPa driven at fr. Figure 2 shows the pressure waveforms obtained at several distances x from the source with the corresponding frequency content.Figure 2
(a) Pressure waveform and (b) frequency content at several distances from the source. Bubble radius isR0r=7 μm, bubble density is Nr=1×1011/m3, driving frequency is fr=20 kHz, and pressure amplitude at the source is p0=20 kPa.
(a) (b)As the wave travels through the bubbly liquid, the strong attenuation due to the bubbles reduces the amplitude and the signal undergoes a strong nonlinear distortion that is made evident by the asymmetry between compression and rarefaction amplitudes and by the creation of higher harmonic components.


## 3.3. Determination of the Nonlinear Parameter of the Bubbly Fluid

The numerical results presented in the previous section revealed the strong nonlinear character of ultrasound at quite high amplitude in (FC-43)-He mixture. In this section, we evaluate the nonlinear parameterβ of this bubbly liquid by applying the second approach described in Section 2 [15, 20]. The acoustic nonlinear parameter, β=1+B/2A, measures the capacity of a medium to distort nonlinearly an acoustic wave, where B/A is a measure of the nonlinearity in the equation of state of a material [19, 21–24]. To this purpose, we use the finite-amplitude technique described in [25]. The amplitude values of the fundamental and second harmonic are calculated in the configuration shown in Section 3.2 (Figure 2) but for a lower amplitude at the source (1800 Pa). This is to avoid the formation of higher harmonic components (from third to higher) that can adversely affect the application of this technique. Figure 3 shows the second harmonic, p2, as a function of the square of the fundamental, p1, up to a distance of x=30.4 cm from the source, as well as a linear fit of the data. The slope of the linear fit is S=3.71×10-4/mPa. The nonlinear parameter β is defined by β=SρBfcBf3/πf [25–27], where ρBf is the density of the bubbly liquid and cBf is the sound speed in the bubbly liquid at frequency f. The β value is thus 1380. This nonlinear parameter value is several orders of magnitude higher than the β value measured in homogeneous Fluorinert FC-43 liquid, which is 7.6 [11].Figure 3
Second harmonic of pressure as a function of the square of the fundamental versus distance from the source. Bubble radius isR0r=7 μm, bubble density is Nr=1×1011/m3, driving frequency is fr=20 kHz, and pressure amplitude at the source is p0=1800 Pa.


## 3.4. Difference Frequency Generation in the Bubbly Fluid

The acoustic nonlinearity of the bubbly medium (FC-43)-He reflected by the very high nonlinear parameterβ obtained in Section 3.3 leads to the prospect of using this medium for generating a low difference frequency signal from two high frequency signals by nonlinear mixing [11, 14, 28]. A cavity filled with the bubbly medium (FC-43)-He is considered. Here, the second approach described in Section 2 is used [29]. The difference frequency fd is chosen to be fr by setting the primary frequencies at the source to f1=800 kHz and f2=820 kHz. The length of the cavity is set to fit (3/4)th of the wavelength at fd in the bubbly liquid. The pressure amplitude at the source p0 is set at 20 kPa. Figure 4(a) displays the ratio of the difference frequency amplitude pd to p0 expressed in percent versus distance x from the source in the resonator. Figure 4(b) represents the maximal value of pd found in the resonator when p0, expressed in kPa, is varied from very small values up to 50 kPa.Figure 4
(a) Ratio ofpd to p0 (%) versus distance from the source in the resonator and (b) maximal value of pd in the resonator versus p0 (in kPa). Bubble radius is R0r=7 μm, bubble density is Nr=1×1011/m3, primary frequencies are f1=800 kHz and f2=820 kHz, and difference frequency is fd=20 kHz. Pressure amplitude at the source is p0=20 kPa in (a) and is varied up to 50 kPa in (b).
(a) (b)Figure4(a) shows a nonhomogeneous distribution in the cavity in terms of amplitude of difference frequency and shows that the point of maximal amplitude is located at 13 mm from the source. Figure 4(b) indicates that at low amplitude no difference frequency is present in the cavity. The difference frequency component is due to the nonlinearity of the medium that is decisive at finite-amplitude pressures. This means that the difference frequency amplitude increases with pressure amplitude at the source. The fitting curve of this component indicates that the response of the system follows a fourth-degree polynomial curve.


## 4. Conclusions

This work shows the high acoustic nonlinearity of a unique medium composed of Fluorinert FC-43 liquid in which a population of tiny helium gas bubbles is evenly distributed. The medium is characterized acoustically by evaluating the sound speed, attenuation parameter, compressibility coefficient, and nonlinear parameter. The effects of this medium on the propagation of ultrasound are presented, as well as its capability for generating difference frequency components. This work suggests the potential use of Fluorinert FC-43 liquid with helium gas bubbles as a nonlinear enhancer in various practical applications, including shock formation in fluids.



---

*Source: 2518168-2017-01-19.xml*