跳转至

第2章 向量演算回顾（Review of Vector Calculus）

作者

本章由 F. Moukalled、L. Mangani 和 M. Darwish 合著，作者署名排序与全书一致。三位作者均为工程流体力学背景的 CFD 学者，并非纯数学家；本章的"向量演算回顾"也明显偏向工程流体力学需要的子集 —— 重点是散度定理（Gauss）、Green 定理、Stokes 定理、Leibniz 积分法则这一组对 N-S 方程组推导与 FVM 离散化至关重要的定理。本章假定读者已具备线性代数与多元微积分基础，因此对证明过程高度压缩，对几何直观与符号约定则给予充分展开。

内容概述

本章是 Ch 3（物理守恒方程推导）的数学预备章节。Ch 3 将基于本章建立的向量与张量记号，从质量、动量、能量三大守恒定律出发推导出 N-S 方程组；Ch 5 又会基于本章的散度定理（式 2.84）将 N-S 方程组写成积分形式，进而展开 FVM 离散化。

本章的逻辑主干可以概括为： 1. §2.1 标量 / 向量 / 张量 / 矩阵的记号约定 —— 标量用斜体、向量用粗体小写罗马字母、张量用粗体希腊字母、矩阵用粗体大写罗马字母；乘积结果用 (scalar)、[vector]、{tensor} 三种括号加以视觉区分。 2. §2.2 向量运算 —— 加法、数乘、点积、叉积、向量投影、单位向量等。 3. §2.3 - §2.4 张量运算 —— 二阶张量的分量表示、转置、迹、对称 / 反对称分解、点积与双点积。 4. §2.5 五大基本定理 —— 梯度定理（线积分）、Green 定理（2D 面-线积分）、Stokes 定理（线-面积分）、散度定理（体-面积分）、Leibniz 积分规则（动控制体积分）。这五个定理是后续守恒定律推导与 FVM 通量展开的"母工具"。 5. §2.6 简短小结 —— 直接指向 Ch 3 守恒方程推导。

作者特意在 §2.1 末尾强调："Throughout this book scalars are designated by lightface italic, vectors by lower boldface Roman, and tensors by boldface Greek letters" —— 这一记号约定在后续 18 章严格保持不变，是 Moukalled 这本教材"可读性 + 跨章一致性"的重要设计。

核心方程与概念

本章 5 节涉及的核心方程与概念，按"信号量"与"应用面"两层列举如下。

一、向量与张量基础（§2.2-§2.4）

  • 向量分量表示 (Eq. 2.1)

$\(\mathbf{v} = u \mathbf{i} + v \mathbf{j} + w \mathbf{k} \tag{2.1}\)$

其中 \((u, v, w)\) 是 \(\mathbf{v}\) 在 \((x, y, z)\) 方向的分量，\((\mathbf{i}, \mathbf{j}, \mathbf{k})\) 是单位基向量。Ch 3 的速度场即采用该符号系统。

  • 向量模 (Eq. 2.10)

$\(\lVert \mathbf{v} \rVert = \sqrt{\mathbf{v} \cdot \mathbf{v}} = \sqrt{u^2 + v^2 + w^2} \tag{2.10}\)$

  • 单位方向向量 (Eq. 2.11)

$\(\mathbf{e}_v = \frac{\mathbf{v}}{\lVert \mathbf{v} \rVert} \tag{2.11}\)$

在 FVM 中，每个控制体界面 \(f\) 都有其面积向量 \(\mathbf{S}_f = S_f \mathbf{n}_f\)，其中 \(\mathbf{n}_f\) 即对应单位方向向量。

  • 点积 (Eq. 2.7-2.9)

$\(\mathbf{v}_1 \cdot \mathbf{v}_2 = \lVert \mathbf{v}_1 \rVert \, \lVert \mathbf{v}_2 \rVert \, \cos(\mathbf{v}_1, \mathbf{v}_2) = u_1 u_2 + v_1 v_2 + w_1 w_2 \tag{2.9}\)$

在 N-S 方程的扩散项中，\(\Gamma \nabla \phi \cdot \mathbf{S}_f\) 即点积运算（标量通量 = 扩散通量向量在面法向的投影）。

  • 叉积 (Eq. 2.12-2.13)

$\(\mathbf{v}_1 \times \mathbf{v}_2 = \begin{vmatrix} \mathbf{i} & \mathbf{j} & \mathbf{k} \\ u_1 & v_1 & w_1 \\ u_2 & v_2 & w_2 \end{vmatrix} \tag{2.12}\)$

叉积的物理意义是"由两向量张成的平行四边形面积向量"。在 FVM 推导中不直接出现，但涡量 \(\boldsymbol{\omega} = \nabla \times \mathbf{v}\)、角动量等概念用到。

  • 张量转置、迹、对称 / 反对称分解 (Eq. 2.40-2.46)

对任意二阶张量 \(\mathbf{s}\)，有

$\(\mathbf{s} = \frac{1}{2}(\mathbf{s} + \mathbf{s}^T) + \frac{1}{2}(\mathbf{s} - \mathbf{s}^T) = \mathbf{s}_{sym} + \mathbf{s}_{asym} \tag{2.46}\)$

这一分解在 Ch 3 推导 N-S 方程的应力张量时被反复使用：\(\boldsymbol{\sigma} = -p \mathbf{I} + \boldsymbol{\tau}\)，其中黏性应力 \(\boldsymbol{\tau}\) 即可分解为对称（黏性耗散）与反对称（涡黏性贡献）两部分。

  • 张量点积与双点积 (Eq. 2.73-2.78)

$\(\mathbf{s} : \nabla \mathbf{v} = s_{xx} \frac{\partial u}{\partial x} + s_{xy} \frac{\partial u}{\partial y} + \cdots + s_{zz} \frac{\partial w}{\partial z} \tag{2.78}\)$

N-S 方程动量方程右侧的"黏性耗散项" \(\nabla \cdot (\mu \nabla \mathbf{u})\) 在严格推导时是双点积的形式，Ch 3 会用到。

二、五大基本定理（§2.5）

  • 梯度定理 (Gradient Theorem) (Eq. 2.79)

$\(\int_C \nabla s \cdot d\mathbf{r} = s(\mathbf{r}(b)) - s(\mathbf{r}(a)) \tag{2.79}\)$

沿闭合曲线积分为零。Ch 3 推导 N-S 方程的瞬态项时，本质上是在用这一定理将"随体导数"展开成"局部导数 + 对流项"。

  • Green 定理 (Eq. 2.80, 2.82)

$\(\oint_C (u \, dx + v \, dy) = \iint_R \left( \frac{\partial v}{\partial x} - \frac{\partial u}{\partial y} \right) dx \, dy \tag{2.80}\)$

向量形式

$\(\oint_C \mathbf{v} \cdot d\mathbf{r} = \iint_R (\nabla \times \mathbf{v}) \cdot d\mathbf{S} \tag{2.82}\)$

2D 情形的线-面积分互换。Ch 3 推导涡量方程时要用到。

  • Stokes 定理 (Eq. 2.83)

$\(\iint_S (\nabla \times \mathbf{v}) \cdot d\mathbf{S} = \oint_C \mathbf{v} \cdot d\mathbf{r} \tag{2.83}\)$

3D 情形的线-面积分互换，旋度场是核心。

  • 散度定理 (Divergence / Gauss Theorem) (Eq. 2.84) —— FVM 的母定理

$\(\int_V \nabla \cdot \mathbf{v} \, dV = \oint_S \mathbf{v} \cdot \mathbf{n} \, dS \tag{2.84}\)$

这是 FVM 一切推导的起点。Ch 3 推导守恒方程（积分形式）、Ch 5 推导 FVM 离散化、Ch 8 推导扩散项面通量、Ch 11 推导对流项面通量，都直接调用本定理。

作者还给出散度定理的两个重要推论：

$\(\int_V \nabla s \, dV = \oint_S s \, d\mathbf{S} \tag{2.85}\)$

$\(\int_V \nabla \cdot \mathbf{s} \, dV = \oint_S \mathbf{s} \cdot \mathbf{n} \, dS \tag{2.86}\)$

式 (2.85) 用于 Ch 3 推导标量方程（如温度方程、组分方程）的守恒形式，式 (2.86) 用于 Ch 3 推导动量方程的守恒形式。

  • Leibniz 积分规则 (Eq. 2.87-2.89) —— 动网格 / 多相流 VOF 的基础

$\(\frac{d}{dt} \int_{a(t)}^{b(t)} f(x,t) \, dx = \int_{a(t)}^{b(t)} \frac{\partial f}{\partial t} dx + f(b(t), t) \frac{db}{dt} - f(a(t), t) \frac{da}{dt} \tag{2.87}\)$

三维形式（动控制体）

$\(\frac{d}{dt} \int_{V(t)} f \, dV = \int_{V(t)} \frac{\partial f}{\partial t} dV + \int_{S(t)} f \, \mathbf{v}_s \cdot \mathbf{n} \, dS \tag{2.88}\)$

式 (2.88) 是 Ch 3 推导一般性守恒方程的核心工具 —— 当控制体 \(V(t)\) 本身随流体运动时（如 Lagrangian 描述、VOF 界面追踪、ALE 动网格），控制体内 \(\phi\) 的总变化率 = 局部时间导数 + 界面对流通量。本书虽然以 Euler 描述为主，但 (2.88) 仍是 Ch 3 §3.3 推导"随体导数 \(D\phi/Dt = \partial \phi / \partial t + \mathbf{u} \cdot \nabla \phi\)"的关键工具。

静止控制体情形 (Eq. 2.89)

$\(\frac{d}{dt} \int_V f \, dV = \int_V \frac{\partial f}{\partial t} dV \tag{2.89}\)$

本书绝大多数离散化都基于此简化形式。

关键结论

  1. 向量 / 张量 / 矩阵的记号约定是全书的语法基础：标量 \(\phi\)、向量 \(\mathbf{v}\)、张量 \(\mathbf{s}\)、矩阵 \(\mathbf{A}\)；乘积结果用 (scalar) / [vector] / {tensor} 区分。这一约定在 Ch 3-19 严格保持，是读者顺畅阅读后续章节的前置条件。
  2. 散度定理 (Gauss 定理) 是 FVM 的数学基石：本书 Ch 5-14 几乎所有面通量展开都直接调用 \(\int_V \nabla \cdot \mathbf{F} \, dV = \oint_S \mathbf{F} \cdot \mathbf{n} \, dS\) 这一恒等式，将控制体上的体积分转化为界面上的面积分。这是 FVM 区别于 FDM（直接对微分形式做差分）的根本所在。
  3. Green 定理与 Stokes 定理的工程价值：Green 定理是 2D 旋度计算与平面流场分析的基础，Stokes 定理是 3D 涡量方程积分形式的基础。这两个定理在 Ch 3 推导涡量输运方程时被反复使用。
  4. Leibniz 积分规则是动控制体分析的起点：固定网格 Euler 描述下可简化为 (2.89)，但理解 (2.88) 的完整形式对学习后续 VOF、ALE、动网格等高级主题至关重要。
  5. 张量的对称 / 反对称分解是流体力学应力张量处理的标准技术：\(\mathbf{s} = \mathbf{s}_{sym} + \mathbf{s}_{asym}\) 在 Ch 3 推导黏性应力 \(\boldsymbol{\tau}\) 的 Stokes 假设（黏性应力张量对称）时被显式使用。

挑战和开放性问题

  1. 张量与矩阵的区分 —— 在 FVM 离散化中，单元刚度矩阵 \(\mathbf{A}\)、界面通量矩阵 \(\mathbf{F}\) 等都是"矩阵"而非"张量"。如何在不同记号下严格区分"张量"（满足坐标变换规则的几何对象）与"矩阵"（仅是数组），在初学者中常引起困惑。本书通过符号约定部分缓解，但未深入讨论。
  2. 曲线坐标系下的向量 / 张量运算 —— 本章全部在 Cartesian 坐标系下展开，但在圆柱 / 球坐标系下，散度、旋度、梯度的具体形式会有附加项（如柱坐标下 \(\nabla \cdot \mathbf{v} = \frac{1}{r} \frac{\partial (r v_r)}{\partial r} + \frac{1}{r} \frac{\partial v_\theta}{\partial \theta} + \frac{\partial v_z}{\partial z}\)）。本书未在本章展开，但 Ch 6（网格）与 Ch 8（扩散项）会涉及。
  3. 高阶张量 —— 流体力学中虽然主要是 0 阶（标量）、1 阶（向量）、2 阶（张量），但在非线性应力、压电材料、复杂介质中会出现 3 阶、4 阶张量。本书未涵盖。
  4. Green 定理与 Stokes 定理的拓扑前提 —— 两者均要求 \(C\) 是简单闭合曲线、\(S\) 是单连通光滑曲面。在工程实际中遇到多连通区域、洞、带尖角边界时，定理的形式需要修正。本书在 Example 7 中演示了半圆+直线段组合的处理，但未深入讨论拓扑边界条件。
  5. Leibniz 规则在多相流 VOF 中的扩展 —— 动界面的法向速度 \(\mathbf{v}_s \cdot \mathbf{n}\) 在 VOF 中常需结合界面曲率、界面张力源项处理，本书未深入展开（VOF 不在本书覆盖范围内）。

个人反思与批判性分析

本章是本书最像教科书式"扫盲"的一章，对本科已修过工程数学的读者而言，5 大定理的回顾略显冗长（特别是 §2.2 的向量加减法与 §2.3 的张量分量表示，作者用了大量篇幅展示元素式推导，这部分对成熟读者可跳过）。但对自学者或跨专业读者（如机械专业转向 CFD 的初学者），这种"从零开始"展开的写法是有价值的。

从教学法的角度，本章有以下几个值得借鉴的写作特点：

  • 符号约定前置 —— 在展开任何公式之前先约定标量 / 向量 / 张量 / 矩阵的记号，避免后续章节因符号歧义引起混淆。这是工程流体力学教材与数学系向量分析教材的重要区别 —— 工程教材更注重"无歧义地传达信息"，数学系教材更注重"严格证明每一个命题"。
  • 5 大定理按"维度递进"组织 —— 1D 梯度定理 → 2D Green 定理 → 3D Stokes 定理 → 3D 散度定理 → 任意维 Leibniz 规则。这种"由低维到高维"展开的逻辑对初学者较友好，读者可以在 1D 情形下先建立"积分 = 函数值之差"的直觉，再推广到 2D 旋度场、3D 散度场。
  • 每个定理配 Example —— Example 7 用 Green 定理计算 \(2y^3 \, dx + 3xy^2 \, dy\) 沿半圆+直线段的线积分（结果 \(-3\pi/8\)），Example 8 用散度定理计算 \(\oint \mathbf{F} \cdot d\mathbf{S}\) 沿长方体表面（结果 39）。这两个例题都很有教学价值：Example 7 演示了"线-面积分互换"的非平凡几何计算；Example 8 演示了"直接积分 vs 散度定理"的计算量差异。

但本章也存在一些不足：

  • 张量的双点积 (Eq. 2.78) 推导过程繁琐，对流体力学的实际应用并不对称 —— 实际 N-S 方程推导中用到的双点积形式 \(\boldsymbol{\tau} : \nabla \mathbf{u} = \sum_{i,j} \tau_{ij} \partial u_i / \partial x_j\) 只是一种特殊形式（两个张量都是 2 阶），更一般的"高阶张量双点积"（如 4 阶弹性张量与 2 阶应变张量的双点积 \(C : \boldsymbol{\varepsilon}\)）未涉及。这对固体力学背景的读者是个遗憾。
  • Stokes 定理的几何直观不够 —— 数学教材（如 Stewart Calculus）会配 3-4 个图示说明"旋度场沿闭合曲线的环量 = 旋度在曲面上的通量"，本书仅给出 1 个示意图（Fig. 2.16），且对边界方向（右手规则）的解释相对简略。
  • Leibniz 规则的工程应用举例不够 —— 公式 (2.87-2.89) 在动网格、VOF、ALE 中有极其重要的应用，但本章仅在文字中提及"对于动控制体……"，没有具体算例。对初学者而言，理解"为什么 \(\frac{d}{dt} \int_{V(t)} \phi \, dV \neq \int_{V(t)} \frac{\partial \phi}{\partial t} dV\)"是非常重要的，本章未给出直接的算例展示这种差异。

总体而言，本章的写作风格"扎实、保守、工程友好"，是 Moukalled 这本教材"中等深度"定位的典型体现 —— 既不是 Ferziger-Peric 那种纯理论（每一节都有详尽证明），也不是 Anderson CFD Basics 那种科普化（几乎不展开公式），而是"在工程读者需要的范围内给到足够深度，超出范围则留白"。

重要参考文献

  • [X1] Aris R. (1962) Vectors, Tensors, and the Basic Equations of Fluid Mechanics. Prentice-Hall. （工程流体力学的经典向量 / 张量入门书，本书写作风格受其影响明显）
  • [X2] Bird R.B., Stewart W.E., Lightfoot E.N. (2002) Transport Phenomena, 2nd Edition. Wiley. （化工传递现象的经典教材，其向量演算与本构方程部分可与本章互补）
  • [X3] Batchelor G.K. (1967) An Introduction to Fluid Dynamics. Cambridge University Press. （理论流体力学经典，本章的张量分解思想与 Batchelor 一脉相承）
  • [X4] Spalding D.B. (1971) A Finite-Difference Procedure for the Solution of the Navier-Stokes Equations Using the Method of Stream Function - Vorticity. Imperial College. （FVM 数值方法原始文献）
  • [X5] Stewart J. (2015) Calculus: Early Transcendentals, 8th Edition. Cengage. （Stokes 定理与 Green 定理的几何直观参考）
  • [X6] Riley K.F., Hobson M.P., Bence S.J. (2006) Mathematical Methods for Physics and Engineering, 3rd Edition. Cambridge University Press. （工程数学标准参考书，对张量部分有更全面的展开）
  • [X7] Hughes W.F., Brighton J.A. (1991) Schaum's Outline of Theory and Problems of Fluid Dynamics, 3rd Edition. McGraw-Hill. （流体力学习题集，例题可作本章补充练习）
  • [X8] Kreyszig E. (2011) Advanced Engineering Mathematics, 10th Edition. Wiley. （梯度定理与散度定理的严格证明参考）