Vorticity generation by shock propagation through bubbles in a gas

We present a new theoretical model of ‘late-time’ phenomena related to the interaction of a planar shock with a local, discrete inhomogeneity in an ambient gas. The term ‘late-time’ applies to the evolution of the inhomogeneity and the flow field after interaction with the incident shock has ceased. Observations of a shock propagating through a bubble or a spherical flame have exhibited or implied the formation of vortex structures and have showed continual distortion of the bubble or flame. Our theory shows that this is due to the generation of long-lived vorticity at the edge of the discrete inhomogeneity. The vorticity interacts with itself through the medium of the fluid, and, depending on the geometry of the discrete inhomogeneity, can roll up into vortex filaments or vortex rings. To verify and amplify this theoretical description, we use numerical solutions of the fluid equations for conservation of mass, momentum, and energy to study the interaction of a weak shock with a cylindrical or spherical bubble. The simulated bubble has either a higher or lower density than the ambient gas. In this way, the calculations provide insights into the effects of both geometry and distortion of the local sound speed. The Mach number of the shock is 1.2, the ambient gas is air, and the pressure is 1 atmosphere. Because of the simple geometry of each bubble, the vorticity generated at the boundary rolls up into a vortex filament pair (cylindrical bubble) or a vortex ring (spherical bubble). The structural features and timescales of the phenomena observed in the calculations agree closely with recent experiments of Haas & Sturtevant, in which helium and Freon bubbles were used to provide the local departures from ambient density. The discussion of results includes a survey of alternative numerical methods, sources of uncertainty in velocities of interfaces or structures, as derived from the laboratory and numerical experiments, and the relationship of our analysis to other theories.

[1]  P. Woodward,et al.  The Piecewise Parabolic Method (PPM) for Gas Dynamical Simulations , 1984 .

[2]  Jay P. Boris,et al.  Convective Cooling of Lightning Channels , 1981 .

[3]  Jay P. Boris,et al.  Vorticity Generation by Asymmetric Energy Deposition in a Gaseous Medium. , 1983 .

[4]  T. Maxworthy,et al.  Some experimental studies of vortex rings , 1977, Journal of Fluid Mechanics.

[5]  S. Zalesak Fully multidimensional flux-corrected transport algorithms for fluids , 1979 .

[6]  Jay P. Boris,et al.  The Lagrangian solution of transient problems in hydrodynamics using a triangular mesh , 1979 .

[7]  H. C. Yee Construction of explicit and implicit symmetric TVD schemes and their applications. [Total Variation Diminishing for fluid dynamics computation] , 1987 .

[8]  N. Didden On the formation of vortex rings: Rolling-up and production of circulation , 1979 .

[9]  J. Peraire,et al.  Finite Element Flux-Corrected Transport (FEM-FCT) for the Euler and Navier-Stokes equations , 1987 .

[10]  J. Haas,et al.  Interaction of weak shock waves with cylindrical and spherical gas inhomogeneities , 1987, Journal of Fluid Mechanics.

[11]  A. Leonard Vortex methods for flow simulation , 1980 .

[12]  J. Boris,et al.  Vortex shedding due to laser ablation , 1984 .

[13]  G. I. Taylor Formation of a Vortex Ring by Giving an Impulse to a Circular Disk and then Dissolving it Away , 1953 .

[14]  Elaine S. Oran,et al.  Surface tension and viscosity with Lagrangian hydrodynamics on a triangular mesh. Memorandum report , 1988 .

[15]  J. Boris,et al.  Numerical simulations of asymmetric mixing in planar shear flows , 1986, Journal of Fluid Mechanics.

[16]  P. Woodward,et al.  The numerical simulation of two-dimensional fluid flow with strong shocks , 1984 .

[17]  G. Rudinger,et al.  Behaviour of small regions of different gases carried in accelerated gas flows , 1960, Journal of Fluid Mechanics.

[18]  R. D. Richtmyer Taylor instability in shock acceleration of compressible fluids , 1960 .