湍流

力学的世纪难题——湍流

周恒

中国航空报

June12,2014

Abstract

人们关心流体的运动是很自然的，因为地球为大气所包围，而地球表面的2/3为水面覆盖。作为科学问题的湍流，是在1883年Reynolds做

了区分层流和湍流这两种不同形态流动的实验后确立的。而自20世纪初以

来，由于工程技术的发展，对认识湍流的规律提出了迫切的要求，从而大

大地推动了湍流的研究。在这100多年中，对湍流的认识的确取得了很大

进展，否则如航空、航天、船舶、动力、水利、化工、海洋工程等工程技

术，以及气象、海洋科学等自然科学都不可能有很大的进展。但另一方面，

人们对湍流的认识又还很不全面，从而制约了这些工程技术和自然科学的

进一步发展，也可能会对21世纪的某些新兴科学技术的形成起到制约作

用。因而在21世纪之始，再一次将这一世纪难题提到科学工作者面前是很

必要的。

1湍流运动的复杂性

湍流运动复杂性的根源在于它是强非线性系统的运动。控制湍流运动的方程:Navier-Stokes(N-S)方程是非线性的。在多数情况下，它的解是不稳定的，从而导致了流动的多次分叉，形成了复杂流态，而方程的非线性又使各种不同尺度的流动耦合起来，无法将它们分别研究。

一个世纪以来，数学家们曾对N-S方程做过大量研究，但由于其非线性带来的困难，正面的成果远不如对其他数学物理方程的研究所得到的多。看起来，进一步对N-S方程的数学性质做研究尽管重要，但依靠这一途径来解决工程技术和自然数学中提出的湍流问题恐怕是不现实的。

物理学家、力学家以及一部分数学家试图从另一途径来解决湍流，即通过直接建立能反映其某些重要特性的模型来认识湍流。例如，在20世纪

1

2对湍流研究的深入2

20年代通过和分子运动论中的分子碰撞及分子自由路程作对比，提出了涡粘性和混合长理论，解决了一些迫切需要解决的工程技术问题。这一理论虽不完整而且要依靠实验来确定其中的参数，因而远不能认为真正解决了湍流问题，但由于其简洁性，至今仍在多种场合被应用。

到20世纪30年代，由于认识到湍流中包含了很多随机的成分，所以提出了湍流的统计理论。但是与分子运动的统计理论不同，这里找不到有关湍流脉动的普适的概率密度分布。所以尽管做了很多努力，并没有能解决湍流的问题。然而，在20世纪40年代，特别是通过苏联Kolmogorov 的工作，在一些合理的假定下，在湍流泳动能谱的惯性区（当雷诺数足够大时会有这一惯性区），的确得到了一些普适的规律。有意思的是，经过50年的相对沉寂后，近来很多人对这一范围的问题再度产生兴趣，提出了不少新的模型，试图找到这一范围速度脉动的所谓标度律，即相距为l的任两点速度差（绝对值）的任一次方的平均值和l的某一次方之间的简单的指数关系。这一方面最成功的是在美华人学者佘振苏等的工作。他们首先通过均匀各向同性湍流的直接数值模拟，发现了其中存在的强间歇现象，以及一些与Kolmogorov假设不一致的现象，在此基础上，他们提出了一种级串（大尺度涡对相对小尺度涡的作用）的统计模型，从而得到了标度律，很好地与已知的实验结果吻合，克服了Kolmogorov理论的不足。这一统计模型的建立，是在用N-S方程做数值模拟的启发下提出的，但后来发现对一些不是由N-S方程控制的其他非线性系统的混沌运动，这一模型也是有效的。说明这一模型的确反映了非线性系统发生的混沌运动的某一共有的性质。其结果对非线性科学可能有比说明湍流更重要的影响，值得大家关注。

2对湍流研究的深入

但是，这一理论目前只对均匀各向同性湍流有效，而且要求雷诺数足够大而存在惯性区，其研究对象从物理机制看实际上是一封闭系统，而绝大多数与工程技术有关的湍流都是非均匀各向同性的、与周围环境有相互作用的开放系统，而且雷诺数又不是大得足以保证惯性区的存在。所以，从20世纪40年代开始，人们就探索另一条路子。对很多工程技术问题，如果能知道与流动有关的平均量就够了，而Reynolds早就推导出了平均量满足的方程（雷诺方程）。可惜其中出现了湍流脉动的二阶关联量，而由

2对湍流研究的深入3

于没有普适的脉动概率密度函数，这些关联量无法事先求出，从而成了新的未知量，使未知量数目多于方程数。

而当人们推导这些二阶关联量满足的方程时，其中又出现了脉动的三阶关联量。它们又是新的未知量，使未知量数仍然多于方程数。这一过程可以不断做下去，方程数目会越来越多和越来越复杂，但未知数总是多于方程数。这就是有名的湍流方程不封闭的问题。人们不得不停止于某阶，并设其中出现的最高阶关联量可以通过低阶的关联量表达出来，从而使方程封闭。这是我国周培源教授开创的研究方向。但是这种做法往往缺乏明确的物理图像，而通常只建立在量纲分析的基础上，因而又带来一些必须通过其他途径（比如将计算结果与实验结果比）来确定的未知参数。这使得从理论上看它很不完美，而从实践上看，这些参数又不具有普适性，对不同类型的流动，在能取得大量实验数据以检验这些参数是否合用之前，其计算结果不能令人放心。但是尽管这一方法有这些缺点，它却是目前唯一能实际提供众多工程技术问题及一些自然科学问题所需实际计算方程的途径。而且它的成果也的确已为很多问题的解决提供了有用的方法，因为这一方法在继续受到人们的关注。如果打个比方的话，近代的绝大多数药物，也都是依靠大量的经验和筛选所产生，而不是事先研究清楚了致病原因及药物治病原理而“设计”出来的。但是为了对付各种疾病，这种在一定的生理、病理、药理等理论指导下通过大量试验而筛选出有用药物的方法，在很长时间内还会是一个主要方法。同样地，大量的工程技术问题，也不能等到“彻底”弄清了湍流本质再去解决。

从20世纪50年代开始，人们已经注意到在剪切湍流中不仅有小尺度的随机脉动，而且存在有一定规律性的大尺度的拟序结构，又称相干结构。到20世纪60年代末，这一发现得到确认。而且这些大尺度结构，正是外部环境对湍流发生影响的中介。以边界层为例，外部流动在近壁面处产生强烈剪切，这一剪切会产生大尺度的结构（人们相信，也有一定证据说明这是流动不稳定的结果）。这些大尺度结构尽管还有一定随机性，但却又有相当的确定性。它们的能量占据了总湍流的大部，而又通过某种目前还不清楚的机制破裂成小尺度的脉动，从而对整个湍流的发生和维持起重要作用。从物理上说，这种剪切湍流是一种非线性开放系统，因为它与外部环境有相互作用（通过剪切力或其他途经如热量交流而产生）。在这种系统中存在着大尺度的非完全随机的运动，似乎是一种普适的规律。而这种大尺度运动的存在，可能是上述模式理论中不能找到普适参数的原因。因为小