本文是《复变函数完全入门》系列的第四篇，覆盖保角变换的定义、性质、儒科夫斯基变换及其在翼型空气动力学中的应用。阅读本文前需掌握解析函数的定义和CR方程。

第四篇：保角变换——从圆到翼型的数学桥梁

一、概念引入——为什么需要保角变换

1.1 问题

势流理论中，圆柱绕流已经有了精确解析解。但工程师关心的是翼型绕流，不是圆柱绕流。

有没有一种方法，能把圆柱绕流”变形”成翼型绕流？

有——保角变换。

核心思想：

在 $\zeta$ 平面（计算平面）上，我们有一个圆——绕流解已知
通过一个解析函数 $z = f(\zeta)$ ，把圆映射为翼型形状
圆柱绕流解同步映射——就得到了翼型绕流解

1.2 历史的五分钟

时间	人物	贡献
1851	黎曼	黎曼映射定理：任何单连通区域都可以保角映射到单位圆
1910	儒科夫斯基（Zhukovsky）	提出儒科夫斯基变换——从圆到翼型的精确映射
1922	库塔（Kutta）	应用儒科夫斯基变换计算翼型升力
1933	西奥多森（Theodorsen）	提出一般翼型的数值保角变换方法

二、保角变换的定义

2.1 定义

一个复变函数 $w = f(z)$ 如果在 $z_0$ 处满足 $f'(z_0) \neq 0$ ，则它在 $z_0$ 处是保角的：

保角性（Angle-preserving）：经过 $f$ 的映射，任意两条曲线在 $z_0$ 处的夹角保持不变
膨胀率：任意方向的伸缩因子都是 $|f'(z_0)|$ ——即各向同性（不因方向而异）

2.2 性质的重要性

保角性的物理含义：经过保角变换后，流线和等势线的夹角仍为90°——势流解在变换前后保持一致。

特性表：

性质	含义	在空气动力学中的意义
保角性	角度不变	流线与等势线的垂直关系不变
各向同性	伸缩因子仅与位置有关	流场结构在变换后保持”形状”
解析函数	变换函数必须解析	只有CR方程满足的变换才可用
$f'(z) \neq 0$	非共形点需特殊处理	翼型后缘就是非共形点→库塔条件

三、儒科夫斯基变换——最重要的保角变换

3.1 定义

儒科夫斯基变换：

\boxed{z = \zeta + \frac{c^2}{\zeta}}

其中 $c$ 是一个实数常数，决定了变换的特性。

3.2 变换的几何效果

令 $\zeta = re^{i\theta}$ ，代入：

z = re^{i\theta} + \frac{c^2}{r}e^{-i\theta} = \left(r+\frac{c^2}{r}\right)\cos\theta + i\left(r-\frac{c^2}{r}\right)\sin\theta

所以：

x = \left(r+\frac{c^2}{r}\right)\cos\theta, \quad y = \left(r-\frac{c^2}{r}\right)\sin\theta

对 $\zeta$ 平面上的一个圆 $r = R$ ，它的像是什么？

如果 $R = c$ ：圆映射为线段（在实轴上从 $-2c$ 到 $2c$ ）
如果 $R > c$ ：圆映射为一个封闭曲线——儒科夫斯基翼型
如果 $R < c$ ：圆映射为另一个封闭曲线（有自交）

3.3 翼型的形成

取 $R > c$ 的圆，圆心在 $(-\mu, 0)$ 处：

$R - c$ 决定了翼型的**厚度**
$\mu$ 决定了翼型的**弯度**
圆上的 $\theta = 0$ 点映射到翼型的后缘

前缘：由圆上 $\theta \approx \pi$ 附近的点映射而来，曲率半径正比于 $(R-c)^2/(R+c)$

后缘角：后缘处的夹角为零——尖锐后缘正是库塔条件要求的几何形式。

四、翼型的升力系数

4.1 库塔条件在后缘的应用

在 $\zeta$ 平面上，我们在圆上加一个点涡（环量 $\Gamma$ ），使得在映射后的翼型后缘处满足库塔条件——速度有限。

由此可以确定环量：

\Gamma = 4\pi V_\infty R\sin(\alpha + \beta)

其中 $\beta$ 与圆心偏移 $\mu$ 相关。

4.2 升力系数

由库塔-儒科夫斯基定理 $L' = \rho V_\infty \Gamma$ ：

L' = 4\pi\rho V_\infty^2 R\sin(\alpha + \beta)

代入升力系数公式 $C_l = L'/(\frac12\rho V_\infty^2 c)$ ：

\boxed{C_l = 2\pi\left[\sin(\alpha + \beta) + \sin\beta\right]}

对于小迎角，简化为：

C_l \approx 2\pi(\alpha + 2\beta) = 2\pi(\alpha - \alpha_{L=0})

这正是薄翼理论结果—— $2\pi$ 斜率在此由保角变换自然导出。

五、其他重要保角变换

5.1 幂函数变换 $w = z^n$

$n=2$ ：把上半平面映射为全平面切除正实轴
$n=\frac12$ ：把角形区域映射为半平面

在势流理论中， $w = z^{1/2}$ 可用于描述90°拐角处的流动。

5.2 对数变换 $w = \ln z$

把环形区域映射为矩形区域。这也是绕流中从物理平面到复势平面的映射。

5.3 Schwarz-Christoffel 变换

可以将多边形内部保角映射到上半平面。这是数值保角变换的基础——用于计算任意形状翼型的绕流。

六、完整概念地图

复变函数理论 → 空气动力学应用

解析函数 f(z)=φ+iψ
    └── 实部=速度势, 虚部=流函数

保角变换 w=f(z)
    │
    ├── 黎曼映射定理: 任何单连通区域→圆
    ├── 儒科夫斯基变换: z=ζ+c²/ζ
    │   ├── 圆→翼型
    │   ├── 环量确定→库塔条件
    │   └── Cl=2π(α-αL=0) 自然导出
    │
    ├── 幂函数: zⁿ → 角形区域
    ├── 对数: ln z → 环形→矩形
    └── Schwarz-Christoffel → 多边形→半平面

七、核心公式速查卡

公式	含义
$z = \zeta + c^2/\zeta$	儒科夫斯基变换
$f'(z) \neq 0 \Rightarrow$ 保角	保角变换条件
$\Gamma = 4\pi V_\infty R\sin(\alpha+\beta)$	库塔条件确定的环量
$C_l = 2\pi(\alpha - \alpha_{L=0})$	翼型升力系数的复变推导

参考文献

Brown, J. W., & Churchill, R. V. (2013). Complex Variables and Applications (9th ed., Chapter 8). McGraw-Hill.
Anderson, J. D. (2010). Fundamentals of Aerodynamics (6th ed., Chapter 4). McGraw-Hill.
Joukowski, N. E. (1910). “Über die Konturen der Tragflächen der Drachenflieger”. Zeitschrift für Flugtechnik und Motorluftschiffahrt, 1, 281-284.
Milne-Thomson, L. M. (1996). Theoretical Hydrodynamics (5th ed., Chapter 9). Dover.
Needham, T. (1997). Visual Complex Analysis (Chapter 8). Oxford University Press.

下一节：复变函数在空气动力学中的应用

用复势统一描述源、汇、偶极子、点涡；构造均匀流+偶极子+点涡的圆柱绕流解；用儒科夫斯基变换得到翼型绕流的精确解——儒科夫斯基变换的全部数学基础。