第3章物理现象的数学描述（Mathematical Description of Physical Phenomena）

作者

本章由 F. Moukalled、L. Mangani 和 M. Darwish 合著。三位作者均有流体力学与传热学的双重学术背景，本章对 N-S 方程组的推导不仅给出标准形式，还强调非牛顿流体、可压缩 / 不可压缩流的简化、量纲分析与无量纲数（Re、Pr、Ma、Pe 等），这些都是工业 CFD 应用中的关键决策点。本章是 Ch 5（积分形式守恒方程）的基础。

内容概述

本章是全书 20 章的"物理理论核心" —— 把流体力学三大守恒定律（质量、动量、能量）从基本原理出发，系统地推导出 N-S 方程组的控制微分方程形式；最后给出一个适用于任意标量 / 向量 / 张量物理量的一般性守恒方程，作为 Ch 5-14 数值推导的统一起点。

本章逻辑主干： 1. §3.1 引言 —— N-S 方程组的历史地位与广泛应用。 2. §3.2 流体流动的分类 —— 牛顿 / 非牛顿流体、可压 / 不可压、稳态 / 非稳态、层流 / 湍流、有旋 / 无旋等；从数学角度又可分为双曲型（超声速、可压缩、瞬态）、抛物型（边界层）、椭圆型（亚声速、不可压、回流），这一分类决定了求解方法的根本差异。 3. §3.3 Eulerian 与 Lagrangian 描述 —— 给出 Lagrangian 物质导数 $D\phi/Dt = \partial \phi/\partial t + \mathbf{v} \cdot \nabla \phi$ (Eq. 3.2)；引入 Reynolds 输运定理 (RTT) (Eq. 3.4) 作为 Lagrangian 物质积分到 Eulerian 控制体积分的桥梁。 4. §3.4 质量守恒 —— 连续性方程 $\partial \rho/\partial t + \nabla \cdot (\rho \mathbf{v}) = 0$ (Eq. 3.22)，可压缩 / 不可压缩简化。 5. §3.5 动量守恒 —— 应力张量 $\boldsymbol{\sigma}$ 的分解 (Eq. 3.40)、Stokes 假设、Newton 流体动量方程 (Eq. 3.49)、N-S 方程的最终形式 (Eq. 3.51)。 6. §3.6 能量守恒 —— 总能量方程 (Eq. 3.61)、Fourier 热传导定律 (Eq. 3.70)、温度形式的能量方程 (Eq. 3.78)，及对不可压流、固体、理想气体的简化 (Eq. 3.79-3.82)。 7. §3.7 一般守恒方程 —— 任意标量 $\phi$ 的统一控制方程 (Eq. 3.93)，其形式为瞬态项 + 对流项 = 扩散项 + 源项。 8. §3.8 量纲分析与无量纲数 —— Reynolds 数 $Re = \rho U L / \mu$、Prandtl 数 $Pr = \mu c_p / k$、Mach 数 $Ma = U/c$、Peclet 数 $Pe = Re \cdot Pr$ 等。 9. §3.9 数学性质 —— 椭圆型、抛物型、双曲型 PDE 的判别方法（系数矩阵特征值判别）。 10. §3.10 闭包。

核心方程与概念

本章的核心方程分布密度大，按"母方程 → 守恒定律 → 简化形式 → 一般化"四层列举。

一、Lagrangian 物质导数（§3.3）

总导数（物质导数 / 实质导数） (Eq. 3.2)

$$\frac{D\phi}{Dt} = \frac{\partial \phi}{\partial t} + \mathbf{v} \cdot \nabla \phi \tag{3.2}$$

其中 $\mathbf{v} = (u, v, w)$。这是从 Lagrangian 描述到 Eulerian 描述的核心桥梁。Ch 13 时间离散、Ch 11 对流离散都隐式用到这一分解。

加速度 (Eq. 3.3)

$$\frac{D\mathbf{v}}{Dt} = \frac{\partial \mathbf{v}}{\partial t} + (\mathbf{v} \cdot \nabla) \mathbf{v} \tag{3.3}$$

注意 $(\mathbf{v} \cdot \nabla) \mathbf{v}$ 是向量，其分量为 $\sum_j v_j \partial v_i / \partial x_j$。N-S 方程的"非线性"完全源自这一项。

Reynolds 输运定理 (RTT) (Eq. 3.4)

$$\frac{dB}{dt}\bigg|_{MV} = \frac{d}{dt} \int_{V(t)} b \rho \, dV + \int_{S(t)} b \rho \mathbf{v}_r \cdot \mathbf{n} \, dS \tag{3.4}$$

其中 $b = dB/dm$ 是 $B$ 的强度量，$\mathbf{v}_r = \mathbf{v} - \mathbf{v}_s$ 是流体相对于控制体表面的相对速度。这是 Ch 5 把 Lagrangian 守恒定律转化为 Eulerian 控制方程的数学基础。

二、质量守恒（§3.4）

连续性方程（一般形式） (Eq. 3.22)

$$\frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{v}) = 0 \tag{3.22}$$

在不可压情形，$\rho = \text{const}$，故

$$\nabla \cdot \mathbf{v} = 0 \tag{3.23}$$

Ch 15（不可压缩流体算法）的核心就是"无显式压力方程的 N-S 求解"。

三、动量守恒（§3.5）

应力张量分解 (Eq. 3.40)

$$\boldsymbol{\sigma} = -p \mathbf{I} + \boldsymbol{\tau} \tag{3.40}$$

其中 $p$ 是热力学压力，$\mathbf{I}$ 是单位张量，$\boldsymbol{\tau}$ 是黏性应力张量。

Stokes 假设 —— 对各向同性牛顿流体

$$\boldsymbol{\tau} = \mu (\nabla \mathbf{v} + \nabla \mathbf{v}^T) - \frac{2}{3} \mu (\nabla \cdot \mathbf{v}) \mathbf{I} \tag{3.45}$$

第二项 $-\frac{2}{3} \mu (\nabla \cdot \mathbf{v}) \mathbf{I}$ 是 Stokes 体积黏性项；在不可压情形 $\nabla \cdot \mathbf{v} = 0$ 时该项消失，化为最常见的对称形式

$$\boldsymbol{\tau} = \mu \left( \nabla \mathbf{v} + \nabla \mathbf{v}^T \right) \tag{3.46}$$

N-S 动量方程（不可压牛顿流体） (Eq. 3.51)

$$\frac{\partial (\rho \mathbf{v})}{\partial t} + \nabla \cdot (\rho \mathbf{v} \mathbf{v}) = -\nabla p + \nabla \cdot (\mu \nabla \mathbf{v}) + \rho \mathbf{g} \tag{3.51}$$

黏性项可显式展开为

$$\nabla \cdot (\mu \nabla \mathbf{v}) = \mu \nabla^2 \mathbf{v} + \nabla \mu \cdot \nabla \mathbf{v}$$

对常黏度情形简化为 $\mu \nabla^2 \mathbf{v}$。

可压缩 N-S 动量方程 (Eq. 3.49)

$$\frac{\partial (\rho \mathbf{v})}{\partial t} + \nabla \cdot (\rho \mathbf{v} \mathbf{v}) = -\nabla p + \nabla \cdot \boldsymbol{\tau} + \rho \mathbf{g}$$

此时黏性项需保留 Stokes 体积黏性项。

四、能量守恒（§3.6）

总能量方程（一般形式） (Eq. 3.61)

$$\frac{\partial (\rho \hat{h})}{\partial t} + \nabla \cdot (\rho \mathbf{v} \hat{h}) = -\nabla \cdot \mathbf{q}_s + \frac{Dp}{Dt} + \boldsymbol{\tau} : \nabla \mathbf{v} + \dot{q}_V \tag{3.61}$$

其中 $\hat{h}$ 是比焓，$\mathbf{q}_s$ 是热通量向量，$\dot{q}_V$ 是体积热源，$\boldsymbol{\tau} : \nabla \mathbf{v}$ 是黏性耗散函数 $\Phi$。

Fourier 热传导定律 (Eq. 3.70)

$$\mathbf{q}_s = -k \nabla T \tag{3.70}$$

其中 $k$ 是热导率。

温度形式的能量方程（不可压流） (Eq. 3.79)

$$\frac{\partial (\rho c_p T)}{\partial t} + \nabla \cdot (\rho c_p \mathbf{v} T) = \nabla \cdot (k \nabla T) + \dot{q}_V + \rho T \frac{Dc_p}{Dt} \tag{3.79}$$

对常物性情形（$c_p, k$ 常数），简化为

$$\frac{\partial T}{\partial t} + \mathbf{v} \cdot \nabla T = \alpha \nabla^2 T + \frac{\dot{q}_V}{\rho c_p}$$

其中 $\alpha = k / (\rho c_p)$ 是热扩散率。该形式与一般守恒方程 (Eq. 3.93) 完美匹配。

五、一般守恒方程（§3.7） —— FVM 各章的母方程

通用 $\phi$ 守恒方程 (Eq. 3.93)

$$\underbrace{\frac{\partial (\rho \phi)}{\partial t}}_{\text{瞬态项}} + \underbrace{\nabla \cdot (\rho \mathbf{v} \phi)}_{\text{对流项}} = \underbrace{\nabla \cdot (\Gamma \nabla \phi)}_{\text{扩散项}} + \underbrace{S_\phi}_{\text{源项}} \tag{3.93}$$

当 $\phi = 1$ 时回到连续性方程；$\phi = \mathbf{v}$ 时回到 N-S 动量方程；$\phi = c_p T$ 时回到能量方程。这就是 Ch 5-14 反复使用的"母方程"。

六、量纲分析与无量纲数（§3.8）

Reynolds 数 (Eq. 3.94)

$$Re = \frac{\rho U L}{\mu} = \frac{U L}{\nu} \tag{3.94}$$

表征惯性力与黏性力之比。$Re \ll 1$ 为 Stokes 流（爬流），$Re \gg 1$ 为惯性主导。

Prandtl 数

$$Pr = \frac{\mu c_p}{k} = \frac{\nu}{\alpha} \tag{3.96}$$

表征动量扩散率与热扩散率之比。$Pr \ll 1$（液态金属）热扩散远快于动量扩散；$Pr \gg 1$（油脂）反之。

Mach 数

$$Ma = \frac{U}{c} = \frac{U}{\sqrt{\gamma R T}} \tag{3.97}$$

表征流速与当地声速之比。$Ma < 0.3$ 通常视为不可压。

Peclet 数

$$Pe = Re \cdot Pr = \frac{UL}{\alpha} \tag{3.98}$$

表征对流热输运与扩散热输运之比，是 Ch 11-12 高分辨率对流格式选择的关键判据（$Pe$ 大 → 数值格式需要捕捉对流方向性）。

Reynolds 数与湍流的关系 —— $Re > Re_{crit}$ 时层流失稳，转变为湍流。临界 $Re$ 因几何不同（圆管、内流、外流）而异。

七、PDE 数学分类（§3.9）

2D 二阶 PDE 一般形式 $A \phi_{xx} + 2B \phi_{xy} + C \phi_{yy} + \cdots = 0$
判别式 $\Delta = B^2 - AC$：$\Delta > 0$ 双曲型，$\Delta = 0$ 抛物型，$\Delta < 0$ 椭圆型
流体力学对应：
不可压 N-S（压力方程）：椭圆型（信息全向传播）
边界层：抛物型（仅顺流方向信息传播）
超声速可压缩：双曲型（沿特征线传播，支持激波间断）

关键结论

Lagrangian 物质导数 $D\phi/Dt = \partial \phi/\partial t + \mathbf{v} \cdot \nabla \phi$ 是 Euler 描述下 N-S 方程非线性的根源：$(\mathbf{v} \cdot \nabla) \mathbf{v}$ 项的存在导致 N-S 方程组无法用线性叠加原理，是湍流、激波、涡旋等复杂流动物理的数学源头。
Reynolds 输运定理 (RTT) 是 Lagrangian 守恒定律到 Eulerian 控制方程的数学桥梁：本书选 Euler 描述为基础（Ch 3.3.1 末段明确给出理由），RTT 把物质积分 $dB/dt|_{MV}$ 转化为控制体积分 + 界面对流通量 (Eq. 3.4)。
N-S 方程组 = 质量 + 动量 + 能量三大守恒定律的 PDE 化：通过 RTT 展开，得到 5 个未知量（$\rho, \mathbf{v}, T$ 或 $\hat{h}$）的 5 个 PDE 方程组（连续性 + 3 个动量 + 能量）；加上状态方程 $p = \rho R T$ 形成闭包。Ch 15-16 的 SIMPLE 系列算法即针对该方程组的数值求解。
一般守恒方程 (Eq. 3.93) 是 Ch 5-14 数值推导的统一起点：$\phi$ 任意（标量、向量分量、张量分量），各项分别在不同章节单独离散化：瞬态项 Ch 13、对流项 Ch 11-12、扩散项 Ch 8、源项 Ch 14。这种"分项离散化"是 FVM 相对 FEM 的优势 —— 可以针对每一项的数学性质选择最合适的数值方法。
量纲分析把"具体工况"抽象为"无量纲参数"：$Re$ 决定流态（层流 / 湍流）、$Ma$ 决定可压性、$Pr$ 决定热-动量耦合方式、$Pe$ 决定对流格式要求。Ch 17 湍流建模、Ch 11-12 对流格式选择都直接依赖这些无量纲数。
PDE 数学分类（椭圆 / 抛物 / 双曲）决定了求解方法的根本差异：椭圆型需全场联立求解（SIMPLE 类的迭代方法），抛物型可顺流推进（边界层方程），双曲型需捕捉特征方向（高分辨率格式、激波捕捉）。本书 Ch 11-12 的 TVD/NVD 格式正是双曲型守恒律的现代数值方法。

挑战和开放性问题

湍流的闭包问题 —— 本章推导的 N-S 方程 (Eq. 3.51) 是层流形式。湍流情形下，$\mathbf{v} \to \bar{\mathbf{v}} + \mathbf{v}'$ 分解后会出现 Reynolds 应力项 $-\rho \overline{\mathbf{v}' \mathbf{v}'}$，这是一个未知的二阶张量，需要额外的湍流模型（如 k-ε, k-ω, RSM）来闭包。Ch 17 会专门处理。本章未深入讨论。
多相流与多组分流的推广 —— 本章只处理单相单组分流体。多相流（VOF、Euler-Euler、Level-Set）和多组分（组分输运方程、化学反应）的控制方程集需要扩展。本书未深入（仅在 §3.7 提到"标量 $\phi$ 可以是化学组分的质量分数"）。
非牛顿流体的本构方程 —— 本章 Stokes 假设 (Eq. 3.45) 仅适用于牛顿流体。血液、聚合物熔体、泥浆等非牛顿流体需要更一般的本构方程（如 Carreau、Cross、Bingham 模型），FVM 离散化中如何处理非牛顿黏度是开放话题。
可压缩流的强间断（激波）捕捉 —— 本章 PDE 分类提到双曲型支持间断，但具体到激波附近的数值稳定性（Carbuncle 现象、HLLC 格式、Riemann 求解器）需要更深入的数学工具，本书未深入展开。
热力学与流体力学的耦合 —— 本章对状态方程、黏性-温度关系（如 Sutherland 公式）、湍流 Prandtl 数变效应等热力学细节展开不够。Ch 16 的可压缩流算法会处理一部分。
非结构网格上的张量运算 —— 本章所有张量运算在 Cartesian 坐标系下展开；在非结构网格（Ch 6）上，应力张量 $\boldsymbol{\sigma}$ 的分量计算需要局部坐标系变换，是 FVM 实现的难点之一，本书未深入展开。

个人反思与批判性分析

本章是全书最"硬核"的一章，70+ 页、100+ 公式、几十个例题。完整读下来需要读者具备本科级工程流体力学（White、Fox-McDonald、Pritchard 级别）+ 高等工程数学（向量演算、张量分析）+ 偏微分方程（椭圆 / 抛物 / 双曲型理论）的基础。

从写作特点看，本章有以下几个值得借鉴之处：

"宏观 → 数学" 的逻辑链条清晰 —— §3.1 引言（应用背景）→ §3.2 分类（物理直觉）→ §3.3 Euler/Lagrangian + RTT（数学基础）→ §3.4-3.6 三大守恒定律（核心推导）→ §3.7 一般化（统一视角）→ §3.8 量纲分析（工程抽象）→ §3.9 PDE 分类（求解方法选择）。这种"先物理再数学后方法"的递进结构非常符合工程读者的学习曲线。
大量例题演示 —— §3.5 中的圆柱绕流 Reynolds 数计算、§3.5 中的边界层厚度估算、§3.6 中的温度场算例等。每个例题都把"推导的方程"应用到"具体物理场景"，帮助读者建立"方程 ↔ 物理"的映射。
应力张量分解的 Stokes 假设详细展开 —— 作者没有简单给出"Newton 流体 $\boldsymbol{\tau} = \mu (\nabla \mathbf{v} + \nabla \mathbf{v}^T)$"，而是从"应力张量的对称性"出发推导，给出体积黏性项 $-\frac{2}{3} \mu (\nabla \cdot \mathbf{v}) \mathbf{I}$ 的来源，再说明在不可压情形该项消失。这种"严格推导 + 简化说明"的双层结构对培养读者的"工程近似思维"非常有价值。

但本章也存在以下不足：

对非牛顿流体的处理过于简略 —— 仅在 §3.2 提到"非牛顿流体剪切应力-剪切速率非线性"，没有给出非牛顿本构方程（Carreau 模型、Cross 模型、Bingham 模型）的任何具体形式。对于工业 CFD 应用（如高分子加工、食品工业、生物流变学），非牛顿本构是必备知识。
量纲分析仅给出"经典 4 个无量纲数"，缺少现代无量纲分析 —— 例如磁流体力学中的 Magnetic Reynolds 数 $Re_m$、Alfvén 数 $Al$、磁 Prandtl 数 $Pr_m$；多相流中的 Weber 数 $We$、毛细管数 $Ca$；燃烧中的 Damköhler 数 $Da$、Karlovitz 数 $Ka$。本书的量纲分析章节显然面向"传统单相不可压流"，对现代工业应用覆盖不全。
RTT 的可压缩性处理 —— 本章的 RTT 推导 (Eq. 3.4) 是"密度 + 强度量 $b$"形式，但在可压缩流中"控制体固定 vs 动"的差异更微妙。作者在 Eq. 3.4 用了相对速度 $\mathbf{v}_r = \mathbf{v} - \mathbf{v}_s$，但没有给出"控制体表面变形"对界面通量积分的具体推导（这在 Ch 13 动网格中会涉及）。对 ALE 网格 / 动边界感兴趣的读者需要补充额外资料。
PDE 分类 (3.9 节) 与 N-S 方程组的"混合型"问题 —— 作者正确指出 N-S 方程组是"混合型"（亚声速段椭圆型 + 超声速段双曲型），但对"跨声速"流动的数值方法（人工黏性、Preconditioning、Entropy 稳定格式）未展开。这一领域在 2010 年后是研究热点（Ch 16 可压缩流部分会触及，但深度有限）。

总体而言，本章的"工程流体力学"色彩浓郁，写作风格"扎实、完整、保守"，与 Ferziger-Peric 的 Computational Methods for Fluid Dynamics（更偏数值理论）、Anderson 的 CFD: The Basics with Applications（更偏科普）形成三角对标。Moukalled 的优势在于"既有理论推导又有工程简化"，对工程类学生与从业者最友好。

重要参考文献

[X1] Currie I.G. (1974) Fundamental Mechanics of Fluids. McGraw-Hill. （N-S 方程推导的经典参考）
[X2] White F.M. (2016) Fluid Mechanics, 8th Edition. McGraw-Hill. （本科流体力学标准教材，与本章互补）
[X3] Batchelor G.K. (1967) An Introduction to Fluid Dynamics. Cambridge University Press. （N-S 方程组的高级推导参考）
[X4] Patankar S.V. (1980) Numerical Heat Transfer and Fluid Flow. Hemisphere. （FVM 与 SIMPLE 算法原始文献）
[X5] Bird R.B., Stewart W.E., Lightfoot E.N. (2002) Transport Phenomena, 2nd Edition. Wiley. （RTT 与守恒定律推导的经典参考）
[X6] Slattery J.C. (1999) Advanced Transport Phenomena. Cambridge University Press. （RTT 严格推导参考）
[X7] Tritton D.J. (1988) Physical Fluid Dynamics, 2nd Edition. Oxford University Press. （PDE 分类与流体力学对应）
[X8] Anderson J.D. (2002) Modern Compressible Flow: With Historical Perspective, 3rd Edition. McGraw-Hill. （可压缩流与无量纲数分析）
[X9] Pope S.B. (2000) Turbulent Flows. Cambridge University Press. （湍流闭包问题与 RANS 推导）
[X10] Spalding D.B. (1972) A Novel Finite-Difference Formulation for Differential Expressions Involving Both First and Second Derivatives. Int. J. Num. Methods Eng., 4: 551-559. （FVM 离散化的开创性文献）

第3章 物理现象的数学描述（Mathematical Description of Physical Phenomena）

作者