A Direct Numerical Simulation (DNS) scheme for solving the temporal-spatial multiscale problem of sympathetic detonation is proposed. According to the physical and chemical properties, the SD process is divided into t...A Direct Numerical Simulation (DNS) scheme for solving the temporal-spatial multiscale problem of sympathetic detonation is proposed. According to the physical and chemical properties, the SD process is divided into two stages in the numerical simulation. Two different grid sizes are used in the two stages to improve calculation accuracy and efficiency. The local Level Set Method is used to accurately track the multi-material interface, and the Harten Lax and van Leer Contact (HLLC) method is used to solve local Riemann problem. Based on the computing method mentioned above, choosing the high evolution Weighted Essentially Non-Oscillatory (WENO) scheme for the spatial discretization of governing equations in alliance with Total Variation Diminishing (TVD) Runge-Kutta for time discretization, the Eulerian code of SD is developed. It solves the data interaction between the two solvers with variable grid size and time step, and realizes the three-dimensional parallel computing of SD. Composition B is taken as the research object. The processes of SD in different stand-off distances are analyzed. The transmission and safety distances of Composition B are given. The test under the same conditions is designed to verify the numerical simulation. The numerical simulation results are in good agreement with the experiment, which proves the accuracy and feasibility of the numerical method.展开更多
基金This work was supported by the National Natural Science Foundation of China[11732003,U1830139].
文摘A Direct Numerical Simulation (DNS) scheme for solving the temporal-spatial multiscale problem of sympathetic detonation is proposed. According to the physical and chemical properties, the SD process is divided into two stages in the numerical simulation. Two different grid sizes are used in the two stages to improve calculation accuracy and efficiency. The local Level Set Method is used to accurately track the multi-material interface, and the Harten Lax and van Leer Contact (HLLC) method is used to solve local Riemann problem. Based on the computing method mentioned above, choosing the high evolution Weighted Essentially Non-Oscillatory (WENO) scheme for the spatial discretization of governing equations in alliance with Total Variation Diminishing (TVD) Runge-Kutta for time discretization, the Eulerian code of SD is developed. It solves the data interaction between the two solvers with variable grid size and time step, and realizes the three-dimensional parallel computing of SD. Composition B is taken as the research object. The processes of SD in different stand-off distances are analyzed. The transmission and safety distances of Composition B are given. The test under the same conditions is designed to verify the numerical simulation. The numerical simulation results are in good agreement with the experiment, which proves the accuracy and feasibility of the numerical method.
