Yu Song Sun^{1} , Sui Hua Zhou^{2} , Xiao Bing Zhang^{3} , Ya Li Xiang^{4}
^{1, 2, 3, 4}Naval University of Engineering, Wuhan, China
^{1}Corresponding author
Vibroengineering PROCEDIA, Vol. 17, 2018, p. 137141.
https://doi.org/10.21595/vp.2018.19793
Received 28 February 2018; accepted 7 March 2018; published 20 April 2018
Copyright © 2018 JVE International Ltd. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Aiming at the characteristics of high speed water entry of heavy projectiles, this article studied the highspeed water impact load of heavy projectiles by the multimaterial ALE method, good agreements were obtained compared with the empirical formula, then analyzed the two parameters of impact velocity and the projectile’s mass which will affect the impact load when changed. Numerical simulation results are as follows: under vertical water entry condition, the acceleration of the projectile is roughly proportional to the square of the impact velocity and inversely proportional to the projectile’s mass during water impact.
Keywords: hemisphericalnose projectile, waterimpact, multimaterial ALE method.
The problem of water entry has been widely concerned for a long time. Research on the water entry phenomena has been motivated by an interest in the load during water impact [13], the shape of cavity [4] and the generation and propagation of the ballistic waves into the fluid [5, 6]. In recent years, with the development of the finite element method (FEM), more and more researchers have begun to use finite element method to solve the problem of water impact aiming at a specific condition. Wang Jian investigated the phenomenon of the rocket sled waterbrake high speed water entry by using the method of fluidsolid coupling. The modalities of water were calculated such as supercavitations, uplift, spatter and pressure, and the velocity, deceleration and resistance coefficient varied with time of waterbrake was obtained [7]. Ma Qingpeng studied the multiphase flow and the movement of a cylinder body with a cone head by using the method of volume of fluid (VOF) and dynamic mesh. The velocity attenuation and induced cavity shape were obtained with an initial velocity of 500 m/s. Good agreements were obtained compared with the theoretical results [8].
Investigations on the highspeed water entry of small projectiles (with a diameter of 1 cm5 cm) are relatively mature [810]. Due to the limitation of experimental conditions, only a few researches have investigated the waterimpact of super heavy projectiles (with a diameter of 40 cm60 cm) at a relatively low speed (40 m/s100 m/s) [11, 12]. While the investigation of high speed (200 m/s300 m/s) water entry of heavy projectiles (with a diameter of 10 cm40 cm) is relatively rare. The method of multimaterial ALE is used in this paper to investigate the impact load of high speed water entry of hemisphericalnosed heavy projectiles.
High speed water entry involving the flow of three phases contains fluid, water and solid. Multimaterial ALE (Arbitrary Lagrangian Eulerian) method is good at investigated it numerically.
ALE method includes mass, momentum and energy conservation equations are as follows [7]:
where $\rho $ is the fluid density, ${v}_{i}$ is the fluid velocity, $\upsilon $ is the hydrodynamic viscosity coefficient. ${\mathrm{\omega}}_{i}$is the convection velocity, ${x}_{i}$ is the Euler coordinate value, ${\sigma}_{ij}$ is the stress strain, ${b}_{i}$ is the body force, $E$ is the total specific energy.
The eightnode hexahedron element is used in the model. The fluid grids are defined as ALE with the material of NULL. The water’s EOS (equation of state) is Gruneisen and the air’s EOS is Linearpolynomial. The solid grids are defined as Lagrangian and the material of the projectile is common structural steel. The projectile was in the air at the initial time. Cause of the symmetry of the model, the waterentry model can be simplified as a quarter models in order to reduce the work of computation. The numerical simulation model is as illustrated in Fig. 1.
The profiles are defined as symmetry planes and the rest of surfaces are defined as nonereflect boundaries in order to simulate the unlimitedsurroundings. Using level set method to tracking the interface of water and air and using operator splitting method to solve the fluid solid coupling problem.
Set the hemisphericalnosed heavy projectile’s diameter is 40 cm, height is 120 cm and weight is 171 kg, with an initial velocity of 200 m/s from a height of 0.1 m above the liquid surface, then the cavity formed during the water entry simulated numerically are shown as Fig. 2 and the curve of the projectile’s acceleration, which represents the impact load, with time is shown as Fig. 3.
Fig. 1. Model of projectile’s waterentry vertically
Fig. 2. Cavity formed during the water entry
Fig. 3. Projectile’s acceleration curve during water entry
During the water entry, the surface water surrounds the cavity is uplifted. A spindly huge cavity formed behind the projectile’s nose.
The projectile’s acceleration is relatively high during the water impact, and then the projectile’s acceleration maintains a relatively stable value after the impact cause that the penetration velocity maintains a relatively stable value during the short time.
During the water impact, the projectile’s velocity can be considered as a constant, May A. considers the cavity formed during early water entry can be approximated as a parabola except the area near the projectile’s nose. Theoretical predictions and experimental measurements predicted that the shape of the cavity can be expressed as below [13]:
where ${d}_{b}$ is the projectile’s diameter, ${C}_{d}$ is the drag coefficient of a cavitation vehicle, for spheres, the correlation of the drag coefficient shows that there are ranges classified by velocity, roughly “subsonic” up to $M=$ 0.5, ${C}_{d}=$ 0.384 [4].
In order to keep the projectile’s velocity stable during water impact, here set the projectile’s weight as 740 kg. Comparisons of the cavity formed during water entry by the heavy projectile are shown as Fig. 4, which showed that numerical results have a good agreement with the analytical results, then the validity of the simulation model is verified.
Fig. 4. Comparisons of the cavities formed during the water entry
From the investigation above, it can be concluded that the maximum value of projectile’s load appears during the early stage (within 0.2 ms after touching the water surface) of water entry. This paper also numerically calculated the acceleration curve during water entry of the projectile with impact velocities of 240 m/s, 280 m/s and 320 m/s, and the comparison of the acceleration curve are shown as Fig. 5, the comparison of the peak acceleration under different impact velocity is as Table 1.
Table 1. Different impact velocities and corresponding peak acceleration
${v}_{0}$ (m/s)

${a}_{max}$ (10^{4}^{}m/s^{2})

${a}_{max}/{v}_{0}^{2}$ (m^{1})

200

2.07

0.52

240

3.04

0.53

280

4.11

0.52

320

5.13

0.50

It can be concluded that the maximum value of the projectile’s acceleration is roughly proportional to the square of the impact velocity.
In order to investigate the difference of different mass projectiles’ impact load, this paper also numerically calculated the acceleration during water entry of the projectiles with different mass of 142 kg, 171 kg, 199 kg and 228 kg, and the comparison of the acceleration curve are shown as Fig. 6, the comparison of the peak acceleration is shown as Table 2.
A conclusion can be derived from Fig. 6 and Table 2 that projectiles’ acceleration is roughly inversely proportional to the projectile’s mass during the water impact.
Fig. 5. Comparison of projectile’s acceleration curves
Fig. 6. Comparisons of different mass projectiles’ acceleration curves
Table 2. The mass of each projectile and its corresponding peak acceleration
$m$ (kg)

${a}_{max}$ (10^{4}^{}m/s^{2})

$m{a}_{max}$ (10^{6}^{}N)

142

2.27

3.2

171

2.07

3.5

199

1.92

3.8

228

1.55

3.5

This article studied the highspeed water entry phenomena of heavy projectiles by the multimaterial ALE method, good agreements were obtained compared with the empirical formula. By changing the parameters of simulation conditions, two conclusions have been obtained.
1) The maximum value of the projectile’s acceleration is proportional to the square of the impact velocity, and the duration time of maximum acceleration will decrease with the increasing of the impact velocities.
2) The maximum value of the projectile’s acceleration is roughly inversely proportional to the projectile’s mass during the water impact.