球形预混火焰的热扩散不稳定性理论研究.docxVIP

中国工程热物理学会 燃烧学 学术会议论文 编号：184008 球形预混火焰的热扩散不稳定性理论研究 蔡雁南1，李挺1，杜明龙2，杨立军2 （1.北京航空航天大学能源与动力工程学院，北京 100191； 2.北京航空航天大学宇航学院，北京 100191） (Tel Email:caiyannan2014@) 摘 要 本文研究了预混可燃气体在球形多孔稳焰器外形成火焰锋面的热扩散不稳定性问题。求解并获得了稳定状态下火焰温度和燃料质量混合比随火焰半径变化的解析解，并通过MATLAB中的bvp4c程序包对球形火焰的稳定性进行了分析，详细探讨了Lewis数和表征球形多孔稳焰器的相关物理参数(孔隙率和多孔壁厚度)对火焰锋面稳定性的影响。 关键词 球形预混火焰；热扩散不稳定；稳定性分析 0 前言 火焰锋面的稳定性问题属于火焰的本征不稳定性，与外界声扰动和燃烧室壁面的相互作用无关。火焰本征不稳定性是火焰动力学的一个经典课题之一[1]。在过去的几十年里，为了揭示火焰锋面稳定性的机理，相关学者们进行了大量的理论和实验研究。预混气体在多孔稳焰器外的一定距离能够形成一个稳定的火焰锋面[2]。但是对于多孔稳焰器，预混气体需要通过多孔壁面扩散至多孔稳焰器外进行燃烧。其质量输运和多孔稳焰器自身的性质（孔隙率和壁厚）对火焰稳定性具有很大的影响，需进行深入研究。由于燃烧过程的复杂性，目前关于多孔稳焰器火焰锋面稳定性的理论研究主要集中在一维层流预混火焰，且关于平面多孔稳焰器预混火焰的研究较多[3,4]，圆柱次之[5]，尚未见到针对球形预混火焰的研究成果。因此本文将对球形多孔稳焰器外球形火焰锋面的稳定性进行理论分析及研究。 1 数学模型 本文建立了球形多孔稳焰器的预混火焰在球坐标系r,θ,φ下的物理模型，原点o设置在球形稳焰器的中心，r轴与稳焰器的半径平行，θ为仰角，φ为方位角。多孔稳焰器的半径为r0，壁厚为h。可燃预混气从多孔壁内流出，在多孔稳焰器外壁面下游形成半径为rf的球形火焰锋面。该模型忽略辐射热损失，化学反应为单步不可逆燃烧反应。稳定性分析基于火焰薄片模型。其最早可以追溯到Zeldovich[6]。火焰被看作是一个包含预热区和反应区的组合，反应区很薄，在大活化能假设下，反应区被看成是一个反应面r=rf，被称为火焰锋面。预热区由对流和扩散作用主导，几乎没有化学反应。燃烧化学反应主要集中在火焰锋面，产生了很大的浓度和温度梯度，导致强烈的质量和热量交换，反应物在反应区内被完全消耗，产生大量的热量。T和Y分别是未燃预混气体中不足组分的温度和质量混合比。图1展示了T和Y随火焰半径的分布情况。根据火焰锋 基金项目：国家自然科学基金青年基面上反应物的质量和能量守恒关系可以得到以下跃变条件： -ρD?Y?r=ρY0S qρD?Y?r+ 其中，*=*r=rf+-*r=rf-表示火焰锋面上某参数的跃变。ρ、Dth、D、λ和Cp分别是预混气体的密度、热扩散系数、质量扩散系数、导热系数和定压比热容，均假设为常数。E是单步不可逆燃烧反应的反应活化能，R是通用气体常数，q是单位质量燃料的放热率，Y0为初始反应物的质量混合比，T0是初始反应物温度，Ta=T0+qY0/Cp是绝热火焰温度，Tf是r=rf处的瞬态火焰温度，SL是层流火焰速度。 图1球形多孔稳焰器火焰的物理模型及T、Y随r的分布曲线 假定层流火焰速度SL为单位速度，球形绝热火焰的热厚度rth=Dth/SL作为单位半径，tth=rth/SL作为单位时间，得出无量纲温度、质量混合比、时间、球形多孔稳焰器壁厚和半径分别为T=T-T0 由于多孔稳焰器的存在，多孔稳焰器外壁面和火焰锋面将燃烧过程划分成三个空间区域I（r0-h≤r≤ I：定义多孔稳焰器孔隙率为μ=V0-VV0，其中V0为材料在自然状态下的体积，V T=0 (3 ?Y?t II： ?T?t+mr ?Y?t III：对于已燃尽气体 T=Tf； 其中m= 无量纲边界条件如下： 当r=r0-h时， T=0 , 当r=r0时， T=0 , T 当r=rf时， T= 当 r=+∞时， T=Tf 在r=rf火焰锋面处，由方程(1)、(2) ?T?r+ 其中γ=Ta- 2模型求解 在线性稳定性分析中，假设火焰的扰动属于谐波扰动，火焰半径、温度和质量混合比可写成稳定状态和小扰动状态下解的叠加形式 rf=rfs T=Ts+ε Y=Ys+ε 其中ε 火焰初始扰动的幅值