简介这份PDF是刘沛清主讲的北京航空航天大学《空气动力学基础》精品课程配套讲义面向航空航天相关专业学生、考研复习者及需要系统梳理空气动力学知识框架的工程技术人员帮助从零建立流体力学与空气动力学的基本概念和分析方法。内容按课程章节展开涵盖绪论、流体的基本属性、流体静力学、流体运动学与动力学基础、不可压缩无粘流体平面势流、粘性流体动力学基础、边界层理论及其近似以及可压缩高速流动基础等核心模块并配有相对飞行原理、风洞实验原理、机翼绕流与边界层分离等应用讲解图文结合便于理解与复习。资源包共1个文件为PDF格式大小19.65MB目录结构与课程章节对应清晰既适合自学入门也可作为课堂教学和考前梳理的重要参考。目前已有674人学习下载适合希望借助北航精品课程体系掌握空气动力学核心知识的读者。1. 空气动力学基础这份北航精品课PDF为什么说钥匙在第零章空气动力学基础这份北航精品课程PDF我拿到手后先翻的不是公式而是第零章绪论里那幅风洞示意图——整套课的钥匙就藏在最前面相对飞行原理。飞行器以速度V在静止空气中飞行与飞行器固定不动、空气反向等速流过两者流场和受力完全等效。后续章节都在回答这个原理引出的问题流体属性与静力学怎么建立压强概念运动学怎么描述流场平面势流怎么给出升力边界层理论怎么解释阻力可压缩流动怎么修正密度变化。适合航空航天专业学生补体系、CFD工程师找理论边界、航空爱好者搞懂机翼升力来源。它解决的核心问题是在设边界条件、调湍流模型之前先把流动机理和适用条件弄对。2. 从三态到流体的静力学把“静止流体不能承受剪力”变成判断依据2.1 固态、液态、气态为什么这句话是整套课的起点课件绪论把物质三态和力学定义放在最前面初看像凑页数实际是整套课的地基。固体在静止状态能承受拉力、压力和剪切力液体和气体在静止状态只能承受压力几乎不能承受拉力和剪切力。这句话直接决定后面所有建模的边界只要流体处于静止或匀速直线运动切应力为零表面力只剩法向压力也就是压强。你在CFD软件里设入口总压还是静压、能不能直接用伯努利方程判断依据就在这句话里。从微观视角更容易记液体和气体的分子间距大、约束弱分子热运动让它们“见缝就钻”宏观表现就是易流动性。所以流体的定义并不是“会流动的东西”而是“静止时不能承受剪切力、一旦有剪切力就持续变形”的介质。这个定义比教材里那句“液体和气体的统称”实用得多——当你判断某种材料能不能按流体处理看的是它对剪切力的响应而不是它分子排列到底像不像液体。2.2 流体静力学压强随深度变化、U形管和连通器静力学部分的核心其实只有一条主线重力场下静止流体内部的压强随深度线性增加。表达式是 (p p_0 \rho g h)其中 (p_0) 是自由面压强(\rho) 是流体密度(h) 是深度。课件里讲的是概念实际做题时容易忽略两个前提第一这个公式要求流体不可压缩且密度均匀如果上部是气体柱、下部是液体柱要分段计算第二压强虽然是标量但计算方向性受力时要区分“压力”和“压强”前者是力、后者是单位面积上的力。U形管压差计是实验室最常见的工具测两点压差的公式是 (\Delta p (\rho_m - \rho) g \Delta h)(\rho_m) 是工作液密度(\rho) 是被测流体密度。当被测流体是气体时(\rho) 通常比 (\rho_m) 小两三个数量级公式简化成 (\Delta p \approx \rho_m g \Delta h)但当被测流体是液体时(\rho) 不能扔否则误差直接到百分之几十。课件里“流体静力学的应用”一节没有展开推导我一般会要求自己把U形管的两种接法——接在壁面静压孔上测静压、接在皮托管上测总压与静压之差——都推一遍后者直接对应飞行速度的测量原理。2.3 量纲与单位最容易翻车的换算以及我的自查方法量纲部分看似废话实际是整套课里我栽过跟头最多的地方。速度用 m/s、压强用 Pa、密度用 kg/m³这三个基本单位一组合动力粘度 (\mu) 的单位是 Pa·s 或 kg/(m·s)运动粘度 (\nu \mu / \rho) 的单位是 m²/s。CFD软件里最常犯的错误是把运动粘度和动力粘度混用空气在标准大气条件下动力粘度约 (1.8 \times 10^{-5}) Pa·s密度约 1.225 kg/m³运动粘度约 (1.5 \times 10^{-5}) m²/s两者数值接近但物理意义完全不同算雷诺数时用错一个结果差一个量级。物理量SI单位常用单位换算关系动力粘度 (\mu)Pa·s kg/(m·s)泊(P)1 P 0.1 Pa·s1 cP 1 mPa·s运动粘度 (\nu)m²/s厘斯(cSt)1 cSt (10^{-6}) m²/s压强 (p)Pabar / atm / mmHg1 bar (10^5) Pa1 atm 101325 Pa密度 (\rho)kg/m³g/cm³1 g/cm³ 1000 kg/m³提示拿到任何公式先做量纲一致性检查。伯努利方程里 (p \frac{1}{2}\rho V^2 \rho g z) 每一项量纲都是 Pa只要有一项对不上要么是推导错误要么是录入错误。我自己的习惯是每学完一章把最终公式的每一项量纲写出来和主干公式对比一次。这个过程能拦住至少一半的低级错误——尤其是从工程图纸或实验结果里抄数据时单位不统一是最常见的坑。3. 运动学与动力学基础流线、连续性方程和伯努利方程怎么串成一套3.1 流线与迹线别混淆“某一时刻”和“某一质点”流体运动学首先解决的是“怎么描述流动”。课件里的关键区分是拉格朗日描述和欧拉描述前者盯着一个流体质点记录它的位置随时间变化后者盯住空间固定点记录不同质点经过该点时的速度、压强和密度。工程上绝大多数问题用欧拉描述因为我们要知道的是流场在空间中的分布而不是某个质点的完整轨迹。由此引出三个容易混淆的概念。流线是某一时刻流场中的一条曲线曲线上每点的切线方向与该点速度方向一致迹线是一个质点在一段时间内走过的轨迹脉线是连续释放的示踪粒子在某一时刻连成的线。定常流场中三者重合非定常流场中都不重合。我见过很多人在分析翼型绕流时把非定常的涡脱落迹线当成瞬时流线来讨论概念一混后面的升力计算全乱。与流线配套的是物质导数概念流场中某物理量对时间的变化率等于当地导数加迁移导数。用速度分量表示就是 ( \frac{D}{Dt} \frac{\partial}{\partial t} u\frac{\partial}{\partial x} v\frac{\partial}{\partial y} w\frac{\partial}{\partial z} )。这个公式看着简单却是从运动学进入动力学的桥梁——加速度不是简单的 (\partial V / \partial t)还要加上流体质点从一处运动到另一处所感受到的变化。忽略迁移加速度是很多初学者把欧拉方程理解成普通牛顿第二定律的根源。3.2 连续性方程与伯努利方程适用条件比公式本身重要连续性方程的本质是质量守恒。定常、不可压缩条件下的形式是 (A_1 V_1 A_2 V_2)即面积与速度成反比。课件里讲的是微分形式 (\nabla \cdot \mathbf{V} 0)意思是速度场的散度为零既没有质量源也没有质量汇。实际使用时要注意这个简化同时要求“定常”和“不可压缩”如果密度变化明显必须用 (\rho_1 A_1 V_1 \rho_2 A_2 V_2)可压缩流动章节会专门处理这个问题。伯努利方程是从欧拉方程沿流线积分得到的前提是定常、无粘、不可压缩、沿同一条流线。这四个条件缺一不可。很多人把它当成“速度大压强小”的万能结论其实忽略了它只适用于无粘流场在边界层内部、激波后部、强旋流区域用伯努利方程结果会和实验差很远。课件把它放在“流体动力学基础”而不是“静力学”里本身就是提醒它描述的是流动过程中机械能的转换关系不是静态的因果关系。3.3 用Python做一段伯努利方程数值验证我学这一章时习惯用一个最简单的变截面管道来验证方程的记忆是否正确。假设空气密度不变忽略高度差和壁面摩擦给定入口状态和出口面积用连续性方程求出口速度再用伯努利方程求出口静压。代码如下import math # 管道参数SI单位 rho 1.225 # 空气密度 kg/m^3 A1 0.10 # 入口面积 m^2 A2 0.05 # 出口面积 m^2 V1 10.0 # 入口速度 m/s p1 101325.0 # 入口静压 Pa近似标准大气压 # 连续性方程rho*A1*V1 rho*A2*V2 # 密度不变时A1*V1 A2*V2 V2 A1 * V1 / A2 # 伯努利方程p1 0.5*rho*V1^2 p2 0.5*rho*V2^2 # 整理得 p2注意这里忽略高度差和粘性损失 p2 p1 0.5 * rho * (V1**2 - V2**2) print(f出口速度 V2 {V2:.2f} m/s) print(f出口静压 p2 {p2:.2f} Pa) print(f压差 p1 - p2 {p1 - p2:.2f} Pa) print(f质量流量 {rho * A1 * V1:.3f} kg/s)运行结果是出口速度 20.00 m/s出口静压 101141.25 Pa入口与出口压差约 183.75 Pa。逻辑说明因为面积缩小一半速度增大一倍动压增大导致静压下降这正是文丘里管测速的原理。参数说明(rho) 用的是标准大气条件下的空气密度如果做水洞实验要改成 (1000) kg/m³(V1) 超过 30 m/s 时需要先检查马赫数是否低于 0.3否则“不可压缩”假设不再成立。这个脚本的价值不在于计算本身而在于把公式变成可验证的数字做题时算出的压强如果是负数通常不是方程错了而是参数组合已经超出模型适用范围。4. 平面势流与粘性动力学位函数、N-S方程和升力来源4.1 不可压缩无粘流假设势流理论能用在哪、不能用在哪第四章是整个课件里最“数学”的一章也是理解升力来源的关键。不可压缩无粘平面势流有两个基本假设流体无粘所以无剪切力流动无旋所以存在速度势函数 (\phi)满足 (\nabla^2 \phi 0)也就是拉普拉斯方程。速度分量可以由势函数求导得到(u \partial \phi / \partial x)(v \partial \phi / \partial y)。无旋条件还引出一个好处伯努利方程在无旋流场中不再局限于同一流线可以在全流场成立这给压强计算带来极大便利。但势流理论有个致命短板无粘流体绕物体流动时物体表面满足法向速度为零却允许切向滑移算出来的阻力恒为零——这就是著名的达朗伯佯谬。课件交代历史背景时一定会提到这一点理想流体力学在18世纪到19世纪发展得很完美却解释不了为什么飞机会受到阻力。所以势流理论的价值不在于给阻力而在于给升力只要边界层很薄、流动没有大范围分离外部流动可以用势流近似机翼升力的主项就是从这个无粘模型里算出来的。4.2 基本解叠加用Python模拟圆柱绕流表面速度势流理论最有用的工具是叠加法。均匀流、点源、点汇、点涡、偶极子都是拉普拉斯方程的基本解把它们叠加起来就得到绕不同外形的流动。经典的例子是均匀流加偶极子模拟无粘圆柱绕流在此基础上再叠加一个点涡圆柱上下表面速度就不对称产生升力。库塔-茹科夫斯基定理 (L \rho V_\infty \Gamma) 把升力与环量直接挂钩这条线就是从势流理论长出来的。import math # 来流参数 V_inf 10.0 # 远前方均匀流速度 m/s R 1.0 # 圆柱半径 m N 16 # 取16个点覆盖0~360度 print(角度(deg) 表面速度(m/s) 表面压强系数Cp) for n in range(N): theta_deg n * 360 / N theta math.radians(theta_deg) # 均匀流 偶极子无粘不可压势流圆柱绕流表面切向速度 # 势函数 phi V_inf * (r R^2 / r) * cos(theta) # 在 r R 处V_theta -2 * V_inf * sin(theta) V_surface abs(-2.0 * V_inf * math.sin(theta)) # 压强系数 Cp 1 - (V_surface / V_inf)^2 Cp 1.0 - (V_surface / V_inf) ** 2 print(f{theta_deg:10.1f} {V_surface:14.2f} {Cp:14.3f})代码逻辑说明在圆柱表面 (rR) 处速度分布只保留切向分量(V_\theta -2V_\infty \sin\theta)。(\theta0^\circ) 和 (180^\circ) 是前后驻点速度为零压强系数 (Cp1)(\theta90^\circ) 和 (270^\circ) 是上下两侧速度达到最大值 (2V_\infty)(Cp-3)。运行这段代码会看到上下表面速度完全对称压强积分后升力为零——这就是达朗伯佯谬的数值体现。叠加一个点涡后上表面速度增大、下表面速度减小压强差产生升力。参数说明(N) 控制采样点数改成 360 可以得到更光滑的压强分布曲线(R) 和 (V_\infty) 在无粘无旋流动中不改变 (Cp) 分布只改变速度的绝对值。4.3 粘性动力学与雷诺数为什么无粘算不出阻力粘性流体动力学基础这一章把问题拉回现实。实际流体都有粘性完整描述流动的是纳维-斯托克斯方程。这个方程比欧拉方程多了粘性项对二维不可压缩流动形式是[ \frac{\partial u}{\partial t} u\frac{\partial u}{\partial x} v\frac{\partial u}{\partial y} -\frac{1}{\rho}\frac{\partial p}{\partial x} \nu \left( \frac{\partial^2 u}{\partial x^2} \frac{\partial^2 u}{\partial y^2} \right) ]方程左边是惯性项右边第一项是压力梯度第二项是粘性扩散项。判断谁主谁次要看无量纲数——雷诺数 (Re \rho V L / \mu)。(Re) 很小粘性项主导流动平滑比如微流控芯片里的液体(Re) 很大惯性项主导粘性只在靠近壁面的薄层里起作用这就是边界层的雏形。课件从“粘性流体动力学基础”过渡到“边界层理论”逻辑上正是用 (Re) 做分水岭。我一般会提醒自己记住几个数量级小型无人机翼弦 0.3 m、速度 20 m/s(Re) 约 (4 \times 10^5)边界层会经历层流到湍流的转捩客机翼弦 3 m、巡航速度 250 m/s(Re) 到了 (10^7) 量级几乎整段边界层都是湍流。不同 (Re) 区间边界层厚度、摩擦阻力系数、分离位置都不一样数值模拟时选湍流模型的标准也由它决定。这就是为什么雷诺数不是“算出来的一个数”而是流动状态的身份证。5. 边界层理论常见问题分离、转捩与网格加密的三个坑5.1 三个典型翻车现场第一个坑是“势流算出来阻力为零实验却测到很大阻力”。现象用上一章的圆柱绕流势流解积分后升力阻力全为零但风洞实验里圆柱明明受到明显的阻力。原因势流假设无粘没有考虑边界层分离实际流动中边界层在圆柱背风面分离形成低压尾迹区前后压力不平衡。解决先用雷诺数判断流动状态当 (Re) 大于临界值时直接放弃无粘模型改用边界层理论估算分离点或者直接上CFD。第二个坑是“把边界层外缘速度当成来流速度直接用”。现象用经验公式算平板摩擦阻力结果和实验数据对不上偏小约 5% 到 10%。原因边界层存在位移厚度 (\delta^*)它把外部无粘流场的有效壁面向外推了一段距离导致边界层外缘的实际速度比远前方来流速度大。解决计算升力、表面压强时先用位移厚度修正有效外形再用修正后的速度作为边界层外缘条件。第三个坑是“CFD里第一层网格高度随便给收敛后结果完全不能用”。现象壁面热流或摩擦阻力严重偏大或者残差迟迟降不下去。原因边界层内部速度梯度大第一层网格高度决定壁面 (y^) 值。如果 (y^) 在几十甚至几百而湍流模型要求解析粘性底层那么近壁速度分布根本解析不出来。解决先估算边界层厚度 (\delta)根据所使用的湍流模型确定目标 (y^)层流模型要求 (y^ \approx 1)壁面函数法可以放宽到 30 到 300再用公式反算第一层网格高度。注意(y^) 的估算公式很多工程软件都内置但前提是来流速度、密度、粘度和壁面到前缘的距离都要取对。很多人输错位置长度 (L)算出来的 (y^) 和实际相差一个量级。5.2 边界层厚度判断层流、转捩与湍流的数量级直觉边界层理论章节里的关键结论是平板零压梯度层流边界层厚度约 (\delta \approx 5x / \sqrt{Re_x})。这个公式要求记住但更要紧的是对数量级有直觉。空气运动粘度约 (1.5 \times 10^{-5}) m²/s来流速度 30 m/s 时不同弦长位置的雷诺数和层流厚度如下位置 (x)雷诺数 (Re_x)层流边界层厚度 (\delta)工程判断0.1 m(2.0 \times 10^5)约 1.1 mm未达转捩临界 (5 \times 10^5)通常仍为层流1 m(2.0 \times 10^6)若保持层流约 3.5 mm实际通常已转捩为湍流湍流边界层更厚湍流边界层厚度比层流厚得多近似公式 (\delta \approx 0.37x / Re_x^{1/5})同样条件下 1 m 处约 5 到 6 mm。这个数量级直觉有什么用当你做翼型网格时第一层网格高度要在 0.01 mm 量级边界层内至少要布置 15 到 20 层网格才能勉强解析近壁速度剖面。如果边界层只有 1 mm而你画的边界层网格总厚度先给了 10 mm那一半网格都在边界层外面完全白费。这是我常说“边界层厚度是玄学但数量级是科学”的原因先用公式估一遍再用软件里的分布函数实测回流两边一对心里就有底。5.3 实验、理论与CFD怎么分工课件绪论里提到空气动力学的研究方法有理论、实验和数值模拟三类很多人以为它们是并列关系实际上三者各有不可替代的位置。实验的优点是真实缺点是风洞有支架干扰、洞壁干扰模型缩比带来雷诺数不匹配理论解干净漂亮但只能解少数规则外形实际机翼都要靠数值方法CFD能给出全场细节但湍流模型本身就是对真实流动的近似转捩位置、分离点预测不准是常事。我的习惯是拿理论解做基准测试。比如用平板边界层公式算一组阻力系数再用同等条件下的CFD结果对比如果偏差超过 10%先怀疑网格和湍流模型而不是怀疑理论公式。反过来实验数据与CFD对不上时优先检查两者是否保证了相同的雷诺数和来流湍流度——风洞来流湍流度会直接改变转捩位置而CFD里默认的自由来流湍流度通常是 1%这个差异足以让边界层分离位置移动几个百分点弦长。三者对照着看才能避免单一方带来的系统性偏差。6. 把课件读薄四个自测题与一个“反查”习惯四份PPT课件最容易被忽略的地方是绪论里的“研究方法”和“发展进程”两页。我当时直接把这两页跳过了后来做项目时才意识到理解“理论、实验、数值模拟”的边界比多背两个公式更有用。所以我把这门课的最后一步设计成四个自测题每看完一章就拿来问自己序号自测题判断标准1给一个翼型截面能否判断前缘驻点、后缘驻点以及上下表面哪边流速大能用伯努利方程反推压强分布并说明下表面压强大于上表面2给定飞行高度和速度能否在一分钟内估算雷诺数并判断边界层状态能写出 (Re \rho V L / \mu)并给出层流或湍流的判断3能否用自己的话解释“相对飞行原理”为什么让风洞实验成立能说清参考系变换并指出地面坐标系与风洞坐标系下流场的等效性4能否指出N-S方程中哪一项是迁移加速度、哪一项是粘性项忽略各有什么代价能正确标出 (u\partial u/\partial x) 和 (\nu\nabla^2 u)并说明简化后果我还养成一个“反查”习惯每学完一个中文术语翻到英文原版教材里找对应表达。比如“边界层位移厚度”对应 displacement thickness它的物理图像是“流线被向外排挤”而不是“壁面变厚”“库塔条件”对应 Kutta condition它解决的是后缘流动光滑脱离的问题中文教材里往往一句话带过英文教材里会画出三种后缘流动的对比图。这个习惯帮我避免了大量“名词都认识、联想不到一起”的尴尬特别是读论文时中文术语和英文术语对不上几乎等于白学。这门课的容量不大但体系非常完整从第零章的相对飞行原理一路推到可压缩高速流动基础正好覆盖亚声速到跨声速分析的主线。我当年学的时候只盯着公式背结果被问“机翼为什么能产生升力”时第一反应答成伯努利效应被人反问“倒过来飞为什么还能产生升力”当场卡住。从那以后我每学一章都强制走一遍自测题和英文反查这个习惯帮我撑过了后来大多数流体力学的实际项目。希望帮到你。本文还有配套的精品资源点击获取