关键词高雷诺数流动强牯性层流理论计算流体力学Nav衙‰方程组.PDF

第19卷第4期 空气动力学学报 V01．19．No4 20叭年12月 ACrA^ERODYN^^皿CASD虹CA Dec．．200I 文章编号：0258-1825(2001)04-0420-07 强粘性层流理论及在粘性流计算中的应用 高 智 (中国科学院力学研究所LHD实验室，北京100080) 摘要：对高雷诺数流动计算，为了解决部分牯性项甚至全部惯性粘性项落人误差与流动物理 尺度的关系问题，本文提出强粘性流动理论。强牯性流中至少有一个粘性项与惯性项同量 阶．理论包古了物理尺度各向相同极限、经典边界层和多层边界层理论为其特例，给出了从经 典边界层向物理尺度各向相同极限演化的尺度规律和牯性惯性诸项变化的量阶关系，阐明了 牯性与惯性力强相互作用将在剪切层的法向以及流向同时“激发”小尺度结构。对牯性流计 算，利用强枯性剪切流尺度律重新标度Ns格式的修正微分方程，给出临界网格尺度与流动 物理尺度和差分格式精度的关系，得到部分粘性项落人误差和计算结果为非物理数值牯性解 的二个判据。并以流场中的边界层、驻点和分离点邻域计算为例说明理论的应用，对强牯性 剪切流计算、证实部分粘性项甚至全部粘性惯性项落人误差的问题值得重视。 关键词：高雷诺数流动；强牯性层流理论；计算流体力学；Nav衙‰方程组 中图分类号：0357文献标识码：A 0引 言 关于高雷诺数流动的Navier-Stokes(NS)方程组计算已有大量文献，近来张涵信等【l。对此作 了评估。结论是在常见的差分格式精度和网格尺度条件下，一些或相当多的粘性项属于计算误 差，NS方程组计算其实相当于计算某种简化的NS方程组。这个结论大体也适用于湍流DNS、 于雷诺数更高，粘性项落入误差的情况可能更严重。因此关于粘性项落人计算误差、特别是计 算结果为非物理数值粘性解的起因、参量条件和解决办法等，成为一个值得研究的课题。为此 本文提出强牯性流动区域的概念，强粘性流区域外的流场中牯性作用比较次要，因此强粘性流 的理论和数值分析结果可以比较完整地解决上述课题，为越来越普及的高船数NS方程组计 算及其改进提供理论依据。 1强粘性层流理论 对高雷诺数流动，我们定义强牯性流动，强粘性流动指流场中至少有一个粘性项与惯性项 同量阶的流动区域；在强粘性流以外的流动区域中粘性项比较次要。强粘性剪切流和强粘性 扩散迎风流动是最主要的强粘性流动。强粘性剪切流的两个重要特例是经典边界层和多层 收稿日期：20014)2．19：修订日期：2001416|93 基金项目：国家自然科学基金资助项妇，中围科学院“九五”重大项日和中国科学院]MD实验审资助项f1 万方数据 第4期 高智：强粘性层流理论及在粘性流计算中的应用421 (三层)边界层的下层【3·…。在流速和坐标变量分别用来流速度U和边界尺度L归一化的条件 下，不可压二维边界层的流向和法向长度尺度却和蚱以及流向和法向速度尺度邯和邯为 ’ Re一；)m】，因此二维不可压强粘性流动的物理尺度可表示为【2] (％，虾；～，印)兰(Re一～，Re一～；Re一～，Re一～) (1) 其中船：譬，Rel，物理尺度指数毗，ny，n。和n。通过分析Ns方程中粘性和惯性诸项的 量阶关系求得。二维不可压缩Ns方程组为 舞+努=o (2a) 尺一m7 “爱+”旁=一意+瓦1、[3。a—u+3a2刊