27.12.2012 Views

Oscillations, Waves, and Interactions - GWDG

Oscillations, Waves, and Interactions - GWDG

Oscillations, Waves, and Interactions - GWDG

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

176 R. Mettin<br />

important quantity characterizing every bubble is the mass of non-condensable gas<br />

inside, defining its equilibrium size. On a longer time scale, this size can change due<br />

to diffusion (see below, Sect. 3.3).<br />

The spherical symmetry renders the problem one-dimensional <strong>and</strong> facilitates calculations<br />

significantly. Indeed, this assumption is surprisingly good in many situations<br />

encountered in acoustic cavitation, <strong>and</strong> mainly violated for larger bubbles <strong>and</strong> for<br />

very close boundaries imposed by other bubbles or objects. Instability of sphericity<br />

is also discussed below (Sect. 3.2).<br />

Spherical gas bubble models mainly vary in the Equation Of State (EOS) for<br />

the liquid <strong>and</strong> the gas. The simplest assumption is an incompressible liquid <strong>and</strong><br />

an ideal gas, resulting in an equation frequently referred to as the Rayleigh-Plesset<br />

equation [4,5]. It can be improved by incorporation of liquid compressibility <strong>and</strong> an<br />

advanced EOS for the gas like a Van-der-Waals model. For the problems considered<br />

here, a sufficient compromise between computational effort <strong>and</strong> accuracy is a model<br />

with slight (first order) compressibility of the liquid, namely the Gilmore model [13]<br />

or the Keller-Miksis model [14]. With an ideal gas law the Keller-Miksis model reads<br />

�<br />

1 − ˙ R<br />

c<br />

�<br />

R ¨ R + 3<br />

2 ˙ R 2<br />

�<br />

pl = p0 + 2σ<br />

R0<br />

�<br />

1 − ˙ R<br />

3c<br />

�<br />

=<br />

�<br />

1 + ˙ R<br />

c<br />

�<br />

pl<br />

ρ<br />

R dpl<br />

+<br />

ρc dt ,<br />

� � �3κ R0<br />

− p0 −<br />

R<br />

2σ 4µ<br />

−<br />

R R ˙ R − pa(t) .<br />

Here, R is the bubble radius, c <strong>and</strong> ρ the sound speed <strong>and</strong> the density of the liquid,<br />

µ <strong>and</strong> σ the viscosity <strong>and</strong> the surface tension, <strong>and</strong> κ the polytropic exponent. The<br />

acoustic pressure is denoted pa(t), while the static pressure is p0. R0 denotes the<br />

bubble equilibrium radius, <strong>and</strong> vapour pressure is neglected.<br />

In the following we focus on the volume oscillations, shape instabilities, <strong>and</strong> rectified<br />

diffusion. Afterwards, acoustic forces are addressed which are induced by a cw<br />

acoustic field on a gas bubble. The acoustic forces finally lead to spatial translation<br />

<strong>and</strong> structure formation of the bubbles. For all calculations, the spherical bubble<br />

model is used <strong>and</strong> thus it is the core of all results presented.<br />

3.1 Bubble oscillations<br />

A spherical gas-filled bubble under static conditions posseses an equilibrium radius<br />

R0 (or an equilibrium volume V0 = 4πR 3 0/3) if the static pressure p0 is larger than<br />

the Blake threshold pressure [4]<br />

pB = p0 + 8σ<br />

9<br />

�<br />

3σ<br />

2[p0 + (2σ/R0)]R 3 0<br />

� 1/2<br />

For p0 < pB the bubble would exp<strong>and</strong> infinitely. 5 A momentary excursion of the<br />

bubble radius from the equilibrium results in damped oscillations around R0, the<br />

5 However, in practice it is normally not possible to sustain a static negative pressure<br />

for a longer time under the condition of bubble expansion, i. e., cavitation. Thus bubble<br />

expansion will be stopped sooner or later in real systems.<br />

.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!