平面自由湍射流拟序结构的大涡模拟的研究.pdfVIP

中国工程热物理学会 燃烧学学术会议 平面自由潴射漉拟序结构的大涡模拟研究。 范奎林”张会强郭印诫王希麟林文漪 ● 她京·清华大学工程力学系(100084) ● ”Tcl：010-6278162l ne‘ E-majl：劬一qu如1in@263 ／ Ⅱ摘要a本文对空间发展的平面不可压缩湍射流拟序结构的非定常演化过程进行 了走涡模拟．采用标志物浓度等值线分布示踪平面射流气相流动，(Re=儿，300’D模拟 结果成动再现了平面射流中Kelvin—Helmholtz不稳定性的发生与展向太尺度涡的卷 起、合并、破碎过程．模拟工作捕获到了射流拟序结构剧烈相互作用地“偶极子”和 “三板子”现象．在拟序鲒构的配对过程中，大尺度旋涡的尾迹混合形成了流向“发 卡”型拟序涡结构．平面射流拟序结构的时空演化过程所呈现出的卡门涡街特扯主要 集中于初始段的后部和过渡段的前端．敷值流场显示结果和前人染线法的实验结果吻 合得很好． B关键词Ⅲ平面射流；拟序结构；大涡模拟；数值模拟；不可压缩流动 1、引言 作为剪切流的一种模型流动．平面射流中的拟序结构一直是拟序结构领域研究的 热点之一IJ．q。但由于作为示踪剂的烟粒子等标志物在空气中具有很快的扩散率，常常 导致欲观察的流动结构与环境介质之间反差(对比度)很小．进而不易分辨剪切层卷 起过程中的涡结构，使得气相射流的实验研究相对来说比较困难，实验成本昂贵：同 时由于实验条件的差异和和观测视角的不同，业已公布的平面射流拟序结构的结果尚 不完整。 ， 而在数值模拟中，人们可以方便地控制初始条件和边界条件．并考察它们对流动 特性的影响。因此近年来数值模拟和试验研究一起，共同成为全面了解流体动力学特 性的有效途径。直接数值模拟(DNs)求解了流动中的所有尺度．是未来踹流数值模拟的 ■ 一个重要发展方向。但遗蟪的是，一方面，目前的计算机性能使得D咐S仍然只能局限 于相当低Rcy∞lds数和规则几何形状的流动问题。另一方面．肼qS在模拟自由湍流流 ■ 动时还有它特殊的困难：第一，在自由湍流流动中．没有明确的、流动受到限制的边 界．因此计算域必须取得足够大，以包含整个流场，这会使计算量成倍增加。第二， 必须恰当处理由于环境流体裹入湍流流动中，而产生的流动的空间发展。大涡模拟(LEs) 的基本思路是对大尺度涡进行直接求解Nav时stokes方程．对小尺度涡通过亚格子模 型进行模拟。对于较高Rcynolds数的湍流流动．太涡模拟是一种模拟湍流流动的现实 可行方法。 基于以上背景，本文开展了平面自由湍射流拟序结构的大涡模拟研究，以期深入 理解简单射流的动力学特征，理解和掌握射流中涡结构的作用，并为以后开展两相、 燃烧过程的全三维LES及脚qs奠定坚实基础． 3撕 ———F———_『_]r1《一—— 2、控制方程及其求解 本文的模拟对象为图1所示的两个半 无限大平板下游的流动。射流 s Re2u060／v=11，300。对直角坐标系 t 下的绝熟、不可压缩平面自由湍射流的质 量连续方程和Navier．scokes方程进行滤 ● 波和无量纲化运算．得到大涡模拟的基本 圈1平面射流搬序结构大涡模拟的示意圈 控制方程：