Next Article in Journal
From Time-Series to Hybrid Models: Advancements in Short-Term Load Forecasting Embracing Smart Grid Paradigm
Previous Article in Journal
Mitigating the Impact of Asymmetric Deformation on Advanced Metrology for Photolithography
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Study on Water Entry of a 3D Torpedo Based on the Improved Smoothed Particle Hydrodynamics Method

1
China Ship Scientific Research Center, Wuxi 210084, China
2
Key Laboratory of Unsteady Aerodynamics and Flow Control, Ministry of Industry and Information Technology, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China
3
School of Physical and Mathematical Sciences, Nanjing Tech University, Nanjing 211816, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2024, 14(11), 4441; https://doi.org/10.3390/app14114441
Submission received: 28 March 2024 / Revised: 15 May 2024 / Accepted: 21 May 2024 / Published: 23 May 2024
(This article belongs to the Section Marine Science and Engineering)

Abstract

:
The water entry of a torpedo is a complex nonlinear problem, involving transient impact, free surface deformation, droplet splashing, and fluid–structure coupling, which poses severe challenges to traditional mesh methods. The meshless smoothed particle hydrodynamics (SPH) method shows unique advantages in capturing the complex features of the water entry of the torpedo at different entry angles. However, it still suffers from some inherent shortcomings, such as low surface discretization accuracy, poor discretization flexibility, and low calculation efficiency. In this study, an improved adaptive SPH algorithm is proposed to investigate the water entry of the torpedo accurately and efficiently. This method integrates meshless point generation and adaptive techniques simultaneously. The numerical results demonstrate that when the torpedo vertically enters the water at different velocities, the induced impact loads acting on the head of the torpedo fluctuate significantly with two peak values in the initial stage and thereafter stabilize in a later stage. The impact load acting on the torpedo, the entry depth of the torpedo, the splash height of the droplets, and the size of the cavity generated around the torpedo increase with the increment in the entry velocity. When the torpedo enters the water at different entry angles under the same initial entry velocity, both the vertical and the horizontal movements of the torpedo are observed, which results in more complex variations in parameters along the x- and y-axes. The findings and the corresponding numerical method in this study can provide a certain basis for the future designs of the entry trajectory and the structural bearing capacity of torpedoes.

1. Introduction

The water entry of a torpedo is a typical example of water entries of structures, relating to large deformation of the free liquid surface. The impact load and the corresponding moving status induced by the water entry of the torpedo may result in strength damage to the torpedo structure or internal device failure. Thus, they should be carefully considered in the design and analysis of a torpedo. Three types of research methods are widely employed in investigating water entry processes of torpedoes, including the experimental method [1], theoretical method [2], and numerical method [3]. With the rapid development of computer technology, the numerical method has become a mainstream research method owing to its essential advantages of low cost and fast calculating speed compared to the other two methods. The traditional mesh-based numerical methods, such as the volume of fluid (VOF) [4,5], the finite volume method (FVM) [6], the boundary element method (BEM) [7], and the constrained interpolation profile (CIP) [8], are widely applied in simulating phenomena of fluid–structure interactions. However, meshes generated by these methods are prone to distortions, leading to inaccurate calculation results when applied to problems involving large deformations of free interfaces, such as the water entry of a torpedo. In recent years, the meshless numerical method has been extensively developed owing to its unique natural advantages in overcoming the problems above [9]. Rather than using meshes as the computational basis, it discretizes the computational domain into particles and then solves various integral equations or partial differential equations (PDEs) on the particles.
The smoothed particle hydrodynamics (SPH) method is a typical meshless numerical method that uses a collection of particles with specific physical properties, such as density, viscosity, mass, velocity, and temperature, to represent fluids [10]. The movement and the changes in the properties of particles are controlled via mass, momentum, and energy conservation equations. The numerical solution to these equations is obtained by employing the kernel approximation and the weighted integration of particle interpolation [11]. The SPH method has an inherent advantage in solving strongly nonlinear problems of the water entry of the torpedo, induced by complex free surface splashing and large deformation. Lots of work has been conducted to improve the numerical accuracy and stability of the SPH method in capturing the moving boundary, the large deformation of the fluid, and the free surface breakage in water entry processes [12,13,14,15,16,17,18,19,20,21,22,23,24,25]. However, there remain certain shortcomings in the SPH method. First, the SPH method performs calculation processes based on a series of discrete particles. Most studies have, first, commonly applied meshes over the computational domain and then extracted particles from the meshes to discretize the computational domain [26]. This existing particle generation method has low surface discretization accuracy and poor discretization flexibility. Second, the SPH method widely adopts a uniform particle distribution over the entire computational domain. To accurately capture the dramatic variations in the flow field around a moving three-dimensional (3D) torpedo, the total number of particles allocated in the computational domain is often very large. It consumes large amounts of computer memory and processing power and thus makes the SPH method less efficient.
This work thoroughly investigates the complex water entry of a 3D torpedo by applying an improved adaptive SPH algorithm. First, a meshless and efficient point distribution method is established for the 3D computational domain of the water entry of the 3D torpedo. Second, based on the adaptive SPH technique reported by Chiron et al. [27], a dynamic refinement and coarsening scheme of particles within the target region is proposed. Third, the improved adaptive SPH algorithm is validated by conducting a simulation of the water entry of a 3D horizontally placed cylinder, Finally, a parametric study on the effects of different entry velocities and angles of the torpedo on the impact loads and the flow fields around the torpedo is conducted using the validated adaptive SPH algorithm.

2. Numerical Method and Model Development

2.1. Mathematical Model Based on the SPH Method

The SPH method is a Lagrangian particle method, wherein the computational domain is composed of a set of discrete particles. Each particle is assigned physical properties such as mass, density, pressure, and velocity. Based on these particles, the time derivatives of various physical quantities can be obtained using the SPH kernel approximation algorithm, which allows the time integration of the physical process of the flow. The kernel approximation of a field function f ( r ) at a spatial position r, f ( r ) , can be expressed as
f ( r ) = Ω f ( r ) W ( r r , h ) d V r
where W ( r r , h ) is the smoothing kernel function [28], h is the smoothing length [29], and V is the volume. Ω represents the support domain of the smoothing kernel function, which is a function of h. A good tradeoff between the computational accuracy and the kernel approximation cost can be realized using the Wendland kernel function [30], which is expressed as
W ( r r , h ) = α D 1 q 2 4 ( 2 q + 1 )
where D is the dimensions of the computational domain and q is the dimensionless variable, expressed as q = (rr*)/h. In addition, αD is a normalization factor which is given as 7/(4πh2) for a two-dimensional (2D) domain and 21/(16πh3) for a 3D domain. As the computational domain is composed of a set of discrete particles with particular physical properties, the continuous format of the kernel approximation, as expressed in Equation (1), should be transformed into a discrete one:
f ( r i ) = j = 1 N f ( r j ) W ( r i r j , h ) V j
where N is the total number of particles located in the support domain of particle i. Vj is the virtual volume of the neighboring particle j, which is the ratio of the mass mj to the density ρj. Consequently, the kernel approximation of the derivative of the field function in the discrete format, f ( r ) , can be derived as
f ( r i ) = j = 1 N f ( r j ) i W ( r i r j , h ) m j ρ j
Equation (4) indicates that the derivative of the field function can be transferred to the derivative of the smooth kernel function, thereby reducing the smoothness requirements for the field function and producing stable results in numerical approximations of PDEs [31]. Thus, the governing equations in the form of weakly compressible SPH are derived as
d ρ i d t = ρ i j = 1 N ( u i u j ) i W ( r i r j , h ) m j ρ j d u i d t = 1 ρ i j = 1 N ( P i + P j ) i W ( r i r j , h ) m j ρ j + g d r i d t = u i P i = ρ 0 c 0 2 γ ρ i ρ 0 γ 1
where u is the velocity. The pressure of particle i, Pi, is expressed in the form of the Tait equation, where ρ0 is the density of the fluid when P = 0 Pa, c0 is the speed of sound [32,33], and γ is a constant [34,35] that usually takes the value of 7 in hydrodynamic problems. The continuity equation employs an asymmetric SPH approximation scheme to minimize errors arising from discontinuities [36], and the momentum equation adopts symmetric SPH approximation schemes to ensure better conservation [37].

2.2. Six-Degree-of-Freedom (6DOF) Motion Algorithm

Assuming that the structure of the torpedo is a rigid body, the structure is discretized into a series of particles which participate in the computation of the flow field [38]. The force per unit mass of each rigid particle is given as
f i = j F P N f i j
where fi denotes the force acting on the rigid body particle i per unit mass, FP represents the fluid particle, and N is the total number of particles within the support domain. fij denotes the force per unit mass exerted by the fluid particle j on the particle i.
For a moving rigid body, the motion equations of the six-degree-of-freedom (6DOF) rigid body can be written as follows:
M d V d t = i B P N m i f i I d ω 0 d t = i B P N m i ( r i R 0 ) × f i
where M denotes the mass of the rigid body, V is the velocity of the center of mass, and BP represents the rigid body particle. ri denotes the position vector of particle i. I denotes the moment of inertia of the rigid body, ω0 is the angular velocity of the rigid body, and R0 is the position vector of the center of mass of the rigid body. The velocity of each mass point on the rigid body can be obtained by
V i = V + ω 0 × ( r i R 0 )

2.3. Spatial Discretization and Adaptive Algorithm

2.3.1. Spatial Discretization

In contrast to the point generation for a curved line boundary in a 2D space, that in a 3D space is considerably more complex. In addition, owing to the meshless property of the algorithm, the generated points over the curved surface lack connecting relationships and thus lack the structure required for the traditional ray method [39]. To address these challenges, as shown in Figure 1, this work firstly applies a mapping method [40] to discretize the curved surfaces into points. Consequently, an improved method of the vector dot product is employed to evaluate the attributes of background Cartesian points. Finally, attribute-based background grids are used to pre-mark the computational domain and thus accelerate the spatial discretization.
The process of generating the discrete points using the mapping method is shown in Figure 2. It comprises the following four steps. First, the boundary curves of the surface are extracted and then discretized into a series of points. Second, the discrete points of the boundary curves are mapped onto an x1x2 2D plane to form a boundary of a 2D domain. Third, the meshless point generation algorithm for the 2D domain is applied to generate a series of discrete points within the domain. Fourth, all the discrete points on the x1x2 plane are mapped back onto the 3D surface to obtain the meshless discrete point set.
The coordinates of the background points in the 3D domain are calculated as
x i , j , k = x s t a r t + ( i 1 ) Δ x , ( 1 i N x ) y i , j , k = y s t a r t + ( j 1 ) Δ y , ( 1 j N y ) z i , j , k = z s t a r t + ( k 1 ) Δ z , ( 1 k N z )
where the subscript start represents the background points at the start; ∆x, ∆y, and ∆z represent the spacings between two neighboring background points along the x, the y, and the z directions, respectively; Nx, Ny, and Nz represent the total number of rows of all background points generated along the x, the y, and the z directions, respectively; and the subscripts i, j, and k denote the orders of the rows with certain background points along the x, the y, and the z directions, respectively.
After generating the background points, the spatial relationship between the background points and the boundary points must be determined. A vector dot product method is introduced in this work to determine the attribute of the background points and thereby accelerate the generation of the discrete points. Considering a computational domain partially composed of a hemispherical surface, as shown in Figure 3, as an example, it is defined that the right side and the left side of the hemispherical surface are the outside and the inside of the computational domain, respectively. Point A is a point on the surface, whose normal vector is n , and point B beside the surface is a background Cartesian point to be evaluated. Point A is strictly selected by following two principles. First, point A is inside the 3D background points where point B is located. Second, point A is the closest surface boundary point to point B.
Considering the vector pointing from point A to point B as A B , the attribute of the background Cartesian point B can be determined as
A B n > 0 , In A B n < 0 , Out
where the indicators In and Out represent situations of point B being inside and outside the computational domain, respectively. As each background Cartesian point needs to be evaluated at least once, the corresponding efficiency is very low in the case of a large number of background Cartesian and boundary points. Therefore, the accelerated ray method developed in one of our previous studies [41] is utilized to improve the efficiency of the evaluation of the background Cartesian points and consequently to accelerate the point generation.

2.3.2. SPH Adaptive Algorithm

  • Dynamic refinement and coarsening algorithm
Following the initial discretization of the overall computational domain, the refinement of particles within the region of interest is performed at the beginning of the numerical calculations to achieve high calculation accuracy. Inspired by previous studies on the refinement algorithm [42,43], a parent particle, within the initially and uniformly discretized computational domain, requiring refinement is split into eight child particles by using a pre-defined cubic refinement pattern as shown in Figure 4. The parent particle is located in the center of the cubic comprising the eight child particles. Thus, the distance between two neighboring child particles, Δdc, is determined by
Δ d c = ε Δ d p
where the subscripts c and p, respectively, represent the child particle and the parent particle, Δdp is the initial distance between neighboring parent particles, and ε is the separation ratio ranging from 0 to 1. In addition, the smooth length of the child particle hc is expressed as
h c = ς h p
where hp is the smooth length of the parent particle and ς denotes the smooth length ratio ranging from 0 to 1. The two refinement parameters, ε = 0.5 and ς = 0.5, are utilized in this work to warrant the same number of neighboring particles before and after refinement.
The physical properties of the child particle, the mass mc and the velocity uc, are derived from the parent particle as
m c = λ c m p
u c = u p
where λc is the refinement coefficient. This is determined based on the refinement pattern and set as 0.125 to achieve a cubic refinement pattern in this work. Consequently, the mass, the energy, the linear momentum, and the angular momentum conservation equations are expressed as
m p = c = 1 8 λ c m c
1 2 m p u p 2 = c = 1 8 m c u c 2
m p u p = c = 1 8 m c u c
r p × m p u p = c = 1 8 r c × m c u c
To further increase the efficiency of the SPH adaptive algorithm, a coarsening algorithm can be reasonably developed to dynamically remove the numerous child particles that move out of the interested region. Inspired by our previous study on the dynamic coarsening algorithm for a 2D computational domain [41], this work further improves the dynamic refining and coarsening algorithm to render it feasible for a 3D computational domain. Once a parent particle enters the region of interest (refinement region), it is split into eight child particles. To conserve all parameters, the parent particle is inactivated and the eight affiliated child particles are simultaneously activated. The inactivated parent particle passively moves with the flow field and all the parameters of the parent particle are not involved in any calculations. Once the parent particle exits the refinement region, it is immediately activated to be involved in the calculations and the eight affiliated child particles are simultaneously removed from the calculation domain. In general, the improved dynamic refinement and coarsening algorithm can simultaneously realize high accuracy, high efficiency, and better conservations during calculations.
  • Adaptive particle refinement (APR) technology
In SPH calculations, two neighboring regions of different resolutions interact with each other through particles in the guard area [41]. The field function f of the particles in the guard area is interpolated from the fluid as
f ( r ) guard   = j = 1 N f ( r j ) W ( r i r j , h ) j = 1 N W ( r i r j , h )
where f ( r ) guard   denotes the approximation of the field function f for particles in the guard area. Interpolations of the field functions are conducted before each time step and the corresponding values are assigned equally to each particle in the guard area using Equation (16). In addition, an enhanced shifting algorithm [41] is also applied to correct the erroneous displacement of particles in the non-guard area and thus further increase the computational accuracy.

3. Results and Discussion

3.1. Numerical Model Validation

A classic study on water entry of a 3D horizontally placed cylinder [44] is conducted to validate the mathematical model and the SPH numerical method developed in this work. As shown in Figure 5, the 3D computational domain contains a large pool, filled with water, and a cylinder, horizontally located upon the pool. The pool is sufficiently large with dimensions of 1 m in length, 0.4 m in width, and 0.6 m in height to ensure that the calculation of the fluid flow is not affected by boundaries. The diameter and the length of the cylinder are 0.05 and 0.2 m, respectively. The densities of the water and the cylinder are 1000 and 1370 kg/m3, respectively. The kinematic viscosity of water is 10−6 m2/s and the initial velocity of the cylinder entering the water is set as 6.22 m/s.
First, grid independence verification is performed by applying four different particle resolutions, including dclddp = 5, 10, 15, and 20. The corresponding time evolutions of the penetration depth y−y0 and the vertical velocity uy of the cylinder entering water achieved from this numerical work and the reported experimental studies [44,45] are summarized in Figure 6. As evident, the numerical variation trends of both parameters gradually converge to the corresponding experimental ones with increasing particle resolution. Considering the numerical results achieved at a particle resolution of dclddp = 20 as an example (Table 1), the errors between the numerical and the reported experimental results at 0.12 s are only 2.11% and 4.15% for the penetration depth y − y0 and the vertical velocity uy, respectively. This comparison result validates that the mathematical model and the SPH numerical method developed in this work can effectively give rise to accurate convergence results.
As calculated, the total number of particles is as high as 1.7 × 107 if the particle resolution of dclddp = 20 is applied for the whole 3D computational domain, which greatly reduces the computational efficiency. Thus, three adaptive refinement areas with three different particle resolutions of dclddp = 20, 10, and 5 are assigned outward from the cylinder to skillfully reduce the total number of particles and thus improve the computational efficiency. In particular, the distance from the boundary of the area of the highest particle resolution to the boundary of the cylinder is controlled at 30Δdp,max to capture the force of the fluid on the cylinder more accurately. Here, Δdp,max is the initial distance between the two neighboring particles within the area of the highest particle resolution, Based on this adaptive technique, the total number of particles is reduced to 1.4 × 106, which is approximately 92% less than that generated under the condition of the uniformly distributed particle resolution of dclddp = 20. Therefore, the adaptive technique significantly reduces the usage of the computer memory and thus improves computational efficiency.
Figure 7 shows the evolutions of the velocities within the middle plane along the z direction, as achieved from this work and the reported experimental work [44]. The results showed significant consistency. In the initial stage of the cylinder entering water, the fluid splashes around the cylinder, forming a thin water curtain, as shown in Figure 7b. As the cylinder continues to sink, an open cavity forms and gradually grows above the cylinder. Simultaneously, the thin water curtain also gradually expands to stand vertically. Owing to the large initial splash velocity of the fluid, the tail of the thin water curtain eventually separates from the thin water curtain, forming a number of small splashed droplets, as shown in Figure 7c. As the cavity continues to expand, the free surface of the fluid forms four obvious cross lines along four corners of the cylinder when observed from the top of the computational domain. Two pairs of cross lines at the left and the right ends of the cylinder shrink inward faster than the other parts of the free surface. Two cross lines at each end approach each other and eventually merge to form a new cross line at the corresponding end, as shown in Figure 7c–f. The two newly formed cross lines continue to approach each other along the axial direction at a faster speed, as shown in Figure 7f–h. The phenomenon observed above is typically referred to as cavity inversion [46]. The remarkable consistency of the results achieved from this work and the reported experimental work proves that the mathematical model and the SPH numerical method developed in this work can accurately and efficiently capture all the physical phenomena of the water entry problems of 3D entities.

3.2. Water Entry of 3D Torpedo

This section systematically investigates the water entry process of a 3D torpedo at different entry velocities and entry angles. The corresponding schematic and detailed parameters are shown in Figure 8. The length, the width, and the height of the pool filled with water are 8, 8, and 4 m, respectively. The length and the diameter of the torpedo are 2.79 and 0.32 m, respectively. The center of mass of the torpedo is located at a distance of 1.18 m from its head. The density of the torpedo is assumed as 1000 kg/m3 and the rotational inertias of the torpedo are set as Ix = Iz = 88.01 kg∙m2 and Iy = 2.31 kg∙m2. The torpedo is initially placed above the water with its head placed tangentially to the surface of the water. In addition, its center of mass is fixed at (x0, y0, z0). At the beginning of the water entry, the entry velocity of the torpedo is set along the same direction as the axis of the torpedo. Therefore, the entry angle β formed between the entry velocity and the water surface is the same as the attitude angle θ of the torpedo formed between the axis of the torpedo and the water surface.

3.2.1. Effect of Entry Velocity

Three different entry velocities of 5, 10, and 15 m/s, at the same entry angle of β = 90°, are applied to investigate the effects of the entry velocity on the water entry of a torpedo. A total of 1.52 × 107 initialization particles are employed as the optimum strategy of particle discretization along with the adaptive refinement performed for the area around the torpedo. Figure 9 and Figure 10 summarize the evolution of the velocity fields within the middle plane of the calculation domain along the z direction and the development processes of the cavities during the water entry of the torpedo at three different entry velocities, respectively. As evident, three stages of evolution are observed.
The first stage is the impact stage. When a torpedo starts to collide with the water surface in a very short time, intensive compression waves are induced and dissipate rapidly to the surroundings. Meanwhile, the head of the torpedo is significantly impacted by the instantaneous impact load, which may result in considerable damage to the torpedo head. Typically, the impact process is accompanied by the exchange of momentum between the torpedo and the surrounding fluid near the impact surface. Consequently, the surrounding fluid, obtaining momentum from the torpedo, is squeezed away from its original location and splashes around at a high velocity, as shown in Figure 9a–c for t = 100, 60, and 40 ms, respectively. In addition, with a gradual increase in the entry velocity from 5 to 15 m/s, the surrounding fluid obtains growing momentum and thus splashes away at an increasing speed.
The second stage is the cavity formation stage. After the impact stage, the torpedo continues to dive downward. Thus, more fluid around the head of the torpedo is squeezed out and gradually forms a growing V-shaped cavity, as shown in Figure 10. In addition, it is also found that a larger entry velocity results in a larger cavity formed at the same depth of the water domain.
The third stage is the cavity closure stage. As the torpedo further enters into the water, the fluid, forming the free surface of the cavity around the torpedo, moves radially toward the torpedo under the hydrostatic force. Consequently, it gradually closes the space formed by the cavity and the torpedo surfaces as shown in Figure 10. For a low entry velocity, the cavity is small and thus is closed earlier. As for the case of a low entry velocity of 5 m/s (Figure 10a), the cavity starts to close from the head to the rear of the torpedo at t = 100 ms, and it is fully closed at t = 370 ms. Subsequently, the free surface of the fluid around the previous cavity continues to be squeezed upward, forming new splashing droplets (Figure 9a and Figure 10a). As shown in Figure 10b, when the entry velocity increases to 10 m/s, the cavity requires approximately 170 ms to close, and the area where the cavity closes is from the middle part instead of the head of the torpedo. In addition, the induced droplets with a series of higher splashing heights are observed (Figure 9b and Figure 10b). When the entry velocity continues to increase to 15 m/s (Figure 10c), the cavity does not close until 202 ms and the induced droplets with a series of highest splashing heights are captured (Figure 9c and Figure 10c). In general, when the entry velocity is low, the free surface of the cavity closes with the torpedo earlier at the front part of the torpedo. In addition, a small quantity of splashing droplets is also generated and they occupy a small space above the water domain. With increasing entry velocity, the cavity size increases significantly, and numerous droplets are simultaneously generated and splash above the water domain.
Figure 11 summarizes the variations in the impact load projecting onto the y-axis Fy, the velocity uy, and the vertical displacement of the center of mass y − y0 of the torpedo with time when the torpedo vertically enters the water at different speeds of 5, 10, and 15 m/s. As shown in Figure 11a, once the torpedo enters the water at a certain entry velocity, it is subjected to an instantaneous impact load along the positive y-direction, Fy, which increases to two peak values in a short period and stabilizes afterward. In addition, it is also observed that a higher entry velocity results in a higher peak and stabilized impact loads. In addition, the gravity force G is calculated as 1.848 kN. Subsequently, the vertical resultant forces and the corresponding accelerations, that occur around the moment of the first peak impact loads, point to the positive direction of the y-axis. Therefore, the absolute values of the velocities decrease at the beginning, as shown in Figure 11b. Although the impact loads along the positive y-direction slightly fluctuate, the regions of Fy < G between two peak loads occur in a very short time. Therefore, the fluctuations in the velocities can be reasonably neglected during the corresponding period. Figure 11b also shows that different stabilized impact loads, induced by different entry velocities, induce different stabilized accelerations and variation trends of velocities. For a low entry velocity of 5 m/s, as the vertical stabilized impact load acting on the torpedo is less than the gravity force, the resultant force along the vertical direction is in the same direction as that of the velocity. Therefore, the torpedo accelerates and the absolute value of velocity of the torpedo linearly increases during the stabilization period. For an entry velocity of 10 m/s, as the vertical stabilized impact load is almost equal to the gravity force, the resultant force and the acceleration along the vertical direction are zeros and thus the corresponding velocity is maintained at a relatively stable value during the stabilization period. Contrary to the findings observed from the case of a low entry velocity of 5 m/s, the large entry velocity of 15 m/s induces a high stabilized impact load, which is greater than the gravity force. Thus, the corresponding resultant force is in the opposite direction from the directions of the gravity force and the velocity. Therefore, the absolute value of the velocity linearly decreases during the stabilization period. Although the variation trends of the velocities are opposite for different entry velocities of 5 and 15 m/s, the displacements of the center of mass are mainly dominated by the initial entry velocity. Therefore, the displacements of the center of mass at the initial entry velocities of 15 and 5 m/s increase the fastest and the slowest, respectively, as shown in Figure 11b.

3.2.2. Effect of Entry Angle

In addition to the case of the entry angle of β = 90° under 10 m/s studied above, two more entry angles of β = 75° and 60° under the same entry velocity are also investigated to determine the effects of the entry angle on the water entry of the torpedo. Figure 12 and Figure 13 summarize the evolution of the velocity fields within the middle plane of the calculation domain along the z direction and the development processes of the cavities during the water entry of the torpedo at different entry angles, respectively.
Once the torpedo collides with water in a very short time, it causes the nearby fluid to splash radially outward with high velocities, as seen from the figures of t = 55 ms in Figure 12a and t = 55 ms in Figure 12b. However, in contrast to the symmetric splash curtain generated by the vertical water entry of the torpedo, the velocity along the x-positive component of the torpedo ux, occurring when the torpedo enters water at angles of β = 75° or 60°, results in a larger impact force on the fluid in front of the torpedo. This results in an asymmetric splash curtain, as shown in Figure 12 and Figure 13. Correspondingly, an asymmetric cavity is also simultaneously generated, which is larger at the front of the torpedo owing to the higher momentum transferred to the fluid there. As the torpedo moves downward, the fluid gradually flows back to contact with the torpedo from the bottom to the upper part of the torpedo. In addition, the cavity closes faster at the back of the torpedo than at the front of the torpedo. Furthermore, the cavity at the back of the torpedo closes faster for a smaller entry velocity.
The variations in the impact loads, the velocities, and the displacements of the center of mass of the torpedo with time when the torpedo enters water at different entry angles of β = 60°, 75°, and 90° under the same entry velocity of 10 m/s are drawn in Figure 14. Compared to the impact load of the x-component of zero at the entry angle of β = 90° depicted in Figure 14a, entry angles of β = 60° and 75° lead to non-zero impact loads of the x-component. These two non-zero horizontal impact loads initially point to the negative direction of the x-axis, impeding the horizontal movement of the torpedo. In addition, a larger impact load of the x-component occurs at the smaller entry angle of β = 60° in the earlier stage of the torpedo entering the water. Correspondingly, the velocity of the x-component at the smallest entry angle of β = 60° is largest and it decreases the most sharply in the earlier stage, as shown in Figure 14c. When the torpedo continues to enter water more deeply in the later stage, the impact loads of the x-component at entry angles of β = 60° and 75° gradually approach zero, as shown in Figure 14a; then, they increase toward the direction of the positive x-axis owing to the momentum introduced by the backflow of the fluid, as observed in Figure 13. Thus, the velocities of the x-component at those two entry angles slightly increase in the later stage of the water entry of the torpedo, as shown in Figure 14c.
The torpedo is initially subjected to fluctuating vertical impact loads along the positive direction of the y-axis at all three entry angles and the larger entry angle yields a higher vertical impact load, as shown in Figure 14b. Correspondingly, the torpedo achieves higher downward velocity for a larger entry angle, as shown in Figure 14d. All those vertical impact loads increase slowly after the fluctuations. Considering the gravity force, the vertical resultant forces at those three entry angles are mostly along the opposite direction with the direction of the vertical velocity during the fluctuation period. Therefore, the absolute values of the velocities decrease first. In the later stage after t = 0.02 s, the variation trends of the vertical velocities at those three entry angles are slightly different, as shown in Figure 14b. The absolute value of the vertical velocity of the torpedo entering the water at β = 90° slowly decreases after t = 0.02 s as the induced vertical impact load is always slightly larger than the gravity force of the torpedo during this period. However, for entry angles of β = 60° and 75°, the corresponding vertical impact loads increase from the loads being smaller than the gravity force to the loads being larger than the gravity force. Thus, the resultant forces at those two entry angles in the later stages first exhibit the same direction with the vertical velocities and then vary to the opposite direction. Correspondingly, the absolute values of the velocities at entry angles of β = 60° and 75° first increase and thereafter decrease in the later stages, as shown in Figure 14d.
Variations in the horizontal displacement x-x0 and the vertical displacement y − y0 of the center of mass of the torpedo with time are shown in Figure 14e,f, respectively. In general, the displacements along the different dimensions are mainly determined by the corresponding initial velocities as the velocities at different dimensions and different entry angles fluctuate within relatively small ranges during such a short period from 0 s to 0.2 s. At the same entry velocity, the smaller entry angle results in a larger absolute value of the horizontal velocity and a smaller absolute value of the vertical velocity according to the principle of velocity orthogonal projection. Therefore, the smaller entry angle ultimately induces a larger absolute value of the horizontal displacement and, simultaneously, a smaller absolute value of the vertical displacement.
Figure 15 shows the variations in the moment about the z-axis mz, the angular velocity ωz, and the variation in the attitude angle of the torpedo θθ0 with time when the torpedo enters the water at different entry angles of 60°, 75°, and 90° under the same entry velocity of 10 m/s. When the torpedo vertically enters water at β = 90°, there is no rotating movement of the torpedo observed, owing to the relatively uniform force acting around the torpedo body. Therefore, the angular velocity and the variation in the attitude angle at β = 90° are always maintained at zero. However, for the torpedo entering the water at the other two entry angles, similar and remarkable fluctuations in different parameters of the torpedo are observed. Owing to the induced complex impact loads acting on the torpedo body at different periods of the water entry, as shown in Figure 14a,b, the moments of the torpedo first decrease to negative values and then increase to positive values and finally decrease again. Affected by the fluctuating moments, the corresponding angular velocities and the variations in the attitude angles exhibit similar variation trends with those of the moments. In general, as the entry angle decreases, the fluctuation in the moment acting on the torpedo becomes more significant, thereby resulting in greater changes in the angular velocity and the attitude angle of the torpedo.

4. Conclusions

This work systematically investigates the complex water entry process of the 3D torpedo by applying an improved adaptive SPH algorithm. The initial flow field is generated by an algorithm of 3D meshless point generation. Based on adaptive particle refinement (APR) technology, the dynamic refinement/coarsening algorithm of particles is improved, avoiding the flow field noise in the dynamic refinement process and maintaining the mass conservation of the flow field. The guard area adopted in APR technology facilitates the realization of the coupling of two regions of different resolutions, avoiding the numerical error caused by the interaction between particles of different smooth lengths. In addition, the enhanced shifting algorithm is used to ensure the uniform distribution of particles during calculations, effectively preventing non-physical oscillations of the flow field and improving the accuracy of local SPH numerical simulations.
Based on the improved SPH method, a parametric study on the water entry of the torpedo is conducted from two aspects: entry velocity and entry angle. When the torpedo vertically enters water at different velocities, all the induced impact loads acting on the head of the torpedo fluctuate significantly with two peak values in the initial stage and thereafter stabilize in the later stage. With the growth in the entry velocity, the impact load acting on the torpedo, the entry depth of the torpedo, the splash height of the droplets, and the size of the cavity generated around the torpedo increase. When the torpedo enters the water at different entry angles under the same initial entry velocity, both the vertical and the horizontal movement of the torpedo are observed, leading to increased and more complex variations in parameters along both the x- and the y-axes. First, the torpedo entering water at angles except β = 90° results in one side of the cavity being larger than another side, forming an asymmetric cavity. Second, with a decrement in the entry angle, the fluctuation in the horizontal impact load becomes more significant, the horizontal velocity increases, and the corresponding horizontal displacement of the center of mass of the torpedo increases, throughout the entire water entry process of the torpedo. On the contrary, opposite conclusions about variations in the fluctuation in the vertical impact load, the absolute value of the vertical velocity, and the absolute value of the vertical displacement of the center of mass of the torpedo with decreasing entry angle are observed. Third, the moment, the angular velocity, and the variation in the attitude angle of the torpedo obliquely entering water initially decrease to negative values, then increases to positive values, and finally decreases again. In addition, the decreasing entry angle increases fluctuations in the moment, the angular velocity, and the attitude angle of the torpedo.
This work systematically investigates the complex nonlinear fluid–structure coupling performance of the water entry of the torpedo, which provides valuable guidance for the future design of the entry trajectory and the structural bearing capacity of the torpedo. However, there are three challenges for future works. Firstly, when a torpedo enters the water, the head will be subjected to a transient peak load during the impact stage, which may damage the structure. This aspect is not fully investigated along with the water entry process in this work. Therefore, it is necessary to develop an SPH-FEM coupling algorithm to solve the elastic deformation problem of the structure in the future, which will give more complicated guidance on the design of the torpedoes. Secondly, when a torpedo enters the water at a high speed, a low-pressure area will be formed at the junction of the torpedo head and the torpedo body, which may cause local cavitation. The existing SPH method lacks an effective cavitation model. Therefore, it is an important research direction to establish an SPH algorithm of a two-phase flow coupled with a cavitation model. Finally, it is also worth studying how to improve the search efficiency of neighboring particles after adaptive refinement.

Author Contributions

Conceptualization, X.C., J.Z. and W.W.; methodology, W.W. and J.Z.; software, X.C. and W.H.; validation, W.H. and W.L.; formal analysis, W.W., W.L. and Y.J.; investigation, W.H. and Y.J.; resources, W.H. and W.W.; data curation, W.H. and X.C.; writing—original draft preparation, X.C., W.W. and Y.J.; writing—review and editing, Y.J. and W.L.; visualization, W.L.; supervision, X.C. and J.Z.; project administration, Y.J.; funding acquisition, J.Z. and Y.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

The authors would like to thank members of the Key Laboratory of Unsteady Aerodynamics and Flow Control for fruitful discussions.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Zhao, G.; Chen, T.; Guo, Z.; Zhang, W. Experimental study on ballistic stability of truncated cone projectile in high-speed oblique water entry. Ocean. Eng. 2023, 283, 115078. [Google Scholar] [CrossRef]
  2. Tassin, A.; Piro, D.J.; Korobkin, A.A.; Maki, K.J.; Cooker, M.J. Two-dimensional water entry and exit of a body whose shape varies in time. J. Fluids Struct. 2013, 40, 317–336. [Google Scholar] [CrossRef]
  3. Wu, Q.; Wang, K.; Kong, D.; Zhang, J.; Liu, T. Numerical simulation of subsonic and transonic water entry with compressibility effect considered. Ocean. Eng. 2023, 281, 114984. [Google Scholar] [CrossRef]
  4. Iranmanesh, A.; Passandideh, F.M. A three-dimensional numerical approach on water entry of a horizontal circular cylinder using the volume of fluid technique. Ocean. Eng. 2017, 130, 557–566. [Google Scholar] [CrossRef]
  5. Nguyen, V.T.; Park, W.G. Enhancement of Navier–Stokes solver based on an improved volume-of-fluid method for complex interfacial-flow simulations. Appl. Ocean. Res. 2018, 72, 92–109. [Google Scholar] [CrossRef]
  6. Wen, X.L.; Ong, M.C.; Yin, G. On the three-dimensional effects of the water entry of wedges. Appl. Ocean. Res. 2023, 138, 103649. [Google Scholar] [CrossRef]
  7. Dong, L.; Ding, W.; Wei, Z.; Liu, P.; Yao, B.; Lian, L. Numerical study on the water entry of two-dimensional airfoils by BEM. Eng. Anal. Bound. Elem. 2023, 151, 83–100. [Google Scholar] [CrossRef]
  8. Zheng, K.Y.; Zhao, X.Z.; Yan, D.M. Numerical simulation of water entry of two-dimensional structures with complex geometry using a CIP-based model. Appl. Ocean. Res. 2021, 106, 102379. [Google Scholar] [CrossRef]
  9. Zhang, G.; Liang, G.; Yang, X.; Zhang, Z. Numerical investigations on water entry and/or exit problems using a multi-resolution Delta-plus-SPH model with TIC. Ocean. Eng. 2024, 292, 116560. [Google Scholar] [CrossRef]
  10. Wang, Z.B.; Chen, R.; Wang, H.; Liao, Q.; Zhu, X.; Li, S.Z. An overview of smoothed particle hydrodynamics for simulating multiphase flow. Appl. Math. Model. 2016, 40, 9625–9655. [Google Scholar] [CrossRef]
  11. Chen, F.G.; Ge, W. A review of smoothed particle hydrodynamics family methods for multiphase flow. Chin. J. Theor. Appl. Mech. 2021, 53, 2357–2373. [Google Scholar]
  12. Oger, G.; Doring, M.; Alessandrini, B.; Ferrant, P. Two-dimensional SPH simulations of wedge water entries. J. Comput. Phys. 2006, 213, 803–822. [Google Scholar] [CrossRef]
  13. Gong, K.; Liu, H.; Wang, B. Water entry of a wedge based on SPH model with an improved boundary treatment. J. Hydrodyn. 2009, 21, 750–757. [Google Scholar] [CrossRef]
  14. Shao, S. Incompressible SPH simulation of water entry of a free-falling object. Int. J. Numer. Methods Fluids 2009, 59, 91–115. [Google Scholar] [CrossRef]
  15. Yang, X.F.; Peng, S.L.; Liu, M.B. Simulation of water entry with smoothed particle hydrodynamics method. Chin. J. Comput. Phys. 2011, 28, 523–528. [Google Scholar]
  16. Panciroli, R.; Abrate, S.; Minak, G.; Zucchelli, A. Hydroelasticity in water-entry problems: Comparison between experimental and SPH results. Compos. Struct. 2012, 94, 532–539. [Google Scholar] [CrossRef]
  17. Koukouvinis, P.K.; Anagnostopoulos, J.S.; Papantonis, D.E. Simulation of 2D wedge impacts on water using the SPH–ALE method. Acta Mech. 2013, 224, 2559–2575. [Google Scholar] [CrossRef]
  18. Farsi, M.; Ghadimi, P. Simulation of 2D symmetry and asymmetry wedge water entry by smoothed particle hydrodynamics method. J. Braz. Soc. Mech. Sci. Eng. 2015, 37, 821–835. [Google Scholar] [CrossRef]
  19. Yan, R. Research on Some Problems in Structure Impact with Water Using SPH Method. Ph.D. Thesis, Northwestern Polytechnical University, Xi’an, China, 2016. [Google Scholar]
  20. Zhou, J.; Xu, S.L. SPH simulation on the behaviors of projectile water entry. Explos. Shock Waves 2016, 36, 326–332. [Google Scholar]
  21. Nair, P.; Tomar, G.A. study of energy transfer during water entry of solids using incompressible SPH simulations. Sādhanā 2017, 42, 517–531. [Google Scholar] [CrossRef]
  22. Sun, P.; Zhang, A.M.; Marrone, S.; Ming, F. An accurate and efficient SPH modeling of the water entry of circular cylinders. Appl. Ocean. Res. 2018, 72, 60–75. [Google Scholar] [CrossRef]
  23. Shao, J.R.; Yang, Y.; Gong, H.F.; Liu, M.B. Numerical simulation of water entry with improved SPH method. Int. J. Comput. Methods 2019, 16, 1846004. [Google Scholar] [CrossRef]
  24. Pei, Z.B.; Yang, X.F.; Wang, X.L. SPH simulation of water entry of a slender body. Unmanned Syst. Technol. 2022, 5, 40–49. [Google Scholar]
  25. Wu, Y.F. Study on Water Entry Process of Slender Body in Wave Environment. Master’s Thesis, Nanjing University of Aeronautics and Astronautics, Nanjing, China, 2022. [Google Scholar]
  26. Goyal, V.K.; Huertas, C.A.; Vasko, T.J. Smooth Particle Hydrodynamics for Bird-Strike Analysis Using LS-DYNA. Am. Trans. Eng. Appl. Sci. 2013, 2, 83–107. [Google Scholar]
  27. Chiron, L.; Oger, G.; De Leffe, M.; Le Touzé, D. Analysis and improvements of Adaptive Particle Refinement (APR) through CPU time, accuracy and robustness considerations. J. Comput. Phys. 2018, 354, 552–575. [Google Scholar] [CrossRef]
  28. Li, D.M.; Xu, Y.N.; Li, L.L.; Lu, H.J.; Ling, B.A.I. Tracking methods for free surface and simulation of a liquid droplet impacting on a solid surface based on SPH. J. Hydrodyn. 2011, 23, 447–456. [Google Scholar] [CrossRef]
  29. Peng, Y.X.; Zhang, A.M.; Ming, F.R. On the comparison of particle regeneration technique and volume adaptive scheme in the compressible flow based on smoothed particle hydrodynamics. J. Hydrodyn. 2021, 376, 113653. [Google Scholar]
  30. Wendland, H. Piecewise polynomial, positive definite and compactly supported radial functions of minimal degree. Adv. Comput. Math. 1995, 4, 389–396. [Google Scholar] [CrossRef]
  31. Liu, X.; Lee, C.K.; Fan, S.C. On using enriched cover function in the Partition-of-unity method for singular boundary-value problems. Comput. Mech. 2002, 29, 212–225. [Google Scholar] [CrossRef]
  32. Ni, X.Y.; Feng, W.B.; Huang, S.C.; Hu, Z.J.; Liu, Y. An SPH Wave-Current Flume using Open Boundary Conditions. J. Hydrodyn. 2020, 32, 12. [Google Scholar] [CrossRef]
  33. He, M.; Gao, X.F.; Xu, W.H. Numerical simulation of wave-current interaction using the SPH method. J. Hydrodyn. 2018, 30, 535–538. [Google Scholar] [CrossRef]
  34. Zhao, Z.X.; Hong, Y.; Gong, Z.X.; Liu, H. Numerical analysis of cavity deformation of oblique water entry using a multi-resolution two-phase SPH method. Ocean. Eng. 2023, 269, 113456. [Google Scholar] [CrossRef]
  35. Rosić, N.M.; Kolarević, M.B.; Savić, L.M.; Đorđević, D.M.; Kapor, R.S. Numerical modelling of supercritical flow in circular conduit bends using SPH method. J. Hydrodyn. 2017, 29, 9. [Google Scholar] [CrossRef]
  36. Monaghan, J.J.; Pongracic, H. Artificial viscosity for particle methods. Appl. Numer. Math. 1985, 1, 187–194. [Google Scholar] [CrossRef]
  37. Colagrossi, A.; Antuono, M.; Le Touzé, D. Theoretical considerations on the free-surface role in the smoothed-particle-hydrodynamics model. Phys. Rev. E 2009, 79, 056701. [Google Scholar] [CrossRef] [PubMed]
  38. Monaghan, J.J.; Kos, A.; Issa, N. Fluid motion generated by impact. J. Waterw. Port Coast. Ocean. Eng. 2003, 129, 250–259. [Google Scholar] [CrossRef]
  39. Hormann, K.; Agathos, A. The point in polygon problem for arbitrary polygons. Comput. Geom. 2001, 20, 131–144. [Google Scholar] [CrossRef]
  40. Pan, W.; Wu, H.; Li, T.R.; Gao, B.Q. Grid generation on free-form surface based on surface flattening. J. Zhejiang Univ. (Eng. Sci.) 2016, 50, 1973–1979. [Google Scholar]
  41. Wu, W. Research on Meshless Point Generation and Algorithm Based on SPH. Master Thesis, Nanjing University of Aeronautics and Astronautics, Nanjing, China, 2023. [Google Scholar]
  42. Feldman, J.; Bonet, J. Dynamic refinement and boundary contact forces in SPH with applications in fluid flow problems. Int. J. Numer. Methods Eng. 2007, 72, 295–324. [Google Scholar] [CrossRef]
  43. Yaidel, R.L.; Dirk, R. Dynamic refinement for fluid flow simulations with SPH. In Proceedings of the International Conference on Particle-Based Methods (PARTICLES)—II International Conference on Particle-Based Methods: Fundamentals and applications (PARTICLES 2011), Barcelona, Spain, 26–28 October 2011. [Google Scholar]
  44. Wei, Z.; Hu, C. An experimental study on water entry of horizontal cylinders. J. Mar. Sci. Technol. 2014, 19, 338–350. [Google Scholar] [CrossRef]
  45. Sun, P.N. Study on SPH Method for the Simulation of Object-Free Surface Interactions. Ph.D. Thesis, Harbin Engineering University, Harbin, China, 2018. [Google Scholar]
  46. Enriquez, O.R.; Peters, I.R.; Gekle, S.; Schmidt, L.E.; Versluis, M.; van der Meer, D.; Lohse, D. Collapse of nonaxisymmetric cavities. Phys. Fluids 2010, 22, 091104. [Google Scholar] [CrossRef]
Figure 1. Flowchart of the spatial discretization of a 3D space.
Figure 1. Flowchart of the spatial discretization of a 3D space.
Applsci 14 04441 g001
Figure 2. Generation of discrete points using the mapping method.
Figure 2. Generation of discrete points using the mapping method.
Applsci 14 04441 g002
Figure 3. Schematic of the vector dot product method used for determining the spatial relationship between boundary point A and background Cartesian point B. n is the normal vector of the surface.
Figure 3. Schematic of the vector dot product method used for determining the spatial relationship between boundary point A and background Cartesian point B. n is the normal vector of the surface.
Applsci 14 04441 g003
Figure 4. Schematic of the cubic refinement pattern for a 3D computational domain. The black and the red dots denote the parent and the child particles, respectively. Δdc and Δdp are the distance between two neighboring child particles and the distance between two neighboring parent particles, respectively.
Figure 4. Schematic of the cubic refinement pattern for a 3D computational domain. The black and the red dots denote the parent and the child particles, respectively. Δdc and Δdp are the distance between two neighboring child particles and the distance between two neighboring parent particles, respectively.
Applsci 14 04441 g004
Figure 5. Three-dimensional computational domain of water entry of a horizontally placed cylinder.
Figure 5. Three-dimensional computational domain of water entry of a horizontally placed cylinder.
Applsci 14 04441 g005
Figure 6. Time evolutions of the penetration depth y − y0 (a) and vertical velocity uy (b) of the cylinder achieved from this numerical work and the reported experimental studies [44,45].
Figure 6. Time evolutions of the penetration depth y − y0 (a) and vertical velocity uy (b) of the cylinder achieved from this numerical work and the reported experimental studies [44,45].
Applsci 14 04441 g006
Figure 7. Evolutions of the velocities, within the middle plane along the z direction, as achieved from this work and the reported experimental work [44].
Figure 7. Evolutions of the velocities, within the middle plane along the z direction, as achieved from this work and the reported experimental work [44].
Applsci 14 04441 g007
Figure 8. Schematic of (a) the torpedo and (b) the calculation domain of the water entry of the torpedo.
Figure 8. Schematic of (a) the torpedo and (b) the calculation domain of the water entry of the torpedo.
Applsci 14 04441 g008
Figure 9. Evolution of velocity fields of the water domain induced by the torpedo vertically entering water at different speeds of (a) 5 m/s, (b) 10 m/s, and (c) 15 m/s.
Figure 9. Evolution of velocity fields of the water domain induced by the torpedo vertically entering water at different speeds of (a) 5 m/s, (b) 10 m/s, and (c) 15 m/s.
Applsci 14 04441 g009
Figure 10. Development of the water cavity induced by the torpedo vertically entering the water at different speeds of (a) 5 m/s, (b) 10 m/s, and (c) 15 m/s. As for the cavity development at a certain entry velocity, the figures within the upper row show the cross-sectional view of the water phase within the middle plane of the computational domain along the z direction and those within the bottom row show the corresponding 3D views of the free surface of the water generated around the torpedo.
Figure 10. Development of the water cavity induced by the torpedo vertically entering the water at different speeds of (a) 5 m/s, (b) 10 m/s, and (c) 15 m/s. As for the cavity development at a certain entry velocity, the figures within the upper row show the cross-sectional view of the water phase within the middle plane of the computational domain along the z direction and those within the bottom row show the corresponding 3D views of the free surface of the water generated around the torpedo.
Applsci 14 04441 g010
Figure 11. Variations in (a) the impact load projecting onto the y-axis Fy and (b) the velocity uy and the vertical displacement of the center of mass y − y0 of the torpedo with time when the torpedo vertically enters the water at different speeds of 5, 10, and 15 m/s.
Figure 11. Variations in (a) the impact load projecting onto the y-axis Fy and (b) the velocity uy and the vertical displacement of the center of mass y − y0 of the torpedo with time when the torpedo vertically enters the water at different speeds of 5, 10, and 15 m/s.
Applsci 14 04441 g011
Figure 12. Evolution of velocity fields of the water domain induced by the torpedo entering the water at different entry angles of (a) β = 75° and (b) β = 60° under the same entry velocity of 10 m/s.
Figure 12. Evolution of velocity fields of the water domain induced by the torpedo entering the water at different entry angles of (a) β = 75° and (b) β = 60° under the same entry velocity of 10 m/s.
Applsci 14 04441 g012
Figure 13. Development of the cavity within the water domain induced by the torpedo entering the water at different entry angles of (a) β = 75° and (b) β = 60° under the same entry velocity of 10 m/s. As for the cavity development at a certain entry angle, the figures within the upper row show the cross-sectional view of the water phase within the middle plane of the computational domain along the z direction and those within the bottom row show the corresponding 3D views of the free surface of the water generated around the torpedo.
Figure 13. Development of the cavity within the water domain induced by the torpedo entering the water at different entry angles of (a) β = 75° and (b) β = 60° under the same entry velocity of 10 m/s. As for the cavity development at a certain entry angle, the figures within the upper row show the cross-sectional view of the water phase within the middle plane of the computational domain along the z direction and those within the bottom row show the corresponding 3D views of the free surface of the water generated around the torpedo.
Applsci 14 04441 g013
Figure 14. Variations in (a) the impact load projecting onto the x-axis Fx, (b) the impact load projecting onto the y-axis Fy, (c) the velocity projecting onto the x-axis ux, (d) the velocity projecting onto the y-axis uy, (e) the horizontal displacement xx0, and (f) the vertical displacement yy0 of the center of mass of the torpedo with time when the torpedo enters water at different entry angles of β = 60°, 75°, and 90° under the same entry velocity of 10 m/s.
Figure 14. Variations in (a) the impact load projecting onto the x-axis Fx, (b) the impact load projecting onto the y-axis Fy, (c) the velocity projecting onto the x-axis ux, (d) the velocity projecting onto the y-axis uy, (e) the horizontal displacement xx0, and (f) the vertical displacement yy0 of the center of mass of the torpedo with time when the torpedo enters water at different entry angles of β = 60°, 75°, and 90° under the same entry velocity of 10 m/s.
Applsci 14 04441 g014
Figure 15. Evolutions of the moment about the z-axis mz, the angular velocity ωz, and the variation in the attitude angle θθ0 of the torpedo when the torpedo enters the water at different entry angles of β = 60°, 75°, and 90° under the same entry velocity of 10 m/s.
Figure 15. Evolutions of the moment about the z-axis mz, the angular velocity ωz, and the variation in the attitude angle θθ0 of the torpedo when the torpedo enters the water at different entry angles of β = 60°, 75°, and 90° under the same entry velocity of 10 m/s.
Applsci 14 04441 g015
Table 1. Errors between the numerical results achieved at a particle resolution of dclddp = 20 and the reported experimental results.
Table 1. Errors between the numerical results achieved at a particle resolution of dclddp = 20 and the reported experimental results.
yy0/muy/m∙s1
Exp. Wei and Hu [44] −0.285\
Sun [45] \−1.421
SPH−0.291−1.362
Error (%)2.114.15
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Cai, X.; Wu, W.; Han, W.; Li, W.; Zhang, J.; Jiao, Y. Study on Water Entry of a 3D Torpedo Based on the Improved Smoothed Particle Hydrodynamics Method. Appl. Sci. 2024, 14, 4441. https://doi.org/10.3390/app14114441

AMA Style

Cai X, Wu W, Han W, Li W, Zhang J, Jiao Y. Study on Water Entry of a 3D Torpedo Based on the Improved Smoothed Particle Hydrodynamics Method. Applied Sciences. 2024; 14(11):4441. https://doi.org/10.3390/app14114441

Chicago/Turabian Style

Cai, Xiaowei, Wei Wu, Wenji Han, Wenjie Li, Jun Zhang, and Yanmei Jiao. 2024. "Study on Water Entry of a 3D Torpedo Based on the Improved Smoothed Particle Hydrodynamics Method" Applied Sciences 14, no. 11: 4441. https://doi.org/10.3390/app14114441

APA Style

Cai, X., Wu, W., Han, W., Li, W., Zhang, J., & Jiao, Y. (2024). Study on Water Entry of a 3D Torpedo Based on the Improved Smoothed Particle Hydrodynamics Method. Applied Sciences, 14(11), 4441. https://doi.org/10.3390/app14114441

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop