平板边界层转捩的数值模拟.pdfVIP

第 6卷 第 2期 航空动力学报 V01．6 No．2 1991焦 4月 JournalofAero6pacePower Apr． l99l 平板边界层转捩的数值模拟 中国科学院工程热物理研究所 蒋洪德 美国弗吉尼亚理工学院机械系 J．G．Moore．J．Moore 【擒要】 对于 自由稿魂弓f起的压力梯度为零的平板边界盛转捷过程的数值模拟，本文提出了一 十单方程jjlf漉模量 fVL^I模型)．对jjlf漉牿性系数和精漉长度尺度靠进行 了謦正以正确引^低雷诺 数效应 ．数值计算表明．该模型在 自由来观湍流度为0．O3～8 的大艳茸内能正确模拟边界盛转捩 全过程 ．包括转捩区起点位置和转捩区长度，边界层 内的平均及脉动参数分布等．本文还提出了上 述转捩过程的时均化五阶段模型． 一 、 前 言 对层流到湍流转捩的深入理解并准确预测仍是当代流体力学面临的难题之一。近年来用 湍流模型来模拟转捩过程的方法 日益受到重视。但通常使用的湍瀛模型是对充分发展的湍流 建立的t因而都不能正确模拟转捩过程。例如 目前广泛使用的k-~模型用于 自由湍流引起的平 板边界层转捩预测时，计算结果与实验数据有很大差别。本文对单方程湍流模型进行了研究， 在综合实验数据的基础上，对湍流粘性系数和湍流长度尺度都进行 了修正。数值计算表明，本 文提出的单方程湍流模型在模拟 自由湍流引起的边界层转按时具有 良好的能力。 二、充分发展湍流的单方程湍流模型 单方程湍流模型使用湍流动能方程 ，它包括端流动能k的对流、扩散、生成和耗散四部 份。日前常用Hassld和 Poreh于1975年 提出的下列方程 ： 孤／知-t-pratt匆一 a[( -t- )孤／匆]／9 十 (a【，／两) 一 [ ·CIC,+ 2~ k／L (1) 在高雷诺数和各向同性湍流的假定下，赭瀛粘性系数 和湍瀛长度尺度L可分别表示为： “= 口 L (2) f坤[1一exp(Ap~ m)] 内层 … I(L )一 ·d 外层 方程 (1)～ (3)构成了充分发展湍流的单方程湍瀛模型 (简称HP模型)，其中的常数 = 0．17， =0．536， 一0．0292， =0．l和 )～ =0．08。图l是 HP模型用于压力梯度为零 的平板边界层转捩的计算结果以及与实验数据 的比较 。可 以看 出计算求得的转捩起始位置 均大大超前于实验结果 、转捩区长度远小于实验值，且转捩后壁面摩擦系数明显低于充分发 展湍流状态值。图2则表明HP模型的计算结果强烈地依赖于平板前缘的初始雷诺数 。 本文于Z990~ 0yJ收到． 58 航 空 动 力 学 报 第 6 卷 图l 来流j褂流度变化时 胛 模型的 图2 计算韧始雷诺敢变化时胛 模型 计算与实验值的比较 的计算与实验值 的比较 总之，高雷诺数单方程湍流模型不能直接用于转捩模拟。 三、湍流粘性一长度尺度修正 (VLM)模型 单方程湍流模型中主要变量是湍流粘性系数 n和湍流长度尺度 L。对 已有实验数据。进 行综合分析后，本文对n和L进行如下修正 ： 1．使用一个分段光滑的函数 修正n，即 = ·C~plc