A spectacular feature of transient cavity collapse in the neighbourhood of a rigid boundary is the formation of a high-speed liquid jet that threads the bubble and ultimately impacts upon the side of the bubble nearest to the boundary. The bubble then evolves into some toroidal form, the flow domain being doubly connected. In this work, the motion of the toroidal bubble is computed by connecting the jet tip to the side of the bubble upon which it impacts. This connection is via a cut introduced into the flow domain and across which the potential is discontinuous, the value of this discontinuity being equal to the circulation in the flow. A boundary integral algorithm is developed to account for this geometry and some example computations are presented. Consideration of the pressure field in the fluid has implications for possible damage mechanisms to structures due to nearby cavity collapse.

The motion of buoyant transient cavities with non-condensible contents is investigated numerically using a boundary-integral method. The bubble contents are described by an adiabatic gas law. Motion is considered in the neighbourhood of a rigid boundary, in an axisymmetric geometry. We investigate whether the non-condensible contents will resist the formation of jets. It is found that jets form upon collapse and, in general, completely penetrate the bubble before it rebounds, but circumstances are identified under which the non-spherical bubble will rebound prior to this occurrence. In these cases the bulk of the jet growth occurs upon rebound. Furthermore, the interaction between the buoyancy force causing jet formation upwards, and the Bjerknes attraction of the rigid boundary causing jet formation towards it, is investigated and general principles discussed which allow the behaviour to be interpreted. The concept of the Kelvin impulse is utilized.

The Lagally theorem is used to obtain an expression for the Bjerknes force acting on a bubble in terms of the singularities of the fluid velocity potential, defined within the bubble by analytic continuation. This expression is applied to transient cavity collapse in the neighbourhood of boundaries, allowing analytical estimates to be made of the Kelvin impulse of the cavity. The known result for collapse near a horizontal rigid boundary is recovered, and the Kelvin impulse of a cavity collapsing in the neighbourhood of a submerged and partially submerged sphere is estimated. A numerical method is developed to deal with more general body shapes and in particular, bodies of revolution. Noting that the direction of the impulse at the end of the collapse phase generally indicates the direction of the liquid jet that may form, the behaviour of transient cavities in these geometries is predicted. In these examples the concept of a zone of attraction is introduced. This is a region around the body, within which the Kelvin impulse at the time of collapse, and consequent jet formation, is expected to be directed towards the body. Outside this zone the converse is true.

Abstract. The study of the propagation of a shock down a tube of slowly varying cross sectional area has proved to be most valuable in understanding the dynamics of shocks. A particular culmination of this work has been the theory of geometrical shock dynamics due to Whitham (1957Whitham ( , 1959. In this theory the motion of a shock may be approximately computed independently of a determination of the flow field behind the shock. In this paper the propagation of a shock down such a tube is reconsidered. It is found that the motion of the shock is described by an infinite sequence of ordinary differential equations. Each equation is coupled to all of its predecessors but only to its immediate successor, a feature which allows the system to be closed by truncation. Of particular relevance is the demonstration that truncation at the first equation in the sequence yields the A-M relation that is the basis for Whitham's highly successful theory. Truncation at the second equation yields the next level of approximation. The equations so obtained are investigated with analytic solutions being found in the strong shock limit for the propagation of cylindrical and spherical shock waves. Implementation of the theory in the numerical scheme of geometrical shock dynamics allows the computation of shock motion in more general geometries. In particular, investigation of shock diffraction by convex corners of large angular deviation successfully yields the observed inflection point in the shock shape near the wall. The theory developed allows account to be taken of non-uniform flow conditions behind the shock. This feature is of particular interest in consideration of underwater blast waves in which case the flow behind the shock decays approximately exponentially. Application of the ideas developed here provides an excellent description of this phenomenon.

/ ~M (~(0,0 o~(t) l OM(~(t), t) n (~(t), t)) ~t 6~(t) + O~(t)., Ox(~q.), t)and equation (175) should readThe reference to the paper by Lighthill should read as set out below. 198:454-470