2018S-92-多舱室内气云爆炸爆轰产物运动二维数值模拟
航空舱内烟气扩散数值仿真研究
航空舱内烟气扩散数值仿真研究随着时代的进步和人们对旅游的需求越来越高，航空业的发展也趋于快速。
然而，随着烟草的流行，烟草制品也被带上了飞机。
虽然现在已经有不吸烟的航班，但是仍有许多人坚持着在航班上吸烟。
烟草带来的不仅仅是健康问题，还可能会对航空运营安全造成影响。
如何将烟气扩散在航空舱内的问题变成了一个需要考虑的问题。
考虑到难以在实验室或飞机上实际进行试验，数值仿真成为了解决该问题的可行方法。
1.数值模拟的基本原理和方法数值模拟是通过计算机技术对模型进行数字化表示，并对其内在运动进行计算的一种方法。
在航空领域，对烟气扩散进行数值模拟，需要考虑气流传输的基本变量、吸附剂物理化学性质、气体的传递过程、热传递和质量传递等因素。
基于以上因素，数值模拟的具体方法按照复杂性程度的不同可以分为以下几类：1）CFD（Computational Fluid Dynamics）数值模拟方法。
该方法主要通过对流动流场的数值计算，来分析舱内的空气流动和烟气扩散特性。
在CFD中，将舱内空气流动视为复杂的流动场，对其进行离散化后通过求解Navier-Stokes方程得出舱内的空气流动特性。
2）统计模型方法。
该方法是基于对化学反应和热力学性质等方面的实验和观测结果进行统计模型的建立。
采用该方法可以在较高的精度下对航空舱内的烟气扩散情况进行模拟，同时也可以获得化学反应的相关数据。
3）热传递和传质模型。
该方法是基于热传递和传质现象的数学模型来解决烟气扩散问题。
通过使用数学方法对烟气扩散的物理现象进行数值求解，来克服烟气扩散实验的困难性。
该方法的优点主要在于可以对烟草烟雾扩散机理的深入理解。
2.航空舱内的烟气扩散源烟气扩散源是航空舱内烟气扩散问题的核心。
当实验中的气流不强或舱内气流受环境因素影响而不稳定时，烟气扩散源也会影响到实验结果。
要保证烟气扩散模拟有效，应对烟气扩散源进行实验测试、定量分析，同时对于模拟不同环境及使用情况下的大量测试需求时，可以采用一般合成方法，把如温度、湿度、自然对流等因素考虑进去并制成程序，可灵活选择模拟环境参数。
大加速度场中二维层流燃过程的数值模拟
第 ８ （０２第 ３ 卷 ２０ ） 期 Ｊ ＵＲ Ｌ ＯＦＣ Ｏ ＮＡ 燃
烧
科
学
与
技
术
ＴＥＣ ＨＮ Ｏ ＬＯ Ｇ Ｙ Ｏ Ｍ ＢＵ Ｓ Ｉ Ｎ Ｓ ＥＮＣＥ Ａ Ｎ Ｄ Ｔ Ｏ ＣＩ
Ｖ 。 ．８（ Ｏ２ Ｎ 。．３ １ ２０ ）
大 加 速 度 场 中二 维 层 流 燃 烧 过 程 的数 值 模 拟
方 向 与 对 称 轴 方 向一 致 时 ， 速 度 场 强 度 和 进 口流 速 等 因素 对 燃 烧 过 程 的 速 度 场 、 度 场 及 燃 烧 稳 定 性 的影 响 。 加 温 模
拟 燃 料 为 Ｃ Ｈ。 空 气 ， 烧 室 壁 面 绝 热 。 算 结 果 表 明 ， 大 加 速 度 场 中 ， 焰 面 附 近 的 温 度 梯 度 所 产 生 的 “ 力 ” 。 和 燃 计 在 火 浮
Ｗ Ｕ ｉ — ｉｎ Ｊｎ ｘａ ｇ ，ＬＩ Ｌｉｎ ｓ ｅ ｇ ，Ｗ ＡＮ Ｇ － ｕ ，ＧＵＯ — ｉ ｇ Ｕ ａ －ｈ ｎ 。 Ｅｎ ｙ Ｙａ ｐｎ
（ ．Ｄｅ ｒ ｍ ｅ １ ｐａ ｔ ｎｔｏｆＴｈｅｒ ａｌＳｃ ｅ ｅ ａ ｄ ｎ ｒ ｍ ｉ ｎｃ ｎ Ｅ ｅ ｇｙ Ｅｎｇｉ ｅ ｉ ｇ 。 ｅ ｎｉ ｅ ｓｔ ｃ ｎｅ ｒ ｎ Ｈ ｂｅｉＵ ｖ ｒ ｉｙ ｏｆＴｅ ｈｎｏｌ ｏｇｙ，
吴 晋 湘 ，刘 联 胜。 ，王 恩 宇 ，郭 雅 萍
（ ．河 北 工 业 大 学 热 能 动 力 系 ， 津 ３ ０ ３ ； ．天 津 大 学 内 燃 机 燃 烧 学 国家 重 点 实 验 室 ， 津 ３ ０ ７ ） １ 天 ０１２ ２ 天 ０ ０ ２
摘
要 ： Ｓ ＭＰ Ｅ 方 法 ， 大 加 速 度 场 中 的 二 维 轴 对 称 预 混 层 流 燃 烧 过 程 进 行 了 数 值 模 拟 ， 析 了 当 加 速 度 场 用 Ｉ Ｌ Ｒ 对 分
PBX炸药二维冲击起爆机理的数值模拟
PBX炸药二维冲击起爆机理的数值模拟
梁增友;黄风雷;张振宇;王辉
【期刊名称】《火炸药学报》
【年(卷),期】2008(031)005
【摘要】建立了基于KIM模型的一个化学反应速率方程模型,用于分析PBX起爆和爆轰行为的损伤作用.运用遗传算法确定了反应速率方程的相关参数,通过与Forest-Fire反应速率模型数值模拟结果对比,验证了所建模型的合理性.将所建反应速率方程模型嵌入有限元程序对PBX炸药二维冲击起爆过程进行数值模拟,数值模拟结果与试验结果吻合较好,可以描述分析受冲击加载造成孔隙率、颗粒尺寸等变化的损伤炸药的冲击起爆过程.
【总页数】4页(P15-18)
【作者】梁增友;黄风雷;张振宇;王辉
【作者单位】中北大学机电工程学院,山西,太原,0300051;北京理工大学爆炸科学与技术国家重点实验室,北京,100081;国防科技大学理学院,湖南,长沙,410073;中北大学机电工程学院,山西,太原,0300051
【正文语种】中文
【中图分类】TJ55;TQ564
【相关文献】
1.PBX炸药细观结构冲击点火的二维数值模拟 [J], 刘群;陈朗;伍俊英;王晨
2.含铝炸药二维冲击起爆的爆轰数值模拟 [J], 丁刚毅;徐更光
3.PBX-9404炸药冲击起爆的一维数值模拟 [J], 曾代朋;张振宇
4.聚能射流冲击起爆屏蔽压装PBX炸药的试验研究 [J], 宋乙丹;陈科全;路中华;陈翔;陈红霞;李兴隆
5.孔隙度对PBX炸药冲击起爆影响的实验研究 [J], 刘海庆; 段卓平; 白志玲; 温丽晶; 欧卓成; 黄风雷
因版权原因，仅展示原文概要，查看原文内容请购买。
储罐内可燃气云爆炸压力数值模拟对比分析研究
全性问题，同时可 将 理 论 解 和 少 量 试 验 结 果 作 为 数值计算结果的验证 。 ［４～６］
常用确定容器内部爆炸载荷数值模拟方法主 要有两种：一种是 将 气 云 爆 炸 的 破 坏 作 用 等 效 为 ＴＮＴ爆炸 的 破 坏 作 用，即 为 ＴＮＴ当 量 法，通 过 模 拟炸药的 爆 炸 压 力 近 似 得 到 可 燃 气 云 的 爆 炸 载 荷；还有一种是基 于 流 体 动 力 学 和 化 学 反 应 动 力 学理 论，实 现 爆 炸 流 场 各 物 理 量 的 模 拟。 文 献 ［７］即 采 用 ＴＮＴ当 量 法，对 １０００ｍ３ 立 式 钢 储 罐 进行了数值模拟，探 讨 了 不 同 因 素 （浓 度、高 径 比 及流固耦 合 等 ）对 爆 炸 载 荷 的 影 响，分 析 了 爆 炸 冲击载荷下储罐的应力。 而 文 献 ［８］则 采 用 流 体 动力 学 （ＣＦＤ）方 法，基 于 Ｆｌｕｅｎｔ软 件 分 析 了 连 通 容器 内 气 体 爆 炸 流 场 中 温 度、压 力、速 度、密 度 和 燃烧速率随时间 的 变 化 规 律，较 为 清 晰 地 反 映 出 气体爆 炸 过 程。 ＴＮＴ当 量 法 和 计 算 流 体 动 力 学 （ＣＦＤ）方法都能 够 比 较 准 确 地 计 算 爆 炸 压 力，但
立式拱顶储罐在我国的大型油库和石油化工 企业得到了广泛 应 用，其 内 部 储 存 油 品 挥 发 出 的 可燃气体与空气 混 合 后，遇 点 火 源 极 易 引 发 爆 炸 事故 。 ［１～３］ 在我国 石 油 石 化 介 质 储 存 中，储 罐 已 逐 步 实 现 大 型 化 ，储 罐 容 量 越 大 ，意 味 着 单 个 储 罐 发生燃烧爆炸事 故 引 起 的 损 失 越 大，因 此 需 要 对 储罐安全性予以更多的重视。合理确定储罐内部 的燃爆载荷及其 变 化 规 律，对 于 减 少 爆 炸 事 故 的 发 生 、储 罐 的 设 计 、研 发 和 后 期 维 护 具 有 重 要 的 指 导意义。试验无疑是科学研究获取相关参数最直 接 有 效 的 方 式 ，但 由 于 爆 炸 具 有 很 强 的 破 坏 作 用 ， 如 果 完 全 通 过 试 验 来 获 取 数 据 ，花 费 巨 大 ，此 外 还 存在不可预测的 安 全 性 问 题。采 用 理 论 模 型 （经 验公式）可 以 直 接 获 得 最 大 超 压，但 过 度 的 简 化 往往使计算结果产生较大的偏差。数值计算为爆 炸问题的研究提 供 了 一 种 新 思 路，能 够 获 得 更 加 全 面 、足 够 精 确 的 计 算 结 果 ，花 费 较 小 又 不 涉 及 安
多舱室内气云爆炸爆轰产物运动二维数值模拟
多舱室内气云爆炸爆轰产物运动二维数值模拟
徐维铮;吴卫国
【期刊名称】《中国舰船研究》
【年(卷),期】2018(013)0z1
【摘要】[目的]为了研究多舱室空间内气云爆炸后爆轰产物的运动规律,[方法]采用五阶WENO有限差分格式,通过耦合求解爆轰产物质量分数输运方程与可压缩欧拉方程,自主编写多舱室内气云爆炸过程二维数值计算程序.利用平面激波绕射矩形台阶试验验证本文数值程序的可靠性与正确性.采用验证的程序开展不同爆炸工况气云爆炸过程数值计算,探讨多舱室内爆炸过程中泄压、后燃烧能量加入对爆轰产物运动过程的影响规律.[结果]结果表明:增大泄压口尺寸,可加剧爆轰产物泄出过程;增大反应速率常数,可加速爆轰产物向相邻舱室的扩张过程.[结论]研究成果可为进一步考虑多组分复杂燃烧效应后爆轰产物扩张过程及毁伤评估提供一定的参考和指导.【总页数】8页(P92-99)
【作者】徐维铮;吴卫国
【作者单位】武汉理工大学高性能舰船技术教育部重点实验室,湖北武汉430063;武汉理工大学交通学院,湖北武汉430063;武汉理工大学高性能舰船技术教育部重点实验室,湖北武汉430063;武汉理工大学交通学院,湖北武汉430063
【正文语种】中文
【中图分类】U674.7
【相关文献】
1.爆轰产物驱动飞片运动数值模拟 [J], 金柯;李平;吴强;金孝刚
2.进气总压对连续旋转爆轰发动机爆轰影响的二维数值模拟 [J], 李宝星;翁春生
3.爆轰产物驱动飞片运动数值模拟研究 [J], 金柯;李平;吴强;金孝刚
4.密闭多舱室内柱形气云爆炸过程二维数值模拟 [J], 徐维铮;吴卫国;况正
5.爆轰产物物态方程(Ⅰ)——爆热、爆轰产物的等熵方程和物态方程 [J], 李银成因版权原因，仅展示原文概要，查看原文内容请购买。
空间爆炸冲击波的数值模拟
空间爆炸冲击波的数值模拟一：无限空间爆炸如图所示，半径为7.0cm的圆柱形TNT装药，质量为5.018Kg从炸药中心单点起爆后在无限空间中传播。
试分析起爆后冲击波的传播及压力分布特性。
二：建模分析1材料模型及参数设置本数值模拟采用的基本材料为TNT炸药，空气。
在无限空气领域中传播。
1）空气空气简化为无粘性理想气体，冲击波的膨胀假设为等熵绝热过程以LS_DYNA中的*MAT_NULL材料模型和线性多项式状态方程*EOS_LINEAR_POLYNOMIAL方程来描述。
*EOS_LINEAR_POLYNOMIAL方程具体表达式为：式中C0、C1、C2、C3、C4、C5、C6是与气体性质有关的常数，C0 = C1 = C2 = C3 = C4 = C6 = 0, ;，、e0及分别为气体的初始密度、密度、初始单位体积内能和绝热指数。
空气材料的模型参数取值：=1.292910-3g/cm3, e0=2.5×105Pa,=1.4。
表 1 空气状态方程参数变量EOSID C0 C1 C2 C3 C4 C5 C6数值 1 -1.0E-6 0.0 0.0 0.0 0.4 0.4 0.02) 炸药以LS_DYNA 中的*MAT_HIGH_EXPLOSIVE_BURN材料模型和*EOS_JWL方程模拟TNT炸药。
*EOS_JWL方程的表达式为：式中，P为压力，V为相对体积，即爆轰产物体积与炸药初始体积之比；E0为炸药的初始比内能，即单位体积内能。
A、B R1、和R2是与炸药性质有关的常数。
炸药模型的各参数取值：密度g/cm3, 爆速D=0.693cm/, 压力PCJ=0.27105MPa，A=3.74105MPa, B=0.0733105MPa,R1=4.15,R2=0.95,=0.3,0=0.07105MPa。
如下表：表2 炸药材料参数变量MID RO D PCJ BETA K G SIGY 数值 1.O 1.63 0.693 0.27 0.0 0.0 0.0 0.0表3 JWL方程参数变量EOSID A B R1 R2 OMEG E0 V数值 1 3.74 0.073 4.15 0.95 0.3 0.07 1.0由于LS_DYNA在爆炸分析中用的基本单位为-ｇ-的单位系统，故表1-3中各参数取值由m-㎏-s单位换算得到。
复杂多舱室内炸药爆炸爆轰产物运动二维数值模拟研究
复杂多舱室内炸药爆炸爆轰产物运动二维数值模拟研究
徐维铮;吴卫国
【期刊名称】《应用力学学报》
【年(卷),期】2019(0)2
【摘要】为了研究复杂多舱室空间内炸药爆炸后爆轰产物运动规律,采用三阶WENO有限差分格式,通过耦合求解两组分可压缩欧拉方程,自主编写了多舱室内炸药爆炸过程二维数值计算程序。
利用平面激波绕射矩形台阶上的竖直平板实验与数值模拟纹影图进行了对比,图中不同时刻的波系结构和台阶、竖直平板处绕射形成的漩涡结构均证明了文中程序开发的可靠性和正确性。
采用验证的程序开展了复杂多舱室内炸药爆炸过程数值计算,探讨了复杂多舱室内炸药爆炸爆轰产物的传播路径及泄压口泄压对爆轰产物运动规律的影响。
研究表明:舱室泄压对复杂多舱室内爆轰产物运动的影响,在爆炸初期较小,在爆炸后期较大;爆炸舱室内部泄压能显著地降低相邻舱室超压载荷强度(40%左右),减弱爆炸对相邻舱室的毁伤效应。
本文的研究可为进一步探讨、考虑燃烧效应后爆轰产物扩张过程及爆炸毁伤评估提供一定的参考和指导。
【总页数】6页(P400-404)
【作者】徐维铮;吴卫国
【作者单位】武汉理工大学高性能舰船技术教育部重点实验室;武汉理工大学交通学院船舶
【正文语种】中文
【中图分类】U674.70;TQ560.1
【相关文献】
1.高能炸药非理想爆轰的二维数值模拟
2.高能炸药非理想爆轰的二维数值模拟
3.多舱室内气云爆炸爆轰产物运动二维数值模拟
4.内爆炸爆轰产物运动过程二维数值模拟研究(二)——泄压空间
5.内爆炸爆轰产物运动过程二维数值模拟研究(一)——密闭空间
因版权原因，仅展示原文概要，查看原文内容请购买。
开敞空间可燃气云爆燃场数值模拟
开敞空间可燃气云爆燃场数值模拟
罗正鸿;詹晓力;丁信伟;王淑兰;阳永荣
【期刊名称】《化工学报》
【年(卷),期】2003(054)003
【摘要】通过在Euler能量守恒方程中加入一表示能量释放速率物理意义的"开关函数"得到了Euler坐标系中描述有限能量释放速率的二维可燃气云爆燃场的数学模型,给出了初值及边界条件的取法,并做了相关的模拟算例.模拟计算得出的爆燃时间-压力曲线与实验结果吻合很好.
【总页数】5页(P383-387)
【作者】罗正鸿;詹晓力;丁信伟;王淑兰;阳永荣
【作者单位】浙江大学化学工程与生物工程学系,浙江,杭州,310027;浙江大学化学工程与生物工程学系,浙江,杭州,310027;大连理工大学化工学院,辽宁,大连,116012;大连理工大学化工学院,辽宁,大连,116012;浙江大学化学工程与生物工程学系,浙江,杭州,310027
【正文语种】中文
【中图分类】TQ016;X928.7
【相关文献】
1.小气量开敞空间可燃气云爆燃实验 [J], 罗正鸿;詹晓力;丁信伟;王淑兰;阳永荣
2.基于小扰动波理论的开敞空间可燃气云爆燃研究 [J], 罗正鸿;詹晓力;丁信伟;王淑兰;阳永荣
3.可燃气体云爆燃场数值模拟 [J], 丁信伟;李志义
4.开敞空间可燃气云爆炸数值模拟研究 [J], 党福辉;董呈杰;孙旭红
5.有障碍物开敞空间可燃气云爆炸超压场的数值模拟 [J], 陈文瑛;柴建设
因版权原因，仅展示原文概要，查看原文内容请购买。
核爆炸冲击波与爆炸产物相互作用的二维数值模拟
在核武器的爆炸过程 中， 核爆炸主冲击波
展研究 ， 对理解爆炸产物 的上升和发展规律具
的背面将形成由稀疏波和次级 冲击波组成的波
系结构［ ， １ 爆炸气体产物的压力 、 ］ 密度和温度等 参数在该波系结构的作用下将发生改变［ ； ］另
界面 出现不规则和涡结构 ， 即爆 炸产物呈现非同一性； 另一方面 ， 地面反射 冲击波将穿越 诱导波系 ， 挤 并
压爆炸产物的下表 面， 物的形状 和参 数发生改变 ， 而加速 爆炸产物的浮力上 升过程并增大其最 大 使产 从 上升高度 ， 继而扩大了放射 性烟云的扩散范围。 关键词 ： 核爆炸放射性烟云 ； 核爆炸 冲击波 ； Ｎ ＷＥ Ｏ格式 ； 产物 ｝ ； 爆炸 火球 真实气体效应
汪 林 ， 毅 陈 琳 程先友 郑 ， ，
（． １防化研究院 ， 京 １２０ ；． 北 ０２５２海军核化安全研究所 ， 北京 １０７ ） ００７
摘要 ： 核爆炸 冲击波与爆炸产物的相互作用 过程是影响放射性烟云扩散 程度 的重 要因素 。为 了阐
明二者相互作用的物理图景 ， 采用 Ｔ Ｄ保持性 质的 Ｒ ｎｅ ｕｔ离散 和 ＷＥ Ｏ格式 对核爆炸 冲击波 Ｖ ｕ ｇ－ ｔ Ｋ ａ Ｎ 在真实气体中的传播 过程进行 了计算 ， 并分析 了其诱 导波 系和反 射波系 与爆炸 产物 的相互 作用规律 。 计算结果表明 ， 在诱导波系的压缩和稀疏作用下 ， 爆炸产 物的压力 和密度值 呈高低 相间分 布 ， 物的下 产
Ｑ （一 Ｅ） （一 Ｆ Ｓ１ ＋ Ｅ 壶 ＋Ｆ 责 一 （ ） ）
不同大气环境下爆炸烟云运动的二维数值计算
边 界条 件 ： 矩 形空 间上 下 边界 为无 滑移 边界 ； 矩形 空 间两 侧边 界 为周 期性 边界 条件 。 计算 结 果如 图 ２所示 。
１２ 计 算模 型 的验 证 ．
用 Ｒ ｙｅ ｈＴ ｙｏ 不稳 定 问 题 进 行 数值 验 ａｌｉ —ａｌ ｇ ｒ
证， 并将 结果 与其 它 公开 发表 的结 果进 行 比较 ，
程 可写 为 ：
理试验Ｌ 】 剖和数值计算 两种途径 ， 由于受试
验条件 、 适用 范 围 、 重复 性 等条 件 的限 制 ， 可 对
烟 云的研究趋 向于数 值计算 方 面。前人 的研究
重 点 在 于均 一 大 气 条 件 烟 云 的运 动或 者 采 用 Ｂ ｕｓ ｅｑ近似模 型 … 进行 讨 论 。均 一 大气 条 ｏ ｓ ｎｓ ｉ 件与真 实 大气 的 分 层 特 征 不 符 ， Ｂ ｕｓ ｅｑ 而 ｏ ｓｉ ｓ ｎ 近似模型较适 用 于浅 对 流运 动 Ｊ 与核 爆 烟 云 ， 的深对 流 运 动 特 征 有 一定 差 异 。本 文 采 用 不可压 缩模 型从数 值计 算 的角 度对 烟云在 标准 大气 和均 一大气 下 的运动 进行 了讨论 ， 出两 指 者 的异 同点 ， 为研 究 核爆 携 带放 射 性粒 子 的烟 云在真 实大气 中的运 动提供基 础 。
向量 之 和 ， ： ７ ＋ 即 ： ７
７ ０。 ７
合 项 目（０ ７ １０ 和 （０ ７ Ａ ６ ４ ８ 资 助 。 １６ ６ ２ ） ２０ Ａ ０ Ａ ０ ） 作 者简 介 ： 晓 丽 （ ９０ ） 女 ． 东人 ， 士 研 究 生 ， 李 １８ 一 ， 山 博
研究方向 ： 核爆 效应 现 象 学 。
云的运动过程 进行 了数值计算研究 ， 并通过元粘不可压缩流体二维 Ｒ ｙ ｉ —ａｌｒ ａｌｇ Ｔｙｏ 不稳定性问题对程序 ｅｈ 进行 了数值验证 ； 通过 比较烟 云在均一和标准分层大气条件下 的运动 ， 得出大气 分层对 爆炸烟 云运动 的 影响主要体现在烟云上升高度上 的结论 ， 为进一步研究爆炸烟云 的运动提供了参考 。
舰艇防火分区内火灾蔓延数值模拟
舰艇防火分区内火灾蔓延数值模拟张光辉1 , 浦金云1 , 营1 , 博2李曹(1. 海军工程大学船舶与动力学院, 湖北武汉430033; 2. 海军装备部武汉局, 湖北武汉430033)摘要: 为了建立多舱室间烟气蔓延模型, 计算舰艇防火分区内各舱室的火灾烟气温度、高度随时间的变化, 并用计算结果指导舰艇防火决策。
首先假设防火分区内起火舱室的烟气沉降符合双层区域模型, 而烟气经过通风口蔓延到相邻舱室时没有卷吸现象发生; 然后运用质量、能量守恒和气体状态方程建立防火分区内多舱室烟气蔓延模型。
应用该模型对某舰艇防火分区的按照t2 发展的火灾场景进行计算, 计算结果比美国的C FA ST 软件计算的结果略高出10% , 主要原因是C FA ST 中考虑了下层冷空气温度变化。
关键词: 防火分区; 双层区域模型; 火灾蔓延中图分类号:U 664 文献标识码:文章编号: 100923443 (2009) 增刊20027206AC om p u te r s i m u la t i o n fo r f i re s p re a d ing in v e s s e l f i re p ro o f i ng s ub a re aZ H A N G G u a n g2h u i1 ,P U J in2y u n1 , L I Y i n g1 ,CA O B o2(1. Sh i p an d Pow e r Co llege, N ava l U n ive r s ity o f E n g i n ee r i n g, W uh a n 430033, C h ina;2. W uh a n B u reau o f N avy E q u i p m en t D ep a r t m en t, W uh a n 430033, C h ina)A b s t r a c t: T o se t up sm o k e sp read i n g m o d e l fo r m u lt i2cab i n i n ve s se l f i rep roo f i n g su b a r ea, sm o k e tem p e ra2 tu re an d h e i gh t t ren d w ith t i m e i n each cab i n w e re com p u tec l, an d th e re su lt s w e re u sed i n f ire f i g h t i n g de2 c is i o n 2m ak i n g. F ir s t ly it w a s suppo sed th a t sm o k e su b s i d i n g i n f ire cab i n acco rd s w ith th e tw o2l aye r z o n e m o de l an d sm o k e sp read i n g th ro u gh th e op en w ill no t w rap a ir. T h en th e law o f m a s s co n ve r s a t i o n,th e law o f en e rgy co n ve r s a t i o n an d th e ga s s ta te equ a t i o n w e re refe r red to e s tab lish th e m o de l i n th e ve s se l f i re2 p roo f su b a r ea. T h e re s u lt s b y C F A ST so f t w a r e w e re 10 p e rcen t l ow e r th a n tho se com p u t ed b y th e m o d e l.T h e rea s o n is th a t C FA ST so f t w a re ta k e s th e ch a n g i n g tem p e ra t u r e i n l ow e r laye r i n t o co n s i de r a t i o n.Ke y w o rd s: f irep roo f i n g su b a r ea; tw o2laye r z o n e m o d e l; f ire sp read i n g目前, 火灾模拟主要包括场模拟和区域模拟, 在区域模拟中, 双层区域模拟应用最广泛, 其思想是将起火房间分为上下两个区域: 上层烟气和下层空气, 如图1 所示, 并假定每个区域内各物理参数均匀一致, 而区域与区域之间的质量、能量没有直接交换, 它们的交换是通过羽流来实现的1 。
舱室火灾轰燃现象重构实验与数值模拟
技 术 的发 展 ， 利 用 适 当 的计 算 机 软 件 对 火 灾 事 故 进 行 重 构， 既 可 节 约 大 量 资金 和 宝 贵 时 间 ， 又 可 以 得 到 调 查 所 需 的火 灾 发 展 中各 种 动 力 学 特 征 参 量 ， 是 目前 火 灾 事 故 调
查 的重 要 辅 助 手段 。
垂 誊 消 防勰 论 研 究囊
舱 室 火 灾 轰 燃 现 象 重 构 实 验 与 数 值 模 拟
陈 晓洪 。 杨 枫。 李 其 修 ，吴 向君 （ 海 军工程 大学 动 力工程 学 院 ， 湖 北 武汉 ４ ３ ０ ０ ３ ３ ）
摘 要： 为 了研 究 船 舶 封 闭舱 室 火 灾 轰 燃 现 象 ， 首 先 用 非
“¨ ¨ － ・¨¨ ¨ 。 ｌ ＿ ¨ ｝ 。 ｌ ｌ ¨ ｌ ¨。 ‘ ｌ ｌ ¨ ＿ － ・ ¨ ¨ Ｉ 一。 Ｉ ｉ ｎ ， ， ． ， ｏ ｌ ＿ ｌ ¨ 。 ｌ ｌ ¨ ｌ ・。 。 …ｉ ｎ ｉ ． ， ＇ ｌ ｌ Ｉ ｌ ¨ ’ ” Ｉ ｉ ｉ ｏ ， ｏ ， ｑ 。 ｌ ｌ ¨ ¨。 Ｉ “ 。 －‘ Ｉ ‘ － 一 … ｌ 。 ｌ ¨ｌ ・ ” ０ ¨ ＂ ‘ ‘ ’ ¨ ¨ Ｉ 一。 ｌ ｌ ¨ 。 Ｉ 。 －
中图 分 类 号 ： Ｘ９ １ ３ ． ４ 。 Ｕ６ ７ ４ ． ７
文献标志码 ： Ａ
文章 编 号 ： １ ０ ０ ９ —０ ０ ２ ９ （ ２ ０ １ ３ ） ０ ４ —０ ３ ５ ４ —０ ４
风控 制 、 火 灾 由发 展 期 进 入最 盛 期 的 标 志 ， 因此 ， 针 对 轰 燃 进行 火 灾 重 构 可 为火 灾 调 查 提供 重 要 的判 断依 据 。 １ 轰 燃 发 生 的 临 近 判 据 由非 线 性 动 力 学 理 论 ， 认 为 上 层 烟 气 温 度 Ｔ 随 时 间
密闭多舱室内柱形气云爆炸过程二维数值模拟
密闭多舱室内柱形气云爆炸过程二维数值模拟徐维铮;吴卫国;况正【摘要】为了研究密闭多舱室空间内爆炸波传播过程及舱室壁面爆炸载荷分布特性,采用三阶WENO有限差分格式,编写了多舱室内爆炸冲击波二维数值计算程序.利用平面激波与矩形方腔相互作用实验验证了本文数值程序的可靠性与正确性.采用验证的程序开展了密闭多舱室内柱形高压气云爆炸波数值计算研究,探讨了多舱室内爆炸波传播路径,不同舱室内部爆炸载荷特性及爆源位置对爆炸载荷特性的影响规律.本文研究可为工程抗爆结构设计及爆炸毁伤评估提供一定的参考和指导.【期刊名称】《舰船科学技术》【年(卷),期】2018(040)011【总页数】9页(P21-29)【关键词】多舱室内爆炸;爆炸波;三阶WENO格式;数值计算;冲量【作者】徐维铮;吴卫国;况正【作者单位】武汉理工大学高性能舰船技术教育部重点实验室,湖北武汉 430063;武汉理工大学交通学院船舶、海洋与结构工程系,湖北武汉 430063;武汉理工大学高性能舰船技术教育部重点实验室,湖北武汉430063;武汉理工大学交通学院船舶、海洋与结构工程系,湖北武汉 430063;武汉理工大学高性能舰船技术教育部重点实验室,湖北武汉 430063;武汉理工大学交通学院船舶、海洋与结构工程系,湖北武汉430063【正文语种】中文【中图分类】O354.50 引言工业可燃性气云爆炸事故大多是由弱点火点燃可燃气体引起的，其传播形式大多是以亚音速传播的爆燃波。
爆燃波的传播过程很复杂，受环境条件和物理因素的影响极大，当气体在传播过程中由于障碍物阻碍发生湍流加强效应或是在初始时刻由高能量直接起爆，其爆炸模式将由爆燃转为爆轰或直接以爆轰模式发生爆炸。
气体发生爆轰以后将在空气中形成强度较大的带有负压区的空气冲击波，对工作人员和周边结构设施将造成较大的损害。
针对空中可燃气云、云雾爆炸模拟，近年来徐胜利等学者进行过相应的数值研究[1 - 4]，他们将可燃气云、云雾等效的简化为一团高压气体，而舱室内部可燃气云爆炸研究较少。
爆震燃烧波在管内传播过程的二维数值模拟
爆震燃烧波在管内传播过程的二维数值模拟
王丁喜;严传俊
【期刊名称】《机械科学与技术》
【年(卷),期】2006(025)004
【摘要】采用三阶TVD迎风格式和Strang-splitting算子分裂法,对氢氧混合物爆震燃烧进行了二维数值计算;获得了爆震燃烧流场的清晰、独特和规则的三波结构;分析了三波结构的构成和成因以及爆震燃烧直接起爆后由强爆震燃烧波衰减为稳定自持的爆震燃烧波的过程.
【总页数】5页(P461-464,475)
【作者】王丁喜;严传俊
【作者单位】西北工业大学,西安,710072;西北工业大学,西安,710072
【正文语种】中文
【中图分类】O38
【相关文献】
1.基于WENO格式的爆震波传播过程一维数值模拟 [J], 陈灯;何立明;闫莉;罗俊
2.利用详细燃烧模型对裂解炉二维模型富氧燃烧过程进行数值模拟 [J], 申东发;王国清;刘俊杰;张利军;周先锋
3.二维会聚波在气体内传播不稳定性的数值模拟 [J], 刘金宏;谭多望;张旭;邹立勇;黄文斌
4.爆震波在爆震燃烧室内传播过程的数值模拟 [J], 何立明;严传俊;范玮;王建峰
5.爆震波点火器管路传播过程数值模拟 [J], 郭红杰;梁国柱
因版权原因，仅展示原文概要，查看原文内容请购买。
PETN、RDX和HMX炸药爆轰参数的数值模拟
PETN、RDX和HMX炸药爆轰参数的数值模拟李德华;程新路;杨向东;吴国栋【期刊名称】《爆炸与冲击》【年(卷),期】2005(025)004【摘要】用吉布斯自由能最小原理,通过解化学平衡方程组,求解PETN、RDX和HMX炸药爆轰产物系统的平衡组分,计算结果与用BKW和LJD方法计算的结果相近.用自编的程序从碳的石墨相、金刚石相、类石墨液相和类金刚石液相4种相态中确定出炸药爆轰产物中游离碳更可能存在的相态,并用此相态计算碳的吉布斯自由能.以WCA状态方程作为爆轰气相产物的物态方程,对PETN、RDX和HMX炸药爆轰参数作了预言,爆轰CJ点的爆速、爆压和爆温的计算结果与实验值吻合得很好.【总页数】5页(P325-329)【作者】李德华;程新路;杨向东;吴国栋【作者单位】四川大学原子分子物理研究所,四川,成都,610065;四川大学原子分子物理研究所,四川,成都,610065;四川大学原子分子物理研究所,四川,成都,610065;中国工程物理研究院流体物理研究所冲击波物理与爆轰物理实验室,四川,绵阳,621900;中国工程物理研究院流体物理研究所冲击波物理与爆轰物理实验室,四川,绵阳,621900【正文语种】中文【中图分类】O381【相关文献】1.某些炸药在CI源中的行为——TNT、RDX、PETN、HMX和特屈儿的正负CI 谱 [J], 徐建中2.RDX/HMX颗粒炸药落锤撞击点火-燃烧机理 [J], 吴艳青;鲍小伟;王明扬;黄风雷;张柱3.α-AlH3对HMX基炸药爆轰参数的影响 [J], 牛磊;王亲会;曹少庭;牛国涛;金大勇;王淑萍4.PETN炸药爆轰参数的数值模拟 [J], 李德华;程新路;杨向东;吴国栋5.糖颗粒对HMX和RDX两种单质炸药非冲击点火的影响 [J], 张钊; 吴艳青因版权原因，仅展示原文概要，查看原文内容请购买。
燃料-空气云雾爆轰的直接引爆实验研究
燃料-空气云雾爆轰的直接引爆实验研究
解立峰;郭学永;果宏;徐晓峰
【期刊名称】《爆炸与冲击》
【年(卷),期】2003(023)001
【摘要】在直径240mm的立式激波管中环氧丙烷(PO)、正己烷、癸烷分别与空气混合进行了直接起爆,测定了不同燃料、不同当量比的云雾直接引爆的临界起爆能.实验发现,当PO-空气混合物、正己烷-空气混合物及癸烷-空气混合物的当量比分别为1.05、1.12及1.15时,其临界起爆能均最小;三种混合物的可爆下限当量比分别为0.47(质量浓度4.23%)、0.75(质量浓度5.29%)及0.89(质量浓度5.8%).实验表明,在本实验条件下,三种燃料与空气混合物的云雾都容易在较大当量比范围内引发爆轰并实现爆轰波稳定传播.
【总页数】3页(P78-80)
【作者】解立峰;郭学永;果宏;徐晓峰
【作者单位】南京理工大学化工学院,江苏,南京,210094;南京理工大学化工学院,江苏,南京,210094;南京理工大学化工学院,江苏,南京,210094;南京理工大学化工学院,江苏,南京,210094
【正文语种】中文
【中图分类】O384
【相关文献】
1.液体碳氢燃料云雾爆轰特性的实验研究 [J], 姚干兵;解立峰;刘家骢
2.碳氢燃料云雾爆轰特性的实验研究 [J], 徐晓峰;解立峰;彭金华;陈瑛;苏华
3.三种燃料-空气混合物爆轰的直接引爆实验研究 [J], 解立峰;刘家骢;果宏
4.激波管内燃料两相云雾爆轰的实验研究 [J], 姚干兵;解立峰;刘家骢
5.燃料空气炸药二次起爆引信与爆轰云雾动态交会浓度检测方法 [J], 吉童安;娄文忠;付胜华;汪金奎;刘伟桐
因版权原因，仅展示原文概要，查看原文内容请购买。
滑移爆轰作用下飞板运动的解析解
滑移爆轰作用下飞板运动的解析解
谭多望;孙承纬
【期刊名称】《高压物理学报》
【年(卷),期】1999(013)002
【摘要】给出了滑移爆轰作用下飞板运动的近似解析解,应用其结果计算炸药滑移爆轰时飞板的运动,同实验结果和二维数值模拟计算结果吻合,表明该近似解析解具有较高的精度.解析解简洁,便于工程应用.
【总页数】7页(P120-126)
【作者】谭多望;孙承纬
【作者单位】中物院流体物理研究所冲击波物理与爆轰物理实验室,成都523信箱,610003;中物院流体物理研究所冲击波物理与爆轰物理实验室,成都523信箱,610003
【正文语种】中文
【中图分类】O383.3
【相关文献】
1.滑移爆轰载荷作用下复板承受的冲击压力的数值计算 [J], 谢飞鸿;罗冠炜;孙伟;廖军生
2.滑移爆轰作用下飞板运动规律的研究——双P—M法 [J], 郭佑雄;赵福兴
3.在两层不同爆速炸药滑移爆轰作用下飞板运动的数值计算研究 [J], 张凯;李晓杰
4.滑移爆轰作用下飞板运动规律的特征线差分法研究 [J], 赵春风;李晓杰;于娜
5.滑移爆轰驱动下飞板运动姿态的连续电阻测试法 [J], 李科斌;董新龙;王永刚;陈翔;李晓杰
因版权原因，仅展示原文概要，查看原文内容请购买。
内爆炸后燃烧反应率时间历程理论预估及能量释放常数确定方法
内爆炸后燃烧反应率时间历程理论预估及能量释放常数确定方法徐维铮; 吴卫国; 况正【期刊名称】《《中国舰船研究》》【年(卷),期】2019(014)004【总页数】8页(P22-29)【关键词】约束空间内爆炸; 后燃烧效应; 反应率模型; 五阶WENO格式; 化学反应动力学; 能量释放常数【作者】徐维铮; 吴卫国; 况正【作者单位】高性能舰船技术教育部重点实验室湖北武汉430063; 武汉理工大学交通学院湖北武汉430063【正文语种】中文【中图分类】U661.40 引言针对约束空间内爆炸问题，爆炸后燃烧效应是需要重点关注的内爆炸物理现象。
大多数炸药（例如TNT、温压炸药、SDF（Shock-Dispered-Fuel）混合型炸药［1］等）的爆轰产物具有负氧性，在膨胀的过程中会与周围空气中氧气发生剧烈的燃烧反应，释放出大量能量，形成爆炸后燃烧效应［2］。
当爆炸发生在约束空间内时，因壁面限制，相比空爆会进一步加剧后燃烧效应，从而对冲击波的壁面反射压力、准静态压力等产生一定的增强效应。
迄今，针对约束空间内后燃烧效应对爆炸载荷的影响规律，国内外学者开展了大量的实验和数值模拟研究。
国内的相关研究主要集中在西安近代化学所和西北核技术研究所。
例如，姬建荣等［3］在小型爆炸容器内，实验测量了在密闭和半密闭爆炸工况下TNT-AL炸药的后燃烧性能，结果表明当铝粉质量分数为20%时，TNT-AL炸药在小型爆炸容器中较易发生后燃烧反应。
研究人员用高速摄影仪观测到了后燃烧反应及其发展过程，发现后燃烧反应在爆炸腔内形成了一定的准静态压力，持续时间达数十毫秒。
钟巍和田宙［4-5］基于理想气体状态方程和能量守恒定律，将准静态压力划分为爆炸压力（等熵假设/等压假设）和化学反应压力2个部分进行计算，推导了考虑爆炸产物发生化学反应后在密闭空间内的准静态压力计算公式；金朋刚等［6］在密闭容器内通过实验测量了含有3种粒度（0，13，130 μm）的铝粉HMX基炸药的准静态压力，结果表明，铝粉可以提高HMX基炸药在密闭空间爆炸的准静态压力，小颗粒铝粉相比大颗粒铝粉更有利于提高准静态压力。
爆轰问题的一个高效二维离散玻尔兹曼模型
爆轰问题的一个高效二维离散玻尔兹曼模型林传栋;许爱国;张广财;李英骏【期刊名称】《凝聚态物理学进展》【年(卷),期】2015(4)3【摘要】本文构建了多松弛时间离散玻尔兹曼模型,并使用该模型模拟爆轰现象。
相对于我们之前的一个模型[Xu A., Lin C., Zhang G., Li Y., Phys. Rev. E 91 (2015) 043306],本模型在模拟有化学反应或无化学反应流体系统时的计算效率更高。
这是因为前者使用了24个离散速度,而本模型只使用16个。
在模拟部分高马赫物理系统时,本模型表现出更高的数值稳定性。
使用该模型,本文分四种情况模拟了爆轰波激发的Richtmyer-Meshkov不稳定性问题。
当爆轰波由反应物传向另一种较轻的不反应的物质时,由于突然失去能量补充,温度急剧下降,在物质界附近将会出现一层高密区域。
【总页数】10页(P102-111)【作者】林传栋;许爱国;张广财;李英骏【作者单位】[1]中国矿业大学(北京)深部岩土力学与地下工程国家重点实验室,北京;;[2]北京应用物理与计算数学研究所计算物理重点实验室,北京;;[3]北京应用物理与计算数学研究所计算物理重点实验室,北京;;[1]中国矿业大学(北京)深部岩土力学与地下工程国家重点实验室,北京【正文语种】中文【中图分类】O35【相关文献】1.无起爆药雷管燃烧转爆轰的二维两相流数学模型 [J], 李梦娇;王雨时;闻泉2.化学反应速率对爆轰特性的影响研究—基于离散Boltzmann模型 [J], 张玉东;许爱国;张广财;祝成民;;;;3.无起爆药雷管燃烧转爆轰的二维两相流数学模型 [J], 李梦娇;王雨时;闻泉;4.爆轰模型熵问题的探讨 [J], 胡绍鸣;李辰芳5.一个高效的离散Boltzmann爆轰模型 [J], 张玉东;许爱国;张广财;祝成民因版权原因，仅展示原文概要，查看原文内容请购买。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
多舱室内气云爆炸爆轰产物运动二维数值模拟
2 2 徐维铮 1， ， 吴卫国 1，
1 武汉理工大学 高性能舰船技术教育部重点实验室， 湖北 武汉 430063 2 武汉理工大学 交通学院， 湖北 武汉 430063
摘
式， ቤተ መጻሕፍቲ ባይዱ过耦合求解爆轰产物质量分数输运方程与可压缩欧拉方程， 自主编写多舱室内气云爆炸过程二维数值计
第 13 卷 增刊 1 2018 年 12 月
中 国 舰 船 研 究 Chinese Journal Ship Research 中 国 舰 of 船 研 究
Vol.13 Supp 1 Dec. 2018 第 13 卷
引用格式： 徐维铮， 吴卫国 . 多舱室内气云爆炸爆轰产物运动二维数值模拟 ［J］ . 中国舰船研究， 2018， 13 （增刊 1） ： 92-99. XU W Z， WU W G. A two-dimensional numerical simulation on dynamic motion of detonation products generated by gas cloud explosion inside the multi-cabin ［J］ . Chinese Journal of Ship Research， 2018， 13 （Supp 1） ： 92-99.
势， 基 于 Fortran 平 台 ， 采 用 五 阶 WENO 有 限 差 分
鉴于 WENO 格式具有精度高、 稳定性好的优
93
格式， 本文拟自主开发多舱室内气云爆炸过程的 二维数值计算程序， 开展几种典型爆炸工况下爆 轰产物的运动规律研究。
1
控制方程
为了计算气云爆轰产物的运动规律， 对爆轰
产物质量分数输运方程与无黏性可压缩欧拉方程 进行耦合求解。考虑到爆轰产物与周围空气燃烧 过程的复杂性， 本文不考虑多组分的复杂燃烧化 学反应以及热传导、 热辐射等效应。针对后燃烧 增强效应， 提出采用简化的反应率模型［14］进行近 似， 初步探讨后燃烧效应对爆轰产物在周围空气 中扩张过程的影响。方程的具体形式如下： ¶U + ¶E x + ¶F = S （1） ¶t ¶x ¶y 其中，
Abstract： ［Objectives］In order to study on the dynamic motion of detonation products generated by gas cloud explosion inside the multi-cabin， ［Methods］ the fifth-order WENO finite difference scheme is adopted to solve the coupling of mass fraction transport equation and compressible Euler equation for the detonation products，and then a two-dimensional numerical simulation program for the gas cloud explosion in the multi-cabin is developed successfully. The reliability and correctness of this numerical simulation program is verified by virtue of the plane shock diffraction rectangular step test，after which the verified program is used to carry out the numerical simulation of the gas cloud explosion under different explosion conditions， so as to explore the influence of venting hole and afterburning energy on the dynamic motion process of the detonation products during the explosion inside the multi-cabin.［Results］The results show that increasing the size of the venting hole can accelerate the emission process of the detonation products； increasing the reaction rate constant can accelerate the expansion of the detonation products to the adjacent cabin.［Conclusions］This study can provide a reference and guidance for further exploring the expansion process of the detonation products generated by the complex multi-component combustion and the evaluation of the explosion fireball damage. Key words：gas cloud explosion in multi-cabin；blast wave；fifth-order WENO scheme；numerical simulation； detonation products
（2）
（3） （4）
舱焊补过程中， 电火花点然了舱室内部的可燃蒸 汽， 引起强烈爆炸， 并殃及泵舱和货油舱发生连续
式中：ρ 为密度； u v 为 x y 方向上的速度分量；
p 为流体压力；E 为单位体积流体的总能量； e为
比内能；λ 为爆轰产物的质量分数； α 为可燃气云
a 为反 的反应率 （初始 α = 1 ， 反应完成后 α = 0 ） ；
文献标志码： A DOI： 10.19693/j.issn.1673-3185. 01186
A two-dimensional numerical simulation on dynamic motion of detonation products generated by gas cloud explosion inside the multi-cabin
要： ［目的］为了研究多舱室空间内气云爆炸后爆轰产物的运动规律， ［方法］采用五阶 WENO 有限差分格
算程序。利用平面激波绕射矩形台阶试验验证本文数值程序的可靠性与正确性。采用验证的程序开展不同爆 炸工况气云爆炸过程数值计算， 探讨多舱室内爆炸过程中泄压、 后燃烧能量加入对爆轰产物运动过程的影响规 律。 ［结果］结果表明： 增大泄压口尺寸， 可加剧爆轰产物泄出过程； 增大反应速率常数， 可加速爆轰产物向相邻 舱室的扩张过程。 ［结论］研究成果可为进一步考虑多组分复杂燃烧效应后爆轰产物扩张过程及毁伤评估提供 一定的参考和指导。 中图分类号 ： U674.7 关键词 ： 多舱室内气云爆炸； 爆炸波； 五阶 WENO 格式； 数值计算； 爆轰产物
基 金 项 目 ：国 家 部 委 基 金 资 助 项 目 ； 国家自然科学基金资助项目 （51409202） ； 中央高校基本科研业务费专项资金资助 作者简介： 徐维铮 （通信作者） ， 男， 1991 年生， 博士生。研究方向： 爆炸波数值计算方法及程序开发研究。 吴卫国， 男， 1960 年生， 教授， 博士生导师。研究方向： 结构动力学， 舰艇抗爆抗冲击。 （2016-YB-016） 收稿日期： 2018 - 01 - 29 网络首发时间： 2018-11-28 16:10
æ ρu ö æ ρv ö æ ρö ç ρu 2 + p ÷ ç ρvu ÷ ç ÷ ç ρu ÷ ç ÷ ç 2 ÷ ç ρv ÷ ÷ ç ρuv ÷ ρv + p ÷ ç ç U=ç çE÷ ÷ E x = ç u( E + p) ÷ F = ç v ( E + p) ÷ ç ç ρλ ÷ ÷ ç ρuλ ÷ ç ρvλ ÷ ç ç ÷ ÷ ç ÷ ç ÷ è ρα ø è ρuα ø è ρvα ø 0 æ ö ç ÷ 0 ç ÷ ç ÷ 0 ç ÷ 1 1 ç ÷ ç S = ρaα 2 p 6 Q ÷ ç af ÷ ç ÷ ç ÷ 0 ç ÷ 1 ÷ 1 ç ç ÷ 2 6 è ρaα p ø E = ρe + 1 ρu 2 + 1 ρv 2 2 2 p = ( γ - 1) ρe
2 2 XU Weizheng1， ， WU Weiguo1， 1 Key Laboratory of High Performance Ship Technology of Ministry of Education， Wuhan University of Technology， Wuhan 430063， China 2 School of Transportation， Wuhan University of Technology， Wuhan 430063， China
增刊 1
徐维铮等： 多舱室内气云爆炸爆轰产物运动二维数值模拟
0
引
言
工业可燃性气云爆炸事故大多由弱点火点燃 可燃气体引起， 以亚音速传播的爆燃波为主。爆 燃波的传播过程复杂， 受环境条件和物理因素 （爆 燃转爆轰物理机制） 的影响极大。在传播过程中， 气体由于障碍物阻碍发生湍流加强效应或是在初 始时刻直接起爆， 其爆炸模式将由爆燃转为爆轰 或直接以爆轰模式发生爆炸。发生爆轰以后， 气 体将在空气中形成强度较大的带有负压区的空气 冲击波， 对工作人员和周边结构设施将造成较大 的危害。针对空中可燃气云和云雾爆炸模拟， 近 年来一些学者进行了相应的数值研究［1-4］， 他们将 可燃气云、 云雾简化为高压气体， 通过数值模拟， 研究了近地面气云爆炸波反射和气云相互作用的 过程。针对管道和单一舱室内部气云爆炸过程， 文献 ［5-10］ 进行了一系列的试验和数值仿真研 爆炸冲击波载荷及传播过程， 较少关注爆炸后爆 轰产物的运动规律。在实际爆炸过程中， 由于爆 轰产物处于高温、 高压状态， 且具有易燃性， 极易 与周围空气中的氧气发生剧烈的化学反应并释放 能量， 进一步增加冲击波的强度或形成高温场和 极强的热辐射。对于多舱室空间， 爆轰产物还将 扩展到其他舱室内部， 进一步扩大毁伤范围。尤 其是在船舶这种典型多域连通空间内非常容易发 生火灾， 而且一旦发生气云爆炸， 轻者货损， 重者 伤人， 甚至导致船舶全损。例如， 1967 年 7 月， 美 国“福莱斯特” 号航母发生火灾爆炸， 造成 134 人 死亡、 62 人受伤、 26 架飞机烧毁， 舰体的 10 层甲 板有 6 层损坏； 1982 年 5 月， “ 大庆 53” 号油轮在机 爆炸， 造成全船沉没， 造成 20 人死亡， 直接经济损 失 达 1 400 万 元 ， 并导致严重的海洋污染。鉴于 爆轰产物的运动规律及后燃烧能量加入的影响。 爆炸属于强间断问题， 其数值模拟对激波捕 捉格式提出了较高的要求。 1994 年， Liu 等［11］提出 了加权本质无振荡 （Weighted Essentially Non-os⁃ NO 格式作为一种典型的高精度激波捕捉格式， 对 此， 有必要探讨和研究多舱室空间内气云爆炸后 究。然而， 这些研究主要集中在气云爆轰形成的