基于幂律褶皱模型对随机填充床内湍流预混火焰的孔隙尺度模拟研究.docVIP

中国工程热物理学会 燃烧学 学术会议论文 编号： 基于幂律褶皱模型对随机填充床内湍流 预混火焰的孔隙尺度模拟研究 姜霖松*，刘宏升，索少一，解茂昭 （大连理工大学 海洋能源利用与节能教育部重点实验室 辽宁 大连 116024） （TelEmail： 1160819157@） 摘 要：采用离散元软件得到了与实际堆积结构相近的3-D几何模型。应用大涡方法结合幂律火焰褶皱模型，模拟了氢气-空气预混气体在多孔介质随机结构下的燃烧过程，并与实验对比验证了模型的有效性。对球形火焰在不同结构内的分布情况，厚度，传播速度及火焰面半径随时间变化情况进行了对比分析，给出了火焰面上Karlovitz数的分布情况，对多孔介质特殊结构内湍流-火焰相互作用情况的定量分析，并根据该结果得到了多种结构内的湍流火焰分区情况。 关键词：随机结构；球形火焰；幂律火焰褶皱模型；湍流-火焰相互作用 0 引言 近年来，随着对多孔介质燃烧技术的认识不断加深，一些学者尝试将此技术应用于内燃机、页岩油开采及燃料电池等相关领域。而对多孔介质材料内预混火焰的研究成为多孔介质燃烧中的研究重点之一。 Schulze等[1]提出了一个粒子求解的CFD模型以对多孔介质堆积床内部的流体流动、传热、传质以及化学反应进行研究。结果表明应用文中的基本模型得到的结果与随机堆积床实验中得到的结果吻合较好，同时发现堆积床粒子的分布对流场，传热传质影响较大，Schulze等通过研究火焰长度对以上影响做出了详细分析。Liu 等[2]以计算颗粒流体动力学（CPFD）为基础，应用欧拉-拉格朗日方法建立了一个三维可压缩循环流化床（CFB）的燃烧模型，为后来的流化床反应器构建奠定了基础。Schwebel等[3]分别将流化堆积床和固定堆积床应用于氧化实验以做出对比，并对差异进行定量分析，结果表明固定堆积床在气体转化率方面具有较大的优势，即便应用相较于流化堆积床更粗的颗粒，其中的气体转化率依然较高，然而在固定堆积床中进行氧化，颗粒的烧结是一个重大的难题。 考虑到封闭腔内多孔介质下燃烧研究十分罕见，且封闭腔内球形火焰又极具特点，下面仅给出一些前人对空腔内球形火焰的研究成果，Tholin et al. [4]对纳秒火花点火器点燃的贫燃预混氢氧球形火焰进行了二维数值模拟，结果表明在初始温度1000K和一个大气压工况下，纳秒火花点火产生的热效应和化学效应是相差无几的，研究中应同时考虑热效应和化学效应。Lawes等[5]利用高速纹影图像推导对湍流燃烧弹中膨胀的异辛烷-空气预混球形火焰的湍流燃烧速度进行了测量，测量的主要工况为U0 = 0.5，1，2，4，6m/s，当量比为0.8，1.0，1.2，1.4，压力为p = 0.1，0.5，1 MPa。并将得到的结果与一些现有的湍流燃烧速度相关的研究结果进行对比，对差异进行了分析和讨论。Bradley等[6]应用内爆技术实现了对更广阔的当量比和压力范围内湍流扩散燃烧中火焰传播速度的测量，同时给出了新的火焰分区. 资助项目：国家自然科学基金面上项目( 资助项目：国家自然科学基金面上项目 1 物理模型与计算方法 1.1 物理模型 图1 计算域示意图 （a） 堆积床整体结构示意图 （b1） 球形空腔计算域 （b2），（b3） 多孔介质球形空腔计算域 在我们以前的研究[7]中，使用离散元软件LIGGGHTS[8]来建立随机堆积床几何模型（图1（a）），通过对真实的堆积过程进行仿真，使模型相比传统的插排或顺排有序模型更接近真实的堆积结构。 如图1所示，本文通过LES（大涡模拟）方法对氢-氧预混燃烧过程进行介尺度模拟，为进一步分析燃烧初期火焰传播的细节，并节省计算时间，从堆积床上切割出半径为10mm的球形空腔，得到图1（b2），（b3）所示计算域，同时将具有相同半径的球形封闭腔（图1（b1））作为对比结构（无多孔介质结构）。 文中初始气体温度、初始壁面温度、初始当量比、初始压力和初始反应进程变量（C）分别设置为300 K、300 K、1、0.1 MPa和0，并通过设置计算域几何中心处半径为1mm的球形高反应进程变量（C＝0.5（21））区域（图1（a））来得到接近实验（准层流状态）的点火源。壁面均采用由Navier-Stokes特征边界条件（NSCBC）给出的等温和无滑移条件，壁面均不参与化学反应，所有组分的质量分数梯度在边界处均为零。 1.2 计算方法 本文中反应区内