CH4—O2混合气中爆燃爆震转捩数值模拟.pdfVIP

维普资讯 翻D 爆 炸 与 冲 击 Vo】．18，No 4 1998年 10月 EXPIOS10N AND SHOCK WAVES oct ，1998 CH一O混合气 中爆燃爆震转捩的数值模拟 高秦荫 黄军涛、， 元明 高 越 信 一 兵 (东北大学 沈阳 110006) O32J 摘要 运用化学流体力学基本理论和两步燃烧反应模型原理，建立了一维封闭体系可燃气爆 燃爆震转变现象的数学模型，利用拉格 朗 日质量坐标变换下的I．ax—Wendroff和 Mecorrnic气动差 分与 Adams化学差分格式 ．求解基本方程 ．成功地完成 了过程 的数埴模拟 ，清楚地说 明了可燃气 中DDT现象 由压缩波到激波达到稳态爆震 的发生机制和火焰带引生爆震波的过程行为 。 l 引 言 建立可燃混合物爆炸发生的临界条件是工业爆炸灾害预报与防护的基础。以拄 ，这一工作 主要通过物理实验进行 ，并曾试图建立获得爆炸临界参数的实验标准 。但是 ，迄今为止的推 荐实验方法都是以等容绝热为条件 ，这与工程实际中常常面临复杂多变的流动传热状况相差 甚远0 后来，开始有人在限定条件下实现可燃混合物中的爆震，实验研究其发生因索和发展 过程 ，丰富了人们对爆震波形成机理和传播结构的基本认识 ，但是对于有工程实际那样复杂的 边界条件和初始条件的实验相似尚难以实现 ]。 可燃混合物 中燃烧波爆燃爆震转捩与过渡 (DDT)特性受系统初始条件和边界条件的显 著影响。正因为这一点+物理实验有许多困难 。随着数值方法和计算技术的进步 ，近年来，DDT 问题的数值分析与模拟研究 日渐活跃 ，涉及的对象有均相混合气 ，也有气一固混合物，研究的问 题有DDT发生条件 ，也有爆震波结构 。已有 的研究预示着数值模拟研究 DDT 问题的前 景，同时也期待在现象本质描述和过程解析计算上有更多的尝试 “ 我们针对CH O混合气燃烧波在一维封闭空问传播的DDT 问题．运用化学反应流动传 热基本理论和两步反应模型原理，建立现象的数学表达，利用拉格朗日质量坐标变换方法，完 成过程的数值解析，为深人研究DDT过程发生的化学反应流动传热本质特征提供基础 。 2 物理模型及其数学表达 2．1 物理模型 讨论 CHrO混合气在一维封闭空间点火形成燃烧波后 ，由爆燃转成爆震的发生过程 ．在 如图l简化条件下 ，考察其化学反应流动参数的基本特征+其假定条件是 ： 国家 自然科学基盎资助项 目 高泰荫：男 1943年 月生、教授 I997一岫一r,i】日收到原稿 ．1997—12-2s收到修改稿 。 维普资讯 爆 炸 与 冲 击 第 l8卷 1_反应装置是一根两端封 闭、边界绝热且光滑的圆柱形 钢管 ，管子的抗压强度足够大，且长度L远远大于直径d。 2．初始 CH一O 气体混合物 的温度、浓度和压力沿径 向 … n §． 始掘度场由两部分组成，即T—T()：当0≤ ≤一 ，． 时，为恒定 的高温区域 ，即点火源，71—71h~COIlst；当 ≤ ^ a(pY )+ 一(puYD ： t 此式描述二步反应模 型 的第二 阶段 ：放 热阶段 。其 中：y，为弟 i种 组分 的质量 日分 比维魇 ．且 yr一 1(N一 4)， 为 i组分的反应率 ·