一维爆轰波不稳定性的数值模拟

洪滔 秦承森

引用本文:
Citation:

一维爆轰波不稳定性的数值模拟

    通讯作者: 洪滔; 

Numerical Simulation of One-Dimensional Instability of Detonation Wave

    Corresponding author: HONG Tao
  • 摘要: 对当量比为1的乙烷与空气的混合气体的一维爆轰不稳定问题进行数值模拟,得到了不同大小的网格对爆轰不稳定问题数值模拟结果的影响。网格大小x由从ZND模型分析得到的导引长度Lin确定,网格大小从0.2Lin变到0.002 5Lin。随着网格变细,没有得到振幅趋于一致的解,每一种网格尺寸得到的解的振幅都互不相同。当网格大小为x=0.01Lin、0.005Lin时,得到有规则的爆轰激波阵面压力的振荡,振荡的模式是峰值一大一小的振荡。网格更细时,爆轰波的振荡在计算范围内由一些有规则的振荡和一些较不规则的振荡组成。但爆轰激波阵面压力振荡的波长最后趋于一致,为91~93 mm,与实验得到的胞格长度88 mm很接近。
  • [1] Fickett W, Davis W C. Detonation [M]. California: University of California Press, 1979: 230-288.
    [2] Lee H I, Stewart D S. Calculation of Linear Detonation Instability: One-Dimensional Instability of Plane Detonation [J]. J Fluid Mech, 1990, 216: 103-132.
    [3] Fickett W, Jacobson J D, Schott G L. Calculated Pulsating One-Dimensional Detonation with Induction-Zone Kinetics [J]. AIAA, 1972, 10: 514-516.
    [4] Abouseif G E, Toony T Y. Theory of Unstable One-Dimensional Detonations [J]. Combustion and Flame, 1982, 45: 68-94.
    [5] Bourlioux A, Majda A J, Roytburd V. Theretical and Numerical Structure for Unstable One-Dimensional Detonations [J]. SIAM J Appl Maths, 1991, 51: 303-343.
    [6] He L, Lee J H S. Dynamical Limit of One-Dimensional Detonations [J]. Phys Fluids, 1995, 7(5): 1151-1158.
    [7] Short M, Quirk J J. On the Nonlinear Stability and Detonability Limit of a Detonation Wave for a Model Three-Step Chain-Branching Reaction [J]. J Fluid Mech, 1997, 339: 89-119.
    [8] Westbrook C K, Dryer F L. Simplified Reaction Mechanisms for the Oxidation of Hydrocarbon [J]. Combust Sci Technol, 1981, 27: 31-43.
    [9] Kiehne T M, Matthews R D, Wilson D E. An Eight Step Kinetics Mechanism for High Temperature Propane Flames [J]. Combust Sci Technol, 1987, 54: 1-23.
    [10] Book D L, Boris J P, Hain K. Flux-Corrected Transport Ⅱ: Generalization of the Method [J]. J Comput Phys, 1975, 18: 248-283.
    [11] Taki S, Fujiwara T. Numerical Simulation of Triple Shock Behavior of Gaseous Detonation [A]. Eighteenth (International) Symposium on Combustion [C]. Pittsburgh, PA: The Combustion Institute, 1981: 1671-1681.
    [12] Tieszen S R, Stamps D W, Westbrook C K, et al. Gaseous Hydrocarbon-Air Detonations [J]. Combust Flame, 1991, 84: 376-390.
    [13] Lefebvre M H, Oran E S, Kailasanath K, et al. The Influence of Heat Capacity and Diluent on Detonation Structure [J]. Combust Flame, 1993, 95: 207-218.
    [14] Bull D C, Elsworth J E, Shuff P J. Detonation Cell Structures in Fuel/Air Mixtures [J]. Combust Flame, 1982, 45: 7-22.
  • [1] 胡湘渝张德良 . 气相爆轰波沿胞格的动力学机理研究. 高压物理学报, 2001, 15(3): 205-214 . doi: 10.11858/gywlxb.2001.03.007
    [2] 王刚张德良刘凯欣 . 气相爆轰波数值模拟中化学反应模型研究. 高压物理学报, 2008, 22(4): 350-356 . doi: 10.11858/gywlxb.2008.04.003
    [3] 董刚范宝春谢波 . 氢气-空气混合物中瞬态爆轰过程的二维数值模拟. 高压物理学报, 2004, 18(1): 40-46 . doi: 10.11858/gywlxb.2004.01.008
    [4] 张焕好陈志华孙晓晖范宝春 . 突扩管爆轰推进及其喷管推力性能数值模拟. 高压物理学报, 2010, 24(6): 415-422 . doi: 10.11858/gywlxb.2010.06.003
    [5] 姚干兵解立峰刘家骢 . 激波管内燃料两相云雾爆轰的实验研究. 高压物理学报, 2007, 21(1): 77-82 . doi: 10.11858/gywlxb.2007.01.013
    [6] 张博白春华 . C2H2-O2-Ar 混合气体爆轰特征参数研究. 高压物理学报, 2013, 27(2): 287-291. doi: 10.11858/gywlxb.2013.02.017
    [7] 陈默白春华刘庆明 . 铝粉-空气混合物爆燃转爆轰过程及爆轰波结构. 高压物理学报, 2012, 26(4): 426-432. doi: 10.11858/gywlxb.2012.04.011
    [8] 陈朗鲁建英冯长根 . 空气中激光支持爆轰波实验及理论分析. 高压物理学报, 2010, 24(2): 154-160 . doi: 10.11858/gywlxb.2010.02.012
    [9] 鲁建英陈朗伍俊英冯长根 . 激光支持等离子体爆轰波温度的实验测量. 高压物理学报, 2009, 23(2): 123-129 . doi: 10.11858/gywlxb.2009.02.008
    [10] 柏劲松李平陈森华廖海东杨礼兵姜洋 . 内爆加载下果冻内外界面不稳定性数值计算. 高压物理学报, 2004, 18(4): 295-301 . doi: 10.11858/gywlxb.2004.04.002
    [11] 袁帅胡海波张崇玉李平 . 爆轰波对碰驱动平面组合飞片的数值模拟. 高压物理学报, 2016, 30(1): 49-54. doi: 10.11858/gywlxb.2016.01.008
    [12] 罗艾民张奇 . 战斗部爆轰产物温度效应数值模拟. 高压物理学报, 2006, 20(1): 45-50 . doi: 10.11858/gywlxb.2006.01.010
    [13] 董玉斌苏林祥陈大年经福谦韩钧万封加波 . 滑移爆轰作用下内爆柱形钢壳层裂的数值模拟. 高压物理学报, 1989, 3(1): 1-10 . doi: 10.11858/gywlxb.1989.01.001
    [14] 曹玉忠卢泽生管怀安张幼平 . 抗爆容器内爆炸流场数值模拟. 高压物理学报, 2001, 15(2): 127-133 . doi: 10.11858/gywlxb.2001.02.009
    [15] 董贺飞赵艳红洪滔 . HMX炸药燃烧转爆轰数值模拟. 高压物理学报, 2012, 26(6): 601-607. doi: 10.11858/gywlxb.2012.06.001
    [16] 向梅黄毅民韩勇饶国宁彭金华 . JO-9159与JB-9014复合药柱爆轰驱动平面飞片实验与数值模拟. 高压物理学报, 2014, 28(3): 379-384. doi: 10.11858/gywlxb.2014.03.018
    [17] 倪小军马宏昊沈兆武蒋耀港李磊 . 流-固界面冲击波入射压力的测试及数值模拟. 高压物理学报, 2013, 27(4): 528-534. doi: 10.11858/gywlxb.2013.04.010
    [18] 段晓瑜郭学永聂建新王秋实焦清介 . RDX基含铝炸药三波点高度的数值模拟. 高压物理学报, 2018, 32(3): 035101-1-035101-8. doi: 10.11858/gywlxb.20170539
    [19] 贾宪振胡毅亭董明荣许学忠刘家骢 . 深水爆炸冲击波作用下圆柱壳动态响应影响因素的数值模拟研究. 高压物理学报, 2008, 22(2): 208-214 . doi: 10.11858/gywlxb.2008.02.016
    [20] 李顺波东兆星齐燕军孙舒 , . 水下爆炸冲击波在含吸收层结构中的传播规律的数值模拟. 高压物理学报, 2009, 23(5): 360-366 . doi: 10.11858/gywlxb.2009.05.007
  • 加载中
计量
  • 文章访问数:  1919
  • 阅读全文浏览量:  11
  • PDF下载量:  803
出版历程
  • 收稿日期:  2002-10-16
  • 录用日期:  2003-07-15
  • 刊出日期:  2003-12-05

一维爆轰波不稳定性的数值模拟

    通讯作者: 洪滔; 
  • 1. 北京应用物理与计算数学研究所,北京 100088

摘要: 对当量比为1的乙烷与空气的混合气体的一维爆轰不稳定问题进行数值模拟,得到了不同大小的网格对爆轰不稳定问题数值模拟结果的影响。网格大小x由从ZND模型分析得到的导引长度Lin确定,网格大小从0.2Lin变到0.002 5Lin。随着网格变细,没有得到振幅趋于一致的解,每一种网格尺寸得到的解的振幅都互不相同。当网格大小为x=0.01Lin、0.005Lin时,得到有规则的爆轰激波阵面压力的振荡,振荡的模式是峰值一大一小的振荡。网格更细时,爆轰波的振荡在计算范围内由一些有规则的振荡和一些较不规则的振荡组成。但爆轰激波阵面压力振荡的波长最后趋于一致,为91~93 mm,与实验得到的胞格长度88 mm很接近。

English Abstract

参考文献 (14)

目录

    /

    返回文章
    返回