四阶精度显式滤波大涡模拟方法研究.pdfVIP

中围T程热物理学会 热机气功热力学 学术会议论文 编号：062047 四阶精度显式滤波大涡模拟方法研究 席光辛宝玉 (西安交通人学能源与动力工程学院陕西两安710049) Emaih 联系电话：029iguang@mail．xjtu．edu．c11 摘要：在传统人涡模拟方法的基础上引入显式滤波技术，有效地控制J，最小解析尺度区域的各种数值 误差。对有限差分格式的守恒特性进行了研究。发展了适用于三维平板通道大涡模拟的四阶精度守恒 程序。分别利用二阶和四阶精度的传统与显式滤波大涡模拟程序模拟平板通道内的湍流流动，将计算 结果与DNS数据进行了对比，得到了合理结果． 关键词：大涡模拟，显式滤波，四阶精度，守恒特性 ’ 1 引言 Simulation，LES)方法通过尺度选择完成对流动尺度的分离， 大涡模拟(LargeEddy 并对各种尺度采取不同的求解策略。传统LES方法将计算网格和离散操作作为“隐式” 滤波操作运用到N．S方程，将流动尺度划分为解析尺度和亚网格尺度。最小解析尺度区 域是截止误差、混淆误差和插值误差最为严重的区域，但其中所包含的运动信息通常被 用来作为弧网格模型建立的依据。显式滤波技术的引入对流动尺度进行重新划分，将该 区域从直接求解范围中去除，可以有效地减小甚至消除各种数值误差。 实践证明，在完成对不可压缩非稳态流动的数值模拟时，有必要确保N．S方程中对 流项的离散形式对动能守恒，以便得到稳定且不受数值耗散影响的解。交错网格系统中 的标准二阶精度有限差分程序具有良好的守恒特性，但其精度太低且要求精密网格划 分。另一方面，同样采用交错网格系统，并非所有的高阶有限差分格式都对动能守恒。 2 显式滤波LES方法 显式滤波技术的运用，将流动尺度划分为解析滤波尺度与亚滤波尺度。由于显式滤 波操作的截止长度～般人于网格宽度，亚滤波尺度又被划分为解析砭滤波尺度与非解析 皿滤波尺度。解析Ⅱ滤波尺度的特征长度介于网格宽度和显式滤波截止长度之间，光滑 滤波器的运用使其可由滤波流场得到重建；非解析弧滤波尺度对解析弧滤波尺度的作用 效果必须通过模型模拟。 以平板通道高度的一半h为参考长度，初始时刻中心线上速度以为参考速度进行 、 无量纲化处理，对应雷诺数Re=uch／v。引入模型后的N-S方程变形为： 重+塑：一重+上盟一互 (1) 研 Ox ReOx Ox ,舐， j 一 式中：“～”表示“隐式”滤波操作，“一”表示显式滤波操作；湍流应力r。=堪“．一历，历．。 为统一所有项的最大波数，对西。西，和乇进行了二次显式滤波。 通常将湍流应力％拆分为两部分进行模拟：。 、 式中：At,=堪“，一蠢，厅，，称为非解析亚滤波尺度应力，反映了求解范围外小尺度的影响 178 作用；色=玩乃一厩彬，称为解析亚滤波尺度应力，反映了求解范围内朱滤波速度场与 滤波速度场间的差异。前者须通过模型模拟，后者可由滤波速度场得到重建。 3 四阶守恒有限差分格式 ． 3．1守恒条件 用以描述不可压缩流体运动的连续方程和动量方程示意性简化为如下形式： (3) 守恒条件：方程中的项可等价改写为散度形式。在满足周期性条件的区域内，连续方程、 动量方程中的压力项和粘性项项自动满足守恒条件。对流项有四种不同的常用表达形 式，散度、