# Research on the Rapid Closing Jet Mechanism of Pistol Shrimp’s Claws Based on Fluid Dynamic Grid

Authors: Wei Wei; Xinyu Quan; Hongchao Cao; Shijie Zhang; Ximing Zhao; Nan Yu; Jie Zhou; Hongxiang Wang; Xuyan Hou
Journal: Mathematical Problems in Engineering
 (2021)

---

## Abstract

The predation behavior of the pistol shrimp is extremely special, and the predation process will produce a huge popping sound, which has caused extensive research by scholars from all over the world. This article carried out a study on the rapid closing jet mechanism of pistol shrimp’s claws. A theoretical model, based on the hydrodynamic characteristics of seawater and the theory of fluid-structure coupling, was proposed for the interaction between the claws and seawater. A simulation model was established using the finite volume software Fluent, and the rapid closing jet mechanism of pistol shrimp’s claws was verified by using fluid dynamic grid. This article studied the influence of different fluid models on the simulation results. The effects of the claws’ closing angular velocity and angular acceleration on the interaction between the claws and seawater were analyzed, which provides a theoretical basis for the development of new underwater kinetic energy weapons.


---

## Body

## 1. Introduction

With the continuous development of science and technology, from shabby wooden fishing boats to battleships with iron bones, from aircraft carrier above water to submarine torpedoes under water, and now many new underwater equipment, the origin of human, the sea, is more and more important. The predation behavior of the pistol shrimp is extremely special, and the predation process will produce a huge popping sound, which has caused extensive research by scholars from all over the world [1–12]. This article carried out a study on the rapid closing jet mechanism of pistol shrimp’s claws. The current research on the pistol shrimp mainly focuses on the following aspects: on the one hand, it focuses on the body structure and evolution process of the pistol shrimp, as well as the basic biological characteristics and living habits [13–16]; on the other hand, it focuses on the shooting of the pistol through the high-speed camera system. The jetting behavior of shrimp’s claws and the accompanying phenomena of pistol shrimp and goby were studied in detail [2, 7, 11, 12]. Another part focuses on the study of the morphological and physiological characteristics of the claws and muscle distribution of the claws [17, 18]; a few scholars studied the evolution process of the claws in detail [9], including the claws joints’ structural differences and evolutionary process [19], and the changes in exercise performance caused by the evolution of the structure of the claws’ joints [20]. The existing problems in the current research are as follows: due to the limitations of the experiment, the flow field has not been analyzed in previous studies. Especially when it comes to the investigation of actual pistol shrimps, there are certain limitations in shrimp handling, experimental environment, and conditions, which inherently limits the applicable measurement technology. High-speed photography becomes difficult to perform because it requires a high frame rate (of the order of 106 fps); lighting and focusing are difficult (animals may move in a less controllable way). The pressure signal recorded by the hydrophone may be oversmoothed or underestimated by the bandwidth of the sensor. In addition, the complexity geometry of shrimps’ claws and the uniqueness of each individual animal hinder systematic and reproducible research. On the other hand, the experimentally replicated claws lack the conditions for cavitation formation; for cavitation to occur, a high-speed moving object (plunger) is required. This kind of plunger is difficult to realize in actual size, and it is difficult to move at the actual closing speed, and the difficulty in experimental technology makes the research results not accurate enough. Existing research methods focus on experimental research, while fluid-structure coupling simulation can improve the efficiency of verification and reduce costs. The fluid-structure coupling simulation of the rapid closing jet mechanism of pistol shrimps’ claws can reduce the limitation of the measurement technology by the uncontrollable movement of live shrimp and it is conducive to systematic and reproducible research. Refer to the reference using CFD method to solve practical engineering problems [21, 22] and the reference on the application of dynamic grid technology in actual engineering problems [23]; based on the cavitation theory and fluid-structure coupling theory, this paper uses the finite volume method software Fluent to construct a model about the seawater contact with the claws and a seawater cavitation model. Carry out a simulation of the claw closure under different working conditions. The effects of different closing angular velocity and different closing angular acceleration on the rapid closing of the claws are analyzed, which provides important support for the study of the rapid closing jet mechanism of pistol shrimps’ claws.


## 2. Materials and Methods

### 2.1. Physical Model

The body of a pistol shrimp is6∼10cm long, with a large claw and a small claw. The large claw is composed of dactylus and pollex. The specific structure is shown in Figure 1. This paper mainly studies the movement of the dactylus of the large claw, so only the model of the large claws needs to be established. At the same time, because the three-dimensional shape of the claws is complicated and difficult to achieve through modeling, only the two-dimensional section of the flow channel is intercepted for modeling. The two-dimensional cross-sectional view of the flow channel and the outline of the large claw are shown in Figures 2(a) and 2(b). MATLAB is used to obtain the boundary contour coordinates of the large claw, and then import to CAD modeling and finally use the model shown in Figure 2(c).Figure 1
Structure of pistol shrimp.Figure 2
(a) Two-dimensional cross section of the channel. (b) Shape of pistol shrimp’s large claw. (c) Two-dimensional model of pistol shrimp’s large claw.
(a)(b)(c)


### 2.2. Hydrodynamic Modeling

The continuity equation for viscous fluid is(1)DρDt+ρ∇⋅v=∂ρ∂t+∇ρv=0.The N-S equation of incompressible viscous fluid is(2)DvDt=∂v∂t+v⋅∇v=f−1p+v∇2v.When the partial pressure of a place in the flow field is lower than the saturated vapor pressure at that place, not only the air dissolved in the water (usually 2%) will escape, but also the water will start to vaporize, forming a lot of cavitation bubbles composed of air and steam in the water. These cavitation bubbles are carried by the water flow and collapse when they reach the high-pressure zone. The whole process of the initiation, growth, and collapse of the cavitation bubbles is called cavitation. The degree of cavitation is usually expressed by the dimensionless cavitation indexσ, which is defined as(3)σ=2P0−PvρV02,where P0 is the environmental pressure, V0 is the exit speed, Pv is the saturated vapor pressure of water, and ρ is the density of water.The cavitation numberσ can be considered as the ratio between the force that inhibits cavitation and the force that promotes cavitation, and it represents whether cavitation appears and the degree of cavitation in the flow field. If the cavitation number when cavitation occurs in water is defined as the initial cavitation number σi, when the cavitation number of the flow field σ>σi, no cavitation occurs; when σ=σi, it is in the critical state of cavitation; when σ<σi, cavitation will occur, and the cavitation area will be enlarged as the cavitation number σ decreases. When σ<1, it is considered that cavitation has occurred, and when σ<0.5, it is considered that stable cavitation has occurred.For an isolated spherical bubble in a static state in the water, it is assumed that the bubble contains only water vapor. In the case of ignoring the gas diffusion in water, the balance condition of the bubble is(4)p=Pv−2SR,where p is the water pressure on the wall around the bubble, S is the surface tension coefficient of water, and R is the spherical bubble radius.In fact, in addition to water vapor in the general bubble, there are some gases originally dissolved in water, which diffuse from the surrounding water into the bubble. At this situation, the static balance equation can be written as(5)p=Pv+Pg−2SR,where Pg is the partial pressure of gas in bubble.IfPv is considered to be a constant, and because the specific heat capacity of the surrounding water is large, the gas mass is very small, and the heat imbalance caused by the diffusion of gas from the water into the bubble will soon be adjusted by the surrounding water. In this way, the temperature of vapor and gas in the bubble can be considered constant. When the change process of R is very fast, it is regarded as an ideal gas adiabatic process, so Pg in equation (5) becomes(6)Pg=Pg0V0V=Pg0R0R3,where Pg0, V0, and R0, respectively, are partial pressure of gas, bubble radius, and bubble volume in an initial state and Pg, V, andR, respectively, are partial pressure of gas, bubble radius, and bubble volume in other states.Substituting equation (6) into (5), the static balance conditions for the adiabatic process of cavitation bubble containing steam and gas can be obtained:(7)p=Pv+P0−Pv+2SRR0R3−2SR.The Rayleigh–Plesset (R-P) equation is an ordinary differential equation that simulates the motion of bubbles in a fluid. At present, most of the cavitation numerical models applied to engineering problems are derived from the R-P equation that only considers the first derivative term of the cavitation bubble radius and the pressure-driven term. This type of model does not consider the influence of the second derivative term, viscosity term, and surface tension term on the growth of bubble radius.The Rayleigh–Plesset (R-P) equation is(8)PB−P∞ρL=RR¨+32R˙2+4ϑLRR˙+2SρLR,where PB is the pressure in the bubble, P∞ is the pressure at infinity from the cavitation bubble, ρL is the density of surrounding liquid, R is the cavitation bubble radius, ϑL is the viscosity of the surrounding liquid, and S is the cavitation surface tension.The second derivative term in the complete R-P equation has an important influence on the growth of the bubble radius and cannot be ignored [24]. The actual fluid also contains noncondensable gas (NCG). The saturated mass concentration of NCG that can be dissolved by pure water at a standard atmospheric pressure is 1.5×103mg/L. This part of the gas acts as a cavitation nucleus in the process of cavitation, which plays an important role in the initiation of cavitation [25].Using multiphase cavitation modeling methods, basic two-phase cavitation models include standard viscous flow equations that use controlled mixture (mixed model) or phase (Eulerian multiphase) transport and conventional turbulence models (k-ε model). In cavitation, liquid-gas mass transfer (evaporation and condensation) is controlled by the vapor transport equation:(9)∂∂tαρv+∇⋅αρvV¯v=Re−Rc,where α is the gas volume fraction, ρv is the gas density, V¯v is the gas velocity, and Re and Rc are the mass transfer source terms, which, respectively, are related to the growth and collapse of the vapor bubble.In formula (9), Re and Rc represent the mass transfer between liquid and gas phases in cavitation. They are based on the Rayleigh-Plesset equation of formula (8), which describes the growth of a single vapor bubble in the liquid. The equation for the vapor volume fraction from equation (9) has a general form:(10)∂∂tαρv+∇⋅αρvV¯v=ρvρlρDαDt.The sources of net quality are as follows:(11)R=ρvρlρdαdt.The vapor volume fraction is related to the number of bubbles per volume of liquid to get the following equation:(12)α=nb4/3πRB31+nb4/3pRB3,where nb is the bubble number density.Substitute equations (11) and (12) into equation (10) to derive the following equation:(13)R=ρvρlρα1−α3RB23Pv−Pρl,RB=α1−α34π1n1/3.where R is the mass transfer rate and RB is the bubble radius.The mass transfer rate in the model is proportional toαv1−αv. In addition, the value of this function is zero when α=0 and α=1 and reaches the maximum value between them. Also, in this model, the only parameter that must be determined is the number of bubbles per unit volume of liquid. If it is assumed that no bubbles are generated or destroyed, the bubble number density will remain unchanged.The final form of the model is as follows:WhenPv>P,(14)Re=ρvρlρα1−α3RB23Pv−Pρl.WhenPv≤P,(15)Rc=ρvρlρα1−α3RB23P−Pvρl.There are three existing turbulence numerical simulation methods: direct numerical simulation (DNS), Reynolds average Navier–Stokes (RANS), and large eddy simulation (LES). DNS is currently limited to small scale and small Re-number simulation because it consumes too much computing resources and is difficult to apply in practical engineering. RANS is the application of statistical theory of turbulence, which is the simulation method commonly used in engineering. The unsteady N-S equation is averaged over time, and the Reynolds stress is assumed based on the theoretical knowledge of turbulence, experimental data, or DNS results; that is, various empirical and semiempirical constitutive relations are assumed, so that the average Reynolds equation of turbulence is closed. According to the different starting points for model processing, turbulence model theory can be divided into two categories: one type introduces the governing equation of the second-order pulsation term to form a second-order moment closed model called the Reynolds stress model; another is the eddy viscosity closed mode based on Boussinesq’s eddy viscosity hypothesis, such as zero-equation model, one-equation model, and two-equation model. Due to the large amount of calculation of the Reynolds stress model, the application occasions are very limited. The eddy viscous closed mode has lower requirements on the computer and at the same time can obtain calculation results that meet engineering requirements.The zero-equation model has a common shortcoming; that is, the turbulence viscosity coefficient only depends on the local flow parameters and has nothing to do with the flow elsewhere, which is inconsistent with experimental observations. On this basis, one-equation model and two-equation model were developed. Among the two-equation model, thek-ε model has certain accuracy and consumes limited computing resources, so it has been widely used. This paper compares the effects of three k-ε models on the experimental results.


### 2.3. Simulation Model

Use the preprocessing software Gambit to mesh the calculation space. Import the model into Gambit and use a rectangular area five times the length and three times the width of the model as the boundary of the computational domain. The area near the claws needs to be focused, so the density of the grid nodes near the claws should be increased. The height of the first cell perpendicular to the body surface is 0.05 mm. This height is chosen to makey+ of most of the cells in contact with the body surface within the effective range of the standard wall function. At the same time, considering the calculation speed, the outer division step of the model is taken as 0.1mm, the boundary and the movable chelate are set after the compliant grid, and the interior is filled with seawater. The model is shown in Figure 3.Figure 3
Grid chart in Gambit.Import the .mesh file exported by Gambit into Fluent, and set the seawater parameters according to the living environment of pistol shrimps, as shown in Table1. Use Fluent for rapid closing simulation, and set the model boundary parameters as shown in Table 2. After consulting the literature and referring to the official Fluent parameter settings, the simulation parameters are shown in Table 3. The movement of the dactylus is controlled by Fluent UDF.Table 1
Seawater parameters.
ParameterModel valueActual valueTemperature (°C)2018–24Salinity5%5%Density (kg/m3)10241024Dynamic viscosity (mPa s)1.1401.140Pressure (Pa)3028201.5∼6.5×105Saturated vapor pressure (Pa)22922291.66Table 2
Model boundary parameters.
Boundary nameBoundary typeDactylusWallPollexWallOutletPressure outlet (302820 Pa)Move1DeformingMove2DeformingSeawaterInteriorTable 3
Preset parameters for CFD simulation.
ParameterPreset valueSolverPressure basedTimeTransientTurbulence modelRNGk-ε modelPressure-velocity couplingPISOTransient formulationFirst-order implicitOther term spatial discretizationFirst-order upwindTime step size1×10−6 (s)Number of time steps500


## 2.1. Physical Model

The body of a pistol shrimp is6∼10cm long, with a large claw and a small claw. The large claw is composed of dactylus and pollex. The specific structure is shown in Figure 1. This paper mainly studies the movement of the dactylus of the large claw, so only the model of the large claws needs to be established. At the same time, because the three-dimensional shape of the claws is complicated and difficult to achieve through modeling, only the two-dimensional section of the flow channel is intercepted for modeling. The two-dimensional cross-sectional view of the flow channel and the outline of the large claw are shown in Figures 2(a) and 2(b). MATLAB is used to obtain the boundary contour coordinates of the large claw, and then import to CAD modeling and finally use the model shown in Figure 2(c).Figure 1
Structure of pistol shrimp.Figure 2
(a) Two-dimensional cross section of the channel. (b) Shape of pistol shrimp’s large claw. (c) Two-dimensional model of pistol shrimp’s large claw.
(a)(b)(c)


## 2.2. Hydrodynamic Modeling

The continuity equation for viscous fluid is(1)DρDt+ρ∇⋅v=∂ρ∂t+∇ρv=0.The N-S equation of incompressible viscous fluid is(2)DvDt=∂v∂t+v⋅∇v=f−1p+v∇2v.When the partial pressure of a place in the flow field is lower than the saturated vapor pressure at that place, not only the air dissolved in the water (usually 2%) will escape, but also the water will start to vaporize, forming a lot of cavitation bubbles composed of air and steam in the water. These cavitation bubbles are carried by the water flow and collapse when they reach the high-pressure zone. The whole process of the initiation, growth, and collapse of the cavitation bubbles is called cavitation. The degree of cavitation is usually expressed by the dimensionless cavitation indexσ, which is defined as(3)σ=2P0−PvρV02,where P0 is the environmental pressure, V0 is the exit speed, Pv is the saturated vapor pressure of water, and ρ is the density of water.The cavitation numberσ can be considered as the ratio between the force that inhibits cavitation and the force that promotes cavitation, and it represents whether cavitation appears and the degree of cavitation in the flow field. If the cavitation number when cavitation occurs in water is defined as the initial cavitation number σi, when the cavitation number of the flow field σ>σi, no cavitation occurs; when σ=σi, it is in the critical state of cavitation; when σ<σi, cavitation will occur, and the cavitation area will be enlarged as the cavitation number σ decreases. When σ<1, it is considered that cavitation has occurred, and when σ<0.5, it is considered that stable cavitation has occurred.For an isolated spherical bubble in a static state in the water, it is assumed that the bubble contains only water vapor. In the case of ignoring the gas diffusion in water, the balance condition of the bubble is(4)p=Pv−2SR,where p is the water pressure on the wall around the bubble, S is the surface tension coefficient of water, and R is the spherical bubble radius.In fact, in addition to water vapor in the general bubble, there are some gases originally dissolved in water, which diffuse from the surrounding water into the bubble. At this situation, the static balance equation can be written as(5)p=Pv+Pg−2SR,where Pg is the partial pressure of gas in bubble.IfPv is considered to be a constant, and because the specific heat capacity of the surrounding water is large, the gas mass is very small, and the heat imbalance caused by the diffusion of gas from the water into the bubble will soon be adjusted by the surrounding water. In this way, the temperature of vapor and gas in the bubble can be considered constant. When the change process of R is very fast, it is regarded as an ideal gas adiabatic process, so Pg in equation (5) becomes(6)Pg=Pg0V0V=Pg0R0R3,where Pg0, V0, and R0, respectively, are partial pressure of gas, bubble radius, and bubble volume in an initial state and Pg, V, andR, respectively, are partial pressure of gas, bubble radius, and bubble volume in other states.Substituting equation (6) into (5), the static balance conditions for the adiabatic process of cavitation bubble containing steam and gas can be obtained:(7)p=Pv+P0−Pv+2SRR0R3−2SR.The Rayleigh–Plesset (R-P) equation is an ordinary differential equation that simulates the motion of bubbles in a fluid. At present, most of the cavitation numerical models applied to engineering problems are derived from the R-P equation that only considers the first derivative term of the cavitation bubble radius and the pressure-driven term. This type of model does not consider the influence of the second derivative term, viscosity term, and surface tension term on the growth of bubble radius.The Rayleigh–Plesset (R-P) equation is(8)PB−P∞ρL=RR¨+32R˙2+4ϑLRR˙+2SρLR,where PB is the pressure in the bubble, P∞ is the pressure at infinity from the cavitation bubble, ρL is the density of surrounding liquid, R is the cavitation bubble radius, ϑL is the viscosity of the surrounding liquid, and S is the cavitation surface tension.The second derivative term in the complete R-P equation has an important influence on the growth of the bubble radius and cannot be ignored [24]. The actual fluid also contains noncondensable gas (NCG). The saturated mass concentration of NCG that can be dissolved by pure water at a standard atmospheric pressure is 1.5×103mg/L. This part of the gas acts as a cavitation nucleus in the process of cavitation, which plays an important role in the initiation of cavitation [25].Using multiphase cavitation modeling methods, basic two-phase cavitation models include standard viscous flow equations that use controlled mixture (mixed model) or phase (Eulerian multiphase) transport and conventional turbulence models (k-ε model). In cavitation, liquid-gas mass transfer (evaporation and condensation) is controlled by the vapor transport equation:(9)∂∂tαρv+∇⋅αρvV¯v=Re−Rc,where α is the gas volume fraction, ρv is the gas density, V¯v is the gas velocity, and Re and Rc are the mass transfer source terms, which, respectively, are related to the growth and collapse of the vapor bubble.In formula (9), Re and Rc represent the mass transfer between liquid and gas phases in cavitation. They are based on the Rayleigh-Plesset equation of formula (8), which describes the growth of a single vapor bubble in the liquid. The equation for the vapor volume fraction from equation (9) has a general form:(10)∂∂tαρv+∇⋅αρvV¯v=ρvρlρDαDt.The sources of net quality are as follows:(11)R=ρvρlρdαdt.The vapor volume fraction is related to the number of bubbles per volume of liquid to get the following equation:(12)α=nb4/3πRB31+nb4/3pRB3,where nb is the bubble number density.Substitute equations (11) and (12) into equation (10) to derive the following equation:(13)R=ρvρlρα1−α3RB23Pv−Pρl,RB=α1−α34π1n1/3.where R is the mass transfer rate and RB is the bubble radius.The mass transfer rate in the model is proportional toαv1−αv. In addition, the value of this function is zero when α=0 and α=1 and reaches the maximum value between them. Also, in this model, the only parameter that must be determined is the number of bubbles per unit volume of liquid. If it is assumed that no bubbles are generated or destroyed, the bubble number density will remain unchanged.The final form of the model is as follows:WhenPv>P,(14)Re=ρvρlρα1−α3RB23Pv−Pρl.WhenPv≤P,(15)Rc=ρvρlρα1−α3RB23P−Pvρl.There are three existing turbulence numerical simulation methods: direct numerical simulation (DNS), Reynolds average Navier–Stokes (RANS), and large eddy simulation (LES). DNS is currently limited to small scale and small Re-number simulation because it consumes too much computing resources and is difficult to apply in practical engineering. RANS is the application of statistical theory of turbulence, which is the simulation method commonly used in engineering. The unsteady N-S equation is averaged over time, and the Reynolds stress is assumed based on the theoretical knowledge of turbulence, experimental data, or DNS results; that is, various empirical and semiempirical constitutive relations are assumed, so that the average Reynolds equation of turbulence is closed. According to the different starting points for model processing, turbulence model theory can be divided into two categories: one type introduces the governing equation of the second-order pulsation term to form a second-order moment closed model called the Reynolds stress model; another is the eddy viscosity closed mode based on Boussinesq’s eddy viscosity hypothesis, such as zero-equation model, one-equation model, and two-equation model. Due to the large amount of calculation of the Reynolds stress model, the application occasions are very limited. The eddy viscous closed mode has lower requirements on the computer and at the same time can obtain calculation results that meet engineering requirements.The zero-equation model has a common shortcoming; that is, the turbulence viscosity coefficient only depends on the local flow parameters and has nothing to do with the flow elsewhere, which is inconsistent with experimental observations. On this basis, one-equation model and two-equation model were developed. Among the two-equation model, thek-ε model has certain accuracy and consumes limited computing resources, so it has been widely used. This paper compares the effects of three k-ε models on the experimental results.


## 2.3. Simulation Model

Use the preprocessing software Gambit to mesh the calculation space. Import the model into Gambit and use a rectangular area five times the length and three times the width of the model as the boundary of the computational domain. The area near the claws needs to be focused, so the density of the grid nodes near the claws should be increased. The height of the first cell perpendicular to the body surface is 0.05 mm. This height is chosen to makey+ of most of the cells in contact with the body surface within the effective range of the standard wall function. At the same time, considering the calculation speed, the outer division step of the model is taken as 0.1mm, the boundary and the movable chelate are set after the compliant grid, and the interior is filled with seawater. The model is shown in Figure 3.Figure 3
Grid chart in Gambit.Import the .mesh file exported by Gambit into Fluent, and set the seawater parameters according to the living environment of pistol shrimps, as shown in Table1. Use Fluent for rapid closing simulation, and set the model boundary parameters as shown in Table 2. After consulting the literature and referring to the official Fluent parameter settings, the simulation parameters are shown in Table 3. The movement of the dactylus is controlled by Fluent UDF.Table 1
Seawater parameters.
ParameterModel valueActual valueTemperature (°C)2018–24Salinity5%5%Density (kg/m3)10241024Dynamic viscosity (mPa s)1.1401.140Pressure (Pa)3028201.5∼6.5×105Saturated vapor pressure (Pa)22922291.66Table 2
Model boundary parameters.
Boundary nameBoundary typeDactylusWallPollexWallOutletPressure outlet (302820 Pa)Move1DeformingMove2DeformingSeawaterInteriorTable 3
Preset parameters for CFD simulation.
ParameterPreset valueSolverPressure basedTimeTransientTurbulence modelRNGk-ε modelPressure-velocity couplingPISOTransient formulationFirst-order implicitOther term spatial discretizationFirst-order upwindTime step size1×10−6 (s)Number of time steps500


## 3. Results and Discussion

The interaction process between the claws and seawater is short, so data postprocessing is needed to select appropriate feature quantities to characterize the interaction process and the result. What can be obtained intuitively are the pressure cloud graph, velocity cloud graph, and bubble volume fraction cloud graph during the closing process, which can reflect the cavitation phenomenon in the closing process in a certain extent. The dynamic response curve of seawater is obtained through data postprocessing, which can show the changes of pressure, velocity, and bubble size at a certain point in the process of interaction between the claws and seawater. Through the comparative analysis under different conditions, the influence of different parameters on the closing speed and cavitation effect of the claws can be concluded.


### 3.1. Simulation Process

Under the condition of closing angular velocity of 3500 rad/s, closing time of 0.5 ms, and closing angle of 90°, the pressure cloud graph, velocity cloud graph, and bubble volume fraction cloud graph are set to monitor, and images of pressure, velocity, and bubble volume change with time are obtained.It can be seen from Figure4 that, due to the rapid movement of the dactylus, a low-pressure area is formed at its tip. At the same time, as the claw closes, the seawater speed increases continuously, and the seawater between the dactylus and the pollex is compressed, and the pressure increases to form a jet. The jet velocity reaches 40 m/s, causing the seawater pressure in front to rapidly decrease below the saturated vapor pressure, and the gas core scattered in the seawater quickly grows into bubbles and moves at high speed in the direction of the jet. At 0.6 ms, the bubbles collapse instantly under the huge pressure of the surrounding seawater, forming huge pressure. The seawater basically returns to normal at 0.7 ms.Figure 4
Changes in pressure of seawater during claw closing.It can be seen from Figures5 and 6 that the claw outlet and bubble formation position are ideal positions to measure the pressure, velocity, and bubble size during the closing process, so mark points 1 and 2 in Figure 7 to analyze the specific numerical changes during the closing process and the law of change.Figure 5
Changes in velocity of seawater during claw closing.Figure 6
Changes in bubble volume during claw closing.Figure 7
Position of point 1 and point 2.


### 3.2. Influence of Turbulence Model on Simulation Accuracy

Using a closing angular velocity of 7000 rad/s, a closing time of 0.25 ms, and a closing angle of 90°, three sets of simulations were performed by changing the turbulence model. The results are shown in Figure8, which is the dynamic response curve of seawater pressure under different turbulence models. It can be seen that, under different turbulence models, the morphological characteristics of the dynamic response curve of seawater are consistent. In the closing stage, the pressure changes under the three models are basically the same, and the overall trend is increasing, with the maximum pressure reaching 2×106Pa. After the cavitation bubbles are produced, the three models appear different. The pressure curve under the Realizable model changes greatly and is unstable, which is manifested by the large difference between the peak and valley values of the pressure. The pressure peak appears earlier under the Standard model, indicating that the bubble bursts earlier.Figure 8
Pressure curves under different turbulence models. (a) Pressure of dactylus. (b) Pressure of point 1. (c) Pressure of point 2.
(a)(b)(c)Figure9 is the velocity change curve of point 1 and point 2 under different turbulence models. It can be found that the dynamic response curve of seawater has the same morphological characteristics. Within 0∼0.4 ms after closing, the velocity changes under the three models are basically consistent; after 0.4 ms, the seawater reverses in the Realizable model and the RNG model, causing the velocity of point1 to increase again, but the Standard model does not appear. The velocity of point 2 does not obviously change because the location is far from the channel.Figure 9
Velocity curves under different turbulence models. (a) Velocity curves of point 1. (b) Velocity curves of point 2.
(a)(b)Figures10(a)–10(c) show the images of cavitation bubbles under different turbulence models. The shapes of cavitation bubbles produced under different turbulence models are similar. The bubble vapor volume fraction generated by the RNG model and the Realizable model is larger, and the bubble vapor volume fraction generated by the Standard model is smaller; the bubble size generated by the RNG model and the Standard model is not much different, and the bubble generated by the Realizable model is larger.Figure 10
The proportion of cavitation bubble under different turbulence models. (a) RNG models. (b) Standard models. (c) Realizable models. (d) The proportion of cavitation bubble curves under different turbulence models.
(a)(b)(c)(d)The same result can be obtained from Figure10(d); that is, the size of the bubbles generated by the RNG model and the Standard model is not much different, and the bubbles generated by the Realizable model are larger. Referring to the simulation settings of other scholars, it is shown that the RNG model is more in line with the actual situation of this research [14, 26, 27].


### 3.3. Effect of Closing Angular Velocity

In order to obtain the relationship between the closing process and the closing speed, three sets of angular velocity-time are designed as shown in Figure11. Keep the closing angle unchanged at 90°, the peak closing angular velocities are 1810 rad/s, 3631 rad/s, and 7300 rad/s, and the closing time is 1 ms, 0.5 ms, and 0.25 ms, respectively. Since the rotation direction of the claws is clockwise, the closing angular velocity is negative.Figure 11
Three types of closing angular velocity of the pistol shrimp’s claw.Import the measurement data under different closing speeds into MATLAB to observe the influence of closing angular speed on pressure more intuitively. As shown in Figures12(a)–12(c), it is the pressure change curve of the dactylus, point 1, and point 2 under different closing angular velocities. It can be found that the dynamic response curve morphology characteristics of seawater under different closing angular velocities are the same but with different peaks. Since the closing angle has a fixed value of 90°, a large closing speed will result in a short closing time, which represents an earlier response; a small closing speed will result in a long closing time, which represents a lag in response. The small pressure is positively correlated with the speed.Figure 12
Pressure curves under different closing angular velocity. (a) Pressure of dactylus. (b) Pressure of point 1. (c) Pressure of point 2. (d) Fitting curve between the minimum of pressure of dactylus and closing angular velocity.
(a)(b)(c)(d)Take the minimum pressure of the dactylus at the jet and the closing angular velocity to perform MATLAB fitting to obtain the relationship between cavitation and the closing angular velocity. The result is shown in Figure12(d).Fitting with exponential function, the fitting results are as follows:(16)P=axb,a=4.147×1026,b=−6.548.The fitting result is good, and the saturated vapor pressure under the model is substituted into the function to obtain the critical closing angular velocity asωz=2466.6rad/s.Figures13(a)and 13(b) show the velocity curves of point 1 and point 2 with time under different closing speeds. It can be found that the velocity curve has the same change trend with the pressure curve. The curve shape characteristics are the same, and the peak values are different. The velocity at point1 and point 2 is positively correlated with the closing velocity.Figure 13
Velocity curves under different closing angular velocity. (a) Velocity of point 1. (b) Velocity of point 2. (c) Fitting curve between the velocity peak of point 1 and closing angular velocity.
(a)(b)(c)In order to summarize the relationship between closure velocity and cavitation, the peak velocity and angular velocity peak at point1 are imported into MATLAB for fitting. The fitting expression is(17)V=p1x2+p2x+p3,p1=3.713×10−6,p2=−0.009242,p3=24.61.The drawn image is shown in Figure13(c), and the fitting result is good. It is known from equation (3) that when the jet velocity V0=24.22m/s, cavitation occurs; when V0=34.26m/s, stable cavitation occurs. Substituting V0=24.22m/s into equation (17), the critical closing angular velocity for cavitation is 2446.8 rad/s, which is basically the same as the result obtained from equation (19).Figure14 is the volume image of cavitation bubbles formed at different closing angular velocities, and Figure 15(a) is the bubble volume fraction curve at different closing angular velocities. The image shows that when the closing angular velocity is small, the velocity of the water jet is small, and the pressure drop cannot meet the requirements of cavitation, so cavitation bubbles cannot be generated. When the angular velocity is large, the velocity of the water jet increases, causing the pressure to drop below the saturated vapor pressure, so the generated cavitation bubbles are larger. The proportion of bubbles is positively correlated with the closing angular velocity, but as the closing speed increases, the growth rate of the proportion of bubbles slows down.Figure 14
Volume of cavitation bubbles formed at different closing speeds. (a) 1810 rad/s. (b) 3631 rad/s. (c) 7300 rad/s.
(a)(b)(c)Figure 15
Curves of bubble size and bubble movement distance under different closing angular velocities. (a) Bubble volume fraction curve. (b) Fitting curve between bubble size and closing angular velocity. (c) Fitting curve between bubble distance and closing angular velocity.
(a)(b)(c)From equations (16) and (17), the critical closing angular velocity is 2456rad/s, and the bubble size and movement distance can be set to zero at this time. The bubble is ellipsoidal at 3631rad/s, and, at 7300rad/s due to the closing angular velocity being too high, the bubble deforms into an irregular shape. Observe that the sizes of the bubble at 3631rad/s and 7300rad/s in Figure 13 are 4mm and 12mm, respectively, and the bubble movement distance is 6mm and 10mm, respectively. The bubble size and closing angular velocity image can be obtained as shown in Figure 15(b), and the fitting result is(18)d=p1x2+p2x+p3,p1=−2.526×10−7,p2=0.004942,p3=−10.61.The image of bubble movement distance and closing angular velocity is shown in Figure15(c) and the fitting result is(19)d2=axb+c,a=−3.168×107,b=−1.9,c=11.44.


### 3.4. Effect of Closing Angular Acceleration

When studying the effect of the closing angular acceleration on the interaction between the claws and seawater, the maximum angular velocity is kept unchanged atωmax=4000rad/s, the closing angle remains unchanged at 90°, and the angular acceleration is α1=0.75×107rad/s2, α2=1×107rad/s2, and α3=1.5×107rad/s2; the angular velocity-time image is shown in Figure 16. Since the rotation direction of the claws is clockwise, the closing angular acceleration is negative.Figure 16
Three types of closing angular acceleration curve.Import the measurement data at different closing speeds into MATLAB to observe the effect of closing angular acceleration on pressure more intuitively. Figure17 is the pressure change curve of the dactylus, point 1, and point 2 under different closing angular accelerations. It can be found that the dynamic response curve of seawater has the same shape characteristics and different peak values under different closing angular accelerations. Since the closing angle is a fixed value of 90°, the closing time under α2=1×107rad/s2 is exactly 0.5ms, and the closing angular acceleration under α3=1.5×107rad/s2 is relatively larger, so the closing time is short, and the response is earlier; when the closing angular acceleration is smaller under α1=0.75×107rad/s2, the closing time is longer, and the response lags. The peaks of pressure are not much different, and they increase slightly with the increase of angular acceleration; the valleys of pressure are exactly the same.Figure 17
Pressure curves under different closing angular acceleration. (a) Pressure of dactylus. (b) Pressure of point 1. (c) Pressure of point 2.
(a)(b)(c)After importing the measurement data into MATLAB, as shown in Figure18, it is the jet velocity of point1 and point 2 with time under different closed angular accelerations. The velocity of point 1 increases again within 0.7∼0.8ms because the velocity of seawater fills the space instantly when the cavitation bubbles collapse. The position difference between point 1 and point 2 is Δd12=6.0mm, indicating that the velocity of the jet and cavitation bubbles moving from point 1 to point 2 is attenuated from about 55m/s to about 20m/s.Figure 18
Velocity curves under different closing angular acceleration. (a) Velocity of point 1. (b) Velocity of point 2.
(a)(b)As shown in Figure19, there are three types of angular acceleration maximum bubble volume cloud images, from top to bottom, α1=0.75×107rad/s2, α2=1×107rad/s2, and α3=1.5×107rad/s2; the time to reach the maximum bubble is 0.67ms, 0.63ms, and 0.59ms, respectively. It can be seen from the image that the cavitation bubble has an ellipsoidal shape, and the farthest movement distance is about d=6.0mm. The section of the cavitation bubble is elliptical, with a major axis of about 5mm and a minor axis of about 3mm. The volume of bubbles formed in the case of α2 is the smallest, and the volume of bubbles formed in the cases of α1 and α3 is approximately the same.Figure 19
Maximum bubble volume cloud image under different closing angular accelerations.Figure20 shows the bubble volume fraction curve under different closing angular accelerations. The bubble volume fraction indicates the bubble volume size. It can be seen that the bubble volume fraction of α1 reaches the maximum at 0.67ms, which is 92%; the bubble volume fraction of α2 reaches the maximum at 0.63ms, which is 78%; the bubble volume fraction of α3 reaches the maximum at 0.59ms, which is 98%. It can be calculated that it takes 0.09ms for the bubble to move to the farthest position in the case of α1; it takes 0.1ms for the bubble to move to the farthest position in the case of α2; it takes 0.09ms for the bubble to move to the farthest position in the case of α3.Figure 20
Bubble volume fraction curve under different closing angular accelerations.The results of the study show the following: (1) The loud popping sound of the pistol shrimp is caused by the collapse of the cavitation bubbles in the high-speed water jet produced by the rapid closing of the pistol shrimp claws. (2) The pressure, velocity, bubble volume fraction, bubble movement distance, and bubble size increase with the closing angular velocity, and there is a minimum closing speed for the cavitation effect. (3) The closing angular acceleration has a very small effect on pressure, velocity, bubble volume fraction, bubble movement distance, and bubble size.


## 3.1. Simulation Process

Under the condition of closing angular velocity of 3500 rad/s, closing time of 0.5 ms, and closing angle of 90°, the pressure cloud graph, velocity cloud graph, and bubble volume fraction cloud graph are set to monitor, and images of pressure, velocity, and bubble volume change with time are obtained.It can be seen from Figure4 that, due to the rapid movement of the dactylus, a low-pressure area is formed at its tip. At the same time, as the claw closes, the seawater speed increases continuously, and the seawater between the dactylus and the pollex is compressed, and the pressure increases to form a jet. The jet velocity reaches 40 m/s, causing the seawater pressure in front to rapidly decrease below the saturated vapor pressure, and the gas core scattered in the seawater quickly grows into bubbles and moves at high speed in the direction of the jet. At 0.6 ms, the bubbles collapse instantly under the huge pressure of the surrounding seawater, forming huge pressure. The seawater basically returns to normal at 0.7 ms.Figure 4
Changes in pressure of seawater during claw closing.It can be seen from Figures5 and 6 that the claw outlet and bubble formation position are ideal positions to measure the pressure, velocity, and bubble size during the closing process, so mark points 1 and 2 in Figure 7 to analyze the specific numerical changes during the closing process and the law of change.Figure 5
Changes in velocity of seawater during claw closing.Figure 6
Changes in bubble volume during claw closing.Figure 7
Position of point 1 and point 2.


## 3.2. Influence of Turbulence Model on Simulation Accuracy

Using a closing angular velocity of 7000 rad/s, a closing time of 0.25 ms, and a closing angle of 90°, three sets of simulations were performed by changing the turbulence model. The results are shown in Figure8, which is the dynamic response curve of seawater pressure under different turbulence models. It can be seen that, under different turbulence models, the morphological characteristics of the dynamic response curve of seawater are consistent. In the closing stage, the pressure changes under the three models are basically the same, and the overall trend is increasing, with the maximum pressure reaching 2×106Pa. After the cavitation bubbles are produced, the three models appear different. The pressure curve under the Realizable model changes greatly and is unstable, which is manifested by the large difference between the peak and valley values of the pressure. The pressure peak appears earlier under the Standard model, indicating that the bubble bursts earlier.Figure 8
Pressure curves under different turbulence models. (a) Pressure of dactylus. (b) Pressure of point 1. (c) Pressure of point 2.
(a)(b)(c)Figure9 is the velocity change curve of point 1 and point 2 under different turbulence models. It can be found that the dynamic response curve of seawater has the same morphological characteristics. Within 0∼0.4 ms after closing, the velocity changes under the three models are basically consistent; after 0.4 ms, the seawater reverses in the Realizable model and the RNG model, causing the velocity of point1 to increase again, but the Standard model does not appear. The velocity of point 2 does not obviously change because the location is far from the channel.Figure 9
Velocity curves under different turbulence models. (a) Velocity curves of point 1. (b) Velocity curves of point 2.
(a)(b)Figures10(a)–10(c) show the images of cavitation bubbles under different turbulence models. The shapes of cavitation bubbles produced under different turbulence models are similar. The bubble vapor volume fraction generated by the RNG model and the Realizable model is larger, and the bubble vapor volume fraction generated by the Standard model is smaller; the bubble size generated by the RNG model and the Standard model is not much different, and the bubble generated by the Realizable model is larger.Figure 10
The proportion of cavitation bubble under different turbulence models. (a) RNG models. (b) Standard models. (c) Realizable models. (d) The proportion of cavitation bubble curves under different turbulence models.
(a)(b)(c)(d)The same result can be obtained from Figure10(d); that is, the size of the bubbles generated by the RNG model and the Standard model is not much different, and the bubbles generated by the Realizable model are larger. Referring to the simulation settings of other scholars, it is shown that the RNG model is more in line with the actual situation of this research [14, 26, 27].


## 3.3. Effect of Closing Angular Velocity

In order to obtain the relationship between the closing process and the closing speed, three sets of angular velocity-time are designed as shown in Figure11. Keep the closing angle unchanged at 90°, the peak closing angular velocities are 1810 rad/s, 3631 rad/s, and 7300 rad/s, and the closing time is 1 ms, 0.5 ms, and 0.25 ms, respectively. Since the rotation direction of the claws is clockwise, the closing angular velocity is negative.Figure 11
Three types of closing angular velocity of the pistol shrimp’s claw.Import the measurement data under different closing speeds into MATLAB to observe the influence of closing angular speed on pressure more intuitively. As shown in Figures12(a)–12(c), it is the pressure change curve of the dactylus, point 1, and point 2 under different closing angular velocities. It can be found that the dynamic response curve morphology characteristics of seawater under different closing angular velocities are the same but with different peaks. Since the closing angle has a fixed value of 90°, a large closing speed will result in a short closing time, which represents an earlier response; a small closing speed will result in a long closing time, which represents a lag in response. The small pressure is positively correlated with the speed.Figure 12
Pressure curves under different closing angular velocity. (a) Pressure of dactylus. (b) Pressure of point 1. (c) Pressure of point 2. (d) Fitting curve between the minimum of pressure of dactylus and closing angular velocity.
(a)(b)(c)(d)Take the minimum pressure of the dactylus at the jet and the closing angular velocity to perform MATLAB fitting to obtain the relationship between cavitation and the closing angular velocity. The result is shown in Figure12(d).Fitting with exponential function, the fitting results are as follows:(16)P=axb,a=4.147×1026,b=−6.548.The fitting result is good, and the saturated vapor pressure under the model is substituted into the function to obtain the critical closing angular velocity asωz=2466.6rad/s.Figures13(a)and 13(b) show the velocity curves of point 1 and point 2 with time under different closing speeds. It can be found that the velocity curve has the same change trend with the pressure curve. The curve shape characteristics are the same, and the peak values are different. The velocity at point1 and point 2 is positively correlated with the closing velocity.Figure 13
Velocity curves under different closing angular velocity. (a) Velocity of point 1. (b) Velocity of point 2. (c) Fitting curve between the velocity peak of point 1 and closing angular velocity.
(a)(b)(c)In order to summarize the relationship between closure velocity and cavitation, the peak velocity and angular velocity peak at point1 are imported into MATLAB for fitting. The fitting expression is(17)V=p1x2+p2x+p3,p1=3.713×10−6,p2=−0.009242,p3=24.61.The drawn image is shown in Figure13(c), and the fitting result is good. It is known from equation (3) that when the jet velocity V0=24.22m/s, cavitation occurs; when V0=34.26m/s, stable cavitation occurs. Substituting V0=24.22m/s into equation (17), the critical closing angular velocity for cavitation is 2446.8 rad/s, which is basically the same as the result obtained from equation (19).Figure14 is the volume image of cavitation bubbles formed at different closing angular velocities, and Figure 15(a) is the bubble volume fraction curve at different closing angular velocities. The image shows that when the closing angular velocity is small, the velocity of the water jet is small, and the pressure drop cannot meet the requirements of cavitation, so cavitation bubbles cannot be generated. When the angular velocity is large, the velocity of the water jet increases, causing the pressure to drop below the saturated vapor pressure, so the generated cavitation bubbles are larger. The proportion of bubbles is positively correlated with the closing angular velocity, but as the closing speed increases, the growth rate of the proportion of bubbles slows down.Figure 14
Volume of cavitation bubbles formed at different closing speeds. (a) 1810 rad/s. (b) 3631 rad/s. (c) 7300 rad/s.
(a)(b)(c)Figure 15
Curves of bubble size and bubble movement distance under different closing angular velocities. (a) Bubble volume fraction curve. (b) Fitting curve between bubble size and closing angular velocity. (c) Fitting curve between bubble distance and closing angular velocity.
(a)(b)(c)From equations (16) and (17), the critical closing angular velocity is 2456rad/s, and the bubble size and movement distance can be set to zero at this time. The bubble is ellipsoidal at 3631rad/s, and, at 7300rad/s due to the closing angular velocity being too high, the bubble deforms into an irregular shape. Observe that the sizes of the bubble at 3631rad/s and 7300rad/s in Figure 13 are 4mm and 12mm, respectively, and the bubble movement distance is 6mm and 10mm, respectively. The bubble size and closing angular velocity image can be obtained as shown in Figure 15(b), and the fitting result is(18)d=p1x2+p2x+p3,p1=−2.526×10−7,p2=0.004942,p3=−10.61.The image of bubble movement distance and closing angular velocity is shown in Figure15(c) and the fitting result is(19)d2=axb+c,a=−3.168×107,b=−1.9,c=11.44.


## 3.4. Effect of Closing Angular Acceleration

When studying the effect of the closing angular acceleration on the interaction between the claws and seawater, the maximum angular velocity is kept unchanged atωmax=4000rad/s, the closing angle remains unchanged at 90°, and the angular acceleration is α1=0.75×107rad/s2, α2=1×107rad/s2, and α3=1.5×107rad/s2; the angular velocity-time image is shown in Figure 16. Since the rotation direction of the claws is clockwise, the closing angular acceleration is negative.Figure 16
Three types of closing angular acceleration curve.Import the measurement data at different closing speeds into MATLAB to observe the effect of closing angular acceleration on pressure more intuitively. Figure17 is the pressure change curve of the dactylus, point 1, and point 2 under different closing angular accelerations. It can be found that the dynamic response curve of seawater has the same shape characteristics and different peak values under different closing angular accelerations. Since the closing angle is a fixed value of 90°, the closing time under α2=1×107rad/s2 is exactly 0.5ms, and the closing angular acceleration under α3=1.5×107rad/s2 is relatively larger, so the closing time is short, and the response is earlier; when the closing angular acceleration is smaller under α1=0.75×107rad/s2, the closing time is longer, and the response lags. The peaks of pressure are not much different, and they increase slightly with the increase of angular acceleration; the valleys of pressure are exactly the same.Figure 17
Pressure curves under different closing angular acceleration. (a) Pressure of dactylus. (b) Pressure of point 1. (c) Pressure of point 2.
(a)(b)(c)After importing the measurement data into MATLAB, as shown in Figure18, it is the jet velocity of point1 and point 2 with time under different closed angular accelerations. The velocity of point 1 increases again within 0.7∼0.8ms because the velocity of seawater fills the space instantly when the cavitation bubbles collapse. The position difference between point 1 and point 2 is Δd12=6.0mm, indicating that the velocity of the jet and cavitation bubbles moving from point 1 to point 2 is attenuated from about 55m/s to about 20m/s.Figure 18
Velocity curves under different closing angular acceleration. (a) Velocity of point 1. (b) Velocity of point 2.
(a)(b)As shown in Figure19, there are three types of angular acceleration maximum bubble volume cloud images, from top to bottom, α1=0.75×107rad/s2, α2=1×107rad/s2, and α3=1.5×107rad/s2; the time to reach the maximum bubble is 0.67ms, 0.63ms, and 0.59ms, respectively. It can be seen from the image that the cavitation bubble has an ellipsoidal shape, and the farthest movement distance is about d=6.0mm. The section of the cavitation bubble is elliptical, with a major axis of about 5mm and a minor axis of about 3mm. The volume of bubbles formed in the case of α2 is the smallest, and the volume of bubbles formed in the cases of α1 and α3 is approximately the same.Figure 19
Maximum bubble volume cloud image under different closing angular accelerations.Figure20 shows the bubble volume fraction curve under different closing angular accelerations. The bubble volume fraction indicates the bubble volume size. It can be seen that the bubble volume fraction of α1 reaches the maximum at 0.67ms, which is 92%; the bubble volume fraction of α2 reaches the maximum at 0.63ms, which is 78%; the bubble volume fraction of α3 reaches the maximum at 0.59ms, which is 98%. It can be calculated that it takes 0.09ms for the bubble to move to the farthest position in the case of α1; it takes 0.1ms for the bubble to move to the farthest position in the case of α2; it takes 0.09ms for the bubble to move to the farthest position in the case of α3.Figure 20
Bubble volume fraction curve under different closing angular accelerations.The results of the study show the following: (1) The loud popping sound of the pistol shrimp is caused by the collapse of the cavitation bubbles in the high-speed water jet produced by the rapid closing of the pistol shrimp claws. (2) The pressure, velocity, bubble volume fraction, bubble movement distance, and bubble size increase with the closing angular velocity, and there is a minimum closing speed for the cavitation effect. (3) The closing angular acceleration has a very small effect on pressure, velocity, bubble volume fraction, bubble movement distance, and bubble size.


## 4. Conclusions

In this paper, firstly, a theoretical modeling of the interaction between the claws and seawater is carried out, and a cavitation model suitable for this subject is derived. Then, based on the theory of fluid-structure coupling and the physical properties of seawater, a simulation model of pistol shrimp claws and seawater was established, and the parameters were matched by combining the parameter matching reference values given by the Fluent database. Finally, taking a certain working condition as an example, the interaction process between the claws and seawater is analyzed, and the pressure, velocity, and bubble volume fraction in the closing process are obtained through the data postprocessing. The influence of different turbulence models in the simulation parameters on the simulation accuracy is analyzed, and the model that best fits the actual situation of the research content of this subject is determined. Then, using the controlled variable method, keeping the other variables the same, the effects of closing angular velocity and angular acceleration on the pressure, jet velocity, and bubble volume fraction in the closing process are studied, and the effects of closing angular velocity and angular acceleration on the pressure, jet velocity, and bubble volume fraction in the closing process are summarized. At the same time, the relationship between bubble movement distance and bubble size, closing angular velocity, and angular acceleration is studied. The mechanism of the closed jet of the claw claws is revealed, which provides a theoretical basis for the development of new underwater kinetic energy weapons.



---

*Source: 9975952-2021-05-15.xml*