流体力学：从概念到计算

FLUIDDYNAMICS流体力学：从概念到计算从物理图像到数值求解的完整链路2026课程导览01流体与流动建立流体的物理图像02控制方程从守恒定律到N-S方程03无量纲与相似用无量纲数刻画流动04CFD数值计算从方程到离散求解SECTION流体与流动认识流体，是一切计算的起点01流体是什么流体在剪切力作用下持续变形，无法保持固定形状。理解基本属性，是建立一切方程的前提。核心性质5项密度单位体积的质量黏度抵抗相对滑动的能力表面张力液体表面收缩趋势压缩性体积随压力变化膨胀性体积随温度变化牛顿流体vs非牛顿流体按剪应力与速度梯度的关系划分：牛顿流体：满足

τ=η·γ̇

线性关系，如空气与水非牛顿流体：不遵循线性规律，如血液、蜂蜜、油漆研究对象液体与气体与固体不同，流体在剪切力作用下持续变形，无法保持固定形状——任意小的切向力都会引发不断形变，这是流体最本质的特征。静止流体的压强规律流体力学分静力学与动力学两大分支，先看平衡态。静压强分布规律水平方向处处相等同一水平面上各点压强一致。竖直方向随深度线性增加压强随深度成正比例变化。任一点压强只与深度和密度有关即p=ρgh，是静力学的基本结论。静止流体阿基米德定律浮力等于排开流体的重量F_b=ρgV，与物体浸没深度无关。解释船舶漂浮排水量等于船体重量时，船舶漂浮于水面。解释人体水中受力人体排开的水量产生向上的支撑浮力。浮力壁面作用力通过压强在平面或曲面上的积分求得。对平面壁：F=∫pdA；对曲面壁需分别求水平与竖直分力后进行矢量合成。积分求解流体的运动如何描述拉格朗日追踪质点，欧拉观察固定点。把物理图像翻译成方程，先选对观察的视角。工程首选欧拉法：直接给出流场中速度、压强的空间分布，满足绝大多数工程问题需求。两种视角，两种描述LagrangianvsEulerian拉格朗日法：追踪单个流体质点轨迹，记录每个时刻的位置与速度，如跟踪随水流漂动的微粒。欧拉法：固定空间点，观察流体流过时物理量随时间的变化，如在河道架设测站记录流速。配套概念：流线与流管Streamline&StreamTube流线：某时刻处处与速度矢量相切的曲线，直观呈现瞬时流场结构。流管：由流线围成，管内外流体互不穿越。SECTION控制方程守恒定律，是流体世界的底层规则02三大守恒定律撑起方程体系守恒型方程基于固定控制体导出，允许参数间断，可处理激波等非连续流动。流体运动的数学描述建立在三条不可违背的物理定律之上：质量守恒、动量守恒、能量守恒。质量守恒→连续性方程微分形式∂ρ/∂t+∇·(ρv)=0，即密度变化率与质量通量散度相平衡不可压缩流体（密度恒定）：简化为

∇·v=0，流入等于流出动量守恒·能量守恒动量守恒体现牛顿第二定律能量守恒体现热力学第一定律封闭条件本构方程应力与变形速率的关系状态方程如完全气体

p=ρRT纳维-斯托克斯方程描述黏性牛顿流体运动的偏微分方程，1827年纳维提出、1845年斯托克斯完善，被视为流体力学的核心。推导基础基于牛顿第二定律将压强、黏性应力与速度场联系起来结构拆解4项对流加速度项

u·∇u—非线性来源压力梯度项

-∇p/ρ—驱动流动黏性项

ν∇²u—耗散能量外力项（如重力）—外部作用非线性意义对流项的非线性意味着即便边界条件简单，流动也可能演化出复杂行为，这正是湍流的数学根源。未解难题三维N-S方程解的存在性与光滑性是克雷数学研究所千禧年七大难题之一，悬赏一百万美元，至今未被正式认证为已解决。从守恒方程到伯努利方程“简化不是妥协，而是让方程走向工程的第一步。”—由守恒方程到伯努利方程伯努利方程能量守恒的直接推论理想流体沿流线定常流动时，能量守恒直接导出

½ρv²+ρgh+p=常量。流速越快处压强越低——这正是机翼升力与喷雾器抽液的原理。同一流线动压+势压+静压泊肃叶公式层流流量的四次方定律水平圆管层流的速度分布与流量

Q=πR⁴Δp/(8ηl)，流量与管半径的四次方成正比。血管或管道内径的微小变化，会带来流量的显著改变。∝R⁴内径敏感SECTION无量纲与相似无量纲数，是流动的指纹03雷诺数：流动状态的判据“一个无量纲数，判定两种流态”“雷诺数相等，流动即相似”定义Re=ρvL/μ：惯性力与黏性力之比两种流态层流黏性力占优→流动平稳有序湍流惯性力占优→微小扰动被放大，流动失稳管道流动的判据Re<约

2000层流Re>约

4000湍流其间过渡区更根本的作用——相似准则几何相似的流动，雷诺数相等则无量纲速度场相同。据此可用缩比模型做风洞或水槽实验，将结果外推到真实尺寸原型，以低成本获得可靠设计依据。马赫数与弗劳德数雷诺数、马赫数、弗劳德数，三者共同构成风洞实验相似律的基础，让模型实验准确预测实物表现。流速接近声速时，黏性让位于压缩效应。三者共同构成风洞实验相似律的基础，让模型实验准确预测实物表现。马赫数：压缩效应的标尺定义Ma=v/c，流速与当地声速之比，无量纲而非速度单位。Ma<0.8亚音速0.8–1.2跨音速>1.2超音速与高超音速应用决定可压缩程度与激波形成，飞机外形设计与发动机工作状态均围绕它展开。弗劳德数：自由液面的判据定义Fr=v/√(gL)，惯性力与重力之比。物理意义用于船舶兴波、水上飞机滑行等自由液面问题。应用惯性力与重力之比，用于自由液面问题。相似律基础风洞实验SECTIONCFD数值计算把方程交给计算机04CFD：把方程交给计算机N-S方程解析解不可得，CFD以数值方法离散求解，揭示流场细节——抽象的偏微分方程，由此变成工程师可用的设计答案。1前处理几何建模、区域离散与网格生成2流场计算按选定模型与格式迭代求解3后处理结果可视化，流场分析与诠释抽象的偏微分方程，由此变成工程师可用的设计答案。解析无门，数值求解解析解不可得，CFD以数值方法离散求解，揭示流场细节第三种研究手段始于20世纪30年代初，与理论分析、实验并列看得更深揭示实验难以观测的内部流场与机理改得更快参数调整快速，适合设计优化迭代省得更多大幅节省风洞试验的人力物力有限体积法：工程的主流选择一套框架，多种物理，天然守恒三种离散方法，有限体积法凭借守恒性成为工程CFD的主流选择三种离散方法离散化路径有限差分法直接在网格节点近似导数，最直观有限元法擅长复杂几何，多用于结构分析有限体积法以守恒性见长，是工程CFD的主流核心思路积分离散将计算域划分为互不重叠的控制体积，在每个控制体积上对守恒方程做积分。把偏微分方程转化为节点未知量的代数方程组，逐点求解。局部守恒⇒全局守恒统一输运形式通用方程对流输运ρuφ+扩散输运−Γ∇φ=源项S对流、扩散与源项三者平衡。φ=1连续性方程φ=u动量方程φ=T能量方程φ取1、速度u、温度T时，同一方程分别退化为上述方程。离散格式与湍流的权衡“格式选择：精度与稳定性不可兼得”对格式与模型的权衡，正是CFD工程实践的核心判断力。格式选择空间离散格式中心差分二阶精度，但对流占优时（网格贝克莱数>2）系数可能为负，引发非物理震荡甚至不收敛一阶迎风强制用上游信息，系数非负、解稳定有界，但数值扩散会抹平陡峭梯度二阶迎风利用上游两个节点线性外推，兼顾精度与稳定性