湍流的数值模拟综述.pdf

湍湍流流的的数数值值模模拟拟综综述述

湍流的数值模拟

⼀、引语

流体的流动形态分为湍流与层流。⽽层流是流体的最简单的⼀种流动状态。流体在管流动时，其质点沿着与管轴平⾏的⽅向

作平滑直线运动。此种流动称为层流或滞流，亦有称为直线流动的。流体的流速在管中⼼处最⼤，其近壁处最⼩。管流体的

平均流速与最⼤流速之⽐等于0.5，根据雷诺实验，当雷诺准数引Re2320时，流体的流动状态为层流。当雷诺数Re2320

时，流体流动状态开始向湍流态转变，湍流是⼀种很复杂的流动状态，是流体⼒学中公认的难题。

⾃从19世纪末O.Reynolds提出湍流的统计理论以来，已经有⼀个多世纪了，经过⼏代科学家的努⼒，湍流研究取得很⼤进

展，但是仍然不能满⾜⼯程应⽤的需要，以⾄于经常有悲观的论调侵袭湍流研究。为什么湍流问题没有圆满地解决会受到如此

关注呢?因为湍流是⾃然界和⼯程中⼗分普遍的流动现豫，对于湍流问题的正确认识和模化直接影响到对⾃然环境的预测和⼯

程的质量。例如，当前影响航天器⽓动⼒和⽓动热预测准确度的主要障碍是缺乏可靠的湍流模型。和其他⼀些⾃然科学的准题

不同，解决湍流问题具有迫切性。

湍流运动的最主要特征是不规则性，这是⼤家公认的。对于湍流不规则性的深⼊认识，是⼀百多年来湍流研究的上要成就之

⼀。早期的科学家认为，像分⼦运动⼀样，湍流是完全不规则运动。类似于分⼦运动产⽣黏性，湍流的耗散可以⽤涡黏系数来

表述。20世纪初，⼀些杰出的流体⼒学家，相继对涡黏系数提出各种流体⼒学的模型，如Taylor(1921年)的涡模

型，Praudtl(1925年)的混合长模型和vonKarman(1930年)相似模型等。当科学家⽤流体⼒学观念(不是分⼦观念)来建⽴湍流耗

散的涡黏模型时，就开始考虑连续介质不规则运动的特点，其中有别于⽓体分⼦不规则运动的最主要特点是运动的多尺度性。

第⼀个提出流体湍流运动中多尺度输运特性的科学家mhardson(1922年)曾描述湍动能的多尺度传输过程如下：“⼤涡包含⼩

涡，并喂予速度；⼩涡包含更⼩的涡，如此继续直到黏性耗散”。多尺度的思想导致产⽣描述多尺度的谱概念和谱分析⽅法，并

最终产⽣了Kolmogorov(1941年)的局部各向同性的通⽤谱(即5/3谱)。

湍流不仅是多尺度的⽽且是有结构的运动。20世纪中叶，⼤量的湍流实验(包括测量和显⽰)发现多尺度的湍流运动存在某种特

殊的运动状态。Townsend(1951年)，Corrsin(1955年)和Lumley(1965年)等从脉动序列的间歇性和空间相关相继推测湍流结构

的可能形态。理论上也提出过各种湍涡的模型：球涡模型，柱涡模型等。早期的湍流结构主要是从运动学上考虑，把旋涡结构

作为湍流统计的样本。我国的周培源教授是近代湍流模式的奠基⼈之⼀，他⾸先提出先解⽅程后平均的统计⽅法，就是说湍涡

必须满⾜Navier—Stokes⽅程(ChouandChou，1995年)。

真实的、可以观察到的湍流结构通过流动显⽰，以及稍后湍流直接数值模拟所证实。典型的例⼦是混合层的Brown—Roshko涡

(1976年)，图1明显地展⽰了混合层中存在规则的⼤涡和分布在⼤涡周围的细⼩湍涡。在边界层、槽道和圆管湍流中也存在各

式各样的⼤涡结构。例如，⽤激光诱导荧光的显⽰⽅法，我们可以在圆管湍流中观察到周向（图2a）和流向⼤涡（图2b）。值

得提出的是，不仅在剪切湍流中有⼤涡结构，简单的均匀各向同性湍流中也存在涡结构。图3展⽰的是各向同性湍流的直接数

值模拟中强涡量等值⾯，它们是管状结构。仔细分析还可以确定管状涡的平均长度约等于各向同性湍流的积分尺度，它们的平

均直径约等于湍流TayLor微尺度，更进⼀步分析可以算出管状涡部的平均速度

场，它们接近于Burgers涡，即有轴向拉伸的柱状涡，在管状涡之间错综复杂地分布着各种尺度的树叉结构。

所有以上发现充分说明：⽆论是简单还是复杂湍流，都存在⼀定的涡结构．⼤尺度结构的发⽣是不规则的，就是说，在长时间

和⼤范围来观察，⼤尺度运动结构发⽣的地点和时划是不确定的．因此在⼤样本统计中我们不可能发现这种结构，这就是为什

么经典的长时间统计未能察觉它们的原因。另⼀⽅⾯，⼤尺度运动结构⼀旦⽣成，它以⼀定的动⼒学规律演化，因此湍流⼤尺

度结构⼜称拟序结构，或相⼲结构。举例来说，在湍流边界层、槽道或圆管湍流的近壁区(5

湍流中⼤涡拟序结构对于湍流⽣成和发展有主宰作⽤，因此抑制或消除⼤涡结构可能抑制整体的湍流强度，甚⾄使流动层流

化。这是近代湍流减阻和降噪的思想(Bushnell等，1989)。湍流是多尺度有结构的不规则流体运动．它指出湍流运动的主要特

征，同时也指出了研