A strong coupling prediction-correction immersed boundary method

In the traditional immersed boundary methods for solving compressible fluid-structure interaction problems, conservation is one of the problems that must be considered. When the coupling boundary moves on the fixed grid, the structure coverage will change, resulting in many dead elements and emerging elements on the fluid grid. In the ghost-cell immersed boundary method, the reconstructed grid can not maintain the strict mass conservation when the dead elements and emerging elements appear. In order to overcome the shortcomings of traditional methods, a strong coupling prediction-correction immersed boundary method for compressible fluid-structure interaction problems was proposed. Firstly, the matrix representation of a general fluid-structure coupling system was described, and a strong coupling Gauss-Seidel iterative scheme of fluid-structure coupling system was derived. Furthermore, a prediction-correction scheme was derived, and a prediction-correction immersed boundary method was proposed. The fluid-structure coupling boundary was regarded as a free surface, and the space originally occupied by the solid was initialized as zero mass elements, allowing the fluid to pass through the coupling boundary freely. For the calculation of fluid, the second-order MUSCL finite volume scheme with the minmod limiter and the AUSM+-up flux based on Zha-Bilgen splitting were used to advance the time step with the third-order Runge-Kutta scheme. In the correction step, the transport process was realized by a set of mass conservation transport rules. The transport algorithm could be summarized as marking the fluid inside the boundary, dividing and moving the fluid in a uniform way according to the marking order, generating a flow pointing to the outside of the boundary, and finally applying a velocity correction near the boundary to ensure the no-slip condition. The marking and transport algorithm avoided the tedious geometric treatment of the cut-cells, and ensured the easy implementation of the algorithm. For the calculation of solids, the first-order difference scheme and the implicit dynamic finite element scheme were used to solve the rigid body and linear elastic body respectively, and the Gauss quadrature was used to obtain the coupling force on the solid surface. The one-dimensional and two-dimensional problems were calculated by the prediction-correction immersed boundary method. In the one-dimensional piston problem, the accuracy, conservation and convergence of the method were investigated by comparing the results with those in the literature. In the two-dimensional shock wave impact problem, the experimental optical schlieren images were compared with those obtained by the numerical simulation, and the deflection history of the plate structure was investigated. The study shows that this method can accurately maintain the mass conservation of the flow field and has the advantage of easy implementation, which is different from the traditional ghost-cell method and the cut-cell method. This method has the first-order convergence accuracy, and can accurately predict the flow field after shock diffraction and the deflection of plate under shock waves. It provides a new idea for the development of fluid-structure coupling algorithms. © 2021, Editorial Staff of EXPLOSION AND SHOCK WAVES. All right reserved.