重现第一架隐形飞机背后的数学原理

# 重制夜鹰1号——隐形物理学 来源：https://www.pramit.gg/post/remaking-echo-1-stealth-physics 苏联数学如何促成了第一架经验最优化的隐形飞机 *事先声明，我绝非该领域的专家；我只是一名普通的软件工程师，但对航空航天和物理学有着超出平均水平的兴趣。所有可视化内容均由我绘制/生成，代码可在本项目GitHub仓库中找到 (https://github.com/pmazumder3927/nighthawk_rcs/tree/main/visualizations)* ## 黑色菱形与反弹的光子 1752308995975_bzq2ct.png 自从我开始钻研工程学以来，我就被那架可以说是最迷人的飞行器——SR-71“黑鸟”所吸引。仅这一架飞机的技术就令人难以置信，如果你不熟悉它，我强烈推荐Real Engineering在YouTube上关于这个主题的视频 (https://www.youtube.com/watch?v=3hYSnyVLmGE&pp=ygUWcmVhbCBlbmdpbmVlcmluZyBzciA3MQ%3D%3D)。 SR-71确实是一架隐身的飞机，但并非以我们今天理解的那种方式。它采用了巧妙的造型和吸波涂料，但其真正的防御手段是速度和高度。然而，后来的飞机，即F-117，则力求完全隐形。为此，它采用了一种完全不同的技术，这种技术由物理学而非工程学开创。 （另注：也强烈推荐关于夜鹰的视频 (https://www.youtube.com/watch?v=1zjcnnx7igc&t=555s)，其中提供了一些更直观的可视化内容） ### 苏联的论文 1971年，彼得·乌菲姆采夫发表了一篇论文¹ (https://www.pramit.gg/post/remaking-echo-1-stealth-physics#user-content-fn-1)，描述了电磁波如何在边缘上散射。这个想法在很大程度上被西方工程师忽视了，他们更关注实际应用而非衍射理论。然而，洛克希德·马丁“臭鼬工厂”的一位工程师丹尼斯·奥弗霍尔瑟意识到，这一理论可以用来制造一架能够最小化雷达散射截面（RCS）的飞机，正是这一概念最终促成了F-117的发展² (https://www.pramit.gg/post/remaking-echo-1-stealth-physics#user-content-fn-2)。 在这篇文章中，我将尝试对隐形的数学原理进行*足够*直观的解读，为下一篇博文中的代码实现做准备。 ### 隐形的物理学 要理解隐形，必须先理解雷达。RCS非正式地指*将雷达信号反射回源的有效面积*。这是什么意思呢？ 假设你有一个金属球，并向它发射雷达波，一部分能量会返回给你，大致随球体面积的增大而增大。大的RCS意味着雷达看到一个明亮的大目标，而小的RCS则产生一个幽灵般的回波。为了让你直观感受RCS数值，并强调这项技术所取得的巨大飞跃，我制作了一张小图表： | 对象/飞机 | RCS (m²) | RCS (dBsm) | 备注 | | :--- | :--- | :--- | :--- | | **小鸟³ (https://www.pramit.gg/post/remaking-echo-1-stealth-physics#user-content-fn-3)** | 0.01 | –20 dBsm | 难以探测 | | **人⁴ (https://www.pramit.gg/post/remaking-echo-1-stealth-physics#user-content-fn-4)** | 1 | 0 dBsm | RCS的基准 | | **F-15“鹰” (1972年)** | ~12 | ~11 dBsm | 传统战斗机 | | **SR-71“黑鸟” (1964年)⁵ (https://www.pramit.gg/post/remaking-echo-1-stealth-physics#user-content-fn-5)** | ~10 | ~10 dBsm | 有一定被动RCS降低 | | **F-117“夜鹰” (1981年)⁶ (https://www.pramit.gg/post/remaking-echo-1-stealth-physics#user-content-fn-6)** | 0.01 – 0.1 | –20 到 –10 dBsm | 拓扑优化 | 现在，让我们深入了解更正式的定义。RCS由以下方程定义： σ = lim(r→∞) 4πr² |E_s / E_i|² 其中 R 是到目标的距离（想象雷达非常远）(en.wikipedia.org)。这个公式可能看起来很吓人，所以我们来分解一下： - 当 R→∞ 时，雷达位于远场，此时电磁波可视为射线。因子 R² 解释了这样一个事实：随着距离增加，散射波在面积为 4πR² 的球面上扩散。将散射场强度 |E⃗_s|² 乘以 R² 本质上是测量该方向上的总散射功率。 - 比值 |E⃗_s|² / |E⃗_i|² 比较了散射波与入射波的强度。如果我们的对象是一个横截面积为 1 m² 的理想反射球体，那么该比值（乘以 4πR²）在精确的背向方向上会等于 1。 - 因此，σ 表现为一个有效面积（以 m² 为单位），代表目标将雷达能量反射回来的强度。如果 σ = 1 m²，目标就像是一个 1 m² 的金属球一样可探测。如果 σ = 0.001 m²，那么雷达看到的就像是一个弹珠大小的东西（是 B-2 有效截面的 10 倍）！ 简而言之，把目标包裹在一个巨大的气泡里，测量通过单位面积反弹回来的功率，用气泡面积 (4πr²) 进行缩放，再与你最初打击的强度相比较——有效面积就出来了。 从这个方程我们还可以看出，即使是很小的改变对探测效果也有巨大影响；探测距离大致与 RCS 的四次方根成正比，因此将 RCS 降低 100 倍，探测距离会降至原来的 1/100 的四次方根 ≈ 56%。那么我们是如何做到这一点的呢？ ### 计算 RCS 要经验性地计算像飞机这样的复杂形状的 RCS，精确解需要求解带有适当边界条件的麦克斯韦方程组；这是因为本质上要求解电磁波如何在形状上散射。这极具挑战性，因为你需要： - 将飞机的表面离散化为 N 个面片 - 建立一个方程组来模拟**所有面片对**之间的相互作用，以计算总散射场。这一效应由格林函数（一个 (3N × 3N) 的线性系统）决定。对于一个有 1000 个顶点的网格 (N=1000)，格林函数评估的规模达到了 3000 * 3000 个方程。在 20 世纪 70 年代，仅仅存储这些方程就需要大约 100MB，更不用说*求解*它们了；当时最好的计算机只有 **1MB** 的内存。 ### 物理光学 (PO) 登场 正如许多伟大工程技巧的诞生一样，物理光学也源于一个简单的技巧：简化假设。 - **表面为理想电导体 (PEC)：** 如果假设表面是理想导体，则表面上的切向电场为零，因此不再需要建模材料特性，并且可以将电磁作用完全视为表面上的积分。 1752312039402_0227on.png - **仅考虑被照射（正面朝向）的面：** 如果一个面片不在雷达波直接照射的那一侧，则完全忽略它。这使计算量减半（或更多），并且不再需要建模波“爬行”到物体另一侧的复杂相互作用或反射。 1752310428680_qyt61x.png - **场通过镜面反射散射：** 反射波总是以与入射角相等的角度从表面反弹，就像镜子一样。这样，你永远不需要建模多次弹跳、衍射、散射或其他任何杂乱的东西。 所有这些假设都牺牲了准确性。但它们极大地简化了计算，而且事实证明（至少在当时），这些简化模型足以大幅降低 RCS。现在，让我们进入数学部分，看看这些假设如何为我们工作： 对于理想导电表面，在 PO 下，感应表面电流密度（以向量形式）由下式给出： J⃗_s(r⃗) = 2 n̂(r⃗) × H⃗_i(r⃗), 其中 n̂(r⃗) 是表面点 r⃗ 处指向外的单位法向量，H⃗_i 是该点的入射磁场。这个公式基本上是说，表面电流的流动方式是为了抵消表面的切向磁场，因子 2 来自于反射场在表面处与入射场的叠加。 现在我们知道表面电流是什么了，就可以通过简单地积分这些电流在表面 S 上的分布，来得到任意方向 k̂_s 上的散射场： E⃗_s(k̂_s) = (j k η) / (4π) ∫_{S_lit} [ k̂_s × ( k̂_s × J⃗_s(r⃗) ) ] e^{j k (k̂_i - k̂_s) · r⃗} dS. 这可能是物理光学中最关键的方程，所以我们来逐项分析： - k = 2π/λ 是波数（λ 为雷达波长）。例如，X 波段雷达在 10 GHz 时，λ ≈ 3 cm，所以 k ≈ 210 m⁻¹。η ≈ 377 Ω 是自由空间的阻抗，是一个联系空气中电场与磁场的常数。k̂_i 是入射波的方向（从雷达到目标），k̂_s 是散射波的方向（从目标返回雷达）。对于单站 RCS（我们通常关心的情况），k̂_s 正好是 k̂_i 的反方向（因为雷达从同一位置发射和接收）。 - 被积函数 k̂_s × (k̂_s × J⃗_s) 是一个向量，它本质上挑选出向雷达方向辐射的表面电流分量。（如果你从天线理论或偶极子物理学中回想起来，电流元会垂直于其方向辐射；这个双叉积有助于对向接收器再辐射场的正确极化进行强制）。 - 指数项 e^{j k (k̂_i - k̂_s) · r⃗} 是一个相位因子。它解释了路径差：如果表面的一块区域相对于另一块区域稍微更远或角度不同，它们的再辐射波将会相位不同地到达。当我们进行积分时，这些相位差会导致干涉，从而在某些方向上抵消或放大散射波。这正如同将许多小的相量（旋转复数的微小箭头）相加，每个小相量对应一个表面面片。积分仅在表面的被照射部分 S_lit 上进行。如果物体的某一部分背对雷达（在阴影中），PO 假设它没有感应电流（因为入射波没有直接照射它）。这是一种近似，因为实际上阴影区域仍然可以通过边缘附近的衍射进行散射，但 PO 忽略了这一点（这是其已知的局限性之一）。 这种方法在频率高且物体尺寸远大于 λ 时快速且相当准确。上述局限性（边缘衍射、爬行波）可能很显著，最初的“Echo 1”程序在认识到其重要性后，又添加了边缘衍射建模。然而在我的项目中，我坚持使用基本的 PO 模型，它应该足以捕捉大部分主要效应。 1752312201386_d3vzw4.png *雷达视角偏离正面的角度越大，PO 的准确性就越低。幸运的是，对于飞机优化而言，这通常是可以接受的。* ### 为什么是刻面？ 1752355473770_quokq2.jpg 此时一个自然的问题是，为什么这架飞机如此棱角分明。需要指出的是，来自 PO 原理优化的一个核心见解是：平坦且有角度的表面非常适合隐身。让我们稍微展开说明： - 记住，PO 的关键见解之一是，一个表面会像镜子反射光一样反射雷达（镜面反射）。因此，任何平坦的表面都可以通过**简单地将其倾斜远离源**，将雷达能量偏转开。 - 曲面：人们可能认为球体接近隐身的最优解（我确实这么认为过）。但如果你想象一个迪斯科球散射光线，**无论你身处何处，曲面的某一部分总会将光线反射回给你**。这正是迪斯科球会闪烁出标志性光芒的原因，雷达上也会发生完全相同的情况。因此，对于隐身而言，球体实际上是最**糟糕的情况**。 - **内部反射**——如果你想象一个完全最优的雷达目标，它会将所有能量完全反射离开来自源的任何角度。因此，在微观层面上，它应该看起来像一组平面镜，全部倾斜远离源（理想情况下是垂直的）。想象一下就像迪斯科球类比中那样用手电筒照射一面镜子。如果镜子倾斜远离你，你发出的光线中不会有太多返回。 1752311938819_bzcici.png *平坦表面产生对齐的电流和可控的散射；弯曲表面产生变化的电流并增加反向散射。* 对此也有具体的数学观点：如果我们观察上面关于 E⃗_s 的积分，可以看出 k̂_i - k̂_s 在一个刻面上是线性变化的，但在曲面上则是非线性变化的。局部曲率会导致反向散射增加，并将能量散射到不需要的方向。 ### 权衡取舍 平坦且有角度的表面对于需要飞行的东西来说并不特别好；没有汽车或飞机被设计成矩形是有原因的。F-117 由于其奇怪的形状以及 inherently unstable (本身是不稳定的)⁷ (https://www.pramit.gg/post/remaking-echo-1-stealth-physics#user-content-fn-7) 而被戏称为“摇摆的地精”，需要计算机系统才能保持在空中。 ### 下一篇 现在我们对数学原理有了较好的理解，我将带你重新创建一个我本周末制作的类似 Echo 1 的程序，使用现代库和一些计算技巧。不过这篇文章已经有点长了，我也挺累了，所以留到下一篇新文章吧。 希望你喜欢！ ### 脚注 1. P. Ya. Ufimtsev, *Method of Edge Waves in the Physical Theory of Diffraction*. FTD 翻译 HT-23-180-71 (1971).↩ (https://www.pramit.gg/post/remaking-echo-1-stealth-physics#user-content-fnref-1) 2. D. Jenkins, *Flying the Black Triangle*, AIAA 2018 — 概述 F-117 开发，引用 Overholser 的见解。↩ (https://www.pramit.gg/post/remaking-echo-1-stealth-physics#user-content-fnref-2) 3. E. Knott, J. Schaeffer, M. Tuley, *Radar Cross Section*, 2e, Artech 2004.↩ (https://www.pramit.gg/post/remaking-echo-1-stealth-physics#user-content-fnref-3) 4. FAA, *Radar Handbook*, 第3章 — 人体目标截面。↩ (https://www.pramit.gg/post/remaking-echo-1-stealth-physics#user-content-fnref-4) 5. B. Rich & L. Janos, *Skunk Works*, Little Brown 1994, 第12章 (SR-71 的 RCS 数据, F-117 开发轶事)。↩ (https://www.pramit.gg/post/remaking-echo-1-stealth-physics#user-content-fnref-5) 6. 美国空军国家博物馆，F-117 资料表 (公开的 RCS 范围估计)。↩ (https://www.pramit.gg/post/remaking-echo-1-stealth-physics#user-content-fnref-6) 7. J. Correll, “History of Stealth,” *Air & Space Forces Magazine*, 2019年9月 — 计算限制、刻面化、边缘衍射修正。↩ (https://www.pramit.gg/post/remaking-echo-1-stealth-physics#user-content-fnref-7)