A coupled FD-SPH method for shock-structure interaction and dynamic fracture propagation modeling

Shock wave propagation and their damage to solid structures, which involve complex multiphase and multiphysics phenomena, are difficult problems to address. In this paper, a three-dimensional graphics processing unit (GPU)-accelerated finite difference-smoothed particle hydrodynamics (FD-SPH) method was developed for the prediction of strongly compressible fluid flow and strong fluid-structure interactions. The conventional mesh-based finite difference method was used to simulate shock wave propagation, while the smoothed particle hydrodynamics method was employed to predict the dynamic behavior of solid materials. In addition, the immersed boundary method (IBM) was implemented in the GPU-accelerated FD-SPH solver for the boundary treatment of the interface between solid and fluid materials. Four numerical cases, namely shock-cylinder obstacle interaction, elastic panel deformation induced by shock wave propagation, the response of stainless steel tubes under internal blast loading, and the damage of blast loading on reinforced concrete slab, were conducted for verification of the GPU-accelerated FD-SPH solver. The numerical results were compared against the available experimental data, which shows that the GPU-accelerated FD-SPH solver is capable of capturing the shock wave propagation and its damage to structures very well.