• ISSN 1001-1455  CN 51-1148/O3
  • EI、Scopus、CA、JST收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊

障碍物对预混火焰特性影响的大涡数值模拟

王公忠 张建华 李登科 陈先锋

张林, 李英华, 张祖根, 李雪梅, 胡昌明, 蔡灵仓. 用于诊断材料微层裂的Asay窗技术[J]. 爆炸与冲击, 2017, 37(4): 692-698. doi: 10.11883/1001-1455(2017)04-0692-07
引用本文: 王公忠, 张建华, 李登科, 陈先锋. 障碍物对预混火焰特性影响的大涡数值模拟[J]. 爆炸与冲击, 2017, 37(1): 68-76. doi: 10.11883/1001-1455(2017)01-0068-09
Zhang Lin, Li Yinghua, Zhang Zugen, Li Xuemei, Hu Changming, Cai Lingcang. Asay window for probing the microspall of materials[J]. Explosion And Shock Waves, 2017, 37(4): 692-698. doi: 10.11883/1001-1455(2017)04-0692-07
Citation: Wang Gongzhong, Zhang Jianhua, Li Dengke, Chen Xianfeng. Large eddy simulation of impacted obstacles' effects on premixed flame's characteristics[J]. Explosion And Shock Waves, 2017, 37(1): 68-76. doi: 10.11883/1001-1455(2017)01-0068-09

障碍物对预混火焰特性影响的大涡数值模拟

doi: 10.11883/1001-1455(2017)01-0068-09
基金项目: 

国家自然科学基金项目 51174153

国家自然科学基金项目 51374164

国家重点研发计划项目 2016YFC0802801

建筑消防工程技术公安部重点实验室开放课题项目 KFKT2014ZD03

详细信息
    通讯作者:

    陈先锋,cxf618@whut.edu.cn

  • 中图分类号: O383

Large eddy simulation of impacted obstacles' effects on premixed flame's characteristics

  • 摘要: 障碍物在预混气体火焰传播过程中对其结构及传播特性造成较大影响,对火焰的加速和爆燃转爆轰过程(deflagration-to-detonation transition, DDT)起到直接的促进作用。通过障碍物条件下可视管道中甲烷/空气预混火焰传播实验,捕获其火焰微观结构变化。采用三维物理模型,采用壁面自适应局部涡黏模型(wall-adapting local eddy-viscosity, WALE)的大涡模拟(large eddy simulation, LES),并用火焰增厚化学反应模型(thickened flame model, TFM)对实验过程进行重现。分析开口管道中预混火焰翻越障碍物后的复杂流场变化,并分析层流向湍流转变过程的特点。揭示了在障碍物影响下预混火焰扰动失稳现象的直接原因,是由障碍物引发的3个气流涡团同时作用而形成Kelvin-Helmholtz不稳定及Rayleigh-Taylor不稳定现象耦合作用所导致。
  • 认识和理解材料在冲击、爆轰或强激光辐照等条件下的动态损伤与破坏具有重要的基础科学意义,在国防和民用工程领域也具有重要的应用价值。几十年来,材料在固态下的动态损伤与破坏得到了广泛关注,特别是层裂问题[1-12]。固体层裂可通过测量粒子速度剖面[13-14]、应力剖面[15]以及冻结回收样品等进行有效观测[1]。但是在强三角波或Taylor波等缺乏支持的强冲击波作用下,材料会发生一种更加极端的破坏。此时,材料在冲击压缩过程中会出现部分甚至完全熔化,当冲击波在材料自由面反射拉伸稀疏波时,由于材料的抗拉强度几乎完全丧失,材料将连续破碎为大量的细小液滴,以云雾形态向外高速喷射。这种在冲击熔化后发生的破碎现象称为“微层裂”[16]。对于微层裂问题,固体层裂实验中常用的VISAR(velocity interferometer system for any reflector)等激光干涉技术的作用非常有限,原因主要有两方面:第一,强冲击下冲击波传出样品几十纳秒以内材料表面便受到严重破坏,之后测试探头将无法获得有效的携带多普勒效应的信号光,从而无法完成对破碎过程的有效诊断,特别是对于厚样品而言(需要微秒量级的有效测量时间);第二,即使在固体层裂实验中,样品自由面速度剖面测量也只能对单次层裂开展,无法诊断样品的多层层裂(第2层及以后的层裂信息都被阻断),而样品的微层裂在形式上可视为连续的层层破碎,所以即使测得了有效的样品自由面速度剖面,也难以获得样品内部的破碎信息。正因为如此,当前的微层裂实验研究主要依靠照相技术(如X射线照相、质子照相、可见光阴影照相等)和回收分析技术[17-20]。通过X射线照相和质子照相,不仅可以获得产物形貌演化图像,还可以提取质量密度分布数据,但是受限于实验室诊断光源品质和动态模糊等问题,照相质量以及密度分布反演精度有待提升。此外,回收实验可能存在粒子的二次破碎,其统计价值也有待商榷。Asay窗是针对固体多层层裂问题建立的诊断技术[21],其基本思想是通过在距样品自由面一定间隙处设置窗口,利用激光干涉测速技术连续测量多层层裂片相继撞击窗口引起的窗口表面光反射信号的改变,反演分析样品的多层层裂特性。L.Signor等[22]、陈永涛等[23]以及本课题组[24]曾尝试利用Asay窗技术对微层裂问题进行诊断,图 1为报道的典型速度剖面[23-24],其中d为样品与窗口的间隙,v为弹速。图 1(a)中,速度剖面存在一定程度的杂乱干扰;而图 1(b)所示的剖面结构特征较为简单,缺乏特征变化。这不禁让人产生疑问:Asay窗技术对于微层裂问题的诊断能力究竟如何,诊断能否有效揭示样品微层裂的内在过程?对此,我们开展了多发实验,测得的速度剖面与图 1相比,具有更加丰富且明显的结构特征;波系作用分析表明,这些特征正是样品微层裂破碎不同特征区的反映。本文中将对上述疑问在一定程度上给予澄清,同时证明优化设计的Asay窗技术能够为微层裂研究提供有力且灵敏的诊断。

    图  1  较早报道的Asay窗诊断样品微层裂时所测的速度剖面
    Figure  1.  Earlier reported velocity profiles measured by Asay window in microspall experiments

    图 2为实验装置结构示意图。实验在冲击波物理与爆轰物理重点实验室的∅30 mm二级轻气炮上进行,由薄飞片对厚样品进行高速碰撞,实现样品在强三角波加载下的微层裂破碎。飞片击靶速度由DPS(Doppler pin system)测量,微层裂由Asay窗诊断。图 2中,Asay窗采用复合结构(由薄、厚两块窗口叠加,并在叠加面上蒸镀约1.0 μm厚铝膜)实现入射激光的反射,实际测量面是有限厚窗口的内界面。在L.Signor等[22]的测试方案中,速度测量面位于Asay窗左端面(碰撞表面),激光的反射通过粘贴在窗口表面的金属膜实现。这种测量方式存在明显缺陷:若膜太薄,易遭受破坏,测量的可靠性无法保证,即使测到信号,其信噪比一般也较低;若膜太厚,则膜本身的响应特性必须予以考虑,给实验解读增添了不利因素。实验装置的关键尺寸设计如下:弹托(Lexan聚碳酸酯)∅30 mm×10 mm,飞片(无氧铜)∅28 mm×0.5 mm,样品(Sn)∅22 mm×2.5 mm,Asay窗(LiF)∅22 mm×21.5 mm(其中薄窗口厚1.5 mm,厚窗口厚20.0 mm)。实验时,样品窗口间隙d可调,其他实验参数保持不变,用以分析间隙d对诊断效果的影响。原则上,如果该间隙过小,样品微层裂产物颗粒撞击Asay窗之前的运动时间较短,不同状态的产物颗粒在时空上的分辨不够明晰,导致产物颗粒的某些力学量分布信息无法测量与分析;反之,如果间隙过大,产物颗粒间距在撞击Asay窗之前已被过度拉开,某些力学量的特征变化在诊断时可能被抹掉或淹没在统计涨落与噪声之中。

    图  2  实验装置结构示意图
    Figure  2.  Experimental design

    共进行了5发实验,实验测量结果见表 1,测得的Asay窗粒子速度剖面见图 3表 1中:df为飞片厚度,ds为样品厚度。

    图  3  实测的Asay窗粒子速度剖面
    Figure  3.  Measured velocity profiles by Asay window for all experiments
    表  1  实验参数及结果
    Table  1.  Experimental parameters and results
    实验编号 df/mm ds/mm d/mm v/(km·s-1)
    Shot 1 0.508 2.522 0.3 3.244
    Shot 2 0.505 2.524 0.5 3.200
    Shot 3 0.503 2.522 1.0 3.254
    Shot 4 0.507 2.525 3.0 3.255
    Shot 5 0.500 2.521 5.0 3.252
    下载: 导出CSV 
    | 显示表格

    表 1所示的5发实验中弹速最高达3.255 km/s,最低为3.200 km/s,平均值为3.241 km/s,表明加载速度得到了很好的控制,可以近似认为5发实验中样品的微层裂过程是一致的。利用自行编写的流体动力学一维模拟程序以及我们早前建立的关于Cu和Sn的物态方程[25],对样品的冲击压缩过程进行模拟。结果显示,Sn样品的初始撞击压力最高达到73 GPa,冲击波传出样品自由面时的压力接近64 GPa,在整个厚度范围内样品压力均高于Sn的冲击熔化压力(约45 GPa[25]),表明Sn样品在冲击压缩过程中已完全进入熔化状态,之后所发生的破碎确属于微层裂。

    诊断的有效时间是实验关注的重要问题。一般而言,引起窗口失效有3种可能:(1)冲击波传出窗口自由面后破坏了自由面的透光性;(2)反射稀疏波与追赶稀疏波相互作用致使窗口破坏;(3)边侧稀疏波相互作用致使反光膜紊乱甚至窗口破坏。图 3显示,5发实验中Asay窗测试的有效时间不足2 μs,由于Asay窗的厚度达到20 mm,依据LiF的Hugoniot数据(c0=5.214 km/s,λ=1.336)[26],估算冲击波的传播时间至少接近3 μs,而Asay窗的直径为22 mm,边侧稀疏波到达窗口中心大约需要1.5 μs,考虑到反光膜的破坏尚需一定的孕育时间,由此可判定实验中Asay窗的测试有效时间主要由边侧稀疏决定。究其根源,主要在于轻气炮的口径有限。

    比较5发实验中Asay窗的测试剖面,发现在Shot 3(d=1.0 mm)、Shot 4(d=3.0 mm)和Shot 5(d=5.0 mm)的速度剖面前端均存在明显的小幅信号,而Shot 1(d=0.3 mm)和Shot 2(d=0.5 mm)的速度剖面前端没有或极不明显。分析认为,这部分信号是由冲击波到达样品自由面时引起的高速微喷粒子撞击Asay窗产生的,属于微喷信号,在此之后幅度较大的信号才属于微层裂信号。Shot 1和Shot 2中没有观测到明显微喷信号的原因在于间隙过小,微喷粒子与微层裂粒子到达Asay窗的时间间隔不够长,导致微喷粒子碰撞Asay窗产生的积分效应尚不明显。在我们的Asay窗诊断中,观测到了微喷信号,充分说明Asay窗诊断是灵敏的,可以利用Asay窗技术对破碎全过程(包含微喷和微层裂)进行诊断,但间隙的选择要适当。参照Shot 3的结果,设置的间隙大小应确保微喷信号持续100 ns以上。

    下面重点分析样品微层裂包含的特征过程,并结合实验结果考察Asay窗的实际诊断能力。图 4为本实验的波系相互作用原理示意图。飞片碰靶后,从O点分别向样品和飞片发出右行和左行冲击波OAOB;左行冲击波在B点反射一束追赶中心稀疏波,至A点追上样品中的冲击波,该冲击波随即衰减为三角波继续向前传播;三角波传至C点到达样品自由面并反射中心稀疏波,该稀疏波与追赶稀疏波相互作用致使样品破碎,即发生微层裂。进一步分析可以发现,实际上样品破碎存在着不同的特征分区。首先,对于最靠近自由面的CD间的区域,由于自由面为零压边界,因此该区域将保持固态,理论上会层裂为“痂片”,其厚度一般为数微米(与冲击波宽度同量级)。在细观不均匀性引起的复杂波系作用下,该痂片一般难以维持其层状结构,将破碎为颗粒。DE间的区域则完全是由一维应力波相互作用导致的破碎,破碎过程表现为一层紧接一层相继发生,但是其破碎边界实际沿着DFEF两条路径传播,并在F点相遇,即DF以及EF间的区域属于两个不同的破碎特征区。样品的其他区域,即EO间的区域,这里称为“残体”,不再发生微层裂,但是由于其处于熔化状态,依据T.de Rességuier等[27]的研究,后继在表面张力及边侧稀疏等作用下,同样会破碎为较大的颗粒。

    图  4  波系相互作用原理示意图
    Figure  4.  Wave interaction

    以上分析表明,理论上样品破碎会产生4个不同的特征区,分别为1个“痂片”特征区、2个微层裂特征区和1个“残体”特征区。D.B.Holtkamp等[19]采用质子照相技术对Sn样品在接触爆轰作用下的微层裂过程进行了研究,清晰地观测到总的微层裂区和“残体”区(见图 5),但未能揭示出微层裂的细致分区。观察图 3发现,Shot 1和Shot 2的速度剖面较完整地展现了样品的破碎过程,其特征点T1~T5将剖面分为4个特征信号段,与图 4所示的各特征区相对应:T1T2段为痂片特征区和微层裂DF特征区碰撞Asay窗产生的信号,T2T3为微层裂EF特征区产生的信号;T3T4为“残体”对应的信号;T4T5则为飞片碰撞Asay窗产生的信号。此结果表明,Asay窗技术能够对主要的特征区进行有效分辨,但是对于痂片特征区则不够灵敏。Asay窗未能有效诊断痂片特征区的原因可能是,痂片特征区太薄,其碰撞Asay窗的作用时间太短,致使作用积分效果不明显。对比图 3中不同间隙的测量结果,Shot 1、Shot 2和Shot 3显示,随着间隙d的增加,后继特征区对应的剖面信号相对前面特征区对应的信号呈现一定的爬升趋势。这一特点应如何理解?目前难以确定,需要通过理论建模和数值模拟深入分析。但是,Shot 5的结果则展示了问题的另一面,即当间隙d增大到一定程度后,除了由于有效测试时间的限制不能实现完整测量以外,更重要的是样品破碎各特征区的边界特征出现模糊,不再清晰,因此正如之前的基本判断,间隙d亦非越大越好,并且存在优化区间。观察本实验结果发现,对于2~3 mm厚的样品,若样品自由面速度在数千米每秒,则优化区间大约在0.5~3.0 mm范围。

    图  5  爆轰加载下Sn样品微层裂的质子照相图像[19]
    Figure  5.  Proton radiography of microspall event of Sn by explosive loading[19]

    Asay窗适用于微层裂问题的诊断,但目前关于其诊断能力、诊断信号特征及其所反映的微层裂内在过程尚缺乏清晰认识。通过将波系分析与实验相结合,考察了Asay窗对于飞片击靶导致微层裂破碎的诊断能力,分析了诊断信号的特征结构。结果表明:样品微层裂破碎存在1个“痂片”特征区、2个微层裂特征区以及1个“残体”特征区。实验结果表明,Asay窗对于样品微层裂的诊断非常灵敏,除了由于“痂片”特征区太薄难以分辨外,其余几个特征区均可以清晰地分辨。但是,诊断效果强烈依赖于Asay窗的优化设计,特别是样品与窗口的间隙尺寸,这些参数不仅影响有效诊断时间,更重要的是影响上述特征分区的辨析。另外,实验还揭示Asay窗可以同时对微喷进行有效诊断,并且同样对Asay窗与样品的间隙具有一定的要求,需要确保微喷信号持续约100 ns以上。

  • 图  1  管道中心轴切面结构图

    Figure  1.  Pipe axis structure in cross-section

    图  2  甲烷/空气预混火焰传播过程纹影图片

    Figure  2.  High-speed schlieren images of premixed methane/air flame propagation

    图  3  LES模型计算三维预混火焰结构时间序列图

    Figure  3.  Sequence diagram three-dimensional premixed flame structure by LES model

    图  4  火焰翻越障碍物时管道内流场结构示意图

    Figure  4.  The pipe flow field structure diagram when the flame climb over the obstacle

    图  5  特征点处压力时程曲线实验与LES模型计算值对比

    Figure  5.  Comparison of pressure histories at characteristic point between experiment and LES

    图  6  预混火焰翻越障碍物过程中管道内的速度矢量及流场变化图

    Figure  6.  Velocity vector and the flow field when premixed flame climbed over obstacle

    图  7  预混火焰法向截面的温度云图及流场结构

    Figure  7.  Temperature contours and the flow field structure of premixed flame normal section

  • [1] Ciccarelli G, Johansen C, Kellenberger M.High-speed flames and DDT in very rough-walled channels[J].Combustion and Flame, 2013, 160(1):204-211. doi: 10.1016/j.combustflame.2012.08.009
    [2] Chen Xianfeng, Zhang Yin, Zhang Ying.Effect of CH4-Air ratios on gas explosion flame microstructure and propagation behaviors[J].Energies, 2012, 5(10):4132-4146. doi: 10.3390/en5104132
    [3] Sklavounos S, Rigas F.Validation of turbulence models in heavy gas dispersion over obstacles[J].Journal of hazardous materials, 2004, 108(1):9-20. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=9c7b9ce98b1236bf2e9888f9b42111d6
    [4] Arntzen B J.Modelling of turbulence and combustion for simulation of gas explosions in complex geometries[J].Journal of Loss Prevention in the Process Industries, 1998, 18(4/5/6):225-237. http://cn.bing.com/academic/profile?id=bd294e5486d8a62e3f3850d02e580e62&encoded=0&v=paper_preview&mkt=zh-cn
    [5] Sarli V D, Benedetto A D, Long E J, et al.Time-resolved particle image velocimetry of dynamic interactions between hydrogen-enriched methane/air premixed flames and toroidal vortex structures[J].International Journal of Hydrogen Energy, 2012, 37(21):16201-16213. doi: 10.1016/j.ijhydene.2012.08.061
    [6] Masri A R, Ibrahim S S, Cadwallader B J.Measurements and large eddy simulation of propagating premixed flames[J].Experimental Thermal and Fluid Science, 2006, 30(7):687-702. doi: 10.1016/j.expthermflusci.2006.01.008
    [7] Gubba S R, Ibrahim S S, Malalasekera W, et al.Measurements and LES calculations of turbulent premixed flame propagation past repeated obstacles[J].Combustion and Flame, 2011, 158(12):2465-2481. doi: 10.1016/j.combustflame.2011.05.008
    [8] Kessler D A, Gamezo V N, Oran E S.Simulations of flame acceleration and deflagration-to-detonation transitions in methane-air systems[J].Combustion and Flame, 2010, 157(11):2063-2077. doi: 10.1016/j.combustflame.2010.04.011
    [9] Gamezo V N, Ogawa T, Oran E S.Flame acceleration and DDT in channels with obstacles:Effect of obstacle spacing[J].Combustion and Flame, 2008, 155(1/2):302-315. http://www.sciencedirect.com/science/article/pii/S0010218008001934
    [10] Ogawa T, Gamezo V N, Oran E S.Flame acceleration and transition to detonation in an array of square obstacles[J].Journal of Loss Prevention in the Process Industries, 2013, 26(2):355-362. doi: 10.1016/j.jlp.2011.12.009
    [11] 任少峰.可燃性气体泄爆动力学机理研究[D].武汉: 武汉理工大学, 2012. http://cdmd.cnki.com.cn/Article/CDMD-10497-1012442364.htm
    [12] Lesieur M.Large-eddy simulations of turbulence[M].Cambridge:Cambridge University Press, 2005.
    [13] 肖华华.管道中氢/空气预混火焰传播动力学实验与数值模拟研究[D].合肥: 中国科学技术大学, 2013. http://www.wanfangdata.com.cn/details/detail.do?_type=degree&id=Y2354259
  • 期刊类型引用(2)

    1. 张凤国,刘军,何安民,赵福祺,王裴. 强冲击加载下延性金属卸载熔化损伤/破碎问题的物理建模及其应用. 物理学报. 2022(24): 284-292 . 百度学术
    2. 贺年丰,任国武,陈永涛,郭昭亮. 爆轰加载下金属锡层裂破碎数值模拟. 爆炸与冲击. 2019(04): 91-98 . 本站查看

    其他类型引用(1)

  • 加载中
图(7)
计量
  • 文章访问数:  4933
  • HTML全文浏览量:  1749
  • PDF下载量:  428
  • 被引次数: 3
出版历程
  • 收稿日期:  2015-05-12
  • 修回日期:  2015-11-18
  • 刊出日期:  2017-01-25

目录

/

返回文章
返回