华院计算 | 纳维、斯托克斯和他们的数学方程

2023-12-19 18:10

本文主要是介绍华院计算 | 纳维、斯托克斯和他们的数学方程,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

8bc314c6e3a2f00ef90b66556620014c.jpeg

6dc04b51f3c9049348616a00b76eb241.jpeg

【一】

1998年,美国投资咨询与矿业金融大咖兰登·克莱(Landon T. Clay,1926–2017)捐助成立了克莱数学研究所。这间非盈利私立研究所在2000年5月24日宣布设置千禧年数学大奖(The Millennium Grand Challenge in Mathematics),给能够解答七大数学难题的每人每项一百万美元的奖金。其中,第四道难题是要证明流体动力学中最重要的“纳维–斯托克斯方程”(Navier-Stokes equation)是否有解且为唯一。这则新闻把一个艰深数学问题及其称谓广为普及,甚至传播到数学领域之外的大庭广众。

纳维-斯托克斯方程是描述粘性不可压缩流体动量守恒的运动方程组。该方程在不同条件下具有不同的形式。它在三维空间的一般向量表达式包括简单的线性质量守恒方程∇⋅𝑣=0和复杂的非线性动量守恒方程

e43e2a16485d3dbd7eb60798b3946009.jpeg

纳维-斯托克斯方程的发展经历了一个漫长的过程:

1775年,欧拉(Leonhard Euler,1707-1783)在论文“流体运动的一般原理”(Principes généraux du mouvement des fluides)中根据理想无粘性不可压缩流体所受外力时的动量变化推导出了一组微分方程式。

1821年,纳维(Claude-Louis Navier,1785-1836)推广了欧拉方程,考虑到流体内部分子之间相互作用力所产生的粘性,特别是在靠边界层的粘性对流体的显著影响,建立了一个不可压缩流体动力平衡和运动基本方程式。

1827年,柯西(Augustin-Louis Cauchy,1789-1857)在欧拉方程中引入流体微团应力张量的概念,推导出一个“柯西动量方程”,准确地描述了流体运动不同于固体运动的基本规律。

1831年,泊松(Siméon D. Poisson,1781-1840)建立了一个可压缩流体的运动方程。

1845年,乔治·斯托克斯(George G. Stokes,1819-1903)在纳维方程的基础上推导出了带常数粘性的不可压缩流体的动量守恒运动方程式,也就是现在著名的“纳维-斯托克斯方程”。

在这近百年的接力研究过程中,不少数学家和物理学家都在多个相关问题上作出过不同形式和不同程度的贡献,包括在非粘性流体数学分析方面的约翰·伯努利(Johann Bernoulli,1667-1748)、达朗贝尔(Jean le Rond d’Alembert,1717-1783)、拉格朗日(Joseph-Louis Lagrange,1736-1813)、拉普拉斯(Pierre-Simon Laplace,1749-1827)、高斯(J. Carl F. Guass,1777-1855)、圣维南(A. Jean Claude B. de Saint-Venant,1797-1886),以及带粘性流体和湍流方面的雷诺(Osborne Reynolds,1842-1912)和泰勒(Geoffrey I. Taylor,1886-1975),等等。

今天,纳维-斯托克斯方程解的存在唯一性和许多特征(例如光滑性)尚待证明。1934年,法国数学家让-勒雷(Jean Leray,1906-1998)证明了方程在某种意义下弱解的存在。目前已知方程有一百多个特解,但一般解只能借助于计算机求数值近似。1999年,美国数学家史蒂夫·斯梅尔(Steve Smale,1930-)在题为“下一个世纪的数学问题”(Mathematical problems for the next century)的文章中也把纳维-斯托克斯方程列为现代数学尚待解决的18个极大难题之一。

【二】

图1 克劳德-路易·纳维塑像
图1 克劳德-路易·纳维塑像

克劳德-路易·纳维(Claude-Louis Navier,1785年2月10日-1836年8月21日)是法国工程师与物理学家,他的主要贡献在力学领域。纳维的父亲克劳德-伯纳德(Claude-Bernard Navier,1756-1793)是个律师,在法国大革命时期是巴黎国民议会(Assemblée Nationale de Paris)一位议员。父亲在1793年去世后,母亲珍妮-玛丽(Jeanne-Marie Pourcher)自己回归老家Chalon-sur-Saône,把纳维留在巴黎由叔公、数学家和土木工程师高西(Émiland Gauthey,1732-1806)负责教育培养。当年,高西在桥梁公路学校(École des Ponts et Chaussées)兼职任教,并于1791年被政府任命为桥梁与公路建筑集团的总督导。叔公是少年纳维心目中的模范榜样,让他立志将来当一个土木工程师。

1802年,17岁的纳维进入巴黎综合理工学院(École Polytechnique)就读。纳维的入学成绩几乎排在最低端,但他后来十分努力,一年后就进入了前十名优等生行列。在理工学院,纳维修读了数学家傅立叶(Jean-Baptiste J. Fourier,1768-1830)主讲的数学分析课程。纳维学习用功、思维敏捷,深得傅立叶教授的喜爱。后来,两人成为终生挚友。1804年,纳维以优等生资格从理工学院毕业。1806年,叔公高西去世,建筑集团邀请纳维去整理他留下的丰富文献记录和研究资料。在接下来的十多年里,纳维编辑了他叔公的手稿,主要是一些传统经验方法在土木工程中的多种应用,其中也加上了自己理论力学研究的一些分析结果。1813年,纳维为叔公出了新版的《桥梁特性》(Traité des ponts)以及为法国著名水力学工程师伯纳德·德·贝利多尔(Bernard-Forest de Bélidor,1698-1761)出了新版的《工程手册》。

1819年,纳维应聘在桥梁公路学校讲授应用力学课程。在那里,纳维对传统的力学课程教学大纲做了一次彻底的改革,增加了大量的物理学和数学分析内容。1820年代初期,纳维与在巴黎综合理工学院任教的数学家柯西和加布里埃尔·拉梅(Gabriel Lamé,1795-1870)一起建立了力学中的弹性基本理论。

1821年,纳维推导出一条新的弹性理论数学公式,第一次获得了具有足够精确度的弹性力计算结果。1823年,纳维受政府委托前往英格兰和苏格兰,去了解链式悬索桥的建设情况。事后他发表了一份著名的调查报告,其中提出了桥梁史上第一个悬索桥理论。这一成绩让他在1824年获选为法国科学院院士。不过,纳维的科学研究和工程实践也并非不无挫折。1824年,纳维在主持荣军院吊桥(Pont des Invalides)设计时,在计算中没有留出足够的安全余量,导致桥梁开裂并最后被拆除。当时指控不断,他因而受到了政府委员会的公开问责,认为他的设计过分依赖于物理数学而不是工程经验。

3452e2a9f45098deb0950037e670756f.jpeg
图2 纳维的《桥梁与道路学院课程总结》

1826年,纳维出版了《桥梁与道路学院课程总结》(Résumé des Leçons Donnees a l'École des Ponts et Chaussées),文中论证了弹性模量是材料的一个基本属性,只和物体的材质有关而与其结构几何形状无关。这项重要成果让他成为结构分析理论的奠基人之一。1828年,纳维和泊松在《化学与物理年鉴》(Annales de Chimie et de Physique)杂志上发生了关于弹性理论的争论。不过,因为两者都基于相同的分子假说,这场辩论无果而终。

1830年7月,法国爆发了“七月革命”。国王查理十世(Charles Philippe X,1757-1836)逃离法国,由波旁王朝(Bourbon Restoration)国王菲利普一世(Louis Philippe I,1773-1850)继任。这事件让爱国者数学家柯西大为震惊。柯西留下了家人,独自离开巴黎到了瑞士,在弗里堡(Fribourg)待了一段短暂时间。在那里,他作了一个重大决定,拒绝宣誓效忠新政权。他因此失去了在巴黎的所有职位,除了不需要宣誓的科学院院士资格。

1831 年,柯西前往意大利都灵(Turin),随后接受了撒丁岛(Kingdom of Sardinia)国王的邀请,担任专门为他设立的理论物理学教授职位。1831年,纳维接替了柯西在巴黎综合理工学院的职位,担任微积分及力学教授。

1836年8月21日,纳维在巴黎去世,享年51岁。他的名字和其他71 位法国名人一起,被镌刻在巴黎埃菲尔铁塔上。

【三】

b1ae54e01a013b33676172bf88780943.jpeg
图3 乔治·斯托克斯肖像(Alice Power绘,英国皇家学会存展)


乔治·斯托克斯(George Gabriel Stoke)于1819年8月13日出生在爱尔兰北部Sligo郡Skreen教区,父亲加布里埃尔(Gabriel Stokes,1762-1834)是教区长,母亲伊丽莎白(Elizabeth Haughton)来自另一位教区长家庭。斯托克斯在家中八个孩子里最小,三个哥哥长大后都成为了牧师。斯托克斯从小接受家教,到1832年13岁时才被送到都柏林(Dublin)一所学校去接受正规教育。在学校里,斯托克斯显示出不凡的几何学观察和推证能力。三年后,斯托克斯到了英国Bristol College学习,该学院院长是位数学家,是斯托克斯哥哥在剑桥大学的同学。院长对斯托克斯的数学天分非常惊讶,力劝他到剑桥深造。

1837年,斯托克斯获得奖学金进入了剑桥大学Pembroke College就读。从此之后,他终其一生从未离开过该学院。在剑桥,斯托克斯得到了数学家威廉·霍普金斯(William Hopkins,1793-1866)的喜爱和特别指导。

1841年,斯托克斯以数学第一名荣誉生(Senior Wrangler)毕业并获Smith奖。学院马上提供了他一个研究员职位(Master of Pembroke College)。在霍普金斯建议和指导下,斯托克斯开始流体动力学方面的研究。斯托克斯在1842年和1843年相继发表了两篇论文“不可压缩流体的稳态运动”(On the steady motion of incompressible fluids)和“关于流体运动的一些研究”(On some cases of fluid motion)。之后,他进一步考虑有黏性流体的运动,并于1845年发表了论文“论有内部摩擦力流体之运动”(On the theories of the internal friction of fluids in motion)。在发表此论文之前,斯托克斯发现纳维早已在1822年写下相同的方程式,不过其结论带有推测性,而斯托克斯则是基于不同的假设条件下用严格的数学理论推导出该方程式。虽然纳维在先,斯托克斯再三考虑之后还是发表了该论文。此举厥功至伟,文中严格建立起来的方程式就是千禧难题“纳维-斯托克斯方程”。

1846年,斯托克斯在英国科学促进会(British Association for the Advancement of Science)作了题为“流体动力学最新研究报告”(Report on recent researches in hydrodynamics)的著名演讲,一时好评如潮,之后各项荣誉接踵而至。

1847年,斯托克斯在数学复变函数积分和无穷级数分析中有一个有趣的发现:他描述了两个独立函数的线性组合曲线(“斯托克斯曲线”)在渐近变化过程中当变量的模越过某个阈值(“斯托克斯界线”)时会出现不连续性,后人称之为“斯托克斯现象”。1847年,斯托克斯参与了对英国Dee Bridge断塌灾难的调查。他被委任为皇家委员会(Royal Commission)成员,对移动机车在大桥上的压力作了详尽计算和分析,确定了主要原因是火车重量超过了桥梁铸铁的承载能力。事后,他还被任命为负责研究风压对建筑结构影响的皇家委员会成员。

1849年,他在《剑桥哲学学会汇刊》发表了一篇题为“衍射动力学理论”(The dynamical theory of diffraction)的长文,证明了光的偏振面是垂直于传播方向的。同年,他还发表了一篇关于大地测量的重要论文“关于地球表面重力的变化”(On the variation of gravity at the surface of the earth),指出地表重力的测量并不依赖于地球内部的性质。同年,他接受了位于伦敦的政府矿业学院物理学教授一职。

1849年,三十而立的斯托克斯获剑桥卢卡斯数学教授(Lucasian Chair of Mathematics)席位,1851年被选为英国皇家学会院士,1854年出任皇家学会秘书,1887-1892年出任皇家学会院长。他是历史上除了牛顿以外唯一曾获英国皇家学会秘书、院长以及卢卡斯数学讲座教授三个头衔的人。1850年,斯托克斯开始尝试考察过滤光(filtered light)。他让阳光穿过蓝色玻璃,然后将光束照射到黄色醌溶液中,以此产生强烈的黄色照明。斯托克斯对不同化合物的溶液进行了相同的实验,但发现只有一些化合物显示出与原始光束不同的颜色。

1851年,斯托克斯对德裔英国天文学家威廉·赫歇尔(William Herschel,1738-1822)在1845年报告的关于奎宁硫酸盐无色透明溶液表面在光照下会产生蓝色闪光的现象产生了极大兴趣。他亲自做了实验,发现蓝色闪光是由奎宁溶液吸收了不可见的紫外线而引起的。他将这种光命名为“荧光”(fluorescence)。后来他还发现,石英棱镜而不是玻璃棱镜可用于研究光谱的紫外线部分,从而确定了玻璃会吸收紫外线。他认为这些材料具有将不可见的紫外线辐射转化为波长较长的可见光的能力。这一材料特性被称为“斯托克斯平移”(Stokes shift),是拉曼散射(Raman scattering)的基础。研究中,他还引进了描述光偏振(polarization)状态的“斯托克斯参数”。翌年,斯托克斯把这项极其重要的研究结果以“论光折射系数的变化”(On the change of the refrangibility of light)为题发表在《英国皇家学会哲学汇刊》(Philosophical Transactions of the Royal Society)。同年,他因此项成果荣获英国皇家学会的Rumford奖章。

1851年,斯托克斯在粘性流体运动的研究工作中计算了球体在粘性介质中下落的终端速度,得到了著名的“斯托克斯定律”(Stokes law)。后来,厘米-克-秒系统中的粘度单位被命名为“斯托克斯”。当年,斯托克斯还探讨了流体内部摩擦力对单摆运动的影响。他的研究成果是流体动力学的一个新里程碑,对一些复杂的自然现象如空中云的飘浮和水中浪的波动等作出了合理的力学解释,并对涉及水在河沟和管道中的流动以及船身在水中的表面阻力等应用设计提供了技术参考。他也曾对声学做过研究,包括风对声音强度影响的分析。

1857年,38岁功成名就的斯托克斯和剑桥Armagh天文台一位天文学家汤姆士·罗宾逊(Thomas R. Robinson,1792-1882)的女儿玛丽·苏珊娜(Mary Susanna Robinson,1825-1899)结了婚。新婚夫妇搬进了Lensfield Cottage,在那里斯托克斯还建立了一个小小实验室。夫妇俩养育有五个子女并在那里度过了余生。婚后,斯托克斯在学校和政府部门担任了较多的行政职务,并且倾向于更为实验性而不是理论性的研究。

1862年,斯托克斯在英国科学促进会作了一个题为“论电光的长光谱”(On the long spectrum of electric light)的著名演讲,报告了从一串圆盘双重反射或透射的光强度变化的实验分析结果。该报告随后发表在《英国皇家学会哲学汇刊》。斯托克斯在化学领域也有过一些研究和贡献。他和英国科学家、英国科学促进会创始人威廉·哈考特(William V. Harcourt,1789-1871)合作,研究了各种玻璃的化学成分与光学性能之间的关系和透明度。他还和德国物理及化学家菲利克斯·霍普-塞勒(E. Felix I. Hoppe-Seyler,1825-1895)合作,发现了血红蛋白具有氧运输功能,并显示了血红蛋白溶液通过气体时所产生的颜色变化。

1869年,斯托克斯被选为英国科学促进协会主席。1877年,斯托克斯荣获德国哥廷根大学颁发高斯奖章。1880年代初期,斯托克斯和英国物理学家威廉·克鲁克斯(William Crookes,1832-1919)合作研究射线。当时他们认为阴极射线(cathode ray)是带负电的粒子,但德国物理学家们则不以为然。1897年3月31日,斯托克斯写信给克鲁克斯,指出德国人所说的阴极射线Kathodenstrahlen“根本就不是射线,而是分子流”。一个月后,英国物理学家约瑟夫·汤姆森(Joseph J. Thomson,1866-1940)宣布从实验中发现了电子(electron)并因之而获得了 1906年的物理学诺贝尔奖。1889年,斯托克斯被英国皇家册封为男爵(Sir)。

1893年,斯托克斯荣获声誉崇高的英国皇家学会科普里奖章(Copley Medal)。物理学家开尔文勋爵(Lord Kelvin,1824-1907)在颁奖词中说:“五十二年前,他就开始了流体运动的主题研究,然后精巧地勾画出流体特性的完全真实本性,并以新颖的视角通过深入到物质的构成中去发现数学问题的极致之美”。记录表明,斯托克斯和开尔文勋爵有长达55年的友谊和合作,期间双方互助共进,各自成果累累。其中一件轶事说,1854年两人在一次讨论中斯托克斯推测可以用太阳光光谱中的Fraunhofer线来阐明太阳的化学成分。但由于他没有正式发表文章,后来1859年德国物理学家古斯塔夫·基尔霍夫(Gustav R. Kirchhoff,1824-1887)和罗伯特·本生(Robert W. Bunsen,1811-1899)发表频谱分析原理时,斯托克斯并不作声,毫无意念去争夺这项成果。
1899年,斯托克斯因担任卢卡斯讲座教授50周年接受了剑桥大学的喜庆,成为“禧年教授”(Professorial Jubilee)。

1901年,斯托克斯获普鲁士科学院颁发Helmholtz 奖章。1902年,斯托克斯获剑桥大学Pembroke College授予最高荣誉。同年,他获挪威 Royal Frederick University(现为University of Oslo)授予荣誉数学博士学位。之前,斯托克斯还获得过英国剑桥、牛津、爱丁堡、格拉斯哥、阿伯丁和爱尔兰都柏林大学的荣誉博士学位。1903年2月1日,斯托克斯在剑桥Lensfield Cottage去世,享年84岁。他被葬于剑桥的Mill Road坟场。为了纪念斯托克斯,后人在他的出生地矗立了一块牌匾。

0599291495beab1cad6ce2cee1be68be.jpeg
图4 斯托克斯出生地纪念牌匾

斯托克斯为人品性谦和,生活简朴,是个虔诚教徒。他毕生对科学与信仰之间的联系非常感兴趣,并积极参与创造论与进化论的各种争论。
1885年,他开始担任维多利亚研究所(Victoria Institute)所长职务直至去世。该研究所探索宗教教义与科学发现之间的关系。

斯托克斯的全部数学和物理学论文经整理后分五卷出版。前三卷是斯托克斯于1880年、1883年和1891年亲自编辑的,后两卷则是在斯托克斯去世后由爱尔兰裔英国物理学和数学家约瑟夫·拉莫尔(Joseph Larmor,1857-1942)于1905年完成的。

c9c0c7e7081da69c454ea5890e5a52d2.jpeg
图5 斯托克斯数学和物理学论文集

回顾十九世纪中叶,剑桥是世界学术界的数学物理学中心,成绩巨大、声誉斐然。斯托克斯和詹姆斯·麦克斯韦(James Clerk Maxwell, 1831-1879)以及开尔文勋爵一起,被誉为是对剑桥数学物理中心做出了最重要学术贡献的“自然哲学家三人组”(The Trio of Natural Philosophers),流芳百世。

b09df1d06030e863f285081c8de6f311.jpeg
图6 乔治·斯托克斯纪念像标题

这篇关于华院计算 | 纳维、斯托克斯和他们的数学方程的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/513195

相关文章

uva 10014 Simple calculations(数学推导)

直接按照题意来推导最后的结果就行了。 开始的时候只做到了第一个推导,第二次没有继续下去。 代码: #include<stdio.h>int main(){int T, n, i;double a, aa, sum, temp, ans;scanf("%d", &T);while(T--){scanf("%d", &n);scanf("%lf", &first);scanf

uva 10025 The ? 1 ? 2 ? ... ? n = k problem(数学)

题意是    ?  1  ?  2  ?  ...  ?  n = k 式子中给k,? 处可以填 + 也可以填 - ,问最小满足条件的n。 e.g k = 12  - 1 + 2 + 3 + 4 + 5 + 6 - 7 = 12 with n = 7。 先给证明,令 S(n) = 1 + 2 + 3 + 4 + 5 + .... + n 暴搜n,搜出当 S(n) >=

uva 11044 Searching for Nessy(小学数学)

题意是给出一个n*m的格子,求出里面有多少个不重合的九宫格。 (rows / 3) * (columns / 3) K.o 代码: #include <stdio.h>int main(){int ncase;scanf("%d", &ncase);while (ncase--){int rows, columns;scanf("%d%d", &rows, &col

【生成模型系列(初级)】嵌入(Embedding)方程——自然语言处理的数学灵魂【通俗理解】

【通俗理解】嵌入(Embedding)方程——自然语言处理的数学灵魂 关键词提炼 #嵌入方程 #自然语言处理 #词向量 #机器学习 #神经网络 #向量空间模型 #Siri #Google翻译 #AlexNet 第一节:嵌入方程的类比与核心概念【尽可能通俗】 嵌入方程可以被看作是自然语言处理中的“翻译机”,它将文本中的单词或短语转换成计算机能够理解的数学形式,即向量。 正如翻译机将一种语言

poj 1113 凸包+简单几何计算

题意: 给N个平面上的点,现在要在离点外L米处建城墙,使得城墙把所有点都包含进去且城墙的长度最短。 解析: 韬哥出的某次训练赛上A出的第一道计算几何,算是大水题吧。 用convexhull算法把凸包求出来,然后加加减减就A了。 计算见下图: 好久没玩画图了啊好开心。 代码: #include <iostream>#include <cstdio>#inclu

uva 1342 欧拉定理(计算几何模板)

题意: 给几个点,把这几个点用直线连起来,求这些直线把平面分成了几个。 解析: 欧拉定理: 顶点数 + 面数 - 边数= 2。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstring>#include <cmath>#inc

uva 11178 计算集合模板题

题意: 求三角形行三个角三等分点射线交出的内三角形坐标。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstring>#include <cmath>#include <stack>#include <vector>#include <

XTU 1237 计算几何

题面: Magic Triangle Problem Description: Huangriq is a respectful acmer in ACM team of XTU because he brought the best place in regional contest in history of XTU. Huangriq works in a big compa

数学建模笔记—— 非线性规划

数学建模笔记—— 非线性规划 非线性规划1. 模型原理1.1 非线性规划的标准型1.2 非线性规划求解的Matlab函数 2. 典型例题3. matlab代码求解3.1 例1 一个简单示例3.2 例2 选址问题1. 第一问 线性规划2. 第二问 非线性规划 非线性规划 非线性规划是一种求解目标函数或约束条件中有一个或几个非线性函数的最优化问题的方法。运筹学的一个重要分支。2

音视频入门基础:WAV专题(10)——FFmpeg源码中计算WAV音频文件每个packet的pts、dts的实现

一、引言 从文章《音视频入门基础:WAV专题(6)——通过FFprobe显示WAV音频文件每个数据包的信息》中我们可以知道,通过FFprobe命令可以打印WAV音频文件每个packet(也称为数据包或多媒体包)的信息,这些信息包含该packet的pts、dts: 打印出来的“pts”实际是AVPacket结构体中的成员变量pts,是以AVStream->time_base为单位的显