基于多介质Riemann问题的流体-固体耦合数值方法及其在爆炸与冲击问题中的应用

doi:10.3969/j.issn.1000-1093.2021.02.012

摘要/Abstract

摘要： 为提高爆炸与冲击下多介质大变形问题的数值模拟精度，针对具有高度非线性特征、通用形式的Mie-Grüneisen状态方程和流体弹塑性本构方程的多种介质间相互作用问题，提出一种通用、健壮的多介质Riemann问题求解方法，以有效提高物质界面上各物理量的计算精度。结合Euler坐标系下具有锐利界面的守恒型多介质流动数值方法，建立一套能够模拟可压缩流体和弹塑性固体在极端物理条件下大变形动力学行为的数值计算体系。对流体-固体耦合、固体-固体耦合Riemann问题，地下强爆炸、空中强爆炸以及高速冲击等问题开展了数值模拟，计算结果和理论、实测数据符合较好，表明该数值方法能够对典型爆炸与冲击等多介质问题进行有效的模拟。

关键词: 爆炸与冲击, 多介质流动, Riemann问题, Mie-Grüneisen状态方程, 流体弹塑性模型

Abstract: In order to improve the accuracy of numerical simulation of multi-medium large deformation under explosive and impact loading, a general and robust multi-medium Riemann solver is proposed to solve the multi-medium interactions for the highly non-linear Mie-Grüneisen equation of state and hydro-elastoplastic constitutive equation. The solver can be used to improve the calculation accuracy of physical quantities at material interface. A numerical system is established to simulate the large deformation of compressible fluid and elastoplastic solid under extreme conditions by using the compressible multi-medium flow scheme. The fluid-solid Riemann problem, solid-solid Riemann problem, underground explosion problem, air blast problem and high-speed impact problem were simulated. The numerical results agree well with the analytical and experimental data, which reveals the capacity of the proposed algorithm in solving multi-medium flows, such as typical explosion and impact problems.

Key words: explosionandimpact, multi-mediumflow, Riemannproblem, Mie-Grüneisenequationofstate, hydro-elastoplasticmodel

中图分类号:

O381

姚成宝，付梅艳，韩峰，闫凯，雷雨. 基于多介质Riemann问题的流体-固体耦合数值方法及其在爆炸与冲击问题中的应用[J]. 兵工学报, 2021, 42(2): 340-355.

YAO Chengbao, FU Meiyan, HAN Feng, YAN Kai, LEI Yu. A Numerical Scheme for Fluid-solid Interactions Based on Multi-medium Riemann Problem and Its Application in Explosion andImpact Problems[J]. Acta Armamentarii, 2021, 42(2): 340-355.

参考文献

［1］ KAPAHI A, SAMBASIVAN S, UDAYKUMAR H. A three-dimensional sharp interface Cartesian grid method for solving high speed multi-material impact, penetration and fragmentation problems[J]. Journal of Computational Physics, 2013, 241: 308-332.
[2］ SAMBASIVAN S, KAPAHI A, UDAYKUMAR H. Simulation of high speed impact, penetration and fragmentation problems on locally refined Cartesian grids[J]. Journal of Computational Physics, 2013, 235: 334-370.
[3］ SAMBASIVAN S, UDAYKUMAR H. A sharp interface method for high-speed multi-material flows: strong shocks and arbitrary materialpairs[J]. International Journal of Computational Fluid Dynamics, 2011, 25(3): 139-162.
[4］ CHEN Q Y, LIU K X. A high-resolution Eulerian method for numerical simulation of shaped charge jet including solid-fluid coexistence and interaction[J]. Computers & Fluids, 2012, 56: 92-101.
[5］ 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]. Journal of Computational Physics, 2010, 229(19): 7503-7519.
[6］陈千一. 层裂与成型装药金属射流问题的数值分析[D]. 北京：北京大学, 2011.
CHEN Q Y. Numerical analysis of spall fracture and shaped charge jet[D]. Beijing: Peking University, 2011. (in Chinese)
[7］丁建许. 爆炸与冲击问题的高精度界面处理及数值研究[D]. 北京：北京理工大学, 2017.
DING J X. The high order inerfacial treatment and numerical investigation of explosion and impact problem[D]. Beijing: Beijing Institute of Technology, 2017. (in Chinese)
[8］王仲琦. 面向对象的爆炸力学 Euler 型多物质数值方法及其应用研究[D]. 北京：北京理工大学, 2000.
WANG Z Q. Object-oriented Euler multi-material numerical met-hods on and its applications[D]. Beijing: Beijing Institute of Technology, 2000. (in Chinese)
[9］ BANKS J R. On exact conservation for the Euler equations with complex equations of state[J]. Communications in Computational Physics, 2010, 8(5): 995.
[10］ KAMM J. Solution of the 1D Riemann problem with a general EOS in ExactPack[C]∥Proceedings of the 4th ASME Conference on Verification and Validation of Simulations. Las Vegas,NV, US：LANL, 2015.
[11］ DESPRES B. A geometrical approach to nonconservative shocks and elastoplastic shocks[J]. Archive for Rational Mechanics and Analysis, 2007, 186(2): 275-308.
[12］ GAVRILYUK S, FAVRIE N, SAUREL R. Modelling wave dynamics of compressible elastic materials[J]. Journal of Computational Physics, 2008, 227(5): 2941-2969.
[13］ KAMM J. FLAG simulations of the elasticity test problem of Gavrilyuk et al[R]. Los Alamos, NV,US: Los Alamos National Laboratory,2014.
[14］ LIN X, BALLMANN J. A Riemann solver and a second-order Godunov method for elasticplastic wave propagation in solids[J]. International Journal of Impact Engineering, 1993, 13(3): 463-478.
[15］ TRANGENSTEIN J, COLELLA P. A higher-order Godunov method for modeling finite deformation in elastic-plastic solids[J]. Communications on Pure and Applied Mathematics, 1991, 44(1): 41-100.
[16］ ABUZYAROV M, BAZHENOV V, KOTOV V. A Godunov-type method in dvnamics of elastoplastic media[J]. Computational Mathematics and Mathematical Physics, 2000, 40(6): 900-913.
[17］ BAZHENOV V, KOTOV V. Modification of Godunov's numerical scheme for solving problems of pulsed loading of soft soils[J]. Journal of Applied Mechanics and Technical Physics, 2002, 43(4): 603-611.
[18］ MENSHOV I, MISCHENKO A, SEREJKIN A. Numerical mo-deling of elastoplastic flows by the Godunov method on moving Eulerian grids[J]. Mathematical Models and Computer Simulations, 2014, 6(2): 127-141.
[19］ TANG H S, HUANG D. A second-order accurate capturing scheme for 1D inviscid flows of gas and water with vacuum zones[J]. Journal of Computational Physics, 1996, 128(2): 301-318.
[20］ LIU T G, XIE W F, KHOO B C. The modified ghost fluid method for coupling of fluid and structure constituted with hydro-elasto-plastic equation of state[J]. SIAM Journal on Scientific Computing, 2008, 30(3): 1105-1130.
[21］ GAO S, LIU T G. 1D exact elastic-perfectly plastic solid Riemann solver and its multimaterial application[J]. Advances in Applied Mathematics and Mechanics, 2017, 9(3): 621-650.
[22］ GAO S, LIU T G, YAO C B. A complete list of exact solutions for one-dimensional elasticperfectly plastic solid Riemann problem without vacuum[J]. Communications in Nonlinear Science and Numerical Simulation, 2018, 63: 205-227.
[23］ BERJAMIN H. Analytical solution to 1D nonlinear elastodynamics with general constitutive laws[J]. Wave Motion, 2017, 74: 35-55.
[24］ HATTORI H. The Riemann problem for thermoelastic materials with phase change[J]. Journal of Differential Equations, 2004, 205(1): 229-252.
[25］ FENG Z, KABOUDIAN A, RONG J. The simulation of compressible multi-fluid multi-solid interactions using the modified ghost method[J]. Computers & Fluids, 2017, 154(5): 12-26.
[26］ FENG Z, RONG J, KABOUDIAN A. The modified ghost method for compressible multimedium interaction with elastic-plastic solid[J]. Communications in Computational Physics, 2017, 22(5): 1258-1285.
[27］ TORO E F. Riemann solvers and numerical methods for fluid dynamics[M]. Berlin, Germany: Springer, 2009: 102-200.
[28］ FEHLBERG E. Low-order classical Runge-Kutta formulas with stepsize control and their application to some heat transfer problems: R-315[R]. Washington, DC,US: NASA, 1969.
[29］ MATHEWS J H, FINK K D. Numerical methods using MATLAB［M］. Upper Saddle River, NJ,US: Pearson Prentice Hall, 2004.
[30］ GUO Y H, LI R, YAO C B. A numerical method on Eulerian grids for two-phase compressible flow[J]. Advances in Applied Mathematics and Mechanics, 2016, 8(2):187-212.
[31］ DI Y N, LI R, TANG T. Level set calculations for incompressible two-phase flows on a dynamically adaptive grid[J]. Journal of Scientific Computing, 2007, 31(1):75-98.
[32］ LIU T G, XIE W F, KHOO B C. The modified ghost fluid method for coupling of fluid and structure constituted with hydro elastoplastic equation of state [J]. SIAM Journal on Scientific Computing, 2008, 30(3): 1105-1130.
[33］郝宝田. 地下核爆炸及其作用[M]. 北京：国防工业出版社, 2002.
HAO B T. Underground explosions and application [M]. Beijing: National Defense Industry Press, 2002. (in Chinese)
[34］姚成宝, 王宏亮, 浦锡锋,等. 空中强爆炸冲击波地面反射规律数值模拟研究[J]. 爆炸与冲击, 2019, 39(11): 112201-1-112201-8.
YAO C B, WANG H L, PU X F, et al. Numerical simulation of intense blast wave reflected on rigid ground[J]. Explosion and Shock Waves, 2019, 39(11): 112201-1-112201-8. (in Chinese)
[35］乔登江. 核爆炸物理概论[M]. 北京：国防工业出版社, 2003:51-55.
QIAO D J. An introduction to nuclear explosion physics［M］. Beijing: National Defense Industry Press, 2003:51-55. (in Chinese)
[36］ GLASSTONE S, DOLAN P J. The effects of nuclear weapons [M]. US: United States Department of Defense and the Energy Research and Development Administration, 1977: 453-501.
[37］江松，王瑞利. 多介质可压缩流体力学的若干测试基准问题[R]. 北京：北京应用物理与计算数学研究所, 2013.
JIANG S, WANG R L. Several benchmark problems of compressible multi-material fluid dynamics[R]. Beijing:Institute of Applied Physics and Computational Mathematics, 2013. (in Chinese)

[1]	罗永晨, 续晗, 张锋, 郑权, 翁春生. 乙烯对煤粉-氧气爆轰波起爆特性影响机制的实验研究[J]. 兵工学报, 2024, 45(3): 754-762.
[2]	李丹一, 程扬帆, 李翔, 王浩, 赵长啸, 沈兆武. Al/PTFE/TiH2三元活性材料与RDX组合装药的爆炸释能特性[J]. 兵工学报, 0, (): 0-0.
[3]	张锋, 罗永晨, 肖博文, 倪晓冬, 郑权, 翁春生, 续晗. 不同物理性质无烟煤粉在乙烯/氧气氛围下爆轰波传播特性实验研究[J]. 兵工学报, 0, (): 0-0.
[4]	夏全志, 吴艳青, 柴传国, 杨昆, 黄风雷. 冲击加载下PBX界面对热点形成和安全性影响[J]. 兵工学报, 0, (): 0-0.
[5]	梁霄, 王瑞利. 基于冲击Hugoniot关系的爆压性质和不确定度量化[J]. 兵工学报, 0, (): 0-0.
[6]	姚奎光，王淑娟，樊星，聂少云，王翔，代晓淦. 不同机械约束下压装PBX炸药反应演化行为[J]. 兵工学报, 2022, 43(8): 1772-1778.
[7]	李淑睿，段卓平，白志玲，张旭，黄风雷. 2，4-二硝基苯甲醚基熔铸含铝炸药冲击起爆特性[J]. 兵工学报, 2022, 43(6): 1288-1294.
[8]	段卓平，白志玲，白孟璟，黄风雷. 强约束固体炸药燃烧裂纹网络反应演化模型[J]. 兵工学报, 2021, 42(11): 2291-2299.
[9]	白志玲，段卓平，黄风雷. 高聚物粘结炸药冲击起爆统计热点反应速率模型[J]. 兵工学报, 2021, 42(11): 2379-2387.
[10]	李淑睿，段卓平，郑保辉，罗观，黄风雷. 24-二硝基苯甲醚基熔铸含铝炸药圆筒试验及爆轰产物状态方程[J]. 兵工学报, 2021, 42(7): 1424-1430.
[11]	魏万里，郑权，鲁江涛，翁春生，武郁文. 中心锥对液态燃料旋转爆轰发动机工作过程与性能的影响[J]. 兵工学报, 2021, 42(6): 1185-1194.
[12]	梁霄，陈江涛，王瑞利. 高维参数不确定爆轰的不确定度量化[J]. 兵工学报, 2020, 41(4): 692-701.
[13]	苗飞超，周霖，张向荣，曹同堂. 点火增长反应速率方程在LS-DYNA软件中嵌入及应用[J]. 兵工学报, 2019, 40(7): 1411-1417.
[14]	李瑞, 李伟兵, 靳洪忠, 王桂林, 洪晓文, 朱建军, 李文彬, 王晓鸣. 基于Jones-Wilkins-Lee状态方程的爆轰波相互作用参数理论分析[J]. 兵工学报, 2019, 40(3): 516-521.
[15]	魏万里，翁春生，武郁文，郑权，李宝星. 氧化剂喷注面积对连续旋转爆轰波传播特性影响的实验研究[J]. 兵工学报, 2018, 39(12): 2345-2353.