1. 从“控制体”视角看流体N-S方程不再是天书我学流体力学的时候也干过一件蠢事把N-S方程当英语单词先抄三遍再背三遍结果考完就忘做题照样卡壳。后来才想明白问题不在于记性而在于我从一开始就选错了理解方式。N-S方程不是需要背诵的结论它完全可以从一个更朴素、更直观的视角推出来这个视角就是“控制体”。控制体是什么你可以把它理解成在流场里画一个虚构的盒子盒子可以很小也可以包围整个机翼、整个管道它不随流体运动固定在你关心的空间里。流体从盒子流进流出我们在盒子表面上记录进出的量再结合盒子内部的积累就能写出各种守恒定律。工程上绝大多数计算管道、泵、风机、飞行器外流场都是站在这种固定盒子的视角去看问题。而流体力学教材里那些让人头大的方程包括N-S方程本质上都是从这个视角出发一步步用守恒定律推出来的。这篇文章适合两类人第一类是刚开始学流体力学、正被N-S方程折磨的学生第二类是工作后想回头把底层逻辑补上的工程师。全文不要求你先背任何公式我带着你从零开始先推导输运方程再用输运方程推出连续性方程和N-S方程。推完之后你会发现N-S方程只是“动量守恒定律在控制体上的展开”想忘记都难。1.1 两种观察方式追着流体跑还是守着空间看流体力学里一直有两种观察方式很多人一开始没分清楚后面全乱套。第一种叫系统视角。系统指的是“一堆固定的流体微团”你跟着这堆流体一起走它变形、移动、旋转你都跟着。就像你骑着一匹马去追踪一条河里的同一个水团这个水团从上游漂到下游你始终盯着它。物理定律天然是在这种视角下成立的牛顿第二定律说的是“这个流体团的动量变化率等于所受合力”而不是“这个固定盒子里的动量变化率”。第二种叫控制体视角。控制体是你预先选定的一块空间区域边界叫控制面。流体穿过控制面进进出出你站在控制体外只关心这块区域内部和边界上的变化。这就像你在河上的一座桥墩旁边圈出一段水面不跟着水跑只记录单位时间内有多少水从上游进来、多少水从下游出去、桥墩这一段里水位涨了还是落了。工程上哪种更实用控制体。因为你要设计的管道、阀门、风洞、机翼都是固定的几何边界你关心的是“流过这个截面有多少流量”“作用在这个物体上的力有多大”。你总不能每算一个力就派个人骑着马追着同一个水团跑遍全流场吧。所以实际计算通常用控制体但物理定律本质上是给系统用的这就产生了一个矛盾定律说系统的事计算要用控制体两边怎么对接对接工具就是输运方程它回答一个问题系统内某个物理量随时间的变化率如何用控制体上的当地导数和对流通量来表示。N-S方程就是从这个对接过程中一步步长出来的。1.2 输运方程到底解决了什么问题先举一个最简单的问题一个浴缸放水水位是涨还是落浴缸就是控制体水龙头流进来的水、排水口流出去的水、浴缸里水量的变化三者关系一眼就能写出来水位变化 流进来的水量 − 流出去的水量。输运方程干的也是这件事只不过把“水量”换成了质量、动量、能量这类一般化的物理量。你可能会觉得这有什么好推的但流体不是刚体它无时无刻不在变形、旋转、穿越边界。同一时刻控制体内可能有流体在加速、有流体在流出、有流体在流入物理量在时间的推进中变化得非常复杂。如果直接写出系统导数和控制体积分之间的关系很容易漏项、错项。输运方程的价值就是给这个关系一个严格的、通用的形式。只要你确定你关心的是哪个量把它代进去各种守恒方程都能出来。所以我的建议很简单别背N-S方程先背输运方程再学会代公式。后面你会看到N-S方程的每一项都可以在输运方程里找到对应位置而不是凭空冒出来的。2. 手把手推导输运方程加了这么多铺垫现在就开始正题。我们要推导的是一个关于任意物理量B的方程。先定义清楚B是一个广延量比如质量、动量、能量这样跟物质的量成正比的总量b是单位质量对应的B专业点叫强度量。比如B是质量b就等于1B是动量b就是速度VB是动能b就是v²/2。B和b的关系通过密度联系B ∫ ρb dV这个积分可以放在系统上也可以放在控制体上区别只在于积分区域不同。我们现在先写成系统的形式B_sys(t) ∫_{sys(t)} ρb dV其中sys(t)表示t时刻系统占据的体积。系统的特点是它永远由同一堆流体微团组成随流场运动但形状和位置都会变。我们要算的是系统量随时间的变化率即 d B_sys / d t。这就是牛顿第二定律左边那种导数。输运方程的目标是把这种“跟随流体”的导数改写成“站在固定控制体上看”的表达式。2.1 最直观的做法对比两个相邻时刻推导的关键在于选t时刻让系统刚好和控制体重合。这只是一个数学上的约定因为物理定律里的系统是任选的此刻它与控制体重合方便我们做转换。t时刻系统占据的体积就是控制体V_CV。但t Δt时刻系统里的流体微团已经跟着流场往前走了它们会有一部分流出原来的控制体同时控制体外面的另一部分流体并没有进入系统。注意系统是那堆固定的流体微团不是控制体内的所有流体。所以t Δt时刻系统的体积可以分成三块仍然留在控制体内的部分从控制面流出去的部分原先在控制体内部、现在依然在那里的那部分里的流体依然属于系统。等一下这里要小心。系统和控制体在t时刻重合但在t Δt时刻系统可能有一部分已经流出控制体边界所以系统体积会偏离控制体。我们推导时考虑的是系统量的变化由“控制体内当地变化”加上“穿越控制面的净流出”共同组成。用数学表达系统量的变化率近似写成d B_sys / dt ≈ [B_sys(t Δt) − B_sys(t)] / Δt其中B_sys(t) ∫_{V_CV} ρb dV因为t时刻系统与控制体重合。B_sys(t Δt)则分成两部分仍然在控制体内的量加上流出控制面的量。经过整理并取Δt → 0的极限就得到输运方程的标准形式d/dt ∫_{sys} ρb dV ∫_{V_CV} ∂(ρb)/∂t dV ∮_{A_CV} ρb (V·n) dA这里V是流体在控制面当地的速度n是控制面外法线方向单位向量V·n表示流体沿外法线方向穿越控制面的速度分量。这个式子就是“系统量的变化率 控制体内当地积累率 控制面上净对流输出率”。为了让你彻底记住我用自己的话翻译一遍系统里某物理量的总量变化等于你站在控制体里看到的当地变化再加上从边界上像传送带一样把物理量送出去的量。传出去的是正的贡献传进来的是负的贡献因为外法线方向是朝外的。就这么简单。2.2 通量项的正负号为什么这么重要很多人在这一步栽跟头问题都出在V·n的符号上。我多说几句。控制面的外法线n是垂直于表面、指向控制体外侧的向量。流体速度V如果顺着n的方向流出V·n就是正的代表物理量正在离开控制体对控制体来说是“损失”所以在方程右侧加号里它贡献为正。流体如果从外部流入V·n是负的代表物理量进入控制体对控制体来说是“获得”所以该项是负贡献。净流出的通量项取整个闭合面积分流出的正贡献和流入的负贡献会自然抵消一部分剩下的就是真正从控制体净流出去的量。我见过很多初学的人因为嫌弃V·n这个点积太数学就自己改成“流出减流入”结果在非平面控制面上算错。建议你无论画什么样的控制体都坚持写ρb(V·n)dA用外法线方向自己判断正负而不是背什么“进为正出为负”的口诀。前者是按定义推后者是凭感觉猜猜一定会出错。再举一个生活化的例子。想象你的手机话费是一个保单你先是这个保单的“系统”每个月你的话费变化等于你充进去的钱减去你花掉的钱。控制体就是你的手机账户余额所在的那个“账面”充值是流入消费是流出。输运方程说的就是你这个“系统里的钱”其实可以理解成你名下某个号码的话费余额的变化等于你最近充值的钱加花掉的钱。只不过在流体里“充值”和“消费”每时每刻都在控制面上发生而且不同位置的速率还不一样。2.3 从系统导数到物质导数输运方程更严谨的数学形式还可以用“物质导数”的概念来理解。我们常看到D/Dt这个符号它就是物质导数意思是“跟随一个流体微团去观察某个量随时间的变化率”。在N-S方程里D/Dt ∂/∂t V·∇前面是当地导数后面是对流导数。如果你把输运方程里的系统导数也写成这种形式会发现它恰恰就是把物质导数、当地导数和对流项串在了一起。∂(ρb)/∂t是控制体内某一点上物理量的当地变化V·∇(ρb)则是流体把物理量从一点搬运到另一点带来的变化。这两者合起来才构成跟随流体时看到的完整变化。这也是为什么N-S方程里总会出现V·∇V这样的项它不是额外的复杂项而是系统在流动过程中“随波逐流”的加速度的一部分。3. 从输运方程推守恒方程直到N-S输运方程本身只是工具真正的重头戏是把它应用到物理定律上。严格来说N-S方程不是单独一个方程而是连续性方程加上动量方程在特定假设下的组合。下面我们就按这个顺序一步一步推。3.1 连续性方程取B质量质量守恒是流体力学最基本的定律它说系统的质量不会凭空产生或消失。也就是说d/dt ∫_{sys} ρ dV 0。对质量这个量b1。把b1代入输运方程得到∫_{V_CV} ∂ρ/∂t dV ∮_{A_CV} ρ(V·n) dA 0把面积分用高斯散度定理转换成体积分这里不展开讲高斯的证明你可以把它理解成“边界穿出的通量等于内部散度的总和”∫_{V_CV} [∂ρ/∂t ∇·(ρV)] dV 0由于控制体可以任意选取要让任意体积内的积分恒等于零只能是括号里的被积函数处处为零∂ρ/∂t ∇·(ρV) 0这就是连续性方程。每一项的物理意义第一项是当地密度随时间的变化率第二项是密度通量的散度说白了就是单位体积内质量“向外发散”的速率。两者相加为零意味着质量既不会多也不会少。对不可压缩流体密度沿路径不变并且通常取为常数。连续性方程退化成∇·V 0也就是速度场的散度为零。这时的直观含义是流体不能在某一点堆积你从这边流进来多少体积就得从那边挤出去多少体积就像地铁车厢门关闭瞬间人挤人但总量不变。这个式子看着很简单却是后面N-S方程简化的重要基础。3.2 动量方程取B动量现在换一个物理量B取动量。系统的总动量是∫ρV dV所以bV。牛顿第二定律说系统内动量随时间的变化率等于作用在系统上的合外力。将输运方程套上这个关系得到d/dt ∫_{sys} ρV dV ∫_{V_CV} ∂(ρV)/∂t dV ∮_{A_CV} ρV (V·n) dA ΣF外力ΣF可以分为体积力和表面力两类。体积力最常见的是重力ρg它是作用在每个流体微团上的。表面力则来自控制面上的应力单位面积上的表面力叫应力向量。表面力可以继续分解为法向应力和切向应力。写成应力张量的形式表面力为∮_A τ_s · n dA。为了从控制体形式的动量方程推出N-S方程需要处理这个应力张量。对于牛顿流体应力包括压强和黏性应力两部分。压强p对任何方向都产生法向的压缩作用黏性应力τ则与流体变形率有关。于是表面力可以写成∮ [-pn τ · n] dA代入前面的输运方程并再次使用高斯散度定理把面积分转换为体积分忽略一些繁琐的推导得到微分形式的动量方程∂(ρV)/∂t ∇·(ρVV) ρg − ∇p ∇·τ这里的∇·(ρVV)是一个二阶张量的散度它看起来不像N-S方程里常见的形式。我们可以用连续性方程把它化简。先把∇·(ρVV)展开为ρ(V·∇)V V∇·(ρV)再用连续性方程∂ρ/∂t ∇·(ρV) 0消去包含∇·(ρV)的项再把∂(ρV)/∂t展开为ρ ∂V/∂t V ∂ρ/∂t。经过整理方程左边变成ρ(∂V/∂t V·∇V)也就是ρDV/Dt这就是我们熟悉的物质导数形式。到这里动量方程已经有了N-S方程的左半部分。右半部分对不可压缩的牛顿流体黏性应力与应变率张量成正比τ μ(∇V (∇V)ᵀ)再经过张量运算化简∇·τ恰好等于μ∇²V。把压强梯度项和重力项放进去ρ DV/Dt ρg − ∇p μ∇²V这就是不可压缩牛顿流体的N-S方程。整个过程没有任何需要背的魔法唯一用到的就是动量守恒加上牛顿流体的本构关系。那条看似复杂的方程本质上说的就是流体微团的加速度乘上密度等于重力、压力差以及黏性力三者的合力。3.3 N-S方程各项的物理意义速记为了防止你推完之后转头就忘我把N-S方程每一项的物理含义整理成一个速查表方程项物理含义类比ρDV/Dt单位体积流体随流运动的惯性力你骑自行车加速时需要克服的惯性ρg单位体积的重力地球对水持续往下拽的力−∇p单位体积的压力梯度力高压区往低压区推动流体的力μ∇²V单位体积的黏性扩散力流体内部试图抹平速度差的“摩擦力”你可以这样记N-S方程就是“惯性力 重力 压力梯度力 黏性力”在流体微团上的牛顿第二定律。以后只要写出这三类力方程自然就出来了不需要从记忆里翻找那一长串偏导符号。4. 实际操作中的避坑指南和问题排查推到这一步你已经拿到了N-S方程。但说实话考试和工程里真正的难点不是把这个方程写出来而是正确使用它。下面我把自己在实际学习和推导中踩过的坑、见过别人踩的坑集中说一下。4.1 六个最常见的坑坑一把控制体当系统。很多人推导时不自觉地用了系统的运动方程却把积分范围写成了固定控制体。比如算一个喷管的受力明明控制体固定不动却跟着流体微团写DV/Dt忽略了进出控制面的动量通量结果力少算了一大截。正确做法是先明确控制体固定在哪再用输运方程把通量项老老实实写上。坑二正负号全靠猜。输运方程里外法线方向n是关键。我见过有人为了省事把所有V·n都当成正速度结果把流入算成了流出得到负的流量还不自知。建议每次写面积分之前在控制面示意图上画一个外法线箭头再判断V和它的夹角。夹角的余弦值超过90度就是负的不会出错。坑三定常流动就直接删掉当地导数。定常流动意味着∂/∂t 0但V·∇项往往还在。比如定常流过弯管流体质点仍然有转弯加速度这就是对流加速度你不能因为“定常”两个字就把它去掉。N-S方程里DV/Dt包含两项一删就错。坑四不可压缩假设被滥用。空气在低速时可以被当作不可压缩流体但一旦马赫数超过0.3密度变化就不能忽略。水在极高压力下密度也会变化。你在推方程时没有考虑能量方程所以所有推导只适用于条件合适的场景。用之前先问一句能不能假设不可压缩牛顿流体的假设还成立吗坑五面对湍流直接用N-S方程。N-S方程本身对层流和湍流瞬时状态都成立但真实湍流场是随机脉动的直接求解每一个瞬时脉动工程上几乎算不动。工程上用的是RANS方程即对N-S方程做时间平均平均后会多出一个雷诺应力项。如果你拿层流牛顿流体的本构关系硬套湍流问题物理上就错了。很多初学者卡在这里不是数学不行而是没搞懂模型适用范围。坑六参考系写错。输运方程和N-S方程都默认在惯性参考系里。如果你在旋转机械、飞机转弯等非惯性系里使用凭空冒出了科里奥利力和离心力方程左边或右边就必须额外加项。否则你算出的流场会在虚假的力作用下错得离谱。4.2 推导完成后的自查清单每次推导完方程我都会做三件事你也可以照做。第一查量纲。量纲是最后一道防火墙。N-S方程左边ρDV/Dt的单位是N/m³右边每一项的单位也都必须是N/m³。如果哪一项出现了Pa·s/m²这种混合单位八成是黏性项错了。量纲不对方程一定错。第二看退化。把速度取零静止流体中N-S方程退化成0 ρg − ∇p也就是静水压强公式。把黏性系数取零退化成欧拉方程对应无黏流动。把∂/∂t取零且速度均匀退化成00。如果退化以后跟已知结论对不上推导里一定藏了问题。第三用简单流动验算。比如均匀流看对流项是否为零一维管流看方程能不能还原成沿流向的动量平衡。我习惯在拿到一个复杂解之前先拿一个最简单的流动去“刺”一下方程刺不破才敢放心用。4.3 一个典型控制体算例收缩管道受力理论说得再多不如动手算一算。我给你一个特别简单的例子帮你建立“控制体→输运方程→受力计算”的完整链路。假设一段水平收缩管道进口截面积A₁出口截面积A₂进口速度V₁出口速度V₂流体不可压缩、定常、忽略壁面摩擦力。我们要算流体给管壁的轴向力。先选控制体取管壁内侧和进出口截面围起来的整个管道内部空间控制面就是入口截面、出口截面和管壁内表面。连续性方程告诉我们A₁V₁ A₂V₂所以V₂ (A₁/A₂)V₁。定常不可压缩时∂/∂t项为零动量方程沿轴向为通过进口进入的动量通量减去出口出来的动量通量等于合外力在轴向的分量。由于出口截面积小、速度大出口动量通量高于进口动量通量合外力的方向必须沿着流向把流体往出口方向加速。这个力就来自管壁对流体施加的压力合力。反过来流体给管壁的力大小相等、方向相反。你不用背任何复杂公式只要写出动量通量差再减去进出口压力差就能得到作用在管壁上的净力。这就是控制体思想的日常用法。4.4 从N-S方程看真实工程问题也许你会问现在计算机都能直接求解N-S方程理解推导还有意义吗我的看法是越依赖数值模拟越需要理解方程来源。因为你在用计算流体力学软件时你选的模型、边界条件、湍流模型、网格分辨率全都是围绕N-S方程展开的。如果你不知道方程两边每一项的物理意义根本没法判断计算结果到底可不可信。比如你在算一个高雷诺数绕流问题如果算出的压力场在某个区域异常你可能需要检查是压力梯度项的问题还是黏性项在网格尺度上被离散误差吃掉了。又比如你做多相流模拟两相间的动量交换本质上还是通过控制体通量来描述追根溯源还是输运方程。工具越来越黑箱但你心里要有一张白纸把方程推一遍白纸上的逻辑才不会丢。我个人在实际操作中最深的体会是输运方程不是流体力学里的一章普通内容它是整套理论的地基。地基不牢你后面每学一个方程都要靠死记死记的东西一定会在关键时刻背叛你。而一旦你真正把控制体、系统、通量这几个概念串起来你会发现空气动力学、水力学、多相流、甚至传热学里的很多方程都可以用同一套逻辑去重新理解。那才是真正“学通了”的感觉。最后再分享一个小技巧如果你手边有N-S方程却怎么都记不住不妨把它的左边改成“ρDV/Dt”右边依次写下“体积力、压力梯度、黏性扩散”然后去想每个力在真实流动里是怎么作用的。比如吹风机吹乒乓球重力往下拉压力差让球往低压区偏黏性让周围空气跟着球走。想一遍这些场景比抄十遍方程都管用。
1. 从“控制体”视角看流体N-S方程不再是天书我学流体力学的时候也干过一件蠢事把N-S方程当英语单词先抄三遍再背三遍结果考完就忘做题照样卡壳。后来才想明白问题不在于记性而在于我从一开始就选错了理解方式。N-S方程不是需要背诵的结论它完全可以从一个更朴素、更直观的视角推出来这个视角就是“控制体”。控制体是什么你可以把它理解成在流场里画一个虚构的盒子盒子可以很小也可以包围整个机翼、整个管道它不随流体运动固定在你关心的空间里。流体从盒子流进流出我们在盒子表面上记录进出的量再结合盒子内部的积累就能写出各种守恒定律。工程上绝大多数计算管道、泵、风机、飞行器外流场都是站在这种固定盒子的视角去看问题。而流体力学教材里那些让人头大的方程包括N-S方程本质上都是从这个视角出发一步步用守恒定律推出来的。这篇文章适合两类人第一类是刚开始学流体力学、正被N-S方程折磨的学生第二类是工作后想回头把底层逻辑补上的工程师。全文不要求你先背任何公式我带着你从零开始先推导输运方程再用输运方程推出连续性方程和N-S方程。推完之后你会发现N-S方程只是“动量守恒定律在控制体上的展开”想忘记都难。1.1 两种观察方式追着流体跑还是守着空间看流体力学里一直有两种观察方式很多人一开始没分清楚后面全乱套。第一种叫系统视角。系统指的是“一堆固定的流体微团”你跟着这堆流体一起走它变形、移动、旋转你都跟着。就像你骑着一匹马去追踪一条河里的同一个水团这个水团从上游漂到下游你始终盯着它。物理定律天然是在这种视角下成立的牛顿第二定律说的是“这个流体团的动量变化率等于所受合力”而不是“这个固定盒子里的动量变化率”。第二种叫控制体视角。控制体是你预先选定的一块空间区域边界叫控制面。流体穿过控制面进进出出你站在控制体外只关心这块区域内部和边界上的变化。这就像你在河上的一座桥墩旁边圈出一段水面不跟着水跑只记录单位时间内有多少水从上游进来、多少水从下游出去、桥墩这一段里水位涨了还是落了。工程上哪种更实用控制体。因为你要设计的管道、阀门、风洞、机翼都是固定的几何边界你关心的是“流过这个截面有多少流量”“作用在这个物体上的力有多大”。你总不能每算一个力就派个人骑着马追着同一个水团跑遍全流场吧。所以实际计算通常用控制体但物理定律本质上是给系统用的这就产生了一个矛盾定律说系统的事计算要用控制体两边怎么对接对接工具就是输运方程它回答一个问题系统内某个物理量随时间的变化率如何用控制体上的当地导数和对流通量来表示。N-S方程就是从这个对接过程中一步步长出来的。1.2 输运方程到底解决了什么问题先举一个最简单的问题一个浴缸放水水位是涨还是落浴缸就是控制体水龙头流进来的水、排水口流出去的水、浴缸里水量的变化三者关系一眼就能写出来水位变化 流进来的水量 − 流出去的水量。输运方程干的也是这件事只不过把“水量”换成了质量、动量、能量这类一般化的物理量。你可能会觉得这有什么好推的但流体不是刚体它无时无刻不在变形、旋转、穿越边界。同一时刻控制体内可能有流体在加速、有流体在流出、有流体在流入物理量在时间的推进中变化得非常复杂。如果直接写出系统导数和控制体积分之间的关系很容易漏项、错项。输运方程的价值就是给这个关系一个严格的、通用的形式。只要你确定你关心的是哪个量把它代进去各种守恒方程都能出来。所以我的建议很简单别背N-S方程先背输运方程再学会代公式。后面你会看到N-S方程的每一项都可以在输运方程里找到对应位置而不是凭空冒出来的。2. 手把手推导输运方程加了这么多铺垫现在就开始正题。我们要推导的是一个关于任意物理量B的方程。先定义清楚B是一个广延量比如质量、动量、能量这样跟物质的量成正比的总量b是单位质量对应的B专业点叫强度量。比如B是质量b就等于1B是动量b就是速度VB是动能b就是v²/2。B和b的关系通过密度联系B ∫ ρb dV这个积分可以放在系统上也可以放在控制体上区别只在于积分区域不同。我们现在先写成系统的形式B_sys(t) ∫_{sys(t)} ρb dV其中sys(t)表示t时刻系统占据的体积。系统的特点是它永远由同一堆流体微团组成随流场运动但形状和位置都会变。我们要算的是系统量随时间的变化率即 d B_sys / d t。这就是牛顿第二定律左边那种导数。输运方程的目标是把这种“跟随流体”的导数改写成“站在固定控制体上看”的表达式。2.1 最直观的做法对比两个相邻时刻推导的关键在于选t时刻让系统刚好和控制体重合。这只是一个数学上的约定因为物理定律里的系统是任选的此刻它与控制体重合方便我们做转换。t时刻系统占据的体积就是控制体V_CV。但t Δt时刻系统里的流体微团已经跟着流场往前走了它们会有一部分流出原来的控制体同时控制体外面的另一部分流体并没有进入系统。注意系统是那堆固定的流体微团不是控制体内的所有流体。所以t Δt时刻系统的体积可以分成三块仍然留在控制体内的部分从控制面流出去的部分原先在控制体内部、现在依然在那里的那部分里的流体依然属于系统。等一下这里要小心。系统和控制体在t时刻重合但在t Δt时刻系统可能有一部分已经流出控制体边界所以系统体积会偏离控制体。我们推导时考虑的是系统量的变化由“控制体内当地变化”加上“穿越控制面的净流出”共同组成。用数学表达系统量的变化率近似写成d B_sys / dt ≈ [B_sys(t Δt) − B_sys(t)] / Δt其中B_sys(t) ∫_{V_CV} ρb dV因为t时刻系统与控制体重合。B_sys(t Δt)则分成两部分仍然在控制体内的量加上流出控制面的量。经过整理并取Δt → 0的极限就得到输运方程的标准形式d/dt ∫_{sys} ρb dV ∫_{V_CV} ∂(ρb)/∂t dV ∮_{A_CV} ρb (V·n) dA这里V是流体在控制面当地的速度n是控制面外法线方向单位向量V·n表示流体沿外法线方向穿越控制面的速度分量。这个式子就是“系统量的变化率 控制体内当地积累率 控制面上净对流输出率”。为了让你彻底记住我用自己的话翻译一遍系统里某物理量的总量变化等于你站在控制体里看到的当地变化再加上从边界上像传送带一样把物理量送出去的量。传出去的是正的贡献传进来的是负的贡献因为外法线方向是朝外的。就这么简单。2.2 通量项的正负号为什么这么重要很多人在这一步栽跟头问题都出在V·n的符号上。我多说几句。控制面的外法线n是垂直于表面、指向控制体外侧的向量。流体速度V如果顺着n的方向流出V·n就是正的代表物理量正在离开控制体对控制体来说是“损失”所以在方程右侧加号里它贡献为正。流体如果从外部流入V·n是负的代表物理量进入控制体对控制体来说是“获得”所以该项是负贡献。净流出的通量项取整个闭合面积分流出的正贡献和流入的负贡献会自然抵消一部分剩下的就是真正从控制体净流出去的量。我见过很多初学的人因为嫌弃V·n这个点积太数学就自己改成“流出减流入”结果在非平面控制面上算错。建议你无论画什么样的控制体都坚持写ρb(V·n)dA用外法线方向自己判断正负而不是背什么“进为正出为负”的口诀。前者是按定义推后者是凭感觉猜猜一定会出错。再举一个生活化的例子。想象你的手机话费是一个保单你先是这个保单的“系统”每个月你的话费变化等于你充进去的钱减去你花掉的钱。控制体就是你的手机账户余额所在的那个“账面”充值是流入消费是流出。输运方程说的就是你这个“系统里的钱”其实可以理解成你名下某个号码的话费余额的变化等于你最近充值的钱加花掉的钱。只不过在流体里“充值”和“消费”每时每刻都在控制面上发生而且不同位置的速率还不一样。2.3 从系统导数到物质导数输运方程更严谨的数学形式还可以用“物质导数”的概念来理解。我们常看到D/Dt这个符号它就是物质导数意思是“跟随一个流体微团去观察某个量随时间的变化率”。在N-S方程里D/Dt ∂/∂t V·∇前面是当地导数后面是对流导数。如果你把输运方程里的系统导数也写成这种形式会发现它恰恰就是把物质导数、当地导数和对流项串在了一起。∂(ρb)/∂t是控制体内某一点上物理量的当地变化V·∇(ρb)则是流体把物理量从一点搬运到另一点带来的变化。这两者合起来才构成跟随流体时看到的完整变化。这也是为什么N-S方程里总会出现V·∇V这样的项它不是额外的复杂项而是系统在流动过程中“随波逐流”的加速度的一部分。3. 从输运方程推守恒方程直到N-S输运方程本身只是工具真正的重头戏是把它应用到物理定律上。严格来说N-S方程不是单独一个方程而是连续性方程加上动量方程在特定假设下的组合。下面我们就按这个顺序一步一步推。3.1 连续性方程取B质量质量守恒是流体力学最基本的定律它说系统的质量不会凭空产生或消失。也就是说d/dt ∫_{sys} ρ dV 0。对质量这个量b1。把b1代入输运方程得到∫_{V_CV} ∂ρ/∂t dV ∮_{A_CV} ρ(V·n) dA 0把面积分用高斯散度定理转换成体积分这里不展开讲高斯的证明你可以把它理解成“边界穿出的通量等于内部散度的总和”∫_{V_CV} [∂ρ/∂t ∇·(ρV)] dV 0由于控制体可以任意选取要让任意体积内的积分恒等于零只能是括号里的被积函数处处为零∂ρ/∂t ∇·(ρV) 0这就是连续性方程。每一项的物理意义第一项是当地密度随时间的变化率第二项是密度通量的散度说白了就是单位体积内质量“向外发散”的速率。两者相加为零意味着质量既不会多也不会少。对不可压缩流体密度沿路径不变并且通常取为常数。连续性方程退化成∇·V 0也就是速度场的散度为零。这时的直观含义是流体不能在某一点堆积你从这边流进来多少体积就得从那边挤出去多少体积就像地铁车厢门关闭瞬间人挤人但总量不变。这个式子看着很简单却是后面N-S方程简化的重要基础。3.2 动量方程取B动量现在换一个物理量B取动量。系统的总动量是∫ρV dV所以bV。牛顿第二定律说系统内动量随时间的变化率等于作用在系统上的合外力。将输运方程套上这个关系得到d/dt ∫_{sys} ρV dV ∫_{V_CV} ∂(ρV)/∂t dV ∮_{A_CV} ρV (V·n) dA ΣF外力ΣF可以分为体积力和表面力两类。体积力最常见的是重力ρg它是作用在每个流体微团上的。表面力则来自控制面上的应力单位面积上的表面力叫应力向量。表面力可以继续分解为法向应力和切向应力。写成应力张量的形式表面力为∮_A τ_s · n dA。为了从控制体形式的动量方程推出N-S方程需要处理这个应力张量。对于牛顿流体应力包括压强和黏性应力两部分。压强p对任何方向都产生法向的压缩作用黏性应力τ则与流体变形率有关。于是表面力可以写成∮ [-pn τ · n] dA代入前面的输运方程并再次使用高斯散度定理把面积分转换为体积分忽略一些繁琐的推导得到微分形式的动量方程∂(ρV)/∂t ∇·(ρVV) ρg − ∇p ∇·τ这里的∇·(ρVV)是一个二阶张量的散度它看起来不像N-S方程里常见的形式。我们可以用连续性方程把它化简。先把∇·(ρVV)展开为ρ(V·∇)V V∇·(ρV)再用连续性方程∂ρ/∂t ∇·(ρV) 0消去包含∇·(ρV)的项再把∂(ρV)/∂t展开为ρ ∂V/∂t V ∂ρ/∂t。经过整理方程左边变成ρ(∂V/∂t V·∇V)也就是ρDV/Dt这就是我们熟悉的物质导数形式。到这里动量方程已经有了N-S方程的左半部分。右半部分对不可压缩的牛顿流体黏性应力与应变率张量成正比τ μ(∇V (∇V)ᵀ)再经过张量运算化简∇·τ恰好等于μ∇²V。把压强梯度项和重力项放进去ρ DV/Dt ρg − ∇p μ∇²V这就是不可压缩牛顿流体的N-S方程。整个过程没有任何需要背的魔法唯一用到的就是动量守恒加上牛顿流体的本构关系。那条看似复杂的方程本质上说的就是流体微团的加速度乘上密度等于重力、压力差以及黏性力三者的合力。3.3 N-S方程各项的物理意义速记为了防止你推完之后转头就忘我把N-S方程每一项的物理含义整理成一个速查表方程项物理含义类比ρDV/Dt单位体积流体随流运动的惯性力你骑自行车加速时需要克服的惯性ρg单位体积的重力地球对水持续往下拽的力−∇p单位体积的压力梯度力高压区往低压区推动流体的力μ∇²V单位体积的黏性扩散力流体内部试图抹平速度差的“摩擦力”你可以这样记N-S方程就是“惯性力 重力 压力梯度力 黏性力”在流体微团上的牛顿第二定律。以后只要写出这三类力方程自然就出来了不需要从记忆里翻找那一长串偏导符号。4. 实际操作中的避坑指南和问题排查推到这一步你已经拿到了N-S方程。但说实话考试和工程里真正的难点不是把这个方程写出来而是正确使用它。下面我把自己在实际学习和推导中踩过的坑、见过别人踩的坑集中说一下。4.1 六个最常见的坑坑一把控制体当系统。很多人推导时不自觉地用了系统的运动方程却把积分范围写成了固定控制体。比如算一个喷管的受力明明控制体固定不动却跟着流体微团写DV/Dt忽略了进出控制面的动量通量结果力少算了一大截。正确做法是先明确控制体固定在哪再用输运方程把通量项老老实实写上。坑二正负号全靠猜。输运方程里外法线方向n是关键。我见过有人为了省事把所有V·n都当成正速度结果把流入算成了流出得到负的流量还不自知。建议每次写面积分之前在控制面示意图上画一个外法线箭头再判断V和它的夹角。夹角的余弦值超过90度就是负的不会出错。坑三定常流动就直接删掉当地导数。定常流动意味着∂/∂t 0但V·∇项往往还在。比如定常流过弯管流体质点仍然有转弯加速度这就是对流加速度你不能因为“定常”两个字就把它去掉。N-S方程里DV/Dt包含两项一删就错。坑四不可压缩假设被滥用。空气在低速时可以被当作不可压缩流体但一旦马赫数超过0.3密度变化就不能忽略。水在极高压力下密度也会变化。你在推方程时没有考虑能量方程所以所有推导只适用于条件合适的场景。用之前先问一句能不能假设不可压缩牛顿流体的假设还成立吗坑五面对湍流直接用N-S方程。N-S方程本身对层流和湍流瞬时状态都成立但真实湍流场是随机脉动的直接求解每一个瞬时脉动工程上几乎算不动。工程上用的是RANS方程即对N-S方程做时间平均平均后会多出一个雷诺应力项。如果你拿层流牛顿流体的本构关系硬套湍流问题物理上就错了。很多初学者卡在这里不是数学不行而是没搞懂模型适用范围。坑六参考系写错。输运方程和N-S方程都默认在惯性参考系里。如果你在旋转机械、飞机转弯等非惯性系里使用凭空冒出了科里奥利力和离心力方程左边或右边就必须额外加项。否则你算出的流场会在虚假的力作用下错得离谱。4.2 推导完成后的自查清单每次推导完方程我都会做三件事你也可以照做。第一查量纲。量纲是最后一道防火墙。N-S方程左边ρDV/Dt的单位是N/m³右边每一项的单位也都必须是N/m³。如果哪一项出现了Pa·s/m²这种混合单位八成是黏性项错了。量纲不对方程一定错。第二看退化。把速度取零静止流体中N-S方程退化成0 ρg − ∇p也就是静水压强公式。把黏性系数取零退化成欧拉方程对应无黏流动。把∂/∂t取零且速度均匀退化成00。如果退化以后跟已知结论对不上推导里一定藏了问题。第三用简单流动验算。比如均匀流看对流项是否为零一维管流看方程能不能还原成沿流向的动量平衡。我习惯在拿到一个复杂解之前先拿一个最简单的流动去“刺”一下方程刺不破才敢放心用。4.3 一个典型控制体算例收缩管道受力理论说得再多不如动手算一算。我给你一个特别简单的例子帮你建立“控制体→输运方程→受力计算”的完整链路。假设一段水平收缩管道进口截面积A₁出口截面积A₂进口速度V₁出口速度V₂流体不可压缩、定常、忽略壁面摩擦力。我们要算流体给管壁的轴向力。先选控制体取管壁内侧和进出口截面围起来的整个管道内部空间控制面就是入口截面、出口截面和管壁内表面。连续性方程告诉我们A₁V₁ A₂V₂所以V₂ (A₁/A₂)V₁。定常不可压缩时∂/∂t项为零动量方程沿轴向为通过进口进入的动量通量减去出口出来的动量通量等于合外力在轴向的分量。由于出口截面积小、速度大出口动量通量高于进口动量通量合外力的方向必须沿着流向把流体往出口方向加速。这个力就来自管壁对流体施加的压力合力。反过来流体给管壁的力大小相等、方向相反。你不用背任何复杂公式只要写出动量通量差再减去进出口压力差就能得到作用在管壁上的净力。这就是控制体思想的日常用法。4.4 从N-S方程看真实工程问题也许你会问现在计算机都能直接求解N-S方程理解推导还有意义吗我的看法是越依赖数值模拟越需要理解方程来源。因为你在用计算流体力学软件时你选的模型、边界条件、湍流模型、网格分辨率全都是围绕N-S方程展开的。如果你不知道方程两边每一项的物理意义根本没法判断计算结果到底可不可信。比如你在算一个高雷诺数绕流问题如果算出的压力场在某个区域异常你可能需要检查是压力梯度项的问题还是黏性项在网格尺度上被离散误差吃掉了。又比如你做多相流模拟两相间的动量交换本质上还是通过控制体通量来描述追根溯源还是输运方程。工具越来越黑箱但你心里要有一张白纸把方程推一遍白纸上的逻辑才不会丢。我个人在实际操作中最深的体会是输运方程不是流体力学里的一章普通内容它是整套理论的地基。地基不牢你后面每学一个方程都要靠死记死记的东西一定会在关键时刻背叛你。而一旦你真正把控制体、系统、通量这几个概念串起来你会发现空气动力学、水力学、多相流、甚至传热学里的很多方程都可以用同一套逻辑去重新理解。那才是真正“学通了”的感觉。最后再分享一个小技巧如果你手边有N-S方程却怎么都记不住不妨把它的左边改成“ρDV/Dt”右边依次写下“体积力、压力梯度、黏性扩散”然后去想每个力在真实流动里是怎么作用的。比如吹风机吹乒乓球重力往下拉压力差让球往低压区偏黏性让周围空气跟着球走。想一遍这些场景比抄十遍方程都管用。
