Experimental Study of a Sphere Bouncing on the Water

Xujian Lyu · Honglu Yun, · Zhaoyu Wei

Abstract In this paper, the flow physics and impact dynamics of a sphere bouncing on a water surface are studied experimentally.During the experiments, high-speed camera photography techniques are used to capture the cavity and free surface evolution when the sphere impacts and skips on the water surface. The influences of the impact velocity (v1) and impact angle (θ1) of the sphere on the bouncing flow physics are also investigated, including the cavitation evolution, motion characteristics, and bounding law. Regulations for the relationship between v1 and θ1 to judge whether the sphere can bounce on the water surface are presented and analyzed by summarizing a large amount of experimental data. In addition, the effect of θ1 on the energy loss of the sphere is also analyzed and discussed. The experiment results show that there is a fitted curve of v1 =17.5θ1 −45.5 determining the relationship between the critical initial velocity and angle whether the sphere bounces on the water surface.

Keywords Water entry · Bouncing sphere · Critical impact angle · High-speed photograph · Splash

1 Introduction

The impact dynamics of structures on the water have been studied for several decades because of their wide engineering applications (Glasheen and Mcmahon 1996; Miloh and Shukron 1991; Shlien 1994; Von Karman 1929), such as the landing of a seaplane and a speedboat gliding on water. Bouncing and skipping a stone and the swift running of insects or basilisk lizards are other interesting natural phenomena on a water surface that attract the attention of scientists (Clanet et al. 2004, Glasheen and Mcmahon 1996). The bouncing mechanism of an elastic sphere on a water surface was studied by Belden et al. (2016) through experiments and numerical simulations. They found that the bouncing dynamics is decided by the ratio of the material shear modulus to hydrodynamic pressure and the wave propagation speed. Hurd et al. (2019) revealed a new mode of bouncing on the water for relatively soft spheres with low impact angles by experiments, which shows good agreement between measured acceleration, number of skipping events,and distanced traveled. Liu and Smith (2014) proposed a new prediction formula for repeated impacts and rebounds,especially over short time scales. Korobkin et al. (2011)summarized and analyzed the mathematical problems of modern hydroelastic mechanics and analyzed some models and techniques. Hewitt et al. (2011) experimentally studied the planning and bouncing of rectangular paddles on a water surface and established a model using shallow water theory.By studying the rebound effect of objects in the shallow layer of a liquid, Hicks and Smith (2011) proposed a model that couples the motion of an object to fluid dynamics.

In addition, projectile and stone skipping is also quite an interesting phenomenon that attracts the attention of researchers. Nishida et al. (2010) studied the dynamic properties of a projectile impinging on a granular medium and derived a formula for determining the critical angle by the density ratio and diameter ratio through experimental techniques. Rosellini et al. (2005) studied the hydrodynamic characteristics of a rotating disk impacting water with multiple bounces. They found that the source of dissipation lies in the dependence of this reaction force on the angle between the water surface and the trajectory of the stone. In addition, they proposed a simple model to measure bounce dissipation. Bocquet (2003) analyzed the physics of stone skipping and found that the maximum number of bounces can be estimated by the deceleration and angular stability of the stone. Faltinsen (2000) reviewed hydroelastic slamming experiments and theoretical studies and noted that hydrodynamics research on slamming problems must be from a structural point of view. Johnson (1998) summarized the theory of ricochet from a liquid surface of a non-rotating sphere and improved the non-rotating projectile portion.Truscott et al. (2014) summarized experimental, theoretical,and numerical studies of an object impacting on the water and introduced several problems to be solved, including the high-fidelity measurements of acceleration, surface stress,and cavity pressure.

The above literature shows that most of the studies addressed aspects of the law of bouncing plate-shaped objects on the water, for instance, the number of bounces.The phenomenon of cavity evolution and kinetic energy dissipation of a symmetric and regularly shaped object(i.e., a sphere) is further from being clearly explained. This forward-looking paper reveals the cavity dynamics and the mechanism of a sphere bouncing on a water surface.Experimental evaluation of the effects of the initial impact angle and the energy dissipation and estimation of the law of the critical impact angle are provided for the bouncing of a sphere on a free water surface.

2 Experimental Methods of High‑Speed Photograph

The experiment was conducted in a water tank located in Nanjing University of Science and Technology, China. During the experiments, the temperature was 23 to 27 °C, and the closed laboratory eliminated the interference of natural light and airflow. The size of the tank was 1500 (L)×400 (W)×600 mm (H), and the water depth was set to 350 mm, as shown in Figure 1. The water tank was made of 10-mm thick optical glass to reduce the optical path loss. The steel sphere used in the experiment has a massm≈1.397 g, a diameterd≈7 mm, and densityρs≈7778.6 kg/m3. The measurement error of these parameters was less than 0.5%. The sphere was launched by an electromagnetic launcher similar to the equipment described in Yun et al.(2020), and the arrangement is shown in Figure 1(a).

Figure 1 Sketch of the experimental setup

A high-speed digital video camera, Phantom VEO410L,with a 105-mm fixed-focus Nikon lens and a resolution of 1280×400 pixels, was used to capture the sphere’s movement and the associated cavity evolution. The horizontal distance between the camera and the electromagnetic launcher was 3.4 m, which allows the microlens to capture the whole flow physics during the impact process, and the camera center was approximately 0.35 m above the ground, which is nearly the same height as the free water surface. For the quality of the flow visualization experiment, an adequate supplementary light source is fundamental. In the experiment, four highfrequency lamps of the type JINBEI EF-200 LED and a rectangular LED light box were used to provide sufficient light intensity. Four high-frequency lamps were used as the foreground lighting arranged in front of the water tank, and the light box was used as the backlighting to provide a uniform white background arranged at the back of the water tank, as shown in Figure 1. On the basis of these strong lightings, the frame rate was set as high as 10 000 fps to capture significantly crisp images to accurately analyze the trajectory of spheres,and the exposure time was decreased to 10 μs.

Seven typical cases with similar initial impact velocities and different impact angles were selected as analysis objects,and the parameters are shown in Table 1. The main parameters include the initial impact velocityv1, initial impact angleθ1, bouncing velocityv2, and bouncing angleθ2; definitions of the variables (i.e.,v1,θ1,v2, andθ2) are shown in Figure 1(a).

3 Results and Discussion on the Sphere Bouncing and Cavity Characteristics

Cavity dynamic characteristics and the kinematics law of a steel sphere bouncing on the water will be discussed in this section.

Table 1 Main parameters of typical cases

3.1 Cavity Formation of a Sphere Bouncing on the Water

To analyze the hydrodynamic characteristics of the sphere,t0is used to indicate the time when the sphere just touches the water surface. Figure 2(a) shows the selected image sequence for case 1 with an initial impact velocity ofv1=58.58 m/s and an initial impact angle ofθ1=1.04°. Because of the very small initial impact angle and the larger initial impact velocity, the center of the sphere was always above the static-free surface throughout the bouncing process, and the sphere was in contact with water for a short period of only 1.5 ms.The cavity is seen to be the shallowest and the splash is weak. Furthermore, the bouncing velocity ofv2=56.81 m/s is quite near the initial impact velocity ofv1=58.58 m/s.The bouncing angle ofθ2=0.94° is near the impact angle,indicating that a slight amount of energy is lost during the water surface impact.

Figure 2 Cavity evolution of a sphere bouncing on the water for respective cases 1 to 4

The typical images of case 2 are shown in Figure 2(b),and the initial impact angle of the sphere was 2.60°. It can be shown that the splash above the free surface is much stronger, and the cavity is relatively deep and large compared with that of case 1. Moreover, the duration of the sphere touching the water is extended slightly to Δt=2.5 ms, compared with Δt=1.5 ms for case 1. Figure 2(c) and (d) show the selected image sequences for case 3 and 4, respectively.Note that the initial impact velocities for the four cases are similar and the initial impact angles are gradually increased.It is seen that the larger the initial impact angle is, the deeper and larger the profile of the cavity formed. Furthermore, it can be seen in case 4 that a smooth and complete water curtain is formed, which almost encloses the entire sphere as it leaves the water surface (Figure 2(d)). Another significant phenomenon seen in cases 1 to 4 is that the bouncing velocity decreases as the initial impact angle increases, although the initial impact velocity is the same in each case.

It is obvious that the flow physics of case 5 is different from those of the other cases, as shown in Figure 3. For case 5,v1is 57.88 m/s andθ1is 5.89°. Notably, the duration that the sphere moved under the free surface is much longer for case 5, reaching 46 ms. Furthermore, the sphere is displaced 675 nm in the water, nearly 100-fold its diameter. A complete cavity, including the smooth water film above the static water surface, encloses the sphere as it moves in the water. The cavity wall behind the sphere finally becomes rough because of the cavity collapse and the effect of capillary waves on the free surface, as described by Grumstrup et al. (2007).

The cavity finally collapses to many small air bubbles,which is possibly associated with the hydrostatic pressure under the free surface. The cavity also causes violent-free surface elevation, where an obvious boundary between the cavity and the surface water film can be noticed. The water film later shrinks into a long water column, possibly associated with the surface tension. The bouncing velocity is justv2=6.22 m/s, much smaller than the initial impact velocity ofv1=57.88 m/s. More than 98.8% of the sphere’s kinetic energy is lost during the impact process, which unboundedly transfers to the energy of the water film and cavity.

Figure 3 Images of cavity evolution for case 5 with v1=57.88 m/s and θ1=5.89°

To clearly analyze the phenomenon and mechanism of a sphere bouncing on a water surface, two experimental cases with different initial impact angles were selected.Figure 4 shows the image sequence of cavity evolution for the sphere in case 6. The initial impact velocity for case 6(v1=57.68 m/s) is almost the same as that of case 5, while the initial impact angle (θ1=6.01°) is slightly larger. The early-stage flow physics for case 5 and case 6 is quite similar,including the cavity and evolution of water film. Nevertheless, obvious differences are noticed in the flow physics and movement of the sphere for the two cases. In case 6, the sphere eventually enters the water after sliding a long distance near the free surface, as shown in Figure 4. The main reason for the difference is that the lift force of the water on the sphere is less than the gravitational and inertial forces associated with the larger initial impact angle for case 6.Eventually, the kinetic energy of the sphere gradually dissipates under the free surface and the sphere continues to move underwater.

Figure 5 shows a selected image sequence of the sphere in case 7 hitting the water with a larger initial impact angle ofθ1=10.23°. The initial impact velocity wasv1=56.78 m/s,close to those of other cases. It can be seen that the sphere enters directly into the water. A perfect and smooth cavity is formed, which expands and is elongated as the sphere descends in the water. The contraction of the tail of the cavity is obvious under the free surface in Figure 5. Moreover,the water curtain above the free surface has a clear tendency to expand outward, accompanied by the forward movement of the sphere. At the later stage, the cavity wall becomes rough when it collapses because of the effects of hydrostatic pressure and surface tension.

3.2 Mechanism of a Sphere Bouncing on the Water

Figure 4 Cavity evolution and movement of a sphere with v1=57.68 m/s and θ1=6.01°for case 6

Figure 5 Evolution of the cavity caused by the sphere with v1=56.78 m/s and θ1=10.23°for case 7

Figure 7 sketches the evolution of the cavity profile and the brief force analysis for case 4, including the whole process of the sphere impacting, gliding, and bouncing on the water surface. The blue dotted lines indicate the profile of the splash, and the red lines are the contact area between the sphere and water. Because of the interaction between the sphere and water in the contact area, the sphere is subjected to a lift force, as shown in Figure 7. During the process of the impact, the external forces acting on the sphere are the weightmg, the buoyancy forceFb, the steady-state drag forceFd, the capillary forceFcdue to surface tension, and the hydrodynamic forceFh, also called the resistive force. On the basis of Newton’s second law, the equation of motion of a sphere bouncing on the water can be written as follows:

Figure 6 Typical images for cases 1–4 with different initial impact angles

Figure 7 Sketches for cavity configurations and force characteristics analysis of a sphere

In the process of the sphere bouncing on the water, the hydrodynamic forceFhis the dominant force, and the direction and magnitude ofFhchange with the direction of movement. The entire force on the sphere can be decomposed into two components: horizontal drag and vertical lift. The drag reduces the velocity, and the lift eventually changes the direction of the sphere. Because of the external forces, the sphere bounces on the free surface.

3.3 Force Characteristics Analysis of a Sphere Bouncing on the Water

In this section, the kinematics and dynamic characteristics of a sphere bouncing on a water surface are analyzed quantitatively. Figure 8 presents the displacement, velocity, and acceleration curves associated with cases 5 and 6 in the horizontal and vertical directions. The analytical method is similar to that described by Epps et al. (2010)and Wei and Hu (2014, 2015). The line (x, 0) indicates the static-free surface, and the point (0, 0) is the impact point of the sphere on the water. The initialθ1andv1for the two cases are quite similar, while the impact phenomena are significantly different.

Displacement curves in cases 5 and 6 show that the sphere trajectories are the same at the initial stage of impact (referring to Figure 8(a)). As the sphere moves forward, its trajectory gradually is deflected upward and eventually sends it rushing out of the free surface in case 5. The obvious difference is that the sphere’s trajectory in case 6 shifts downward and is gradually moved away from the static-free surface. The velocity variations of spheres along thex- andy-directions are shown in Figure 8(b) and(c). It is found that the tendencies for the two cases are similar along thex-direction (Figure 8(b)), while the tendencies along they-direction gradually differ over time,as shown in Figure 8(c). It can be seen that the vertical velocity of the sphere in case 5 changes from negative to positive att=11 ms, which indicates that the sphere reaches the deepest point in the cavity and begins to move upward at this time. In contrast, the vertical velocity of the sphere in case 6 almost levels off after a rapid decrease in the early stage of impact.

Figure 8 Variation of parameters for cases 5 and 6

The acceleration curves in thex- andy-directions in Figure 8(d) almost overlap in the early stage. Nevertheless,slight differences are noticed along they-direction after impact, as shown in detail in Figure 8(d). The acceleration in case 5 is positive, which causes the sphere to ascend gradually until rushing out of the water. In contrast, the acceleration in case 6 is negative, which leads the sphere to move downward throughout the process of impacting the water.

3.4 Statistics of a Successful Bounce

To reveal the bouncing law of the sphere on a free water surface, numerous experiments have been conducted. A schematic diagram of the statistical distribution ofθ1andv1during the experiments is shown in Figure 9. The black hollow squares stand for the cases in which the sphere bounces on the water surface, and the red solid triangles are the cases without bouncing. A conspicuous boundary is found to differentiate the phenomena of bouncing and non-bouncing. The boundary line is shown with a blue dotted line in Figure 9 and can be fitted as follows:

Figure 9 Distribution of the initial impact angle (θ1) and initial impact velocity (v1)

For the process of a sphere bouncing on the water, another significant characterization parameter is the bouncing angleθ2. The distribution ofθ2andθ1is shown in Figure 10. Experimental data points of the bouncing cases can be distinguished through the durationt, which is the period from the sphere impacting on the water to leaving the static-free surface. These black solid squares replace the cases in which the durationtis smaller than 10 ms, and the red hollow triangles replace the cases withtlarger than 10 ms. Note that the blue dotted line corresponds toθ1=θ2. It is clearly seen that the statistical data points are mainly distributed near the blue dotted line, and the distribution is more concentrated for the cases oft≤10 ms. For cases with the spheres gliding on the water for a longer time(t>10 ms), the bouncing angleθ2tend to be smaller than the initial impact angleθ1.

As mentioned above, the energy is transmitted from the sphere to water during the bouncing process. Figure 11 shows the statistical data points for the variation in kinetic energy loss percentageηwith bouncing angle. The loss of kinetic energy is more obvious at the moment for the cases with larger initial impact angles. Furthermore, the percentage of kinetic energy loss is below 50% when the impact angle is smaller than 3°.The energy attenuation rate corresponding to different angles is between two curves for all the cases, as shown in Figure 11,and the two curves can be fitted using the following equations:

Figure 10 Distribution of initial impact angle (θ1) and bouncing angle distributions (θ2)

Figure 11 Variation of the percentage of kinetic energy loss (η) with bouncing angle (θ1)

4 Conclusions

In this work, the cavity configurations and kinematic characteristics of a sphere bouncing on water are experimentally investigated. Through continuous improvement of experimental methods, high-resolution images of the cavity evolution of a sphere bouncing on the water are obtained, and they show the evolution process of cavitation more intuitively and facilitate the analysis of complex flow phenomena during the process. The motion and force characteristics of the sphere in typical cases are extracted and analyzed by an advanced image processing technique. More than 200 cases of a sphere impacting on a water surface have been considered, and the effect of initial impact velocity and initial impact angle on the phenomena of sphere bouncing has been investigated. The results indicate that a critical curve ofv1=17.5θ1−45.5 shows the relationship between the initial velocity and angle and determines whether the sphere bounces on the water surface. In addition, for the impacting cases with the same initial velocity, the energy attenuation of the sphere is closely associated with the initial impact angle during bouncing, and more than half of the kinetic energy is lost for an impact angle ofθ1> 3°.

FundingSupported by the Fundamental Research Funds for the Central Universities (30918012201) and the Fund of the State Key Laboratory (6142604190302).

Open AccessThis article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source,provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.