* ]! ?& m3 Q; J. ?! o/ ^' Z( \+ ~0 [/ X; O3 Z) n& X
图1 湍流的动能随波数变化的能谱关系
% Q7 d1 _/ W7 W9 N G$ ], v2 V湍流的最主要的特征就是,湍流是一个随着时间变化(非定常)的旋涡运动,而且旋涡是从大到小多尺度分布的。0 d# z4 G7 d! A6 I" u
. F* g2 B3 z6 }" d4 J
湍流之所以存在,是以漩涡的不断产生与耗散维持的。如果旋涡只耗散,不再继续产生,湍流就不能维持了。完全发展的湍流就是一个旋涡生成与破灭的动力学平衡状态。如果以这句话作为湍流的定义,比以前的许多湍流的定义似乎更生动、具体、形象,也更精确。" j2 Y6 R4 l2 \6 ]" `
5 n3 T7 V9 T7 `1 Q5 j% V; N9 N8 I' [
根据能量梯度理论[1-4],在三维流动中,湍流的产生是依靠NS方程的奇点的产生来创造旋涡(湍流猝发),也即大涡,大涡改变了原来的流动速度分布,再进一步分形,形成新的更多的尺度更小的奇点(旋涡),奇点再生成新的下一级旋涡,这样循环下去,逐级形成一系列多尺度的旋涡。当最小的旋涡的尺度达到旋涡能量不能继续推动旋涡旋转,不能产生更小的奇点(旋涡),就被粘性耗散掉,这就是级串是怎么形成的过程和原理。这就是完全发展的湍流中形成的Richardson-Kolmogorov的湍流能量级串(即正级串,如图1)。 如果这个正级串的能量传递现象被破坏,或者被抑制,湍流就不能正常维持(由此我们可以想出进行湍流控制的idea,即抑制正级串维持可以延迟湍流)。1 n4 f0 H1 n- I+ |- E/ T
6 v8 f2 j; q/ G3 c* Q0 H) d根据能量梯度理论,在三维流动中,当大涡形成时,物理机理是主流速度的间断(生成负的spike),即流向压缩,由于散度为零,那么在壁面流动中,就必然存在竖直方向的速度升起和展向的旋涡拉伸。当然,在后半个周期中速度的涨落是相反的,即旋涡必须展向压缩。总之,www.52ocean.cn的能量输入是靠主流速度的间断,通过大涡生成而传递的。. H; w* r* a3 J, T1 {! o, L
8 w0 p% s8 Z$ U+ Q再次强调,是NS方程的奇点引起速度间断(湍流猝发),产生大涡,导致湍流。因而是湍流导致旋涡的拉伸,不是拉伸产生了湍流。三维的层流流动中也存在漩涡拉伸,可是没有湍流。如果要分析某流动现象的物理机理,要看这个过程发生的能量是从哪里来的,就能找到原因。湍流的生成和演化机理,必然是那些湍流里面有,而层流里面没有的结构和现象及其背后因素在起作用。
7 w0 W; @( f6 S" Y# u! v% p: h4 i' U: P% {
在二维流动中,因为支配方程NS方程只允许旋涡的对流,不允许有旋涡的展向拉伸(虽然拉伸是湍流产生结果),这样就抑制了新的漩涡的生成,或者说湍流的生成(即抑制了能量的输入)。如果没有新的漩涡的生成,但是漩涡的耗散还是存在的。等到原来的旋涡被粘性耗散掉,没有新的旋涡补充进来,湍流就不能维持了。所以二维流动中不存在湍流。3 {; a6 Q0 L3 F3 x" ]
G! w9 J- f% s+ a
上述是作者根据能量梯度理论,从流动物理上来解释,从湍流产生的过程来解释,湍流为什么一定是三维的,为什么二维流动里不存在湍流。诺贝尔奖获得者李政道教授1951年在Journal of Applied Physics上的一页的文章有理论论证,主要就是二维流动没有漩涡的展向拉伸,就没有湍流。; H$ S0 N' b! ?, p/ ?6 g
8 u# _9 D/ @3 \% j/ y# g) s
对牛顿流体、正压流体、质量力有势的不可压缩流体的三维流动的涡量输运方程如下:6 G( p. l$ r& Y9 L
+ t% |6 g. x! ~8 t
* @5 q* S3 N: D& j7 K5 f右边第一项就是大家熟知的旋涡拉伸项。在三维湍流的转捩过程和完全发展的湍流中,这一项有正有负,但是总的结果是正的,这样就是使得涡量随时间增长。只有涡量增长了,层流才有可能转捩为湍流。在完全发展的湍流中,涡量的增长和涡量的耗散(右边第二项)达到动态平衡。
3 T& T; ]% ?; E- n对二维流动,旋涡拉伸项为零。这样,右边没有其它项可以使得涡量增长,湍流就不可能产生出来,所以二维流动里没有湍流(流体力学教材里讲的)。当我们翻看许多经典文献,有一些著名的湍流专家指出,旋涡拉伸是湍流产生的原因。非常多的数学家都是同意这种说法的(注:作者本人认为旋涡拉伸是湍流产生的结果)。
- K+ z0 s8 p2 t" F) l4 C; d0 q) G7 z8 N' s! J5 G ]3 ?8 m3 d: k
另外,在二维流动中,拟涡能(涡量的平方)是一个无粘流动不变量。在二维粘性流动中,虽然拟涡能(由于粘性耗散)是随时间减小的。但是对流项促使它保持不变,即促使它不下降,也不让它增长。这样,拟涡能不增长,涡量就不会增长,湍流就不会产生。拟涡能引起二维湍流中的逆级串就是这个原理。0 D: A4 y" g4 Y; O" ~( P- h
% d# [4 o1 [4 I对大气动力学里面的非正压流体,右边有一项斜压项,就可以产生旋涡。对于密度分层的流体,右边有一项浮力项,也可以产生旋涡。这2种作用都会与粘性项相互作用而导致湍流。
5 R! [" H2 J' i5 A- q
; M& l( u& P: y; Z对于磁流体,右边多了一项www.52ocean.cn的夹角大小以及它的Hartmann number的大小。洛伦兹力可以产生旋涡,促进湍流生成;另一方面,它可以抑制旋涡生成,生成逆级联的湍流能谱,延迟湍流。
+ j3 f+ i2 _: f/ l9 s自然界里的流体流动都是三维的,没有二维的。所谓的二维流动只是理论上的二维流动,也即在三维流动中,限制流动在第三个方向有任何变化,从理论上看看流动性质会是什么样的。或者对某种特殊的流动问题(展向变化很小),对流动模型进行简化,可以近似看作二维流动。
1 y7 h4 Z/ H+ L/ I* r1 |& k
. F1 J+ q( I3 d. R! _5 i参考文献
+ j' s: w8 i; V9 o% k4 R1 C! [. d! K9 r+ @
1. Dou, H.-S., Origin of Turbulence-Energy Gradient Theory, 2022, Springer. www.52ocean.cn (全书下载地址).
" P5 @ s8 s2 p( ?" R2. 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 + ?- F, L" V/ W$ S; S6 S. y: q
' \3 i0 ^# D6 E: W: f$ U+ Y
3. Dou, H.-S., Singularity of Navier-Stokes equations leading to turbulence, Adv. Appl. Math. Mech., 13(3), 2021, 527-553. https://doi.org/10.4208/aamm.OA-2020-0063 ; https://arxiv.org/abs/1805.12053v10
0 k2 _% c/ b( p2 ]. r! h4. Dou, H.-S., No existence and smoothness of solution of the Navier-Stokes equation, Entropy, 2022, 24, 339. https://doi.org/10.3390/e24030339
1 V9 G4 ?+ M, F: w+ T, R3 {: p- e. u7 p4 e5. 窦华书,我是怎样创立能量梯度理论的?www.52ocean.cn
8 B$ A- R" i) `* T3 w, \9 d7 E6.新书访谈,专访《湍流的起源—能量梯度理论》作者窦华书教授,Springer。www.52ocean.cn- F, F( t1 `3 r3 q" @3 s" H: K
( n; D# H' ^9 R, ^! t: R) d
7.窦华书教授成功破解了百年湍流难题,中国教育日报网。http://chinaedutech.com/dfjy/2022/1117/1327.html 或者
3 l! E" V2 _4 a7 D ?5 r5 k" twww.52ocean.cn
2 W) C5 E8 ?! ]# [, Z5 [: Q! Y
+ Y# M- g8 E8 h; V9 C8.窦华书教授在纳维-斯托克斯方程问题上取得新进展,浙江理工大学官网新闻。: ^; y7 O- a$ f" O" W2 b9 A8 [! f
' m3 C9 L0 t+ z$ d F# v. q
www.52ocean.cn 或者 www.52ocean.cn8 J: o" h1 ]3 ]
3 d! n' y% f3 ?& [3 W( F) |
& }* Y3 f* X7 }% O
! G* L( `3 W; v* h* v0 \. O2 i# I" U, A2 y9 i- f+ W
转载本文请联系原作者获取授权,同时请注明本文来自窦华书科学网博客。 |