1. Introduction
A prevailing myth is that water feels like concrete if one jumps onto it from a great enough height. Although this may seem like an oversimplification, the statement is somewhat truthful. The impact force felt at the time of penetrating a quiescent water surface can be very high (Thompson Reference Thompson1928; Von Karman Reference Von Karman1929; Watanabe Reference Watanabe1933; Shiffman & Spencer Reference Shiffman and Spencer1945; May Reference May1975; Grady Reference Grady1979; Moghisi & Squire Reference Moghisi and Squire1981; Korobkin & Pukhnachov Reference Korobkin and Pukhnachov1988), much higher than the subsequent sustained underwater drag. One such example is shown in figure 1(a), where a 50 mm sphere dropped from 0.72 m above the free surface results in an impulse with a peak impact acceleration of 
$\sim$8 g whereas the underwater acceleration is close to a constant value of 
$\sim$2 g (1c), indicating that at the moment of surface penetration the drag coefficient 
$C_d$ is four times the steady state underwater free fall. This initially high impact force is primarily due to the large rate of change of momentum of the added fluid mass (May Reference May1975; Wang, Lugni & Faltinsen Reference Wang, Lugni and Faltinsen2015; Wang, Faltinsen & Lugni Reference Wang, Faltinsen and Lugni2019), which is the highest during a submergence depth of 10 %–20 % of the radius for spheres (figure 1a) (Shiffman & Spencer Reference Shiffman and Spencer1945; Moghisi & Squire Reference Moghisi and Squire1981). Reducing this peak impact force is of significant interest because it presents structural failure risk to impinging bodies like aircraft landing on water, water landing spacecraft, underwater missiles, divers, base jumpers, etc. (Kornhauser Reference Kornhauser1964; May Reference May1975; Seddon & Moatamedi Reference Seddon and Moatamedi2006; Guillet et al. Reference Guillet, Mouchet, Belayachi, Fay, Colturi, Lundstam, Hosoi, Clanet and Cohen2020). Previous studies have shown that impact forces can be reduced not only through object geometry (McGehee, Hathaway & Vaughan Reference McGehee, Hathaway and Vaughan1959; Thompson Reference Thompson1965; Li & Sigimura Reference Li and Sigimura1967; May Reference May1970; Qi et al. Reference Qi, Feng, Xu, Zhang and Li2016; Sharker et al. Reference Sharker, Holekamp, Mansoor, Fish and Truscott2019; Güzel & Korkmaz Reference Güzel and Korkmaz2020), but also by modifying the near-surface region via, for example, aeration (Elhimer et al. Reference Elhimer, Jacques, Alaoui and Gabillet2017) or liquid jet-induced acceleration (Speirs et al. Reference Speirs, Belden, Pan, Holekamp, Badlissi, Jones and Truscott2019a). An interesting extension to the idea of free surface modification is to launch a precursory object to agitate the free surface before entry. Such a concept has been proposed in popular culture (e.g. Mythbusters, Hollywood movies), yet has not received careful scientific investigation.
Two different impact phenomena with their impact acceleration data plotted. (a) Free-surface water impact of a 50 mm Vero plastic sphere (
$U_o=3.76\ \textrm {m s}^{-1}$). (b) The same sphere impacting at the same speed onto a cavity already created by a 10 mm sphere (3.38 m s
$^{-1}$). (c) An accelerometer embedded in the Vero plastic sphere reveals peaks of 
$\sim$8 g for the case in (a) and 
$\sim$4.8 g for the case in (b). The reduction in peak initial impact acceleration between the two indicates that the pre-formed cavity reduces the impact force on the trailing sphere.
Here, we present the findings from an experiment investigating the consecutive water entry of two spheres, where the spheres are axially aligned and vertically separated (figure 1b). Only a few other studies have looked into multiple sphere entry. Yang et al. (Reference Yang, Vaikuntanathan, Terzis, Cheng, Weigand and Helmig2019) studied the water entry of a linear array of magnetic spheres focusing on the cavity shape of the array. Yun, Lyu & Wei (Reference Yun, Lyu and Wei2019) showed that oblique two sphere entry resulted in contact with the cavity wall and eventual collision of the two spheres, but neither study focuses on the hydrodynamic forces involved. We measure the forces of water entry through an accelerometer embedded in the upper sphere which provides time-resolved measurements of the impact force. The lower sphere hits the water and creates a cavity through which the upper sphere falls, which can result in a reduced impact force on the upper sphere. Figure 1(c) presents an example, where the upper sphere experiences a 
$\sim$40 % reduction in impact acceleration compared to the case where the same sphere impacts the quiescent free surface at the same velocity but without a cavity in front. We propose a modified version of the non-dimensional parameter called the ‘Matryoshka’ (
$Mt$) number (Hurd et al. Reference Hurd, Fanning, Pan, Mabey, Bodily, Hacking, Speirs and Truscott2015; Speirs et al. Reference Speirs, Pan, Belden and Truscott2018) based on the cavity characteristics and the vertical spacing between the two spheres, which allows us to build an experimental regime diagram correlating different cavity conditions with the upper sphere impact force reduction and results in the observation of four distinct classes of consecutive two-sphere water entry behaviour.
2. Experimental methods
Figure 2(a) illustrates the experimental set-up used for this study. Two spheres of diameter 
$d_1$ and 
$d_2$ are placed on two vertically separated axially aligned platforms held above a glass water tank filled to a height of 1.2 m. The platforms are kept parallel to the water surface with the help of a clamped string and pulley mechanism. When the platforms are released at the same moment, the two spheres free fall and impact the water surface in tandem. The lower sphere (
$d_1$) impacts the water surface first, creating a cavity. The trailing upper sphere impacts on different stages of the cavity in front depending on the spacing between the two spheres, and the time difference of impact 
${\rm \Delta} t$. The impact acceleration for the upper sphere is recorded by an accelerometer embedded in an inertial measurement unit (IMU) placed inside the upper sphere. The impact events are captured with a Photron SA-3 high-speed camera at 1000 f.p.s. The videos taken were used to calculate different important parameters such as the cavity opening diameter, pinch-off time and to differentiate between different modes of consecutive two-sphere entry. The spheres are kept at heights of 
$h_2$ and 
$h_1$ from the water surface, the spacing (
${\rm \Delta} h = (h_2-h_1)$) between the spheres varies from 0.07 to 1.24 m. For every two-sphere impact combination tried, a quiescent upper sphere drop from the same drop height 
$h_2$ was also done to measure the peak quiescent impact acceleration for comparison with the two-sphere peak impact acceleration.
A schematic diagram of the experimental set-up for consecutive two-sphere water entry is shown in (a). The spacing 
$(h_2-h_1) = {\rm \Delta} h$ is varied to attain different modes of two-sphere impact. (b) A quiescent drop where the upper sphere impacts the undisturbed water pool from height 
$h_2 = 1.27$ m and creates a cavity. (c) The data taken by the accelerometer inside the sphere show the acceleration at different stages of the quiescent sphere impact and entry, with the blue, green and red coloured dashed boxes (b) and lines (c) indicating the time of peak acceleration, steady state underwater drag and starting point of cavity pinch-off. (d) The measured peak impact accelerations experienced by the upper sphere for the consecutive two-sphere impacts studied herein, presented with respect to the variation of the time difference 
${\rm \Delta} t$ between impacts. Variation in 
${\rm \Delta} h$ results in different 
${\rm \Delta} t$. Different markers represent the different diameter (
$d_1$) lower spheres (that travel in front of the upper sphere and impact the quiescent water pool first) used in the study.
The upper sphere is a 3D-printed Vero plastic sphere with a fixed diameter of 
$d_2=50$ mm. The sphere consists of two separable parts for helping in removing the IMU between experiments. The two parts of the sphere are pressed together with the help of butyl-ethylene tape to seal off and prevent water from entering. The contact area of the two parts is located one third the diameter from the top of the sphere, minimizing its influence on water entry events. Weights are inserted in the upper sphere to make it bottom heavy resulting in a density of 2290 kg m
$^{-3}$ for the sphere as a whole. Five different diameter steel spheres (
$d_1=10\text {--}38$ mm, density 7800 kg m
$^{-3}$) are used as the lower sphere, sprayed with Cytonix WX-2100 coating to make them hydrophobic, resulting in a surface contact angle of 117
$^{\circ }$ and increased the roughness of the spheres to 
$R_z =50.2 \pm 21.4\ \mathrm {\mu }$m (95 % confidence). The Vero plastic upper sphere has a hydrophilic surface of wetting angle 
$\theta = 80 \pm 8^{\circ }$ and surface roughness 
$R_{z} = 7.2 \pm 1.2\ \mathrm {\mu }$m (95 % confidence) such that it forms a cavity during water entry (see Table 1 in Appendix for more information).
Parameter space for consecutive two-sphere water entry.
The upper sphere housed an IMU built in house, with two three-axis accelerometers, one gyroscope and one magnetometer embedded. The two accelerometers on-board were one low range and one high range. The low range accelerometer has a measurement range of 
$\pm$16 g, it is a MPU-9250 motion tracking device manufactured by Invensen Inc. The high range accelerometer is a chip called H3LIS331DL produced by STMicroelectronics, and was set to a maximum range of 
$\pm$100 g. Both would register data for any drop event, whenever possible the data from the low accelerometer are reported, because the high accelerometer is more prone to noise. Data from both are comparable for cases where acceleration values were within 
$\pm$16 g. The accelerometer sampling rate is limited to 1000 Hz. The root sum square of the acceleration values from the three axes is calculated and reported as the total acceleration.
Figure 2(b) shows a typical quiescent upper sphere impact event, with the acceleration output from the IMU shown in figure 2(c). The sphere impacts the free surface at 0 ms, and a sudden increment of acceleration is registered. The peak acceleration is reached soon after, shown by the black dot and the time marked by the blue dashed line. This impact pulse lasts fleetingly until 
$\sim$8 ms, after which the sphere travels downwards with an air cavity in its wake until 97 ms, when the cavity pinches-off and divides into two halves. The sphere travels downwards after pinch-off with an air bubble attached to it. The bubble experiences pressure perturbations, which are registered in the accelerometer as oscillations after 97 ms (red dashed box and line in figure 2b,c). The quiescent peak values are compared with the peak accelerations measured from consecutive two-sphere impact drops (figure 2d) to calculate the change in impact acceleration.
2.1. Uncertainty
The uncertainty relating to the measurements are calculated and the uncertainty bands in the figures represent 95 % confidence interval of the measurement (Coleman & Steele Reference Coleman and Steele2018). The method of calculating uncertainties for peak acceleration values and non-dimensional ‘Matryoshka’ number is explained in detail in appendix A.1. The uncertainty bands are placed on acceleration plots whenever possible.
3. Results and discussion
3.1. Scaling analysis and formulation of ‘Matryoshka’ number
When an object impacts a water pool it displaces some of the water with air and accelerates fluid downwards as the object falls through the pool, leaving an air filled cavity in its wake (Truscott, Epps & Belden Reference Truscott, Epps and Belden2014). Creating a cavity in front of an impacting object can be conducive to reducing its impact acceleration, as evident from figure 1(c). The state of the cavity over time indicates the local liquid flow field surrounding the cavity (Truscott, Epps & Techet Reference Truscott, Epps and Techet2012; Mansoor et al. Reference Mansoor, Marston, Vakarelski and Thoroddsen2014), which may help explain the change in impact acceleration for any trailing object. Thus, understanding cavity creation and evolution is paramount for determining why and how an air cavity may reduce impact force.
Objects with rough and hydrophobic surfaces almost always create cavities even at very low impact velocities (Duez et al. Reference Duez, Ybert, Clanet and Bocquet2007; Zhao, Chen & Wang Reference Zhao, Chen and Wang2014; Speirs et al. Reference Speirs, Mansoor, Belden and Truscott2019b). Figure 3 shows such a case where a 38 mm hydrophobic sphere creates an axisymmetric cavity at an impact velocity of 
$\sim$3.2 m s
$^{-1}$. The cavity elongates with the downward moving sphere, until the point when hydrostatic pressure forces the cavity to seal near the cavity mid-point at 92 ms. This sealing event is popularly referred to as ‘deep-seal’ pinch-off. After pinch-off the cavity divides into two parts, a pulsating air bubble attached to the downwards moving sphere and the upper bowl shaped distortion in the free surface retreating upwards creating a high-speed axisymmetric ‘Worthington’ jet (Worthington & Cole Reference Worthington and Cole1897; Gekle & Gordillo Reference Gekle and Gordillo2010), as seen in figure 2(a) from 92 to 115 ms. Different impact velocities and sphere sizes result in different cavity behaviours which can be classified by cavity seal type (Aristoff & Bush Reference Aristoff and Bush2009; Speirs et al. Reference Speirs, Mansoor, Belden and Truscott2019b). For example, the cavity shown in figure 3 is referred to as deep-seal cavity because of the characteristic mid depth deep-seal pinch-off.
A typical lower sphere impact on quiescent water surface without a trailing upper sphere. A 38 mm hydrophobic steel sphere impacts a quiescent water pool (
$U_o\sim 3.2$ m s
$^{-1}$) creating a subsurface air cavity (0–68 ms). The cavity elongates in time with the downwards moving sphere (30–92 ms) and finally collapses on itself in a deep-seal (92 ms) pinch-off. After the deep seal, the cavity is divided in two parts, a cavity bubble attached with the sphere moving downwards (97 ms) and the upper bowl shaped distortion of the free surface, which eventually create a Worthington jet (97–115 ms).
Increasing the impact velocity for the same sphere sizes results in surface seal (Aristoff & Bush Reference Aristoff and Bush2009; Mansoor et al. Reference Mansoor, Marston, Vakarelski and Thoroddsen2014; Speirs et al. Reference Speirs, Mansoor, Belden and Truscott2019b), which is distinguished by the splash crown sealing above the free surface, and the resulting detachment and pull away of the cavity below the free surface. In the context of consecutive two-sphere water entry, one might expect the pinch-off (or seal) event from the cavity of the first sphere to affect the dynamics of the trailing sphere. For two axially aligned, vertically separated spheres (upper sphere diameter 
$d_2$, lower sphere diameter 
$d_1$) as shown in figure 2(a), varying the spacing 
${\rm \Delta} h = (h_2-h_1)$ between the two spheres will result in the upper sphere interacting with the cavity either before or after pinch-off, which we anticipate will lead to different sphere–cavity interaction modes. Taking the pinch-off time as a characteristic time scale, we propose to characterize consecutive two-sphere water entry with a modification of the non-dimensional parameter known as the ‘Matryoshka’ number 
$Mt$. This term has been used in prior research to describe successive cavity formation from multi-droplet impacts using droplet frequency and cavity formation time as the fundamental time scales (Hurd et al. Reference Hurd, Fanning, Pan, Mabey, Bodily, Hacking, Speirs and Truscott2015; Speirs et al. Reference Speirs, Pan, Belden and Truscott2018).In a physical sense, 
$Mt$ can be considered a ratio of the time to completion of a single event to the consecutive initiation of the same event by the second sphere. Here, we define
\begin{equation} Mt = \frac{{\rm \Delta} t}{t_p}. \end{equation}
where 
${\rm \Delta} t = |t_2 - t_1|$ is the time difference between the two spheres passing the free surface (figure 4a), and 
$t_p$ is the pinch-off time of the first cavity. Thus, 
$Mt$ parameterizes the state of the cavity formed by the first sphere at the time when the second sphere interacts with it. The case 
$Mt<1$ indicates that the first cavity has not gone through pinch-off, which results in cases where the upper sphere interacts with an elongating cavity. The cavity opening diameter 
$d_c$ (itself a function of time and the lower sphere diameter Duclaux et al. Reference Duclaux, Caillé, Duez, Ybert, Bocquet and Clanet2007; Aristoff & Bush Reference Aristoff and Bush2009) is measured at time 
$t_2$. If 
$d_c < d_2$, the upper sphere falls on the cavity when impacting the water pool, and we name these ‘on cavity’ cases (figure 4b). When 
$d_c > d_2$, then the upper sphere falls through the cavity opening, which we call the ‘inside cavity’ case, with 
$d_c/d_2 = 1$ working as the transition between the two cavity cases predicted for 
$Mt < 1$. For 
$Mt>1$, the trailing sphere interacts with the upper detached portion of the cavity, either falling through (
$d_c > d_2$) or falling on (
$d_c < d_2$) the bowl shaped retreating free surface; or falling through a Worthington jet resulting from cavity pinch-off at higher 
$Mt$ (figure 4b). We call these cases ‘bowl’ and ‘on jet’, respectively. Figure 5(a) shows snapshots of these predicted modes of consecutive two-sphere water entry from actual experiments, with the modes observed at various 
$Mt$. The experimental map of the consecutive two-sphere entry cases tested herein (figure 5b) illustrates that the ratio 
$d_c/d_2=1$ does work as a separation criterion between the on cavity and inside cavity cases. Note that figure 5(b) also denotes the existence of two different on bowl modes (on cavity bowl and inside cavity bowl as expected from figure 4b). These two can also mostly be separated using the 
$d_c/d_2 =1$ condition. The transition from on cavity to bowl mode happens at 
$Mt=1$, as predicted, denoted by the magenta dashed line in figure 5(a), but the transition from inside cavity to inside bowl cases happens earlier (
$Mt=0.65$, vertical black dashed line on figure 5b). This discrepancy in transition 
$Mt$ between the two pre-pinch-off cavity modes is surprising, but can be explained analytically, as shown in later discussion in § 3.2.
The formulation of the non-dimensional ‘Matryoshka’ number and predicting different modes of two-sphere water entry using cavity pinch-off time information. (a) When two spheres vertically separated by a distance 
${\rm \Delta} h$ are released from their rest position as shown in figure 2(a), the variation in spacing will result in the upper sphere impacting various stages of the cavity created by the lower sphere. Taking the characteristic deep-seal pinch-off time as 
$t_p$, a non-dimensional time parameter called ‘Matryoshka’ number can be defined. (b) The case 
$Mt = 1$ predicts the boundary between different modes of consecutive two-sphere impact. For 
$Mt<1$, the upper sphere can interact with a still growing cavity, yielding pre-pinch-off cavity cases: ‘on cavity’ and ‘inside cavity’. Beyond 
$Mt>1$, the upper sphere should interact with the collapsed cavity, resulting either in two different ‘bowl’ modes based on the cavity opening diameter 
$d_c$, or with the Worthington jet (‘on jet’ cases).
(a) Snapshots of five different two-sphere entry modes observed through experiments. (b) The experimental map of 
$Mt$ and ratio of cavity to the upper sphere diameter 
$d_c/d_2$. The different coloured markers denote different two-sphere entry modes, with red, blue and green denoting inside cavity, on cavity and on jet cases; black and magenta represent cases belonging to the two different bowl modes: inside cavity bowl and on cavity bowl, respectively. The dotted yellow line at 
$d_c/d_2 =1$ marks the separation between the inside cavity and on cavity cases, as shown in figure 4(b). The transitions between the modes are represented by the vertical dashed black line at 
$Mt = 0.65$ between pre-pinch-off inside cavity and post-pinch-off inside bowl modes, the vertical dashed magenta line at 
$Mt = 1$ between on cavity and on cavity bowl cases, and the vertical dashed green line at 
$Mt = 1.43$ between post-pinch-off bowl and on jet modes. The emergence of these values as separation criteria is explained in § 3.2. (c) The experimental value of 
$\beta$ for different sphere diameters over the Froude (
$Fr$) number range tried herein for 169 different single-sphere drop experiments. The black dashed line shows the averaged value of 
$\beta$ to be 
$2.037\pm 0.097$ (95 % confidence with vertical error bar).
To determine 
$Mt$ as formulated in (3.1) requires experimental measurements of 
${\rm \Delta} t$ and 
$t_p$. To make our results more broadly useful to other researchers, we seek to redefine 
$Mt$ in terms of independent variables. Prior research has shown that the deep-seal pinch-off time can be written as 
$t_p = \beta \sqrt {d_s/(2g)}$, (Glasheen & McMahon Reference Glasheen and McMahon1996; Duclaux et al. Reference Duclaux, Caillé, Duez, Ybert, Bocquet and Clanet2007; Truscott & Techet Reference Truscott and Techet2009), where 
$\beta$ is an experimental constant with different values ranging from 1.72 to 2.285 proposed in the literature (Duclaux et al. Reference Duclaux, Caillé, Duez, Ybert, Bocquet and Clanet2007; Bergmann et al. Reference Bergmann, Van Der Meer, Gekle, Van Der Bos and Lohse2009; Marston, Vakarelski & Thoroddsen Reference Marston, Vakarelski and Thoroddsen2012). This equation works well for low Froude number (
$Fr$) impacts, where 
$Fr = U_o^{2}/(gd_s)$, 
$U_o$ is the impact velocity, 
$d_s$ is the sphere diameter and 
$g$ is the gravitational constant. In this study, experimentally calculated 
$\beta = 2.03\pm 0.0974$ (95 % confidence) is used (figure 5c). For higher 
$Fr$ where surface seals are expected to happen, this equation overestimates the pinch-off time marginally (figure 16a), but the lack of a good consensus in the literature about the surface seal time and the scale of the experiments carried out in this paper (
$Fr<400$) makes this equation a good approximation. The time spacing 
${\rm \Delta} t$ can be calculated in terms of sphere drop heights using the free-fall relationship between the two spheres. If the upper sphere takes time 
${\rm \Delta} t = (t_2-t_1)$ to pass the free-surface line, then 
${\rm \Delta} t = (\sqrt {2(h_1+{\rm \Delta} h)/g}-\sqrt {2(h_1)/g})$. Substituting this in (3.1) and using the pinch-off time relation 
$t_p = \beta \sqrt {d_1/2g}$ yields,
\begin{equation} Mt = \frac{2(-\sqrt{h_1}+\sqrt{h_1+{\rm \Delta} h})}{\beta \sqrt{d_1}} , \end{equation}which is as an alternative interpretation for the definition of the Mt number. Not only is (3.2) easier to use for prediction than (3.1), but it also explains how the lower sphere diameter and the height difference between the two spheres are required to adequately explain the phenomenon (figure 16b).
3.2. Consecutive sphere entry: different modes
As explained in the previous section, 
$Mt<1$ results in pre-pinch-off cavity cases (on cavity and inside cavity) depending on the ratio 
$d_{c}/d_{2}$ (figures 4(b), 5(b)), and 
$Mt>1$ indicates the post-pinch-off cases (bowl and on jet). Figures 6 and 7 present time-series image sequences of all two-sphere modes including two different bowl cases (figure 7a,b), with their dynamic acceleration response plotted with the image sequences (figures 6(c), 7(d)). The accelerations at impact are reduced for both on cavity and inside cavity cases compared to the quiescent case peak. In the bowl cases, the initial impact pulse has a higher peak than the quiescent case, indicating a higher impact force experienced during free-surface entry. For the on jet case, the peak acceleration value is significantly smaller than the quiescent peak. The time of the peak accelerations for the two-sphere cases happen later than the quiescent case (see figures 6(a–c), 7(a–d)), since the sphere interacts with modified free-surface conditions.
Time series for two example pre-pinch-off cavity cases in consecutive sphere entry. A hydrophobic sphere (
$d_1 = 10$ mm (a) and 38 mm (b)) impacts the quiescent free surface and creates a cavity, followed by a 50 mm sphere with embedded IMU in free fall from a height of 
$h_2=0.72$ m, where 
${\rm \Delta} h= 0.14$ and 0.11 respectively for (a) and (b). Acceleration of both cases is shown in (c). The time the upper sphere passes the free-surface line is considered as 
$t = 0$ ms for each of the cases, yielding 
$Mt\sim 0.82$ and 
${\sim }0.35$ for the on cavity (a) and the inside cavity (b) cases, respectively. The blue and red dashed boxes and lines in (a–c) denote the time of peak acceleration felt by the upper sphere for on cavity and inside cavity cases. The peak acceleration in both cases are significantly less than the quiescent value of 8 g. The uncertainty bands are as shown. See supplemental movies S1 and S2 available at https://doi.org/10.1017/jfm.2020.1165.
A hydrophobic sphere (
$d_1 = 10$ mm (a) and 38 mm (b,c)) impacts the quiescent free surface and creates a cavity, followed by the 50 mm upper sphere in free fall from a height of 
$h_2=0.72$ m, where the cases are: (a) inside cavity bowl is 
${\rm \Delta} h=0.22$ and 
$Mt\sim 1.26$, (b) on cavity bowl is 
${\rm \Delta} h= 0.3$ and 
$Mt\sim 1.05$ and (c) on jet is 
${\rm \Delta} h= 0.54$ and 
$Mt\sim 2.20$. (d) Acceleration of each of the cases in (a–c) with peak acceleration time for all marked with magenta, black and green dashed boxes (a–c) and lines (d). The black dashed-dotted horizontal line at 
$a\sim 8$ g shows the peak acceleration value for the quiescent case where the upper sphere dropped from the same height (0.72 m) impacts a quiescent free surface without a cavity in front of it (figure 6c). The bowl cases show significant increase in acceleration than the quiescent value of 8 g, but the on jet case shows dramatic reduction in peak impact acceleration, which never rises beyond the level of free-fall acceleration through water. See supplemental movies S3–S5.
Figure 8 shows a regime diagram where reduction in acceleration for varying 
$Mt$ is plotted for the range of experimental conditions. Reduction in acceleration is computed as 
$1-a/a_q$, where 
$a$ is the measured peak acceleration of the trailing sphere in a two-sphere water entry, and 
$a_q$ is the peak acceleration of the same sphere impacting quiescent water from the same drop height (
$h_2$). Similar to the standalone cases presented in figures 6 and 7, on cavity (
$0 < Mt < 1$) and inside cavity cases (
$0 < Mt < 0.65$) experience a notable reduction in impact acceleration, with a downward linear trend in reduction values present for both as 
$Mt$ increases. Alternatively, the bowl cases experience an increase in impact acceleration, evident from the negative reduction in acceleration in figure 8 up to 
$Mt\sim 1.43$. The on jet cases experience dramatic reduction in impact acceleration in 
$1.43 < Mt < 4$ range, with the mode extending to 
$Mt\sim 6$.
Percentage reduction in the initial impact force for the upper sphere versus Mt for different two-sphere water entry modes. Here, 
$a$ is the peak impact acceleration for the second of the two-sphere impact cases, and 
$a_q$ is the corresponding quiescent impact peak acceleration of the same sphere dropped from the same height. Cavity impact types are the same as figure 4. Significant force reductions are achieved in the inside cavity (
$0 < Mt < 0.65$), on cavity (
$0 < Mt < 1$) and on jet (
$1.43 < Mt < 4$) cases. The transition from inside cavity to on bowl occurs at 
$Mt\sim 0.65$ (black dashed line) and on cavity to on bowl mode at 
$Mt\sim 1$ (magenta dashed line), with both bowl modes showing a significant increase in impact acceleration. The inside cavity bowl and on cavity bowl cases transition to on jet cases at 
$Mt\sim 1.43$ indicated by the dashed vertical green line. A variation of the on jet cases occurs when the splash domeover at the free surface suppresses the Worthington jet, creating a small jet-like water column at the surface we term the ‘on surface jet’, marked with cyan coloured markers. The horizontal black dashed-dot line at 
$1-a/a_q = 0$ in the 
$y$-axis indicates 
$a = a_q$, while the vertical green dashed line represents the experimentally found transition from both on bowl modes to both on jet modes. Uncertainty bands are marked with 95 % confidence with more information available in appendix A.1. Note: the 
$y$-axis values above and below 0 are scaled differently.
We can gain a mechanistic understanding of how the peak accelerations are reduced for some regimes (on cavity, inside cavity, on jet) and dramatically increased for others (on bowl, inside bowl) by considering the contributions to the large force during the first few moments of impact. For inertially dominated regimes, this peak force stems from the rate of change of momentum of the added fluid mass (Shiffman & Spencer Reference Shiffman and Spencer1945; May Reference May1975). For impact of a sphere on a quiescent free surface, Shiffman & Spencer (Reference Shiffman and Spencer1945) derived an expression for this force invoking conservation of momentum during water impact such that 
$MU_o = (M+m)U$, where 
$M$ is the sphere mass, 
$m$ is the added mass of the fluid, 
$U_o$ the impact velocity and 
$U$ the instantaneous velocity of the sphere and added fluid mass after impact. Defining 
$b = s/R$, where 
$R$ is the sphere radius and 
$s$ the submergence below the initial free-surface plane, the vertical force on the sphere is given as
\begin{equation} F_z ={-}M\frac{\textrm{d} U}{\textrm{d} t} = \frac{U_o^2}{R}\frac{1}{ \left(1+\dfrac{m}{M}\right)^3}\frac{\textrm{d} m}{\textrm{d} b} \end{equation}
(see appendix A.2 for detailed derivation). Using the typical definition of the vertical force 
$F_z = \rho {\rm \pi}R^2 U_o^2 C_d/2$, where 
$C_d$ is the impact drag coefficient, we get,
\begin{equation} C_d = \frac{2}{\rho{\rm \pi} R^3}\frac{1}{ \left(1+\dfrac{m}{M}\right)^3}\frac{\textrm{d} m}{\textrm{d} b},\end{equation}
which Shiffman & Spencer (Reference Shiffman and Spencer1945) report for values of 
$m/M$. Thus, the vertical force during early water entry has a first-order dependence on the rate of change of added mass with submergence, 
$\textrm {d} m/\textrm {d} b$. In the on cavity and inside cavity cases, the leading sphere opens a hole on the surface such that the first contact between the water and the sphere is below the free surface and at some higher azimuthal angle on the sphere (see schematic in figure 9). In the quiescent sphere impact, the value of 
$\textrm {d} m/\textrm {d} b$ peaks at small submergence of 
$b\approx 0.1$ (Shiffman & Spencer Reference Shiffman and Spencer1945; May Reference May1975). Therefore, one might expect that for on cavity or inside cavity cases for which the bottom of the sphere does not initially contact the water, the value of 
$\textrm {d} m/\textrm {d} b$ and thus the impact force would be smaller than in the quiescent case, as is observed. While the on cavity or inside cavity cases occur prior to first cavity pinch-off, for the on bowl and inside bowl cases the cavity has pinched off and the fluid has significant upward momentum as the bowl retracts. We can rationalize the effect of this to first order by again considering conservation of momentum before and after impact 
$MU_o - mU_m= (M+m)U$, where 
$U_m$ is the average velocity of the added fluid mass just before impact (assumed to be in the direction opposite of 
$U_o$). Following the same derivation for impact on a quiescent free surface from Shiffman & Spencer (Reference Shiffman and Spencer1945), we arrive at the expression,
\begin{equation} F_z = \frac{\left(U_o-\dfrac{m}{M}U_m\right)}{\left(1+\dfrac{m}{M}\right)^3} \frac{\left(U_o+U_m\right)}{R}\frac{\textrm{d} m}{\textrm{d} b}\end{equation}
(see derivation in appendix A.2). Now, the vertical force depends not only on 
$\textrm {d} m/\textrm {d} b$, but also increases with 
$U_m$ assuming 
$m/M$ is sufficiently small. Thus, one would expect the force to increase relative to the quiescent case due to the upward momentum of the retracting bowl that must be reversed during impact. Furthermore, in most on bowl and inside bowl cases, the shape of the surface at the moment of water impact is contoured to the sphere rather than flat (see schematic in figure 9). Therefore, the rate of change of added mass with submergence, 
$\textrm {d} m/\textrm {d} b$, would be larger than impact of a sphere with a flat free surface. Intuitively, one can think of a disk impacting a flat free surface, for which 
$\textrm {d} m/\textrm {d} b$ would go to infinity if the surface were not first deformed by an over-pressurized air layer (May Reference May1975). These two effects – the upward momentum of the fluid prior to impact, and larger rate of change of added mass during impact – combine to make the impact force larger for on bowl and inside bowl cases relative to the quiescent case.
The interaction of the trailing sphere with different stages of the cavity in front leads to varying impact force experienced by the trailing sphere. The impact force depends on the rate of change of added mass 
$\textrm {d} m/\textrm {d} b$, and the upward velocity 
$U_m$ of the added fluid mass. Notional increase/decrease of these values relative to the quiescent case are shown for all but the on jet case (the physics behind which is explained in later discussions).
While the trends predicted by this first-order analysis explain the observations of force increase for the bowl cases and force reduction for on and inside cavity cases, the experimental data show a wide range in the peak acceleration in these regimes (see figure 8). The impact acceleration increases with increasing 
$Mt$ for inside cavity cases, which can be understood qualitatively by considering the cavity evolution relative to the second sphere impact. For small 
$Mt$, the cavity is open and cylinder-like at the moment of impact, resulting in water contact near the equator of the sphere where 
$\textrm {d} m/\textrm {d} b$ is small. As 
$Mt$ increases, sphere impact happens closer to the time of pinch-off and the cavity is more conical, resulting in sphere wetting closer to the bottom of the sphere where 
$\textrm {d} m/\textrm {d} b$ is larger. It is also possible that the fluid has some upward momentum in the vicinity of the impact location as the cavity evolves toward pinch-off (Truscott et al. Reference Truscott, Epps and Techet2012; Mansoor et al. Reference Mansoor, Marston, Vakarelski and Thoroddsen2014), contributing to larger fluid momentum that must be reversed during impact. For the on cavity regime, the cavity diameter at the surface does not vary much in time. Therefore, it seems unlikely that a significant difference in 
$\textrm {d} m/\textrm {d} b$ exists between the on cavity cases. A possible source of the increase in acceleration is the fluid momentum in the vicinity of the cavity opening. For small 
$Mt$, the cavity is expanding and the added fluid momentum is outward. With increasing 
$Mt$, the added fluid momentum becomes neutral as the cavity stops expanding, and thus the second sphere needs to provide a larger change of momentum for the added fluid mass. In general, as 
$Mt$ increases and the two-sphere impacts transition between the regimes depicted in figure 9, the relative contributions from 
$\textrm {d} m/\textrm {d} b$ and 
$U_m$ change. Detailed modelling of the sphere impact with the various cavity states is beyond the scope of this work, but the framework presented here lends insight into how the peak force changes with 
$Mt$ for the observed two-sphere regimes.
At first glance, we would expect both on cavity and inside cavity cases to transition to on bowl and inside bowl cases respectively when 
$Mt\sim 1$. This is true for the on-cavity-to-on-bowl transition. Yet, when 
$d_c > d_2$ and 
$Mt<1$, pinch-off can occur while the upper sphere falls through the upper portion of the cavity leading to an in bowl impact. This reduces the 
$Mt$ value for the in-cavity-to-in-bowl transition, which we can predict by taking into account the time 
$t_{HP}$ for the upper sphere to fall from the original free-surface level to the pinch-off point. We find that the limiting pinch-off time for this to occur is,
\begin{equation} t_p = (t_2-t_1)+t_{HP}, \end{equation}
Dividing both sides of (3.6) by 
$t_p$ and rearranging we get,
\begin{equation} Mt_c = 1-\frac{t_{HP}}{t_p}. \end{equation}
From (3.7) we see that 
$Mt_c$ is the critical value where in-cavity-to-in bowl transition would occur, and 
$Mt_c$ is smaller than 1. The time 
$t_{HP}$ can be estimated using the relation 
$t_{HP} = (u_{2p}-u_2)/{g}$, where 
$u_{2p}$ is the velocity with which the upper sphere impacts the pinch-off point and 
$u_2$ is the initial velocity with which the upper sphere passes the free surface (
$u_2 = \sqrt {2gh_2}$). The value 
$u_{2p}$ can be calculated from the relation 
$u_{2p} = \sqrt {u_2^2+2gH_p}$, where 
$H_p$ is the pinch-off depth; 
$H_p$ can be estimated using the empirical relation 
$H_p = 0.57d_{1}\sqrt {Fr}$ (Duclaux et al. Reference Duclaux, Caillé, Duez, Ybert, Bocquet and Clanet2007), where 
$d_1$ is the lower sphere diameter. Then, (3.7) can be rewritten as,
\begin{equation} Mt_c = 1-\frac{\sqrt{u_2^2+2gH_p}-u_2}{gt_p} = 1-\frac{2(\sqrt{h_2+0.57d_{1}Fr^{1/2}}-\sqrt{h_2})}{\beta\sqrt{d_1}}, \end{equation}
where 
$t_p$ is written as 
$t_p = \beta \sqrt {d_1/(2g)}$. Equation (3.8) reveals the transition 
$Mt_c$ at which the upper sphere falls onto the pinch-off singularity from independent parameters 
$h_1,d_1$ and 
$h_2$. Using these parameters from 33 different two-sphere cases, this 
$Mt_c$ is determined to be 
$0.6539\pm 0.022$ (95 % confidence). Plotting this line as the separation criterion for inside cavity and bowl cases in figure 8 illustrates that the experimental data agree with this 
$Mt_c = 0.65$ separation line quite well.
Depending on the cavity opening diameter 
$d_c$, two different on bowl modes are possible: on cavity bowl (
$d_c < d_2$) and inside cavity bowl (
$d_c>d_2$). They transition from on cavity and inside cavity modes at 
$Mt\sim 1$ and 
$\sim 0.65$ respectively as discussed earlier. In both, as discussed with the help of (3.5), the coupled effects of upward added mass velocity 
$U_m$ and the contoured curvature of the bowl surface increasing the 
$\textrm {d} m/\textrm {d} b$ value leads to dramatically higher peak impact accelerations than a normal quiescent case (figure 7(a), 0–9 ms, 7(b), 10–20 ms). Figure 8 shows that the bowl cases almost always have higher impact acceleration; in some instances impact forces almost quadruple that of the quiescent case value (
${\approx }427\,\%$ at 
$Mt\sim 0.96$). Thus, this 
$0.65 < Mt < 1.43$ range of Matryoshka values (
$Mt\sim 1.43$ being the transition from bowl to jet mode, discussed in the next paragraph) is a range to avoid if one wishes to achieve any sort of reduction in impact force. The inside cavity bowl cases show much higher accelerations than their on cavity bowl counterparts, which can be attributed to the upward region of the bowl being much larger for inside cavity bowl cases and hence the average 
$U_{m}$ and 
$\textrm {d} m/\textrm {d} b$ are greater over the frontal area of the sphere.
The bowl eventually forms an upward moving axial jet called a Worthington jet (Worthington & Cole Reference Worthington and Cole1897; Worthington Reference Worthington1908; Gekle & Gordillo Reference Gekle and Gordillo2010) coming out of the base of the distorted bowl-shaped free surface following the collapse of the cavity after pinch-off (
$Mt\geq 1.43)$. Here, the on jet cases are characterized by the upper sphere dramatically passing through the Worthington jet into the water pool (figure 7c). Finding the 
$Mt$ at which the on bowl cases would transition to on jet cases is a complicated task, since theoretical consideration of the bowl shaped free-surface distortion that results in the Worthington jet is scarce in the literature. To estimate the 
$Mt$ where we could expect the transition from bowl to jet mode to happen, we experimentally observed 20 different quiescent drops that result in Worthington jets and calculated the time 
$t_j$ after pinch-off (
$t_p$) for the underwater bowl to transform into a fully formed jet above the free surface, assuming the jet to be fully formed when the upwards moving sub-surface bowl has fully vanished. If the upper sphere passes the free-surface line at the exact moment when the bowl vanishes, we can consider the upper sphere free-surface passing time 
$t_2$ as the onset of on jet cases. Thus, we can write
\begin{equation} (t_2-t_1) = t_p+t_j \end{equation}
to define the transition from bowl to on jet (
$t_1$ is the time of the lower sphere impacting the free surface). Now dividing both sides of (3.9) by 
$t_p$, we end up with,
\begin{equation} Mt = 1+\frac{t_j}{t_p} . \end{equation}
Thus (3.10) finds the 
$Mt$ value of the onset of the on jet mode. From the 20 observations, the fractional term 
$t_j/t_p$ is found to be a constant value of 
$0.43\pm 0.06$ for varied 
$Fr$ impacts, resulting in the 
$Mt$ value of 
$1.4345\pm 0.06$ (standard uncertainty, 95 % confidence, see appendix A.3, figure 17). The value of 
$Mt \sim 1.43$ can be used as the separation criterion between on bowl and on jet modes, and the experimental data presented in figure 8 follow this separation reasonably well.
The acceleration reduction for these 
$Mt\geq 1.43$ cases are also dramatic, with the highest reduction of up to 
$\sim$78 % observed experimentally (
$Mt=2.67$, figure 8). The time steps and the acceleration plot presented in figure 7(c) provide an explanation of how the reduction in impact acceleration manifests. The jet starts wetting the trailing upper sphere long before it has reached the free surface (80 ms before free-surface impact, not shown in figure 7(d) on jet case, see figure 18). This drawn out collision with the narrow axial jet results in a reduction in momentum of the upper sphere over a longer period of time, and also results in partial wetting of the upper sphere by the jet. At the time of impact, the water-enveloped upper sphere does not abruptly go through an air–water interface like the quiescent case does. Instead, it enters the water pool through a thick jet base and is enveloped by the jet to such extent that it resembles a sphere inside a thick jet entering the pool. Speirs et al. (Reference Speirs, Belden, Pan, Holekamp, Badlissi, Jones and Truscott2019a) showed that, for a similar kind of water entry of a sphere inside a free-falling jet, the maximum impact acceleration is a function of 
$Fr$, and can be modelled using a modification to the force balance equation,
\begin{equation} \frac{a_q}{g} = \frac{3}{4}{\frac{\rho}{\rho_s}C_{d_{max}}Fr} ,\end{equation}
where 
$\rho /\rho _s$ is the density ratio of water and the sphere (
${\sim }2.29$ in our study), and 
$C_{d_{max}}$ is the impact drag coefficient. Shiffman & Spencer (Reference Shiffman and Spencer1945) showed that, for large density ratios, the value of 
$C_{d_{max}}$ can be taken as 
$\sim$1 during the initial moment of entry (when 
$b\sim 0.1$–0.2 per discussion of (3.3)) of a single sphere into a quiescent pool. Speirs et al. (Reference Speirs, Belden, Pan, Holekamp, Badlissi, Jones and Truscott2019a) argued that when a sphere enters the pool travelling within a jet, this 
$C_{d_{max}}\sim 1$ value decreases, and (3.11) becomes,
\begin{equation} \frac{a}{g} = \frac{3}{4}{\frac{\rho}{\rho_s}(\alpha C_{d_{max}})Fr}, \end{equation}
where they reported 
$\alpha =0.25$, reducing the impact drag coefficient in sphere in a jet cases to 
$C_{d}^{new} = 0.25C_{d_{max}}$. Despite the opposing jet direction between Speirs et al. (Reference Speirs, Belden, Pan, Holekamp, Badlissi, Jones and Truscott2019a) and the Worthington jets in this study, plotting the measured impact acceleration values as a function of 
$Fr$ in figure 10 reveals a similar reduction in impact drag coefficient for the on jet cases, albeit only in the 
$1.43\leq Mt<4$ range. We experimentally observe 
$\alpha \sim 0.4$, yielding 
$C_{d}^{new} = 0.4C_{d_{max}}$ (figure 10a,b) for 
$Mt<4$ on jet cases, with the Froude of impact calculated from 
$Fr = {u_{imp}^2/(gd_2)}$, where 
$u_{imp}$ is the reduced impact velocity resulting from the sustained collision of the upper sphere with the Worthington jet (inset, figure 10a). At 
$Mt>4$, the on jet 
$C_d$ values are closer to one (figure 10b), which can be attributed to the retreating nature of the Worthington jet as 
$Mt$ increases. A greater 
$Mt$ indicates a larger time difference between the jet creation and the upper sphere passing the free surface. Whereas thick jet diameters and prominent jet bases are typical of 
$Mt<4$ cases (figure 11a), as the time between the two-sphere free-surface impact is increased, the Worthington jet changes from an upward rising jet to a descending jet, with the jet peak starting to thin, and the jet base eventually vanishing beyond 
$Mt \sim 4$ (figure 11b).
The maximum impact accelerations (a) and the associated impact drag coefficients 
$C_d$ (b) in on jet cases. Prolonged collision with the upward moving jet results in a reduced trailing sphere impact velocity 
$u_{imp}$, shown for several cases in the inset in (a). The dashed line in (a) shows the theoretical impact acceleration using this reduced 
$u_{imp}$ in (3.11), and the dotted line represents (3.12) with 
$\alpha = 0.4$. The squares plot actual measured impact accelerations in different on jet cases, the green squares are for 
$Mt<4$ cases, and the red squares are for 
$Mt >4$. In (b), the impact drag coefficients are plotted for both 
$Mt<4$ and 
$Mt>4$, with the dashed line showing quiescent impact drag coefficient of 
$1$, and the dash-dotted line representing 
$C_d^{new} \sim 0.4C_d$. Uncertainty bands are marked with 95 % confidence interval.
Two different on jet cases happening at either side of 
$Mt\sim 4$ resulting in contrasting impact accelerations. In (a), 
$Mt = 2.67$, and the upper sphere falls through a fully formed Worthington jet. The lower sphere diameter is 38 mm. The crash through the jet is violent, the sphere gets entirely enveloped by the jet and falls through a thick jet base with a reduced impact velocity of 4.36 m s
$^{-1}$. The acceleration reduction here is prominent, 
$\sim$78 % for this particular case. (b) An on jet case of 
$Mt = 4.50$. The lower sphere (
$d_1\sim 23.75$ mm) creates a cavity that evolves into the Worthington jet. Larger 
$Mt$ results in a thinner jet and a significantly calmer impact event when the trailing sphere passes the free surface with a reduced impact velocity 4.89 m s
$^{-1}$. The jet base is non-existent, the sphere falls through a retracting water column. The reduction in impact acceleration here is a mere 
$\sim$7 %. There is stark difference in jet diameter 
$d_j$ at the jet base also; 
$d_j$ is much bigger in (a) than in (b).
We can understand how the jet base diameter affects impact acceleration by again considering the rate of change of added mass 
$dm/db$. As the sphere moves through the jet, some amount of the jet fluid passes around the sphere, but presumably some amount is also accelerated downward by the sphere. As the sphere descends, the jet diameter increases smoothly to the value 
$d_j$ (measured at the top of the jet base when the sphere is one diameter above the base as shown in figure 11). One might expect the fluid below the flat water surface to be accelerated by the downward moving jet fluid prior to the sphere entering the water pool. This effect would reduce the value of 
$\textrm {d} m/\textrm {d} b$ during water entry of the sphere, thus reducing the impact force and resultant acceleration. As 
$d_j \longrightarrow d_2$ and a larger amount of fluid mass exists in the jet base, 
$\textrm {d} m/\textrm {d} b$ would be smaller at the moment of water entry, and one would expect the impact acceleration reduction to be greater. This notion is supported by the data in figure 12, which show acceleration reduction increasing with increasing 
$d_j/d_2$.
The reduction in impact acceleration in on jet cases are linearly dependent on the jet base diameter prior to the upper sphere passing the free-surface line. Increasing jet diameter results in higher reduction, with the jet diameter monotonically decreasing over increasing 
$Mt$, as shown in the inset. The different coloured markers follow the same convention presented in figure 10. The blue dashed line is fitted through least-square regression giving the equation: 
$1-a/a_q = 1.89(d_j/d_2)-0.1306$.
When the second sphere falls through a thick Worthington jet a chaotic but pretty watery mess can occur, as shown in figure 13. As the second sphere falls through the jet, it forces it into a thin film, prompting intriguing fingering patterns, buckling and droplets at the top of the sphere as the thin film envelops the sphere over the apex. A time series showing this violent yet beautiful splash phenomenon can also be seen in figure 11(a).
A 50 mm steel sphere falling through the Worthington jet originating from the collapse of a cavity created by a 38 mm hydrophobic steel sphere impacting a quiescent water pool at 2.09 m s
$^{-1}$. The beautiful fluid patterns observed at the surface of the sphere originate from the jet water films following the sphere curvature and crashing together at the top of the sphere.
In some of the on jet cases, the Worthington jet created from pinch-off is suppressed by the splash crown domeover at the free surface. This results in a disturbed, jet-like water column at the free surface, through which the upper sphere falls. Here, we consider these cases as ‘on surface jet’ cases, since their formation mechanism is different from that of the on jet cases (see appendix A.4, figure 19). These cases show reduction in acceleration similar to the on jet cases as illustrated in figure 8.
4. Conclusion
The impulse force felt by any object at the initial moment of water impact can be very high (figure 1c), and may prove to be catastrophic for water landing craft or missiles, and fatal even for thrill seeking bungee jumpers (Von Karman Reference Von Karman1929; McGehee et al. Reference McGehee, Hathaway and Vaughan1959; Kornhauser Reference Kornhauser1964; Thompson Reference Thompson1965). Herein, we have shown through a canonical sphere impact study that the initial impact force can be greatly reduced by first launching another object in front of the body of interest. The force of impact is reduced by the cavity of the first object providing less initial water impact area, lower relative velocities or upward jets that wet and decelerate the trailing body. However, if the object encounters the collapsing upward cavity, the upper sphere may experience a larger force of impact than if the leading object were not present at all. A non-dimensional number called the Matryoshka number 
$Mt$ is defined to classify two-sphere consecutive water entry behaviour based on the object size and the cavity pinch-off time. Experimental results show that for 
$0.2< Mt < 0.65$ and 
$1.43 < Mt < 4$, significant reduction of the impact acceleration of the trailing sphere is achieved. In the interim range of 
$0.65 < Mt < 1.43$, a sudden rise in impact acceleration is observed, and must be avoided if trying to avoid catastrophic failures. This 
$Mt$ formulation can potentially be used to predict interactions in any multi-object water entry system, and the regime diagram proposed in figure 8 does well to aid in making predictions for size differences and timing.
Supplementary movies
Supplementary movies are available at https://doi.org/10.1017/jfm.2020.1165.
Funding
R.R., N.S., J.B. and T.T.T. acknowledge funding from the Office of Naval Research, Navy Undersea Research Program (Grant no. N000141812334), monitored by Ms M. Medeiros. A.K. is a JSPS Overseas Research Fellow.
Declaration of interests
The authors report no conflict of interest.
Appendix A
A.1. Uncertainty analysis
To calculate the uncertainty relating to the variables measured for our two-sphere study we first look at our reported results. From figures 6(c), 7(d) and 8 in the main text it is evident that we need to consider the uncertainty relating to measured acceleration peaks and the reported 
$Mt$ numbers by propagating uncertainty through the acceleration reduction equation (
$1-a/a_q$) and the 
$Mt$ number formulation (3.1).
First, we propagate the uncertainty through 
$Mt$ formulation,
\begin{align} Mt &= \frac{t_2-t_1}{t_p} \nonumber\\ &= \frac{t_2-t_1}{\beta\sqrt{\dfrac{d_1}{2g}}}, \end{align}
where 
$d_1$ is the sphere diameter and 
$g$ is the gravitational constant. Applying the Taylor series method (TSM) approach to propagate uncertainty,
\begin{equation} u_{Mt}^2 = \left(\frac{\partial Mt}{\partial (t_2-t_1)}\right)^2u_{t_2-t_1}^2+ \left(\frac{\partial Mt}{\partial \beta}\right)^2u_{\beta}^2 + \left(\frac{\partial Mt}{\partial d_{1}}\right)^2u_{d_{1}}^2 .\end{equation} Now, the partial derivatives of 
$Mt$ with respect to the variables are,
\begin{equation} \frac{\partial Mt}{\partial (t_2-t_1)} = \frac{1}{\beta \sqrt{\dfrac{d_1}{2g}}},\quad \frac{\partial Mt}{\partial \beta} ={-}\frac{(t_2-t_1)}{\beta^2\sqrt{\dfrac{d_1}{2g}}},\quad \frac{\partial Mt}{\partial R} ={-}\left(\frac{2}{d}\right)^{1.5}\frac{(t_2-t_1)\sqrt{g}}{\beta}. \end{equation} The values 
$(t_2-t_1)$, sphere diameter 
$d_1$ and 
$\beta$ are experimentally measured values all of which has uncertainty associated with them. These uncertainties can be reported as,
\begin{equation} u_{(t_2-t_1)}=0.05\times 10^{{-}3},\quad u_{\beta}=0.0497,\quad u_{d_1} = 0.02\times 10^{{-}3}. \end{equation}Using these and substituting the partial derivatives of (A2) by (A3) we get,
\begin{align} u_{Mt}^2 &= \left(\frac{\sqrt{2g}}{2.03\sqrt{d_1}}0.5\times 10^{{-}3}\right)^2 + \left(\frac{(t_2-t_1)}{2.03^2\sqrt{\dfrac{d_1}{2g}}}0.0497\right)^2 \nonumber\\ &\quad +\left(\frac{(t_2-t_1)\sqrt{g}}{2.03\left(\dfrac{d_1}{2}\right)^{3/2}}0.02\times 10^{{-}3}\right)^2 . \end{align}
Thus (A5) gives the uncertainty for individual 
$Mt$ numbers explored in the study.
We can propagate uncertainty in the same manner for the acceleration reduction equation,
\begin{equation} R = 1-\frac{a}{a_q}, \end{equation}
where 
$R$ denotes the amount of reduction experienced for a two-sphere case, 
$a_q$ is the quiescent acceleration and 
$a$ is the two-sphere acceleration for any particular case. Propagating uncertainty through TSM,
\begin{equation} u_{R}^2 = \left(\frac{\partial R}{\partial a}\right)^2u_{a}^2 + \left(\frac{\partial R}{\partial a}_q\right)^2u_{a_q}^2 .\end{equation} Now, the individual uncertainties associated with the acceleration measurement come from the device uncertainty and the low sampling rate of the accelerometer (1000 Hz). The uncertainties relating to measured accelerations (
$u_{a_q}$ and 
$u_a$) are one sided, assuming the low sampling rate of the accelerometer fails to capture the exact highest peak of any acceleration event. Applying these uncertainties coupled with device uncertainty itself, we see that we would have two-sided uncertainty, where 
$u_{a_q}$ will drive the upper uncertainty for any data point since a bigger 
$a_q$ will increase the reduction value in (A6), and 
$u_a$ drives the lower uncertainty since a greater impact acceleration 
$a$ will reduce the amount of acceleration reduction. Taking these into account, we can rewrite (A7) as a two-sided uncertainty where,
\begin{equation} u_{up}^2 = \left(\frac{\partial R}{\partial a}\right)^2\times(u_{device})^2 + \left(\frac{\partial R}{\partial a_q}\right)^2\times(\sqrt{u_{device}^2+u_{a_q}^2})^2 \end{equation}
\begin{equation}u_{down}^2 = \left(\frac{\partial R}{\partial a}\right)^2\times(\sqrt{u_{device}^2+u_{a}^2})^2 + \left(\frac{\partial R}{\partial a_q}\right)^2\times (u_{device})^2. \end{equation}The partial derivative terms from (A8) and (A9) can be written as,
\begin{equation} \frac{\partial R}{\partial a} ={-}\frac{1}{a_q},\quad \frac{\partial R}{\partial a_q} = \frac{a}{a_q^2}. \end{equation} Device uncertainty 
$u_{device}$ can be calculated from the datasheets of the accelerometers used (MPU 9250A manufactured by Invensen Inc. and H3LIS331DL produced by STMicroelectronics). For calculating 
$u_a$ and 
$u_{a_q}$, we employed a different method. Based on the assumption that the low sampling rate of 1000 Hz is not enough to capture the highest point of impact acceleration, we assumed that the highest acceleration point for any two-sphere impact event would reside between the two highest points physically measured by the accelerometer. Then we fit two straight lines comprising of the two maximum points and a few other points on both sides of the true peak acceleration point. The lines would intersect in between the two measured highest acceleration point, and the difference between this intersection peak and measured peak will be our acceleration uncertainty 
$u_{a_q}$ and 
$u_a$. Figure 14(a) illustrates this idea of estimating the uncertainty relating to the true peak acceleration. Figure 14(b) shows the application of the uncertainty bands calculated this way to experimentally measured quiescent peak acceleration data. The measured value underestimates the peak values when compared to theoretical peak values (calculated using (3.11) with 
$C_d \sim 1$) with increasing 
$Fr$, but falls within the uncertainty bounds.
(a) Peak estimation using the line intersection method. The true peak should lie between the two highest measured peak values. Fitting two straight lines to these two data points results in an intersection point (denoted by the green circle) higher in magnitude than both the highest points and is in between the points also. The difference between the measured peak (blue circle in the plot) and the estimated peak (green circle) can be taken as the uncertainty value for this two-sphere impact event. (b) Measured quiescent acceleration data compared with theoretical estimations calculated using (3.11). At lower 
$Fr$, the measured acceleration follows (3.11). With increasing 
$Fr$, underestimation of peak acceleration may happen, but it happens mostly within the uncertainty bounds.
Using this approach, we estimated the peak acceleration values for each of the data points presented in figure 8. This also gives us an idea of how much uncertainty we are dealing with for the individual acceleration curves reported in figures 6(c) and 7(d). Applying the calculated uncertainty values for 
$u_{a_q}$ and 
$u_a$ and using (A8) and A9 we measured two-sided uncertainty for all the reported cases from figure 8. Figure 8 in the main text plots all the data points with their uncertainty bands presented.
A.2. Impact force dependence on added mass
The analysis here follows from Shiffman & Spencer (Reference Shiffman and Spencer1945) and May (Reference May1975), who essentially report the same analysis. For impact of a sphere on a quiescent free surface, their analysis starts with the notion that momentum is conserved during a water impact such that
\begin{equation} MU_o = \left(M+m\right)U, \end{equation}
where 
$M$ is the sphere mass, 
$m$ is the added mass of the fluid, 
$U_o$ the impact velocity and 
$U$ the instantaneous velocity of the sphere and added fluid mass after impact. As the force due to added mass dominates forces in the early stages of impact for the regime of our study, this force is given by
\begin{equation} F_z ={-}M \frac{\textrm{d} U}{\textrm{d} t}. \end{equation}
Solving (A11) for 
$U$ and taking the derivative with respect to time gives
\begin{equation} \frac{\textrm{d} U}{\textrm{d} t} ={-}U_o\frac{M}{ \left(m+M\right)^2}\frac{\textrm{d} m}{\textrm{d} t}. \end{equation}
The time derivative can be converted to a submergence-dependent derivative by first defining 
$b = s/R$, where 
$R$ is the sphere radius and 
$s$ the submergence below the initial free-surface plane (figure 15). Derivatives are related by
\begin{equation} \frac{\textrm{d}()}{\textrm{d} t} = \frac{\textrm{d}()}{\textrm{d} b}\frac{\textrm{d} b}{\textrm{d} s} \frac{\textrm{d} s}{\textrm{d} t} = \frac{U}{R}\frac{\textrm{d}()}{\textrm{d} b}. \end{equation}Combining (A11)–(A14) and simplifying gives
\begin{equation} F_z = \frac{U_o^2}{R}\frac{1}{ \left(1+\dfrac{m}{M}\right)^3}\frac{\textrm{d} m}{\textrm{d} b}. \end{equation}
Using the typical definition 
$F_z = \rho {\rm \pi}R^2 C_d/2$, where 
$C_d$ is the impact drag coefficient, gives
\begin{equation} C_d = \frac{2}{\rho{\rm \pi} R^3}\frac{1}{ \left(1+\dfrac{m}{M}\right)^3}\frac{\textrm{d} m}{\textrm{d} b}, \end{equation}
which Shiffman & Spencer (Reference Shiffman and Spencer1945) report for values of 
$m/M$. We can see from (A16) that the vertical force varies linearly with the rate of change of added mass with submergence.
Schematic diagram of a sphere penetrating the free surface.
(a) The non-dimensional pinch-off time 
$t_p u_1/d_1$ can be scaled with 
$\sqrt {2}Fr^{1/2}$, which agrees with previous such reports in the literature (Aristoff et al. Reference Aristoff, Truscott, Techet and Bush2010). (b) The relation between 
$Mt$ and non-dimensional spacing between the spheres show that, for the same spacing between different diameter spheres, 
$Mt$ will be different, leading to different reduction of impact forces. The on surface jets are cases which can be treated as extension to on jet cases.
The time the upwards moving bowl takes to evolve into a fully fledged Worthington jet from the moment of pinch-off is invariant with the impact parameters (a), resulting in 
$Mt = 1.43$ as the criterion for transition from bowl modes to on jet modes (b).
Acceleration plot for different two-sphere water entry cases. This plot shows the acceleration felt by the free-falling sphere, going back 100 ms before free-surface impact. In the on jet case shown by the green line, the sphere meets the upward moving Worthington jet significantly ahead in time, 
$\approx$80 ms before impact in the case presented here, the green circle in the plot denotes the approximate time when the sphere meets the jet.
A consecutive two-sphere impact event with 
$Mt=1.99$. A 23.75 mm sphere impacts the free surface at a velocity of 
$\sim$3.6 m s
$^{-1}$ creating a cavity. The cavity domes over at the surface (29 ms) that suppresses the axial jet originating from pinch-off (70–119 ms). This results in a smaller jet at the free surface through which the upper sphere falls (129–146 ms). Similar cases are termed as on surface jet cases.
For the on bowl and inside bowl cases, the cavity has pinched off and the fluid has significant upward momentum as the bowl retracts. We can rationalize this effect to first order by again considering conservation of momentum before and after impact
\begin{equation} MU_o - mU_m= \left(M+m\right)U, \end{equation}
where 
$U_m$ is the average velocity of the added fluid mass just before impact (assumed to be in the direction opposite of 
$U_o$). We follow the same derivation for impact on a quiescent free surface as above. Solving (A15) for 
$U$ and differentiating with respect to time gives
\begin{equation} \frac{\textrm{d} U}{\textrm{d} t} ={-}(U_o+U_m)\frac{M}{\left(m+M\right)^2}\frac{\textrm{d} m}{\textrm{d} t}. \end{equation}Applying (A14) yields
\begin{equation} \frac{\textrm{d} U}{\textrm{d} t} ={-}\frac{M(MU_o - mU_m)}{ \left(m+M\right)^3} \frac{(U_o+U_m)}{R}\frac{\textrm{d} m}{\textrm{d} b}. \end{equation}Finally, inserting (A19) into (A12) and simplifying gives
\begin{equation} F_z = \frac{\left(U_o-\dfrac{m}{M}U_m\right)}{\left(1+\dfrac{m}{M}\right)^3} \frac{\left(U_o+U_m\right)}{R}\frac{\textrm{d} m}{\textrm{d} b}. \end{equation}
Now the vertical force depends not only on 
$\textrm {d} m/\textrm {d} b$, but also has a first-order dependence on 
$U_m$.
A.3. Transition from on bowl to on jet mode
Figure 17(a,b) shows that the non-dimensional time of the Worthington bowl evolving into a fully fledged Worthington jet does not vary with increasing 
$Fr$. Employing the mean 
${t_j}/{t_p}$ value in (3.10) gives us the transition 
$Mt \sim 1.43$, which acts as an excellent separation criterion between the bowl mode and the on jet mode (figure 8).
A.4. Surface jet cases and expansion of reduction regime
For a few of the higher 
$Fr$ cases, the splash crown domeover at the surface suppresses the Worthington jet. This suppression leads to the creation of a not so high jet-like water column at the free surface. The upper sphere falls through this jet into the water pool. One such example is given in figure 19, where an 
$Mt=1.99$ impact event is presented.
Open access
