惰性气体和水蒸气对长直空间燃气爆炸超压及其振荡的抑制作用

刘洋 李展 方秦 王森佩 陈力

刘洋, 李展, 方秦, 王森佩, 陈力. 惰性气体和水蒸气对长直空间燃气爆炸超压及其振荡的抑制作用[J]. 高压物理学报, 2021, 35(5): 055201. doi: 10.11858/gywlxb.20200654
引用本文: 刘洋, 李展, 方秦, 王森佩, 陈力. 惰性气体和水蒸气对长直空间燃气爆炸超压及其振荡的抑制作用[J]. 高压物理学报, 2021, 35(5): 055201. doi: 10.11858/gywlxb.20200654
LIU Yang, LI Zhan, FANG Qin, WANG Senpei, CHEN Li. Inert Gas and Water Vapor Suppressing Overpressure and Its Oscillation of Gas Explosion in Long Straight Space[J]. Chinese Journal of High Pressure Physics, 2021, 35(5): 055201. doi: 10.11858/gywlxb.20200654
Citation: LIU Yang, LI Zhan, FANG Qin, WANG Senpei, CHEN Li. Inert Gas and Water Vapor Suppressing Overpressure and Its Oscillation of Gas Explosion in Long Straight Space[J]. Chinese Journal of High Pressure Physics, 2021, 35(5): 055201. doi: 10.11858/gywlxb.20200654

惰性气体和水蒸气对长直空间燃气爆炸超压及其振荡的抑制作用

doi: 10.11858/gywlxb.20200654
基金项目: 江苏省自然科学基金(BK20190571);国家自然科学基金(52008392)
详细信息
    作者简介:

    刘 洋(1996-),男,硕士研究生,主要从事燃气爆炸荷载研究. E-mail:20145035@cqu.edu.cn

    通讯作者:

    李 展(1990-),男,博士,讲师,主要从事燃气爆炸灾害效应研究. E-mail:lz.9008@163.com

  • 中图分类号: O381; TD712

Inert Gas and Water Vapor Suppressing Overpressure and Its Oscillation of Gas Explosion in Long Straight Space

  • 摘要: 长直空间燃气爆炸超压及其振荡将对人员和结构安全产生不利影响。为减轻燃气爆炸危害,基于CFD软件FLACS建立了长直管道空间燃气爆炸数值模型,并对模型进行了验证。利用已验证的数值模型,研究了添加不同体积分数CO2、N2和水蒸气的化学当量比CH4/空气混合气体的爆炸,讨论了惰性气体和水蒸气的体积分数对爆炸超压及其振荡的影响,并对比了3种气体的抑爆效果。结果表明:CO2、水蒸气和N2的体积分数每增加10%,密闭管道气体爆炸的最终超压将分别下降81、47、65 kPa,尾端泄爆管道分别下降24、25、20 kPa,3种气体的体积分数分别为25%、26%、30%时,爆炸被完全抑制;CO2、水蒸气和N2均能有效抑制爆炸超压的振荡,压力振幅和压力振荡频率均随添加气体体积分数的增加而减小;CO2对爆炸超压及其振荡的抑制效果最好,水蒸气次之,N2最弱,这与3种气体的物理特性及其抑爆机理的差异有关。

     

  • 长直空间具有横截面尺寸小、长度尺寸远大于截面尺寸、泄爆面积有限等特点,常见于煤矿巷道、隧道洞库、管涵暗渠、城市综合管廊等。在煤矿采挖、隧道施工以及油气管道输运过程中存在燃气泄漏的可能,泄漏的燃气在长直空间内与空气混合后遇点火源极易引发爆炸,进而造成人员伤亡和财产损失。近年来相关事故时有发生,例如:2013年青岛排水暗渠油气爆炸,2015年成都洛带古镇隧道瓦斯爆炸,2016年重庆金山沟煤矿瓦斯爆炸等。燃气爆炸产生的压力波是造成人员及结构损害的主要原因之一[1-2],有学者发现特定条件下长直空间燃气爆炸存在明显的超压振荡现象[3-8],而超压振荡将引起荷载峰值及升压速率的明显升高[9-10]、流场湍流程度增强[11],并可能造成结构高温蠕变[3]以及引发结构共振[12],从而进一步加剧爆炸对人员和结构的危害。考虑到长直空间燃气爆炸事故可能造成的严重后果,开展爆炸超压及其振荡抑制措施研究十分必要。

    Xing等[12]利用长度为2 m的圆管进行了甲烷/空气爆炸实验,发现CH4的体积分数为9.5%~12.0%时,密闭管道爆炸出现了剧烈的压力振荡;朱传杰等[13-14]通过实验和数值模拟(AutoReaGas)研究了不同尺寸密闭方管内CH4/空气爆炸超压的振荡特性,认为超压的振荡是压力波在管道两端产生反射的结果;韦世豪等[15-16]进行了小尺寸密闭管道内汽油蒸气/空气爆炸实验,认为前驱波和压缩波在容器内的反射叠加是产生压力振荡的主要原因,而火焰和压力波耦合所形成的正反馈现象进一步加剧了振荡。研究人员在长直尾端泄爆空间内同样观察到了压力振荡现象[6-8];Xing等[12]发现随着泄爆面积的增加,压力振荡相比于密闭管道逐渐衰减,但另外出现了亥姆霍兹振荡;李毅[17]、路长等[18]通过实验发现气体流入流出的交替变化对压力的振荡具有影响;王亚磊等[19]在长0.5 m的尾端管道内开展了燃气爆炸实验,结果表明管内气体的快速泄放和回流造成了压力的振荡。

    目前应用较为广泛的抑爆措施主要包括水系抑爆及惰性气体抑爆,二者分别通过在可燃气体中加入一定的水雾/水蒸气和惰性气体对爆炸产生抑制。用于抑爆的惰性气体主要为氮气(N2)和二氧化碳(CO2[20-21],实验发现随着N2和CO2浓度的增加,惰性气体对甲烷(CH4)爆炸的抑制作用增强[22-23],CO2和N2对甲烷/空气混合气体的抑爆浓度分别为22%~23%和30%~37%[24-26]。李成兵[27]、Xie等[28]认为CO2的比热容更高以及CO2作为燃烧产物能够从化学角度对燃烧产生抑制是CO2抑爆性能优于N2的主要原因。水系抑爆主要包括水雾抑爆和水蒸气抑爆,秦文茜等[29]、康泉胜等[30]分别通过实验研究了超细水雾对瓦斯和丙烷爆炸的抑制效果,发现添加超细水雾能够有效降低爆炸超压及其升压速率,并使爆炸下限大幅提高,随着超细水雾添加量的增加,抑制效果增强。高旭亮[31]利用Fluent软件研究了超细水雾对CH4爆炸的抑制效果,得到了与上述实验相同的结论,此外还发现当施加量达到一定程度时,超细水雾甚至会完全抑制CH4的爆炸。刘丹等[32]通过实验研究了水蒸气对CH4爆炸超压及爆炸极限的影响,发现水蒸气同样能够对爆炸产生抑制,且随着水蒸气浓度的增加,爆炸超压及升压速率降低,爆炸极限范围缩小。

    研究人员对长直空间燃气爆炸超压的振荡现象进行了研究,但大多关心超压的振荡特性及其产生原因,鲜有关于超压振荡抑制措施的研究报道。虽然添加N2、CO2和水蒸气能够在一定程度上对燃气爆炸超压及其升压速率产生抑制,但其能否抑制超压的振荡仍有待深入探讨。本研究拟通过FLACS数值模拟软件,建立长直密闭和长直尾端泄爆管道模型,并利用已有实验数据对模型进行验证;基于已验证的数值模型,以化学当量比甲烷/空气混合气体为研究对象,分析不同浓度CO2、N2及水蒸气对长直空间燃气爆炸超压及升压时间的抑制作用,对比3种气体对超压振荡的抑制效果,探究不同气体抑爆效果差异的原因。

    甲烷/空气混合气体的爆炸是一个快速释能的化学过程。FLACS描述可压缩流体流动的数学模型包括质量守恒方程、动量守恒方程以及焓的输运方程[33]

    质量守恒方程:

    t(βvρ)+xj(βjvρuj)=˙mV
    (1)

    动量守恒方程:

    t(βvρui)+xj(βjρuiuj)=βvpxi+pxi(βjσij)+Fo,i+Fw,i+βv(ρρ0)gi
    (2)

    焓的输运方程:

    t(βvρh)+xj(βjρujh)=xj(βjμeffσhhxj)+βvdpdt+˙QV
    (3)

    式中:t为时间;βv为体积孔隙率;ρ为密度;ρ0为初始密度;m为质量;β为火焰模型中的转换系数;V为比容;h为比焓;σh为焓对应的Prandtl数,取0.7; μeff为有效黏度;Q为热量流率;p为压强;uiuj分别为i方向和j方向的平均速度;σij为应力张量分量;gii方向的重力加速度;Fw,iFo,i分别为容器壁和亚网格阻塞引起的流动阻力。

    FLACS的湍流模型为标准k-ε 模型[33],属于涡黏模型中的两方程模型,湍动能k的运输方程为

    (ρk)t+(ρkui)xi=xj[(μ+μtσk)kxj]+Gk+GbρεYM+Sk
    (4)

    耗散率 ε 的运输方程为

    (ρε)t+(ρεui)xi=xj[(μ+μtσε)εxj]+C1εεk(Gk+C3εGb)C2ερε2k+Sε
    (5)

    式中:μ为动力黏度;μt为湍动黏度;Gk是由平均速度梯度引起的湍动能k的产生项;Gb是由浮力引起的湍动能k的产生项;YM代表可压湍流中脉动扩张的贡献;σk为湍动能k对应的Prandtl数,取为1.0;σε为耗散率ε对应的Prandtl数,取为1.3;SkSε为用户定义的源项;C1εC2εC3ε为经验常数,FLACS根据Launder等[34]的研究结果,分别取1.44、1.92、0.80。

    FLACS燃烧模型由火焰模型和燃烧速度模型组成[33]。火焰模型采用β模型,通过定义扩散系数D和反应速率RF确定燃烧标准和反应速率在数值火焰区域的空间分布。燃烧速度取准层流燃烧速度和湍流燃烧速度中的较大值。

    建模过程中对FLACS管道数值模型做以下简化:(1)模型管道壁面为刚性,不考虑其变形和损坏;(2)假设爆炸过程为绝热爆炸,不考虑由管道壁面散热引起的能量损失;(3)忽略模型管道壁上的各种可关闭预留孔口、传感器以及管道连接部位等细节,忽略管道壁对压力波的阻碍作用,认为容器内壁光滑。

    考虑到长直密闭和尾端泄爆管道内燃气爆炸均可能存在超压的振荡现象,建立如图1所示的完全密闭和尾端泄爆两种管道模型,两种管道内部尺寸均为0.8 m × 0.8 m × 30 m,尾端泄爆管道泄爆口尺寸为0.8 m × 0.8 m,泄爆压力为0 kPa。点火点设置在管道端部中央,分别在每节管道侧壁距点火点1、4、7、···、28 m每隔3 m处布置压力测点,依次编为测点1~测点10。流速尺度(CFLV)和声速尺度(CFLC)分别取为0.5和5.0,初始温度为20 ℃,初始压强100 kPa,相对湍流强度和湍流长度尺度均为零,其余设置均采用默认值。

    图  1  管道模型及其网格划分
    Figure  1.  Tube model and its meshing

    密闭模型各边界均设置为EULER边界,尾端泄爆管道气体泄放到开敞空间,故将除z轴负方向(地面)外的其他方向均设置为PLAN_WAVE边界。计算域及网格尺寸设置如下:将管道内部设置为核心区,核心区分别采用0.10、0.05、0.03 m 3种尺寸的网格。密闭管道核心区以外0.2 m范围为拉伸区,尾端泄爆管道核心区外x轴正负方向3 m,y轴正方向15 m,y轴负方向0.5 m,z轴正方向10 m,z轴负方向0.2 m为拉伸区,拉伸系数为1.2。

    研究人员对FLACS模拟惰性气体抑爆的可靠性进行了验证,陈晓坤等[35]利用FLACS软件对多组分可燃气体添加CO2后的爆炸参数进行了模拟研究,并将模拟结果与已有实验数据进行对比,发现不同配比条件下的最大爆炸超压模拟值与实验值的相对误差均小于16%,说明利用 FLACS 模拟多组分可燃性气体爆炸可行。罗振敏等[36]在内径88 mm、长1600 mm的尾端泄爆管道内进行了CO2抑制甲烷/空气爆炸实验,并通过FLACS软件对实验进行了模拟,结果表明,除CO2添加量为18%时误差较大外,模拟得到的火焰传播时间和火焰传播速度均与实验值基本吻合且变化趋势一致,认为误差主要来源于模型的理想环境假设以及FLACS对小尺寸管道模拟的偏差。张印等[37]利用FLACS软件对20 L密闭容器内添加CO2的CH4/空气预混爆炸进行了数值模拟,发现最大爆炸压力的模拟值与理论计算值的相对误差在6.5%~10.5%之间,爆炸温度的模拟值与计算值的相对误差在3%~13%之间。

    为验证FLACS模拟大尺寸长直管道燃气爆炸的准确性,对Li等[38]开展的长直空间燃气爆炸实验进行了模拟,模拟工况见图2。该实验管道尺寸为0.8 m × 0.8 m × 30 m,由10节长度为3 m的钢制管道组成。工况A和工况B均在管道端部充入3 m长的9.5%甲烷/空气混合气体,点火位置位于管道端部中央,点火能量100 mJ。工况A中管道完全密闭,工况C中第4、7、10节管道顶部泄爆口(尺寸0.6 m × 0.6 m)打开,泄爆压力0 kPa。需要指出的是,工况A的数值模型与完全密闭管道模型(图1(a))相同;工况C的数值模型除泄爆口数量与分布尾端泄爆管道模型(图1(b))不同外,其他设置均与之相同。

    图  2  模型验证工况[38]
    Figure  2.  Cases of model verification[38]

    工况A和工况C的模拟结果与实验值的对比如图3所示,以靠近点火端的测点1为例,可以看到密闭和泄爆条件下FLACS均能够正确地模拟爆炸超压随时间的变化规律。对于密闭管道(工况A),网格尺寸为0.10 m时模拟值与实验值最接近,第1个超压峰值(65 kPa)与实验值完全相同,第2个超压峰值(110 kPa)的相对误差为14.6%。对于泄爆管道(工况C),网格尺寸为0.05和0.03 m时误差较小且接近,超压峰值的误差分别为9.5%、11.9%。造成模拟值与实验值差异的主要原因:实验管道壁为散热良好的钢制粗糙壁面,而数值模型做了光滑、绝热的假设。模拟结果表明,数值模型基本可靠,模拟误差在可接受范围内。综合考虑结果的准确性和计算效率,后续密闭管道和泄爆管道分别采用0.10和0.05 m尺寸的网格进行计算。

    图  3  模型验证(测点1)
    Figure  3.  Numerical model verification (pressure sensor 1)

    为研究添加惰性气体及水蒸气对长直空间内燃气爆炸超压及其振荡的抑制效果,利用FLACS软件对添加了不同体积分数N2、CO2和水蒸气的化学当量比CH4/空气混合气体进行数值模拟,添加气体的体积分数如表1所示。充气体积如图4所示,在密闭管道内充满混合气体,在尾端泄爆管道点火端充入3 m长的混合气体。由于密闭管道尾端[14]和泄爆管道点火端[19]的爆炸超压及其振荡相较管道其他位置更明显,如未特别说明,本研究所讨论的密闭管道内的超压时程曲线均取自管道尾端测点(测点10),所讨论的尾端泄爆管道内的超压时程曲线均取自管道点火端测点(测点1)。

    表  1  添加气体的体积分数(密闭/泄爆)
    Table  1.  Volume fraction of added gas (closed/end-vented) % 
    CO2N2Water vapor (H2O) CO2N2Water vapor (H2O)
    000 202520
    454 243024
    8108 253126
    121512 27
    162016
    下载: 导出CSV 
    | 显示表格
    图  4  混合气体体积示意图
    Figure  4.  Diagram of premixed gas volume
    2.1.1   CO2的抑爆效果

    通过模拟发现,密闭管道内CO2的体积分数超过24%时不能点火,爆炸被完全抑制,与前人实验结果[24-26]非常接近。图5为添加不同体积分数CO2预混气体的爆炸超压时程曲线。可以发现,随着CO2体积分数的增加,超压的振荡现象逐渐减弱,振荡持续时间逐渐变短,当CO2添加量达到16%时,振荡完全消失。CO2对超压振荡的抑制可归因于其对混合气体燃烧速率的降低作用[27, 39]:首先,添加CO2稀释了甲烷及氧气,使反应速率降低,燃烧产生的热量减少;其次,CO2将吸收燃烧产生的热量使反应区温度降低,抑制链式反应的进行;最后,CO2作为燃烧产物之一将与燃烧区自由基反应,使自由基浓度降低,爆炸链传递速度减缓。密闭管道内超压的振荡主要是压力波来回反射传播的结果[13-16],燃烧速率的减缓将造成压力波强度(大小及升压速率)下降,使管内压力能够及时分布均匀,所以最终未表现出压力振荡现象。

    图  5  不同CO2体积分数条件下的超压时程曲线
    Figure  5.  Overpressure-time history curves of CO2 with different volume fractions

    图6所示,当CO2体积分数(φ)从零逐渐上升到24%时,升压时间(达到压力峰值的时间)从1.28 s逐渐增加到19.60 s,最终爆炸超压由785 kPa呈线性降低到590 kPa,即CO2的体积分数每增加10%,升压时间平均增加7.63 s,爆炸超压下降81 kPa。升压时间的增加是燃烧速率降低的结果,爆炸超压的降低则与添加CO2后管内燃料质量的降低有关。需要指出的是,未添加CO2时,超压峰值(897 kPa)受压力振荡的影响高于最终爆炸超压,但当CO2的体积分数达到4%时,受CO2抑爆作用的影响,超压的振荡减弱并不再使超压峰值升高,此后管内最终压力即为爆炸最大超压。

    图  6  最终超压及升压时间随CO2体积分数的变化
    Figure  6.  Change of final overpressure and pressure rising time with CO2 volume fraction

    超压振幅随CO2体积分数的变化见图7。可以看到:未添加CO2时压力振幅存在一个明显的上升阶段(0~0.96 s),该阶段点火后气体燃烧产生的压力波在管道内的往复反射形成了超压的振荡;燃烧与压力波的耦合作用使压力振幅不断升高[15-16],并在0.96 s达到最大值(337 kPa);0.96 s之后压力振幅进入下降阶段,此时燃料被大量消耗,预混气体浓度减小,燃烧速率降低,燃烧对压力波的增强逐渐小于压力波的摩擦损耗,振幅呈指数函数下降。需要指出的是,数值模型假设管道壁绝热,所以不考虑散热造成的压力振幅降低。随着CO2添加量的增加,超压振幅逐渐减小,振幅随时间的变化由升高和降低两个阶段逐渐变为仅有降低阶段,这是因为燃烧产生的压力波强度随着CO2体积分数的升高而降低,当CO2的体积分数达到一定值(约8%),燃烧对压力波的增强作用将小于气体摩擦损耗对压力波的衰减作用,振幅上升阶段消失。

    图  7  不同CO2体积分数下超压振幅随时间的变化
    Figure  7.  Variation of overpressure amplitude of different volume fractions of CO2 with time

    图8为不同CO2体积分数条件下超压振荡频率随时间的变化。可以看到,振荡频率随时间逐渐升高,并最终趋于一定值。此外,可以发现当CO2添加量增加时,超压振荡频率的上升速率和最大超压振荡频率降低。超压振荡归因于压力波的往复反射传播,在管道长度不变的情况下,超压振荡频率的大小主要取决于压力波的传播速度。在管道模型的绝热假设条件下,随着燃烧的进行,管内温度逐渐上升,并在燃烧结束后维持在一定值,而压力波传播速度与气体温度的1/2次方正相关[40],故超压振荡频率表现为先上升后维持的变化规律。另一方面,随着CO2体积分数的增加,甲烷的质量减小,燃烧产生的热量减少,造成管内温度降低,故最大超压振荡频率随CO2体积分数的增加而降低。

    图  8  不同CO2体积分数下超压振荡频率随时间的变化
    Figure  8.  Variation of overpressure oscillation frequency of different volume fractions of CO2 with time
    2.1.2   N2的抑爆效果

    图9为添加不同体积分数N2预混气体的爆炸超压时程曲线。模拟发现:N2的体积分数等于30%时,管内压力在10 s内均无明显变化;当N2的体积分数超过30%时,无法点火,可以认为模拟得到的N2完全抑爆浓度为30%,与已有实验结果[24-26]基本一致。图10显示了最终爆炸超压及升压时间随N2体积分数的变化。当N2的体积分数由零增加到25%时,升压时间从1.28 s逐渐增长到11.58 s,最终爆炸超压由785 kPa线性降低到667 kPa,即N2的体积分数每增加10%,升压时间平均增加4.12 s,爆炸超压下降47 kPa。随着N2体积分数的升高,超压的振荡逐渐减弱,振荡持续时间逐渐变短,N2的体积分数为20%时,振荡被完全抑制。

    图  9  不同N2体积分数条件下的超压时程曲线
    Figure  9.  Overpressure-time history curves of N2 with different volume fractions
    图  10  最终超压及升压时间随N2体积分数的变化
    Figure  10.  Change of final overpressure and pressure rising time with N2 volume fraction

    图11为超压振幅随N2体积分数的变化。随着N2添加量的增加,超压振幅逐渐降低,N2的体积分数达到10%时,振幅上升阶段消失,高于CO2(8%)。图12为超压振荡频率随N2体积分数的变化。随着N2添加量的增加,管内温度降低,压力波传播速度降低,超压振荡频率减小。N2主要通过稀释甲烷和氧气浓度抑制爆炸,其次是吸收燃烧产生的热量,相比CO2缺少了化学反应抑制一项,故N2对爆炸超压及其振荡的抑制作用弱于CO2[27, 39]

    图  11  不同N2体积分数下超压振幅随时间的变化
    Figure  11.  Variation of overpressure amplitude of different volume fractions of N2 with time
    图  12  不同N2体积分数下超压振荡频率随时间的变化
    Figure  12.  Variation of overpressure oscillation frequency of different volume fractions of N2 with time
    2.1.3   水蒸气的抑爆效果

    图13为不同水蒸气体积分数下甲烷/空气混合气体的爆炸超压时程曲线。可见,当水蒸气的体积分数达到16%时,超压振荡被完全抑制。此外,当水蒸气的体积分数达到26%时,管内压力在16 s内几乎未上升;当水蒸气的体积分数超过26%时,混合气体无法点火,由此可知模拟得到的水蒸气完全抑爆浓度为26%。

    图  13  不同水蒸气体积分数条件下的超压时程曲线
    Figure  13.  Overpressure-time history curves of N2 with different volume fractions

    图14所示,当水蒸气的体积分数由零增加到24%时,升压时间从1.28 s逐渐增长到13.79 s,最终爆炸超压由785 kPa线性降低到629 kPa,即水蒸气的体积分数每增加10%,升压时间平均增加5.21 s,爆炸超压下降65 kPa。

    图  14  最终超压及升压时间随水蒸气体积分数的变化
    Figure  14.  Change of final overpressure and pressure rising time with water vapor volume fraction

    图15图16分别显示了超压振幅和超压振荡频率随水蒸气体积分数的变化。与添加CO2和N2类似,当水蒸气的体积分数约为8%时,超压振幅的上升阶段消失,随着水蒸气体积分数的继续增大,最大超压振幅减小,超压振荡频率降低。

    图  15  不同水蒸气体积分数下超压振幅随时间的变化
    Figure  15.  Variation of overpressure amplitude of different volume fractions of water vapor with time
    图  16  不同水蒸气体积分数下超压振荡频率随时间的变化
    Figure  16.  Variation of overpressure oscillation frequency of different volume fractions of water vapor with time
    2.2.1   CO2的抑爆效果

    图17为不同CO2体积分数条件下混合气体爆炸超压时程曲线。可以发现,当尾端泄爆管道内CO2的体积分数达到24%时,混合气体仅发生缓慢燃烧,管内压力不发生明显变化;当CO2的体积分数为25%时,气体将不能点火,可以认为,尾端泄爆管道CO2的完全抑爆浓度为24%,与密闭管道相同。

    图  17  不同CO2体积分数条件下的超压时程曲线
    Figure  17.  Overpressure-time history curves of CO2 with different volume fractions

    图18所示,随着CO2的体积分数由零增长到24%,升压时间从0.24 s逐渐增长到0.71 s,超压峰值由57 kPa接近线性降低到0 kPa,即CO2的体积分数每增加10%,爆炸超压下降约24 kPa。

    图  18  最终超压和升压时间随CO2体积分数的变化
    Figure  18.  Change of final overpressure and pressure rising time with CO2 volume fraction

    图19为爆炸超压振幅随时间的变化。可见,不同CO2体积分数条件下混合气体的压力振幅均随时间呈指数函数下降。泄爆管道内压力的振荡类似于阻尼振动[41],振荡可归因于管内气体过度泄放引起气体往复的泄出与回流[19],因此可以认为,压力振幅的下降与振荡过程中气体动能被耗散有关。此外,随着CO2的体积分数由零增加到20%,最大压力振幅由90 kPa减小为7 kPa,这是CO2导致的混合气体燃烧速率降低、爆炸超压和管内气体流速减小的结果。

    图  19  不同CO2体积分数下超压振幅随时间的变化
    Figure  19.  Variation of overpressure amplitude of different volume fractions of CO2 with time

    图20给出了压力振荡频率随时间的变化。可以发现,振荡频率随CO2体积分数的增加而降低,可解释为:CO2的体积分数越大,燃烧速率越低,管内气体通过燃烧获取的速度越小,气体往复运动一个周期所需时间更长,最终导致超压振荡频率降低。不考虑由超压时程曲线得到压力振荡频率的误差,可以发现,各工况下的压力振荡频率几乎不随时间而明显上升或下降,如果认为尾端泄爆管道内压力的振荡为阻尼振动,则压力振荡频率的不变可能与阻尼振动的等时性[41]有关。

    图  20  不同CO2体积分数下超压振荡频率随时间的变化
    Figure  20.  Variation of overpressure oscillation frequency of different volume fractions of CO2 with time
    2.2.2   N2的抑爆效果

    图21为不同N2体积分数下混合气体爆炸超压时程曲线。可以发现,当N2的体积分数达到30%时,混合气体仅发生缓慢燃烧,管内压力无明显变化;当N2的体积分数为31%时,爆炸被完全抑制,说明尾端泄爆管道内N2的完全抑爆浓度与密闭管道相同,均为30%。

    图  21  不同N2体积分数条件下的超压时程曲线
    Figure  21.  Overpressure-time history curves of N2 with different volume fractions

    图22所示,当N2的体积分数由零增加到30%,升压时间由0.24 s增长到0.79 s,超压峰值由57 kPa接近线性降低到0 kPa,即CO2的体积分数每增加10%,爆炸超压下降20 kPa。

    图  22  最终超压及升压时间随N2体积分数的变化
    Figure  22.  Change of final overpressure and pressure rising time with N2 volume fraction

    压力振幅和压力振荡频率随时间的变化分别如图23图24所示。可以看到,随着N2体积分数的升高,超压振幅和振荡频率均逐渐减小,说明N2对尾端泄爆管道内爆炸超压的振荡效应具有一定的抑制作用。

    图  23  不同N2体积分数下超压振幅随时间的变化
    Figure  23.  Variation of overpressure amplitude of different volume fractions of N2 with time
    图  24  不同N2体积分数下超压振荡频率随时间的变化
    Figure  24.  Variation of overpressure oscillation frequency of different volume fractions of N2 with time
    2.2.3   水蒸气的抑爆效果

    图25为不同水蒸气体积分数下甲烷/空气混合气体的爆炸超压时程曲线。可以发现,尾端泄爆管道内水蒸气完全抑制爆炸的体积分数与密闭管道相同,均为26%。

    图  25  不同水蒸气体积分数条件下的超压时程曲线
    Figure  25.  Overpressure-time history curves of water vapor with different volume fractions

    图26可以看出,随着水蒸气体积分数由零增加到24%,升压时间从0.24 s逐渐增长到0.68 s,最终爆炸超压由57 kPa接近线性降低到0 kPa,即水蒸气的体积分数每增加10%,爆炸超压下降约25 kPa。

    图  26  最终超压及升压时间随水蒸气体积分数的变化
    Figure  26.  Change of final overpressure and pressure rising time with water vapor volume fraction

    图27图28分别显示了不同水蒸气体积分数下超压振幅和振荡频率随时间的变化。可以看出,振幅和频率均随着水蒸气含量的增加而减小,说明添加水蒸气同样能够有效抑制尾端泄爆管道内甲烷爆炸超压的振荡效应。

    图  27  不同水蒸气体积分数下超压振幅随时间的变化
    Figure  27.  Variation of overpressure amplitude of different volume fractions of water vapor with time
    图  28  不同水蒸气体积分数下超压振荡频率随时间的变化
    Figure  28.  Variation of overpressure oscillation frequency of different volume fractions of water vapor with time

    由2.1、2.2节可知,CO2、水蒸气和N2的完全抑爆浓度分别为25%、26%、30%,且泄爆条件对完全抑爆浓度没有影响。密闭条件下完全抑制超压振荡时CO2、水蒸气和N2的体积分数分别为16%、16%、20%,尾端泄爆条件下压力振荡被完全抑制时的体积分数则与完全抑爆浓度相同,为25%、26%、30%。图29图30分别为添加CO2、水蒸气和N2对长直空间内甲烷爆炸超压抑制作用的对比。可以发现,当添加气体浓度相同时,密闭和尾端泄爆管道均是添加CO2的爆炸超压及其振幅下降最明显,升压时间增加最长,添加水蒸气次之,添加N2最差,说明CO2对爆炸超压及其振荡的抑制作用最强,N2最弱。

    图  29  密闭管道内CO2、N2和水蒸气对爆炸的抑制作用对比
    Figure  29.  Comparison of suppression of CO2, N2 and water vapor on the explosion in the closed tube
    图  30  尾端泄爆管道内CO2、N2和水蒸气对爆炸的抑制作用对比
    Figure  30.  Comparison of suppression of CO2, N2 and water vapor on the explosion in the end-vented tube

    3种气体抑爆效果的差异可解释为[27, 39]:N2仅能从物理角度(稀释、吸热)对超压及其振荡产生抑制,而CO2与水蒸气作为甲烷爆炸的燃烧产物,除物理抑爆外,还能够与燃烧区自由基反应,降低自由基浓度,抑制爆炸链传递速度,进而从化学角度对爆炸产生抑制,故抑爆效果优于N2。水蒸气和CO2抑爆效果的差异可能与二者比热容不同及参与链式反应的种类和程度不同有关。

    利用FLACS软件建立了长直(密闭/尾端泄爆)空间管道模型,并对模型进行了验证。基于已验证的模型研究了添加CO2、N2和水蒸气的化学当量比CH4/空气混合气体的爆炸情况,讨论了添加不同浓度气体对爆炸超压及其振荡的抑制作用,并对比了3种气体的抑爆效果。

    (1)管道内气体爆炸超压随添加气体体积分数的增加而降低,CO2、水蒸气和N2的体积分数每增加10%,密闭管道气体爆炸的最终超压分别下降81、65、47 kPa,尾端泄爆管道超压峰值分别下降24、25、20 kPa,3种气体的完全抑爆浓度分别为25%、26%和30%(体积分数)。

    (2)CO2、水蒸气和N2对长直空间燃气爆炸超压的振荡具有抑制作用,压力振幅和压力振荡频率均随添加气体体积分数的增加而减小,密闭管道内3种气体完全抑制超压振荡时的体积分数分别为16%、16%、20%,尾端泄爆管道分别为25%、26%、30%。

    (3)CO2对爆炸超压及其振荡的抑制效果最好,水蒸气次之,N2最弱。这是因为N2仅能从物理角度对爆炸产生抑制,而CO2与水蒸气作为甲烷爆炸的燃烧产物,还能通过参与链式反应而抑制爆炸。水蒸气和CO2抑爆效果的差异可能与二者比热容以及参与链式反应的种类和程度不同有关。

  • 图  管道模型及其网格划分

    Figure  1.  Tube model and its meshing

    图  模型验证工况[38]

    Figure  2.  Cases of model verification[38]

    图  模型验证(测点1)

    Figure  3.  Numerical model verification (pressure sensor 1)

    图  混合气体体积示意图

    Figure  4.  Diagram of premixed gas volume

    图  不同CO2体积分数条件下的超压时程曲线

    Figure  5.  Overpressure-time history curves of CO2 with different volume fractions

    图  最终超压及升压时间随CO2体积分数的变化

    Figure  6.  Change of final overpressure and pressure rising time with CO2 volume fraction

    图  不同CO2体积分数下超压振幅随时间的变化

    Figure  7.  Variation of overpressure amplitude of different volume fractions of CO2 with time

    图  不同CO2体积分数下超压振荡频率随时间的变化

    Figure  8.  Variation of overpressure oscillation frequency of different volume fractions of CO2 with time

    图  不同N2体积分数条件下的超压时程曲线

    Figure  9.  Overpressure-time history curves of N2 with different volume fractions

    图  10  最终超压及升压时间随N2体积分数的变化

    Figure  10.  Change of final overpressure and pressure rising time with N2 volume fraction

    图  11  不同N2体积分数下超压振幅随时间的变化

    Figure  11.  Variation of overpressure amplitude of different volume fractions of N2 with time

    图  12  不同N2体积分数下超压振荡频率随时间的变化

    Figure  12.  Variation of overpressure oscillation frequency of different volume fractions of N2 with time

    图  13  不同水蒸气体积分数条件下的超压时程曲线

    Figure  13.  Overpressure-time history curves of N2 with different volume fractions

    图  14  最终超压及升压时间随水蒸气体积分数的变化

    Figure  14.  Change of final overpressure and pressure rising time with water vapor volume fraction

    图  15  不同水蒸气体积分数下超压振幅随时间的变化

    Figure  15.  Variation of overpressure amplitude of different volume fractions of water vapor with time

    图  16  不同水蒸气体积分数下超压振荡频率随时间的变化

    Figure  16.  Variation of overpressure oscillation frequency of different volume fractions of water vapor with time

    图  17  不同CO2体积分数条件下的超压时程曲线

    Figure  17.  Overpressure-time history curves of CO2 with different volume fractions

    图  18  最终超压和升压时间随CO2体积分数的变化

    Figure  18.  Change of final overpressure and pressure rising time with CO2 volume fraction

    图  19  不同CO2体积分数下超压振幅随时间的变化

    Figure  19.  Variation of overpressure amplitude of different volume fractions of CO2 with time

    图  20  不同CO2体积分数下超压振荡频率随时间的变化

    Figure  20.  Variation of overpressure oscillation frequency of different volume fractions of CO2 with time

    图  21  不同N2体积分数条件下的超压时程曲线

    Figure  21.  Overpressure-time history curves of N2 with different volume fractions

    图  22  最终超压及升压时间随N2体积分数的变化

    Figure  22.  Change of final overpressure and pressure rising time with N2 volume fraction

    图  23  不同N2体积分数下超压振幅随时间的变化

    Figure  23.  Variation of overpressure amplitude of different volume fractions of N2 with time

    图  24  不同N2体积分数下超压振荡频率随时间的变化

    Figure  24.  Variation of overpressure oscillation frequency of different volume fractions of N2 with time

    图  25  不同水蒸气体积分数条件下的超压时程曲线

    Figure  25.  Overpressure-time history curves of water vapor with different volume fractions

    图  26  最终超压及升压时间随水蒸气体积分数的变化

    Figure  26.  Change of final overpressure and pressure rising time with water vapor volume fraction

    图  27  不同水蒸气体积分数下超压振幅随时间的变化

    Figure  27.  Variation of overpressure amplitude of different volume fractions of water vapor with time

    图  28  不同水蒸气体积分数下超压振荡频率随时间的变化

    Figure  28.  Variation of overpressure oscillation frequency of different volume fractions of water vapor with time

    图  29  密闭管道内CO2、N2和水蒸气对爆炸的抑制作用对比

    Figure  29.  Comparison of suppression of CO2, N2 and water vapor on the explosion in the closed tube

    图  30  尾端泄爆管道内CO2、N2和水蒸气对爆炸的抑制作用对比

    Figure  30.  Comparison of suppression of CO2, N2 and water vapor on the explosion in the end-vented tube

    表  1  添加气体的体积分数(密闭/泄爆)

    Table  1.   Volume fraction of added gas (closed/end-vented) % 

    CO2N2Water vapor (H2O) CO2N2Water vapor (H2O)
    000 202520
    454 243024
    8108 253126
    121512 27
    162016
    下载: 导出CSV
  • [1] 李建国, 郭昭胜, 张永生, 等. 室内燃气爆炸对混凝土结构的破坏及有限元分析 [J]. 消防科学与技术, 2018, 37(4): 563–566. doi: 10.3969/j.issn.1009-0029.2018.04.041

    LI J G, GUO Z S, ZHANG Y S, et al. The concrete structural failure caused by internal gas explosion and its finite element analysis [J]. Fire Science and Technology, 2018, 37(4): 563–566. doi: 10.3969/j.issn.1009-0029.2018.04.041
    [2] 闫秋实, 刘晶波, 伍俊. 典型地铁车站内爆炸致人员伤亡区域的预测研究 [J]. 工程力学, 2012, 29(2): 81–88.

    YAN Q S, LIU J B, WU J. Estimation of casualty areas in subway station subjected to terrorist bomb [J]. Engineering Mechanics, 2012, 29(2): 81–88.
    [3] 王波, 杜扬, 齐圣, 等. 油气爆炸在细长密闭管道内的振荡传播特性 [J]. 振动与冲击, 2017, 36(17): 97–103, 126.

    WANG B, DU Y, QI S, et al. Oscillation propagation characteristics of gasoline-air mixture explosion in elongated closed tubes [J]. Journal of Vibration and Shock, 2017, 36(17): 97–103, 126.
    [4] XIAO H H, MAKAROV D, SUN J H, et al. Experimental and numerical investigation of premixed flame propagation with distorted tulip shape in a closed duct [J]. Combustion and Flame, 2012, 159(4): 1523–1538.
    [5] LI H W, GUO J, YANG F Q, et al. Explosion venting of hydrogen-air mixtures from a duct to a vented vessel [J]. International Journal of Hydrogen Energy, 2018, 43(24): 11307–11313.
    [6] YAO Z F, DENG H X, ZHAO W L, et al. Experimental study on explosion characteristics of premixed syngas/air mixture with different ignition positions and opening ratios [J]. Fuel, 2020, 279: 118426.
    [7] LV P F, ZHANG J X, JU M H, et al. Influence of branch pipes on the deflagration characteristics of methane in confined space [J]. Process Safety Progress, 2020, 40(1): e12152.
    [8] 李国庆, 杜扬, 王波, 等. 点火位置对管道内油气泄压爆炸超压特性影响 [J]. 振动与冲击, 2017, 36(24): 204–212.

    LI G Q, DU Y, WANG B, et al. Effects of ignition position on overpressure characteristics of vented gasoline-air mixture explosion in a pipe [J]. Journal of Vibration and Shock, 2017, 36(24): 204–212.
    [9] PHYLAKTOU H, ANDREWS G E. Gas explosions in long closed vessels [J]. Combustion Science and Technology, 1991, 77(1/2/3): 27–39.
    [10] 郑立刚, 朱小超, 于水军, 等. 浓度和点火位置对氢气-空气预混气爆燃特性影响 [J]. 化工学报, 2019, 70(1): 408–416.

    ZHENG L G, ZHU X C, YU S J, et al. Effect of concentration and ignition position on characteristics of premixed hydrogen-air deflagration [J]. CIESC Journal, 2019, 70(1): 408–416.
    [11] SEARBY G. Acoustic instability in premixed flames [J]. Combustion Science and Technology, 1992, 81(4/5/6): 221–231.
    [12] XING H D, XU Q M, SONG X Z, et al. The effects of vent area and ignition position on pressure oscillations in a large L/D ratio duct [J]. Process Safety and Environmental Protection, 2020, 135: 166–170.
    [13] 朱传杰, 林柏泉, 江丙友, 等. 瓦斯爆炸在封闭管道内冲击振荡特征的数值模拟 [J]. 振动与冲击, 2012, 31(16): 8–12, 17. doi: 10.3969/j.issn.1000-3835.2012.16.002

    ZHU C J, LIN B Q, JIANG B Y, et al. Numerical simulation on oscillation and shock of gas explosion in a closed end pipe [J]. Journal of Vibration and Shock, 2012, 31(16): 8–12, 17. doi: 10.3969/j.issn.1000-3835.2012.16.002
    [14] ZHU C J, LIN B Q, JIANG B Y, et al. Numerical simulation of blast wave oscillation effects on a premixed methane/air explosion in closed-end ducts [J]. Journal of Loss Prevention in the Process Industries, 2013, 26(4): 851–861.
    [15] 韦世豪, 杜扬, 王世茂, 等. 不同形状受限空间内油气爆燃特性的实验研究 [J]. 中国安全生产科学技术, 2017, 13(5): 41–47.

    WEI S H, DU Y, WANG S M, et al. Experimental study on deflagration characteristics of gasoline-air mixture in confined space with different shapes [J]. Journal of Safety Science and Technology, 2017, 13(5): 41–47.
    [16] 韦世豪, 杜扬, 王世茂, 等. 非均匀容积式受限空间油气爆炸超压与火焰特征 [J]. 后勤工程学院学报, 2017, 33(4): 30–36. doi: 10.3969/j.issn.1672-7843.2017.04.006

    WEI S H, DU Y, WANG S M, et al. Overpressure and flame characteristics of gasoline-vapor explosion in the non-uniform volumetric confined space [J]. Journal of Logistical Engineering University, 2017, 33(4): 30–36. doi: 10.3969/j.issn.1672-7843.2017.04.006
    [17] 李毅. 管道中氢-空和甲烷-空预混火焰传播与压力振荡研究 [D]. 焦作: 河南理工大学, 2015.
    [18] 路长, 李毅, 潘荣锟. 管道氢气-空气预混气体爆炸特征的试验研究 [J]. 安全与环境学报, 2016, 16(3): 38–42.

    LU C, LI Y, PAN R K, et al. Experimental study on explosion tendency of hydrogen-air premixed gases in the duct [J]. Journal of Safety and Environment, 2016, 16(3): 38–42.
    [19] 王亚磊, 郑立刚, 于水军, 等. 约束端面对管内甲烷爆炸特性的影响 [J]. 爆炸与冲击, 2019, 39(9): 139–148.

    WANG Y L, ZHENG L G, YU S J, et al. Effect of vented end faces on characteristics of methane explosion in duct [J]. Explosion and Shock Waves, 2019, 39(9): 139–148.
    [20] WANG Z R, NI L, LIU X, et al. Effects of N2/CO2 on explosion characteristics of methane and air mixture [J]. Journal of Loss Prevention in the Process Industries, 2014, 31: 10–15. doi: 10.1016/j.jlp.2014.06.004
    [21] LI M, XIAO Y, DENG J, et al. Effect of CO2 on explosion limits of flammable gases in goafs [J]. Mining Science and Technology, 2010, 20(2): 193–197.
    [22] LI M H, XU J C, WANG C J, et al. Thermal and kinetics mechanism of explosion mitigation of methane-air mixture by N2/CO2 in a closed compartment [J]. Fuel, 2019, 255(1): 115747.
    [23] 张迎新, 吴强, 刘传海, 等. 惰性气体N2/CO2抑制瓦斯爆炸实验研究 [J]. 爆炸与冲击, 2017, 37(5): 906–912.

    ZHANG Y X, WU Q, LIU C H, et al. Experimental study on coal mine gas explosion suppression with inert gas N2/CO2 [J]. Explosion and Shock Waves, 2017, 37(5): 906–912.
    [24] 李凌飞. 甲烷爆炸特性及其抑爆技术研究[D]. 太原: 中北大学, 2012.
    [25] 任韶然. 惰性及特种可燃气体对甲烷爆炸特性的影响实验及分析 [J]. 天然气工业, 2013, 33(10): 110–115. doi: 10.3787/j.issn.1000-0976.2013.10.019

    REN S R. An experimental study of effects of insert and special flammable gases on methane’s explosion characteristics [J]. Natural Gas Industry, 2013, 33(10): 110–115. doi: 10.3787/j.issn.1000-0976.2013.10.019
    [26] 王华, 葛岭梅, 邓军. 惰性气体抑制矿井瓦斯爆炸的实验研究 [J]. 矿业安全与环保, 2008, 35(1): 4–7. doi: 10.3969/j.issn.1008-4495.2008.01.002

    WANG H, GE L M, DENG J. Experimental study of using inert gas to suppress mine gas explosion [J]. Mining Safety & Environmental Protection, 2008, 35(1): 4–7. doi: 10.3969/j.issn.1008-4495.2008.01.002
    [27] 李成兵. N2/CO2/H2O抑制甲烷爆炸化学动力学机理分析 [J]. 中国安全科学学报, 2010, 20(8): 88–92. doi: 10.3969/j.issn.1003-3033.2010.08.014

    LI C B. Chemical kinetics mechanism analysis of N2/CO2/H2O suppressing methane explosion [J]. China Safety Science Journal, 2010, 20(8): 88–92. doi: 10.3969/j.issn.1003-3033.2010.08.014
    [28] XIE Y L, WANG J H, ZHANG M, et al. Experimental and numerical study on laminar flame characteristics of methane oxy-fuel mixtures highly diluted with CO2 [J]. Energy and Fuels, 2013, 27: 6231–6237. doi: 10.1021/ef401220h
    [29] 秦文茜, 王喜世, 谷睿, 等. 超细水雾作用下瓦斯的爆炸压力及升压速率 [J]. 燃烧科学与技术, 2012, 18(1): 90–95.

    QIN W Q, WANG X S, GU R, et al. Methane explosion overpressure and overpressure rise rate with suppression by ultra-fine water mist [J]. Journal of Combustion Science and Technology, 2012, 18(1): 90–95.
    [30] 康泉胜, 李振明, 王睿, 等. 超细水雾对管内丙烷爆炸火焰抑制效果的试验研究 [J]. 安全与环境学报, 2015, 15(6): 111–114.

    KANG Q S, LI Z M, WANG R, et al. Experimental study on the inhibition effects of the ultra-fine water mist on the propane explosion flame in the tube [J]. Journal of Safety and Environment, 2015, 15(6): 111–114.
    [31] 高旭亮. 超细水雾抑制甲烷爆炸实验与数值模拟 [D]. 大连: 大连理工大学, 2014.
    [32] 刘丹, 司荣军, 李润之. 环境湿度对瓦斯爆炸特性的影响 [J]. 高压物理学报, 2015, 29(4): 307–312. doi: 10.11858/gywlxb.2015.04.011

    LIU D, SI R J, LI R Z. Ambient humidity on explosion characteristics of methane/air mixture [J]. Chinese Journal of High Pressure Physics, 2015, 29(4): 307–312. doi: 10.11858/gywlxb.2015.04.011
    [33] GEXCON A S. FLACS V10. 8 user’s manual [Z]. Bergen, Norway: GEXCON AS, 2017.
    [34] LAUNDER B E, SPALDING D B. The numerical computation of turbulent flows [J]. Computer Methods in Applied Mechanics and Engineering, 1974, 3(2): 269–289.
    [35] 陈晓坤, 丁园月, 程方明, 等. CO2对矿井多组分可燃性气体抑爆特性的影响 [J]. 煤炭科学技术, 2015, 43(3): 43–47.

    CHEN X K, DING Y Y, CHENG F M, et al. Influence of CO2 on explosion suppression characteristics of multicomponent flammable gas in coal mine [J]. Coal Science and Technology, 2015, 43(3): 43–47.
    [36] 罗振敏, 王涛, 程方明, 等. 小尺寸管道内二氧化碳抑制甲烷爆炸效果的实验及数值模拟 [J]. 爆炸与冲击, 2015, 35(3): 393–400. doi: 10.11883/1001-1455-(2015)03-0393-08

    LUO Z M, WANG T, CHENG F M, et al. Experimental and numerical studies on the suppression of methane explosion using CO2 in a mini vessel [J]. Explosion and Shock Waves, 2015, 35(3): 393–400. doi: 10.11883/1001-1455-(2015)03-0393-08
    [37] 张印, 赵东风, 刘义. 基于FLACS的CH4/CO2/air混合气爆炸参数分析 [J]. 中国安全生产科学技术, 2016, 12(9): 36–40.

    ZHANG Y, ZHAO D F, LIU Y. Analysis on explosion parameters of CH4/CO2/air mixed gas based on FLACS [J]. Journal of Safety Science and Technology, 2016, 12(9): 36–40.
    [38] LI Z, CHEN L, YAN H C, et al. Gas explosions of methane-air mixtures in a large-scale tube [J]. Fuel, 2021, 285: 119239.
    [39] 杨春丽, 刘艳, 胡玢, 等. 氮气和水蒸气对瓦斯爆炸基元反应的影响及抑爆机理分析 [J]. 高压物理学报, 2017, 31(3): 301–308. doi: 10.11858/gywlxb.2017.03.012

    YANG C L, LIU Y, HU F, et al. Effect of nitrogen an water vapor on methane-air mixture explosion elementary reaction suppression mechanism [J]. Chinese Journal of High Pressure Physics, 2017, 31(3): 301–308. doi: 10.11858/gywlxb.2017.03.012
    [40] 袁渭兰, 吕卫民, 贾忠湖. 气体动力学[M]. 北京: 科学出版社, 2013: 55.
    [41] 卓士创, 董慎行. 弹簧振子在定常干摩擦阻尼作用下的振动 [J]. 大学物理, 2001(6): 17–21, 33. doi: 10.3969/j.issn.1000-0712.2001.06.005

    ZHUO S C, DONG S X. The oscillation of a spring oscillator under constant dry frictional damping [J]. College Physics, 2001(6): 17–21, 33. doi: 10.3969/j.issn.1000-0712.2001.06.005
  • 加载中
图(30) / 表(1)
计量
  • 文章访问数:  2003
  • HTML全文浏览量:  1118
  • PDF下载量:  30
出版历程
  • 收稿日期:  2020-12-11
  • 修回日期:  2021-01-13

目录

/

返回文章
返回