《湍流的数值模拟方法进展..pdf》由会员分享，可在线阅读，更多相关《湍流的数值模拟方法进展..pdf（18页珍藏版）》请在三一文库上搜索。
1、高等计算流体力学课程作业 湍流的数值模拟方法进展 1 概述 自然环境和工程装置中的流动常常是湍流，模拟任何实际过程首先遇到的就 是湍流问题， 而湍流问题本身又是流体力学理论上的难题。对于某些简单的均匀 时均流场，如果湍流脉动是各向均匀及各向同性的，可以用经典的统计理论来分 析，但实际上的湍流往往是不均匀的，给理论分析带来了极大困难。 湍流是空间上不规则和时间上无秩序的一种非线性的流体运动，表现出非常 复杂的流动状态，主要表现在湍流流动的随机性、有旋性、统计性。传统计算流 体力学中描述湍流的基础是Navier-Stokes（N-S）方程，根据 N-S 方程中对湍流 处理尺度的不同，湍流数值模拟方
2、法主要分为：直接数值模拟（DNS） 、雷诺平 均方法（ RANS）和大涡模拟（ LES） 。 直接数值模拟可以获得湍流场的精确信息，是研究湍流机理的有效手段， 但 现有的计算资源往往难以满足对高雷诺数流动模拟的需要，从而限制了它的应用 范围。雷诺平均方法可以计算高雷诺数的复杂流动，但给出的是平均运动结果， 不能反映流场脉动的细节信息。 大涡模拟基于湍动能传输机制， 直接计算大尺度 涡的运动，小尺度涡运动对大尺度涡的影响则通过建立模型体现出来，既可以得 到比雷诺平均方法更多的诸如大尺度涡结构和性质等的动态信息，又比直接数值 模拟节省计算量，从而得到了越来越广泛的发展和应用。 2 雷诺平均方法 (
3、RANS) 雷诺平均模拟（ RANS）即应用湍流统计理论，将非定常的 N - S 方程对时 间作平均，求解工程中需要的时均量。利用湍流模式理论，对Reynolds应力做 出各种假设，即假设各种经验的和半经验的本构关系，从而使湍流的平均 Reynolds方程封闭。 2.1 控制方程 对非定常的 N - S 方程作时间演算， 并采用 Boussinesp 假设， 得到 Reynolds 方程 2 1ij iii ji jjjjj u uuuu p ufv txxxxx =0 i i u x 式中，附加应力可记为 ijij pu u，称为雷诺应力。 这种方法只计算大尺度平均流动，而所有湍流脉动对平均
4、流动的影响，体现 到雷诺应力 ij 中。由于雷诺应力在控制方程中的出现，造成了方程不封闭，为 使方程组封闭，必须建立湍流模型。 2.2 湍流模型 目前工程计算中常用的湍流模型从对模式处理的出发点不同，可以将湍流模 式理论分类成两大类：一类引入二阶脉动项的控制方程而形成二阶矩封闭模型， 或称为雷诺应力模型，另一类是基于Boussinesq 的涡粘性假设的涡粘性封闭模 式，如零方程模型，一方程模型和二方程模型。 2.2.1雷诺应力模型 雷诺应力模型 (RSM)从 Reynolds应力满足的方程出发，直接建立以 i j u u 为 因变量的偏微分方程，将方程右端未知的项（生成项，扩散项，耗散项等）用
5、 平均流动的物理量和湍流的特征尺度表示出来，并通过模化封闭。 封闭目标是雷 诺应力输运方程 : ijijji kikjkijijij kkk u uu uu u uu uu uD txxx （7） 式中 ij 是雷诺应力再分配项， ij D是雷诺应力扩散项， ij 是雷诺应力耗 散项。 典型的平均流动的变量是平均速度和平均温度的空间导数。这种模式理论， 由于保留了 Reynolds应力所满足的方程，如果模拟的好，可以较好地反映 Reynolds应力随空间和时间的变化规律， 因而可以较好地反映湍流运动规律。因 此，二阶矩模式是一种较高级的模式，但是，由于保留了Reynolds应力的方程， 加上平
6、均运动的方程整个方程组总计15 个方程，应用这样一个庞大的方程组来 解决实际工程问题，计算量很大，极大地限制了二阶矩模式的应用。 2.2.2 涡粘性模型 涡粘性模型在工程湍流问题中得到广泛应用。这是由 Boussinesq仿照分子粘 性的思路提出的，即设Reynolds应力为， , 22 () 33 ijTi jj ik kijij u uUUUk 这里 1 2 ij ku u 是湍动能， T称为涡粘性系数，这是最早提出的基准涡粘性模式，即假设 雷诺应力与平均速度应变率成线性关系，当平均速度应变率确定后，六个雷诺应力只需要通 过确定一个涡粘性系数 T就可完全确定，且涡粘性系数各向同性，可以通过
7、附加的湍流量 来模化，比如湍动能k，耗散率，比耗散率w 以及其它湍流量/k，/ 2/3 kl， kq，根据引入的湍流量的不同，可以得到不同的涡粘性模式，比如常见的k，k-w 模式，以及后来不断得到发展的k，q-w，k-l 等模式，涡粘性系数可以分别表示为 / 2 kC T k C T kC T 2 q C T lkC T 雷诺平均方法的优点为 (1) 对计算机的要求较低，同时可以得到符合工程要求的计算结果。 (2)一旦给定合理的Reynolds应力模型，可以很容易地从RANS 方程解出湍 流的统计量，所需要的计算资源小。 (3)几乎能对所有雷诺数范围的工程问题求解，并得出一些有用的结果。 其不
8、足之处在于： (1) 对不同类型的湍流， 需要采用不同的 Reynolds应力模型， 甚至对于同一 类型的问题，对应于不同的边界条件需要修改模型的常数。 (2) 由于不区分旋涡的大小和方向性，对旋涡的运动学和动力学问题考虑不 足，不能用来对流体流动的机理进行描述。 (3) 对于非定常流动、大分离流动、逆压力梯度数值模拟等问题，受湍流模 型条件的限制，很难得到满意的计算结果。 (4)严重依赖流场形状和边界条件，普适性差，计算很大程度上依赖于经验。 2.3 常用的湍流模型 常用的湍流模型有： 零方程模型： C-S模型，由 Cebeci-Smith给出；B-L 模型，由 Baldwin-Lomax
9、给出。 一方程模型： 来源由两种， 一种从经验和量纲分析出发，针对简单流动逐步 发展起来，如 Spalart-Allmaras(S-A)模型；另一种由二方程模型简化而来，如 Baldwin-Barth(B-B) 模型。 二方程模型：应用比较广泛的两方程模型有Jones与 Launder提出的标准 k-e 模型，以及 k-omega模型。 1、零方程模型 上世纪 30 年代发展的一系列湍流的半经验理论，如 Prandtl 的混合长度理论、 Taylor 的涡量输运理论、 von Karman 的相似性理论等， 本质上即是零方程湍流模 型。零方程模型直接建立雷诺应力与平均速度之间的代数关系，由于不
10、涉及代数 关系故称为零方程模型： m u u v y 其中 m称为涡粘系数， 与分子的运动粘性系数 有相同的量级。对于一般的三维 的情况，上式可写为： 2 2 3 ijmijij u vSK K 为单位质量的湍流脉动动能。为了发展上述方法，需要建立 m与平均速度之 间的关系。 1925 年，普朗特提出混合长度理论，认为存在这样的长度l，在此长 度内流体质点运动是自由的，l称为混合长度。由于湍流漩涡的作用，到达新位 置后他会低于当地周围的平均速度，此即流向脉动速度 10 ()()uU yU y， 显然， 此速度差取决于当地的平均速度梯度Uy 与微团沿y向跳动的距离l，即： U ul y 此l表示
11、在此距离内微团沿y向脉动时基本不丧失其原有速度。实际测量表明， 虽然一般情况下流向的脉动速度的均方根值大于法向值，但他们有相同的量级， 因此有： U vl y 所以有： 2 uu uvl yy 由此可算出涡粘性系数为： 2 m u l y 由此可见，若假设l不随速度变化， 则可得出湍流切应力与平均速度平方成比例， 这与实验结果是一致的。 混合长度理论已成功用于研究多种湍流剪切流，如流管、边界层和各种湍流 剪切流。 目前应用最广泛的零方程模型是Baldwim-Lomax 模型，该模型对湍流边界层 的内层和外层采用不同的混合长度假设，在流体分离不严重的流场计算中结果较 好。但是实际上，零方程湍流模
12、型仅适用于局部平衡状态的湍流流动。 2， 、一方程模型 一方程模型一般求解湍流动能或涡粘性系数的输运方程，精度较好，鲁棒性 也比较好， B-B 模型和 S-A 模型是典型的单方程模型。特别是S-A 模型，从经 验和量纲分析出发得出了涡粘性系数的输运方程，采用大量的实验结果标定模型 系数，具有良好的鲁棒性和计算准确性，目前已经被集成在各种商业软件和科学 计算的代码中，在航空航天领域的空气动力学计算中得到了十分广泛的应用。 S-A模型常被认为介于 B-L 代数模型和两方程模型之间。 由于其容错功能好， 处理复杂流动的能力强，已得到广泛应用。与B-L 模型相比，其湍流涡粘场是 连续的。且容错性好，
13、计算量少。 该湍流的原理是建立在一个附加的涡粘输运方 程的解决上。方程中包含对流项，扩散项和源项，以非守恒形式建立。S-A 模型 不同于其他一些单方程模型， 是直接根据经验和量纲分析，从简单流动开始， 直 接得到最终的控制方程。 该模型具有一些很好的特点， 相对于两方程模型计算量 小和稳定性好， 同时又有较高的精度。 由于模型方程的因变量函数在对数律区内 与到壁面的距离成线性关系， 所以可以使用相对与低雷诺数模型较粗的网格。另 外，模型是非当地型的，方程中没有诸如y+这类当地型的项在内，所以在有多 个物理面的复杂流场中不需要特殊处理，使用方便。 3、两方程模型 上世纪 70 年代，Launde
14、r发展的 k-模型被称为标准 k-模型，它求解湍流 动能 k 及湍流动能耗散率的输运方程，能够反映一定的湍流物理量的输运特性， 是两方程湍流模型的先驱性工作。之后研究人员又发展了重整化群k-(RNG k-)模型、可实现性 k-模型等，进一步强化了k-系列模型的计算性能。另外 一个系列的两方程模型为-k模型系列，其中比较有代表性的有标准-k模型 和 SST-k模型。一般来说， k-模型对高 Re 数充分发展的湍流模拟结果较好， 而-k模型改进了 k-模型对受壁面影响湍流模拟的缺陷， 对壁面附近的湍流模 拟精度较高。 （1）k-模型 k-模型是分别引入关于湍动能k 和耗散率的方程： () ()()
15、 e kkb kkkk kk u kGG txxx 12 () ()()() e kk kkk k uc Gc txxxk 其中： 222 2()2()() kt uvuv G xyyx () tt bxy tt TT Ggg xy et 2 t k C 模型中各通用常数据计算经验可取为： 12 0.09,1.44,1.92,1,1.3 k Ccc 标准 K- 模型的特点： 可用于边界层型流动和分离流； 近壁需修正或在计算边界上用壁函数（半经 验公式）作边界条件；属于涡粘模型； 方程模化不确定因素多，可靠性差；模 型常数通用性差；不能模拟强各向异性流（如矩形槽道中的二次流）；不能计入 涡量的影
16、响。 除此之外还有各种改进的k模型，比较著名的是 RNGk模型和带旋流 修正的k模型。 （2）k-模型 标准-k模型是基于 Wilcox-k模型，它是为考虑低雷诺数、可压缩性和 剪切流传播而修改的。 Wilcox-k模型预测了自由剪切流传播速率，像尾流、混 合流动、平板绕流、 圆柱绕流和放射状喷射， 因而可以应用于墙壁束缚流动和自 由剪切流动。标准k模型的一个变形是SST-k模型。 SST-k模型由 Menter 发展，以便使得在广泛的领域中可以独立于-k模 型，使得在近壁自由流中-k模型有广泛的应用范围和精度。为了达到此目的， -k模型变成了-k公式。SST-k模型和标准-k模型相似，但有以
17、下改进： （1）SST-k模型是由标准的-k模型和变形的-k模型分别乘上一个混合函 数相加得到的，在近壁面混合函数将为1，此时启用标准-k模型，在远壁面， 混合函数将为 0，此时启用变形的-k模型。 （2）SST-k模型合并了来源于方 程中的交叉扩散。（3） 湍流粘度考虑到了湍流剪应力的传播。 （4）模型常量不同。 这些改进使得SST-k模型比标准-k模型在在广泛的流动中有更高的精度和 可信性。 SST-k模型的方程为： k ()()() ikk ijj k kkuGY txxx ()()() i ijj uGYD txxx 式中： k G 由层流速度梯度而产生的湍流动能； k和 K 和的扩散
18、率； k和 K 和的扩散率； k YY和 K 和的发散项； D 正交发散项。 4、其他模型 其他形式的湍流模型涡粘系数输运（SA）模型 (3 方程)，雷诺应力模型（ 2 阶矩模型）、雷诺应力模型方程（ 7 方程模型）。 一阶矩模型在工程湍流计算中获得了很大的成功，但它们存在一些本质上的 缺陷，即这些模型均是基于Boussinesq线性各向同性的假设， 导致雷诺正应力在 三个方向上的分量相等， 这与很多实际的湍流流动矛盾。因此，一阶矩模型对强 逆压梯度下的流动、 强分离流动、 二次流、存在旋转和曲率效应的复杂湍流等预 测精度较差，需要进行相应的修正。 二阶矩模型， 即雷诺应力输运模型， 通过求解
19、雷诺应力各个分量的输运方程 来封闭雷诺应力项， 可以考虑湍流的各向异性及历史效应，理论上具有一阶矩所 不能及的模拟复杂流动的能力。 我国周培源教授首次建立了雷诺应力的输运方程 组，1951年 Rotta在这个基础上发展了完整的雷诺应力模型。他们的工作是最早 的奠基性工作。 Launder、Reece和 Rodi 对二阶矩模型进行了标定，建立了著名 的 LRR 二阶矩封闭模型。后来很多研究者又提出了多种形式的二阶矩模型。不 同二阶矩模型之间的区别在于扩散性、压力应变率关联项和耗散项的具体模化 形式，其中最关键的是压力。 应变率关联项的模化， 但到目前为止对这一项的模 化还是不成熟。 尽管二阶矩模
20、型模拟复杂湍流流动理论上具有较大的优势，但它 需要求解 6 个雷诺应力的强非线性方程及附加的湍流动能耗散率的方程，鲁棒性 较差，计算量较大， 而且实际流场中的计算精度并不不尽如意，因此在很大程度 上限制了二阶矩模型在工程中的应用。后来Rodi 提出把雷诺应力输运方程简化 为代数应力模型 (Algebraic Stress Model ，ASM)的思想。假设雷诺应力的输运正 比于湍流动能 k 的输运，带入压力 应变率关联项和湍流动能耗散率的模型，从 而得到代数应力模型。 ASM 模型不考虑雷诺应力的时间和空间导数，比较合理 地对二阶矩模型进行了简化。 介于一般意义上的一阶矩和二阶矩模型之间，另外
21、 重要的一类湍流模型即为非线性涡粘性湍流模型。尽管它的推导过程与代数应力 模型不同，但在表达形式上完全相同。Pope指出虽然非线性涡粘性模型和代数 应力模型在推导时所基于的出发点不同，但他们在数学上是等价的。 非线性涡粘 湍流模型的基本思想是改进Bousincsq假设的线性应力应变本构关系，采用非 线性的多阶表达式。 早在 20世纪 70 年代，Lumley 和 Pope就已经给出雷诺应力 的通用非线性表达形式。 非线性模型的二阶项可以反映雷诺应力的各向异性，三 阶项可以反映流线弯曲及旋转效应等。 2.4 各种湍流模型的特点 各类模型基于粗略的假设、类比、量纲分析，无可靠物理基础，因此需引进
22、经验系数。 1、零方程模型不能反映输运效应，计算量最小， 一般适用于边界层型流动， 引进各种修正可扩大适用范围； 2、K 方程模型特征长度不易确定，应用较少； 3、 方程模化不确定因素多，可靠性差； 4、标准 K 模型近壁需修正，且不能模拟二次流； 5、非线性 K 模型能反映各向异性，璧面 的仍有奇异； 6、 涡粘模型不能反应各向异性和松弛效应； 7、二阶矩模型适用范围较广，计算量较大，模型常数的通用性仍差； 8、SA（3）模型近壁无奇异性，可模拟流场变化较剧烈和曲率较大湍流， 但仍具有涡粘模型特点。 多数模型不能完全满足真实性条件，需要改进。 脉动结构信息多的模型， 应 用面较广，但模拟的对
23、象越多，不确定的因素就越多，计算量越大。现在还没有 一个模型能满意预测所有湍流，所以，选模型时应综合考虑流动类型、计算量与 精度等因素。 3 大涡模拟（ LES） 湍流大涡数值模拟（ LES）是有别于直接数值模拟和雷诺平均模式的一种数 值模拟手段。利用次网格尺度模型模拟小尺度湍流运动对大尺度湍流运动的影响 即直接数值模拟大尺度湍流运动，将 N-S 方程在一个小空间域内进行平均（或 称之为滤波），以使从流场中去掉小尺度涡，导出大涡所满足的方程。 3.1 基本思想 湍流运动是由许多尺度不同的旋涡组成的。那些大旋涡对于平均流动有比较 明显的影响，而那些小旋涡通过非线性作用对大尺度运动产生影响。大量的
24、质量、 热量、动量、能量交换是通过大涡实现的，而小涡的作用表现为耗散。流场的形 状，阻碍物的存在，对大旋涡有比较大的影响，使它具有更明显的各向异性。小 旋涡则不然， 它们有更多的共性， 更接近各向同性， 因而较易于建立有普遍意义 的模型。基于上述物理基础， LES 把包括脉动运动在内的湍流瞬时运动量通过滤 波分解成大尺度运动和小尺度运动两部分。大尺度通过数值求解运动微分方程直 接计算出来，小尺度运动对大尺度运动的影响在运动方程中表现为类似于雷诺应 力一样的应力项， 该应力称为亚格子雷诺应力， 通过建立模型来模拟。 即实现大 涡数值模拟， 首先要把小尺度脉动过滤掉， 然后再导出大尺度运动的控制方
25、程和 小尺度运动的封闭方程。 3.2 滤波函数 大涡模拟首先要流动变量划分成大尺度量和小尺度量，这一过程称之为滤 波。滤波运算相当于在一定区间内按一定条件对函数进行加权平均，其目的是滤 掉高波数而只保留低波数， 截断波数的最大波长由滤波函数的特征尺度决定。目 前较为常用的滤波函数主要有以下三种：Deardorff 的盒式 (BOX)滤波函数、富 氏截断滤波函数和高斯 (Gauss) 滤波函数。 不可压常粘性系数的湍流运动控制方程为N-S 方程: j ij ij ji i x S x P x uu t u )2( 1 式中： S 拉伸率张量，表达式为：2/)/( ijjiij xuxuS；分子粘
26、性 系数；流体密度。设将变量 i u分解为方程 (11)中 i u和次网格变量 (模化变量 ) i u ， 即 iii uuu， i u可以采用 Leonard 提出的算式表示为 : （11） 式中)(xxG称为过滤函数，显然G(x)满足 3.3 控制方程 将过滤函数作用与N-S 方程的各项，得到过滤后的湍流控制方程组： 由于无法同时求解出变量 i u和 jiu u，所以将 jiu u分解成 ijijij u uu u， ij 即称为次网格剪切应力张量(亦称为亚格子应力 )。 由此动量方程又可写成： xdxuxxGxu ii )()()( 1)(dxxG ()(2) 1ijij i jij u
27、 uS uP txxx ()(2) 1ijijij i jijj uuS uP txxxx 式中 ij 代表了小涡对大涡的影响。 3.4 常用的亚格子模型 目前， 在大涡模拟中经常广泛采用的亚格子模型有标准的Smagorinsky模型、 动态涡粘性模型、动态混合模型、尺度相似模型、梯度模型、选择函数模型等。 其中 Smagorinsky模型被广泛应用。 3.4.1 亚格子涡粘和涡扩散模型 不可压缩湍流的亚格子涡粘和涡扩散模型采用分子粘性和分子热扩散形式，即 kkijij t ij S 3 1 2 i ti x T 以 上 公 式 中 t 和 t 分 别 称 为 亚 格 子 涡 粘 系 数 和
28、亚 格 子 涡 扩 散 系 数 ； (1/ 2) (/)(/) ijijji Suxux是可接尺度的变形率张量。式( 14)第 2 项是为了满足 不可压缩的连续方程，当 ij S收缩时 ( ij S=0) 等式两边可以相等。 将亚格子应力的涡粘模型公式( 14) 代入到 ( 13) 式中，变形得 )() 3 ( i j i i t i kk ii i j i x u x u x p xx u u t u 0 i x ui 3.4.2 Smagorinsky模型 Smagorinsky模型是由 Smagorinsky于 1963 年提出来的， 该模型是第一个亚 格子模型。广泛用于大涡模拟中的涡粘
29、模型认为亚格子应力的表达式如下： ijTkkijij S2 3 1 式中 (1/ 2) (/)(/) ijijji Suxux 是可接尺度的变形率张量， T是涡粘 系数。 1963年 Smagorinsky定义了涡粘系数： ) 式中 2/1 )2( ijij SSS是变形率张量的大小，是过滤尺度， CS无量纲参数， 称为 Smagorinsky系数。 3.4.3 动态亚格子模式 1991 年，Germano提出了动态亚格子模式，该模式以Smagorinsky 模式 为基本模型， 但克服了 Smagorinsky 模式的部分缺陷。 动力模型实际上是动态确 定亚格子涡粘模型的系数。动力模型需要对湍
30、流场做两次过滤，一次是细过滤， 细过滤后再做一次粗过滤。 通过在网格尺度和检验滤波器尺度条件下计算得到的 应力差来确定应力模型系数， 使模型系数成为空间和时间的函数，从而避免了在 模拟过程中对系数进行调节。 因此比 Smagorinsky 模式所采用的固定系数值更加 合理。 3.4.4 相似性模式 1980 年 Bardina 提出了尺度相似模式。该模式假定从大尺度脉动到小尺度 脉动的动量输运主要由大尺度脉动中的最小尺度脉动来产生，并且过滤后的最小 尺度脉度速度和过滤掉的小尺度脉动速度相似。通过二次过滤和相似性假定可以 导出亚格子应力表达式。 采用这种模式能正确预测墙壁面附近的渐近特性，但预
31、测各向不均匀的室内空气复杂流动准确性较差。 3.4.5 混合模式 混合模式是将尺度相似模式和Smagorinsky 模式叠加来确定亚格子应力。 这 种模式既有和实际亚格子应力良好的相关性，又有足够的湍动能耗散。 SC ST 2 )( 3.5 大涡模拟的特点 (1)能够描述小尺度湍流流动，但是计算量远小于DNS，在科学研究和工程 应用上都显示出良好的发展前景。 (2)用非均匀网格能够使网格数达到最少，节省计算资源，同时又能够保证 足够的计算精度。 (3)网格尺度比湍流尺度大，可以模拟湍流发展过程的一些细节。 (4)相较于 RANS 方法，LES 可以模拟更多的湍流大尺度运动，LES 所用的 湍流
32、亚网格应力模型受边界的几何形状和流动类别的影响小，比RANS 方法所 用的 Reynolds应力更具普适性。 其不足之处在于： (1)小涡模型网格节点的划分极密集，需要庞大的计算机存储能力; (2) 大量数据处理和非线性偏微分方程的求解需要高速数值处理能力; (3)仅用于比较简单的剪切流运动及管流。 (4) 由于实际湍流极其复杂，数值模拟仍需要非常可观的计算时间和实验 经费。 4 直接数值模拟(DNS) 湍流直接数值模拟 (DNS)就是不用任何湍流模型，直接求解完整的三维非定 常的 N - S 方程组，计算包括脉动在内的湍流所有瞬时运动量在三维流场中的时 间演变。 4.1 控制方程 用非定常的
33、 N - S 方程对湍流进行直接计算，控制方程以张量形式给出 : 2 1 iii ji jiji uuup ufv txxxx =0 i j u x 4.2 常用数值方法 由于最小尺度的涡在时间与空间上都变化很快，为能模拟湍流中的小尺度结 构，具有非常高精度的数值方法是必不可少的。 4.2.1 谱方法或伪谱方法 谱方法或伪谱方法是目前直接数值模拟用得最多的方法，其主要思路为， 将 所有未知函数在空间上用特征函数展开，成为以下形式： 123 (, ) m n pmnp mnp Vx tatxxx（3） 其中 m，n与p，都是已知的正交完备的特征函数族。在具有周期性或统 计均匀性的空间方向一般都采
34、用Fourier 级数展开，这是精度与效率最高的特征 函数族。在其它情形，较多选用Chebyshev多项式展开，它实质上是在非均匀网 格上的 Fourier 展开。此外，也有用Legendre， Jacobi， Hermite 或 Laguerre等 函数展开，但它们无快速变换算法可用。如将上述展开式代入N-S 方程组，就 得到一组 mnp at 所满足的常微分方程组，对时间的微分可用通常的有限差分法 求解。 在用谱方法计算非线性项例如V的 Fourier 系数时，常用伪谱法代替直 接求卷积。伪谱法实质上是谱方法与配置法的结合，具体做法是先将两量用 Fourier 反变换回到物理空间，再在物理
35、空间离散的配置点上计算两量的乘积， 最后又通过离散 Fourier 变换回到谱空间。在有了快速Fourier 变换(FFT)算法以 后，伪谱法的计算速度高于直接求两Fourier 级数的卷积。但出现的新间题是存 在“混淆误差”，即在做两个量的卷积计算时会将本应落在截断范围以外的高波 数分量混进来， 引起数值误差。 严重时可使整个计算不正确甚至不稳定，但在多 数情形下并不严重， 且有一些标准的办法可用来减少混淆误差，但这将使计算工 作量增加。 4.2.2高阶有限差分法 高阶有限差分法的基本思想是利用离散点上函数值 i f 的线性组合来逼近离 散点上的导数值。设 i F 为函数 j fx 的差分逼
36、近式，则 jjj Ff 式中系数 j 由差分逼近式的精度确定，将导数的逼近式代入控制流动的N - S 方程，就得到流动数值模拟的差分方程。差分离散方程必须满足相容性和稳 定性。 4.3 直接数值模拟的特点 (1)直接数值求解 N-S 方程组，不需要任何湍流模型，因此不包含任何人为 假设或经验常数。 (2)由于直接对 N - S 方程模拟，故不存在封闭性问题，原则上可以求解所有 湍流问题。 (3)能提供每一瞬时三维流场内任何物理量（如速度和压力）的时间和空间 演变过程，其中包括许多迄今还无法用实验测量的量。 (4)采用数量巨大的计算网格和高精度流体力学计算方法，完全模拟湍流流 场中从最大尺度到最
37、小尺度的流动结构，描写湍流中各种尺度的涡结构的时间演 变，辅以计算机图形显示，可获得湍流结构的清晰与生动的流动显示。 DNS 的主要不足之处在于：要求用非常大的计算机内存容量与机时耗费。 据 Kim ，Moin &Moser 研究，即使模拟 Re 仅为 3300 的槽流，所用的网点数 N 就约达到了 6 2 10 ，在向量计算机上进行了250 h。 5 不同方法的计算效率比较 LES， DNS，RANS 三种方法中 DNS 的计算量最大， LES 的计算量介于另 外两者之间，而 RANS 的计算量最小。影响计算量的因素有三个：网格数量、 流场的时间积分长度 （与计算时间长度有关） 和最小旋涡的
38、时间积分长度 （与时 间步长有关），其中网格数量是重要因素。 直接数值模拟为了得到湍流问题足够精确的解，要求能够数值求解所有旋涡 的运动，因此要求网格的尺度和最小旋涡的尺度相当，即使采用子域技术， 其网 格规模也是巨大的。 为了求解各个尺度旋涡的运动， 要求每个方向上网格节点的 数量与 3 4 Re 成比例，考虑一个三维问题，网格节点的数量与 9 4 Re 成比例。一般 的估计如下 :湍流中包含许多尺度不同的涡，为能模拟最小涡的运动，计算网格 的分辨率应足以分辨最小尺度的涡， 后者以 Kolmogorov 定义的内尺度 3 =v（） 为代表。而计算区域的尺寸应足以容纳最大尺度的涡，最大涡的尺度
39、为L。因此 在一个空间方向上的网点数目至少应与 L 同量阶，而根据统计理论知道这个比 值 3 43 2 L L RR或 于是整个三维空间所需的网点总数至少为 3 9 9 4 2 L L NRR或 此数字也正是按非线性动力系统理论所估计的湍流的吸引子维数的上确界。 计算所需的内存容量应与此数成正比。另一方面计算的时间步长应小于最小涡的 时间尺度 u，而总的计算时间应大于最大涡的特征时间 L u，因此需要计算 的步数应不少于 L.如假设每一时间步长的计算工作量，即使按最低限估计， 与 N 成正比，则总的计算工作量至少也要正比于 3 L R或 6 R 。假如对每一时间步的每 一网点需执行 100 条
40、机器指令，则对一个 5 10 L R的湍流问题，就需执行总共约 17 10条指令。这意味着在一个计算速度为每秒一亿次的超级计算机上也要运行约 30 年。如此巨大的计算工作量即使对当今世界上最大的计算机也是不可接受的。 据 Kim ，Moin &Moser 研 3 究，即使模拟 Re 仅为 3300 的槽流，所用的网点数 N 就约达到了 6 2 10，在向量计算机上进行了250 h。在现有的计算机能力的限 制下，即使在少数拥有世界最大的超级计算机的科学大国，目前也只能计算中等 以下雷诺数且有简单几何边界的湍流流动。 大涡模拟与直接数值模拟相比节省很大的计算量。湍流大涡数值模拟将湍流 的大尺度脉动
41、和小尺度脉动分开，对大尺度结构进行直接数值模拟，通过建立亚 格子尺度 (亚格子尺度 ) 模型来模拟小尺度脉动的作用。理想的湍流直接数值模 拟需要包含所有尺度的湍流脉动，一般最小的脉动尺度等于Kolmogorov 耗散尺 度，流动的最大尺度L 由流动的几何条件确定。直接数值模拟的一维网格数 应为: DNS NL， 而大涡数值模拟的一维网格数为: LES NL可以节省网格 数 3333 1 DNSLESDNS NNN，如果过滤尺度等于2 倍柯氏耗散尺度 的话，就可以比 DNS 节省 87. 5 %的网格。这里我们可以看到完全的湍流直接数 值模拟中， 绝大部分的计算量花费在耗散尺度中，对于高雷诺数流动， 这是很不 经济的计 4 算。 雷诺时均方程法先将湍流中的物理量如速度、浓度等分成扰动量及平均量， 再利用对控制方程作时间平均， 同时采用湍流模型仿真湍流的效应，因此大大降 低了计算量，但其结果受湍流模型的影响很大。 结束语 湍流是一个世纪性难题， 近一个世纪以来对其的研究不断深入。在此过程中 对湍流的模拟方法也得到了很大发展。对湍流数值模拟的研究已经成为计算流体 力学中一个热点和前沿问题。 在未来，对这方面的研究还将不断发展。对湍流复 杂流动现象会有更深刻的认识。