流体力学劝退第一人非N-S方程莫属。我当年上课时老师用一种平静而坚定的语气说这个方程大家记住。然后一整墙的偏微分符号就挂在了黑板上。期中考试默写加解释全班过半靠考前突击。可突击完除了会写那几个符号真到做项目、调CFD仿真软件的时候根本不知道每一项在干什么。后来我换了个思路不从方程本身入手而是回到它的源头——从控制体的视角先把输运定理推明白再把动量守恒一步步走通。结果发现N-S方程根本不是背出来的它是顺着物理直觉自然长出来的。今天这篇就带你重新走一遍这条路顺便把推导过程中我踩过的坑都摊开说一说。1. 从“背公式”到“看机制”先建立两种视角1.1 系统 vs 控制体一个像跟拍一个像监控流体力学里有两种看待流体的方式这是整个推导的地基也是绝大多数人没搞懂的地方。第一种叫系统System也叫物质体。它的定义是始终锁定同一批流体粒子跟着它们走。想象你手里拿着摄像机跟拍一位跑步的人无论他跑到哪里、摆出什么姿势镜头都跟着他。系统的边界是物质边界可以变形、移动但系统里包含的物质永远不会换——不增不减。第二种叫控制体Control Volume它是一个固定的空间区域边界叫控制面Control Surface。流体粒子可以自由进出我们只关心“在某个空间区域里发生什么”。这就像小区门口的监控摄像头位置固定记录进出的人流——有人进来、有人出去但它不关心某人昨天住哪、明天去哪。两种视角各有利弊。系统视角天然契合物理定律因为牛顿第二定律、质量守恒这些定律本来就是针对“同一团物质”写的。但它不实用工程上你没法给每一团流体粒子装GPS况且我们关心的往往是泵、管道、涡轮、翼型周围这些固定空间里的流动而不是某一团粒子的旅程。控制体视角更符合工程习惯但它的麻烦在于——物理定律不是直接作用在控制体上的因为控制体里的物质一直在换。那怎么办答案就是输运定理。1.2 输运定理两个视角之间的“翻译器”输运定理要解决的问题用一句话说清楚跟随物质运动时某个物理量B的变化率等于控制体内B的变化率加上通过控制面净流出的B的速率。这不是一个需要背的孤立公式它本质上就是一道会计平衡账变化总量 内部累积 净流出。我常用一个极其接地气的类比一个小区的人数变化。系统视角是给每个居民装定位算这群人总人数的变化控制体视角是在小区大门安装闸机数进进出出的人流。如果小区里没有新生儿、没有死亡没有源项那么“这群老住户的总人数变化”正好等于“闸机口出去的人减去进来的人”。输运定理把这笔账精确写成了数学公式。为什么要费这么大劲因为流体力学的三大守恒律——质量、动量、能量——在最底层都是针对系统写的。而我们工程上需要的是控制体方程。输运定理正是把“系统的物理定律”翻译成“控制体的可用形式”的那座桥。理解这一点后面的推导就有方向感了先是物理学再是数学翻译最后才是方程。1.3 广延量与强度量先给物理量分个类展开推导之前一个基本概念必须先说清楚广延量Extensive Quantity和强度量Intensive Quantity。广延量的大小和物质总量成正比。一桶水有2千克两桶同样温度的水混一起有4千克加量就加这就是广延量。常见的广延量有质量、动量、动能、内能、熵。而强度量是单位质量的广延量它不随物质多少改变比如温度、密度、速度、单位质量的内能。关键联系是如果B是一个广延量那么单位质量的B记为b就有[ B \int_{\text{物质体}} \rho b , dV ]为什么要引入b而不直接用B因为在积分形式下物理量在场中的分布可能不均匀。比如动量不同位置速度不同必须用速度场V乘密度再积分。选好b之后输运定理就可以统一处理所有守恒量。这套“分类学”看着不起眼但它决定了整个推导的通用性——你推一次输运定理可以套出质量、动量、能量一整套方程。2. 手把手推导输运定理一页草稿纸搞定2.1 推导的起点让系统与控制体在t时刻重合这是全篇最巧妙的一个几何操作也是你亲手推导时最容易卡住的地方。设某一时刻为t。我们让系统在这个时刻恰好占据控制体CV。也就是说t时刻系统边界和控制面完全重合。这一招的意义在于t时刻系统内的物理量B可以直接写成控制体内的积分[ B_{sys}(t) B_{CV}(t) \int_{CV} \rho b , dV ]但过了Δt之后系统随着流体运动走了它不再和控制体重合——一部分原来的物质流出控制面控制体内又涌进了新物质。正因为t时刻两者重合我们才能把“系统的变化率”和“控制体的变化率”干净地联系起来。如果没有这一步重合后面根本没有可比对的对象。2.2 拆解Δt内的变化控制体积累 流出 − 流入现在看tΔt时刻系统内的B怎么表达。t时刻系统里的那团物质在Δt后分成了几部分一部分仍然留在控制体内一部分已经流出控制面。 而控制体内新流入的物质并不属于这个系统它是别人的“系统”。所以系统在tΔt时的B可以这么拆[ B_{sys}(t\Delta t) B_{CV}(t\Delta t) B_{out}(t\Delta t) - B_{in}(t\Delta t) ]这里面B_{out}是Δt内流出控制面的那部分BB_{in}是Δt内流入控制面的那部分B。注意系统在tΔt时不再包含流入的物质所以B_{in}要减去系统在tΔt时大于控制体多了流出的部分所以B_{out}要加上。把t时刻的B_{sys}(t)B_{CV}(t)移到左边除以Δt再取Δt→0极限[ \frac{dB_{sys}}{dt} \frac{\partial B_{CV}}{\partial t} \lim_{\Delta t \to 0} \frac{B_{out}(\Delta t)}{\Delta t} - \lim_{\Delta t \to 0} \frac{B_{in}(\Delta t)}{\Delta t} ]第一项很好理解控制体内B随时间的变化率写成偏导数因为在固定空间观察。后两项就是穿过控制面的净通量。2.3 通量项的几何表达V·n出场关键来了怎么把一个Δ t内穿过面积dA的B表达出来考虑控制面上一个微元面dA它的外法线单位向量记为n。流体相对控制面的速度为V。那么在Δt内穿过dA的流体体积是一个以dA为底、以VΔt为斜高的“斜柱体”。这个斜柱体的体积等于[ (V \cdot n) , \Delta t , dA ]看V·n的正负如果流速方向和外法线方向一致往外流V·n0这是流出反之往内流时V·n0自动带上负号。这一个点积符号解决了所有方向问题不需要分“出口面”“入口面”分别讨论全部统一成一个面积分。这块体积对应的B为[ \rho b , (V \cdot n) , \Delta t , dA ]除以Δt再对整个控制面积分就是后两项的极限合并[ \int_{CS} \rho b , (V \cdot n) , dA ]所以雷诺输运定理的积分形式出来了[ \frac{dB_{sys}}{dt} \frac{\partial}{\partial t} \int_{CV} \rho b , dV \int_{CS} \rho b , (V \cdot n) , dA ]我在第一次推导时总犯一个错误把“流出”的项写成正号、“流入”的项写成负号然后在面积分里又搞一次方向判断结果正负号翻倍。其实只要记住一个原则面积分里统一用外法线nV·n自然区分进出前面的正负号不要动。这个约定能帮你省掉至少一半的符号错误。2.4 从积分形式到微分形式高斯定理救场积分形式很有用CFD里的有限体积法直接用它。但N-S方程更常见的写法是微分形式需要再做一步转换。把通量面积分用高斯散度定理改写[ \int_{CS} \rho b , (V \cdot n) , dA \int_{CV} \nabla \cdot (\rho b V) , dV ]代回去输运定理变成[ \frac{dB_{sys}}{dt} \int_{CV} \left[ \frac{\partial (\rho b)}{\partial t} \nabla \cdot (\rho b V) \right] , dV ]因为控制体是任意取的想让这个积分式对任意控制体都成立唯一的可能是被积函数处处为零。这就是局部化引理——一个从积分形式跳到微分形式的“万能跳板”。于是[ \frac{\partial (\rho b)}{\partial t} \nabla \cdot (\rho b V) 0 ]这个式子堪称流体力学领域的“母方程”。它本身没有任何物理内容只是一个数学模板。质量、动量、能量方程全都长这样区别只是各自带不带源项、b具体是什么。掌握它等于握住了一把通用钥匙。2.5 第一站验证质量守恒方程拿到模板第一个要复现的当然是最简单的质量守恒。质量对应的广延量是Bm那么b1每千克物质含1千克质量。代入“母方程”得到[ \frac{\partial \rho}{\partial t} \nabla \cdot (\rho V) 0 ]这就是连续性方程。物理含义清晰某点密度随时间增大的速度加上质量通量从这个点向外发散的速度合计为零。如果没有源没有化学反应生成质量那密度变化只能靠通量来补。工程上常见的简化场景恒定流动∂ρ/∂t 0得∇·(ρV)0一维管道中就是ρ₁V₁A₁ρ₂V₂A₂流量守恒不可压缩流动ρ恒定更进一步退化成∇·V0意思是流场速度无散度像橡皮泥一样体积不可压缩。我自己验证这一步时特意用一维管道算了流量结果和教材完全一致。那种“公式对上了”的感觉就是理解在生根的信号。3. 从输运定理走向N-S方程关键一跃3.1 动量守恒的骨架输运定理 牛顿第二定律质量守恒只是热身接下来才进入正题动量方程。动量的广延量是BmV所以bV。动量守恒的物理内核就是牛顿第二定律一个物质系统的动量变化率等于它受到的合力。[ \frac{dP_{sys}}{dt} \sum F ]把输运定理套上去左边变成[ \int_{CV} \left[ \frac{\partial (\rho V)}{\partial t} \nabla \cdot (\rho V \otimes V) \right] dV \sum F ]这里V⊗V是并矢积张量积因为动量是矢量对流出去的必须是二阶结构才能同时表达“多少动量往哪个方向流”。不少人第一次看到V⊗V会懵其实拆开来看运输量包含了速度分量与流动方向是比标量通量更丰富的量。右边合力F要拆成两类体积力作用在流体内每一个质点上如重力、电磁力用密度表示就是ρg积分得∫_CV ρg dV表面力作用在控制面上包括压力和粘性应力。用应力张量σ描述单位面积上的力是σ·n积分得∫_CS σ·n dA。表面力再用地散度定理一转换[ \int_{CS} \sigma \cdot n , dA \int_{CV} \nabla \cdot \sigma , dV ]于是动量方程积分形式变身为[ \frac{\partial (\rho V)}{\partial t} \nabla \cdot (\rho V \otimes V) \rho g \nabla \cdot \sigma ]这就是柯西运动方程对任何连续介质都成立——还不只是流体固体变形也常用它。它和N-S方程之间只差一个“本构关系”了。3.2 应力张量压强不见了来了个更大的角色很多初学者一听到应力张量就头大。我帮你拆掉这个心理门槛。应力张量σ描述的是控制面上任意一点、任意朝向的面元上单位面积受到的力。流体微团表面有无数个朝向的面每个面受力有三个分量光靠一个标量压强根本不够。什么时候压强够用静止流体、无粘流体。一旦有粘性微生物表面还受切向的摩擦力于是需要一个二阶张量九个分量把各方向受力全装进去。理解σ的关键在于它天然包含两项[ \sigma -pI \tau ]-pI是压强贡献方向永远垂直表面向内压应力τ是偏应力张量来自流体的变形。在牛顿流体中偏应力与变形率线性相关[ \tau \mu \left( \nabla V (\nabla V)^T \right) \lambda (\nabla \cdot V) I ]μ是动力粘性系数λ是第二粘度在斯托克斯假设下取-2μ/3。对于不可压缩流体∇·V0第二项直接消失式子大幅简化[ \tau \mu \left( \nabla V (\nabla V)^T \right) ]3.3 化简的关键手笔用连续性方程收拾左边从柯西方程到N-S方程中间这一步是很多人“背公式”时完全忽略的但它才是理解N-S方程“形状”的钥匙。柯西方程的左边是∂(ρV)/∂t ∇·(ρV⊗V)。直接用链式法则展开第i个分量左边 ρ∂V_i/∂t V_i∂ρ/∂t V_i∂(ρV_j)/∂x_j ρV_j∂V_i/∂x_j组合一下把带∂V_i/∂t和∂V_i/∂x_j的项挑出来剩下的项放到一起左边 ρ(∂V_i/∂t V_j∂V_i/∂x_j) V_i[∂ρ/∂t ∂(ρV_j)/∂x_j]看仔细方括号里是什么那不是连续性方程的左边吗质量守恒告诉我们它等于零。于是整个“多余部分”干干净净地消失左边 ρ(∂V_i/∂t V_j∂V_i/∂x_j)写回矢量形式就是[ \rho \left( \frac{\partial V}{\partial t} (V \cdot \nabla) V \right) ]这一步的深刻含义是质量守恒律“嵌入”了动量方程。所以N-S方程左边那个看似很特别的括号∂V/∂t (V·∇)V其实是物质加速度跟随流体粒子一起运动时感受到的加速度而不是某个固定观测点的变化率。这就是为什么它是Fma里的那个“a”。3.4 加上牛顿流体本构N-S方程诞生不可压缩牛顿流体的本构关系代回柯西方程∇·σ ∇·(-pI τ)其中-∇p项直接出来。τ的散度呢[ \frac{\partial \tau_{ij}}{\partial x_j} \mu \frac{\partial^2 u_i}{\partial x_j \partial x_j} \mu \frac{\partial}{\partial x_i} \left( \frac{\partial u_j}{\partial x_j} \right) ]第一个交叉项正是μ∇²V。第二个括号里是∇·V不可压缩条件下为零。所以[ \rho \left( \frac{\partial V}{\partial t} (V \cdot \nabla) V \right) - \nabla p \mu \nabla^2 V \rho g ]这就是大名鼎鼎的N-S方程。整个过程走下来每个符号都有明确来历没有一个是需要“硬背”的。从输运定理出发加上牛顿第二定律、应力分解、牛顿本构、连续性方程一气呵成。3.5 N-S方程每一项的物理直觉把方程当成Fma读每一项都有对应的受力来源项名称物理含义生活类比ρ(∂V/∂t)局部加速度固定位置上速度随时间的变化站在桥上看同一位置水流忽快忽慢ρ(V·∇)V对流加速度粒子搬运到新位置引起的速度变化从平静河湾漂到急流中央被加速-∇p压力梯度力高压指向低压的力风从高压区吹向低压区μ∇²V粘性力速度分布弯曲导致的耗散力搅动糖浆时周围的糖浆被带动ρg体积力外部力场重力水往低处流这套表格是我在课堂上讲N-S时的保留节目。它让方程不是一个符号堆砌而是一张“力的全家福”。下次看到方程哪怕不看推导你也应该能在脑子里把每一项“翻译”成一股力、一种流动现象。4. 推导中的翻车现场与排查技巧4.1 通量项符号翻车流出到底加还是减这是新手翻车率最高的一处。问题几乎都出在同一件事上既在公式里写了面积分又在脑子里做了一次“流出加、流入减”的判断结果符号被处理了两遍。正确的做法是面积分里统一写∫_CS ρb(V·n)dA前面别再加额外的正负号。V·n为正时流出逆然为正贡献为负时流入自动贡献为负。你先入为主地“手动”减流入就会在入口上搞出双重减号。我的习惯是凡涉及方向判断立即在草稿纸上画出控制体左边一个入口、右边一个出口标出外法线n的箭头。符号判断只依靠V·n的数值正负不做任何直觉正负。这个小习惯救过我无数次。4.2 把(V·∇)和∇·V搞混看起来像物理上是两个东西(V·∇)V和V(∇·V)是完全不同的两项。前者是“速度场作为微分算子去作用于自己的速度场”是N-S方程中典型非线性项后者是“速度乘以速度场的散度”只在密度变化、体积胀缩时才出现。不幸的是它们在张量运算中长得很像缩写得不好甚至分不清。一个可靠的辨认方法V·∇是一个标量算子作用在矢量上仍得到矢量∇·V则是一个标量。你在N-S方程左边看到的、产生非线性效应的那个是(V·∇)V它对应流体粒子经过空间不均匀流场时产生的“搬运加速度”不是体积膨胀率。考试和编程里把这两项写反方程物理性质会彻底变坏。4.3 高斯定理的隐藏前提间断面会让你翻车高斯定理要求被积函数在区域内连续可微。当流场出现激波、自由面等间断结构时微分形式的N-S方程在间断面上不再成立只有积分形式输运定理仍然有效。这就是为什么CFD里所有主流求解器比如有限体积法都在积分形式下做离散而不是硬刚微分方程。实际工程提示跨音速流、气液两相流这类含强间断的场景如果你在校核方程时非要用微分形式去理解必须额外引入Rankine-Hugoniot跳变条件。而用积分形式写瞬时尚则完全不用操心间断。这也是我在学了CFD后才真正体会到“输运定理的积分形式才是流体力学的底层原生形态”。4.4 不可压缩假设什么时候可以偷懒、什么时候不能不可压处理后的N-S方程精简很多密度常数、∇·V0、粘性项退化为μ∇²V。但“不可压缩”不等于“密度不变”那么简单。它的严格含义是流动中压强、温度变化引起的密度相对变化可以忽略。低速液体流动、马赫数低于0.3的气体流动放心用不可压形式。但高速气流、燃烧、激波、高超声速流动密度剧烈变化必须用完整可压缩N-S方程并且要额外补一个能量方程和状态方程来闭合。初学者最典型的问题是把水和低速空气当成永不可压却在跨声速喷管问题上套不可压模型结果数值噪音暴涨物理结果错得离谱。4.5 总在漏掉边界条件即使完整推导出了N-S方程你还需要边界条件才能真正求解。无滑移壁面速度与壁面一致、入口速度/压力条件、出口通量条件、自由表面条件——每一样都直接影响解的形态。我见过不少学完推导的人以为有了方程就有了一切结果跑到仿真软件里眼巴巴盯着残差看却不知边界条件已经给了物理上荒谬的设定。记住方程描述规律边界条件定义场景两者合起来才是完整问题。这个意识比方程本身早点建立更好。5. 把输运定理装进自己的知识体系学习建议与心得5.1 以输运定理为核心建立“守恒律控制表”到了这一步掌握最有力的方式是把之前所有的B、b代换整理成一个控制表守恒定律广延量B强度量b方程形式质量m1∂ρ/∂t∇·(ρV)0动量mVVρ(∂V/∂t(V·∇)V)∇·σρg能量m(e|V|²/2)e|V|²/2∂E/∂t∇·((Ep)V)...这等于告诉你整个流体力学理论树其实是从一个输运定理“长”出来的。质量、动量、能量守恒并不是三个不相干方程而是同一个数学模板在不同物理量上投影的结果。有了这张表你连“为什么能量方程那么复杂”都能猜到大概——因为b里带了内能、动能还做功传热。5.2 学习路径建议先积分、后微分、再上CFD很多人一上课就直奔微分形式N-S方程这是最大的误解。正确顺序是先在积分形式下理解输运定理理解“系统变化率 控制体积累 净流出”这句会计平衡账学会把积分形式转微分形式高斯定理 局部化引理再给特定守恒律配b和源项最后才谈CFD离散理解有限体积法为什么在每条网格边界计算通量。我做CFD项目时有一个体会有限体积法程序里的每一个通量计算单元本质上都在帮你“手算”输运定理中那个∫_CS ρb(V·n)dA项。如果你搞懂了输运定理看CFD代码会特别亲切——完全不像看天书。5.3 一点个人经验推导两遍胜过默写十遍最后分享一个我带人时反复验证过的经验。凡是能把输运定理自己在草稿纸上完整推两遍——第一遍照着推第二遍合上笔记推——的人之后看N-S方程就像看Fma一样自然。那些总想跳过推导直接做题的人通常在各类测验、面试、项目里反复栽跟头因为稍微变形一下题目就识别不出来。具体操作建议第一遍照着笔记或本文流程推理解每一步第二遍合上笔记从“设系统与控制体在t时刻重合”开始独立推完输运定理再推连续性方程和柯西方程第三遍可选加入牛顿本构推到不可压N-S方程。三遍下来你对N-S方程的“熟悉感”会从“背出来了”变成“本该如此”。真到用的时候——写个脚本、调个仿真器、分析实验数据——你会发现自己完全不需要去“回忆”方程因为整个推导图景已经在脑子里长成了本能。顺带说一句别小看控制体的选取技巧。做工业计算时控制体选得好不好直接影响方程的简洁程度。压气机内部流动选随体控制体管道流量计算选固定控制体两相流界面问题选变形控制体——这些都是课本不细讲、但实际项目里天天要用的经验。等你把输运定理真正用熟了这些“怪招”自然会从直觉里长出来。 
长河君 · 2026/10/9 23:47:02 · 阅读 1.2k 从输运定理到N-S方程：告别死记硬背的流体力学推导
流体力学劝退第一人非N-S方程莫属。我当年上课时老师用一种平静而坚定的语气说这个方程大家记住。然后一整墙的偏微分符号就挂在了黑板上。期中考试默写加解释全班过半靠考前突击。可突击完除了会写那几个符号真到做项目、调CFD仿真软件的时候根本不知道每一项在干什么。后来我换了个思路不从方程本身入手而是回到它的源头——从控制体的视角先把输运定理推明白再把动量守恒一步步走通。结果发现N-S方程根本不是背出来的它是顺着物理直觉自然长出来的。今天这篇就带你重新走一遍这条路顺便把推导过程中我踩过的坑都摊开说一说。1. 从“背公式”到“看机制”先建立两种视角1.1 系统 vs 控制体一个像跟拍一个像监控流体力学里有两种看待流体的方式这是整个推导的地基也是绝大多数人没搞懂的地方。第一种叫系统System也叫物质体。它的定义是始终锁定同一批流体粒子跟着它们走。想象你手里拿着摄像机跟拍一位跑步的人无论他跑到哪里、摆出什么姿势镜头都跟着他。系统的边界是物质边界可以变形、移动但系统里包含的物质永远不会换——不增不减。第二种叫控制体Control Volume它是一个固定的空间区域边界叫控制面Control Surface。流体粒子可以自由进出我们只关心“在某个空间区域里发生什么”。这就像小区门口的监控摄像头位置固定记录进出的人流——有人进来、有人出去但它不关心某人昨天住哪、明天去哪。两种视角各有利弊。系统视角天然契合物理定律因为牛顿第二定律、质量守恒这些定律本来就是针对“同一团物质”写的。但它不实用工程上你没法给每一团流体粒子装GPS况且我们关心的往往是泵、管道、涡轮、翼型周围这些固定空间里的流动而不是某一团粒子的旅程。控制体视角更符合工程习惯但它的麻烦在于——物理定律不是直接作用在控制体上的因为控制体里的物质一直在换。那怎么办答案就是输运定理。1.2 输运定理两个视角之间的“翻译器”输运定理要解决的问题用一句话说清楚跟随物质运动时某个物理量B的变化率等于控制体内B的变化率加上通过控制面净流出的B的速率。这不是一个需要背的孤立公式它本质上就是一道会计平衡账变化总量 内部累积 净流出。我常用一个极其接地气的类比一个小区的人数变化。系统视角是给每个居民装定位算这群人总人数的变化控制体视角是在小区大门安装闸机数进进出出的人流。如果小区里没有新生儿、没有死亡没有源项那么“这群老住户的总人数变化”正好等于“闸机口出去的人减去进来的人”。输运定理把这笔账精确写成了数学公式。为什么要费这么大劲因为流体力学的三大守恒律——质量、动量、能量——在最底层都是针对系统写的。而我们工程上需要的是控制体方程。输运定理正是把“系统的物理定律”翻译成“控制体的可用形式”的那座桥。理解这一点后面的推导就有方向感了先是物理学再是数学翻译最后才是方程。1.3 广延量与强度量先给物理量分个类展开推导之前一个基本概念必须先说清楚广延量Extensive Quantity和强度量Intensive Quantity。广延量的大小和物质总量成正比。一桶水有2千克两桶同样温度的水混一起有4千克加量就加这就是广延量。常见的广延量有质量、动量、动能、内能、熵。而强度量是单位质量的广延量它不随物质多少改变比如温度、密度、速度、单位质量的内能。关键联系是如果B是一个广延量那么单位质量的B记为b就有[ B \int_{\text{物质体}} \rho b , dV ]为什么要引入b而不直接用B因为在积分形式下物理量在场中的分布可能不均匀。比如动量不同位置速度不同必须用速度场V乘密度再积分。选好b之后输运定理就可以统一处理所有守恒量。这套“分类学”看着不起眼但它决定了整个推导的通用性——你推一次输运定理可以套出质量、动量、能量一整套方程。2. 手把手推导输运定理一页草稿纸搞定2.1 推导的起点让系统与控制体在t时刻重合这是全篇最巧妙的一个几何操作也是你亲手推导时最容易卡住的地方。设某一时刻为t。我们让系统在这个时刻恰好占据控制体CV。也就是说t时刻系统边界和控制面完全重合。这一招的意义在于t时刻系统内的物理量B可以直接写成控制体内的积分[ B_{sys}(t) B_{CV}(t) \int_{CV} \rho b , dV ]但过了Δt之后系统随着流体运动走了它不再和控制体重合——一部分原来的物质流出控制面控制体内又涌进了新物质。正因为t时刻两者重合我们才能把“系统的变化率”和“控制体的变化率”干净地联系起来。如果没有这一步重合后面根本没有可比对的对象。2.2 拆解Δt内的变化控制体积累 流出 − 流入现在看tΔt时刻系统内的B怎么表达。t时刻系统里的那团物质在Δt后分成了几部分一部分仍然留在控制体内一部分已经流出控制面。 而控制体内新流入的物质并不属于这个系统它是别人的“系统”。所以系统在tΔt时的B可以这么拆[ B_{sys}(t\Delta t) B_{CV}(t\Delta t) B_{out}(t\Delta t) - B_{in}(t\Delta t) ]这里面B_{out}是Δt内流出控制面的那部分BB_{in}是Δt内流入控制面的那部分B。注意系统在tΔt时不再包含流入的物质所以B_{in}要减去系统在tΔt时大于控制体多了流出的部分所以B_{out}要加上。把t时刻的B_{sys}(t)B_{CV}(t)移到左边除以Δt再取Δt→0极限[ \frac{dB_{sys}}{dt} \frac{\partial B_{CV}}{\partial t} \lim_{\Delta t \to 0} \frac{B_{out}(\Delta t)}{\Delta t} - \lim_{\Delta t \to 0} \frac{B_{in}(\Delta t)}{\Delta t} ]第一项很好理解控制体内B随时间的变化率写成偏导数因为在固定空间观察。后两项就是穿过控制面的净通量。2.3 通量项的几何表达V·n出场关键来了怎么把一个Δ t内穿过面积dA的B表达出来考虑控制面上一个微元面dA它的外法线单位向量记为n。流体相对控制面的速度为V。那么在Δt内穿过dA的流体体积是一个以dA为底、以VΔt为斜高的“斜柱体”。这个斜柱体的体积等于[ (V \cdot n) , \Delta t , dA ]看V·n的正负如果流速方向和外法线方向一致往外流V·n0这是流出反之往内流时V·n0自动带上负号。这一个点积符号解决了所有方向问题不需要分“出口面”“入口面”分别讨论全部统一成一个面积分。这块体积对应的B为[ \rho b , (V \cdot n) , \Delta t , dA ]除以Δt再对整个控制面积分就是后两项的极限合并[ \int_{CS} \rho b , (V \cdot n) , dA ]所以雷诺输运定理的积分形式出来了[ \frac{dB_{sys}}{dt} \frac{\partial}{\partial t} \int_{CV} \rho b , dV \int_{CS} \rho b , (V \cdot n) , dA ]我在第一次推导时总犯一个错误把“流出”的项写成正号、“流入”的项写成负号然后在面积分里又搞一次方向判断结果正负号翻倍。其实只要记住一个原则面积分里统一用外法线nV·n自然区分进出前面的正负号不要动。这个约定能帮你省掉至少一半的符号错误。2.4 从积分形式到微分形式高斯定理救场积分形式很有用CFD里的有限体积法直接用它。但N-S方程更常见的写法是微分形式需要再做一步转换。把通量面积分用高斯散度定理改写[ \int_{CS} \rho b , (V \cdot n) , dA \int_{CV} \nabla \cdot (\rho b V) , dV ]代回去输运定理变成[ \frac{dB_{sys}}{dt} \int_{CV} \left[ \frac{\partial (\rho b)}{\partial t} \nabla \cdot (\rho b V) \right] , dV ]因为控制体是任意取的想让这个积分式对任意控制体都成立唯一的可能是被积函数处处为零。这就是局部化引理——一个从积分形式跳到微分形式的“万能跳板”。于是[ \frac{\partial (\rho b)}{\partial t} \nabla \cdot (\rho b V) 0 ]这个式子堪称流体力学领域的“母方程”。它本身没有任何物理内容只是一个数学模板。质量、动量、能量方程全都长这样区别只是各自带不带源项、b具体是什么。掌握它等于握住了一把通用钥匙。2.5 第一站验证质量守恒方程拿到模板第一个要复现的当然是最简单的质量守恒。质量对应的广延量是Bm那么b1每千克物质含1千克质量。代入“母方程”得到[ \frac{\partial \rho}{\partial t} \nabla \cdot (\rho V) 0 ]这就是连续性方程。物理含义清晰某点密度随时间增大的速度加上质量通量从这个点向外发散的速度合计为零。如果没有源没有化学反应生成质量那密度变化只能靠通量来补。工程上常见的简化场景恒定流动∂ρ/∂t 0得∇·(ρV)0一维管道中就是ρ₁V₁A₁ρ₂V₂A₂流量守恒不可压缩流动ρ恒定更进一步退化成∇·V0意思是流场速度无散度像橡皮泥一样体积不可压缩。我自己验证这一步时特意用一维管道算了流量结果和教材完全一致。那种“公式对上了”的感觉就是理解在生根的信号。3. 从输运定理走向N-S方程关键一跃3.1 动量守恒的骨架输运定理 牛顿第二定律质量守恒只是热身接下来才进入正题动量方程。动量的广延量是BmV所以bV。动量守恒的物理内核就是牛顿第二定律一个物质系统的动量变化率等于它受到的合力。[ \frac{dP_{sys}}{dt} \sum F ]把输运定理套上去左边变成[ \int_{CV} \left[ \frac{\partial (\rho V)}{\partial t} \nabla \cdot (\rho V \otimes V) \right] dV \sum F ]这里V⊗V是并矢积张量积因为动量是矢量对流出去的必须是二阶结构才能同时表达“多少动量往哪个方向流”。不少人第一次看到V⊗V会懵其实拆开来看运输量包含了速度分量与流动方向是比标量通量更丰富的量。右边合力F要拆成两类体积力作用在流体内每一个质点上如重力、电磁力用密度表示就是ρg积分得∫_CV ρg dV表面力作用在控制面上包括压力和粘性应力。用应力张量σ描述单位面积上的力是σ·n积分得∫_CS σ·n dA。表面力再用地散度定理一转换[ \int_{CS} \sigma \cdot n , dA \int_{CV} \nabla \cdot \sigma , dV ]于是动量方程积分形式变身为[ \frac{\partial (\rho V)}{\partial t} \nabla \cdot (\rho V \otimes V) \rho g \nabla \cdot \sigma ]这就是柯西运动方程对任何连续介质都成立——还不只是流体固体变形也常用它。它和N-S方程之间只差一个“本构关系”了。3.2 应力张量压强不见了来了个更大的角色很多初学者一听到应力张量就头大。我帮你拆掉这个心理门槛。应力张量σ描述的是控制面上任意一点、任意朝向的面元上单位面积受到的力。流体微团表面有无数个朝向的面每个面受力有三个分量光靠一个标量压强根本不够。什么时候压强够用静止流体、无粘流体。一旦有粘性微生物表面还受切向的摩擦力于是需要一个二阶张量九个分量把各方向受力全装进去。理解σ的关键在于它天然包含两项[ \sigma -pI \tau ]-pI是压强贡献方向永远垂直表面向内压应力τ是偏应力张量来自流体的变形。在牛顿流体中偏应力与变形率线性相关[ \tau \mu \left( \nabla V (\nabla V)^T \right) \lambda (\nabla \cdot V) I ]μ是动力粘性系数λ是第二粘度在斯托克斯假设下取-2μ/3。对于不可压缩流体∇·V0第二项直接消失式子大幅简化[ \tau \mu \left( \nabla V (\nabla V)^T \right) ]3.3 化简的关键手笔用连续性方程收拾左边从柯西方程到N-S方程中间这一步是很多人“背公式”时完全忽略的但它才是理解N-S方程“形状”的钥匙。柯西方程的左边是∂(ρV)/∂t ∇·(ρV⊗V)。直接用链式法则展开第i个分量左边 ρ∂V_i/∂t V_i∂ρ/∂t V_i∂(ρV_j)/∂x_j ρV_j∂V_i/∂x_j组合一下把带∂V_i/∂t和∂V_i/∂x_j的项挑出来剩下的项放到一起左边 ρ(∂V_i/∂t V_j∂V_i/∂x_j) V_i[∂ρ/∂t ∂(ρV_j)/∂x_j]看仔细方括号里是什么那不是连续性方程的左边吗质量守恒告诉我们它等于零。于是整个“多余部分”干干净净地消失左边 ρ(∂V_i/∂t V_j∂V_i/∂x_j)写回矢量形式就是[ \rho \left( \frac{\partial V}{\partial t} (V \cdot \nabla) V \right) ]这一步的深刻含义是质量守恒律“嵌入”了动量方程。所以N-S方程左边那个看似很特别的括号∂V/∂t (V·∇)V其实是物质加速度跟随流体粒子一起运动时感受到的加速度而不是某个固定观测点的变化率。这就是为什么它是Fma里的那个“a”。3.4 加上牛顿流体本构N-S方程诞生不可压缩牛顿流体的本构关系代回柯西方程∇·σ ∇·(-pI τ)其中-∇p项直接出来。τ的散度呢[ \frac{\partial \tau_{ij}}{\partial x_j} \mu \frac{\partial^2 u_i}{\partial x_j \partial x_j} \mu \frac{\partial}{\partial x_i} \left( \frac{\partial u_j}{\partial x_j} \right) ]第一个交叉项正是μ∇²V。第二个括号里是∇·V不可压缩条件下为零。所以[ \rho \left( \frac{\partial V}{\partial t} (V \cdot \nabla) V \right) - \nabla p \mu \nabla^2 V \rho g ]这就是大名鼎鼎的N-S方程。整个过程走下来每个符号都有明确来历没有一个是需要“硬背”的。从输运定理出发加上牛顿第二定律、应力分解、牛顿本构、连续性方程一气呵成。3.5 N-S方程每一项的物理直觉把方程当成Fma读每一项都有对应的受力来源项名称物理含义生活类比ρ(∂V/∂t)局部加速度固定位置上速度随时间的变化站在桥上看同一位置水流忽快忽慢ρ(V·∇)V对流加速度粒子搬运到新位置引起的速度变化从平静河湾漂到急流中央被加速-∇p压力梯度力高压指向低压的力风从高压区吹向低压区μ∇²V粘性力速度分布弯曲导致的耗散力搅动糖浆时周围的糖浆被带动ρg体积力外部力场重力水往低处流这套表格是我在课堂上讲N-S时的保留节目。它让方程不是一个符号堆砌而是一张“力的全家福”。下次看到方程哪怕不看推导你也应该能在脑子里把每一项“翻译”成一股力、一种流动现象。4. 推导中的翻车现场与排查技巧4.1 通量项符号翻车流出到底加还是减这是新手翻车率最高的一处。问题几乎都出在同一件事上既在公式里写了面积分又在脑子里做了一次“流出加、流入减”的判断结果符号被处理了两遍。正确的做法是面积分里统一写∫_CS ρb(V·n)dA前面别再加额外的正负号。V·n为正时流出逆然为正贡献为负时流入自动贡献为负。你先入为主地“手动”减流入就会在入口上搞出双重减号。我的习惯是凡涉及方向判断立即在草稿纸上画出控制体左边一个入口、右边一个出口标出外法线n的箭头。符号判断只依靠V·n的数值正负不做任何直觉正负。这个小习惯救过我无数次。4.2 把(V·∇)和∇·V搞混看起来像物理上是两个东西(V·∇)V和V(∇·V)是完全不同的两项。前者是“速度场作为微分算子去作用于自己的速度场”是N-S方程中典型非线性项后者是“速度乘以速度场的散度”只在密度变化、体积胀缩时才出现。不幸的是它们在张量运算中长得很像缩写得不好甚至分不清。一个可靠的辨认方法V·∇是一个标量算子作用在矢量上仍得到矢量∇·V则是一个标量。你在N-S方程左边看到的、产生非线性效应的那个是(V·∇)V它对应流体粒子经过空间不均匀流场时产生的“搬运加速度”不是体积膨胀率。考试和编程里把这两项写反方程物理性质会彻底变坏。4.3 高斯定理的隐藏前提间断面会让你翻车高斯定理要求被积函数在区域内连续可微。当流场出现激波、自由面等间断结构时微分形式的N-S方程在间断面上不再成立只有积分形式输运定理仍然有效。这就是为什么CFD里所有主流求解器比如有限体积法都在积分形式下做离散而不是硬刚微分方程。实际工程提示跨音速流、气液两相流这类含强间断的场景如果你在校核方程时非要用微分形式去理解必须额外引入Rankine-Hugoniot跳变条件。而用积分形式写瞬时尚则完全不用操心间断。这也是我在学了CFD后才真正体会到“输运定理的积分形式才是流体力学的底层原生形态”。4.4 不可压缩假设什么时候可以偷懒、什么时候不能不可压处理后的N-S方程精简很多密度常数、∇·V0、粘性项退化为μ∇²V。但“不可压缩”不等于“密度不变”那么简单。它的严格含义是流动中压强、温度变化引起的密度相对变化可以忽略。低速液体流动、马赫数低于0.3的气体流动放心用不可压形式。但高速气流、燃烧、激波、高超声速流动密度剧烈变化必须用完整可压缩N-S方程并且要额外补一个能量方程和状态方程来闭合。初学者最典型的问题是把水和低速空气当成永不可压却在跨声速喷管问题上套不可压模型结果数值噪音暴涨物理结果错得离谱。4.5 总在漏掉边界条件即使完整推导出了N-S方程你还需要边界条件才能真正求解。无滑移壁面速度与壁面一致、入口速度/压力条件、出口通量条件、自由表面条件——每一样都直接影响解的形态。我见过不少学完推导的人以为有了方程就有了一切结果跑到仿真软件里眼巴巴盯着残差看却不知边界条件已经给了物理上荒谬的设定。记住方程描述规律边界条件定义场景两者合起来才是完整问题。这个意识比方程本身早点建立更好。5. 把输运定理装进自己的知识体系学习建议与心得5.1 以输运定理为核心建立“守恒律控制表”到了这一步掌握最有力的方式是把之前所有的B、b代换整理成一个控制表守恒定律广延量B强度量b方程形式质量m1∂ρ/∂t∇·(ρV)0动量mVVρ(∂V/∂t(V·∇)V)∇·σρg能量m(e|V|²/2)e|V|²/2∂E/∂t∇·((Ep)V)...这等于告诉你整个流体力学理论树其实是从一个输运定理“长”出来的。质量、动量、能量守恒并不是三个不相干方程而是同一个数学模板在不同物理量上投影的结果。有了这张表你连“为什么能量方程那么复杂”都能猜到大概——因为b里带了内能、动能还做功传热。5.2 学习路径建议先积分、后微分、再上CFD很多人一上课就直奔微分形式N-S方程这是最大的误解。正确顺序是先在积分形式下理解输运定理理解“系统变化率 控制体积累 净流出”这句会计平衡账学会把积分形式转微分形式高斯定理 局部化引理再给特定守恒律配b和源项最后才谈CFD离散理解有限体积法为什么在每条网格边界计算通量。我做CFD项目时有一个体会有限体积法程序里的每一个通量计算单元本质上都在帮你“手算”输运定理中那个∫_CS ρb(V·n)dA项。如果你搞懂了输运定理看CFD代码会特别亲切——完全不像看天书。5.3 一点个人经验推导两遍胜过默写十遍最后分享一个我带人时反复验证过的经验。凡是能把输运定理自己在草稿纸上完整推两遍——第一遍照着推第二遍合上笔记推——的人之后看N-S方程就像看Fma一样自然。那些总想跳过推导直接做题的人通常在各类测验、面试、项目里反复栽跟头因为稍微变形一下题目就识别不出来。具体操作建议第一遍照着笔记或本文流程推理解每一步第二遍合上笔记从“设系统与控制体在t时刻重合”开始独立推完输运定理再推连续性方程和柯西方程第三遍可选加入牛顿本构推到不可压N-S方程。三遍下来你对N-S方程的“熟悉感”会从“背出来了”变成“本该如此”。真到用的时候——写个脚本、调个仿真器、分析实验数据——你会发现自己完全不需要去“回忆”方程因为整个推导图景已经在脑子里长成了本能。顺带说一句别小看控制体的选取技巧。做工业计算时控制体选得好不好直接影响方程的简洁程度。压气机内部流动选随体控制体管道流量计算选固定控制体两相流界面问题选变形控制体——这些都是课本不细讲、但实际项目里天天要用的经验。等你把输运定理真正用熟了这些“怪招”自然会从直觉里长出来。