开个新坑：把Pytorch当科学计算库用

我就是来作死的。不知道为什么Taichi-Lang好像越来越不热门了，赶紧趁体量少的时候切到Pytorch，反正能帮忙并行计算就是好库（）

2026-02-02

在比我强得多的VSCode+Agent+GPT 5.2 Codex帮助下，起码一个起草的原型给出来了，一个使用相场建模了表面张力和别的两相流模拟，准备debug了，希望不要被坑死

2026-04-23

初步完成。他能跑了！！！

2026-04-24

准备开始自学，叫Gemini给了个初始的目录。肯定是错漏百出，毕竟LLM真的很傻逼，但我起码有个方向。

# 《计算流体力学与非线性偏微分方程：从物理建模到大规模稀疏求解》

## 预备篇：计算数学语言、量纲分析与离散思维启蒙
*为刚接触连续介质与数值计算的大二读者铺设平稳的认知斜坡。本章不引入新物理，仅统一符号、建立量纲直觉、明确浮点运算边界，并演示“连续公式如何安全落地为离散代码”。*
*   **P.1 符号约定与索引规范**
    *   向量/矩阵/张量的粗体与指标记法对照表；自由指标与哑指标的运算纪律。
    *   偏导数记号 $\partial_i$、$\nabla$、$D/Dt$ 的上下文切换规则。
*   **P.2 量纲分析与无量纲化阶梯**
    *   物理量的量纲矩阵与 Buckingham $\pi$ 定理的构造性证明。
    *   特征尺度选取策略：如何将任意 PDE 转化为无量纲形式，并自然浮现 $Re, Ma, Pr, Kn$ 等控制参数。
    *   无量纲化对数值条件数与迭代收敛性的隐性影响。
*   **P.3 浮点算术与离散误差源启蒙**
    *   IEEE 754 标准下的机器精度 $\epsilon_{mach}$、舍入误差累积与灾难性抵消（Catastrophic Cancellation）。
    *   条件数 $\kappa(A)$ 的几何直觉：为什么“数学上可解”不等于“计算机上可算”。
    *   从连续极限到离散网格的映射检查表：一致性、守恒性、有界性的初步定义。

---

## 导论：多尺度物理建模的数学抽象与连续介质假设
*本章奠定全书的理论基准，明确流体力学中不同物理尺度对应的控制方程及其基本假设。所有后续离散格式的物理合法性均溯源至此。*
*   **0.1 物质的微观离散性与分子动力学（MD）**
    *   薛定谔方程的降维：Born-Oppenheimer近似与哈密顿动力学。
    *   相空间（Phase Space）高维灾难与决定论动力学的计算极限。
*   **0.2 介观尺度的统计流体力学**
    *   从多体相空间到单体分布函数 $f(\mathbf{x}, \mathbf{v}, t)$ 的数学投影。
    *   分布函数的物理意义：位置与速度联合空间中的概率密度演化。
    *   Liouville定理、分子混沌假设（Stosszahlansatz）与玻尔兹曼（Boltzmann）积分-微分方程的推导。
*   **0.3 宏观尺度的连续统逼近**
    *   Knudsen数（$Kn$）的定义及其作为尺度划分的无量纲判据。
    *   连续介质假定（Continuum Hypothesis）：将离散粒子场粗粒化（Coarse-graining）为时空连续且中心极限定理成立的宏观场变量（密度、速度、温度）。
    *   统计力学向热力学的过渡：流体微团（Fluid Element）的数学定义。
*   **0.4 导论逻辑闭环与学习路径导引**
    *   微观→介观→宏观的投影算子链条总结。
    *   何时必须用离散粒子？何时可安全使用 PDE？$Kn$ 阈值的工程判据。
    *   顶刊对接提示：多尺度耦合论文（如 DSMC-NS 混合算法）的引言阅读范式。

---

## 第一卷：连续介质力学的数学重构与场论基础
*承接微积分与线性代数，建立描述流体运动变形规律的经典场论与张量数学体系。本卷所有推导均以“坐标无关性”为第一原则，确保后续离散格式不依赖特定网格拓扑。*

**第 1 章 张量代数与向量微积分进阶**
*   **1.1 坐标无关性与张量指标记法**
    *   爱因斯坦求和约定、Kronecker Delta 与 Levi-Civita 符号的代数恒等式。
    *   标量（零阶）、向量（一阶）与二阶张量的坐标变换协变性证明。
    *   二阶张量的特征值不变量（Principal Invariants）及其物理几何意义。
*   **1.2 经典微分算子及其场论本质**
    *   梯度算子（$\nabla$）：标量场的方向导数最大化与驱动势的数学定义。
    *   散度算子（$\nabla \cdot$）：向量空间中控制体积内通量的局部守恒度量。
    *   旋度算子（$\nabla \times$）：向量场局部旋转特性的表征。
    *   拉普拉斯算子（$\nabla^2$ 或 $\Delta$）：标量/向量场局部各向同性扩散的数学抽象。
*   **1.3 积分恒等式与全局-局部的转换**
    *   高斯散度定理（Gauss's Divergence Theorem）的严格定义与证明。
    *   斯托克斯定理（Stokes' Theorem）与闭合曲线积分。
    *   多维空间中的分部积分法则（Integration by Parts）与格林公式（Green's Identities）的推演。
*   **1.4 自解释桥梁：从张量恒等式到守恒律模板**
    *   为什么所有守恒律最终都呈现 $\partial_t (\cdot) + \nabla \cdot (\cdot) = \nabla \cdot (\cdot) + S$ 形式？
    *   指标记法手算演练：将 $\nabla \cdot (\rho \mathbf{u} \otimes \mathbf{u})$ 展开为分量形式并还原。

**第 2 章 流体运动学与变形几何学**
*   **2.1 观测参考系：欧拉描述与拉格朗日描述**
    *   物质坐标与空间坐标的映射雅可比（Jacobian）行列式。
    *   随体导数（Material Derivative, $D/Dt = \partial/\partial t + \mathbf{u} \cdot \nabla$）的链式法则推导；对流项（Advective term）非线性的数学起源。
*   **2.2 雷诺输运定理（Reynolds Transport Theorem）**
    *   随流控制体（Material Volume）积分随时间演化的 Leibniz 法则三维推广。
*   **2.3 速度梯度张量的加性分解**
    *   对称部分：应变率张量（Rate of Strain Tensor）与流体微团的线性拉伸/剪切变形。
    *   反对称部分：旋转张量（Vorticity Tensor）与涡量场的关系。
    *   运动学的无散流体条件：$\nabla \cdot \mathbf{u} = 0$ 等价于局部体积守恒的纯数学证明。
*   **2.4 自解释桥梁：运动学量的离散友好性评估**
    *   为什么涡量 $\boldsymbol{\omega} = \nabla \times \mathbf{u}$ 在交错网格上更易保持无散？
    *   拉格朗日粒子追踪与欧拉场插值的误差传递路径。

**第 3 章 宏观守恒定律的闭合系统推导**
*   **3.1 质量守恒：连续性方程**
    *   基于雷诺输运定理推导积分形式与微分形式的质量守恒（$\frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{u}) = 0$）。
*   **3.2 动量守恒：柯西运动方程与 Navier-Stokes 方程**
    *   微团队列表面应力张量（Cauchy Stress Tensor）与体积力。
    *   热力学平衡假设下的本构方程（Constitutive Equation）：广义牛顿流体定律。
    *   层流动力粘度 $\mu$ 与体积粘度 $\mu_v$ 的引入，斯托克斯假设（Stokes' Hypothesis）。
    *   Navier-Stokes (N-S) 动量方程的完整张量表达形式。
*   **3.3 能量守恒：总能、动能与内能方程**
    *   热力学第一定律的控制体积展开。
    *   热传导本构模型：傅里叶定律（Fourier's Law）与热导率 $k$。
    *   粘性耗散函数（Viscous Dissipation Function, $\Phi$）的张量推导：动能向内能的不可逆转化。
*   **3.4 方程组的闭合（Closure Problem）**
    *   热力学状态假设与量热状态方程（Caloric Equation of State, EoS）。
    *   完全气体方程与液体泰特方程（Tait Equation）的代数结构分析。
    *   6张场变量（1密+3速+1压+1温）对6组偏微分方程的系统闭合证明。
*   **3.5 自解释桥梁：守恒律的“强形式”与“积分形式”等价性检查**
    *   何时必须退回积分形式？（间断、非光滑边界、非结构网格）
    *   本构方程的线性假设边界：牛顿流体 vs 非牛顿流体的判据。

**第 4 章 偏微分方程的适定性与数学分类**
*   **4.1 阿达马适定性条件（Hadamard's Well-posedness）**
    *   解的存在性、唯一性，以及解对初边值条件的连续依赖性（稳定微扰准则）。
*   **4.2 二阶线性 PDE 的算子分类与物理特征**
    *   判别式理论初步。
    *   椭圆型（Elliptic）：如泊松方程。边值问题（BVP），无穷大信息传播速度，经典极值原理（Maximum Principle）的连续时间雏形。
    *   抛物型（Parabolic）：如热传导方程。初边值问题（IBVP），耗散结构，频谱的高阶衰减特性。
    *   双曲型（Hyperbolic）：如波动方程/对流方程。初试演化问题（IVP），特征线理论（Method of Characteristics），有限波速与间断平移。
*   **4.3 自解释桥梁：PDE 类型如何决定离散策略？**
    *   椭圆型→全局耦合/迭代求解；抛物型→时间步进/隐式耗散；双曲型→迎风/特征分解/限制器。
    *   混合类型方程（如不可压 N-S）的解耦哲学预告。

---

## 第二卷：空间离散矩阵化与逼近哲学的分歧
*连续偏微分方程在计算机架构下的降维实现。本卷详细论述三种截然不同的空间离散进路，它们各自根植于不同的数学土壤：泰勒展开分析、高斯积分守恒律，以及泛函空间变分原理。每章末尾提供“格式选择决策树”。*

**第 5 章 有限差分法（FDM）：微积分的代数退化**
*   **5.1 离散网格拓扑与泰勒展开定理**
    *   笛卡尔网格节点索引与连续函数的多阶 Taylor 展开式。
*   **5.2 差分算子与截断误差（Truncation Error）分析**
    *   前向、后向与中心差分格式的一阶导数模板（Stencil）。
    *   二阶导数中心差分的矩阵形式构建与大 O 符号（$\mathcal{O}(\Delta x^p)$）精度界定。
*   **5.3 修正方程（Modified Equation）与数值耗散/色散**
    *   将代数差分反向还原为修正的 PDE，揭示奇、偶数阶导数残留导致的“数值耗散”与“数值色散”。
*   **5.4 紧致差分格式（Compact FDM）**
    *   隐式导数方程系统：如何用较小模板宽度换取极高空间精度（Padé 逼近）。
*   **5.5 边界条件的差分离散**
    *   Dirichlet 与 Neumann 边界处理，虚拟网格（Ghost points）技术。
*   **5.6 自解释桥梁：FDM 的矩阵装配流水线**
    *   从 Stencil 到稀疏矩阵 $A$ 的索引映射规则（CSR 格式预备）。
    *   顶刊对接提示：高阶紧致格式在声学/湍流 DNS 中的应用场景。

**第 6 章 有限体积法（FVM）：守恒本源的积分忠诚**
*   **6.1 逆向思维：以积分形式代替微分形式**
    *   控制体积全局通量守恒性证明（Telescoping property）。FVM 在非结构多面体网格的绝对优势。
*   **6.2 高斯散度定理在网格面上的投影**
    *   体心（Cell-centered）数据的多边形面积向计算与梯度重构（格林-高斯方法）。
*   **6.3 通量插值（Flux Interpolation）与信息传递对流**
    *   迎风格式（Upwind scheme）的物理因果律；中心差分在纯对流方程下的虚假破缺。
    *   非正交交叉扩散项（Cross-diffusion）的补偿策略。
*   **6.4 自解释桥梁：FVM 的通量平衡检查表**
    *   如何验证代码是否真正守恒？（全局残差求和应为机器零）
    *   面心值重构的精度阶梯：一阶迎风 → 二阶线性 → MUSCL/WENO 预告。

**第 7 章 有限元法（FEM）：泛函拓扑与变分原理**
*   **7.1 变分法漫论：微积分的最优几何推广**
    *   泛函导数与欧拉-拉格朗日方程（Euler-Lagrange Equation）。
*   **7.2 加权残量法（Weighted Residual Method）**
    *   测试函数（Weighting function）的引入与逼近函数的残差正交化。
*   **7.3 伽辽金投影（Galerkin Projection）与弱形式（Weak Form）**
    *   降型约束：从强形式向弱形式转移的自洽性，能量范数内积的基础。
*   **7.4 局部基函数与单元装配**
    *   线性帽子函数（Hat function），整体刚度矩阵（Stiffness Matrix）的组装逻辑与其对称正定特征分析。
*   **7.5 自解释桥梁：FDM/FVM/FEM 的统一视角**
    *   在均匀网格上，二阶 FDM、线性 FVM 与线性 FEM 的代数等价性证明。
    *   何时选谁？几何复杂度、守恒要求、高阶精度需求的决策矩阵。

---

## 第三卷：时间的推演与数值稳定性、唯一性的审判
*在通过半离散化将 PDE 转化为 ODE 组之后，本卷系统解析时间积分的代数构造，确立稳定性极限，并首次在代数层面探讨高阶精度与解的唯一可解性。所有时间格式均与第二卷的空间离散显式耦合演示。*

**第 8 章 常微分方程的数值积分与二阶精度构造**
*   **8.1 初值问题的时间步进离散**
    *   局部截断误差与全局累积误差的映射关系。
*   **8.2 显式与隐式的一阶格式**
    *   前向欧拉（Forward Euler）的条件限制与后向欧拉（Backward Euler）的无条件隐式耗散。
*   **8.3 跨越时间的中轴：二阶精度（Second-order Accuracy）体系**
    *   梯形规则与 Crank-Nicolson 格式（CN格式）：时间对恒性（Time-Reversible）与其在抛物型方程中极难维持的虚假振荡阻尼缺陷。
    *   多步法之向后微分公式（BDF2）：牺牲对称性以换取对高频波成分强力压制的优越二阶刚性求解器。
    *   单步多级法：龙格-库塔方法（Runge-Kutta Methods）的 Butcher 表阵架构。
*   **8.4 自解释桥梁：半离散化（Method of Lines）实战映射**
    *   将 $\partial_t u = \mathcal{L}(u)$ 转化为 $\frac{d\mathbf{U}}{dt} = \mathbf{F}(\mathbf{U})$ 的完整代码结构模板。
    *   显式 vs 隐式的计算成本-稳定性权衡曲线。

**第 9 章 稳定性判定机制与系统唯一可解性**
*   **9.1 复平面稳定域与刚性系统问题**
    *   Dahlquist 测试方程特征值衰减分析。绝对稳定区间理论，A-稳定与 L-稳定的代数定义。
*   **9.2 频域波前分析：冯·诺依曼（Von Neumann）稳定性定理**
    *   复指数傅里叶分解，放大因子 $G(k \Delta x, \Delta t)$ 推导，稳定性严苛判据 $\max_k |G| \le 1$。
*   **9.3 CFL 条件与依赖域极限**
    *   Courant Number 对连续波传播数值依赖域的限制物理直觉。
*   **9.4 解的存在域：唯一稳定解（Unique Solvability）的代数定论**
    *   当隐式格式生成系统方程 $A\mathbf{x} = \mathbf{b}$ 时，如何通过 Lax-Milgram 引理（离散形式）或对角占优法则证明离散矩阵 $A$ 严格非奇异（严格正定阵）。
    *   证明唯一稳定解不仅是数值计算不至于“死机”的前提，更是原偏微分方程适定性在网格层面的绝对倒影。
*   **9.5 Lax 等价定理**
    *   一致相容（Consistency） + 数值稳定（Stability） $\implies$ 严格收敛（Convergence）。
*   **9.6 自解释桥梁：稳定性分析的标准化工作流**
    *   手写 Von Neumann 分析步骤清单；CFL 数自动计算代码片段。
    *   顶刊对接提示：如何在新格式论文中规范呈现稳定域与收敛阶测试。

---

## 第四卷：线性解耦系统的代数寻优与跨尺度引擎
*对大规模复杂PDE，直接耦合求解效率极低。只有通过科学解耦与克雷洛夫空间投影，才能征服千万级自由度的刚度怪兽。本卷将第三卷生成的线性系统转化为可并行求解的代数对象。*

**第 10 章 稀疏系统的结构属性与迭代基础**
*   **10.1 矩阵的几何谱特性**
    *   对称正定（SPD）、谱半径 $\rho(A)$、条件数 $\kappa(A)$（误差放大本原）的代数结构。
*   **10.2 不动点弛豫法（Relaxation）的频率偏好**
    *   雅可比（Jacobi）、高斯-赛德尔（Gauss-Seidel）与 SOR 迭代法。
    *   弛豫方法致命缺陷：迅速消灭高频锯齿，却对全局平缓的“低频长波误差”引发停滞综合征。
*   **10.3 自解释桥梁：稀疏矩阵存储与访存优化**
    *   CSR/CSC/COO 格式图解；为什么稀疏矩阵乘法（SpMV）是内存带宽受限型操作？
    *   缓存友好型循环重排与 OpenMP 并行基础。

**第 11 章 基于正交投影的克雷洛夫子空间法（Krylov Subspace）**
*   **11.1 Krylov 幂空间定义**
    *   扩张 $\mathcal{K}_m(A, r_0) = \text{span}\{r_0, Ar_0, A^2 r_0, \dots \}$。
*   **11.2 针对对称正定阵：共轭梯度法（CG）**
    *   优化视角：求解 $A\mathbf{x} = \mathbf{b}$ 等价于二次泛函最低点爬降。共轭方向保证不走回头路的几何正交原理。
*   **11.3 针对非对称对流阵：广义最小残量法（GMRES）**
    *   Arnoldi 过程推导正交基底并在 Hessenberg 矩阵中解取极小化。
*   **11.4 矩阵预处理（Preconditioning）机制**
    *   利用 $M^{-1}A\mathbf{x} = M^{-1}\mathbf{b}$ 强行聚敛特征值。ILU分解阈值截断策略。
*   **11.5 自解释桥梁：Krylov 方法的收敛诊断与重启策略**
    *   残差历史曲线解读；GMRES(m) 重启参数的工程选取法则。
    *   顶刊对接提示：PETSc/Trilinos 求解器接口配置逻辑。

**第 12 章 算法皇冠：多重网格方法（Multigrid Methods, MG）**
*   *保持计算复杂度 $O(N)$（网格数翻倍，时长仅翻倍）的不可思议引擎。*
*   **12.1 频率搬移原理（Frequency Shifting Theory）**
    *   在细网格上的低频长波停滞误差，下放至粗网格中即伪装为高频锯齿，从而被平滑器完美阻断。
*   **12.2 几何多重网格（GMG）**
    *   限制投影（Restriction）、插值外延（Prolongation）与 V-Cycle/W-Cycle。
*   **12.3 代数多重网格（AMG）**
    *   抛弃几何网格，纯依靠系数强弱关联抽象出图论代数粗结构。工业不可压缩流底层的核心主脑。
*   **12.4 自解释桥梁：多重网格的装配与调试指南**
    *   平滑器选择（Jacobi/GS/ILU）与粗化策略的匹配表。
    *   为什么 AMG 是黑盒？强连接阈值 $\theta$ 的敏感性分析。

---

## 第五卷：介观统计演化、渐近展开与格子玻尔兹曼方法（LBM）溯源
*绕开基于大量折中常数的 NS 模型，由“碰撞-迁移”体系切入多孔介质流体力学的介观降维本质。本卷首次系统展示如何通过多尺度渐近展开，从同一 kinetic 方程严格导出不同宏观流体模型。*

**第 13 章 平衡态分布、正交基展开与 Chapman-Enskog 渐近阶梯**
*   **13.1 麦克斯韦分布向离散空间的泰勒截断**。
*   **13.2 速度空间的高斯-厄米特正交（Gauss-Hermite Quadrature）**
    *   权函数导出宏观矩（积分密度与动量）。证明有限速度方向（如 D2Q9, D3Q27）的权重和即可精密还原宏观通量，强制约束网格必须为结构方阵。
*   **13.3 Chapman-Enskog 多尺度展开：从玻尔兹曼到宏观方程族的严格推导**
    *   展开参数 $\epsilon$ 的物理意义（Knudsen 数阶数）与多时间尺度假设 $t_0, t_1, t_2$。
    *   **零阶展开 $O(\epsilon^0)$：欧拉方程（Euler Equations）的涌现**。无粘、绝热极限下的动量与能量守恒；沿流线积分导出伯努利方程（Bernoulli Equation）的适用边界与数学前提。
    *   **一阶展开 $O(\epsilon^1)$：Navier-Stokes 方程的粘性闭合**。非平衡分布函数 $f^{(1)}$ 的矩积分如何精确生成牛顿粘性应力张量与傅里叶热流；输运系数（$\mu, k$）的微观表达式。
    *   **二阶展开 $O(\epsilon^2)$：Burnett 方程与超非平衡修正**。高阶应力/热流项的引入；Knudsen 层内的滑移边界与温度跳跃；Burnett 方程的线性不稳定性缺陷与正则化策略（R13/R26）。
    *   渐近一致性检查：为什么 LBM 天然恢复 NS 而非 Burnett？截断阶数与离散速度集的匹配法则。
*   **13.4 自解释桥梁：渐近展开的代数操作模板**
    *   分布函数展开 $f = f^{(0)} + \epsilon f^{(1)} + \epsilon^2 f^{(2)}$ 的逐阶代入与矩投影手算演示。
    *   顶刊对接提示：如何阅读含高阶滑移边界或稀薄气体效应的 PoF/JCP 论文。

**第 14 章 离散演化与宏观动力学涌现**
*   **14.1 LBM 迁移-碰撞算子（Stream-Collide Algorithm）**
    *   BGK 单弛豫时间模型，完美契合 GPU 的无偏微分纯代数映射。
*   **14.2 粘度反推的 Chapman-Enskog 展开**
    *   证明截断误差化身为“负粘度”，揭示粘度闭环公式 $\nu = c_s^2 (\tau - 0.5) \Delta t$ 的物理本源。
*   **14.3 全局部反弹边界格式（Bounce-Back Scheme）**。
*   **14.4 自解释桥梁：LBM 与宏观求解器的性能-精度对照**
    *   何时 LBM 碾压 FVM？（复杂孔隙、多相界面、GPU 众核）
    *   何时 LBM 失效？（高 $Ma$、高 $Re$ 湍流、强压缩性）

---

## 第六卷：工程求解全域解耦：压力、对流与间断捕获
*直面非线性方程组的计算噩梦，利用算子分裂体系将极为复杂的完全耦合组，优雅地砍散为一系列高效稳定的线性子问题。本卷将可压与不可压体系并列对照，揭示其底层数学结构的同源性。*

**第 15 章 非线性矩阵的微分解降术**
*   **15.1 次线性逼近：皮卡迭代（Picard Iteration）及其滞后场局限。**
*   **15.2 牛顿-拉夫逊法（Newton-Raphson Method）**
    *   多维雅可比矩阵（Jacobian）刚度矩阵组装，超二次收敛与初始刚性。
*   **15.3 工业巨核：JFNK (Jacobian-Free Newton-Krylov)**
    *   应对多物理全耦合方程无法写出解析雅可比的全维灾难。利用方向导数差分 $\frac{F(\mathbf{x} + \epsilon \mathbf{v}) - F(\mathbf{x})}{\epsilon}$ 逼近 $J\mathbf{v}$ 乘积，直接内嵌 GMRES。
*   **15.4 自解释桥梁：非线性迭代的收敛盆地与步长控制**
    *   线搜索（Line Search）与信赖域（Trust Region）的几何直觉。
    *   初始猜测构造策略：从零场、低 $Re$ 解或粗网格插值启动。

**第 16 章 不可压缩流的绝对壁垒：算子分裂法（Operator Splitting/Fractional Step）**
*   *压力失去状态方程，完全退化为迫使 $\nabla \cdot \mathbf{u} = 0$ 的约束乘子。*
*   **16.1 Stokes 问题与同位网格棋盘振荡的交错网格（MAC）救赎。**
*   **16.2 算子分裂方法（Operator Splitting Method）深探**
    *   不再试图全耦合。将流体更新裂解为“纯对流扩散的粘性步”与“施加散度约束的无粘投影步”。
    *   **一阶 Lie 分裂（Lie Splitting）**：Chorin-Temam 投影法推导。求解预测动量、调用泊松方程（AMG解压）及最终无散度映射修正。
    *   **二阶 Strang 分裂（Strang Splitting）**：通过“半步对流-全步扩散-半步对流”的几何对称排序，严格证明如何在完全分离方程独立求解的前提下，最终恢复整体时间推进的**二阶精度（Second-order Accuracy）**。
*   **16.3 商业软件中的 SIMPLE 算法阵列。**
*   **16.4 自解释桥梁：压力泊松方程的离散相容性检查**
    *   为什么 $\nabla \cdot \mathbf{u}^* = \nabla^2 p$ 的右端项必须满足离散零均值条件？
    *   边界压力梯度的隐式处理与质量守恒残差监控。

**第 17 章 可压缩动力学与激波的双曲之战**
*   **17.1 黎曼问题与熵条件（Oleinik Entropy Condition）。**
*   **17.2 Godunov 定理与二阶对称矩阵在跃变处的数值灾难。**
*   **17.3 非线性自适应通量限制（Limiters）**
    *   依靠网格探测在间断处降一阶防止震荡，在平滑处弹回高阶精度的 TVD 与 MUSCL 重构理论。
*   **17.4 WENO 高阶加权近似黎曼求解。**
*   **17.5 自解释桥梁：可压/不可压求解器的架构分岔点**
    *   状态方程的角色切换；全耦合密度基求解器 vs 压力基求解器的适用 $Ma$ 阈值。
    *   顶刊对接提示：激波-边界层干扰（SBLI）论文的数值格式披露规范。

---

## 第七卷：重铸当代物理顶刊黑科技：极值、能量与结构保持（Structure-Preserving）
*这是现代计算数学的皇冠。单纯减小泰勒截断误差只是初级工巧；真正的“道”在于让离散代数体系如同连续物理定律一般绝对忠诚。只有在离散层面严格复现宏观不等式法则，偏微分程序才能跨越刚性，万古长存。本卷所有算法均提供“物理不等式→离散代数约束→线性解耦”的完整推导链。*

**第 18 章 浓度场的生死线：极值原理与有界性保持（Maximum Principle & Bound-Preserving）**
*   *体积分数、化学浓度如果在离散演化中越界（$<0$ 或 $>1$），将导致状态方程输出复数并瞬间摧毁程序。这比简单的稳定性要求更为苛刻（即 MPP, Maximum Principle Preserving）。*
*   **18.1 经典极值原理向全离散系统的等价判定**
    *   单调矩阵（$M$-Matrix）理论推演：为什么隐式扩散算子天然拥有保界性特性？
    *   时间步限制下的对流溢出（Overshoot / Undershoot）分析。
*   **18.2 构造不等式方法（Inequality Construction Method）**
    *   一种严密的代数钳制逻辑：将目标节点离散层面的更新值，严格证明可展开为包围节点上一时间步值的**全凸组合（Convex Combination，即所有组合权重系数 $\ge 0$ 且和为 1）**。
    *   通过此不等式的构造性放缩限制，证明数值解的上/当下界绝不会超越前续迭代。
*   **18.3 截断限制器与算法保界修正（Algorithmic Limiting）**
    *   拒绝野蛮的后处理“剪裁（Clipping）”（它破坏质量守恒）。引入通量矫正传输技术（FCT）进行自适应的有界性重构。
*   **18.4 自解释桥梁：MPP 格式的调试与验证清单**
    *   如何构造违反极值原理的测试用例？凸组合权重的运行时监控代码。

**第 19 章 复杂软物质界面：Cahn-Hilliard 相场建模的非线性刚性势阱**
*   *模拟油水不溶边界时避开几何重构噩梦的弥散界面模型论。*
*   **19.1 吉布斯自由能泛函与双势阱（Double-Well Potential）结构**
    *   $F(\phi) = \frac{1}{4}(\phi^2-1)^2$ 的惩罚拔河拉锯战。
*   **19.2 $H^{-1}$ 梯度流与质量守恒的四阶高频耗散PDE推演。**
    *   动能张力：化学势如何回馈 Navier-Stokes，引发马兰戈尼破裂与毛细凝聚涡旋。
*   **19.3 自解释桥梁：四阶方程的降阶离散策略**
    *   引入辅助变量 $\mu = \delta F/\delta \phi$ 将四阶 PDE 拆分为两个耦合二阶 PDE 的代数等价性。
    *   界面厚度参数 $\epsilon$ 与网格分辨率 $\Delta x$ 的匹配法则（$\Delta x \le \epsilon/3$）。

**第 20 章 热力学第二定律的离散法庭：能量稳定（Energy Stability）算法群**
*   *自由能必须随时间单调递减（$\frac{dE}{dt} \le 0$）。若强行应用显式泰勒展开项，差分系统的代数能量将在某个时间步意外增高，从而引发灾难性爆炸。*
*   **20.1 凸分裂方法（Convex Splitting Method）**
    *   Eyre大神的绝妙思维放缩：由于双阱势非凸而难以求极值，我们强制将其拆分为一个“绝对凸的抛物线项”（必须使用全隐式处理以释放稳定耗散阻力）以及一个“上凸的凹项”（使用显式外推以降低求解难度）。
    *   通过差分法引理，数学上严格保证该离散组合所求数值解的能量随系统迭代不可逆递减（条件/无条件稳定）。
*   **20.2 人工稳定项方法（Artificial Stabilization Method）**
    *   规避复杂凸分裂分解难度的捷径构造法。
    *   在隐式-显式（IMEX）框架中强行加入一个满足一阶/二阶相容原理的纯数值稳定项，例如：$S(\phi^{n+1} - \phi^n)$（其中常数 $S$ 即为无穷大惩罚参数）。依靠该项引入的离散高频强阻尼，抵消显式非线性带来的能量反跳。
*   **20.3 SAV 与 IEQ 黑魔法技术的革命框架**
    *   **不变能量二次化（IEQ）**与**标量辅助变量法（SAV, Scalar Auxiliary Variable Method）**。
    *   当代顶刊大杀器核心奥义：定义一个人造标量变量 $R(t) = \sqrt{E(\phi) + C}$，利用微分链式法则彻底重写自由能相关的高次多项式项。
    *   **线性解耦与高效性（Linear Decoupling & Efficiency）体现**：使得原极为庞大、需要复杂牛顿非线性内迭代的四阶微分组，在每一次时间层仅退化蜕变为“求解常系数线性泊松方程叠加极其廉价的标量代数更新方程”。
    *   **最终成就**：实现**绝对无条件的能量极小化稳定（Unconditionally Energy Stable）**、全线性降维、而且完美保持系统系数矩阵处于 SPD 域从而确认每一时间步拥有唯一稳定解（Unique Solvability）。标志着超刚性问题解算正式进入计算自由时代。
*   **20.4 自解释桥梁：能量稳定格式的通用设计范式**
    *   从连续能量律 $\frac{dE}{dt} = -\mathcal{D} \le 0$ 到离散能量差 $E^{n+1}-E^n = -\Delta t \mathcal{D}_h$ 的构造步骤。
    *   顶刊对接提示：CMAME/JCP 中“Unconditionally Energy Stable”声明的数学审查清单。

**第 21 章 几何界面追踪的降维超脱：散步阈值动力学（Threshold Dynamics / MBO Method）**
*   *除了复杂的非线性四阶相场，还有更疯狂的高效追踪方法吗？*
*   **21.1 平均曲率流（Mean Curvature Flow）的纯几何本质**。
*   **21.2 物理算子分裂（Operator Splitting）在前端算法的高能演绎：MBO 方法**
    *   利用算子分裂哲学的极致，不解非线性势阱。而是将动态界面的平滑问题硬生生劈开。
    *   第一步（扩散算子）：只需无脑演化一个最经典的常数系数隐式热方程（平滑界面）。
    *   第二步（阈值截断算子 Thresholding）：应用 Heaviside 阶跃函数暴力二值化（重新削出尖锐界面）。
*   **21.3 绝境突围：不可思议的网格逼近**
    *   证明这“简单的漫步与砍削”循环竟然严格一致逼近极其困难的高维平均曲率流；在处理极其大量微反应多胞（如细胞动力学、金相结晶晶界演化）时，彻底碾压传统 Level-set 与 CH 相场的运算阻力边界，达到降维式模拟极速。
*   **21.4 自解释桥梁：MBO 与相场方法的精度-成本权衡**
    *   阈值步长的几何误差分析；何时 MBO 会丢失拓扑细节？
    *   多相扩展与体积守恒修正策略。

---

## 终章：从离散代码到顶刊论文的跨越
*全书知识闭环。将前七卷的数学工具、离散格式、求解器与结构保持哲学熔铸为一套可复现、可验证、可发表的计算工作流。*
*   **T.1 数值实验的标准化范式**
    *   网格收敛性研究（GCI）与阶数验证；制造解法（MMS）构造任意复杂 PDE 的精确测试源项。
    *   残差监控、守恒量漂移检查与能量/极值演化曲线的学术绘图规范。
*   **T.2 开源架构对接与高性能计算入口**
    *   从手写 MATLAB/Python 原型到 C++/Fortran 生产代码的迁移路径。
    *   MPI 区域分解、GPU CUDA 核函数映射、PETSc/Hypre 求解器调用的最小可行模板。
*   **T.3 顶刊阅读与复现指南**
    *   如何拆解 JCP/PoF/CMAME 论文的“格式声明→稳定性证明→基准测试→工程应用”四段式结构。
    *   常见“数学严谨但工程失效”的陷阱识别；如何向审稿人证明你的离散格式真正 Structure-Preserving。
*   **T.4 全书逻辑总览与进阶路线图**
    *   物理建模→场论语言→空间离散→时间积分→线性求解→非线性解耦→结构保持的完整知识图谱。
    *   致读者：计算流体力学不是调用黑盒，而是用代数忠诚复刻自然法则。愿您推开大规模方程演算架构的新世界大门。

*(全书完。全卷每一处离散公式的映射、每一个不等式界点的推演、每一段渐近展开的阶数匹配，皆以“直觉引导→严格推导→离散实现→稳定性闭环→顶刊对接”为唯一编排逻辑。确保零基础大二读者可沿平滑阶梯直达前沿研究。)*