不可压湍流反应流直接模拟和RANS二阶矩燃烧模型检验.pdfVIP

中国工程热物理学会 燃烧学 ·学术会议论文 编号：074005 不可压湍流反应流直接模拟和 RANS二阶矩燃烧模型检验 王方h2，周力行1，许春晓1，黄勇2， 1清华大学工程力学系，j匕京100084；2北京航空航天大学热能工程系，北京lo0083 010鲰协lg@buaa．edu．cn 摘要：用谱方法对三维槽道不可压湍流反应流动进行了直接模拟，并用直接模拟数据库对关联量精 确输运方程中各项的贡献进行统计，发现产生项和耗散项起重要作用。用直接模拟数据库对RANs模 拟中的二阶矩(SoM)湍流燃烧模型进行检验，发现扩散项和产生项的模拟在大部分区域都和精确值 接近，但在近壁区需要修正；耗散项的模拟有待改进。代数二阶矩模型(As0M)的模拟值和精确值 符合很好。 关键词直接模拟，湍流反应流动，二阶矩燃烧模型 1．前言 ’湍流反应流动和燃烧的直接数值模拟(DNS)不需要任何湍流模型和化学反应模型， 能够求解所有尺度的湍流结构，因此是基础研究的有力工具，最近几年方兴未艾。虽 然计算技术突飞猛进，但是DNS对计算存储和时间需求很大，尚不能直接解决：【程问 题，目前主要研究以下几方面：①湍流燃烧的瞬时细观流动结构和火焰结构…；②燃 后一个方向是目前研究的热点。不过目前直接模拟的研究多数局限于PDF输运方程模 型、条件矩模型和火焰面模型。 目前，RANS模拟仍然是湍流燃烧工程预报的主要方法。然而广泛应用的 结果出入很大。条件矩模型和火焰面模型的前景很好，但是仍然有待改进。PDF输运 方程模型中的化学反应项是精确的，不需要封闭模型，可以与详细化学反应动力学结 合，但是它的小尺度混合模型有待改进，而且计算耗费很大。为此周力行等提出了二 阶矩(SOM)湍流化学反应模型，已经成功应用于甲烷空气射流燃烧和甲烷空气旋流燃 烧模型和设定的PDF模型，又比PDF输运方程和条件矩模型简单。但是SoM模型的封 闭假设需要进一步的理论检验。 本文使用谱方法对三维不可压湍流反应流动进行直接数值模拟，并用DNS数据库 检验RANS模拟的SoM湍流燃烧模型。 2．直接模拟的控制方程和数值方法 对有组分扩散、传热及化学反应的不可压反应流动，忽略密度的变化，模拟用的 246 瞬态万程为： 垒：o (1) ^ 徼i 鲁M筹一言+毒cy譬+詈，， c2， 昙(z)+毒(啪)=岳(。筹)一w。／p (3) 昙(K)+丢(啪)=丢(。等)_叫p (钔 昙(D+毒(垆)=毒(专考)+嵋蜴饵力、 (5) 所模拟的是充分发展的湍流反应流动，计算域和坐标系统如图1所示，x，y，z 分别表示流向、法向和展向坐标，槽道法向高度为2H，槽道的流向长度是12．6H，槽 流向和展向采用128个平均分布的网格节点，法向采用壁面加密的129个网格节点， 同时采用三阶精度的时间分裂格式。时间步长为H／U。的O．Ol倍。 图1计算域示意图 模拟中，上壁面处反应物1的浓度质量分数Y1取为1．O，反应物2的质量分数Y2 壁面处速度取为0。模拟工况见表1。其中工况1’2，3分别是不同的指数前因子和活 3．湍流燃烧RANs模拟的二阶矩精确输运方程及其封闭模型 首先给出精确的二阶矩输运方程。如果瞬时的化学反应率为： K=Bp2‰‰exp(一言)