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

*e 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. *is 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. *is article studied the influence of different fluid models on the simulation results. *e 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.


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. e 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][2][3][4][5][6][7][8][9][10][11][12]. is article carried out a study on the rapid closing jet mechanism of pistol shrimp's claws. e 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][14][15][16]; on the other hand, it focuses on the shooting of the pistol through the high-speed camera system. e 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]. e 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). e 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. is 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. e 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. e 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.

Physical Model.
e body of a pistol shrimp is 6 ∼ 10 cm long, with a large claw and a small claw. e large claw is composed of dactylus and pollex. e specific structure is shown in Figure 1. is 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. e 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).

e N-S equation of incompressible viscous fluid is
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. ese cavitation bubbles are carried by the water flow and collapse when they reach the high-pressure zone. e whole process of the initiation, growth, and collapse of the cavitation bubbles is called cavitation. e degree of cavitation is usually expressed by the dimensionless cavitation index σ, which is defined as where P 0 is the environmental pressure, V 0 is the exit speed, P v is the saturated vapor pressure of water, and ρ is the density of water. e 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 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 where P g is the partial pressure of gas in bubble.
If P v 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 2 Mathematical Problems in Engineering Mathematical Problems in Engineering regarded as an ideal gas adiabatic process, so P g in equation (5) becomes where P g0 , V 0 , and R 0 , respectively, are partial pressure of gas, bubble radius, and bubble volume in an initial state and P g , V, and R, 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: e 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. is 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.
e Rayleigh-Plesset (R-P) equation is where P B 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. e 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]. e actual fluid also contains noncondensable gas (NCG). e saturated mass concentration of NCG that can be dissolved by pure water at a standard atmospheric pressure is 1.5 × 10 3 mg/L. is 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: where α is the gas volume fraction, ρ v is the gas density, V v is the gas velocity, and R e and R c are the mass transfer source terms, which, respectively, are related to the growth and collapse of the vapor bubble. In formula (9), R e and R c represent the mass transfer between liquid and gas phases in cavitation. ey are based on the Rayleigh-Plesset equation of formula (8), which describes the growth of a single vapor bubble in the liquid. e equation for the vapor volume fraction from equation (9) has a general form: e sources of net quality are as follows: e vapor volume fraction is related to the number of bubbles per volume of liquid to get the following equation: where n b is the bubble number density. Substitute equations (11) and (12) into equation (10) to derive the following equation: where R is the mass transfer rate and R B is the bubble radius. e mass transfer rate in the model is proportional to 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. e final form of the model is as follows: When P v ≤ P, ere 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 Renumber 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. e 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. e eddy viscous closed mode has lower requirements on the computer and at the same time can obtain calculation results that meet engineering requirements.
e 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, the k-ε model has certain accuracy and consumes limited computing resources, so it has been widely used.
is paper compares the effects of three k-ε models on the experimental results.

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. e area near the claws needs to be focused, so the density of the grid nodes near the claws should be increased. e height of the first cell perpendicular to the body surface is 0.05 mm. is height is chosen to make y+ 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.1 mm, the boundary and the movable chelate are set after the compliant grid, and the interior is filled with seawater. e model is shown in Figure 3. 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 Table 1. 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. e movement of the dactylus is controlled by Fluent UDF.

Results and Discussion
e 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. e 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. rough the comparative analysis under different conditions, the influence of different parameters on the closing speed and cavitation effect of the claws can be concluded.

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 Figure 4 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. e 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. e seawater basically returns to normal at 0.7 ms.
It can be seen from Figures 5 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.

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. e results are shown in Figure 8, 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 × 10 6 Pa. After the cavitation bubbles are produced, the three models appear different. e 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.
e pressure peak appears earlier under the Standard model, indicating that the bubble bursts earlier.

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 Figure 11. 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. 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 Figures 12(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. e small pressure is positively correlated with the speed.
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. e result is shown in Figure 12(d).     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. e fitting expression is e drawn image is shown in Figure 13(c), and the fitting result is good. It is known from equation (3) (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). Figure 14 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. e 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. e 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.     Mathematical Problems in Engineering From equations (16) and (17), the critical closing angular velocity is 2456 rad/s, and the bubble size and movement distance can be set to zero at this time. e bubble is ellipsoidal at 3631 rad/s, and, at 7300 rad/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 3631 rad/s and 7300 rad/s in Figure 13 are 4 mm and 12 mm, respectively, and the bubble movement distance is 6 mm and 10 mm, respectively. e bubble size and closing angular velocity image can be obtained as shown in Figure 15(b), and the fitting result is    Mathematical Problems in Engineering

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 � 4000 rad/s, the closing angle remains unchanged at 90°, and the angular acceleration is α 1 � 0.75 × 10 7 rad/s 2 , α 2 � 1 × 10 7 rad/s 2 , and α 3 � 1.5 × 10 7 rad/s 2 ; 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. Import the measurement data at different closing speeds into MATLAB to observe the effect of closing angular acceleration on pressure more intuitively. Figure 17 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 × 10 7 rad/s 2 is exactly 0.5 ms, and the closing angular acceleration under α 3 � 1.5 × 10 7 rad/s 2 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 × 10 7 rad/s 2 , the closing time is longer, and the response lags. e 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.     After importing the measurement data into MATLAB, as shown in Figure 18, it is the jet velocity of point1 and point 2 with time under different closed angular accelerations. e velocity of point 1 increases again within 0.7 ∼ 0.8 ms because the velocity of seawater fills the space instantly when the cavitation bubbles collapse. e position difference between point 1 and point 2 is Δd 12 � 6.0 mm, indicating that the velocity of the jet and cavitation bubbles moving from point 1 to point 2 is attenuated from about 55 m/s to about 20 m/s. As shown in Figure 19, there are three types of angular acceleration maximum bubble volume cloud images, from top to bottom, α 1 � 0.75 × 10 7 rad/s 2 , α 2 � 1 × 10 7 rad/s 2 , and α 3 � 1.5 × 10 7 rad/s 2 ; the time to reach the maximum bubble is 0.67 ms, 0.63 ms, and 0.59 ms, 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.0 mm. e section of the cavitation bubble is elliptical, with a major axis of about 5 mm and a minor axis of about 3 mm. e 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 20 shows the bubble volume fraction curve under different closing angular accelerations. e bubble volume fraction indicates the bubble volume size. It can be seen that the bubble volume fraction of α 1 reaches the maximum at 0.67 ms, which is 92%; the bubble volume fraction of α 2 reaches the maximum at 0.63 ms, which is 78%; the bubble volume fraction of α 3 reaches the maximum at 0.59 ms, which is 98%. It can be calculated that it takes 0.09 ms for the bubble to move to the farthest position in the case of α 1 ; it takes 0.1 ms for the bubble to move to the farthest position in the case of α 2 ; it takes 0.09 ms for the bubble to move to the farthest position in the case of α 3 .
e results of the study show the following: (1) e 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) e 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) e closing angular acceleration has a very small effect on pressure, velocity, bubble volume fraction, bubble movement distance, and bubble size.

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. en, 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. e 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. en, 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. e 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.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.