Nonlinear Modes of Liquid Drops as Solitary Waves

The nolinear hydrodynamic equations of the surface of a liquid drop are shown to be directly connected to Korteweg de Vries (KdV, MKdV) systems, giving traveling solutions that are cnoidal waves. They generate multiscale patterns ranging from small harmonic oscillations (linearized model), to nonlinear oscillations, up through solitary waves. These non-axis-symmetric localized shapes are also described by a KdV Hamiltonian system. Recently such ``rotons'' were observed experimentally when the shape oscillations of a droplet became nonlinear. The results apply to drop-like systems from cluster formation to stellar models, including hyperdeformed nuclei and fission.Comment: 11 pages RevTex, 1 figure p

Secondary Bjerknes forces between two bubbles and the phenomenon of acoustic streamers

The translational velocities of two spherical gas bubbles oscillating in water, which is irradiated by a high-intensity acoustic wave field, are calculated. The two bubbles are assumed to be located far enough apart so that shape oscillations can be neglected. Viscous effects are included owing to the small size of the bubbles. An asymptotic solution is obtained that accounts for the viscous drag on each bubble, for large Re based on the radial part of the motion, in a form similar to the leading-order prediction by Levich (1962), C-D = 48/ReT; Re-T --> infinity based on the translational velocity. In this context the translational velocity of each bubble, which is a direct measure of the secondary Bjerknes force between the two bubbles, is evaluated asymptotically and calculated numerically for sound intensities as large as the Blake threshold. Two cases are examined. First, two bubbles of unequal size with radii on the order of 100 pm are subjected to a sound wave with amplitude P-A < 1.0 bar and forcing frequency w(f) =0.51w(10), so that the second harmonic falls within the range defined by the eigenfrequencies of the two bubbles, w(10)<2w(f)<w(20). It is shown that their translational velocity changes sign, becoming repulsive as PA increases from 0.05 to 0.1 bar due to the growing second harmonic, 2w(f), of the forcing frequency. However, as the amplitude of sound further increases, P-A approximate to 0.5 bar, the two bubbles attract each other due to the growth of even higher harmonics that fall outside the range defined by the eigenfrequencies of the two bubbles. Second, the case of much smaller bubbles is examined, radii on the order of 10 mum, driven well below resonance, w(f)/2pi = 20 kHz, at very large sound intensities, P-A approximate to 1 bar. Numerical simulations show that the forces between the two bubbles tend to be attractive, except for a narrow region of bubble size corresponding to a nonlinear resonance related to the Blake threshold. As the distance between them decreases, the region of repulsion is shifted, indicating sign inversion of their mutual force. Extensive numerical simulations indicate the formation of bubble pairs with constant average inter-bubble distance, consisting of bubbles with equilibrium radii determined by the primary and secondary resonance frequencies for small and moderate sound amplitudes or by the Blake threshold for large sound amplitudes. It is conjectured that in experiments where 'acoustic streamers' are observed, which are filamentary structures consisting of bubbles that are aligned and move rapidly in a cavitating fluid at nearly constant distances from each other, bubbles with size determined by the Blake threshold are predominant because those with size determined by linear resonance are larger and therefore become unstable due to shape oscillations

Spherical capsules in three-dimensional unbounded Stokes flows: effect of the membrane constitutive law and onset of buckling

The dynamic response of an initially spherical capsule subject to different externally imposed flows is examined. The neo-Hookcan and Skalak et al. (Biophys. J., vol. 13 (1973), pp. 245-264) constitutive laws are used for the description of the membrane mechanics, assuming negligible bending resistance. The viscosity ratio between the interior and exterior fluids of the capsule is taken to be unity and creeping-flow conditions are assumed to prevail. The capillary number epsilon is the basic dimensionless number of the problem, which measures the relative importance of viscous and elastic forces. The boundary-element method is used with bi-cubic B-splines as basis functions in order to discretize the capsule surface by a structured mesh. This guarantees continuity of second derivatives with respect to the position of the Lagrangian particles used for tracking the location of the interface at each time step and improves the accuracy of the method. For simple shear flow and hyperbolic flow, an interval in 8 is identified within which stable equilibrium shapes are obtained. For smaller values of E, steady shapes are briefly captured, but they soon become unstable owing to the development of compressive tensions in the membrane near the equator that cause the capsule to buckle. The post-buckling state of the capsule is conjectured to exhibit small folds around the equator similar to those reported by Walter et al. Colloid Polymer Sci. Vol. 278 (2001), pp. 123-132 for polysiloxane microcapsules. For large values of 6, beyond the interval of stability, the membrane has two tips along the direction of elongation where the deformation is most severe, and no equilibrium shapes could be identified. For both regions outside the interval of stability, the membrane model is not appropriate and bending resistance is essential to obtain realistic capsule shapes. This pattern persists for the two constitutive laws that were used, with the Skalak et al. law producing a wider stability interval than the neo-Hookean law owing to its strain hardening nature