Numerical simulation of a complex moving rigid body under the impingement of a shock wave in 3D

Advances in Aerodynamics

Yan Jiang¹  Shihao Liu¹  Mengping Zhang¹  Ziqiang Cheng²  Jianfang Lu³  Shuhai Zhang⁴ 

[1] School of Mathematical Sciences, University of Science and Technology of China, 230026, Hefei, Anhui, China;School of Mathematics, Hefei University of Technology, 230000, Hefei, Anhui, China;School of Mathematics, South China University of Technology, 510641, Guangzhou, Guangdong, China;State Key Laboratory of Aerodynamics, China Aerodynamics Research and Development Center, 621000, Mianyang, Sichuan, China;

关键词: Inverse Lax-Wendroff procedure;    High order accuracy;    Complex moving boundary treatment;    Compressible inviscid fluid;    Compressible viscous fluid;    3D;   

DOI  :  10.1186/s42774-021-00096-5

来源: Springer

【 摘 要 】

In this paper, we take a numerical simulation of a complex moving rigid body under the impingement of a shock wave in three-dimensional space. Both compressible inviscid fluid and viscous fluid are considered with suitable boundary conditions. We develop a high order numerical boundary treatment for the complex moving geometries based on finite difference methods on fixed Cartesian meshes. The method is an extension of the inverse Lax-Wendroff (ILW) procedure in our works (Cheng et al., Appl Math Mech (Engl Ed) 42: 841–854, 2021; Liu et al.) for 2D problems. Different from the 2D case, the local coordinate rotation in 3D required in the ILW procedure is not unique. We give a theoretical analysis to show that the boundary treatment is independent of the choice of the rotation, ensuring the method is feasible and valid. Both translation and rotation of the body are taken into account in this paper. In particular, we reformulate the material derivative for inviscid fluid on the moving boundary with no-penetration condition, which plays a key role in the proposed algorithm. Numerical simulations on the cylinder and sphere are given, demonstrating the good performance of our numerical boundary treatments.

【 授权许可】

CC BY