众所周知，湍流是著名的物理学世纪难题。著名物理学家诺贝尔奖获得者费曼说过,湍流是至今没有得到解决的经典物理学的最后一个问题。自从Navier-Stokes方程（1821-1845）建立以来的200年间，自从Reynolds进行的著名的层流-湍流转捩实验（1883）以来的140年来，众多位科学先辈做了大量的理论、实验和计算研究工作，对于湍流是怎么产生的问题都百思不得其解。在湍流研究方面做出巨大贡献和努力的著名科学家包括普朗特、泰勒、von Karman、Kolmogorov，周培源、林家翘、John von Neumann等人，还包括著名诺贝尔奖获得者Rayleigh、朗道、海森堡、费曼、Chandrasekhar、Onsager等人。在二十世纪，经典物理学被得以终结，相对论和量子力学被得到普遍理论认可和已经获得广泛实际应用的情况下，湍流问题还没有得到解决，湍流问题真的是太难了。

偏偏有人不怕难，为了彻底解决湍流问题，窦华书经过30多年的艰苦努力，建立了能量梯度理论，针对不可压缩流体，进一步发展了关于流体力学中的湍流及流动稳定性的公理、推论和定理，形成了一个崭新的关于湍流及流动稳定性的理论体系，开创了一个新的学派。根据能量梯度理论和流体流动的动量方程即Navier-Stokes方程，窦华书在其发表的学术论文及专著中给出了若干公理、推论和定理，见下面 [1-8]。这些科学定律得出的条件是流动受到有限扰动，而不是线性扰动；这里面的主要区别是有限扰动能够改变流场中的机械能分布并诱发奇点，而线性扰动不能，这是导致湍流产生的流动之所以失稳的关键机理。这些流体力学的公理、推论和定理，在其他教科书上是没有的，是窦华书首次建立的 [1]。正是基于这些研究基础和结论，窦华书在国际上首次成功地解决了湍流的百年难题（湍流是怎么产生的），发现了湍流是由Navier-Stokes方程的奇异性所引起（图1-5）。2000年，美国Clay数学所宣布了7个千禧年大奖难题，Navier-Stokes方程的存在性及光滑性问题为其中之一，窦华书分别采用能量梯度理论和泊松方程分析方法，在国际上首次给出了正确的答案，即对转捩流动和湍流流动，Navier-Stokes方程不存在全局定义域上的光滑解 [1-8]。

窦华书的理论与所有能够得到的大量实验数据获得了一致性；推翻了过去100多年间前人得出的许多有关湍流及流动稳定性的已有结论。原有的那些定理、结论，互相矛盾，与实验不一致，没有能够解决湍流问题。湍流问题的真正解决必然推翻那些互相不一致的定理和结论。作者认为，一个正确的物理学的准则，必然既是必要条件也是充分条件。一个所谓的必要条件，还有进一步的必要条件，另外的必要条件（没有充分条件），都不能拿来作为物理学的准则，可是这些被写在了教科书中，并被用来分析湍流的发生机理。作者认为，湍流问题没有得到解决与这些误导是有很大关系的。一个没有充分条件的必要条件是没有多大意义的。但是，过去前人走过的弯路可以理解，毕竟在黑暗中摸索非常艰难，同时前人为后人趟了一次雷，让我们知道了那些路走不通。

通过作者的研究还发现，Navier-Stokes 方程与Euler方程是性质完全不同的方程，湍流是Navier-Stokes方程控制的，由非线性和粘性共同相互作用导致的复杂流动，缺一不可。在研究过程中发现，研究Euler方程对研究湍流没有看到有任何帮助（反而有时候会有误导），这在发现了湍流生成的机理后，对此理解更加深刻。湍流的产生绝对不是“NS方程中对流非线性项的作用超过了粘性耗散性的影响而产生的”；这句话是过去100多年最具有误导性的一句话，它出现在几乎所有的相关的教科书及学术文献中。

公理是自然界能够被普遍认可的科学规律，是不需要证明的。定理是根据已有的科学公理和定理进行证明的科学规律。推论是根据已有的公理和定理进一步逻辑推断得到的结论。下面是窦华书建立的公理、推论和定理 [1]:

1.公理 (Axiom)

公理5.1 在实验室系统中，如果总机械能沿流动方向的梯度小于零，则流动是可能的。

公理5.2 在实验室系统中，如果沿着垂直于流线的方向的总机械能梯度为零，则流动是稳定的，否则可能不稳定。

公理5.3 在实验室系统中，如果沿着流线的总机械能梯度为零，则流速为零。

公理5.4 在实验室系统中，在剪切驱动的平行流中，如果对流体单元上的剪切应力做功为零，则流动将立刻停止，速度为零。

上面的公理5.1，5.2，5.3适用于压力驱动流动（pipe flow, channel flow, Dean flow, etc）；公理5.4适用于剪切驱动流动(plane Couette flow, etc)。这两种情况的区别就是看看有没有外部对流体做功。因为能量和做功都不是伽利略不变量，参照系选择在实验室系统。

2. 推论 (Corollary)

推论5.1 湍流中猝发（burst）的发生伴随着压力峰值。

推论8.1 粘性在平行流中仅起稳定作用（粘度越大越稳定）。

推论8.2 粘性在圆周流动中仅起稳定作用（粘度越大越稳定）。

推论8.3 势流中不会产生湍流。

3. 定理 (Theorem)

定理4.1 当系统没有功输入或输出时，具有拐点的速度分布是不稳定的。

定理5.1 速度剖面上的拐点是Navier-Stokes方程的奇点，此处速度发生间断。

定理5.2 在高雷诺数（对转捩流动和湍流流动）下，Navier-Stokes方程不存在光滑的且物理上合理的解。（此定理在数学上精确地回答了千禧年大奖难题）。

定理5.3 无粘流动中不会产生湍流 （是因为不会产生奇点）。

定理5.4 在剪切驱动流（平行流）中，速度剖面上存在Navier-Stokes方程的奇点，在该奇点处速度梯度为零，且速度发生间断。

定理6.1 在剪切驱动流（非平行流）中，速度剖面上存在Navier-Stokes方程的奇点，在该奇点处 ，且速度发生间断。 

定理7.1 湍流转捩/湍流产生的充分必要条件是流场中能量梯度函数K (K=法向能量梯度/流向能量损失)为无穷大。

定理7.2 对于压力驱动流动，湍流转捩的充分必要条件是速度剖面上存在速度拐点。

定理7.3 对于剪切驱动流动，湍流转捩的充分必要条件是速度剖面上存在零速度梯度。

定理7.4 湍流产生/湍流转捩的充分必要条件是流场中出现Navier-Stokes方程的奇点（速度发生间断）(此定理在物理学上成功地解决了百年湍流难题)。

定理8.1 势流（无粘和无旋）是稳定的。

定理8.2 无粘有旋流动是不稳定的。

定理8.3 自由涡是稳定的。

定理8.4 强制涡是不稳定的。

定理7.1的证明:

充分性：如果K无穷大，则K的分母为零，即没有能量消耗，则u=0，发生间断，湍流产生。

必要性：如果湍流产生，则发生间断u=0，即能量消耗为零，分母为零，则K无穷大。

定理7.4的证明:

充分性：如果速度发生间断u=0，有了奇点，则会产生涨落、产生新涡量，即湍流产生。

必要性：如果湍流产生，则必须有新涡量，必须有涨落，必须间断u=0，发生奇点。

定理7.1和7.4是等价的；其他定理的证明参见[1]。定理4.1和5.1适用于压力驱动流动。定理5.2是在压力驱动流动条件下得到的，此定理也适用于剪切驱动流动。湍流只能产生在粘性流体中，是因为粘性的作用能够制造奇点。没有粘性就没有奇点，所以没有粘性，就没有湍流。如果要更好地理解这些公理、推论和定理，请参考专著[1]的上下文。奇点的理论速度为零，与驻点速度为零是不一样的，驻点是沿流线速度逐渐变化到u=0。奇点是沿流线速度突然断崖式变化到u=0。

能够敢于在国际同行和权威面前，公开主张自己的这些研究成果，对科学发现的可靠性、正确性以及创新性，自然是有100%的自信和底气的，毕竟作者做的是从0到1的原始创新工作（比如，如果自己都搞不准的东西，怎么会可能说服别人来相信呢？）。从2022年5月开始，作者在科学网和知乎网以及微信公众号（流体空间、风流知音等，声振之家、飞机设计视界，也有转载）精心写了20多篇的长篇科普文章，就是为了多和同行交流，促进湍流科学的研究进展。真诚欢迎各学科的同行，从学术方面提出意见，批评指正。

参考文献

1. Dou, H.-S., Origin of Turbulence-Energy Gradient Theory, 2022, Springer. https://link.springer.com/book/10.1007/978-981-19-0087-7 (全书下载地址).

2. Dou, H.-S., Energy Gradient Theory of Hydrodynamic Instability, The Third International Conference on Nonlinear Science, Singapore, 30 June-2 July, 2004. 链接如下： https://arxiv.org/abs/nlin/0501049

3.科学网-海森堡的第二个问题终于有了答案，窦华书的博文。

https://blog.sciencenet.cn/blog-3057857-1361491.html

4.新书访谈，专访《湍流的起源—能量梯度理论》作者窦华书教授。https://mp.weixin.qq.com/s/unFxknnohDUwp11150OfZw

5.窦华书教授成功破解了百年湍流难题，中国教育日报网。http://chinaedutech.com/dfjy/2022/1117/1327.html

6.窦华书教授在纳维-斯托克斯方程问题上取得新进展，浙江理工大学官网新闻。https://news.zstu.edu.cn/info/1033/41169.htm

7. 湍流是怎样产生的? 最新研究进展！https://blog.sciencenet.cn/blog-3057857-1341235.html

8.千禧年大奖难题之一纳维-斯托克斯方程的解的存在性与光滑性的证明,

https://blog.sciencenet.cn/blog-3057857-1337452.html

图1 边界层转捩过程中湍流斑的形成（Gad-El-Hak et al. 1981）。实验表明湍流斑是局部湍流区，它由局部低速区（奇点）导致的湍流猝发所产生。

图2 边界层转捩中的流动可视化。左边 Experiment of Hama & Nutant (1963); 右边 simulation of Zang & Hussaini (1987)。奇点诱导的猝发清晰可见。猝发引起了流体内部旋涡的产生。

图3 边界层湍流中的速度亏损（奇点）导致的猝发现象（Kline et al. 1967, experiment）。

图4 平面Couette流动中的湍流猝发和转捩（DNS结果，奇点引起），t=10,20,140,160 （Cherubini and De Palma，2013）

图5 平面Poiseuille流动中由奇点导致的湍流转捩（LES结果，Schlatter et al 2006）。Colour plot of the streamwise velocity u for the one-, three- and five-spike stages in the peak plane (y = Ly/2) with superimposed contours of negative λ2. Left: temporal evolution of data obtained from the spatial LES at a fixed location x. Right: temporal simulation. 

根据能量梯度理论的预测，奇点位于中心面两侧正负z/h=0.6的位置，实验结果已经证明最初产生的奇点位于正负z/h=0.6的位置（Nishioka et al. 1975）。关于奇点的分布，图5的LES (文中也包括DNS结果）数值模拟预测结果与能量梯度理论结果基本一致。从上图可以看出流道纵剖面上空间和时间分布的奇点的演化过程，随后，这些奇点就会演化成为转捩流动和湍流中的大尺度流动结构。实验得到的湍流的大尺度拟序结构，正是这些奇点演化的外在表现。窦华书根据能量梯度理论和NS方程，提出了湍流是NS方程的奇点产生、演化和发展的结果，完全发展的湍流实际上由流场中大量的奇点所构成，这些理论分析结果与所得到的大量实验数据一致。