logo

Chinese Science Bulletin, Volume 60, Issue 10: 882-898(2015) https://doi.org/10.1360/N972014-00936

Progress in high-resolution numerical simulation of explosion mechanics

More info
  • AcceptedDec 16, 2014
  • PublishedMar 30, 2015

Abstract

nitride semiconductor material, microwave power device, power electronic device High-resolution numerical simulation of explosion problems is of significant theoretical value owing to the numerous potential applications in both national defense and civilian projects, such as ammunition design, damage assessment, target protection, and industrial gas explosion prevention. In this paper, we review the research history of explosion mechanics and present the highest-order numerical simulation in decades. Building upon the author's research work in recent years, we specifically discuss the development of high-order positivity-preserving schemes, high-order boundary treatment, an adaptive mesh method, and multi-material interface treatment. We extended the WENO and RKDG methods and constructed conservative positivity-preserving high-order schemes. The problems of negative density and pressure were solved through the simulation of complex detonation propagation. The detonation wave front structure and flow characteristics were captured clearly. In this paper, a high-order inverse Lax-Wendroff method on a Cartesian mesh is presented to treat the complex boundary. The main principle is that the first-order normal derivative can be obtained by the inverse Lax-Wendroff method, while all the higher-order normal derivatives can simply be obtained by fifth-order WENO-type extrapolation. The efficient implementation in the detonation process further indicates that the inverse Lax-Wendroff method is stable and robust when applied to a complex-geometry boundary. An a posteriori error estimate, which is often used as a refinement criterion, is applied to solve the Euler equations. To achieve mesh refinement, WENO and Hermite interpolations are used to prolong the solutions from coarse to fine grid in space and time, respectively. For mesh merging, the direct point value replacement method is employed to update the solutions of coarse grids. For this study, a local level-set tracking method was developed to treat the multi- material interface, and the velocity field in the computational domain was modified in order to effectively avoid errors caused by velocity field discontinuity when the interface position was solved. The method presented in this paper combines the GFM (Ghost Fluid Method) with the RGFM (Real Ghost Fluid Method) to deal with the material interface, thereby avoiding non-physical solutions when the GFM is used to deal with a strong discontinuity, and the CPU running time is increased when RGFM is employed to deal with a weak discontinuity in solving the Riemann problem repeatedly. Finally, the future of high-order numerical simulation of explosion mechanics is discussed.


References

[1] Harten A, Osher S. Uniformly high order accurate nonoscillatory schemes I. SIAM J Numer Anal,1987, 24: 279-309. Google Scholar

[2] Liu X D, Osher S, Tony C. Weighted essentially non-oscillatory schemes. J Comput Phys, 1994, 115: 200-212. Google Scholar

[3] Jiang G S, Shu C W. Efficient implementation of weighted ENO schemes. J Comput Phys, 1996, 126: 202-228. Google Scholar

[4] Zhang H X. The exploration of the spatial oscillations in finite difference solutions for navier-stokes shocks (in Chinese). Acta Aerodyn Sin, 1984, 1:12-19 [张涵信. 差分计算中激波上、下游出现波动的探讨. 空气动力学学报, 1984, 1: 12-19]. Google Scholar

[5] Zhang H X. Non-oscillatory dissipation and non-free-parameter difference scheme (in Chinese). Acta Aerodyn Sin, 1988, 6: 143-165 [张 涵信. 无波动、无自由参数的耗散差分格式. 空气动力学学报, 1988, 6: 143-165]. Google Scholar

[6] Shen M Y, Li H D, Liu Q S. Analytic discrete WENO-FCTscheme (in Chinese). Acta Aerodyn Sin, 1998, 16: 56-63 [沈孟育, 李海东, 刘秋生. 用解析离散法构造WENO-FCT格式. 空气动力学报, 1998, 16: 56-63]. Google Scholar

[7] Ren Y X, Liu Q S, Wang S P, et al. A high order accurate non-oscillating finite volume scheme using spline interpolation for solving hyperbolic conservation laws (in Chinese). Acta Aerodyn Sin, 1996, 14: 281-287 [任玉新, 刘秋生, 王少平, 等. 求解双曲型守恒律的高 精度、无波动样条逼近有限体积方法. 空气动力学学报, 1996, 14: 281-287]. Google Scholar

[8] Zong W G, Deng X G, Zhang H X. Double weighted essentially non-oscillatory shock-capturing schemes (in Chinese). Acta Aerodyn Sin,2003, 21: 218-224 [宗文刚, 邓小刚, 张涵信. 双重加权实质无波动激波捕捉格式. 空气动力学学报, 2003, 21: 218-224]. Google Scholar

[9] Zhu J, Zhao N. A kind of MWENO scheme and its applications (in Chinese). Acta Aerodyn Sin, 2005, 23: 330-333 [朱君, 赵宁. 一种 MWENO格式的构造与应用. 空气动力学学报, 2005, 23: 330-333]. Google Scholar

[10] Einfeldt B, Munz C D, Roe P L, et al. On Godunov-type methods near low densities. J Comput Phys, 1991, 92: 273-295. Google Scholar

[11] Linde T, Roe P L. Robust Euler codes. In: Proceedings of the 13th Computational Fluid Dynamics Conference, Snowmass Village, 1997, AIAA-97-2098. Google Scholar

[12] Dou H S, Khoo B C. Effect of initial disturbance on the detonation front structure of a narrow duct. Shock Waves, 2010, 20: 163-173. Google Scholar

[13] Wang C, Ma T B, Lu J. Influence of obstacle disturbance in a duct on explosion characteristics of coal gas. Sci China Phys Mech Astrom,2010, 53: 269-278. Google Scholar

[14] Cockburn B, Shu C W. The Runge-Kutta discontinuous Galerkin method for conservation laws V: Multidimensional systems. J Comput Phys, 1998, 141: 199-224. Google Scholar

[15] Shu C W, Osher S, Efficient implementation of essentially non-oscillatory shock-capturing schemes. J Comput Phys, 1988, 77: 439-471. Google Scholar

[16] Zhang X X, Shu C W, On positivity preserving high order discontinuous Galerkin schemes for compressible Euler equations on rectangular meshes. J Comput Phys, 2010, 229: 8918-8934. Google Scholar

[17] Zhang X X, Shu C W. Positivity-preserving high order discontinuous Galerkin schemes for compressible Euler equations with source terms. J Comput Phys, 2011, 230: 1238-1248. Google Scholar

[18] Zhang R, Zhang M, Shu C W, High order positivity-preserving finite volume WENO schemes for a hierarchical size-structured population model. J Comput Appl Math, 2011, 236: 937-949. Google Scholar

[19] Zhang X X, Shu C W, Positivity preserving high order finite difference WENO schemes for compressible Euler equations. J Comput Phys,2012, 231: 2245-2258. Google Scholar

[20] Cheng J, Shu C W, Positivity-preserving Lagrangian scheme for multi-material compressible flow. J Comput Phys, 2014, 257: 143-168. Google Scholar

[21] Wang C, Zhang X X, Shu C W, et al. Robust high order discontinuous Galerkin schemes for two-dimensional gaseous detonations. J Comput Phys, 2012, 231: 653-665. Google Scholar

[22] Nobuyuki T, Asahara M, Eto K, et al. Numerical simulation of spinning detonation in square tube. Shock Waves, 2008, 18: 329-344. Google Scholar

[23] Wang C, Bi Y, Han W, et, al. Large-scale parallel computing for 3D gaseous detonation. Procedia Eng, 2013, 61: 276-283. Google Scholar

[24] Hanana M, Lefebvre M H. Pressure profiles in detonation cells with rectangular and diagonal structures. Shock Waves, 2001, 11: 77-88. Google Scholar

[25] Moretti, G. Importance of boundary conditions in the numerical treatment of hyperbolic equations. Phys Fluids, 1969, 12: II 13-20. Google Scholar

[26] Peskin C S. Flow patterns around heart values: A numerical method. J Comput Phys, 1972, 10: 252-271. Google Scholar

[27] Farooq M A, Skøien A A, Müller B. Cartesian grid method for the compressible Euler equations using simplified ghost points treatments at embedded boundaries. Comput Fluids, 2013, 82: 50-62. Google Scholar

[28] Kreiss H O, Petersson N A, Yström J. Difference approximations for the second order wave equation. SIAM J Numer Anal, 2002, 40:1940-1967. Google Scholar

[29] Sjögreen B, Petersson N A. A Cartesian embedded boundary method for hyperbolic conservation laws. Commun Comput Phys, 2007, 2:1199-1219. Google Scholar

[30] Appelö D, Petersson N A. A fourth-order accurate embedded boundary method for the wave equation. SIAM J Sci Comput, 2012, 34:2982-3008. Google Scholar

[31] Huang L, Shu C W, Zhang M. Numerical boundary conditions for the fast sweeping high order WENO methods for solving the Eikonal equation. J Comput Math, 2008, 26: 336-346. Google Scholar

[32] Xiong T, Zhang M, Zhang Y T, et al. Fifth order fast sweeping WENO scheme for static Hamilton-Jacobi equations with accurate boundary treatment. J Sci Comput, 2010, 45: 514-536. Google Scholar

[33] Goldberg M, Tadmor E. Scheme-independent stability criteria for difference approximations of hyperbolic initial-boundary value problems I. Math Comput, 1978, 32: 1097-1107. Google Scholar

[34] Tan S R, Shu C W. Inverse Lax-Wendroff procedure for numerical boundary conditions of conservation laws. J Comput Phys, 2010, 229:8144-8166. Google Scholar

[35] Filbet F, Yang C. Inverse Lax-Wendroff method for boundary conditions of Boltzmann type models. J Comput Phys, 2013, 245: 43-61. Google Scholar

[36] Tan S R, Wang C, Shu C W, et al. Efficient implementation of high order inverse Lax-Wendroff boundary treatment for conservation laws. J Comput Phys, 2012, 231: 2510-2517. Google Scholar

[37] Vilar F, Shu C W. Development and stability analysis of the inverse Lax-Wendroff boundary treatment for central compact schemes. Esaim-Math Model Num, 2015, 49: 39-67. Google Scholar

[38] Wang R L, Liu Q, Wu Z H, et al. Advances in the adaptive grid methods on the area of science and engineering application (in Chinese). Sci Technol Rev, 2012, 30: 72-79 [王瑞利, 刘全, 吴子辉, 等. 科学与工程应用中自适应网格方法进展. 科学导报, 2012, 30: 72-79]. Google Scholar

[39] Berger M J, Oliger J. Adaptive mesh refinement for hyperbolic partial differential equations. J Comput Phys, 1984, 53: 484-512. Google Scholar

[40] Norman M L. The impact of AMR in numerical astrophysics and cosmology. Adaptive mesh refinement-theory and applications. In: Proceedings of the Chicago Workshop on Adaptive Mesh Refinement Methods, 2005, 2: 413-430. Google Scholar

[41] Trebotich D, Straalen B V, Graves D, et al. Performance of embedded boundary methods for CFD with complex geometry. J Comput Phys, 2008, 125: 012083. Google Scholar

[42] Khokhlov A M. Fully threaded tree algorithms for adaptive mesh fluid dynamics simulations. J Comput Phys, 1998, 143: 543-519. Google Scholar

[43] Ji H, Lien F S, Yee E. A new adaptive mesh refinement data structure with an application to detonation. J Comput Phys, 2010, 229:8981-8993. Google Scholar

[44] Berger M J, Saltzman J S. AMR on the CM-2. Appl Numer Math, 1994, 14: 239-253. Google Scholar

[45] Colella P, Crutchfield W Y. A parallel adaptive mesh refinement algorithm on the C-90. In Proceedings of the Energy Research Power Users Symposium, 1994, 9: 11-12. Google Scholar

[46] Rendleman C A, Beckner V E, Lijewski M, et al. Parallelization of structured hierarchical adaptive, mesh refinement algorithms. Comput Visual Sci 2000, 3: 137-147. Google Scholar

[47] Hendrickson B, Devine K. Dynamic load balancing in computational mechanics. Comput Method Appl M. 2000, 184: 485-500. Google Scholar

[48] Henshaw W D, Schwendeman D W. Parallel computation of three-dimensional flows using overlapping grids with adaptive mesh refinement. J Comput Phys, 2008, 227: 7469-7502. Google Scholar

[49] Teyssier R, Fromang S, Dormy E. Kinematic dynamos using constrained transport with high order Godunov schemes and adaptive mesh refinement. J Comput Phys, 2006, 218: 44-67. Google Scholar

[50] Fromang S, Hennebelle P, Teyssier R. A high order Godunov scheme with constrained transport and adaptive mesh refinement for astrophysical magnetohydrodynamics. Astron Astrophys, 2006, 457: 314-384. Google Scholar

[51] Baeza A, Mulet P. Adaptive mesh refinement techniques for high-order shock capturing schemes for multi-dimensional hydrodynamic simulations. Int J Numer Meth Fl, 2006, 52: 455-471. Google Scholar

[52] Shen C P, Qiu J M, Christlieb A. Adaptive mesh refinement on high order finite difference WENO scheme for multi-scale simulations. J Comput Phys, 2011, 230: 3780-3802. Google Scholar

[53] Dumbser M. High order space-time adaptive ADER-WENO finite volume schemes for non-conservative hyperbolic systems. Comput Method Appl M, 2014, 268: 395-387. Google Scholar

[54] Valorani M, Di Giacinto M. Performance assessment of an adaptive mesh refinement technique for detonation waves. In: Proceedings of 13th International Conference on Numerical Methods in Fluid Dynamics. Berlin: Springer Berlin Heidelberg, 1993. 300-304.. Google Scholar

[55] Valorani M, Di Giacinto M. Numerical simulation of high speed reactive flows with adaptive mesh refinement. In: Proceedings of the 4th International Conference on Hyperbolic Problems. Braunschweig: Friedr Vieweg & Sohn. 1993. 565-577. Google Scholar

[56] Quirk J J. A parallel adaptive grid algorithm for computational shock hydrodynamics. Appl Numer Math, 1996, 20: 427-453. Google Scholar

[57] Sharpe G J, Falle S A E G. Two-dimensional numerical simulations of idealized detonations. P Roy Soc Lond A, 2000, 456: 2081-2100. Google Scholar

[58] Loth E, Sivier S, Baum J. Adaptive unstructured finite element method for two-dimensional detonation simulations. Shock Waves, 1998,8: 47-53. Google Scholar

[59] Deiterding R. Detonation simulation with the AMROC framework. Lect Notes Comput Sci, 2005, 372: 916-927. Google Scholar

[60] Deiterding R. Parallel adaptive fluid-structure interaction simulation of explosions impacting on building structures. Comput Fluids, 2013, 88: 719-729. Google Scholar

[61] Ravindran P, Lu F K. Multi-dimensional adaptive simulation of shock-induced detonation. In: Proceedings of the 39th AIAA Fluid Dynamics Conference. Gyeongju, 2009. No. 3D-3. Google Scholar

[62] Cao J, Zhou G, Wang C, et al. A hybrid adaptive finite difference method powered by a posteriori error estimation technique. J Comput Appl Mathes, 2014, 259: 117-128.. Google Scholar

[63] Amsden A A. The partical-in cell method for the calculations of the dynamics of compressible fluids. Technical Report LA-3466, Los Alamos Scientific Laboratory, 1966. Google Scholar

[64] Gentry R A, Martin R E, Daly B J. An Eulerian differencing method for unsteady compressible flow problems. J Comput Phys, 1966, 1: 87-118. Google Scholar

[65] Hardlow F H, Welch J E. Numerical calculation of time-dependent viscous incompressible flow of fluid with free surface. Phys Fluids, 1965, 8: 21-82. Google Scholar

[66] Amsden A A, Harlow F H. The SMAC method: A numerical technique for calculating incompressible fluid flow. Technical Report LA-4370, Los Alamos Scientific Laboratory, 1970. Google Scholar

[67] Tome M F, McKee S. GENSMAC: A computational marker and cell method for free surface flows in general domains. J Comput Phys, 1994, 110: 171-186. Google Scholar

[68] Sulsky D, Chen Z, Schreyer H L. A particle method for history-dependent materials. Comput Methods Appl Mech Eng, 1994, 118: 179-196. Google Scholar

[69] Zhang X, Sze K Y, Ma S. An explicit material point finite element method for hyper-velocity impact. Int J Numer Meth Eng, 2006, 66: 689-706. Google Scholar

[70] Lian Y P, Zhang X, Liu Y. Coupling of finite element method with material point method by local multi-mesh contact method. Comput Meth Appl Mech Eng, 2011, 47: 3482-3494. Google Scholar

[71] Hirt C W, Nichols B D. Volume of fluid (VOF) method for the dynamics of free boundaries. J Comput Phys, 1981, 39: 201-225. Google Scholar

[72] Noh W F, Woodward P. SLIC(simple line interface calculation). In: Proceedings. of the 15th international conference on numerical methods in fluid dynamics. Springer Lecture Notes in Physics, 1976. 59: 330-340. Google Scholar

[73] Youngs D L. Time-dependent multi-material flow with large fluid distortion. In: Morton K W, Baines M J, eds. Numerical Methods for Fluid Dynamics. Pittsburgh: Academic press, 1982. 273-285. Google Scholar

[74] Ashgriz N, Poo J Y. FLAIR: Flux line-segment model for advection and interface reconstruction. J Comput Phys, 1991, 93: 449-468. Google Scholar

[75] Feng Q J, He C J, Zhang M, et al. Calculation of mixed cells in simulating high-speed collision problems with two-dimensional Eulerian hydrodynamic method (in Chinese). Chin J Comput Phys, 2003, 20: 147-152 [冯其京, 何长江, 张敏, 等. 二维高速碰撞问题欧拉数值 模拟的混合网格计算. 计算物理, 2003, 20: 147-152]. Google Scholar

[76] He C J, Yu Z L, Feng Q J. 3D Eulerian numerical simulation method of high speed impact (in Chinese). Explos Shock Waves, 1999, 19: 216-221 [何长江, 于志鲁, 冯其京. 高速碰撞的三维欧拉数值模拟方法. 爆炸与冲击, 1999, 19: 216-221]. Google Scholar

[77] Poo J Y , Ashgriz N. A computational method for determining curvatures. J Comput Phys, 1989, 84: 483-491. Google Scholar

[78] Dyadechko V, Shashkov M. Moment-of-fluid interface reconstruction. Technical Report LA-UR-05-7571, Los Alamos National Laboratory,2005. Google Scholar

[79] Jia Z P. A multi-material arbitrary lagrangian-eulerian method based on MOF interface reconstruction (in Chinese). Chin J Comput Phys, 2010, 27: 353-360 [贾祖朋. 基于MOF界面重构的多物质ALE方法. 计算物理, 2010, 27: 353-360]. Google Scholar

[80] Feng Q J, Liu J, Wang Y J, et al. Interface reconstruction based on volume and centroid-of-volume in two dimensional axisymmetry geometry (in Chinese). Sci China Phys Mech Astron, 2014, 44: 203-211 [冯其京, 刘军, 王言金, 等. 二维轴对称坐标下基于介质形心 和体积的界面重构方法. 中国科学: 物理学 力学 天文学, 2014, 44: 203-211]. Google Scholar

[81] Ning J G, Chen L W. Fuzzy interface treatment in Eulerian method. Sci China Ser E, 2004, 47: 550-568. Google Scholar

[82] Wang X, Ma T B, Ning J G. A pseudo arc-length method for numerical simulation of shock waves. Chin Phys Lett, 2014, 31: 030201. Google Scholar

[83] Ning J G, Ma T B, Fei G L. Multi-material Eulerian method and parallel computation for 3D explosion and impact problems. Int J Comp Meth, 2014, 11: 1350079. Google Scholar

[84] Ning J G, Ma T B, Lin G H. A grid generator for 3-D explosion simulations using the staircase boundary approach in Cartesian coordinates based on STL models. Adv Eng Softw, 2014, 67: 148-155. Google Scholar

[85] Osher S, Sethian J A. Fronts propagating with curvature dependent speed: Algorithms based on Hamilton-Jacobi formulations. J Comput Phys, 1988, 79: 12-49. Google Scholar

[86] Enright D, Fedkiw R, Ferziger J, et al. A hybrid particle level set method for improved interface capturing. J Comput Phys, 2002, 183: 92-116. Google Scholar

[87] Fedkiw R P, Aslam T, Merriman B, et al. A non-oscillatory eulerian approach to interfaces in multi-material flows (the ghost fluid method). J Comput Phys, 1999, 152: 457-492. Google Scholar

[88] Liu T G, Khoo B C, Yeo K S, et al. Underwater shock-free surface-structure interaction. Int J Numer Method Eng, 2003, 58: 609-630. Google Scholar

[89] Liu T G, Khoo B C, Xie W F. The modified ghost fluid method as applied to extreme fluid-structure interaction in the presence of cavitation. Comm Comput Phys, 2006, 1: 898-919. Google Scholar

[90] Wang C W, Liu T G, Khoo B C. A real ghost fluid method for the simulation of multimedium compressible flow. SIAM J Sci Comput, 2006, 28: 278-302. Google Scholar

[91] Qiu J X, Liu T G, Khoo B C. Runge-Kutta discontinuous Galerkin methods for compressible two-medium flow simulations: one-dimensional case. J Comput Phys, 2007, 222: 353-373. Google Scholar

[92] Qiu J X, Liu T G, Khoo B C. Simulations of compressible two-medium flow by Runge-Kutta discontinuous Galerkin methods with the ghost fluid method. Comm Comput Phys, 2008, 3: 479-504. Google Scholar

[93] Kim K H, Yoh J J. A particle level-set based Eulerian method for multi-material detonation simulation of high explosive and metal confinements. P Combust Inst, 2013, 34: 2025-2033. Google Scholar

[94] Tran L B, Udaykumar H S . A particle-level set-based sharp interface Cartesian grid method for impact, penetration, and void collapse. J Comput Phys, 2004, 193: 469-510. Google Scholar

[95] Wang J T, Liu K X, Zhang D L. An improved CE/SE schemes for multi-material elastic-plastic flows and its applications. Comput Fluids, 2009, 38: 544-551. Google Scholar

[96] Chen Q Y, Wang J T, Liu K X. Improved CE/SE scheme with particle level set method for numerical simulation of spall fracture due to high-velocity impact. J Comput Phys, 2010, 229: 7503-7519. Google Scholar

[97] Chen Q Y, Liu K X. A high-resolution Eulerian method for numerical simulation of shaped charge jet including solid-fluid coexistence and interaction. Comput Fluids, 2012, 56: 92-101. Google Scholar

[98] Chen S Y, Wu Z W, Liu K X, et al. Atomic diffusion behavior in Cu-Al explosive welding process. J Appl Phys, 2013, 113: 044901. Google Scholar

[99] Zhao H T, Wang C. High resolution numerical simulation of air explosion (in Chinese). Acta Armamentarii, 2013, 34: 1536-1546 [赵海 涛, 王成. 空中爆炸的高精度数值模拟研究. 兵工学报, 2013, 34: 1536-1546]. Google Scholar

[100] Zhao H T, Wang C, Ning J G. High resolution numerical simulation of compressible multi-medium flow (in Chinese). Chin J High Pressure Phys, 2013, 27: 261-267 [赵海涛, 王成, 宁建国. 可压缩多介质问题的高精度数值模拟. 高压物理学报, 2013, 27: 261-267]. Google Scholar

[101] Falcovitz J, Birman A. A singularities tracking conservation laws scheme for compressible duct flows. J Comput Phys, 1994, 115: 431-439. Google Scholar

[102] Hilditch J, Colella P. A front tracking method for compressible flames in one dimension. SIAM J Sci Comput, 1995, 16: 755-772. Google Scholar

[103] Glimm J, Graham M J, Grove J. Front tracking in two and three dimensions. J Comput Math, 1998, 35: 1-12. Google Scholar

[104] Du J, Fix B, Glimm J, et al. A simple package for front tracking. J Comput Phys, 2006, 213: 613-628. Google Scholar

[105] Qian J, Tryggvason G, Law C K. A front tracking method for the motion of premixed flames. J Comput Phys, 1998, 144: 52-69. Google Scholar

[106] Esmaeeli A, Tryggvason G. A front tracking method for computations of boiling in complex geometries. Int J Multiphas Flow, 2004, 30: 1037-1050. Google Scholar

[107] Mao D K. A shock tracking technique based on conservation in one space dimension. SIAM J Numer Anal, 1995, 32: 1677- 1703. Google Scholar

[108] Mao D K. Towards front-tracking based on conservation in tow space dimensions II, tracking discontinuities in capturing fashion. J Comput Phys, 2007, 226: 1550-1588. Google Scholar

[109] Wen W Z. A method of marker on cell line (MOCL) for tracking interfaces in Eulerian code (in Chinese). Acta Armamentarii, 1999, 20:212-216 [温万治. 欧拉程序中描述分界面的网格线示踪点法(MOCL). 兵工学报, 1999, 20: 212-216]. Google Scholar

[110] Wang C, Ning J G, Ma T B. Numerical simulation of detonation and multi-material interface tracking. Comput Mater Con, 2011, 22: 73-96. Google Scholar

Copyright 2020 Science China Press Co., Ltd. 《中国科学》杂志社有限责任公司 版权所有

京ICP备18024590号-1