注:此篇文章是在知乎网上回答网友的问题“流体一定是从压强大的地方流向压强小的地方,对吗?依据是什么?”而写的,以促进流体力学和湍流问题研究在国内外华人学术界的研究进展。 ! d8 f4 l8 `: {
* l J+ f1 q6 ~$ f" _' Q
图1 巴西伊泰普水电站--世界七大奇迹之一 0 h# v3 {$ a/ _ P2 u" ^
N5 B8 p3 R) j6 `% Q
(一)流体一定是从压强大的地方流向压强小的地方,对吗?依据是什么?4 |* t0 |- B1 e) s. a
流体不一定是从压强大的地方流向压强小的地方,一般与压强的大小没有必然关系,因为压强大小不是决定流体流动的准则。流动要根据能量守恒定律,要看总的机械能大小和分布来决定,还有外部做功以及能量损失来决定。8 X8 }. R/ ^4 ~/ \8 Y. o3 x7 n: Y
学过流体力学的人都知道,一般来说,在截面积逐渐缩小的管道里面是压强逐渐减小,速度增大,这时是流体从压强大的地方流向压强小的地方。而在截面积逐渐增大的管道里面是压强逐渐增大,速度减小,这时是流体从压强小的地方流向压强大的地方。但是,大家要明白,在这两种情况下,入口的机械能肯定是一定大于出口的机械能的。否则,流体一定是向回倒流的。) p( ]9 t6 Y+ f1 I. r
(1)对不可压缩粘性流体,如果外部没有对流体做功,流体一定是从总的机械能大的地方流向总的机械能小的地方。即:3 c, d+ M, P* H4 ^0 y5 c- b0 F
如果E1>E2, 则流体从1流向2,其中 E为总的机械能。E1-E2=H,就是沿程能量损失了。
: H; |/ R6 c, I* r2 s$ Z; B7 l例如水管子里面的水,进水口处E1, 大于出水口处E2, E1>E2, 这样水才能从入口流向出口。如果发生入口的E1小于出口处E2的情况, E1<E2,这时水就会从出口向入口处倒流了。这里我没有说水管子是等直径的管道,还是逐渐缩小,或者逐渐扩大的管道,都没有关系。
( |: B1 W; l( l) m" a(2)对不可压缩粘性流体,如外部对流体做功了,上述是不一定成立的。要根据热力学第一定律(能量守恒定律)来定。
3 y z" w. V4 p7 p% q如果E1+W>E2, 则流体从1流向2,其中 E为总的机械能, W为外部对流体做的功。E1+W-E2=H, 就是沿程能量损失了。
8 v- ?- G% u) W0 w' O8 _例如水泵入口E1小于水泵出口E2, E1<E2,但水泵叶轮对水做功了(W), 这样 E1+W>E2, 水才能从水泵入口1流向水泵出口2。
: S. s/ `' I" c如果水泵的扬程不够(W小了), E1+W<E2, 水泵是不能把水从水泵入口1送到水泵出口2的。相反,你用水桶提的一桶水倒入水泵出口处,都通过水泵进行倒流而流向水泵入口处了,做过水泵实验的同学都见过吧?在农村用水泵浇过地的同学也可能见过的。
! Y7 y/ X/ }+ |再例如水力发电站用水轮机(图1),水轮机入口E1大于水轮机出口E2, E1>E2,水对水轮机转轮做了功W,即水轮机转轮对水做了负功(-W), 这样 E1-W>E2。E1-W-E2=H,H就是由于水的粘性引起的水力损失(图1)。
. {# k2 A: R) `7 X( m6 O/ f注: 总的机械能=压力能+动能+位能(势能)。! \! R6 z! C: R3 Z- O# k
(3)流动速度与能量损失的相关性
; T4 y3 T4 I5 Z; I# o6 z通过上面两种情况的理论分析,对于流体为什么会流动,我们可以总结出一个推论(Corollary): 流体流动存在的必要及充分条件是沿程能量损失(H)不为零。而且发现,沿程流动损失越大,速度越大。如果沿程流动损失为零,则速度为零。这个推论,恰恰是湍流生成的物理机理的关键所在,湍流产生的秘密就在这里 [1]。# ^8 a; r. ?/ t7 ^. V: h+ o
没有能量损失(能量差),就不能产生速度。有了能量损失(能量差),就必然产生速度。有了速度,就必然产生能量损失。没有速度,就不能产生能量损失。速度产生的必要及充分条件是存在能量梯度或有外部做功。速度和能量损失(能量梯度及做功)互为其单值函数。
. j0 N# }* B8 ~" M1 B; H+ Z8 k说到这里,这里再次顺便通俗易懂地简要讲一下湍流是怎么产生的。当大量流体质点在流动时,每一个质点都要消耗一定能量(有一个能量梯度),才产生一个速度(相互平衡)。在扰动影响下,如果流体内部突然有一个流体质点上的能量损失(能量梯度)变为零,则这个质点的速度就突变为零。这时周围的质点还是继续走的,这样此处立刻引起了混乱,引起了湍流旋涡,湍流就产生了,这就是湍流的起源,而且这种说法与实验结果是一致的。
: B# F% I# L$ M6 x作者得出的这个湍流产生的机理对wall bounded flow and free shear flow 都是一样的,是通用的,湍流由一个唯一的机理支配,因为控制方程是同一个NS方程。而层流到湍流转捩的机理由NS方程的奇点出现的临界条件唯一地确定。这一点与国际上的最近10年的许多文章的看法是不一样的,他们认为在这两种类型的流动中,支配湍流产生和结构维持的机理是不同的。本文作者认为这种看法是错误的,在剪切湍流中,无论是wall bounded flow还是 free shear flow,湍流产生的机理应该是唯一的。举一个例子,我们在用水壶烧水,水沸腾开了的准则是水的温度达到了100度。你用一个大铁锅去烧水,水开了的准则还是水的温度达到了100度,临界点的准则是不会变的。不管你是改变了几何形状,还是改变了边界条件,由同一个数学物理方程所支配的产生转捩的临界条件的准则是不会改变的。
( k, N; w) g, G3 @许多湍流专家认为湍流是速度剖面的剪切通过壁面边界层里内外层的各种不稳定性作用导致了速度的脉动,而且还认为湍流转捩和湍流维持的机理是不同的。还有专家认为wall bounded flow里湍流的产生是由壁面作用生成的,在早些年代,还有的专家认为边界层湍流是由于边界层外部流动产生的。这些研究结果虽然发表在了国际流体力学顶级期刊上,但是没有坚实的理论支撑,也没有实验数据的支持。本文作者提出了能量梯度理论,发现湍流产生是速度剖面上同一条流线上的来流的动能通过奇点转换(非线性失稳)成了速度的脉动,这个奇点的产生正是湍流产生的临界条件,这是根据非定常的三维NS方程精确推导出来的,而且理论与实验是一致的。
6 c+ J$ W' I$ A作者上面讲的奇点,一定不在壁面上,是在流体内部,是在流体内部速度理论上瞬时为零的点,是速度间断点。反之不一定成立,流体内部速度为零的点不一定是奇点,要看有没有速度间断。
' F; E1 F: J0 u7 e5 H上述这是作者给研究生和本科生上课讲的部分内容。如果想进一步理解,请参考下面文献[1]的第二章,第五章及第九章273-275页内容。6 O7 k( w5 Z' j" Z
(二)流体流动及失稳的基本问题:, E0 B6 x5 D! p5 L& M. c/ a
从回答上述问题,我们可以进一步讨论两个重要的基础问题,就是(a)流体为什么会流动?(b)流动为什么会失稳?现代流体力学并没有回答这2个最最基础的问题。这只能说牛顿力学或者说至少现代流体力学还不完备。如果这2个问题早一些时间回答了,可能湍流的问题早就解决了。
. E/ P- f4 j. M牛顿力学的三大定律:(a)惯性定律;(b)F=ma;(c)作用力与反作用力。牛顿力学定律奠定了经典力学的基础,成为现代科学的基石。飞机飞行,火箭发射,航天飞船,水泵设计,桥梁建筑,水利工程,等等都是依靠牛顿力学的基本原理。
4 u, M: k9 Y4 c% [) V- A可能有的同学会问,牛顿力学诞生于17世纪(牛顿,公元1643-1727),中国的都江堰工程由李冰父子在战国时期(公元前256-251年)修建,赵州桥由李春于隋朝时期修建(公元605-618),那时还没有牛顿,这些工程怎么会利用牛顿力学。那时这几位古代贤人虽然没有从实践中提炼出牛顿力学,但是在他们设计的工程中实际上已经应用上了牛顿力学的原理,例如作用力与反作用力的原理。中国古代很早(隋朝)就发明了火药,春节元宵期间,民间放的“钻天炮”就是今天火箭的雏形,是典型的牛顿第三定律的应用范例。到了21世纪的今天,可能人们仍然需要从实验中提炼自然科学的规律,例如湍流的产生原理。" a- e- O. u% o' R" H, a3 D
需要指出的是,对于(a)流体为什么会流动?(b)流动为什么会失稳?这2个基础科学问题,至今没有人给出什么定律或者公理。可能,我们需要牛顿第四定律,甚至牛顿第五定律。
a, G/ k! N; j% B. a为了顺利解决湍流的问题,作者苦苦追寻30多年,针对不可压缩流体,作者在文献[1]中的第五章给出了下面4个公理(Axiom)。2 }7 \8 G& h. t
(1)公理 5.1: 在实验室系统中,如果总机械能沿流动方向的梯度小于零,则流动是可能的。
& i c* \: C3 f' F! g: A(2)公理 5.2: 在实验室系统中,如果沿着垂直于流线的方向的总机械能梯度为零,则流动是稳定的,否则可能不稳定。3 m4 S) C' |6 ]; k
(3)公理5.3: 在实验室系统中,如果沿着流线的总机械能梯度为零,则流速为零。* L* `5 u( J2 O$ ~4 n! d
(4)公理 5.4: 在实验室系统中,在剪切驱动的平行流中,如果对流体单元上的剪切应力做功为零,则流动将立刻停止,速度为零。
- U# V) C8 [' E/ H7 O- S7 g上面的公理 5.1,5.2,5.3适用于压力驱动流动(pipe flow, channel flow, Dean flow, etc);公理5.4适用于剪切驱动流动 (plane Couette flow, etc)。这两种情况的区别就是看看有没有外部对流体做功。有外部做功和没有外部做功是完全不同的。因为外部对流体做功可以改变流体运动的方向,可以加速或延迟流动的失稳,可以控制湍流的发展。
. {, a5 [: Y" t: @# h为什么要强调“在实验室系统中”,因为能量和做功都不是伽利略不变量,它们随着坐标系转换就不成立了。什么是公理,就是自然界被普遍认可的科学规律。公理与定理的区别是,公理是不需要证明的,也不可能证明。公理可以用实例、用计算和实验结果进行验证,但是公理不允许存在反例。定理是需要根据已有的科学规律和定理去进行逻辑证明的,定理也不允许存在反例,如果存在一个反例,那就是错的。% ?$ n2 C9 U% S" R8 O
我们有了上面这4个公理,流体为什么会流动及流动方向如何,流动为什么会失稳,湍流为什么能够产生,根据速度场的理论分析结果,立刻就能解决了。而且,根据这4个公理得到的定理是普遍适用的。正是基于上述4个公理,“湍流是怎么产生的”这个百年难题终于被成功地解决了[1-6]。$ ?* Q. F% A: K$ g9 F9 j; y/ t
首先解决了湍流专家们最最急切关心的几个问题:(a)湍流怎么产生的;(b)湍流大涡怎么出现的;(c)湍流怎么从平均流动获得能量的。其次理论解释了湍流的其他若干问题,如间歇性、拟序结构、湍流自持过程、重复性、NS正则性、波涡干扰、尾迹稳定性、转捩的扰动标度律,等等 [1-6]。
) Y, ~! `* c/ S5 m' y6 E作者利用所建立的能量梯度理论得到的2个最突出的重大贡献是[1,6]:. H( f) \7 S8 J0 V/ s
(1)定理1:湍流转捩或湍流产生的必要既充分条件是流场中出现Navier-Stokes方程的奇点(速度发生间断)(此定理在物理学上解决了百年湍流难题)。
, r( V( E! t6 A& y(2)定理2:对转捩流动和湍流流动,Navier-Stokes方程不存在全局定义域上的光滑解(此定理在数学上精确地回答了千禧年大奖难题)。
% ` l5 @, q$ a作者提出的能量梯度理论是唯一的一个根据NS方程预测湍流产生和解释湍流维持的理论,也是唯一的一个与实验数据符合一致的理论。在国际上,首次解决了湍流世纪难题; 首次对千禧年大奖难题--Navier-Stokes方程的存在性及光滑性,给出了答案[1-6]。
$ i5 o+ ^" l8 w+ G# u在基础研究中和工程设计研究中,所有的湍流数值计算,所有的湍流modelling的可靠性及精确性,都应该基于对上面两个定理的理解,基于根据湍流生成原理的数值处理。
: a1 ^4 o7 f f9 m* y3 D4 O5 q科学真理终究是科学真理,要经得起实验验证和时间的考验。
- U, j# F1 j+ F2 T x( Q上述这是作者给博士研究生上课讲的部分内容,也是作者创立的能量梯度理论的部分内容,在其他教科书上是没有的。! v" r$ t: x- M. B8 A9 B |4 C
参考文献
! t' E6 Y r* v) A6 g1. Dou, H.-S., Origin of Turbulence-Energy Gradient Theory, 2022, Springer. www.52ocean.cn (全书下载地址).* I- e6 X7 O5 X( j
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% |5 {2 V: g9 g1 t
3. 科学网-海森堡的第二个问题终于有了答案 - 窦华书的博文。
- c* ]2 E( K% Q' o: awww.52ocean.cn7 u* B5 w5 w' B {: }/ c
4. 新书访谈,专访《湍流的起源—能量梯度理论》作者窦华书教授。 www.52ocean.cn* g N" C2 P' S
5. 窦华书教授成功破解了百年湍流难题,中国教育日报网。http://chinaedutech.com/dfjy/2022/1117/1327.html
# m1 j* k; U, d6.窦华书教授在纳维-斯托克斯方程问题上取得新进展,浙江理工大学官网新闻。 www.52ocean.cn. d" D" Z9 P0 F1 h+ |
3 T2 `" e7 Q% X9 r; j) c1 b* {9 t& G7 b- Q* A5 e! `
2 Y2 [; d( ]/ n7 N7 D) Q3 f, ]
4 e/ Y/ }$ J3 ^2 k 转载本文请联系原作者获取授权,同时请注明本文来自窦华书科学网博客。 |