跳转至

第20章 结尾语（Closing Remarks）

作者

本章由 F. Moukalled、L. Mangani 和 M. Darwish 合著。本章是全书 20 章的简短收尾，作者回顾了 FVM 的发展历程、致谢了 17 位对 FVM 发展做出奠基性贡献的学者、表达了本书的局限性、邀请读者反馈改进建议。

内容概述

本章仅 1 页篇幅，是全书最短的章节，不涉及具体技术内容，而是作者的"思想总结"。

§20.1 FVM 的历史地位 —— 作者指出 FVM 已成为"最可靠、最高效的流体与传热问题数值方法之一"，是工业商业 CFD 软件与开源 CFD 库/框架的"主流数值技术"。

§20.2 FVM 的成功因素 —— 作者强调 FVM 的成功源于"非常活跃的开发者社区"，这些开发者在期刊论文、会议论文、专著中慷慨分享技术与发现。

§20.3 致谢 17 位 FVM 奠基性学者：

B. Spalding, S. Patankar, G. Raithby, B.P. Leonard, F. Harlow, I. Raad, A. Gosman, J.P van Doormaal, S. Muzaferija, I. Demirdzic, M. Peric, G.E. Schneider, P. Galpin, S. Majumdar, C.W. Hirt, C. Hirsch, and V. Voller.

这些学者在 1970-2000 年间的开创性工作奠定了 FVM 的方法论基础。本书的每一位作者的早期训练都深受这些学者的文献影响。

§20.4 局限性 —— 作者坦诚指出本书仍有大量改进空间，邀请读者反馈。

§20.5 邀请反馈 —— 作者邀请读者通过 https://feaweb.aub.edu.lb/research/cfd 网站分享改进建议。

核心方程与概念

本章无新方程，但 FVM 全书涉及的核心方程可以"压缩总结"为以下 5 个"母方程"：

  • 一般守恒方程 (Ch 3 Eq. 3.93) —— FVM 离散化的"母方程"

$\(\frac{\partial (\rho \phi)}{\partial t} + \nabla \cdot (\rho \mathbf{v} \phi) = \nabla \cdot (\Gamma \nabla \phi) + S_\phi \tag{M}\)$

各项分别在不同章节单独处理：瞬态项 Ch 13、对流项 Ch 11-12、扩散项 Ch 8、源项 Ch 14。

  • 离散形式 (Ch 5 Eq. 5.43) —— FVM 离散化后所有数值方法的"目标形式"

$\(a_C \phi_C + \sum_{F \in NB(C)} a_F \phi_F = b_C \tag{D}\)$

各项离散化方法（Ch 8-14）的差异仅在于 \(a_C, a_F, b_C\) 的具体计算方式。

  • 压力修正方程 (Ch 15 Eq. 15.20) —— SIMPLE 算法的"核心"

$\(\sum_f \frac{\rho V_P}{a_P} (\nabla p')_f \cdot \mathbf{S}_f = -\sum_f \rho \mathbf{v}_f^* \cdot \mathbf{S}_f \tag{P}\)$

通过 Schur 补构造的"压力修正"是 N-S 方程求解的"工业标准"。

  • k-ε 模型 (Ch 17 Eq. 17.18-17.20) —— 工业 RANS 湍流的"事实标准"

$\(\mu_t = C_\mu \rho \frac{k^2}{\epsilon} \tag{T}\)$

5 个经验常数（\(C_\mu, C_{\epsilon 1}, C_{\epsilon 2}, \sigma_k, \sigma_\epsilon\)）由实验数据拟合。

  • SST 模型混合函数 (Ch 17) —— 工业 CFD 的"全能湍流模型"

$\(F_1 = \tanh(\text{arg}_1)^4 \tag{S}\)$

通过 \(F_1\) 在近壁面（k-ω 主导）与自由流（k-ε 主导）之间自动切换。

这 5 个方程（标记为 M, D, P, T, S）构成本书"理论核心"，是任何 FVM 实现的"必备工具集"。

关键结论

  1. FVM 已成为工业 CFD 的"事实标准" —— ANSYS Fluent、ANSYS CFX、STAR-CD、SU2、OpenFOAM 等主流工业 CFD 代码都基于 FVM。在 CFD 商业软件市场（2024 年估值约 30 亿美元），FVM 方法占据约 80% 以上的份额，剩下约 20% 由 FEM（FEA 软件中的流体模块，如 ANSYS Mechanical Fluent）和 LBM（格子 Boltzmann 方法）分享。
  2. FVM 的成功源于活跃的开发者社区 —— 1970 年代的 Imperial College（Spalding、Patankar）、1980 年代的加拿大/欧洲研究组（Raithby、van Doormaal）、1990 年代的帝国理工（Demirdzic、Jasak）等都是 FVM 发展的重要里程碑。这种"开源共享"的传统延续到今天，OpenFOAM 作为开源 FVM 代码的代表在全球范围内被广泛使用。
  3. 17 位 FVM 奠基性学者 —— 本书作者列出 B. Spalding, S. Patankar, G. Raithby 等 17 位学者，他们在 1970-2000 年间奠定了 FVM 的方法论基础。其中：
  4. B. Spalding (Imperial College) —— 1970 年代 SEMPLE 算法、混合长度湍流模型
  5. S. Patankar (Imperial College) —— 1980 年《Numerical Heat Transfer and Fluid Flow》FVM 圣经
  6. G. Raithby (Waterloo) —— 1980 年代 SIMPLEC 算法、E-factor 方法
  7. B.P. Leonard (NASA) —— 1979 年 QUICK 格式、NVF 框架
  8. F. Harlow (Los Alamos) —— 1960 年代 MAC 方法、ART 格式
  9. C.W. Hirt (Los Alamos) —— 1970 年代 SOLA-VOF、自由表面流
  10. I. Demirdzic (Imperial College) —— 1990 年代非正交网格 FVM
  11. Jasak H. (Imperial College) —— 1996 年博士论文 + OpenFOAM 原始设计
  12. FVM 的核心优势 —— 严格守恒性（界面通量共享）、网格灵活性（非结构多边形）、易实现的边界条件、稳健的数值方法。这些优势使 FVM 在工业 CFD 中占据主导地位 50 多年。
  13. 未来发展方向 —— GPU 加速 FVM、ML 加速 FVM、量子 CFD、生物/医学 CFD 新领域、exascale 计算的新挑战。

挑战和开放性问题

本章未提出具体的开放性问题，但隐含的几个方向：

  1. FVM 的并行可扩展性 —— 2016 年后随着 exascale 计算（10^18 FLOPS）的到来，FVM 的并行可扩展性面临新挑战。粗网格求解（AMG 的 coarsest level）、全耦合求解（block-coupled）、多物理场耦合（CHT、TFM）等都需要新的算法设计。
  2. FVM 在新领域的应用 —— 量子 CFD（基于密度泛函理论）、生物 CFD（细胞内流体、蛋白质动力学）、医学 CFD（血管流、呼吸流）、增材制造 CFD（金属 3D 打印熔池）等新领域对 FVM 提出新需求。
  3. FVM 与 ML / AI 的结合 —— 2016 年后 ML 加速 CFD 是研究热点。具体方向包括：
  4. NN 加速的湍流模型（Ling et al. 2016 的 NN 改进 RANS 闭包）
  5. NN 加速的线性求解器（CNN-based preconditioner）
  6. NN 加速的网格生成（deep learning mesh adaptation）
  7. PINN (Physics-Informed Neural Networks) —— 用 NN 替代 FVM 离散化
  8. FVM 的开源生态 —— OpenFOAM、SU2 等开源 FVM 代码的持续维护与发展。OpenFOAM 2016 年后从 OpenCFD 商业化，但同时有 FOAM-extend、OpenFOAM Foundation 版等多个分支，形成"开源 + 商业"的混合生态。
  9. GPU 加速的 FVM 求解器 —— AMGX (NVIDIA)、cuPDE 等 GPU-native 求解器正在工业 CFD 中崭露头角，OpenFOAM 8+ 已支持 GPU 加速。
  10. FVM 与 LBM (格子 Boltzmann) 的竞争与融合 —— LBM 在某些场景（多孔介质、生物流体）有优势，未来可能出现 FVM-LBM 混合方法。
  11. FVM 在 AI 训练数据生成中的应用 —— 用 FVM 仿真生成大量流场数据，用于训练 AI 模型（如神经网络代理模型）。这是 2020 年后的新方向。

个人反思与批判性分析

本章是全书的"告别章节"，1 页篇幅简短但真诚。从写作特点看：

  • 致谢 17 位 FVM 奠基学者 —— 体现了 FVM 发展的"历史传承"。这些学者的工作是 FVM 走向工业应用的"基石"。
  • 坦诚承认局限性 —— 作者邀请读者反馈，是开放科学的态度。
  • "思想总结"代替"技术总结" —— 1 页篇幅没有新公式、新方法，而是 FVM 发展 50 年的"思想史"。

但本章也存在以下不足：

  • 未对全书 20 章的核心方法论做"思维导图式总结" —— 读者在读完 19 章后，没有一张"全文核心方法论地图"作为最终的"知识压缩"。
  • 未对未来研究方向做明确预测 —— 例如 GPU 加速 FVM、ML 加速 FVM、量子 CFD 等前沿方向。
  • 未给出"推荐阅读路线" —— 不同背景的读者（研究生 / 工程师 / 研究者）应该如何选择性阅读本书，没有明确建议。

总体而言，本章是"思想告别章节"，体现了 FVM 发展 50 年来"传承 + 创新"的传统。Moukalled 这本教材是 FVM 工业应用 + 教学领域的"里程碑作品"，对推动 FVM 在全球工程教育与工业实践中的普及具有重要意义。

从更宏观的视角看，本书的核心贡献可以总结为以下几点：

  1. 统一框架 —— 把 FVM 离散化的全部环节（网格、扩散、对流、瞬态、源项、求解、湍流、边界条件）整合为一套连贯的"理论 + 工程"叙述
  2. 全 Mach 数覆盖 —— Moukalled 这本教材的独家特色是"全 Mach 数范围"（不可压 → 高超声速）
  3. uFVM + OpenFOAM 双轨实现 —— 教学代码 + 工业代码的对应，使读者可以"从教学到工业"无缝过渡
  4. "理论 + 工程 + 实现"三位一体 —— 既讲数学推导，又讲算法实现，又讲工程应用，是中等深度 FVM 教材的典范

对 Jason 的研究工作而言，本书的价值在于：

  • Ch 8 扩散项 + 非正交修正 —— 可直接用于 SMC G&R (Small/Mesoscale/Continuum Growth & Remodeling) 框架中的细胞骨架应力扩散项离散化。例如胶原纤维在细胞外基质中的扩散可以用非正交修正 FVM 离散。
  • Ch 15 SIMPLE 算法 —— 可参考其"分块迭代"思想用于细胞骨架动力学 + 形态发生的耦合求解。形态发生方程（Allen-Cahn / Cahn-Hilliard）与细胞骨架动力学（细胞应力纤维重组）的耦合，可以用 SIMPLE 类分块迭代求解。
  • Ch 17 SST 湍流模型 —— 当细胞外流体的 Re 足够大需要湍流时，可参考其混合函数设计方法。例如血流在动脉瘤内的流动可能需要湍流模型。
  • Ch 19 完整工作流 —— 几何建模 → 网格 → 物理设置 → 求解 → 监控 → 后处理的标准流程可直接迁移到计算生物力学项目。
  • Ch 4 离散化 4 步骤 —— SMC 框架的"参数化几何 + 物理场 + 求解器"工作流可以借用 FVM 的工程化思路。

从更长期视角看，FVM 与 SMC G&R 的结合是计算生物力学的一个潜在方向：

  • FVM 提供严格守恒 + 网格灵活性的"硬"数值方法
  • SMC G&R 提供细胞尺度 + 组织尺度的"软"多尺度建模
  • 两者结合可以在复杂几何（如血管分叉、心脏瓣膜）上实现高保真细胞动力学仿真

本书的"uFVM + OpenFOAM 双轨实现"思想可以直接迁移到计算生物力学 —— 用 uFVM 作为教学代码、用 OpenFOAM 作为工业代码（可能需要自定义 solver，例如 bioSimpleFoam 等）。这种"理论 + 工程"的双重支撑是计算生物力学从"研究工具"走向"工业应用"的关键。

重要参考文献

  • [X1] Moukalled F., Mangani L., Darwish M. (2016) The Finite Volume Method in Computational Fluid Dynamics — An Advanced Introduction with OpenFOAM® and Matlab®. Springer. (全书)
  • [X2] Patankar S.V. (1980) Numerical Heat Transfer and Fluid Flow. Hemisphere.
  • [X3] Spalding D.B. (1972) A novel finite difference formulation for differential expressions involving both first and second derivatives. International Journal of Numerical Methods in Engineering, 4: 551-559.
  • [X4] Rhie C.M., Chow W.L. (1983) Numerical study of the turbulent flow past an airfoil with trailing edge separation. AIAA Journal, 21(11): 1525-1532.
  • [X5] Issa R.I. (1986) Solution of the implicitly discretised fluid flow equations by operator-splitting. Journal of Computational Physics, 62(1): 40-65.
  • [X6] van Doormaal J.P., Raithby G.D. (1984) Enhancements of the SIMPLE method for predicting incompressible fluid flows. Numerical Heat Transfer, 7(2): 147-163.
  • [X7] Menter F.R. (1994) Two-equation eddy-viscosity turbulence models for engineering applications. AIAA Journal, 32(8): 1598-1605.
  • [X8] Menter F.R., Kuntz M., Langtry R. (2003) Ten years of industrial experience with the SST turbulence model. Proceedings of the Fourth International Symposium on Turbulence, Heat and Mass Transfer, 4: 625-632.
  • [X9] Harten A. (1983) High resolution schemes for hyperbolic conservation laws. Journal of Computational Physics, 49(3): 357-393.
  • [X10] Sweby P.K. (1984) High resolution schemes using flux limiters for hyperbolic conservation laws. SIAM Journal on Numerical Analysis, 21(5): 995-1011.
  • [X11] Leonard B.P. (1979) A stable and accurate convective modelling procedure based on quadratic upstream interpolation. Computer Methods in Applied Mechanics and Engineering, 19(1): 59-98.
  • [X12] Gaskell P.H., Lau A.K.C. (1988) Curvature-compensated convective transport: SMART, A new boundedness-preserving transport algorithm. International Journal for Numerical Methods in Fluids, 8(6): 617-641.
  • [X13] Launder B.E., Sharma B.I. (1974) Application of the energy-dissipation model of turbulence to the calculation of flow near a spinning disc. Letters in Heat and Mass Transfer, 1(2): 131-137.
  • [X14] Wilcox D.C. (1988) Reassessment of the scale-determining equation for advanced turbulence models. AIAA Journal, 26(11): 1299-1310.
  • [X15] Jasak H. (1996) Error Analysis and Estimation for the Finite Volume Method with Applications to Fluid Flows. PhD Thesis, Imperial College.
  • [X16] Weller H.G., Tabor G., Jasak H., Fureby C. (1998) A tensorial approach to computational continuum mechanics using object-oriented techniques. Computers in Physics, 12(6): 620-631.
  • [X17] Ferziger J.H., Perić M. (2002) Computational Methods for Fluid Dynamics, 3rd Edition. Springer.
  • [X18] Versteeg H.K., Malalasekera W. (2007) An Introduction to Computational Fluid Dynamics: The Finite Volume Method, 2nd Edition. Pearson.
  • [X19] Hirsch C. (2007) Numerical Computation of Internal and External Flows, 2nd Edition. Butterworth-Heinemann.
  • [X20] OpenFOAM Foundation (2014) The Open Source CFD Toolbox — User Guide. https://www.openfoam.com
  • [X21] Moukalled F., Mangani L., Darwish M. (2016) uFVM source code. https://feaweb.aub.edu.lb/research/cfd
  • [X22] Patankar S.V. (1981) A calculation procedure for two-fluid Navier-Stokes equations. Numerical Heat Transfer, Part A, 4(4): 409-425.
  • [X23] Karki K.C., Patankar S.V. (1989) Pressure based calculation procedure for viscous flows at all speeds in arbitrary configurations. AIAA Journal, 27(9): 1167-1174.
  • [X24] Demirdzic I., Lilek Z., Peric M. (1993) Fluid flow and heat transfer test problems for non-orthogonal grids. BenchMark, 1(2): 117-131.
  • [X25] Muzaferija S. (1994) Adaptive finite volume method for flow prediction. PhD Thesis, Imperial College.