固定界面Level Set方法在爆轰波阵面传播中的应用

doi:10.3969/j.issn.1000-1093.2019.10.013

兵工学报 ›› 2019, Vol. 40 ›› Issue (10): 2080-2087.doi: 10.3969/j.issn.1000-1093.2019.10.013

固定界面Level Set方法在爆轰波阵面传播中的应用

张莉^1,2,3，吴开腾^1,2,3

（1.内江师范学院数学与信息科学学院，四川内江 641110；2.内江师范学院四川省高等学校数值仿真重点实验室，四川内江 641110；3.内江师范学院四川省数据恢复重点实验室，四川内江 641110）

收稿日期:2018-11-19 修回日期:2018-11-19 上线日期:2019-12-18
通讯作者: 吴开腾（1964—），男，教授，博士 E-mail:wukaiteng@263.net
作者简介:张莉（1981—），女，教授，博士。E-mail: lizhang0212@126.com
基金资助:
国家自然科学基金项目（11502121）

The Application of Fixed-interface Level Set Method in the Calculation of Detonation Wavefront Propagation

ZHANG Li^1,2，3， WU Kaiteng^1,2,3

（1.School of Maths and Informations Science, Neijiang Normal University, Neijiang 641110, Sichuan, China;2.Numerical Simulation Key Laboratory of Sichuan Province, Neijiang Normal University, Neijiang 641110, Sichuan, China；3.Data Recovery Key Laboratory of Sichuan Province, Neijiang Normal University, Neijiang 641110, Sichuan, China）

Received:2018-11-19 Revised:2018-11-19 Online:2019-12-18

摘要/Abstract

摘要： 针对重新初始化造成的界面位置偏离问题，提出一种固定界面Level Set方法。该方法的基本思想是在重新初始化相邻两次迭代过程中，通过构建保持界面不动满足的条件，从而推导出光滑参数的计算公式。与传统Level Set方法相比较，固定界面Level Set方法将固定的光滑参数修正为变化的光滑参数。数值模拟结果表明：固定界面Level Set方法在重新初始化过程中能够很好地固定界面位置，从而提高Level Set方法的质量守恒性；固定界面Level Set方法能够有效地避免由于误差累积导致的爆轰波阵面偏离现象，具有较高的计算精度和分辨率，适用于求解爆轰波阵面传播问题。

关键词: LevelSet方法, 固定界面, 光滑参数, 爆轰波阵面, 传播

Abstract: A fixed-interface Level Set method is proposed for the deviation of interface position caused by reinitialization. The basic idea of the fixed-interface Level Set method is to construct a condition that keeps the interface immobile in the process of reinitializing the two adjacent iterations, thus deducing the calculation formula for the smooth parameters. The proposed method is used to correct the fixed smooth parameters to the changed smooth parameters. The numerically simulated results show that, compared with the traditional Level Set method, the fixed-interface Level Set method can fix the interface position well in the process of reinitialization, and the mass conservation is improved. The fixed-interface Level Set method effectively avoides the deviation of detonation wavefront caused by error accumulation. It is suitable for solving the propagation problem of detonation wavefront. Key

Key words: LevelSetmethod, fixed-interface, smoothparameter, detonationwavefront, propagation

中图分类号:

O382⁺.1

张莉，吴开腾. 固定界面Level Set方法在爆轰波阵面传播中的应用[J]. 兵工学报, 2019, 40(10): 2080-2087.

ZHANG Li， WU Kaiteng. The Application of Fixed-interface Level Set Method in the Calculation of Detonation Wavefront Propagation[J]. Acta Armamentarii, 2019, 40(10): 2080-2087.

参考文献

［1］ MADER C L, FOREST C A. Two-dimensional homogeneous and heterogeneous detonation wave propagation［C］∥Proceedings of the 6th International Symposium on Detonation. Arlington, VA,US： Office of Naval Research, Department of the Navy, 1976：405- 413.
［2］ HALLQUIST J O. User's manul for DYNA2D［R］. San Francisco，CA，US: Lawrence Livermore National Laboratory, 1984: 18756-18763.
［3］ BDZIL J B, STEWART D S. Time-dependent two-dimensional detonation: the interaction of edge rarefactions with finite-length reaction zones［J］. Journal of Fluid Mechanics, 2006, 171:1-26.
［4］ BDZIL J B, STEWART D S. Modeling two-dimensional detonations with detonation shock dynamics［J］. Physics of Fluids A Fluid Dynamics, 1989, 1(7): 1261-1267.
［5］ BDZIL J, FICKETT W. DSD technology: a detonation reactive Huygens code: LA-12235-MS［R］. Los Alamos, NM，US: Los Alamos National Laboratory,1992.
［6］ ASLAMA T D. Level set methods applied to modeling detonation shock dynamics［M］. Boston, MA, US: Academic Press Professional, Inc.，1996.
［7］符松, 王智平, 张兆顺，等. 近水面水下爆炸的数值研究［J］. 力学学报, 1995, 27(3):267-276.
FU S, WANG Z P, ZHANG Z S, et al. Numerial study of underwater explosion near air-water surface［J］. Chinese Journal of Theoretical and Applied Mechanics, 1995, 27(3):267-276. (in Chinese)
［8］陈森华, 张旭, 柏劲松. 贴体坐标系中多维非理想爆轰波阵面传播计算的位标函数法［J］. 高压物理学报, 2000, 14(4):280-284.
CHEN S H, ZHANG X, BAI J S. LS method in body-fitted coordinate system for multi-dimendional non-ideal detonation propagation［J］. Journal of High Pressure Physics, 2000, 14(4):280-284. (in Chinese)
［9］文尚刚. 三维爆轰波传播的LS方法研究［D］. 绵阳：中国工程物理研究院, 2001:107-127.
WEN S G. Study on LS method of three-dimensional detonation wave propagation［D］. Mianyang: China Institute of Engineering Physics, 2001:107-127. (in Chinese)

［10］钟敏, 陈森华, 李平. 贴体坐标系中非理想爆轰波阵面传播计算的一个问题［J］. 高压物理学报, 2006, 20(4):397-402.
ZHONG M, CHEN S H, LI P. An issue in the numerical simulation for non-ideal detonation propagation in body-fitted coordinate systems［J］. Journal of High Pressure Physics, 2006, 20(4):397-402. (in Chinese)
［11］柏劲松, 李平, 钟敏,等. 以DSD理论和LS方法为基础的程序燃烧法［J］. 爆炸与冲击, 2008, 28(5): 401-406.
BAI J S, LI P, ZHONG M, et al. Aprogram burn based on detonation shock dynamics and level set methods［J］. Explosion and Shock Waves, 2008, 28(5): 401-406. (in Chinese)
［12］ WIGGINS D J R, FALLS S A E G. Two-dimensional detonation waves in Type Ia supernovae［J］. Monthly Notices of the Royal Astronomical Society, 2018, 287(3):575-582.
［13］ CLAVIN P. Nonlinear dynamics of shock and detonation waves in gases［J］. Combustion Science & Technology, 2018, 189(5):747-775.
［14］ CICCARELLI G, LI Q, METROW C.The three-dimensional structure of a detonation wave propagating in a round tube with orifice plates［J］. Shock Waves, 2018,28(5):1019-1030.
［15］ AVDEEV K A, AKSENOV V S, BORISOV A A, et al. Calculation of shock wave propagation in water containing reactive gas bubbles［J］. Russian Journal of Physical Chemistry B, 2017, 11(2): 261-271.
［16］马天宝, 陈建良, 宁建国. 爆轰波强间断问题的伪弧长算法及其人为解验证［J］. 爆炸与冲击, 2018, 38(2):271-278.
MA T B, CHEN J L, NING J G. Pseudo arc length algorithm and its artificial solution verification for detonation wave intensity discontinuity problem［J］. Explosion and Shock Waves, 2018, 38(2):271-278.(in Chinese)
［17］ PENG D, MERRIMAN B, OSHER S,et al. A PDE-based fast local level set method［J］. Journal of Computational Physics, 1999, 155(2):410-438.
［18］ SUSSMAN M, SMERKA P, OSHER S. A level set approach for computing solutions to incompressible two-phase flow［J］. Journal of Computational Physics, 1994, 114(1):146-159.
［19］张莉, 吴开腾, 宁建国. 一种固定界面的Level Set重新初始化方法［J］. 兵工学报, 2016, 37(增刊2):185-190.
ZHANG L, WU K T, NING J G. A fixed interface re-initialization Level Set method［J］. Acta Armamentarii, 2016, 37（S2）: 185-190.(in Chinese)
［20］ BDZIL J B. Steady-state two-dimensional detonation［J］. Journal of Fluid Mechanics, 2006, 108:195-226.

第40卷
第10期2019 年10月兵工学报ACTA
ARMAMENTARIIVol.40No.10Oct.2019

[1]	励明君，姜毅，马立琦，潘霄. 基于神经网络算法的发射场坪承载能力预测方法[J]. 兵工学报, 2022, 43(5): 982-991.
[2]	孟豪龙，翁春生，武郁文，郑权，肖强，王放，白桥栋. 环形燃烧室中凹腔对C₂H₄/Air旋转爆轰流场影响的数值模拟[J]. 兵工学报, 2022, 43(5): 1063-1074.
[3]	范超, 王鼎, 杨宾, 尹洁昕. 信号传播速度未知下水下多基地声纳定位算法[J]. 兵工学报, 2022, 43(3): 637-652.
[4]	吴明亮，郑权，续晗，翁春生，白桥栋. 氢气占比对氢气-煤油-空气旋转爆轰波传播特性的影响[J]. 兵工学报, 2022, 43(1): 86-97.
[5]	李宝星, 王中, 许桂阳，翁春生, 赵凤起. 煤油燃料旋转爆轰波起爆与传播特性实验研究[J]. 兵工学报, 2020, 41(7): 1339-1346.
[6]	高婕，牛刚，杜太行，雷正伟. 超短波多径信号相关干涉仪测向算法[J]. 兵工学报, 2020, 41(11): 2252-2259.
[7]	赵星宇，白春华，姚箭，孙彬峰. 燃料空气炸药爆轰产物JWL状态方程参数计算[J]. 兵工学报, 2020, 41(10): 1921-1929.
[8]	李军宝，李伟兵，汪衡，袁书强，王晓鸣，洪晓文，徐赫阳. 爆炸载荷下铝粉与橡胶复合材料中的冲击波传播特性[J]. 兵工学报, 2020, 41(10): 2001-2007.
[9]	林逸雪，狄长安，狄长春，宫新宇，计寒. 基于激波传播路径的弹头斜入射双三角阵定位模型[J]. 兵工学报, 2019, 40(2): 251-256.
[10]	何耀民，何华锋，徐永壮，苏敬，王依繁. 基于反向传播神经网络的海杂波参数估计[J]. 兵工学报, 2019, 40(12): 2473-2481.
[11]	邱文昊，连光耀，杨金鹏，黄考利. 基于多影响因子和重要度的故障样本优选[J]. 兵工学报, 2019, 40(12): 2551-2559.
[12]	魏万里，翁春生，武郁文，郑权，李宝星. 氧化剂喷注面积对连续旋转爆轰波传播特性影响的实验研究[J]. 兵工学报, 2018, 39(12): 2345-2353.
[13]	潘曦，佟颖，王华阳. 球谐域传播算子快速声定向算法[J]. 兵工学报, 2018, 39(10): 1936-1944.
[14]	高玉波，张伟, 李达诚，宜晨虹, 汤铁钢. 基于反向传播神经网络的陶瓷损伤参数反演分析[J]. 兵工学报, 2018, 39(1): 146-152.
[15]	黄大荣, 陈长沙，孙国玺，赵玲，米波. 复杂装备轴承多重故障的线性判别分析与反向传播神经网络协作诊断方法[J]. 兵工学报, 2017, 38(8): 1649-1657.

固定界面Level Set方法在爆轰波阵面传播中的应用

The Application of Fixed-interface Level Set Method in the Calculation of Detonation Wavefront Propagation

PDF

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 15

编辑推荐

Metrics

本文评价