首页/文章/ 详情

雷诺数的前世今生

1月前浏览526
       

自从人们发现湍流无规则的现象以后,便有无数的流体力学家们为了降服它而前赴后继。奈何湍流如同会七十二变的孙猴子一样闹腾,时至今日也无法完全收服。《西游记》的总策划如来佛祖为了让孙大圣乖乖的护送唐僧去西天取经,就安排观音菩萨下界给他脑袋上套了一个“紧箍”。有了紧箍的束缚,孙猴子似乎老实了许多,总算给取经路上添了个好帮手。


 


而面对更加上蹿下跳的湍流,流体力学家们便决定给它也套上一个紧箍——雷诺数,期望让湍流变得有迹可循。那么作为湍流的“紧箍”,雷诺数究竟是怎么来的呢?它和湍流之间又有哪些秘密?《西游记》开篇有言:“混沌未分天地乱,茫茫渺渺无人见。自从盘 古破鸿蒙,开辟从兹清浊辨。”本着公众 号一直坚持的钢铁直男的探索精神,我们也从雷诺数的起源说起......


01

雷诺数的起源


最早提出以无量纲数理论来预测流体流动的是N-S组合中的S——乔治·斯托克斯。1851年,他曾试图计算忽略惯性项的球体阻力,并导出了作用在球形物体上的摩擦力(雷诺数很小时也即阻力)的表达式,即斯托克斯定律。


 


人们基于斯托克斯的工作发明了落球粘度计——让一个已知大小和密度的球体通过液体下降,通过测量小球的终端速度,便可利用斯托克斯定律计算流体的粘度。


“形而上者谓之道,形而下者谓之器”,比起落球粘度计这个具体的器物,斯托克斯的思想显然更为重要——斯托克斯定律很好的启发了后来的流体工作者使用无量纲数来预测流动现象。


02

经典的雷诺实验


用无量纲数来预测流态是个很好的想法,但还未来得及进行数学实施。另一边的雷诺则迫不及待的做起了实验,用染色的墨水向人们展示流态当圆管中的流速较低时,流体平稳地沿轴向运动,而无其它方向的流动,流体的层次分明,称为层流。随着流速的增大,原本平稳的流动便逐渐演化为杂乱无章的流动,即为湍流。


 


和蔼可亲的雷诺堪称流体力学界的太白金星,他通过实验向人们展示了湍流这只“泼猴”的百般变化并验证了流体的状态与管道直径、流速、流体的密度和粘度等参数相关。可是,湍流触发的机理究竟是什么?何时会产生?湍流的规律又是什么?如何驯服“一天闹腾24小时”的湍流,不仅困扰着整个天界,更是让雷诺夜不能寐。


既然不易驯服,那就招安吧。雷诺实验5年后,大彻大悟的雷诺不再纠结于湍流的复杂运动,而是将湍流分解成平均速度和脉动速度,并将流体力学的葵花宝典——N-S方程推导成雷诺方程,开启了一段快乐时光。此时的湍流就像被封了“弼马温”的猴子一般,着实老实了一段时间。


03

克莱因的天才直觉


流体力学家们可以用这种化繁为简的方式演绎湍流,可是内功深厚的数学家却觉得它总有一天会大闹天宫,必须斩草除根方得罢休。众多数学家们都想从数学层面上解释湍流触发的机理,其中尤以星光璀璨的哥廷根学派和慕尼黑学派最甚。


作为哥廷根学派的太上老君,克莱因手握着可以容纳整个宇宙的紫金葫芦——克莱因瓶,在数学界和物理学界挥斥方遒。不过众人皆知,紫金葫芦虽法力无穷,但在足智多谋的孙大圣面前一样也哑了火。


 


虽然紫金葫芦无法收服湍流,不过太上老君还是够老辣。关于雷诺实验,克莱因仍然有自己独特的见解。一次课后和他的助手索末菲闲谈时,克莱因便提到,为了解释管道中湍流运动的出现,我们可以认为,当管道中的流速超过临界速度时,层流变成了一种“不稳定”的运动形式,但是,不稳定发生的原因尚不清楚。因此,湍流的触发机理或许可以转化为一个稳定性问题。


言者无心,听者有意。作为克莱因助手的索末菲当即认为,或许可以以此为出发点,从理论上求解层流到湍流的过渡。然而十多年的时光飞逝,索末菲娶妻生子,并从哥廷根大学去到了亚琛工大,又来到他大放异彩的慕尼黑大学,并开启了他统领理论物理的时代。可是关于湍流触发的问题,依然没有答案。


04

雷诺数的命名


无论索末菲在量子力学的领域有多么大的成就,湍流触发机理始终像他无法得到的“床前明月光”一样,时刻萦绕在他的心头。到了1908年,恰逢第四届国际数学大会的邀请,索末菲将计就计,将他对流动稳定性和湍流触发机理的思考写成文章,邀请整个学界的同仁陪他一起纠缠。


这篇文章在历史上第一次明确以“雷诺数”命名了前人提到的某个神奇的无量纲数。于是流体力学中最重要的无量纲数终于正式登场:“如果ρ是密度,μ是流体的粘度,U是横截面上的平均速度,b是管道的直径,那么R=ρUb/μ是一个纯数,我们称之为雷诺数。


另外,在文章中,索末菲还写下了著名的奥尔-索末菲方程,并明确应用经典的小扰动方法来解决湍流触发的问题。


 


所有看到奥尔-索末菲方程的人都拍案叫绝,可是却没有人能解出来,这就像有人造出了紧箍,但没有人会念咒语果然,奥尔-索末菲方程在接下来十几年的时间困扰着整个数学和物理界。


05

特征雷诺数的纠缠之路


本着“我们解决不了的问题,就留给后面的人来做”的科学思想,索末菲把方程的求解丢给了自己的博士生,另一个量子力学的未来之星海森堡——就是那个曾经感慨,见到上帝之后一定要问两个问题,一个是“相对论”,另一个就是“湍流”的物理学家。可见海森堡被湍流虐的有多惨。


回到海森堡的博士课题上来,奥尔-索末菲方程岂是那么容易求解的。海森堡使用的小参数摄动法几乎在一开始就碰到了无法逾越的鸿沟,即解的收敛性问题和奇点问题。博士的生涯总是过得很快,眼看博士答辩的期限将至,海森堡还是没能找到答案。不过凭借天才的直觉,海森堡猜到了解的大概样子,并拼凑出了一个临界稳定雷诺数的解答最后在导师索末菲的保护下顺(蒙)利(混)毕(过)业(关)了。


 


尽管海森堡的工作在这一领域具有很重要的建设意义。但其在数学分析上使用了一定的限制和近似,使得海森堡的博士论文并未得到整个学界的完全认可。其后二十年,人们都没有找到驾驭特征雷诺数的“紧箍咒”。


06

来自东方的神秘力量


直到1941年,来自东土大唐的林家翘携带着东方的神秘力量来到加州理工大学攻读博士学位。在导师冯卡门的建议下,林家翘也选择了奥尔-索末菲方程作为自己的博士课题。在博士论文中,林家翘使用精妙的渐进方法解出了特征方程,并同时给出了对称型和边界层型两类速度剖面下的特征雷诺数,并基于稳定性理论给出了流体稳定的“拇指图”曲线。整个过程逻辑严谨,堪称此类问题求解的典范。


 


1944年,林家翘的博士论文一经问世,便引发了整个学界的关注,也终于结束了特征雷诺数在数学求解层面的纠缠。更令人难以置信的是,林家翘的结果居然和海森堡当年凭借直觉推测出来的结果基本一致。天才果然就是这么豪横!


07

高雷诺数就是湍流吗?


按照林家翘的理论推导,圆管流动的临界雷诺数为5772,这个值和小伙伴们在流体力学课本上学到的2300偏差不小,这也说明了湍流本身的复杂性使得理论很难完全覆盖实际的流动情况。另外,随着人们对于湍流的理解不断加深,大家明白影响转捩的不只是雷诺数,还有扰动这个重要的因素。


于是,很多人都试图通过消除扰动来获取临界雷诺数,这个值在1980年甚至达到了107这个量级。因此,我们有理由相信,临界雷诺数可以无穷大,对应零扰动的流动。


 


如果把稳定性理论应用到泊肃叶流动,我们可以得到类似上图所示的拇指图。纵坐标是扰动,横坐标是雷诺数,可以看到当扰动足够小的时候,需要非常高的雷诺数才能触发流态的变化。


08

湍流和雷诺数之间的密语


即便有了雷诺数这个“紧箍”,加上众多大神们不知如何推导出来的“紧箍咒”,湍流世界似乎依然难以捉摸。而关于高雷诺数的湍流运动,其复杂程度大概不亚于齐天大圣大闹天宫般的混乱吧。幸得仙界还有另外一群老神仙,仍凭你猴头上蹿下跳,我自岿然不动,恪守着“万变不离统计”的思想,将复杂的湍流和雷诺数联系在了一起。


 


还记得我们在前面的文章中提到过的“流场中的大涡尺寸L0是由边界条件决定的,而小涡的尺寸η是由雷诺数决定的”比如上图所示的方块绕流,决定涡旋大尺度的一定是物体的尺寸,而非雷诺数,即无论来流的速度大小,尾流中最大涡旋的尺度基本上都与物体尺寸保持相同的量级。关于雷诺数和小涡尺度之间的关系,则要重提我们之前一直给大家念叨的湍流的能级串理论:湍流中遍布着大大小小的涡系结构,能量就是在这些涡系中相互传递;一般而言,大涡生成能量,依次传递给小涡,并最终耗散。而对于一个平衡的系统,比如完全发展的湍流,大涡生成的能量不多不少的被小涡全部耗散掉。


09

雷诺数密语的数学表达


现象和理论有了,老神仙们追寻的则是更加极致的数学表达:对于充分发展的湍流,可以认为决定小尺度湍流现象(耗散尺度η)的参数只有流体的运动粘性ν(决定传递)和能量耗散率ε(决定耗散)。而此时,耗散率ε又等于大涡传递过来的能量。于是,雷诺数与大涡尺度L0和小涡尺度η之间的关系便可以推导出来:


 


可以看到,对于给定边界的流动,即大涡尺度L0一定。雷诺数越高,则耗散尺度η越小,即流场中的小涡尺度越小。因此,也可以认为,雷诺数越高,流场中从大涡到小涡之间的过渡越复杂,涡系结构越丰富。


1974年刊登于《流体力学》期刊中的一篇经典文章,便用流动显示技术把这一机理进行了生动形象的描绘。如下图所示,随着雷诺数的增大,涡旋的整体形状和尺度基本一致,然而射流掺混的涡系结构却越来越丰富。


 


带着紧箍的孙大圣陪着唐僧一路西行并取得真经,佛祖也顺势取下了大圣头上的紧箍,而获得佛法感悟的孙猴子,即使没有了紧箍,也能够造福人间。而对于更加闹腾的湍流,雷诺数这个“紧箍”在相对低速的工况时或许还能够勉强应对,对于高速流动,甚至是超音速流动,湍流则变得更加的“嚣张”。幸好,流体力学的老神仙们又想出了马赫数这个“佛法”来感化高速下的流动。敬请期待下一篇“突破音障的LBM”,带给你不一样的马赫数体验。


-END-


来源:CFD仿真区
湍流理论管道
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2024-10-13
最近编辑:1月前
濮小川CFD
硕士 心不唤物,物不至!
获赞 16粉丝 42文章 103课程 0
点赞
收藏
作者推荐

流体江湖最耀眼的“过客”

前段时间在网上看到一个有趣的问题:“为什么前苏联以及俄罗斯的数学那么强,但许多数学家却不为人所熟知?”有一个抖机灵的答案挺有意思:数学强是因为研究数学在前苏联不会犯政治错误,不为人知是因为诺贝尔奖中没有数学。虽然这个答案道出了一些要素,不过自从伯努利的两个儿子以及大数学家欧拉入驻俄国皇家科学院的那刻起,数学的基因便深植于俄国的血脉之中。 如果说欧美的科学家在学术界的地位是武林的泰山北斗—少林和武当,那么苏俄的科学家便是神秘而强大的“逍遥派”。前苏联鼎盛时期,号称有一万名数学家,甚至可以跟整个西方数学界抗衡。而莫斯科大学汇聚了众多有名的苏联科学家,组成了赫赫有名的莫斯科学派。而我们提到过无数次的老神仙柯尔莫果洛夫,无疑是莫斯科学派中的翘楚。01逍遥派的集大成者作为逍遥派的集大成者,柯尔莫果洛夫的数学研究范围非常广泛,他在不同时期对基础数学、数理逻辑、实变函数论、微分方程、概率论、数理统计、信息论、泛函分析、拓朴学等众多数学分支均有所涉猎,而且在他耕耘过的许多领域都取得了基础而重要的成果。 柯尔莫果洛夫的重要成就还包括他将数学广泛的应用到了其它专业领域,比如物理、化学、生物、地质、冶金、信息技术等学科。在流体力学领域,柯尔莫果洛夫也跨界双修,并成为一代宗师,而其最重要的成就莫过于他创立的K41理论,首次揭开了湍流空间分布特性的面纱,成为湍流星空中的一轮皓月。02平凡的出身不平凡研究数学在欧洲曾经是一项贵族运动,大家熟知的很多数学家都是含着金汤匙出生,从小就高唱着“小小少年,没有烦恼,万事都有老爹罩”,一边在学术界挥斥方遒,另一边也在各种科研机构觥筹交错,形成了各自的圈子。相比之下,柯尔莫果洛夫的出身则显得平平无奇。1903年4月25日,他出生于莫斯科东南部的小城坦波夫(Tambov)。当年的俄国还发生了一件大事,即布尔什维克党成立。所以柯尔莫果洛夫成长于沙皇俄国的末期,那是一个政治动荡而又思想迸发的年代,而这种艰苦又不安定的环境也似乎更加能够激发人们的意志和创造力。回顾历史,我们也会发现,许多大家都是这种环境下成长为让世人仰望的参天大树。 柯尔莫果洛夫的父亲据称是一名博学的园艺师,但因参加反对沙皇的革命而被流放,后在内战中失踪。而其母亲在当时虽说家庭条件不差,但却在生下他不久后不幸离世。柯尔莫果洛夫的名字来源于他的外祖父,而他从小由姨妈养大。虽然无双亲照料,但姨妈如母亲般的悉心照料也抚慰了他幼小的心灵,并开始了其不平凡的求学经历。和所有豪横的天才一样,柯尔莫果洛夫在幼年时期便展现出对于数字超乎常人的理解,他在五、六岁时便发现奇数之和的规律:1=12,1+3=22,1+3+5=32…,将其发表在家庭课堂的期刊《春燕》上,同时他还负责编写这份期刊的趣味数学部分。1910年,柯尔莫果洛夫随姨妈移居莫斯科,进入前列普曼文法学校学习,开启了更广泛的数学启蒙。03初入逍遥派1920年,柯尔莫果洛夫经过短暂的工作之后,进入莫斯科大学主修物理数学。或许是天才少年的青春实在难以安放,除了主修的物理数学之外,柯尔莫果洛夫又去了门捷列夫化工学院学习冶金,大学期间他还对历史感兴趣并将统计方法引入历史研究中。几年的大学生涯让天才少年渐渐的看清了自己的志向,并决心投身数学的研究。1925年,柯尔莫果洛夫毕业后被著名的数学分析大师鲁金看中,并开始与师兄辛钦一起从事概率论的研究。 此时正值鲁金学派的鼎盛时期,柯尔莫果洛夫就如同进了逍遥派的琅嬛玉洞一样,感受着武学的博大精深,不断演习着各路内功心法,并揣摩如何将这些心法融会贯通,再建立新的经典。1929年,颇具少侠风范的柯尔莫果洛夫获得博士学位,并携带着论文《概率论与测度论的一般理论》横空出世,首次给出了概率论公理结构,使概率论从自然哲学领域纳入严格的数学殿堂之内,终于概率论被确认为数学的一个重要分支,并为随机过程理论的创立提供了必要基础。这样的成果或许是普通人一生都难以企及的高度,然而对于彼时的柯少侠来说,却仅仅是他在学术上的起点。04游历四方的少侠博士毕业后,年仅26岁的柯少侠成为莫斯科大学数学研究所的助理研究员,开启了他成为逍遥派一代宗师的道路。如同所有身怀绝技的少侠一样,彼时的柯尔莫果洛夫也希望在学成之后去外面游历江湖,会会各路大神,一方面可以增长见识,一方面也可以检验自己的功夫是否到家。1930年1月,柯少侠开始了他人生第一次出国访学,问道于当时的数学圣地哥廷根大学与巴黎高师。不过,柯少侠也并非喜欢独来独往的江湖剑客,而是邀请了自己一生的“好基友”亚历山德洛夫和自己一同闯荡江湖。 10个月的学术访问让柯少侠和好基友看到了更加磅礴浩瀚的江湖。他们与库兰特研究极限定理,和韦尔演绎直觉逻辑,和兰道探讨函数理论,也结识了众多江湖前辈,如希尔伯特和波莱尔等人。不过即便是在外闯荡江湖,柯少侠的研究工作也没有停止。1931年,他发表了开创性著作《关于概率论的分析方法》,并在同年晋升为莫斯科大学教授,年仅 28岁。05成就一代宗师1933年,三十而立的柯尔莫果洛夫成为莫斯科大学数学研究所所长,正式执掌逍遥派。当了领导的柯大侠不仅没有躺在功劳簿上面喝茶看报,反而变得更加勤奋。柯大侠在逍遥派内部创立了概率论、数理统计、数理逻辑等教研室,开启了其创造性研究的新阶段。这个阶段持续了相当长的时期,柯大侠的研究也呈现爆发式的增长,他不仅在概率论、随机过程、数理统计、拓扑学、泛函分析等分支成就斐然,还将其数学方法成功应用到众多其它学科。柯大侠创造的大量构造方法和基本理论已成为教科书和专著中的经典内容。 身怀绝技且大公无私的顶级武林高手在江湖上常常被人们敬称为“大侠”,而柯大侠在学术圈也有一个响亮的名号“大数学家”。“好基友”亚历山德罗夫曾诙谐地说:“数学天才有敏捷型与迟缓型两种,柯尔莫果洛夫属于前者,希尔伯特属于后者”。柯大侠之所以能够享有如此高的盛誉,除了他广泛的涉猎和极高的学术成就,还在于他具有发现重要数学概念的能力和直觉。而事实上柯大侠对现实问题的洞察力更为敏捷,他能够将一些模糊的粗线条纳进有序而完备的系统中,并给出精妙的答案。柯大侠对于湍流的研究便是如此。06纷繁的流体江湖还记得我们在前面的文章中给大家介绍过流体江湖的纷争:流体剑宗的侠客们通过“走一步,假设一步”的做法将NS方程推广到湍流问题的数值求解,搭建起了湍流计算从数学到工程的桥梁,也开启了长达百年的流体剑宗的演绎。 但是,面对纷繁复杂的湍流世界,人们似乎依然一无所知。直到科学家们认识到湍流中存在的大涡生小涡的能级串现象,而关于湍流能级串最具代表性的描述便是1922年理查德森的那首小诗。 07柯大侠的另辟蹊径有了湍流能级串的定性认识,柯大侠敏锐而深刻的意识到,湍流的定量描述或许也可以统计。于是他另辟蹊径,使出了“万物皆可统计”的小无相功。柯大侠首先使出了“小无相功”的第一招:如果小涡的尺度足够小,那么它无法直接感受到各向异性的大涡,因此小尺度涡的运动特征独立于大涡,此尺度内两点间速度差的概率密度函数与空间位置无关——也就是局部均匀各向同性。 这也意味着,能级串中各级传递特征相似。且由于此尺度范围内粘性几乎不起作用,无能量耗散,因此各级能量传递率也相同,并等于最终的能量耗散率ε——这便是“柯尔莫果洛夫的第一相似性假设”。柯大侠紧接着使出了“小无相功”的第二招:对于尺度为G的流动结构,如果η<<G<<L0(其中η为耗散尺度,L0为宏观尺度),那么此尺度范围的涡不仅不受大尺度各向异性的涡的影响,也不受耗散尺度的涡的影响。因此该尺度内的含能仅取决于能量传递速率ε,与粘性也无关——即“柯尔莫果洛夫第二相似性假设”。 最后,柯大侠使出了一招平平无奇的量纲分析,在1941年提出了湍流世界最著名的-5/3幂律,并在众多学者的实验中得到了验证。当然随着湍流理论的不断发展,人们也发现了柯大侠的K41理论并非完美无缺,不过这并不妨碍它成为湍流研究史上最耀眼的一章。08荫泽众人相对于柯大侠在赫赫生辉的一生中取得的无数成就比起来,或许他在教育上的投入和付出更加让人敬佩,也引领了更多的后来人。柯尔莫果洛夫进行科学研究的特点是:几乎在他关心的所有领域,首先搭起框架,接着和他的学生们一起研究,达到深入完备的程度,并吸引大量研究人员加入,形成新的科学研究分支。莫斯科学派之所以在武林中有着神秘而强大的地位,柯大侠的引领功不可没。 对于学生,柯尔莫果洛夫则为他们创造了严谨又神圣的研究气氛。他非常擅于和学生们交往,发现适合每个人的问题和任务,并激发他们的创造力。柯大侠一生桃李满天下,比如在我们熟知的湍流研究领域,柯大侠的两位弟 子,奥布霍夫和莫宁就分别掌管着大气和海洋的湍流研究。 柯尔莫果洛夫Andrey Nikolaevich Kolmogorov作为上世纪最伟大的科学家之一,柯尔莫果洛夫还非常关心基础教育工作。他不仅亲自担任前苏联《大百科全书》数学部的主任与编委,还主持了全国奥林匹克竞赛,并为中学生编写杂志与丛书。侠之大者,为国为民,而柯大侠的行为无愧“大数学家”的称号,也无愧他一生获得的无数荣誉。后记POSTSCRIPT作为赫赫有名的“大数学家”,柯大侠行走江湖时的惊鸿一瞥便将流体力学的侠客们引入了一个新的时代。而笔者在撰写此文时,查阅了一些资料,却发现许多论述并不一致。诚然,在西方主导科学与文化输出的时代,许多成就斐然的苏联学者并未像西方学者那样为人所熟知。同时笔者也感受到,似乎每个人在谈到柯尔莫果洛夫时,都能在自己熟知的领域里发现他的身影,或许我们看到的都只是他浩瀚成就中的冰山一角。最后,谨以此文致敬流体江湖最耀眼的“过客”。-END-来源:CFD仿真区

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习 福利任务 兑换礼品
下载APP
联系我们
帮助与反馈