运动壁面槽道流动的直接数值模拟-应用数学和力学.pdf

摇 , 31 1 Applied Mathematics and Mechanics摇 摇 2010年1月15 日出版 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 Vol.31,No.1,Jan.15,2010 文章编号:1000鄄0887(2010)01鄄0091鄄11 訫 应用数学和力学编委会,ISSN 1000鄄0887 运动壁面槽道流动的直接数值模拟* 葛铭纬,摇 许春晓,摇 崔桂香 (清华大学 航天航空学院,北京 100084) (林建忠推荐) 摘要:摇 采用谱方法,在曲线坐标系下对不可压缩Newton流体的N鄄S方程进行求解,采用定义在物 理空间中的流动物理量以避免使用协变、逆变形式的控制方程.在计算空间采用Fourier鄄Chebyshev 谱方法进行空间离散,时间推进采用高精度时间分裂法.为了减小时间分裂带来的误差,采用了高 精度的压力边界条件.与其他求解协变、逆变形式控制方程的谱方法相比,该方法在保持谱精度的 同时减小了计算量.首先通过静止波形壁面和行波壁面槽道湍流的直接数值模拟,对该数值方法 进行了验证;其次,作为初步应用,利用该方法研究了槽道湍流中周期振动凹坑所产生的流动结 构. 关摇 键摇 词:摇 谱方法;摇 运动壁面;摇 槽道;摇 湍流 中图分类号:摇 O357.5摇 摇 摇 文献标志码:摇 A DOI:10.3879/ j.issn.1000鄄0887.2010.01.010 引摇 摇 言 湍流的直接数值模拟是湍流基础研究的重要手段,由于计算机的限制,湍流的直接数值模 拟大多是对中低Reynolds数下简单边界的湍流流动.迄今为止,人们通过直接数值模拟获得 [1] 了许多对湍流机理的重要认识,而由于实验条件的限制,有些发现很难利用实验得到 .对于 壁湍流,大多数研究都是针对静止壁面的,近年来,复杂的运动壁面在湍流控制的应用中逐渐 [2] [3] 得到了人们的重视,例如智能蒙皮 、柔性壁面 在湍流减阻中的应用等.但是,目前关于复 杂运动壁面湍流的研究还非常缺少. 谱方法具有高效和高精度的特点,因此被广泛应用于典型湍流的直接数值模拟,例如: [4] Fourier鄄Galerkin方法是各向同性湍流直接接数值模拟的常规方法 ,Fourier鄄Chebyshev方法在 [5] [6] 槽道湍流直接数值模拟中发挥了重要作用 .对于复杂运动壁面槽道中的流动,Carlson等 在正交坐标系下利用伪谱方法对协变、逆变形式的控制方程进行了求解.但是,Carlson等在坐 [7] 标变换的方程中遗漏了一些重要的附加项,Luo等 指出了该错误并通过张量运算在随时间 变化的曲线坐标系下得到了N鄄S方程完整的协变、逆变形式.尽管协变、逆变形式的N鄄S方程 对于谱方法的应用非常方便,但是因为协变、逆变所产生的大量附加项使得计算量大大增加, * 收稿日期:摇 2009鄄05鄄13;修订日期:摇 2009鄄11鄄26 基金项目:摇 国家自然科学基金资助项目 作者简介:摇 葛铭纬(1984—),男,山东人,博士生(E鄄mail:gmw06@mails.tsinghua.edu.cn); 许春晓,教授,博士(联系人.Tel:+86鄄10E鄄mail:xucx@tsinghua.edu.cn). 9