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.
In this paper, an opposing jet is applied to the space shuttle arc leading edge and the lifting body cone leading edge in the hypersonic condition of the rarefied flow field. The DSMC numerical algorithm is used to simulate and analyze the underlying working physics of the opposing jet. The results provide a reference for designing hypersonic vehicles in near space that reduce drag and protect against heat.
Fish-like robots have been widely used in intelligent surveillance and investigation because of their high swimming efficiency and low traveling noise. Numerical simulations are usually selected to simulate the movement modes and hydrodynamic characteristics of fish-like robots during design and manufacture. However, the body-fitted grid method traditionally utilized in numerical simulations often has difficulty dealing with moving solid boundaries. In this work, the immersed boundary method, superior in handling the moving boundary conditions, is employed to simulate the movement of a fish-like robot swimming in high Reynolds number flows in combination with the RANS turbulence model. The numerical method is first validated using a fluid flowing over a square block, and the corresponding results are in good agreement with the ones reported in reference. Then, the swing of the fish-like robot under three different Reynolds numbers is studied. The lift coefficient and the drag coefficient of the fish-like robot decrease with increasing the Reynolds number. This paper provides remarkable support for future designs and applications of fish-like robots.
An adaptive particle refinement (APR) algorithm has been developed for the smoothed particle hydrodynamics (SPH) method to augment the resolution of the region of interest to achieve high accuracy and simultaneously reduce the cost of computational resources. It is widely applied in the field of fluid-controlling problems involving large interface deformations, such as the two-phase flow and fluid–structure interaction because this algorithm can capture the interface with high accuracy. Nonetheless, existing APR algorithms widely encounter computational dispersion issues at the interface of regions of different particle resolutions. Moreover, traditional shifting algorithms applied in the APR processes also have difficulties in dealing with particles with different smooth lengths. In this work, an algorithm for fast particle generation was first developed based on the accelerated ray method, which accelerates the discretization of the flow field into particles. Then, a dynamic refinement/coarsening algorithm based on the APR algorithm is proposed to solve the computational dispersion problem that occurs at the refinement/coarsening interfaces. In addition, the shifting algorithm was improved in this work to ensure the particles are always well distributed during numerical calculations and, thus, can efficiently facilitate the adaptive particle refinement/coarsening processes. Comparative analysis indicates that the robust algorithms developed for the SPH method in this work can lead to more precise and reasonable flow fields compared with the conventional SPH adaptive methods.
In natural environments, fish almost always swim in groups. Investigating the coupled mechanism of biomimetic fish exhibiting autonomous swimming capabilities advances our understanding of fish schooling phenomena and simultaneously aids in refining the structural and formation configurations of underwater robotic vehicles. This work innovatively develops an algorithm based on the Direct-Forcing Immersed Boundary Method (DF-IBM) and implements it in an efficient, modular software program written in C++. The program accelerates the calculation process by using a multigrid method. Validation against a benchmark case of flow around a cylinder, with comparison to data from the existing literature, verifies the program’s precision with discrepancies of less than 3.6%. Based on this algorithm, the paper analyzes the incompressible viscous flow during the movement of parallel-aligned biomimetic fish. It uncovers the interaction between the fish’s motion and the surrounding flow field and also reveals the hydrodynamic mechanisms of the group motion of the parallel-aligned biomimetic fish. The flow field under varying spacing and phases between the parallel-aligned biomimetic fish proves that the interaction between the flow fields induced by the two fish bodies becomes increasingly significant when decreasing the lateral spacing from to . Notably, an initial lateral convergence of the fish bodies is observed, followed by a sideways swimming pattern at a particular pitch angle, accompanied by a decrement in their forward swimming velocity as they approach each other. Additionally, this study compares flow field alterations in parallel-aligned biomimetic fish with identical lateral spacing but opposing flapping phases. The findings indicate that, irrespective of the phase, the fish exhibit an initial convergence followed by a sideways motion at a specific pitch angle. However, due to disparities in the tail’s flow field, a larger pitch angle is generated when the fish swim in unison. All the findings above will provide a solid theoretical foundation for the design and optimization of underwater robotic vehicles.
During the water entry process of a trans-domain morphing aircraft, significant impact forces are generated when the aircraft hits the water surface, which will potentially cause the deformation of the cabin structure and might damage the structure or onboard devices. Thus, it is necessary to investigate the water entry process of the cabin structure. This paper analyses changes in fluid loads and the corresponding structural responses during the water entry process. Firstly, the numerical model is established for the water entry process and the modeling method is validated by comparing the results to the experimental data. An empirical formula is developed to correlate the impact loads with the water entry velocities. Then, fluid–structure interaction analysis of the water entry process is performed using a two-way coupling approach. The relationship between structural deformation and the water entry process is then investigated. The results are compared with those without considering the structural deformation. The empirical formula is then modified to reflect the effects of the deformation. The results show that structural deformation will disperse the impact load, which represents different responses compared to the rigid cabin structure.