简介：这份来自北京航空航天大学精品课程的空气动力学基础教学课件，以PDF格式呈现，共1个文件，压缩包大小19.65MB，适合航空航天类专业学生、教师及相关工程技术人员系统学习与参考。课件内容涵盖绪论、流体的基本属性、流体静力学、流体运动学、不可压缩无粘流体平面势流、粘性流体动力学基础、边界层理论及其近似、可压缩高速流动基础等核心章节，并对空气动力学的研究方法、分类、发展进程及应用场景做了系统梳理，结构完整，便于按章节循序渐进地学习。目前已有675人学习使用。学习者可通过本课件建立空气动力学整体知识框架，理解飞行器绕流现象与受力分析的基本原理，为后续从事航空航天工程实践或深入研究打下扎实基础。
1. 一份北航精品课程PDF，凭什么值得从头啃一遍
拿到《空气动力学基础(北航精品课程)-刘沛清.pdf》，很多人第一反应是“又是一堆公式”。但我见过太多做CFD的人，能跑通算例却说不清升力从哪来，残差曲线一乱只能靠调松弛因子碰运气——这是基础没补好。这门课是北航飞行器类专业的核心课，刘沛清老师把流体力学筛了一遍，留下的都是直接服务于飞行器受力的内容。它适合刚入行的气动仿真工程师、在考研或读研的学生、做总体设计需要估算气动数据的人。它不教你花哨技巧，每一页都在回答“受力从哪来、什么条件下成立、算错了会怎样”。这篇笔记写的是我怎么用它搭起一套可复现的学习框架，以及踩过的坑。
2. 空气动力学基础的骨架：控制方程与无量纲参数
这份PDF的目录通常先讲大气环境、流体属性，然后进入流动控制方程，再按“不可压缩、可压缩、粘性”三大块展开。如果你跳着读，很容易被符号淹没。我的建议是先抓骨架：控制方程决定物理模型，无量纲参数决定相似性，流动分界决定简化方向。这三件事不立住，后面读翼型、读边界层、读激波都会像是在背字典。
2.1 从NS方程到欧拉方程：简化条件决定适用范围
N-S方程是空气动力学的地板，但直接从三维非定常N-S方程推起会劝退很多人。刘沛清老师的讲义做法是先写出守恒形式，然后一步步删除项，告诉你每一项在什么条件下可以忽略。比如不可压缩、无粘、无热传导情况下，N-S方程退化为欧拉方程；再叠加无旋条件，变成速度势方程，这就是全位势方程的来源。这里要记住的不是方程的形式，而是“删掉哪一项，对应什么假设”。
我一般看N-S方程会先看它的四个来源：质量守恒、动量守恒（三个方向）、能量守恒。低速、无热源流动中，能量方程可以和动量方程解耦，剩下连续性和动量方程就够用了。不可压缩时密度视为常数，连续性方程里密度对时间的偏导直接清零；动量方程里的粘性应力项在边界层外很小，但在物体表面附近又必须保留。所以同一个流动在不同位置可以用不同方程，这是初学者最容易困惑的地方——为什么“无粘”还能描述真实机翼外面的流场？因为粘性的作用局限在边界层和无粘区域边界的相互作用里。
在做实际工程估算时，我会先把流动归入某一类：低速不可压无粘、低速有粘、高速无粘、高速有粘带激波。这个分类直接决定选什么公式。比如估算翼型升力用势流理论足够，但估算阻力必须回到边界层；估算超音速激波位置要用斜激波关系式，而激波后的边界层问题又得回到N-S。PDF里大量篇幅在讲这些简化条件，背后的目的就是让你拿到一个问题时，能判断该用哪一层工具。
| 方程/模型 | 速度适用范围 | 忽略内容 | 典型应用 |
|---|---|---|---|
| 拉普拉斯方程 | 低速 Ma<0.3 无粘 | 粘性、旋涡、压缩性 | 势流法求升力、压力分布 |
| Euler方程 | 跨音速/超音速无粘 | 粘性、热传导 | 激波位置、波阻估算 |
| 边界层方程 | 近壁区粘性主导 | 法向压力梯度近似为零 | 摩擦阻力、转捩判断 |
| 完整N-S方程 | 全速域 | 无 | CFD底层、分离流/湍流直接模拟 |
这个表是我从PDF的公式脉络里抽出来的。它不是为了让你背方程，而是告诉你：当你跑CFD时，求解器内部默认用哪一套方程，你就能解释为什么某些工况会发散。比如在分离流大攻角下，有些基于欧拉方程加边界层耦合的工具会给出乐观的升力，因为把分离后的压力分布算得太“干净”了。
2.2 三个绕不开的无量纲参数：雷诺数、马赫数、升力系数
无量纲参数是空气动力学从“玄学”变科学的底气。三个参数必须形成条件反射：雷诺数决定流动状态，马赫数决定压缩性，升力系数把形状和工况解耦。刘老师的讲义里所有实验数据都以这三个参数为主轴，你只要摸清它们的量级和物理意义，就能看懂大多数图表坐标轴。
雷诺数定义为Re = ρVL/μ，其中ρ是密度，V是特征速度，L是特征长度，μ是动力粘度。它表示惯性力与粘性力的比值。Re低的时候流动是“粘着的”，容易层流；Re高的惯性力主导，扰动容易放大成湍流。实际工程里，小型无人机翼弦0.2m、速度15m/s，Re大概20万左右；民航客机翼弦3m、速度250m/s，Re超过5000万。两个Re差了三个数量级，同样一块翼型的最大升力系数和失速特性完全不同。很多风洞数据在Re百万以下测得，直接用到飞机上必然翻车。
马赫数Ma = V/a，a是当地声速。Ma<0.3可以认为密度不变，这是不可压缩近似的标准门槛。从0.3到0.8是亚跨音速，翼面上出现局部超声速区但主流仍是亚声速，必须考虑压缩性修正；Ma>1是超声速，存在激波。课程里对Ma的讨论不只是分类，还会给出一套等熵关系式：T0/T = 1 + (γ-1)/2 · Ma²，p0/p = (T0/T)^(γ/(γ-1))。这两个公式在计算风洞试验条件、推算飞行器驻点温度时特别常用。比如飞行器Ma=2时，驻点温度按γ=1.4算，大约是静温的1.8倍，如果不做热防护设计，结构强度会直接出问题。
升力系数Cl = L/(q∞·S)，其中q∞是动压，S是参考面积。它把升力从速度中解耦出来，方便对比不同翼型。薄翼型线性范围内Cl≈2π(α-α0)，α0是零升攻角。这个线性关系在真实流体中只在一定攻角内成立，失速后Cl急降。学习时要把Cl-α曲线从讲义中抽出来，反复对照：附着流阶段斜率大约等于2π（单位弧度），非线性段是分离的开始，失速点是最大Cl对应的攻角。这三个阶段对应三种物理行为，也对应三种CFD收敛难度。
2.3 不可压缩与可压缩流动的分界在哪里
分界不是一条严格的线，而是一个工程决策点。通常以Ma=0.3为界。在这个Ma下，密度变化最多5%，工程上当作不可压缩处理，误差可以接受。超过0.3后，密度变化带来的流管面积和速度之间的耦合不能忽略，需要使用可压缩质量方程、动量方程和能量方程的联立形式。这个分界在PDF里反复出现，也在求解器里对应两套压力修正算法：不可压缩用压力泊松方程，通过压力修正保证连续；可压缩用密度基方法，直接解守恒变量。
在低速部分，课程会先引入流函数和速度势，建立无旋流的基本解：均匀流、点源、点涡、偶极子。叠加这些基本解，可以得到圆柱绕流的势流解，但没有升力——这就是著名的达朗贝尔佯谬。加入点涡后，通过库塔-茹可夫斯基定理L' = ρ∞·V∞·Γ得到升力，Γ是环量。升力从“环量”来，这是理解升力的正确入口，而不是单纯靠“上表面快下表面慢”那套直觉。
进入可压缩部分后，课程会从一维等熵流讲起，再讲正激波、斜激波和膨胀波。初学者最容易走神的地方是：总压和静压的区别。总压是把气流等熵滞止到零速后得到的压力，它是流动“有效做功能力”的量度；静压就是当地真实气体分子撞击壁面的压力。皮托管测的是总压，静压孔测的是静压。跨过激波时，总压下降，静压上升，熵增加。这些概念在飞机设计里直接决定发动机进气道的总压恢复系数，也决定气动加热程度。学习时不要只记公式，要在PDF的每个例题里标出：哪些量是滞止量，哪些是当地量。
3. 把PDF读成技能：按模块推进的学习路线与笔记方法
PDF课件和教材不一样。教材有完整的前后推导，课件往往把过程压到只剩关键结论。直接从头翻课件，很容易在导数、旋度、斯托克斯定理之间迷路。我的建议是，不要按页序读，把它当成一个“技能树”来铺设：先建立受力直觉，再按低速、可压缩、粘性三个模块推进，最后把每个结论整理成可检索的知识卡片。
3.1 第一步：建立坐标系、受力定义和速度三角
打开PDF第一遍，先找到翼型受力图，把坐标系、速度三角形、升力和阻力方向定义抄一遍。注意升力永远垂直于来流方向，而不是垂直于翼弦。很多初学者画图时随手一画，把升力画成和机翼表面垂直，后面所有力矩分析全部错位。正确做法是：来流速度V∞水平向右，攻角α是翼弦与来流的夹角，升力垂直V∞，阻力平行V∞，气动力矩以某参考点为正。这张图几乎决定后面所有公式的方向和符号。
然后去PDF里找“压力中心”和“焦点”的定义。压力中心是气动力合力与弦线的交点，它随攻角移动；焦点是气动力矩系数不随攻角变化的参考点，也叫气动中心。你做纵向静稳定性分析时，只需要看焦点相对于重心的位置：焦点在重心之后，扰动攻角产生的低头力矩会恢复原状，即静稳定；焦点在前，扰动会发散。这个判断不用算压力中心，算焦点就够了。
我建议你把这三个参考点（重心、焦点、压力中心）画在一张图上，标出典型攻角下各自位置，再用一个简单算例验证稳定性。比如一个常规布局飞机，焦点约在平均气动弦25%处，重心在30%处，那么焦点在重心前？实际上重心在前才是稳定。画图时别搞反。这一步过关了，后面读升力线斜率、下洗流、机翼-尾翼干扰才有坐标感。
3.2 第二步：按低速、可压缩、粘性三个模块推进
学习顺序不要按PDF目录从前到后闷头啃。我推荐倒着安排，因为低速无粘模块是基础，粘性模块解释阻力和分离，可压缩模块最后加宽带压缩性和激波。先低速无粘，再把无粘结果和真实结果的差距用边界层补上，最后处理高速效应。
低速无粘部分，重中之重是势流理论和环量。手推至少一遍圆柱绕流加点涡的叠加：均匀流加偶极子加点涡，通过速度势求表面速度，再用伯努利求压力，对压力沿圆柱表面积分，得到升力为ρ∞V∞Γ，阻力为零。这个推导串起速度势、流函数、叠加原理、伯努利、库塔-茹可夫斯基定理五件事。做完这一遍，你会突然看懂翼型升力线理论里的“平板替换”是怎么回事。
粘性部分，核心是边界层概念。刘沛清老师在讲义里会对边界层厚度、位移厚度、动量厚度做严格定义。位移厚度是边界层内因速度亏损而把外流“推开”的等效距离，它影响翼型压力分布；动量厚度直接进入阻力公式。边界层是否分离，主要看逆压梯度大小和雷诺数。一个常见误用是“只要Re大，边界层就薄，阻力就小”，实际上Re大往往是湍流边界层，摩擦阻力系数反而相对下降，但转捩位置提前，可能带来更大的分离趋势。这里要理解层流和湍流在速度分布上的差别。
可压缩部分，不要急于去背激波公式。先掌握一个一维变截面管流问题：缩放喷管里面积变化如何影响速度和压力。然后理解正激波的基本关系，比如马赫数越大，总压损失越大。用PDF里的表格或图，查一道“马赫2气流过正激波后马赫数约为0.577，静压比约为4.5”的题。这个数量关系能让你知道超音速进气道设计为什么难：损失太大。
每个模块学完，找一道例题，不看答案，按PDF给的步骤手算。手算不是要你数值精确，而是强迫自己停下来看单位和适用条件。算完后用讲义后面的参考数据对照，误差在一个数量级内即可。
3.3 第三步：用算例和可视化把公式变成直觉
公式背得再熟，不如看一次流场图。PDF里会有翼面压力分布图，升力系数曲线，也可能有油流或PIV实验照片。你要做的事情是：把公式结果和这些可视化对应起来。比如看到上表面吸力峰附近等压线密集，就应该想到那里速度高、压力低、局部静压最低。看到失速状态下吸力峰崩溃，就应该想到边界层分离、环量损失。
我常用的验证方法是做一个手算小表格：给定翼型弦长0.15m，速度30m/s，空气密度1.225kg/m³，粘度1.8e-5Pa·s，Re≈30.6万。按照PDF里的经验值查这个大攻角下可能的Cl范围，比如18万Re的小型无人机在攻角12°时Cl可能在1.1到1.3之间。然后手算升力L=Cl·0.5·ρ·V²·S，得到单位展向升力大约多少牛。这个数字再去和CFD结果对比，差得远了就说明计算域或湍流模型有问题。这个“手算预判”的习惯，能让你的数值模拟不再是黑匣子。
可视化方面，PDF里如果有实验影子图或纹影图，先看激波前后流线的折转角，再对照斜激波理论预测的折转角。一般偏差都来自粘性边界层对等效楔角的影响。别怕偏差，怕的是你完全没有预测值。
3.4 笔记怎么记：参数、公式、适用条件三栏卡片
我当年是吃过亏的：拿着PPT从头抄了一遍，合上文档什么也没留下。后来改成卡片法，每个知识点写一张卡片，分三栏：参数定义与单位、公式与各符号含义、适用条件与反例。比如“马赫数”卡片：定义是V/a；公式Ma=V/a，其中a=sqrt(γRT)；适用条件是无所谓，但用途是判断压缩性；反例是有人把Ma<1全叫亚音速，其实跨音速已经在翼面局部超音速。
公式卡片的重点不是公式本身，而是“此公式在什么假设下才成立”。例如伯努利方程p + 0.5ρV² + ρgh = const只沿同一条流线成立，且要求定常、无粘、不可压。如果你在边界层里用伯努利，就是自找麻烦。把这个“适用条件”写得越严格，后面复习越省力。
我会把PDF里的所有相关曲线重新画一遍草图：升力系数随攻角曲线、极曲线、压力系数分布、激波前后参数表。画图不用精美，关键在于自己动手标注每个拐点的物理含义。失速点对应什么？极曲线上最小阻力点对应什么攻角？这些一旦画过，就不会忘。
4. 避坑：自学空气动力学最容易翻车的五个概念
下面这条是通读讲义、同时做了多轮工程估算后总结的。以下每一个坑我都亲身踩过，有的在考试里扣分，有的在仿真结果里绕了几天才发现。按“现象、原因、解决”三步写出来，希望能帮你少走这些弯路。
4.1 压力中心与焦点的混淆：稳定性算错一整片
现象：做飞机纵向静稳定性估算时，有人用压力中心位置随攻角变化来判定稳定，算出焦点在重心之后，却把结论写反；或者拿两个不同攻角下的压力中心说事，得出“压力中心随攻角变化所以不稳定”的错误结论。
原因：压力中心是气动合力作用点，它受攻角影响明显。纵向静稳定性看的是“攻角扰动后，气动力矩是否产生恢复效应”，这需要用焦点（气动中心）来定义，而不是压力中心。焦点被定义为力矩系数斜率dCm/dα为常数的参考点，实际上就是矩心变化不影响斜率的那一点。判定静稳定，看焦点与重心的相对位置，焦点在前重心在后则不稳定。
解决：把所有力矩分析都换成焦点。从讲义里找“焦点定义和推导力矩系数”的那一节，用一次代数运算验证：如果把力矩系数对攻角求导，就会发现焦点处的导数与参考点无关。计算时用一个固定参考点（比如平均气动弦前缘），求出dCm/dα和dCL/dα，再用公式x_focus = x_ref - (dCm/dα)/(dCL/dα)（无因次化）求出焦点，然后对比重心位置。这一步不能省。
4.2 把伯努利方程当作一切升力的解释
现象：解释升力时只讲“上表面气流快压力低，下面慢压力高”，然后试图用伯努利反推上表面流速，误差很大。还有人用这个去解释平板大攻角分离流，算出的升力方向和大小都对不上。
原因：伯努利是能量方程沿流线的体现，它把速度场和压力场联系起来，但不解释为什么上表面流速会被加快。升力的根源是流体偏转带来的动量通量变率，等价于环绕翼型存在净环量，可以用库塔-茹可夫斯基定理表示。只有在附着流、无粘、定常的条件下，伯努利才是一个便捷的计算工具，把它当物理因果就反了。
解决：遇到升力问题，先画控制体积，用动量定理写出升力等于控制面上下压力积分或动量通量积分。然后用自由涡系（起动涡+附着涡）说明环量来源。最后再用伯努利去算压力分布，得到和库塔-茹可夫斯基一致的结果。这一套流程走完，你才能对“升力从哪来”有完整认识。
4.3 忽略雷诺数相似：风洞数据和飞行数据对不上
现象：从文献里查一个翼型的最大升力系数，直接用于自己的飞行器设计，结果试飞时发现失速速度比设计值高，升力不足。
原因：文献数据对应某个雷诺数，可能是Re=1e6或3e6，而你的实际工况只有Re=2e5。低雷诺数下边界层更容易分离，最大升力系数明显下降，失速攻角提前。雷诺数相似是风洞实验的第一前提，少了这个，数据只能是装饰。
解决：看数据时先确认Re范围。如果是小尺寸低速无人机，就不要采用大型飞机的Re数据。如果没有对应Re的数据，可以降低攻角安全裕度，或者用转捩带提前固定转捩点来模拟高Re。CFD也是一样，湍流模型和转捩模型对Re敏感，边界条件里输入速度、弦长、粘度后，先算一下Re是否和参考工况一致。这个自查两分钟，能省掉一整晚的重新计算。
4.4 把摩擦阻力与压差阻力混在一起调参
现象：模型阻力偏大，有人认为是表面太粗糙，把CFD壁面粗糙度调到0，结果阻力几乎没变，甚至分离区更大。
原因：气动阻力分解为摩擦阻力和压差阻力。粗糙度主要影响湍流摩擦阻力，但压差阻力由逆压梯度引起的分离区大小决定。如果机翼已经在大攻角下有分离，压差阻力占总阻力大头，抛光表面当然没效果。
解决：用CFD后处理把阻力拆开看。从几何表面做积分，摩擦阻力用壁面剪切应力面积分，压差阻力用表压沿法向投影积分。如果压差阻力占比超过70%，你的调整方向是翼型弯度、前缘半径、加装涡流发生器或调整攻角，而不是改粗糙度。PDF里型阻极性和极曲线的部分，会告诉你最小阻力点在哪里，尽量在那个攻角附近飞行，是最直接的调整。
4.5 总压、静压、动压：测出来的是什么
现象：用皮托管测机翼表面压力，把读数直接当静压算压力系数，结果压力系数偏差极大。
原因：皮托管对着来流中心孔感受到的是滞止压力，即总压，侧面孔才测静压。如果探针安装角度不正，总压管还可能感受到不完全滞止，读数异常。把总压当静压，相当于动压被重复计入，速度越高压差越大。
解决：先用一个三孔或五孔探针的校准表，确认偏航角和俯仰角对测压的影响。手算时记住：总压p0 = 静压p + 动压q。在低速不可压缩流中，q = 0.5ρV²；在可压缩流中，需用等熵关系修正，q = p0·(1 - (p/p0)^((γ-1)/γ))·γ/(γ-1)（实际上是总静压比公式的变形）。无论做实验还是CFD，都要明确每个输出的压力是总压还是静压。CFD里总压常标为p0，静压标为p，别混着用。
5. 一个值得养成的验证习惯：动压与马赫数的快速手算
这套课程给我留下最深的影响，不是背下了哪个公式，而是养成了“开工前先手算一个量级”的习惯。现在我拿到一个气动问题，无论跑不跑CFD，都会先算动压、马赫数和雷诺数，再把预期的升力、阻力数量级写在纸上。
比如一个飞行器在海平面飞行，速度137m/s，其中弦长1.2m，标准大气密度1.225kg/m³，粘度1.79e-5Pa·s。那么动压 q = 0.5×1.225×137² ≈ 11500Pa，马赫数Ma = 137/340 ≈ 0.4，已经超过0.3，不可压缩近似需要谨慎使用，至少要做普朗特-格劳厄特修正。雷诺数Re ≈ 1.225×137×1.2/1.79e-5 ≈ 1.12e7。如果预期Cl=0.8，单位面积升力就是 0.8×11500 ≈ 9200Pa，对应升力比重量超出多少，马上能判断设计是否合理。
在可压缩修正上，有一个足够用的技巧：当Ma在0.3到0.7之间，用普朗特-格劳厄特修正Cp = Cp_incompressible / sqrt(1 - Ma∞²)评估翼型表面压力系数变化。这样即使不做完整CFD，也能判断压缩性会不会让吸力峰提前达到临界马赫数。这个细节我曾在一次方案讨论中用来质疑一个高翼载无人机的巡航速度设计，避免了一次绿皮初样阶段的返工。
超音速工况下，我会先查一维等熵流表或直接用公式算总压比、静压比。Ma=2.0时静压比p/p0约0.128，驻点温度比T0/T约1.8。把这个记在心里，看到风洞报告里的压力探针读数就能快速排除一个量级误差。每一次这样的手算都在验证你是否真理解了PDF里的物理模型，而不是把CFD结果当黑匣子吞下去。
我自己的教训是，有次为了赶节点，跳过手算直接起了一个跨音速算例，结果求解器把速度入口压力给成表压而非绝对压力，导致全流场密度差了一倍。这种错，靠残差曲线看不出来，靠手算一分钟就能发现。现在我只信两条：先定参数量级，再动仿真软件。希望帮到你。
本文还有配套的精品资源，点击获取
