A numerical method for the dynamics of non-spherical cavitation bubbles

A boundary integral numerical method for the dynamics of nonspherical cavitation bubbles in inviscid incompressible liquids is described. Only surface values of the velocity potential and its first derivatives are involved. The problem of solving the Laplace equation in the entire domain occupied by the liquid is thus avoided. The collapse of a bubble in the vicinity of a solid wall and the collapse of three bubbles with collinear centers are considered

Oscillations of a gas pocket on a liquid-covered solid surface

The dynamic response of a gas bubble entrapped in a cavity on the surface of a submerged solid subject to an acoustic field is investigated in the linear approximation. We derive semi-analytical expressions for the resonance frequency, damping and interface shape of the bubble. For the liquid phase, we consider two limit cases: potential flow and unsteady Stokes flow. The oscillation frequency and interface shape are found to depend on two dimensionless parameters: the ratio of the gas stiffness to the surface tension stiffness, and the Ohnesorge number, representing the relative importance of viscous forces. We perform a parametric study and show, among others, that an increase in the gas pressure or a decrease in the surface tension leads to an increase in the resonance frequency until an asymptotic value is reached

Phase Diagrams for Sonoluminescing Bubbles

Sound driven gas bubbles in water can emit light pulses. This phenomenon is called sonoluminescence (SL). Two different phases of single bubble SL have been proposed: diffusively stable and diffusively unstable SL. We present phase diagrams in the gas concentration vs forcing pressure state space and also in the ambient radius vs gas concentration and vs forcing pressure state spaces. These phase diagrams are based on the thresholds for energy focusing in the bubble and two kinds of instabilities, namely (i) shape instabilities and (ii) diffusive instabilities. Stable SL only occurs in a tiny parameter window of large forcing pressure amplitude $P_a \sim 1.2 - 1.5$atm and low gas concentration of less than $0.4\%$ of the saturation. The upper concentration threshold becomes smaller with increasing forcing. Our results quantitatively agree with experimental results of Putterman's UCLA group on argon, but not on air. However, air bubbles and other gas mixtures can also successfully be treated in this approach if in addition (iii) chemical instabilities are considered. -- All statements are based on the Rayleigh-Plesset ODE approximation of the bubble dynamics, extended in an adiabatic approximation to include mass diffusion effects. This approximation is the only way to explore considerable portions of parameter space, as solving the full PDEs is numerically too expensive. Therefore, we checked the adiabatic approximation by comparison with the full numerical solution of the advection diffusion PDE and find good agreement.Comment: Phys. Fluids, in press; latex; 46 pages, 16 eps-figures, small figures tarred and gzipped and uuencoded; large ones replaced by dummies; full version can by obtained from: http://staff-www.uni-marburg.de/~lohse

On the Classical Theory of the Electron

A classical theory of the electron, proposed by one of us several years ago and based on finite-difference equations, is discussed by considering the three possible following cases: radiating electron, absorbing electron and nonradiating, nonabsorbing electron. In particular the so-called transmission laws necessary to determine, in conjunction with the dynamical equations, the motion of a charged particle corresponding to given initial values of position and velocity are critically reconsidered. The general characteristics of the one-dimensional motion in the non-relativistic approximation are discussed in detail. It is found that in the case of the radiating electron the particle position tends asimptotically to the point of stable equilibrium. The present theory is, therefore, free from the unphysical phenomenon of runaway solutions. These general results are illustrated by studying the motion of a particle under the action of a restoring elastic force and under the action of purely time-dependent forces

Bubble Shape Oscillations and the Onset of Sonoluminescence

An air bubble trapped in water by an oscillating acoustic field undergoes either radial or nonspherical pulsations depending on the strength of the forcing pressure. Two different instability mechanisms (the Rayleigh--Taylor instability and parametric instability) cause deviations from sphericity. Distinguishing these mechanisms allows explanation of many features of recent experiments on sonoluminescence, and suggests methods for finding sonoluminescence in different parameter regimes.Comment: Phys. Rev. Lett., in pres

Heat transfer mechanisms in bubbly Rayleigh-Benard convection

The heat transfer mechanism in Rayleigh-Benard convection in a liquid with a mean temperature close to its boiling point is studied through numerical simulations with point-like vapor bubbles, which are allowed to grow or shrink through evaporation and condensation and which act back on the flow both thermally and mechanically. It is shown that the effect of the bubbles is strongly dependent on the ratio of the sensible heat to the latent heat as embodied in the Jacob number Ja. For very small Ja the bubbles stabilize the flow by absorbing heat in the warmer regions and releasing it in the colder regions. With an increase in Ja, the added buoyancy due to the bubble growth destabilizes the flow with respect to single-phase convection and considerably increases the Nusselt number.Comment: 11 pages, 14 figure

An Alternative Method to Deduce Bubble Dynamics in Single Bubble Sonoluminescence Experiments

In this paper we present an experimental approach that allows to deduce the important dynamical parameters of single sonoluminescing bubbles (pressure amplitude, ambient radius, radius-time curve) The technique is based on a few previously confirmed theoretical assumptions and requires the knowledge of quantities such as the amplitude of the electric excitation and the phase of the flashes in the acoustic period. These quantities are easily measurable by a digital oscilloscope, avoiding the cost of expensive lasers, or ultrafast cameras of previous methods. We show the technique on a particular example and compare the results with conventional Mie scattering. We find that within the experimental uncertainties these two techniques provide similar results.Comment: 8 pages, 5 figures, submitted to Phys. Rev.

Orthogonal, solenoidal, three-dimensional vector fields for no-slip boundary conditions

Viscous fluid dynamical calculations require no-slip boundary conditions. Numerical calculations of turbulence, as well as theoretical turbulence closure techniques, often depend upon a spectral decomposition of the flow fields. However, such calculations have been limited to two-dimensional situations. Here we present a method that yields orthogonal decompositions of incompressible, three-dimensional flow fields and apply it to periodic cylindrical and spherical no-slip boundaries.Comment: 16 pages, 2 three-part figure

Cerebellar control of cortico-striatal LTD

Purpose: Recent anatomical studies showed the presence of cerebellar and basal ganglia connections. It is thus conceivable that the cerebellum may influence the striatal synaptic transmission in general, and synaptic plasticity in particular. Methods: In the present neurophysiological investigation in brain slices, we studied striatal long-term depression (LTD), a crucial form of synaptic plasticity involved in motor learning after cerebellar lesions in rats. Results: Striatal LTD was fully abolished in the left striatum of rats with right hemicerebellectomy recorded 3 and 7 days following surgery, when the motor deficits were at their peak. Fifteen days after the hemicerebellectomy, rats had partially compensated their motor deficits and high-frequency stimulation of excitatory synapses in the left striatum was able to induce a stable LTD. Striatal plasticity was conversely normal ipsilaterally to cerebellar lesions, as well as in the right and left striatum of sham-operated animals. Conclusions: These data show that the cerebellum controls striatal synaptic plasticity, supporting the notion that the two structures operate in conjunction during motor learning