第 3 章 · 不可压无粘流：最简单的完整案例

这一节把方程砍到最简：丢掉粘性（无粘）、丢掉压缩性（不可压），得到拉普拉斯方程。它是唯一能同时给出"解析解 + 数值解"并排比较的场合，是检验你离散代码的第一个考场。

从方程到拉普拉斯

不可压无粘无旋流的连续方程与无旋条件：

速度势 $\phi$ 满足拉普拉斯方程。它的数学性质是椭圆：任何一点的变化瞬间影响全场，没有特征线、没有波。所以求解必须全场联立（或大量迭代），不能像双曲问题那样一路推进。

基本解叠加：解析一端

直匀流、点源、点汇、偶极子、点涡——五个基本解各满足拉普拉斯方程，而拉普拉斯方程是线性的，所以基本解可以任意叠加。圆柱绕流 = 直匀流 + 偶极子；升力圆柱 = 再加一个点涡。这就是经典势流理论的全部玩法。

数值一端：面板法与网格法

• 源面板法：把物面切成 N 段，每段布置强度未知的源，令物面法向速度为零，得到 N 元线性方程组。这是工业软件 XFOIL 类方法的原型。
• 网格差分法：在矩形网格上用二阶中心差分离散 $\nabla^2\phi$，得到大型稀疏线性方程组，用高斯-赛德尔或 SOR 迭代解。

这个五点格式是 CFD 里最重要的一个式子，务必亲手推一遍泰勒展开，确认截断误差是 $O(\Delta x^2)$。

老师的提醒

• 无粘流里圆柱阻力为零（达朗贝尔佯谬），这不是代码错了，是模型（丢掉粘性）的错。数值结果违反物理直觉时，先检查模型假设。
• 中心差分对拉普拉斯方程好用，但对对流主导的问题会发散——同一个格式不能到处套用。

自检

• 为什么势流方程是椭圆型？这对求解方式意味着什么？
• 手推五点差分格式的截断误差阶数。
• 源面板法的线性方程组是怎么来的？

更深一层：势流理论的繁荣与它的阿喀琉斯之踵

势流理论在 19 世纪末达到了惊人繁荣：环量定理、镜像法、保角变换、基本解叠加——整个无粘不可压无旋流的理论大厦优雅到近乎完美，然后撞上了达朗贝尔佯谬：理想流体中任意物体运动不受阻力。这显然荒谬（ otherwise 船不用帆划），根源在于无粘假设强迫流场前后对称、压力积分相消。理论越美，佯谬越刺眼——历史上正是这个佯谬逼出了普朗特的边界层理论（第 4 节伏笔）：粘性只在很薄的一层里起作用，但它决定了尾流与压差阻力的存在。

这段历史的工程启示：当你发现模型结论与直觉剧烈冲突时，先别怀疑计算，怀疑模型假设。把假设一条条列出来（无粘？定常？无旋？），问每一款能不能松——达朗贝尔佯谬的解决史就是"放松无粘假设的局部化"的历史。

推导细读：源面板法的线性方程组是怎么长出来的

把物面离散成 $N$ 段直线面板，每段布置常强度源 $\lambda_j$，控制点取各段中点。第 $i$ 段控制点处的法向速度 = 来流法向分量 + 所有源段的贡献：

影响系数 $A_{ij}$ 有解析表达式（直线段源对其本身法向的贡献是 $\frac12$——自诱导项），于是 $N$ 元线性方程组一次解出全部源强，切向速度、压力系数随之而得。整个过程没有迭代、没有网格体积、没有边界层——这是"边界元"思想的最纯形态：只离散边界，维数降一维。数值实验：$N = 20$ 段的圆柱表面 $C_p$ 分布与解析解 $1 - 4\sin^2\theta$ 的最大偏差已在 2% 以内，且随 $N$ 二阶收敛——亲手跑一遍这段 50 行的小代码，对"离散收敛"的体感胜过十页讲义。

在真实世界里：势流至今活在哪里

不要以为势流是古董。机翼设计的第一轮快速扫描（上千个外形的气动筛选）仍用二维面板法 + 边界层耦合（XFOIL 的内核）；帆船船体、风电机组的初步水动力/气动力估算用三维面元法；甚至心脏内的血流初步分析也有势流近似的身影。它的价值在于毫秒级的评估速度：先把设计空间粗筛掉 95%，剩下的 5% 交给 RANS 精算——多保真度设计的金字塔，塔基就是势流。

常见误区与辨析

• "无旋 ⟹ 无粘"：因果反了。是无粘 + 从静止启动（开尔文定理）⟹ 保持无旋；反之无旋流放回粘性流体也只是局部近似。
• "势流解出负压力系数就错了"：$C_p < 0$ 完全合法（低于来流静压），负到 $C_p < -\frac{1}{2}\gamma M_\infty^2$ 附近才会触及"汽化/空化"或超音速的物理红线。
• "面板法解不收敛"：多半是控制点与节点重合导致自诱导奇异处理不当，或角点处法向不唯一——离散细节问题，不是方法问题。

与其他节的联系

本节的椭圆型方程（拉普拉斯）与第 6 节压力泊松方程同族——它们的数值处理（全场耦合、迭代解法）一脉相承；而"无粘模型的失效点"直接通向第 4 节有限差分与第 5 节的稳定性问题：一旦引入对流项（双曲性），数值方法的性格就完全变了。

专题深入：达朗贝尔佯谬的三种"修法"与边界层的诞生

达朗贝尔佯谬（理想流体中物体零阻力）的三种修法构成流体力学的一段思想史。修法一（库塔-茹科夫斯基）：加环量产生升力——但阻力仍为零，佯谬没全解；修法二（普朗特边界层）：粘性只在薄层里起作用，物面摩擦与尾流压差给出阻力——佯谬解决；修法三（湍流）：高雷诺数下流动的非定常性（涡脱落的平均阻力）——现代理解。佯谬的解决不是"加粘性项"这么粗暴，而是"粘性在局部扮演全局角色"——边界层与尾流的微小粘性效应决定了全局阻力。这段历史是"模型假设的局部失效产生全局后果"的教科书标本。

数学细节展开：保角变换——二维势流的"万能地图"

二维势流的终极工具是保角变换：解析函数 $w = f(z)$ 把复杂边界映成简单边界（圆柱 → 翼型），而拉普拉斯方程在共形变换下不变——复杂区域的势流 = 简单区域的势流 + 坐标变换。茹科夫斯基变换 $\zeta = \frac12(z + a^2/z)$ 把圆映成翼型是这套工具的巅峰：圆的绕流解直接翻译成翼型的升力、压力分布与失速特性。"用复变函数解流体力学"是 19 世纪末数学物理的辉煌时刻——它在数值方法兴起前是翼型设计的主力，今天仍是理解升力理论的最短路径。

工程案例细节：XFOIL 的工作流解剖

XFOIL（二维翼型分析的事实标准）的工作流值得解剖：面板法（势流解——毫秒级）+ 边界层积分（粘性修正——位移厚度反馈给势流）+ 转捩预测（e^N 方法）。三步循环迭代——势流给边界层提供边缘速度、边界层给势流提供修正后的等效外形——粘性-无粘耦合迭代是这个方法的全部智慧。它在一秒内给出翼型的升阻极曲线，是机翼初步设计的第一工具。理解它的架构，就理解了"快速评估"在工程设计中的位置：先把设计空间粗筛掉 95%，剩下的交给 RANS 精算。

常见误区与辨析（补充）

• "升力来自伯努利原理（上表面流速快压力低）"：伯努利是记账不是因果——升力的因果链是库塔条件 → 粘性产生起动涡 → 总环量守恒 → 附着环量 → 库塔-茹科夫斯基。流速差与压力差都是同一物理的两面。
• "面板法能算粘性效应"：不能——粘性修正靠边界层耦合（位移厚度反馈），不在面板法内部。