首页/文章/ 详情

纳维点亮N-S方程之光,照亮CFD前行之路

3小时前浏览7


导读:春节在家,早上醒来,突然想写《一个工程师的救赎》。这个故事,关于纳维在流体力学领域的执着坚守,还有那份不被理解的孤独与坚持。此刻我满心期待将这个故事讲述给大家,让更多人了解科学背后的动人历程。
值此新春佳节,我衷心祝愿每一位读者朋友新年快乐,万事顺遂!2025 年-期待我精心撰稿的计算流体力学书籍能早日与大家见面(同步课程会在仿真秀官网发布),它就像久别重逢的老友,为大家揭开流体力学的神秘面纱。

一、纳维就像被黏性困住的普罗米修斯

1826年的深秋,巴黎塞纳河畔的一座大桥上挤满了愤怒的市民。因为这座大桥是法国科学院的的天才院士纳维设计的,但桥基在涡旋的撕扯下不断剥落碎石,这座桥建成还不到三年眼瞅着就要塌了!人群愤怒又嘲笑工程师纳维的无能。法兰西科学院公开谴责纳维“用数学谋杀了一座桥”,人们画漫画讽刺纳维:“他举着写满方程的盾牌,被塞纳河女神用浪涛击碎” 。工程师纳维因桥梁事故遭学界排挤,人们不再信任他。
图1.塞纳河上的桥(已不是纳维建造的了)
纳维是严格按照当时的流体力学知识进行设计的,还加了安全阈值,按理说大桥应该稳如磐石。他实在想不通问题出在哪里,他决定把问题搞清楚——“哪怕死,也要死的明明白白”。
工程师纳维开始追根溯源,当时的流体力学课本教导学生流体符合欧拉方程,但现实证明严格按照欧拉方程设计的桥梁不堪一击,自己却背负了骂名。

图2 流体力学中的欧拉方程
从工程的角度出发,他意识到流体力学中的欧拉方程是有问题的,他开始进入长时间沉浸式的思考。经过长达四年的思考,工程师纳维认识到“流体内存在一种隐秘的摩擦,它不屈服于欧拉的傲慢……。”他把这种摩擦定义为流体的黏性,他开始尝试修正欧拉方程!
1830年法国七月革命期间,暴民烧毁了纳维的实验室,他蜷缩在圣母院钟楼里,用炭笔在《论流体运动定律》手稿边缘写下:“每一座倒塌的桥梁,都是流体写给人类的情书——只是我们总误读它的笔迹。黏性不是流体的缺陷,而是它对爱人的执着。” 纳维从工程的角度出发在欧拉方程的右边加上了黏性项,从此后世每一个学习流体力学的学子都记住了工程师纳维的名字。然而这位伟大的工程师生前并没有被世人原谅,他背负着骂名,遭受着排挤,郁郁而终。
1845年,剑桥大学的斯托克斯在图书馆发现纳维残破的手稿,他通过严谨的连续介质理论修正分子假设,最终完成N-S方程,从此这个千年公式的完整形式问世,他们两人的名字也因这个公式而流芳千古。
图3 完整形式的N-S方程,在欧拉方程的基础上加了粘性项
在致皇家学会的信中,斯托克斯写道:“纳维就像被黏性困住的普罗米修斯,而我不过拾起了他掉落的火种。”
法国人后知后觉,开始原谅与感念纳维。他们在法国的巴黎铁塔上镌刻上纳维的名字,以纪念与感怀这位为人类工程与科学做成伟大贡献的工程师。

图4 法国埃菲尔铁塔上镌刻有纳维的名字

图5 纳维肖像

二、我所理解的CFD的基本求解思想

如今,纳维成为了科学史上的一座丰碑,激励着无数科研人员和工程师。在现代工程领域,从航空航天到汽车制造,从船舶设计到建筑施工,N - S 方程无处不在,它就像一把万能 钥 匙,为工程师们打开了理解流体世界的大门。在航空发动机的设计中,工程师们运用 N - S 方程精确计算气流在发动机内部的流动,从而优化发动机的性能,提高燃油效率,降低污染物排放。在汽车风洞试验中,借助 N - S 方程模拟汽车周围的气流,减少风阻,提升汽车的稳定性和速度。而这些成就的背后,都离不开纳维当初的坚持与探索。

提到对流体世界的探索,就不得不说到 CFD(计算流体动力学)。CFD 的基本求解思想,是将描述流体运动的 N - S 方程进行离散化处理,转化为计算机能够处理的代数方程组。通过数值计算的方法,求解这些方程组,从而得到流场中各个物理量的分布,如速度、压力、温度等。那为什么我们称N-S方程为非线性偏微分方程?

1、线性和非线性

首先说明什么是线性系统。比如两个同学称体重,两个人的总重量是他们各自体重之和,他们的体重不会相互影响,即整体等于部分之和,这是线性系统的第一个关键特性;线性系统的第二个特性是,原因与结果成正比,想象一下你用脚踢球的情形,我们忽略空气阻力,如果你用F的力来踢球,球运动的直线距离是x,如果你用2F的力来踢球,球运动的直线距离应该是2x。即原因和结果成正比。满足这两个特性(整体等于部分之和,原因和结果成正比)的系统就是线性系统。

然而自然界的很多事情都比踢球复杂得多。当系统的各个部分相互干扰、合作或竞争时,就会发生非线性的相互作用[11]。大部分日常活动都是非线性的,比如你同时听两首最喜欢的歌,你不会得到双倍的快乐。你左眼看英语,右眼看高数,你也不会得到双倍的知识。

非线性让世界变得丰富多彩、美妙而复杂,还常常不可预测,而CFD恰恰是研究非线性。

2、什么是(非)线性微分方程呢?

简单得说,线性微分方程是指未知函数及其各阶导数都是一次幂,且未知函数及其各阶导数互不掺混的微分方程,否则称其为非线性微分方程。线性微分方程需要满足以下性质:

(1) 只能出现自变量、函数本身,以及函数的任何阶次的导函数;

(2)函数本身跟所有的导函数之间除了加减之外,不可以有任何其他运算(比如乘除、平方、开根号等);

(3)函数本身跟函数本身,各阶导数本身跟各阶导数本身,都不可以有加减之外的运算;

(4)不允许对函数本身、各阶导数做任何形式的复合运算,例如siny、cosy、tany、lny、y^2、y^3、y^x、exp^y。

若一个微分方程不满足上述要求,则是非线性微分方程。

我们举例说明,例如

因为方程中 y、y、y"、y"的最高次幂都是 1,且没有出现y、y、y'、y"(未知函数及其各阶导数)之间的相互掺混。虽然函数y和导数y、y"前面出现了 x、sin(x)、In(x)等自变量的初等函数,但都是允许的,所以以上方程都是线性微分方程。

反之: 

因为方程中的最高次幂要么不全是 1,要么出现了(未知函数y及其各阶导数)的相互掺混,所以以上方程都是非线性微分方程。

我们经常说N-S方程是非线性偏微分方程,学完上面的知识,我们自己来判定一下。

N-S方程

方程中出现了自变量x、y、z、t,自变量的个数有4个,所以是偏微分方程。同时方程中出现了函数u与其对应导数的乘积,比如函数u和其一阶导数乘在了一起,所以是非线性。综合起来,就叫作非线性偏微分方程。

三、我的2025新年愿望

去年,我在仿真秀官网和APP原创的视频课程传热学理论课程86讲》、《流体力学考研》、《工程热力学考研》和《计算流体力学湍流模型大串讲》,深受用户认可和喜爱。其中《工程热力学考研》课程播放量达到230000+次。因此,我荣获仿真秀年度优秀讲师称号。


面向本科生群体,我主要讲授传热学、工程热力学以及流体力学这三门基础且重要的课程,为本科生构建扎实的专业知识框架。而对于研究生和博士生这些更高层次的学习者,我则聚焦于计算流体力学、湍流模型大串讲以及 OpenFOAM 初级入门教程等进阶课程,助力他们深入探索专业前沿知识,提升科研能力。

田东Joshua_Fluent_流体基础_暖通_仿真优质内容创作者-仿真秀

2025 年,期望能与仿真秀共同打造更多高质量的学术与教学成果我还计划将 CFD 理论方面的知识整理出版成书,并将其录制成配套课程,发布在仿真秀官网和APP,以便更广泛地传播专业知识,帮助更多学习者深入理解 CFD 理论。其次,我打算深入学习模型预测控制,不断提升自己在该领域的专业能力,为后续的科研与教学工作注入新的活力。

(完)


来源:仿真秀App
OpenFOAM非线性湍流航空航天船舶汽车建筑理论控制试验
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2025-02-06
最近编辑:3小时前
仿真圈
技术圈粉 知识付费 学习强国
获赞 10420粉丝 21869文章 3640课程 230
点赞
收藏
未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈