大风起兮云飞扬 1

大风起兮云飞扬

—漫话流动显示及纳斯方程

王 振 东

古代诗词:以流动显示来抒发情思

大风起兮云飞扬,

威加海内兮归故乡,

安得猛士兮守四方!

这是汉高祖刘邦(公元前247~前195)在击破英布军以后,回长安时,途经他的故乡沛(今江苏徐州市沛县),设宴招待家乡的故交父老,酒酣时自己击筑(古代乐器)而歌,所作慷慨豪情的《大风歌》。

《史记:高祖本纪》:“高祖(刘邦)还归,过沛、留。置酒沛宫,悉召故人父老子弟纵酒,发沛中儿得百二十人,教之歌。酒酣,高祖击筑,自为歌诗曰:大风起兮云飞扬,威加海内兮归故乡,安得猛士兮守四方!令儿皆和习之。高祖乃起舞,慷慨伤怀,泣数行下”,正是记载了这段历史。刘邦短短三句,洋洋自得,气壮山河,但并没有被胜利冲昏头脑,最后一句流露出了居安思危的忧患意识。

刘邦在这里是以“云飞扬”流动显示大气运动的物理图像,来抒发衣锦还乡、荣归故里的壮志豪情。这是历史上有名的一则典故,“大风歌”或“大风诗”的来历。之后直至现代,不少人皆仿此“歌大风、唱大风”,以表示慷慨悲歌、治国安邦的豪情壮怀。如: 汉武帝刘彻(前156~前87)也有—首以风吹白云飞,表达情感的诗《秋风辞》 秋风起兮白云飞,草木黄落兮雁南归。

兰有秀兮菊有芳,携佳人兮不能忘。

泛楼舡兮济汾河,横中流兮扬素波。

箫鼓鸣兮发棹歌,欢乐极兮哀情多。

少壮几时兮奈老何。

唐太宗李世民(599~649)《辛武功庆善宫》诗

共乐还乡宴,欢比大风诗。

《过旧宅二首》之二

八表文同轨,无劳歌大风

李白《登广武古战场怀古》诗

按剑清八极,归酣歌大风

林宽《歌风台》诗

蒿棘空存百尺基,酒酣曾唱大风词

王德贞《奉和圣制过温汤》诗

停舆兴睿览,还举大风篇

直到近代也有类似的大风诗,如:

董必武(1885~1975)《感时杂咏》诗

欲守四方歌大风,飞鸟未尽先藏弓。

朱德(1886~1976)《赠友人》诗

北华收复赖群雄,猛士如云唱大风。

陈毅(1901~1972)《莱芜大捷》诗

鲁中霁雪明飞帜,渤海洪波唱大风。

现在以云来显示大气的流动,己很常见。如在电视台的气象预报节目中,人们常能看到由云显示千姿百态流动图案的卫星云图,所显示大气中所发生的动力过程。

也有古诗用风叶和船只所显示的流体运动,来形象、生动地比喻和描述远行在外人的行迹和旅途。如宋代诗人范成大(1126~1193)的五言律诗《道中》

月冷吟蛩草,湖平宿鹭沙。客愁无锦字,乡信有灯花。

踪迹随风叶,程途犯斗槎。君看枝上鹊,薄暮亦还家。

程途是指旅程途中,槎(chá)亦做查、楂,系水中木筏意,犯斗槎是指远行所乘的船只。

古代诗人还常以杨絮、柳絮以及虫类拉的丝(亦名游丝、晴丝),所显示的空气流动情况(风、对流或布朗运动),来抒发各种各样的情思,如:

韩愈(768~824)《晚雪》诗

杨花榆荚无才思,唯解漫天作雪飞。

以及《次同冠峡》诗

落英千尺堕,游丝百丈飘。

周紫芝《踏莎行》词

情似游丝,人如飞絮,泪珠阁定空相觑。

范成大《碧瓦》诗

无风杨柳漫天絮,不雨棠梨满地花。

以及《初夏二首》诗

晴丝千尺挽韶光,百舌无声燕子忙。

韶光是美好的时光,这里指春天。诗人想象春末夏初的游丝是在恋惜时光,想把春天挽留住。

石矛《绝句》诗

来时万缕弄轻黄,去日飞毬满路旁。

我比杨花更飘荡,杨花只是—春忙。

以杨花比喻自己奔波游宦,道出了深沉的乡思旅愁。

苏轼《水龙吟·和章质夫杨花韵》词

似花还似非花,也无人惜从教坠,抛家旁路,思量却是,无情有思。

将杨花比作缠绵衰感的思妇。

文天祥《过零丁洋》诗

山河破碎风抛絮,身世飘摇雨打萍

把杨花比作日益沦丧的国土。

各种各样的流动显示方法

流动显示是在力求不改变流体运动性质的前提下,用图像显示流体运动的方法,其任务

是使流体不可见的流动特征,成为可见的。俗话说“百闻不如—见”,人们通过流动显示看到了流场的特征,从而可进一步研究探索和应用流体运动的规律。

西方一些人认为,意大利文艺复兴时期的艺术家和科学家达•芬奇(Da.Vinci,1452~1519),是第一个运用流动显示的方法,来叙述涡旋构图的人。但比起运用流动显示的图像,来描述峡江水流涡旋的运动特征,和抒情言志的我国古代诗人,达•芬奇却要落后好几个世纪了。

首先应用流动显示方法,对现代流体力学发展做出重要贡献,当推英国科学家雷诺(O.Reynolds,1842~1912)。他在1883年,将苯胺染液注入长的水平管道水流中做示踪剂,从而可以看出管中水的流动状态。当流速小时,苯胺染液形成一根纤细的直线与管轴平行,表示流动是稳定的和有规则的流动,称为层流;当流速慢慢地增加,达到某一数值时,流动形式突然发生变化,那根苯胺染液细线受到激烈的扰动,苯胺染液迅速地散布于整个管内,表示流动己十分紊乱,称为湍流。这一试验明确提出了两种不同的流动状态,及其转捩的概念,还提出了后来被称为“雷诺数”的这一十分重要的无量纲参数。至今湍流研究的历史,一般都公认从1883年雷诺这个经典的流动实验算起。

德国科学家普朗特(L.Prandtl,1875~1953),1904年用在水中撤放粒子的方法,获得了水沿薄平板运动的画面。由于画面上粒子留下的轨迹正比于流动的速度,在靠近壁面有一薄层,其中速度比离壁面较远处的速度明显较小,且有大的速度梯度。正是对这一流动显示画面的观察和分折,使他提出了边界层的概念,指出在远离壁面处,可不计黏性,能应用理想流体力学的研究结果;而在物体表面附近的薄层中,由于有很大的速度梯度,从而产生很大的剪切力,不能忽略黏性。这一基于流动显示的新观点,使得可利用边界层很薄的特点,使问题的数学处理大为简化,至今它仍是黏性流体力学最重要的基础理论之一。

20世纪50年代,有人提出了氢气泡显示技术:用很细的金属丝放在水中作为阴极,通电后在金属丝上形成的氢气泡随水流走,而成为显示流场的示踪粒子。克拉茵(Kline)等1967年首先用氢气泡显示技术,发现了近壁湍流的相干结构(Coherent Structure,也有译为拟序结构)。这是一种大尺度的涡旋运动,它在将平均运动动能转变为湍流动能的过程中,作了大部分贡献。后来经许多人用更精确、先进的实验手段(热线热膜测速、激光测速以及数据自动采集、图像处理技术等)进行重复,使实验越做越精确。不但对壁湍流,而且对自由剪切湍流也发现了相干结构,到20世纪80年代,湍流相干结构已为国际流体力学界公认,并认为这是对湍流生成、维持、演化起主要作用的结构。这一由流动显示所发现的相干结构,被认为是对湍流认识上的一次革命,是在湍流研究上的一次重大进展。80年代之后至今,关于湍流相干结构及其控制的研究,一直湍流研究的热点课题。

由以上三个例子可见,流动显示是了解流体运动特性,并深入探索其物理机制的一种直观、有效的手段。它能发现新的流动现象,如层流和湍流两种流动状态及其转捩、涡旋、分离、激波、边界层、壁湍流相干结构等;据了解,流动显示技术己在许多实际问题的研究中,发挥了很大的作用,如三角翼和双三角翼的前缘主涡、二次涡和尾涡的形成和发展,钝物体尾迹的涡旋结构,以及多体干扰等。

上面提到的流动显示方法,,主要只涉及到示踪法。示踪法是在流体中加入某些示踪物质,通过对加入物质踪迹观察得到流体运动的图像。由于所加示踪物质的不同,又可分为用途不一的烟迹(含烟丝)法、染色线法、空气泡和氢气泡法、氦气泡法、激光-荧光法、蒸汽屏法等。当然,在流体中加入了示踪粒子,就又存在粒子的跟随性问题。

除示踪法外,流动显示的方法还有光学方法和表面涂料显迹法。光学方法又分阴影法、

纹影法和干涉法。前两者利用了光通过非均匀流场不同部位时的折射效应,后者通过扰动光和未扰动光的相互干涉得到干涉条纹图,从而进一步可得到流动参数的定量结果。表面涂料显迹法是在物面上涂以薄层物质,以其与流动相互作用时,产生一定的可见图像,从而可定性或定量的推断物面附近的流动特性。按所涂物质的不同,还可分为油流(荧光油流)、丝线(荧光微丝)、染料、升华、相变涂层、液晶、感温漆等方法。

流动显示技术目前发展相当快,特别是与计算机图像处理技术相结合,使传统的流动显示方法得到很大的改进。计算机数据的采集与处理,可对显示结果进行深度的加工分析,以获得更清晰的流动图像,以及有关流动参数的分布。

多种流动显示方法的联合使用,又可得到更丰富的流动信息。随着光学技术和计算机技术的发展,激光全息术、光学层析术、散斑、粒子成像测速(PIV—Particale Image Velocimetry)、激光诱导荧光(LIF—Laser Induce Fluorescent)等方法也己出现并在发展完善之中,为实现瞬时、高分辨率和定量化的空间流动显示展现了美好的前景。

数值模拟、实验检验和世纪数学难题

要弄清流动显示对流体力学的研究能有多大的作用,还需要从流体力学的研究现状来说起。

力学是以实验为基础的科学,流体力学更是建立在实验的基础之上。在流体力学中,绝大多数重要的概念和原理都源于实验,例如:大气压强,流体的可压缩性,黏性剪应力,层流,湍流,雷诺数,卡门涡,二次流,附加质量,激波,孤立波,湍剪切流的相干结构,声障现象等;又如,完全气体的状态方程,连续性方程,能量守恒原理,达西定律,托里拆利原理,伯努利原理等。

瑞士数学家、力学家欧拉(Euler,L. 1707~1783)于1755年,建立了理想流体的动力学方程组,现称为欧拉方程组。法国力学家、工程师纳维(Navier,C.L.M.H. 1785~1836)于1821年,以及英国力学家、数学家斯托克斯(Stokes,G.G. 1819~1903)于1845年,分别对黏性不可压缩流体建立了动力学方程组,现称为纳维—斯托克斯方程组。在无黏性的情况下,纳维―斯托克斯方程组可简化为欧拉方程组。现在人们对于自然界、国防和各种工程技术中的流体力学问题,都在用纳维―斯托克斯方程组进行分析、计算和研究。纳维―斯托克斯方程组(亦可简称为:纳斯方程),现被公认是描述流体运动规律的流体力学基本方程组。

对于纳维—斯托克斯方程组,经过150多年的研究,仅在—些简化的特殊情况下,找到不多的准确解。由于纳维—斯托克斯方程组光滑解的存在性问题,至今尚没有在数学上解决,且这个问题又关系到人类的生产、生活、军事和对大自然的认识,极其重要,所以克莱数学促进会(Clay Mathematics Institute )于2000年5月24日在法国巴黎的法兰西学院,将其发布为新千年数学大奖悬赏的7个世纪数学难题之—,悬赏奖金高达一百万美元。克莱数学促进会发布的7个世纪数学难题是:P与NP问题、黎曼(Riemann)假设、庞加莱(Poincaré)猜想、霍奇(Hodge)猜想、贝尔什和斯威尔顿(Birch及Swinnerton-Dyer)猜想、纳维―斯托克斯方程、杨―米尔斯(Yang-Mills)理论。比纳维—斯托克斯方程组简单得多的欧拉方程组,解的存在性的问题也尚未得到证明,只是它不属于悬赏奖励的问题内容。

在学习微分方程理论時,我们知道:

(1)如果某物理问题的微分方程,被证明其解不仅存在而且唯一时,则无论用何种方法找到这个微分方程的解,可以认为这就是该物理问题方程的解。

(2)当某物理问题的微分方程,被证明解是存在的,但却不见得唯一时,则如用—种方法

找到了解,还必须研究解的稳定性问题,只有证明了所找到的解是稳定的,才能认为这个解有可能代表实际存在的物理现象。

(3)如果某物理问题的微分方程,解的存在性尚还不能被证明,若用某种近似方法(如渐近方法或差分法、有限元法等各种数值方法)找到了“解”,则难以肯定它是否真是代表实际存在的物理现象的解。

不幸的是,流体力学中所遇到的欧拉方程组和纳维—斯托克斯方程组,正好都属于第三种情况。

如果经过数学家的努力,解决了悬赏的问题,纳维—斯托克斯方程组解的存在性问题得到了证明,这自然是皆大欢喜的好事。可是关于纳维—斯托克斯方程组解的存在性问题的悬赏,也还包括给出其解不存在的证明。如果是后者获奖,那问题就大了。当然也有这种可能,经过仔细研究后认为纳维—斯托克斯方程组应做出某些修正和改进,才能使解存在。如是这样,流体力学教科书就需要改写了。

可是,大量的自然界、国防和各种工程实际中的流体力学问题需要解决,并不能等你弄清方程组解的存在性后再说。人们只能在用理论分析、数值计算、物理实验相结合的方法,研究、解决所遇到的流体力学问题。

这三种方法各有优缺点。实验方法的优点是能直接解决生产中的复杂问题,能发现流动中的新现象和新原理,其结果可作为捡验其他方法是否正确的依据;缺点是对不同情况需做不同的实验,且所需人力、财力、物力较多,花费大。分析方法的优点是可明确给出各物理量与流动参数之间的变化关系,普适性较好;缺点是数学上的困难很大,能获得的分析解(包括近似的分析解)的数量有限。数值计算方法的优点是可对分析法无法求解的问题,求得其数值解,且花费相对较小;缺点是对复杂而又缺乏完善数学模型的问题,仍无能为力。分析解及数值解都是建立在具有—定假设条件的运动方程组之上的,其结果仍都应受到物理实验结果的捡验。由于纳维—斯托克斯方程组解的存在性问题至今尚未解决,就更难以肯定数值方法找到的解,是否代表真实的流体运动。所以,数值摸拟与物理实验的本质差别并未消失,数值模拟尚不能替代物理实验,数值摸拟的结果必须用物理实验来捡验其正确性。

由于计算机和数值计算技术的快速发展,出于科学研究和生产实际的需要,对于流体力学问题进行大规模数值模拟,现己很常见,国内已有几种功能较强的计算流体动力学的商品软件(如 FLUENT, STAR—CD, TASC flow,PHOENICS 等)在应用,且已使用并行计算机进行大规模数值模拟。但所得到的数值模拟结果,仍须用物理实验来检验其正确性。而作物理实验又需要投入更多的人力、财力、物力的支持,所以巧妙地构思、设计小规模、精细的物理实验,以较少的花费来捡验大规模数值模拟的正确性,就显得十分重要。

流动显示方法和技术,正是我们在流体力学研究中,能达到上述目的的重要实验方法和技术,它不仅能提出新的观念、新的研究模型,揭示流体运动规律,也能为流体力学计算提供可靠的流动条件(如边界层转捩点、激波位置、涡核位置、尾迹宽度等),和对数值模拟的结果进行检验。

附录:新千年数学大奖悬赏的7个世纪数学难题

Notices of the AMS(美国数学会(AMS)的会刊)在克莱数学促进会发布7个世纪数学难题后,曾为悬赏问题准备了如下的简介:

P和NP问题:一个问题称为是P的,如果它可以通过运行多项式次(即运行时间至多是输入量大小的多项式函数)的一种算法获得解决;一个问题称为是NP的,如果所提出的解答,可以用多项式次算法来检验。P等于NP吗?

Riemann假设:黎曼ζ函数的每个非平凡零点,有等于1/2的实部。

Poincaré猜想:任何单连通闭3维流形同胚于3维球。

Hodge猜想:任何霍奇类关于一个非奇异复射影代数簇,都是某些代数闭链类的有理线性组合。

Birch 及Swinnerton–Dyer猜想:对于建立在有理数域上的每一条椭圆曲线,它在1处的L函数变为零的阶,等于该曲线上有理点的阿贝尔群的秩。

Navier–Stokes方程组:(在适当的边界及初始条件下)对3维纳维–斯托克斯方程组,证明或反证其光滑解的存在性。

Yang–Mills理论:证明量子杨–米尔斯场存在,并存在一个质量间隙。

参考文献

1、王振东,但见流沫生千涡 —漫话流体中的旋涡,科学,2004,56(2):59~61

2、王振东,峡江漱石水多漩 —漫谈流体中的旋涡,力学与实践,1995,17(2):68~70

3、王振东,踪迹随风叶,程途犯斗搓 —漫谈流动显示及应用,力学与实践,1996,18(2):

75~77

4、 王振东、姜楠,新千年数学大奖问题 ——证明纳维—斯托克斯方程组光滑解的存在性,

力学与实践,2003,25(3):72—73

5、王振东,关于流体力学方法论问题,力学与实践,2004,26(2):83~85

6、王振东、武际可,力学诗趣,天津:南开大学出版社,1998

7、曹道衡选注,两汉诗选,北京:中华书局,2005

8、缪斌、霍松林等编,宋诗鉴赏辞典,上海:上海辞书出版社,1987

9、徐成志、王光汉等编,常用典故词典,上海:上海辞书出版社,1985

10、Allyn Jackson ,Million-dollar Mathematics Prizes Announced ,Notices of the AMS,2000,47

(8):877~879

11、Wenjei Yang ,Handbook of Flow Visualization ,Hemispere Publishing Corporation,1980

12、Smits AJ,Lim TT , Flow Visualization ,Imperial College Press 2000

关键词 流动显示 , 数值模拟 , 实验检验 , 纳维‐斯托克斯方程

附图

云显示的美国2005年卡特里娜飓风

尘物显示的龙卷风

液晶显示汽轮机叶片附近的湍流

烟迹显示的流动图

染色法显示的卡门涡街

(原刊登于《自然杂志》2006年28卷5期)

大风起兮云飞扬

—漫话流动显示及纳斯方程

王 振 东

古代诗词:以流动显示来抒发情思

大风起兮云飞扬,

威加海内兮归故乡,

安得猛士兮守四方!

这是汉高祖刘邦(公元前247~前195)在击破英布军以后,回长安时,途经他的故乡沛(今江苏徐州市沛县),设宴招待家乡的故交父老,酒酣时自己击筑(古代乐器)而歌,所作慷慨豪情的《大风歌》。

《史记:高祖本纪》:“高祖(刘邦)还归,过沛、留。置酒沛宫,悉召故人父老子弟纵酒,发沛中儿得百二十人,教之歌。酒酣,高祖击筑,自为歌诗曰:大风起兮云飞扬,威加海内兮归故乡,安得猛士兮守四方!令儿皆和习之。高祖乃起舞,慷慨伤怀,泣数行下”,正是记载了这段历史。刘邦短短三句,洋洋自得,气壮山河,但并没有被胜利冲昏头脑,最后一句流露出了居安思危的忧患意识。

刘邦在这里是以“云飞扬”流动显示大气运动的物理图像,来抒发衣锦还乡、荣归故里的壮志豪情。这是历史上有名的一则典故,“大风歌”或“大风诗”的来历。之后直至现代,不少人皆仿此“歌大风、唱大风”,以表示慷慨悲歌、治国安邦的豪情壮怀。如: 汉武帝刘彻(前156~前87)也有—首以风吹白云飞,表达情感的诗《秋风辞》 秋风起兮白云飞,草木黄落兮雁南归。

兰有秀兮菊有芳,携佳人兮不能忘。

泛楼舡兮济汾河,横中流兮扬素波。

箫鼓鸣兮发棹歌,欢乐极兮哀情多。

少壮几时兮奈老何。

唐太宗李世民(599~649)《辛武功庆善宫》诗

共乐还乡宴,欢比大风诗。

《过旧宅二首》之二

八表文同轨,无劳歌大风

李白《登广武古战场怀古》诗

按剑清八极,归酣歌大风

林宽《歌风台》诗

蒿棘空存百尺基,酒酣曾唱大风词

王德贞《奉和圣制过温汤》诗

停舆兴睿览,还举大风篇

直到近代也有类似的大风诗,如:

董必武(1885~1975)《感时杂咏》诗

欲守四方歌大风,飞鸟未尽先藏弓。

朱德(1886~1976)《赠友人》诗

北华收复赖群雄,猛士如云唱大风。

陈毅(1901~1972)《莱芜大捷》诗

鲁中霁雪明飞帜,渤海洪波唱大风。

现在以云来显示大气的流动,己很常见。如在电视台的气象预报节目中,人们常能看到由云显示千姿百态流动图案的卫星云图,所显示大气中所发生的动力过程。

也有古诗用风叶和船只所显示的流体运动,来形象、生动地比喻和描述远行在外人的行迹和旅途。如宋代诗人范成大(1126~1193)的五言律诗《道中》

月冷吟蛩草,湖平宿鹭沙。客愁无锦字,乡信有灯花。

踪迹随风叶,程途犯斗槎。君看枝上鹊,薄暮亦还家。

程途是指旅程途中,槎(chá)亦做查、楂,系水中木筏意,犯斗槎是指远行所乘的船只。

古代诗人还常以杨絮、柳絮以及虫类拉的丝(亦名游丝、晴丝),所显示的空气流动情况(风、对流或布朗运动),来抒发各种各样的情思,如:

韩愈(768~824)《晚雪》诗

杨花榆荚无才思,唯解漫天作雪飞。

以及《次同冠峡》诗

落英千尺堕,游丝百丈飘。

周紫芝《踏莎行》词

情似游丝,人如飞絮,泪珠阁定空相觑。

范成大《碧瓦》诗

无风杨柳漫天絮,不雨棠梨满地花。

以及《初夏二首》诗

晴丝千尺挽韶光,百舌无声燕子忙。

韶光是美好的时光,这里指春天。诗人想象春末夏初的游丝是在恋惜时光,想把春天挽留住。

石矛《绝句》诗

来时万缕弄轻黄,去日飞毬满路旁。

我比杨花更飘荡,杨花只是—春忙。

以杨花比喻自己奔波游宦,道出了深沉的乡思旅愁。

苏轼《水龙吟·和章质夫杨花韵》词

似花还似非花,也无人惜从教坠,抛家旁路,思量却是,无情有思。

将杨花比作缠绵衰感的思妇。

文天祥《过零丁洋》诗

山河破碎风抛絮,身世飘摇雨打萍

把杨花比作日益沦丧的国土。

各种各样的流动显示方法

流动显示是在力求不改变流体运动性质的前提下,用图像显示流体运动的方法,其任务

是使流体不可见的流动特征,成为可见的。俗话说“百闻不如—见”,人们通过流动显示看到了流场的特征,从而可进一步研究探索和应用流体运动的规律。

西方一些人认为,意大利文艺复兴时期的艺术家和科学家达•芬奇(Da.Vinci,1452~1519),是第一个运用流动显示的方法,来叙述涡旋构图的人。但比起运用流动显示的图像,来描述峡江水流涡旋的运动特征,和抒情言志的我国古代诗人,达•芬奇却要落后好几个世纪了。

首先应用流动显示方法,对现代流体力学发展做出重要贡献,当推英国科学家雷诺(O.Reynolds,1842~1912)。他在1883年,将苯胺染液注入长的水平管道水流中做示踪剂,从而可以看出管中水的流动状态。当流速小时,苯胺染液形成一根纤细的直线与管轴平行,表示流动是稳定的和有规则的流动,称为层流;当流速慢慢地增加,达到某一数值时,流动形式突然发生变化,那根苯胺染液细线受到激烈的扰动,苯胺染液迅速地散布于整个管内,表示流动己十分紊乱,称为湍流。这一试验明确提出了两种不同的流动状态,及其转捩的概念,还提出了后来被称为“雷诺数”的这一十分重要的无量纲参数。至今湍流研究的历史,一般都公认从1883年雷诺这个经典的流动实验算起。

德国科学家普朗特(L.Prandtl,1875~1953),1904年用在水中撤放粒子的方法,获得了水沿薄平板运动的画面。由于画面上粒子留下的轨迹正比于流动的速度,在靠近壁面有一薄层,其中速度比离壁面较远处的速度明显较小,且有大的速度梯度。正是对这一流动显示画面的观察和分折,使他提出了边界层的概念,指出在远离壁面处,可不计黏性,能应用理想流体力学的研究结果;而在物体表面附近的薄层中,由于有很大的速度梯度,从而产生很大的剪切力,不能忽略黏性。这一基于流动显示的新观点,使得可利用边界层很薄的特点,使问题的数学处理大为简化,至今它仍是黏性流体力学最重要的基础理论之一。

20世纪50年代,有人提出了氢气泡显示技术:用很细的金属丝放在水中作为阴极,通电后在金属丝上形成的氢气泡随水流走,而成为显示流场的示踪粒子。克拉茵(Kline)等1967年首先用氢气泡显示技术,发现了近壁湍流的相干结构(Coherent Structure,也有译为拟序结构)。这是一种大尺度的涡旋运动,它在将平均运动动能转变为湍流动能的过程中,作了大部分贡献。后来经许多人用更精确、先进的实验手段(热线热膜测速、激光测速以及数据自动采集、图像处理技术等)进行重复,使实验越做越精确。不但对壁湍流,而且对自由剪切湍流也发现了相干结构,到20世纪80年代,湍流相干结构已为国际流体力学界公认,并认为这是对湍流生成、维持、演化起主要作用的结构。这一由流动显示所发现的相干结构,被认为是对湍流认识上的一次革命,是在湍流研究上的一次重大进展。80年代之后至今,关于湍流相干结构及其控制的研究,一直湍流研究的热点课题。

由以上三个例子可见,流动显示是了解流体运动特性,并深入探索其物理机制的一种直观、有效的手段。它能发现新的流动现象,如层流和湍流两种流动状态及其转捩、涡旋、分离、激波、边界层、壁湍流相干结构等;据了解,流动显示技术己在许多实际问题的研究中,发挥了很大的作用,如三角翼和双三角翼的前缘主涡、二次涡和尾涡的形成和发展,钝物体尾迹的涡旋结构,以及多体干扰等。

上面提到的流动显示方法,,主要只涉及到示踪法。示踪法是在流体中加入某些示踪物质,通过对加入物质踪迹观察得到流体运动的图像。由于所加示踪物质的不同,又可分为用途不一的烟迹(含烟丝)法、染色线法、空气泡和氢气泡法、氦气泡法、激光-荧光法、蒸汽屏法等。当然,在流体中加入了示踪粒子,就又存在粒子的跟随性问题。

除示踪法外,流动显示的方法还有光学方法和表面涂料显迹法。光学方法又分阴影法、

纹影法和干涉法。前两者利用了光通过非均匀流场不同部位时的折射效应,后者通过扰动光和未扰动光的相互干涉得到干涉条纹图,从而进一步可得到流动参数的定量结果。表面涂料显迹法是在物面上涂以薄层物质,以其与流动相互作用时,产生一定的可见图像,从而可定性或定量的推断物面附近的流动特性。按所涂物质的不同,还可分为油流(荧光油流)、丝线(荧光微丝)、染料、升华、相变涂层、液晶、感温漆等方法。

流动显示技术目前发展相当快,特别是与计算机图像处理技术相结合,使传统的流动显示方法得到很大的改进。计算机数据的采集与处理,可对显示结果进行深度的加工分析,以获得更清晰的流动图像,以及有关流动参数的分布。

多种流动显示方法的联合使用,又可得到更丰富的流动信息。随着光学技术和计算机技术的发展,激光全息术、光学层析术、散斑、粒子成像测速(PIV—Particale Image Velocimetry)、激光诱导荧光(LIF—Laser Induce Fluorescent)等方法也己出现并在发展完善之中,为实现瞬时、高分辨率和定量化的空间流动显示展现了美好的前景。

数值模拟、实验检验和世纪数学难题

要弄清流动显示对流体力学的研究能有多大的作用,还需要从流体力学的研究现状来说起。

力学是以实验为基础的科学,流体力学更是建立在实验的基础之上。在流体力学中,绝大多数重要的概念和原理都源于实验,例如:大气压强,流体的可压缩性,黏性剪应力,层流,湍流,雷诺数,卡门涡,二次流,附加质量,激波,孤立波,湍剪切流的相干结构,声障现象等;又如,完全气体的状态方程,连续性方程,能量守恒原理,达西定律,托里拆利原理,伯努利原理等。

瑞士数学家、力学家欧拉(Euler,L. 1707~1783)于1755年,建立了理想流体的动力学方程组,现称为欧拉方程组。法国力学家、工程师纳维(Navier,C.L.M.H. 1785~1836)于1821年,以及英国力学家、数学家斯托克斯(Stokes,G.G. 1819~1903)于1845年,分别对黏性不可压缩流体建立了动力学方程组,现称为纳维—斯托克斯方程组。在无黏性的情况下,纳维―斯托克斯方程组可简化为欧拉方程组。现在人们对于自然界、国防和各种工程技术中的流体力学问题,都在用纳维―斯托克斯方程组进行分析、计算和研究。纳维―斯托克斯方程组(亦可简称为:纳斯方程),现被公认是描述流体运动规律的流体力学基本方程组。

对于纳维—斯托克斯方程组,经过150多年的研究,仅在—些简化的特殊情况下,找到不多的准确解。由于纳维—斯托克斯方程组光滑解的存在性问题,至今尚没有在数学上解决,且这个问题又关系到人类的生产、生活、军事和对大自然的认识,极其重要,所以克莱数学促进会(Clay Mathematics Institute )于2000年5月24日在法国巴黎的法兰西学院,将其发布为新千年数学大奖悬赏的7个世纪数学难题之—,悬赏奖金高达一百万美元。克莱数学促进会发布的7个世纪数学难题是:P与NP问题、黎曼(Riemann)假设、庞加莱(Poincaré)猜想、霍奇(Hodge)猜想、贝尔什和斯威尔顿(Birch及Swinnerton-Dyer)猜想、纳维―斯托克斯方程、杨―米尔斯(Yang-Mills)理论。比纳维—斯托克斯方程组简单得多的欧拉方程组,解的存在性的问题也尚未得到证明,只是它不属于悬赏奖励的问题内容。

在学习微分方程理论時,我们知道:

(1)如果某物理问题的微分方程,被证明其解不仅存在而且唯一时,则无论用何种方法找到这个微分方程的解,可以认为这就是该物理问题方程的解。

(2)当某物理问题的微分方程,被证明解是存在的,但却不见得唯一时,则如用—种方法

找到了解,还必须研究解的稳定性问题,只有证明了所找到的解是稳定的,才能认为这个解有可能代表实际存在的物理现象。

(3)如果某物理问题的微分方程,解的存在性尚还不能被证明,若用某种近似方法(如渐近方法或差分法、有限元法等各种数值方法)找到了“解”,则难以肯定它是否真是代表实际存在的物理现象的解。

不幸的是,流体力学中所遇到的欧拉方程组和纳维—斯托克斯方程组,正好都属于第三种情况。

如果经过数学家的努力,解决了悬赏的问题,纳维—斯托克斯方程组解的存在性问题得到了证明,这自然是皆大欢喜的好事。可是关于纳维—斯托克斯方程组解的存在性问题的悬赏,也还包括给出其解不存在的证明。如果是后者获奖,那问题就大了。当然也有这种可能,经过仔细研究后认为纳维—斯托克斯方程组应做出某些修正和改进,才能使解存在。如是这样,流体力学教科书就需要改写了。

可是,大量的自然界、国防和各种工程实际中的流体力学问题需要解决,并不能等你弄清方程组解的存在性后再说。人们只能在用理论分析、数值计算、物理实验相结合的方法,研究、解决所遇到的流体力学问题。

这三种方法各有优缺点。实验方法的优点是能直接解决生产中的复杂问题,能发现流动中的新现象和新原理,其结果可作为捡验其他方法是否正确的依据;缺点是对不同情况需做不同的实验,且所需人力、财力、物力较多,花费大。分析方法的优点是可明确给出各物理量与流动参数之间的变化关系,普适性较好;缺点是数学上的困难很大,能获得的分析解(包括近似的分析解)的数量有限。数值计算方法的优点是可对分析法无法求解的问题,求得其数值解,且花费相对较小;缺点是对复杂而又缺乏完善数学模型的问题,仍无能为力。分析解及数值解都是建立在具有—定假设条件的运动方程组之上的,其结果仍都应受到物理实验结果的捡验。由于纳维—斯托克斯方程组解的存在性问题至今尚未解决,就更难以肯定数值方法找到的解,是否代表真实的流体运动。所以,数值摸拟与物理实验的本质差别并未消失,数值模拟尚不能替代物理实验,数值摸拟的结果必须用物理实验来捡验其正确性。

由于计算机和数值计算技术的快速发展,出于科学研究和生产实际的需要,对于流体力学问题进行大规模数值模拟,现己很常见,国内已有几种功能较强的计算流体动力学的商品软件(如 FLUENT, STAR—CD, TASC flow,PHOENICS 等)在应用,且已使用并行计算机进行大规模数值模拟。但所得到的数值模拟结果,仍须用物理实验来检验其正确性。而作物理实验又需要投入更多的人力、财力、物力的支持,所以巧妙地构思、设计小规模、精细的物理实验,以较少的花费来捡验大规模数值模拟的正确性,就显得十分重要。

流动显示方法和技术,正是我们在流体力学研究中,能达到上述目的的重要实验方法和技术,它不仅能提出新的观念、新的研究模型,揭示流体运动规律,也能为流体力学计算提供可靠的流动条件(如边界层转捩点、激波位置、涡核位置、尾迹宽度等),和对数值模拟的结果进行检验。

附录:新千年数学大奖悬赏的7个世纪数学难题

Notices of the AMS(美国数学会(AMS)的会刊)在克莱数学促进会发布7个世纪数学难题后,曾为悬赏问题准备了如下的简介:

P和NP问题:一个问题称为是P的,如果它可以通过运行多项式次(即运行时间至多是输入量大小的多项式函数)的一种算法获得解决;一个问题称为是NP的,如果所提出的解答,可以用多项式次算法来检验。P等于NP吗?

Riemann假设:黎曼ζ函数的每个非平凡零点,有等于1/2的实部。

Poincaré猜想:任何单连通闭3维流形同胚于3维球。

Hodge猜想:任何霍奇类关于一个非奇异复射影代数簇,都是某些代数闭链类的有理线性组合。

Birch 及Swinnerton–Dyer猜想:对于建立在有理数域上的每一条椭圆曲线,它在1处的L函数变为零的阶,等于该曲线上有理点的阿贝尔群的秩。

Navier–Stokes方程组:(在适当的边界及初始条件下)对3维纳维–斯托克斯方程组,证明或反证其光滑解的存在性。

Yang–Mills理论:证明量子杨–米尔斯场存在,并存在一个质量间隙。

参考文献

1、王振东,但见流沫生千涡 —漫话流体中的旋涡,科学,2004,56(2):59~61

2、王振东,峡江漱石水多漩 —漫谈流体中的旋涡,力学与实践,1995,17(2):68~70

3、王振东,踪迹随风叶,程途犯斗搓 —漫谈流动显示及应用,力学与实践,1996,18(2):

75~77

4、 王振东、姜楠,新千年数学大奖问题 ——证明纳维—斯托克斯方程组光滑解的存在性,

力学与实践,2003,25(3):72—73

5、王振东,关于流体力学方法论问题,力学与实践,2004,26(2):83~85

6、王振东、武际可,力学诗趣,天津:南开大学出版社,1998

7、曹道衡选注,两汉诗选,北京:中华书局,2005

8、缪斌、霍松林等编,宋诗鉴赏辞典,上海:上海辞书出版社,1987

9、徐成志、王光汉等编,常用典故词典,上海:上海辞书出版社,1985

10、Allyn Jackson ,Million-dollar Mathematics Prizes Announced ,Notices of the AMS,2000,47

(8):877~879

11、Wenjei Yang ,Handbook of Flow Visualization ,Hemispere Publishing Corporation,1980

12、Smits AJ,Lim TT , Flow Visualization ,Imperial College Press 2000

关键词 流动显示 , 数值模拟 , 实验检验 , 纳维‐斯托克斯方程

附图

云显示的美国2005年卡特里娜飓风

尘物显示的龙卷风

液晶显示汽轮机叶片附近的湍流

烟迹显示的流动图

染色法显示的卡门涡街

(原刊登于《自然杂志》2006年28卷5期)


相关文章

  • 大风起兮云飞扬--刘邦[大风歌]赏析
  • 英国的历史学家约瑟.汤恩比在评论世界历史人物时说:"人类历史上最有远见.对后世影响最大的有两位政治人物,一位是开创罗马帝国的恺撒,另一位便是创建大汉文明的汉太祖刘邦." 这个评价显然是够高的了.但平心而论,这不是无原则的 ...查看


  • 大风起兮云飞扬 2
  • 大风起兮云风扬 --论传统文化的命运 从鸦片战争开始,古老的中国第一次遇到了来自西方资本主义侵略势力的冲击,东西方两种社会制度的碰撞,引发了中国传统文化与西方近代资本主义文化的冲突.历史第一次向中国人提出中国文化向何处去的问题.随着时代的前 ...查看


  • 大风起兮云飞扬
  • 经过紧张的筹备,四川博物院藏张大千书画作品在甘肃省博物馆展出了.四川是张大千的故乡,甘肃是他曾经留恋徘徊的艺术之乡.70余年之后,张大千的作品在他曾经为艺术奋斗过的地方展出,不觉让人慨叹岁月易逝,人间沧桑. 大风起兮云飞扬. 威加海内兮归故 ...查看


  • 关于风的诗句:[大风歌]
  • 出自西汉诗人刘邦的<大风歌> 大风起兮云飞扬,威加海内兮归故乡,安得猛士兮守四方. 赏析 刘邦得以战胜项羽,是依靠许多支军队的协同作战.这些军队,有的是他的盟军,本无统属关系;有的虽然原是他的部属,但由于在战争中实力迅速增强,已 ...查看


  • 戏说[大风歌]
  • 此诗作于汉十二年(公元前195年),其时刘邦征讨英布叛军,胜局已定,返回长安时途经故乡沛县,与父老.故人.诸母置酒高会.据<汉书・高帝纪第一下>记载,刘邦击筑自歌,令儿皆和习之,"乃起舞,慷慨伤怀,泣数行下" ...查看


  • 水利工程施工现场扬尘污染防治方案
  • 河南省前坪水库主体工程第一批项目 主体施工4标 (合同编号:QPSK-ZTSG-04) 施工现场扬尘污染防治方案 批准: 审核: 编制: 河南省水利第二工程局前坪水库 主体工程第一批项目主体施工4标项目部 二零一六年七月十二日 施工现场扬尘 ...查看


  • 鸿门宴资料
  • 刘邦入关约法三章 秦王子婴降.沛公入秦宫,意欲留居之.樊哙.张良谏,乃封秦重宝财物府库,还军霸上.召诸县父老豪杰曰:"吾与父老约法三章耳:杀人者死,伤人及盗抵罪."秦人大喜,争持牛羊酒食献飨军士.沛公不受.秦人又益喜,唯 ...查看


  • 北京五月天气众生相:裙角飞扬 湿冷寒凉
  • 2016-05-12 16:04 千龙网 导语 5月11日,北京室外气温接近30℃,12日的最低气温又跌至9℃;昨日还是阳光明媚酷热难当,一夜阵雨又把市民们拉回了清冷的春日.有市民不禁感叹,冬天的衣服刚刚收起,这一降温冷的叫人无处躲藏.图像 ...查看


  • 神话故事:智搏"大风"
  • 羿的第四个目标是去东方的青丘之泽,捕杀一只叫"大风"的猛禽. "大风"就是"大凤",也就是孔雀.这只大风,性极凶悍,身体特大,一展双翅能遮住半边天,双翅一肩动就刮起了飓风,大树被连 ...查看


热门内容