逆向阴燃波的结构和稳定性：热损失效应 袁枋平.doc

中国工程热物理学会 燃烧学 学术会议论文 编号：114232 逆向阴燃波的结构和稳定性：热损失效应 袁枋平， 卢占斌 （上海大学应用数学和力学研究所，上海市延长路149号，200072） （电话：0215633-3707，E-mail:zblu@shu.edu.cn）[6]的实验阐明了逆向阴燃波传播的物理机理，指出受迫流动的空气对于逆向阴燃波的传播具有双重作用。此后Schult et al[7]提出了与上述实验结果相符的数学模型，用以解释平面逆向阴燃波的结构和传播机理。1998年，Olson et al[8]和Zik et al[9]两个研究组在微重力实验中分别独立发现薄固体燃料的逆向阴燃在特定来流条件下会表现出指进不稳定性。除此之外，在Lu et al关于逆向阴燃稳定性的一项工作中[10]，数值模拟结果显示当来流空气流量趋近熄灭极限时，在某些参数条件下平面阴燃波会发生断裂（fragmentation）失稳，失稳过程中伴随着急剧的局部升温，且升温幅度最高可达700摄氏度，足可以触发向火焰燃烧的转化。 文献[10]中采用了一个绝热模型，即没有考虑到阴燃过程中热损失效应的影响。本文我们将研究非绝热逆向阴燃波，也就是考虑阴燃向外界的对流热损失效应，相应的阴 基金项目：国家自然科学基金（No和上海市科委启明星计划（No.09QA1402300）资助。 燃结构示意图如图1所示。我们将采用数值计算的方法研究热损失效应对逆向阴燃的结构和稳定性的影响。 图1：非绝热逆向阴燃结构示意图 1 数学建模与控制方程 考虑如图1所示的一个无穷长的薄固体可燃孔隙材料的阴燃，空气由左向右流过固体燃料的表面。采用如下的单步不可逆Arrhenius化学反应模型： （1） 其中是基于质量的各组分的反应配比系数，是单位质量的固体燃料完全反应后释放的热量。这里我们采用文献[10]中所采用的简化方式，即假设,由此可以得到简化后的非绝热逆向阴燃的控制方程如下： （2） 分别对应于能量守恒方程以及固体燃料和氧气的组分守恒方程。这里，为反应混合物的温度，为固体燃料的质量百分比，为氧气在气体混合物中的质量浓度，为固体（包括固体燃料和固体反应产物）的密度，为气体混合物的密度，为比热，为热传导系数，为氧气在气体混合物中的质量扩散率，和分别为气体和固体的质量流量，且二者均为常数，为固体燃料的厚度，为空气的对流传热系数，为无穷远处外围空气的温度，为反应速率，定义为： （3） 其中是一个反应常数，为反应的活化温度。 我们选择如下特征长度、特征时间，以及特征温度尺度： （4） 其中，是在无穷远处来流空气中氧气的浓度。定义以下无量纲变量： （5） 由此可以得到如下无量纲形式的控制方程： （6） 其中。除此之外，我们还引进了三个重要的无量纲参数： （7） 其中为Damkohler数，为Lewis数，而为无量纲化的对流热损失系数。进一步假设对流热损失系数，其中是一个常比例系数，则可以变化为。 与无量纲化方程（6）对应的方向边界条件为： （8） 2一维定常解 在本节中我们将通过数值方法来模拟一维逆向阴燃波（绝热和非绝热）的结构。在数值方法上，对方程组（6）和边界条件（8）在空间上采用四阶中心差分离散，时间上采用一阶Euler格式，并通过迭代法求解。为了能更直观的表现出阴燃波的温度，所以在这里不对温度无量纲化。控制域内的边界条件改写为： （9） 物理参数的选取如下表所示： 表1 物理参数取值 物理参数 取值 C 1.1 J/g/K ? 4.6×10-4 W/cm/K Ta 8200 K B 2.53×105 g/cm3/s D 0.453 cm2/s T0 300 K ?s 0.0265 g/cm3 ?g 1.0×10-3 g/cm3 yO0 0.235 ? 0.41 Q 1500J/g k 0.03 图2给出了逆向阴燃波随来流空气流量变化时的响应曲线，分别以阴燃的最大温度和固体燃料的质量通量为代表。图3分别给出了位于富燃料分支和贫燃料分支上两个代表性时的温度剖面。很明显，对于绝热情况，反应层后面具有恒定的阴燃反应温度；然而当考虑对流热损失时，反应层后的温度剖面开始逐渐衰减。 3 线性稳定性分析 本节中将通过数