简介《空气动力学基础》是北京航空航天大学北京市精品课程的配套PDF课件面向航空航天、力学及相关专业的本科生、考研复习者与自学读者帮助建立从流体属性到高速可压缩流动的系统知识框架。内容按绪论、流体基本属性与静力学、流体运动学与动力学、不可压缩无粘平面势流、粘性流体动力学、边界层理论及可压缩高速流动等模块展开并穿插飞行器绕流、风洞实验、边界层分离等典型图例便于对照课堂进度梳理概念、公式与工程应用。整包仅含1个PDF文件大小约19.65MB页面以课程讲义形式组织目录层次清晰适合按章检索和打印批注。目前已有686人学习下载可作为课程预习、期末复盘和考研专业课复习的参考材料。1. 六章骨架先摆出来空气动力学不是“算升力”而是逐层拆假设很多人翻开一份空气动力学讲义第一反应是直奔机翼升力那一页结果卡在物质三态、连续介质、量纲单位这几页后面的势流、边界层、激波全成了断点。这门课真正的骨架是一条“拆假设”的路线先假设空气连续、无粘、不可压得到一套干净的势流解再把“无粘”拆掉换来粘性方程和边界层最后把“不可压”拆掉换来马赫数、总温总压和激波。讲义里的六章就是这三刀切出来的——绪论给坐标系第二三章建理想模型第四章拆粘性第五章处理边界层第六章进高速。它适合三类人正在上这门课要交作业和考试的、准备复试被追问基础概念的、做无人机或气动仿真需要补理论底子的。2. 流体属性与静力学从三态定义到欧拉平衡方程的复习重点这一章看起来最“水”实际上后面所有公式的合法性都建立在这里。三态定义决定了“流体为什么只能承受压力”连续介质假设决定了 N-S 方程能不能写出来量纲决定了你推导完的式子敢不敢用。2.1 三态对比与流体的严格定义讲义把固、液、气三态的区别落在两个维度上有没有固定形状、有没有固定体积以及静止状态下能承受哪种应力。这三条一列出来后面“流体静止时只有法向压强”就是自然推论而不是凭空规定。物态固定形状固定体积静止时承受拉力静止时承受压力静止时承受剪力固体有有能能能液体无有几乎不能能不能气体无无几乎不能能不能流体是液体和气体的统称核心特征只有一句易流动性静止状态下不能承受剪应力。这句话直接推出静力学的全部前提——既然静止时没有剪应力那么作用在微元体上的表面力只有法向压强压强就成了静力学里唯一的表面力变量。提示判断某类物质能不能当流体处理看的不是“它像不像水”而是“静止时能不能承受剪应力”。沥青、蜂蜜、熔融金属都算流体只是粘度大。讲义还特别强调了物质形态由分子热运动和分子间作用力共同决定。气体分子间距大、热运动占主导所以既没有固定体积也没有固定形状液体分子间作用力还够强能锁住体积但锁不住形状。这段定性描述在考试里常以选择或简答出现背结论不如记住“分子间距 热运动”这个因果链。2.2 连续介质假设与流体的四个基本属性连续介质假设是整门课的第一块地基把由大量分子组成的空气看成连续充满空间、物理量处处有定义的介质。它成立的条件是分子平均自由程远小于流动的特征尺度工程上用克努森数衡量Kn λ/L一般认为 Kn 0.01 时连续介质假设成立。讲义没有展开这一条但它是后面所有微分方程的前提复习时值得单独记一笔。在这块地基上流体有四个反复出现的属性每一个都对应后面一章的主线属性定义式或度量常温空气典型值后续用在哪密度 ρ单位体积质量1.225 kg/m³海平面标准大气连续方程、伯努利、动量方程可压缩性体积模量 K -V·dp/dV或声速 a √(dp/dρ)a ≈ 340.3 m/s第六章可压缩流的分区依据粘性 μ牛顿内摩擦定律 τ μ·du/dy1.789×10⁻⁵ Pa·s第四章 N-S 方程、边界层表面张力 σ气液界面单位长度上的力空气与水的界面约 0.072 N/m主要是液体问题气体通常忽略可压缩性这一栏要特别注意气体“容易被压缩”不等于“任何流动都必须按可压缩算”。判据是马赫数Ma 0.3 时密度相对变化小于 5%按不可压处理误差可接受。这条界线是第三章和第六章的分水岭也是很多计算题的隐含条件。2.3 静压强的两个特性与欧拉平衡微分方程静压强有两个必须背准的特性第一方向沿作用面的内法线也就是垂直指向受压面第二同一点的压强大小与作用面方位无关各个方向相等。第二条成立的前提正是“静止流体没有剪应力”所以它不是数学结论而是 2.1 节定义的直接后果。把这两个特性和质量力一起作用在微元六面体上就得到欧拉平衡微分方程dp ρ(X dx Y dy Z dz)其中 X、Y、Z 是单位质量质量力在三个方向的分量。在重力场里只有 Z -g方程退化成 dp -ρg dz积分得到 p p₀ ρgh这就是连通器、测压管、U 形管差压计全部题目的来源。等压面是这一节的另一个高频考点等压面上 dp 0所以质量力与等压面正交。重力场中等压面是水平面这一条能秒杀一批“斜置容器里液面怎么画”的题。压强还有三套口径必须分清绝对压强以完全真空为零点物理方程里用的都是它表压绝对压强减去当地大气压压力表读数真空度当地大气压减去绝对压强只在低于大气压时存在。注意伯努利方程、状态方程里代入的压强必须是同一套口径。表压和绝对压强混着代量级上差 10⁵ Pa结果会离谱到看不出来错在哪。2.4 量纲与单位MLT 和 FLT 两套体系别混用量纲分析在讲义里只占一节但它是检查推导的最省力工具。同一物理量可以挂在两套基本量纲体系下MLT质量、长度、时间和 FLT力、长度、时间两者通过 F MLT⁻² 互换。物理量MLT 量纲FLT 量纲SI 单位密度 ρML⁻³FL⁻⁴T²kg/m³压强 pML⁻¹T⁻²FL⁻²Pa动力粘度 μML⁻¹T⁻¹FL⁻²TPa·s运动粘度 νL²T⁻¹L²T⁻¹m²/s力 FMLT⁻²FN复习时养成一个习惯任何一条推出来的公式先看两边量纲是否齐次。量纲不齐一定是抄漏了系数或者漏了一项这比重新推一遍快得多。量纲分析的另一个用途是导出无量纲数雷诺数、马赫数、欧拉数都是这么来的第四章讲相似准则时会再遇到。3. 运动学与无粘势流三组方程撑起的“理想空气”主线第三章是整门课最“漂亮”的一章也是考试计算题最集中的一章。它建立在两个假设上——无粘、不可压——在假设之内流动可以用一个标量函数描述叠加原理随便用一旦出了这个假设就要回到第四章。3.1 拉格朗日与欧拉两种描述法怎么选描述流体运动有两套语言。拉格朗日法跟踪每一个流体质点记号是 (x₀, y₀, z₀, t)关心“这个质点去哪儿了”欧拉法盯住空间点记号是 (x, y, z, t)关心“这个位置上此刻的速度是多少”。空气动力学几乎全部用欧拉法因为实验测量风洞里的探针、仿真里的网格点天然是欧拉视角。配套的三个几何概念必须分清流线同一时刻线上每一点的切线方向就是该点速度方向迹线同一个流体质点在一段时间内走过的轨迹脉线连续流过同一空间点的流体质点连成的线。定常流动时三者重合非定常时一般不重合。考题里“定常”两个字就是提示看到它就可以放心把流线当迹线用。欧拉法带来的核心运算叫质点导数也叫随体导数DV/Dt ∂V/∂t (V·∇)V右边第一项是当地加速度固定空间点上速度随时间的变化第二项是迁移加速度质点从一个位置挪到另一个位置因速度场空间不均匀而产生的变化。很多同学在算加速度时只写 ∂V/∂t把迁移项漏掉这是运动学部分最常见的失分点。3.2 连续方程、动量方程与伯努利方程三组方程里连续方程管“质量守恒”动量方程管“力与加速度”伯努利方程是动量方程在特定条件下的积分。连续方程的微分形式是 ∂ρ/∂t ∇·(ρV) 0不可压条件下密度是常数退化成 ∇·V 0也就是速度场无源无汇。定常管流的形式更常用ρ₁V₁A₁ ρ₂V₂A₂不可压时简化为 V₁A₁ V₂A₂。动量方程在无粘条件下的形式就是欧拉方程加上粘性项之后就是第四章的 N-S 方程。考试里更常出现的是它的积分形式控制体上的合外力等于流出与流入的动量通量之差。伯努利方程是沿流线积分欧拉方程的产物形式是p ½ρV² ρgz 常数它的四个前提必须一字不差地记住定常、无粘、不可压、沿同一条流线。第四条最容易被忽略——两条不同流线上的常数一般不同只有在无旋流场里常数才全场统一。这一点在第五章还会作为易错点单独展开。伯努利最直接的工程应用是皮托管测速。把总压孔和静压孔测到的压差代进去速度就是 V √(2(p₀ - p)/ρ)。同样的思路反过来用就是风洞实验里通过皮托管排测量速度场的方法。3.3 平面势流速度势、流函数与四种基本解无粘不可压再加一条“无旋”流场就变得极其规整。无旋意味着速度的旋度为零于是存在一个标量函数 φ使 V ∇φ这个 φ 叫速度势不可压意味着 ∇·V 0代进去就得到拉普拉斯方程 ∇²φ 0。不可压平面流动还另外存在一个流函数 ψ等 ψ 线就是流线。拉普拉斯方程是线性的所以解可以叠加这就是势流理论能玩出花来的原因。四种基本解搭起了所有复杂流动基本解速度势 φ流函数 ψ速度场特征均匀流UxUy处处速度相同源汇(m/2π)·ln r(m/2π)·θ径向辐射强度 m点涡(Γ/2π)·θ-(Γ/2π)·ln r周向旋转环量 Γ偶极子(M/2π)·cosθ/r-(M/2π)·sinθ/r源汇无限靠近的极限实际流动靠叠加得到均匀流加偶极子就是圆柱绕流均匀流加偶极子再加点涡就是有环量的圆柱绕流后者能算出升力升力大小是 ρUΓ这就是库塔-儒可夫斯基定理。这条定理是整章的高潮也是后面解释机翼升力为什么和环量挂钩的起点。下面这段代码把圆柱绕流的表面压力系数算出来可以直观看到势流解的几个标志性数字import numpy as np # 无粘不可压圆柱绕流半径 R1计算表面速度与压力系数 R 1.0 U 1.0 # 远前方来流速度m/s rho 1.225 # 海平面标准大气密度kg/m^3 theta np.linspace(0, 2*np.pi, 9) # 从 0 到 360 度取 9 个点 Vt -2 * U * np.sin(theta) # 势流解圆柱表面切向速度 Vθ -2U·sinθ Cp 1 - (Vt / U) ** 2 # 伯努利Cp 1 - (V/U)^2 for t, v, c in zip(theta, Vt, Cp): print(ftheta{np.degrees(t):6.1f} deg V{v:.2f} m/s Cp{c:.2f})这段代码只用了两个公式圆柱表面切向速度 Vθ -2U·sinθ以及沿流线的伯努利关系 Cp 1 - (V/U)²。运行结果里三个数字值得记牢θ 0°前驻点速度为零Cp 1θ 90°顶部速度达到 -2UCp -3θ 180°后驻点速度又回到零Cp 1。前后驻点压力完全对称这就是达朗贝尔佯谬——势流解算出来的圆柱阻力是零。真实圆柱当然有阻力差值全部来自粘性和分离这正是第四章存在的理由。提示如果结果里 Cp 的最小值不是 -3先检查是不是把 Vθ 的系数 2 漏了再检查 Cp 公式里有没有用错参考速度。4. 粘性、边界层与可压缩流把假设一层层拆掉的三个模块前三章的全部结论都建立在“无粘、不可压”上。这一章做两件事先把粘性放回来看方程变成什么样、壁面附近怎么处理再把可压缩性放回来看马赫数一高哪些关系会变。4.1 牛顿内摩擦定律与两类粘度粘性的定量描述就是牛顿内摩擦定律两层流体之间的剪应力与速度梯度成正比τ μ·du/dy比例系数 μ 叫动力粘度单位 Pa·s。满足这个线性关系的是牛顿流体空气和水都算血液、高分子溶液这类不满足的属于非牛顿流体。工程计算里更常用的是运动粘度 ν μ/ρ单位 m²/s因为它直接出现在雷诺数里。常温空气 μ ≈ 1.789×10⁻⁵ Pa·sν ≈ 1.46×10⁻⁵ m²/s。温度对粘度的影响方向相反气体温度升高分子热运动加剧动量交换增强粘度上升液体温度升高分子间作用力减弱粘度下降。这一条在判断“高空低温环境下粘度怎么变”时很有用。4.2 N-S 方程、雷诺数与相似准则把牛顿内摩擦定律代入动量方程就得到纳维-斯托克斯方程简称 N-S 方程。它的物理含义很直白左边是单位体积流体的惯性力右边依次是压力梯度力、粘性力和质量力。这门课不要求解 N-S 方程但要求看得懂每一项代表什么、哪一项在什么条件下能扔。扔项的判据是量级比较而量级比较的核心参数是雷诺数Re ρVL/μ VL/νRe 的物理意义是惯性力与粘性力之比。Re 很小粘性主导流动规整、粘性影响遍布全场Re 很大惯性主导粘性只在紧贴壁面的一薄层里起作用其余区域可以按无粘处理。这条判断是边界层理论的全部出发点也是风洞实验相似准则的第一条。相似准则这一块要记的是一组无量纲数Re粘性相似、Ma可压缩相似、Pr传热相似、St非定常相似。做缩比模型实验时如果所有相关无量纲数都与实物相同流动就相似但实际做不到全部相同只能抓住主要矛盾做部分相似这也是风洞数据外推的误差来源。4.3 边界层理论三个厚度与分离判据普朗特的做法是把流场分成两块远离壁面按无粘势流算紧贴壁面划出一层薄薄的粘性区叫边界层。边界层厚度 δ 随距前缘距离 x 的增长满足 δ/L 量级为 1/√ReRe 越大边界层越薄所以高雷诺数下把粘性限制在薄层里处理是合理的。对边界层方程做量级比较会得到两个关键结论。第一法向动量方程里压力沿壁面法向不变即 ∂p/∂y 0边界层内的压力由外部势流“压”进来壁面压力等于同位置外流的势流压力。第二流向动量方程里粘性项保留、其余项简化形成抛物线型的边界层方程。第二条结论决定了边界层要从前往后积分不能反过来算。三个厚度必须分清厚度定义物理意义工程用途边界层厚度 δu/U 0.99 处到壁面的距离粘性影响的粗略范围划分粘性与无粘区位移厚度 δ*∫(1 - u/U)dy因速度亏损导致流线向外偏移的量修正有效外形≈ 0.3δ 量级动量厚度 θ∫(u/U)(1 - u/U)dy动量通量的亏损量计算摩擦阻力、判断分离形状因子 H δ*/θ 用来判断边界层状态层流大约 2.6湍流大约 1.4中途减小说明转捩发生。平板边界层的转捩雷诺数常取 Re_x ≈ 5×10⁵。边界层分离的判据要记两条外部流动出现逆压梯度即 dp/dx 0壁面剪应力降到零即 ∂u/∂y|_{y0} 0。两条同时满足的位置就是分离点分离点之后流体不再贴壁形成回流区。机翼失速、汽车尾流、扩压器效率下降本质上都是分离在不同场景下的表现。这一块用“玄学”形容不算夸张因为分离点位置对雷诺数、壁面粗糙度、来流湍流度都敏感实验里经常对不上。注意分离的判据是逆压梯度加上壁面剪应力为零只看压力升高不看剪应力或者只看剪应力不看压力梯度都会把分离点判错位置。4.4 可压缩流声速、马赫数与总参数把“不可压”这条假设拆掉最先要补的就是声速和马赫数。声速是微小扰动在介质中的传播速度对完全气体有 a √(γRT)其中 γ 是比热比空气取 1.4R 是气体常数空气约 287 J/(kg·K)。马赫数是速度与当地声速之比Ma V/a注意分母是当地声速不是来流声速。马赫数分区是考试里的送分题也是判断题的高发区马赫数范围流动类型主要特征Ma 0.3不可压流密度变化小于 5%可按不可压处理0.3 ~ 0.8亚声速流密度变化明显但全场无激波0.8 ~ 1.2跨声速流局部先出现激波阻力急剧上升1.2 ~ 5超声速流出现明显的激波与膨胀波系Ma 5高超声速流气动加热、真实气体效应不可忽略等熵流动的总静参数关系是这一节的骨架T₀/T 1 (γ-1)/2 · Ma² p₀/p (T₀/T)^(γ/(γ-1)) ρ₀/ρ (T₀/T)^(1/(γ-1))总温 T₀ 是气流绝热滞止到速度为零时达到的温度总压 p₀ 是等熵滞止时达到的压强。下面这段代码把常用马赫数下的总静比算出来方便对照记忆# 等熵流总静参数比输入马赫数与比热比输出 T0/T、p0/p、rho0/rho gamma 1.4 def isentropic_ratios(Ma, gamma1.4): T0_T 1 0.5 * (gamma - 1) * Ma ** 2 # 总温比 p0_p T0_T ** (gamma / (gamma - 1)) # 总压比 rho0_rho T0_T ** (1 / (gamma - 1)) # 总密度比 return T0_T, p0_p, rho0_rho for Ma in (0.3, 0.8, 1.0, 2.0, 5.0): t, p, r isentropic_ratios(Ma) print(fMa{Ma:4.1f} T0/T{t:6.3f} p0/p{p:7.3f} rho0/rho{r:7.3f})结果里几个数值得背下来Ma 0.3 时 p₀/p ≈ 1.064密度变化只有约 6%这就是 0.3 这条不可压界线的来源Ma 1 时 T₀/T 1.2p₀/p ≈ 1.893密度比约 1.577Ma 2 时 p₀/p ≈ 7.82。把这些比值记住做题时可以先估数量级再动笔。临界状态是等熵关系里的特殊点指当地马赫数刚好等于 1 的状态用星号表示。对 γ 1.4 的空气临界参数与总参数之比是固定值T*/T₀ ≈ 0.833p*/p₀ ≈ 0.528ρ*/ρ₀ ≈ 0.634。这三个数在喷管计算里反复出现属于必须条件反射的程度。最后是激波。正激波是最简单的一种气流穿过它速度从超声速降到亚声速压强、温度、密度、熵全部上升总压下降但总温不变。总压下降这一点是判断激波存在与否的硬指标也是第五章要展开的易错点。膨胀波则相反气流经过后速度增加、压强温度下降。超声速气流绕凸角是膨胀波绕凹角是斜激波这两条规则能解释大部分超声速流动图画。5. 易漏点与易错点五个最常在考场上翻车的细节从第二章到第四章的内容连起来看真正的难点不在于公式多而在于“每个公式都有前提前提经常被默默省略”。下面五条是梳理下来最容易漏、也最容易在判卷时被扣分的地方。5.1 连续介质假设有适用边界易漏点把空气当成连续介质当成天经地义不去问它在什么条件下成立。为什么漏讲义在绪论里只给了一句定义后面所有方程都默认成立学生自然把它当成公理而不是假设。怎么补记住克努森数 Kn λ/L 这个判据Kn 0.01 才适用连续介质。高空、真空、微通道这几类场景要主动检查 Kn一旦超限N-S 方程和伯努利方程都不再成立得换稀薄气体动力学的处理方式。5.2 伯努利方程的四个前提被默默省略易漏点只要看到“p ½ρV²”就当伯努利用不管流动是不是定常、有没有粘性、是不是同一条流线。为什么漏伯努利的形式太简洁做题时容易形成条件反射看到压强和速度就往里代。怎么补每次动笔前在草稿纸角落写四个词——定常、无粘、不可压、同一条流线。四条里缺一条就不能直接用。特别是第四条有旋流场里不同流线的常数不同跨流线比较必须补上环量或者重新积分。5.3 激波前后总温总压的“一变一不变”易漏点认为激波前后的总压和总温都守恒或者都损失。为什么漏等熵流里总压总温都守恒激波本身是绝热过程容易把“绝热”和“等熵”混为一谈。怎么补抓住“绝热但不等熵”这个关键。正激波是绝热过程没有热交换所以总温不变但过程有熵增所以总压必然下降。做题时先判总压总压降了就是激波没降就是等熵压缩。5.4 位移厚度不是“速度为零的那一层”易漏点把位移厚度理解成边界层里速度为零的贴壁薄层厚度或者干脆当成边界层厚度的一部分。为什么漏名字里带“厚度”直觉上以为是某种几何尺寸。怎么补位移厚度的定义是 ∫(1 - u/U)dy本质是速度亏损沿法向的积分代表外部无粘流动被“推出去”的距离物理上相当于给物体加厚了一层。它和边界层厚度的关系是 δ* ≈ 0.3δ 量级不是同一个东西。5.5 相似准则里的特征长度取哪个易漏点算雷诺数时随便取一个长度算完发现和别人对不上。为什么漏雷诺数公式 Re VL/ν 里 L 没有规定取什么教材里不同例子取的东西不一样。怎么补特征长度必须跟研究对象匹配。平板边界层取距前缘距离 x机翼取弦长 c圆管取直径 d管道取水力直径。特征速度同样有讲究外流取来流速度内流取截面平均速度。写答案时把 L 和 V 的取法写明比结果本身更能说明你懂。6. 复习清单按“能不能自己复述”来查漏补缺看完整本讲义之后检验标准不是“看懂了”而是“合上书能不能把这条线讲出来”。下面这张表按章节把必须能回答的问题列出来每一行都可以当作一次自测。章节必须能回答的问题关联公式绪论连续介质假设成立的条件是什么相对飞行原理为什么成立Kn λ/L流体属性流体静止时为什么不能承受剪应力可压缩性用什么衡量τ μ·du/dy静力学静压强的两个特性怎么来的表压和绝对压强差在哪p p₀ ρgh运动学当地加速度和迁移加速度分别由什么引起DV/Dt ∂V/∂t (V·∇)V伯努利方程成立需要哪四个条件常数在什么情况下全场统一p ½ρV² ρgz C势流速度势和流函数存在的条件分别是什么四种基本解怎么叠加出圆柱绕流∇²φ 0L ρUΓ粘性雷诺数的物理意义是什么大 Re 和小 Re 的流场差别在哪Re ρVL/μ边界层位移厚度和动量厚度的物理意义分别是什么分离的两个判据δ*、θ、H δ*/θ可压缩Ma 0.3 这条不可压界线从哪来临界参数是多少T*/T₀ ≈ 0.833p*/p₀ ≈ 0.528具体怎么用这张表遮住右列只看中间的“必须能回答的问题”每一条先不用公式、用一段话说清楚说不出来再翻左列对应的章节。凡是要么说不清、要么说到一半要翻书的标记成重点第二天再看一遍。这种“先复述再回看”的方式比自己骗自己地重读一遍强得多。还有三个快速自检的动作值得养成。第一看到任何一个公式先问它的前提条件把条件写在公式旁边第二量纲检查两边量纲不齐一定有问题第三把结果代回极端情况验证比如 Ma 取 0 时等熵关系应该退化为 1圆柱绕流的前后驻点应该对称。我自己最开始学的时候总觉得第四章的边界层和第六章的可压缩是两块孤立的难点背公式背得挺熟一做题就发现不知道该用哪一套。后来才发现问题出在第二章的假设没吃透——每个公式都是一个“在这个假设下”的结论假设一变结论就废。从那以后我每次做题前都强制走一遍“这条流是定常吗、无粘吗、不可压吗、同一条流线吗”的清单四问过了再动笔错题率降得很明显。这份讲义的体系性正好适合这么用把六章当成一条拆假设的链条来复习比零散地背公式省力得多。希望帮到你。本文还有配套的精品资源点击获取