4.3 位势流基础

4.3 位势流基础

本节摘要：当流动无旋（流体微团不自转）且不可压时，速度场可以用一个标量函数——速度势——的梯度表示，整套运动学塌缩成一个拉普拉斯方程。势流方法的回报巨大：伯努利方程升级为全场通行，叠加原理成立，绕流问题从此可以"拼积木"。代价是粘性与分离被请出场外——这层欠账由第 5 章偿还。

前两节靠几何极端简单才求出精确解。势流方法换一条路：不追求方程全解，先在运动学上做减法。

无旋与速度势：减法的起点

流体微团在运动中可能自转（想象顺流而下的树叶打着旋）。数学上用旋度度量每一点的自转强度。若流场中处处旋度为零，称无旋流动。无旋是个强约束，它带来一个漂亮结论：存在标量函数 φ，使速度等于 φ 的梯度（速度沿哪个方向、变化多快，都由 φ 的坡度给出）。这个 φ 叫速度势——"势"字的含义与重力势能同源：顺着势的坡往下走就是顺着流线走。

再叠加不可压（速度散度为零），立刻得到速度势满足拉普拉斯方程——数学物理中研究最透的方程之一，调和函数的一切性质（极值原理、唯一性、叠加性）全部继承。运动学的两条约束（无旋 + 不可压）把求解流场变成了求一个调和函数，动力学问题暂时退场，这就是势流方法轻快的原因。

配一个镜像工具：流函数 ψ。不可压二维流动自动存在流函数，它有两个用途——ψ 相等的曲线就是流线；任意两条流线之间 ψ 的差值等于其间通过的体积流量。速度势管"快慢"，流函数管"路径"，两者正交网格铺满流场，是可视化与解析的共同底座。

伯努利的全场通行版

第 2 章的伯努利只沿流线成立，跨流线要额外条件。势流把条件补齐了：无旋流动中，p + ½ρV² + ρgz = 常数 在整个连通的流场内成立。这一步升级的工程价值立竿见影：只要解出速度场（拉普拉斯方程，线性好解），全场压强分布一次伯努利全给——升力正是从压强分布里长出来的（第 6 章）。

无旋条件的实际合理性也要说清：远离物体的主流区通常均匀来流，天然无旋；粘性在物面上产生的旋量，在雷诺数较大时被限制在边界层内（第 5 章的核心图像）。所以"主流无旋 + 边界层有旋"是一幅自洽的合作分工图，势流负责前者。

势流的适用体检表

飞机巡航时前三条都过；直升机桨-涡干扰、失速后的机翼则全线告警。用势流前先做体检，跟用任何简化模型一样。

动手算：圆柱表面最高流速

势流解出的圆柱绕流（下一节拼装）在表面给出的速度为 u = 2V·sinθ（V 为来流速度，θ 从前驻点量起）。取来流 30 m/s：最侧面 θ = 90° 处表面流速 60 m/s——两倍来流。用全场伯努利算那里表压：p − p∞ = ½ρ(V² − u²) = ½×1.225×(900 − 3600) ≈ −1654 Pa，比来流静压低约 1.6 kPa。两倍来流速度、四倍动压的表面吸力峰值，这个"2 倍"系数在所有钝体与翼型上反复出现（翼型最厚处附近表面流速同样可达来流的 1.5~2 倍），值得作为量级直觉收藏。

流函数的实用三招

流函数 ψ 不只是对称美的装饰，工程上有三招直接用法。第一招：画流线不用解方程——给定 ψ 的表达式，令其等于一串常数，曲线族即流线，可视化一步完成。第二招：流量心算——两条流线之间的体积流量就是 ψ 值之差乘以单位厚度，做通道设计的估算时比积分快得多。第三招：构造满足壁面条件的解——把物面设为某条 ψ 等值线，自动保证无穿透；镜像法（把物体对地面做个"镜像源"，对称面的 ψ 为常数）就是这一招的直接应用，地面效应、风洞壁面修正都靠它。

速度势与流函数在不可压二维无旋流动里满足柯西-黎曼关系，构成正交复势——这门"复分析走进流体"的老手艺至今仍是快速估算绕流的利器：单位阶跃、保角变换，把复杂外形"烫"成圆柱来解。

问题：真实流动怎么可能严格无旋？

的确不可能处处无旋——壁面边界层里旋度高度集中。势流的合法性来自"旋度在哪里"而不是"有没有"：高雷诺数下旋度被粘性锁在薄层与尾流里，主流区基本干净，势流解在主流区的误差很小。对升力问题（压差主导、分离不严重），这套"主流势流 + 薄层修正"的组合拳精度惊人；对大分离流动（失速、钝体），则必须换方法。判断旋度的"领土面积"，就是 4.3 节那张体检表的物理版。

问题：伯努利"全场通行"会不会跟"沿流线"版本矛盾？

不矛盾，只是通行条件升级。沿流线版本不要求无旋；全场通行版本要求无旋（或更弱的"机械能沿流线守恒且各流线同值"，即无旋流场）。换句话说全场版本是更苛刻条件下的特例，特例给出了更强的结论——跨流线也能比压强。翼型压力分布计算之所以能一步从速度场跳到全场压强，用的正是这个加强版。

问题：三维流动还有流函数吗？

部分有。三维轴对称流动（火箭、旋成体）存在斯托克斯流函数，两条流面之间的流量同样守恒；一般三维流动则没有单个标量能同时管流线与流量，速度势（若无旋）仍可用。这就是为什么势流方法在二维玩具问题里无所不能、在三维工程问题上要靠面元法把物体表面离散成无数“局部二维”来拼——降维是三维势流的生存策略。

问题：速度势与流函数哪个更常用于工程？

看任务。求压强与载荷用速度势（伯努利要吃速度模长）；画流动图案、定流量分配用流函数。二维翼型的经典手算几乎都在流函数与复势里进行，三维面元法则全部工作在速度势上。两者由柯西-黎曼关系绑定，知道一个即可恢复另一个——工程上更像“同一工具的两个握法”，而不是两门学问。

本节要点

• 无旋 + 不可压 → 速度势满足拉普拉斯方程，运动学与动力学解耦。
• 流函数：等值线即流线，差值即流量，与速度势构成正交网格。
• 伯努利全场版只在无旋时通行，势流因此能一步从速度场跳到压强场。
• 适用体检：低马赫、高雷诺、薄边界层三条件，缺一先三思。
• 圆柱侧点流速两倍于来流、吸力四倍动压，是绕流的通用量级锚点。

运动学工具齐备。下一节请出四种"基本解"，示范积木拼装——绕流场的建造术。