简介《空气动力学基础》是北京航空航天大学精品课程配套讲义由刘沛清老师主讲面向航空航天工程及相关专业学生系统梳理流体力学与空气动力学的核心知识体系。内容从绪论出发依次涵盖流体的基本属性、流体静力学与运动学、不可压缩无粘流体平面势流、粘性流体动力学基础、边界层理论及其近似、可压缩高速流动基础等模块同时引入相对飞行原理、风洞实验原理并介绍了飞机绕流、边界层分离等现象兼顾理论推导与工程应用。这份PDF讲义共1个文件压缩包大小约19.65MB内容层次分明、图表丰富适合用作课程学习笔记、考研复习提纲或自学入门的系统参考资料。目前已有674人学习浏览对于希望建立空气动力学完整知识框架、理解飞行器绕流与气动力的读者而言具有较高的参考价值。1. 空气动力学基础为什么北航这门课值得你啃 PDF提到空气动力学基础从业者第一反应往往是“伯努利方程、附面层、激波”但真正吃透的人很少。不少搞飞行器设计、CFD 仿真、风电叶片、乃至赛车尾翼的工程师工作几年后发现性能瓶颈都卡在基础概念上——升力到底是靠“流速差”还是“环量”附面层分离为什么单用雷诺数判断总出错这时候回头找一门能系统讲清楚的课北航刘沛清老师的《空气动力学基础》是绕不开的经典。这份 PDF 不是简单 PPT 堆砌它能把“会背公式”升级成“会算流动”。本文按一线应用视角拆解这门课程的知识骨架、可复现的算例流程、关键参数设置以及自学过程中最容易翻车的几个坑。适合准备进入飞行器设计、CFD 仿真、风电气动设计的新手也想帮工作几年的熟手把零散经验接回理论体系。2. 从流场描述到控制方程先搞懂这门课的底层语言空气动力学的问题本质上都在回答“给定外形和来流求受力与热流”。刘沛清老师的课程前三分之一就是在建立这套描述语言——从流线、迹线、流管这些可视化工具一路推到欧拉方程与纳维-斯托克斯方程。很多人读 PDF 时着急跳到翼型升力跳过这一部分结果后面看激波边界层干扰的图像时完全看不懂。2.1 欧拉观点与拉格朗日观点CFD 后处理里最容易混淆的视角我在做 CFD 时见过太多新人在后处理软件里问“为什么这个粒子轨迹和流线不一样”根源就是没分清欧拉观点和拉格朗日观点。刘老师在课里明确说流线Streamline是某一瞬时速度场的切线族属于欧拉描述迹线Pathline是一个流体微团随时间走过的轨迹属于拉格朗日描述。只有在定常流动中两者才重合。提示后处理软件里Paraview 画 Stream Tracer 要看“当前瞬时场”而 FLUENT 的 Pathlines粒子追踪数值上就是迹线。如果计算是瞬态的两线必然分离。判断一个近似是否可用标准是 Strouhal 数St fL/U。当特征频率 f 趋近于 0定常时St→0流线与迹线才等价。这个参数在旋翼、涡轮机械的仿真里必须检查如果稳态假设下 St 过高流线图会误导你对分离区的判断。实际项目中我一般先确认旋翼转速对应的叶尖频率是否远低于来流特征频率否则就直接切到瞬态计算别省那点算力。2.2 控制方程的无量纲化三个必须记死的特征数CFD 工程师的判断力很大程度来自对无量纲数的直觉。刘老师的课把方程无量纲化的目标讲得极清楚消除单位影响让相似流动归并。实际操作中你至少要把下面三个特征数刻进概念里无量纲数定义物理意义常见量级参考雷诺数 ReρUL/μ惯性力/粘性力民用飞机机翼 1e6~1e7小型无人机 1e5马赫数 MaU/a流速/当地声速巡航客机 0.78~0.85低速风洞 0.3克努森数 Knλ/L分子自由程/特征长度标准大气低空 1e-6高空稀薄流 0.01 以上实战里我有个习惯拿到一个新几何模型先算 Re 和 Ma判断流动是层流、湍流、还是可压缩区域。Re 决定了附面层状态层流转捩在工程上通常取 Re≈5e5Ma 决定能不能用不可压缩假设。刘老师在课上强调的“连续介质假设失效”条件对应 Kn0.01在高空高超声速飞行器设计中是选 DSMC 还是 NS 求解器的门限这个在一般工程教程里很少讲透。2.3 从欧拉方程到位函数低速计算的最短路如果只做低速不可压缩流动Ma0.3连续方程简化为 ∇·V0再配合无旋假设∇×V0速度场可以写成 V∇φ整个问题塌缩成一个拉普拉斯方程 ∇²φ0。这就是刘老师课里反复强调的位函数方法——它让手动估算和早期 CFD 都成了可能。我一般会用它做二维翼型升力的快速估算在升力线理论之前先用叠加法均匀流源汇偶极子算对称翼型再用库塔条件定环量。这样算出的升力系数在小迎角5°下有不错的量级误差约 5%~10%用来预判 CFD 初场很省计算资源。把位函数方程离散成有限差分格式在网格上跑一遍高斯-赛德尔迭代不需要任何商业软件一张二维方格网就能验证“绕过圆柱有升力库塔-儒可夫斯基定理”这个反直觉结论——这个验证演示了环量产生的本质。3. 附面层与粘性流动把“边界层分离”钉死在工程直觉里附面层这一章是整个课程的枢纽也是工程里吃掉最多 debug 时间的地方。无粘位势流动算不出阻力达朗伯悖论只有引入粘性边界层才能解释压差阻力和摩擦阻力。刘老师的讲法是从普朗特 1904 年的思想出发高雷诺数流动中粘性只在壁面附近薄层内重要把外流看成有滑移的无粘流把层内看成剪切流两者通过边界层边缘速度 Ue 衔接。3.1 层流附面层动量积分用手算也能求出平板摩擦阻力工程里有时不需要立刻开大涡模拟而是先用动量积分关系式判断趋势。这是刘老师在课程中单列的重点von Kármán 动量积分方程dθ/dx (2H) θ/Ue · dUe/dx Cf/2其中 θ 是动量厚度H 为形状因子位移厚度 δ*/动量厚度 θCf 为壁面摩擦系数。对零压力梯度平板取层流自相似剖面Blasius 解得到 δ/x ≈ 5.0/√Re_xCf 0.664/√Re_x。这段我经常用来给年轻同事做“手算前理解”机翼前缘雷诺数 Re_x 对应位置如果 Re_x 从 1e5 变成 1e6层流附面层厚度减小 √10 倍摩擦阻力系数也减小 √10 倍。这个量级直觉在后处理核验 CFD 残差时很管用——如果你看到平板摩擦阻力比理论值大一个量级先怀疑网格 y 或者湍流模型而不是怀疑物理模型。3.2 转捩判断很多 CFD 算不准升阻比的元凶层流向湍流转捩的判据工程上最常用的是临界雷诺数。但直接拿 Re_cr5e5 套所有流动就是最常见的翻车现场平板自然转捩 Re_x≈5e5 到 1e6有逆压梯度的翼型可能提前到 3e5而表面粗糙或自由流湍流度高可能掉到 1e5 以下。刘老师的课给出两个稳定的工程抓手动量厚度雷诺数 Re_θ ≈ 200~400自然转捩形状因子 H 达到 2.6 附近时层流剖面趋于分离实际算翼型时我一般先看失速迎角附近的压力分布。如果层流分离泡出现Re_θ 判据会被局部逆压梯度破坏这时需要用 e^N 方法配合线性稳定性理论。在商业 CFD 里就是打开转捩模型给自由流湍流度一般取 0.5%~2%和转捩雷诺数这样翼型升力系数的攻角曲线才不会比实验提前失速。3.3 分离与失速为什么大迎角下升力突然掉附面层分离是逆压梯度造成的压力沿流向增加dP/dx0近壁流体动能被不断消耗最终在壁面速度梯度降为零并形成回流。失速的本质不是来流“冲角太大”而是分离区扩展导致环量衰减。刘老师课里用圆柱绕流的分离点位置讲了几种状态亚临界 Re约 1e5分离点约在从前驻点算起 80° 的位置超临界 Re 约在 120°分离点后移阻力系数出现“阻力危机”。实战中我用这个来理解汽车外流场车顶后缘的分离位置决定尾部低压区大小直接影响气动阻力。想延迟分离的常见手段是控制逆压梯度翼型上表面后段设计成“缓压恢复”汽车尾部加小鸭尾或增加尾部上翘角本质就是给分离区补一个压力恢复缓冲。这些设计反馈在刘老师的“分离控制”一节里都有物理依据不是单纯造型美学。4. 可压缩流动与激波高亚声速和超声速绕不开的修正民航客机巡航 Ma≈0.8军机超声速火箭更是高超声速全覆盖。刘老师的可压缩流动章节提供了从亚声速到高超声速的统一框架。工程上要抓的核心当地 Ma1 时会发生什么激波两侧参数如何突变。4.1 正激波关系式记住总压损失与熵增的非单调性超音速气流通过正激波必然减速到亚音速M21总压下降最厉害。刘老师给出了完整关系式组工程上最关心的是总压恢复比p02/p01 [ (γ1)M1² / ((γ-1)M1²2) ]^(γ/(γ-1)) · [ (γ1) / (2γM1²-(γ-1)) ]^(1/(γ-1))当 M11.5总压恢复约 0.93M12.0约 0.72M13.0约 0.33。这个衰减趋势决定了一台超声速进气道在 Ma2.5 后设计斜率要压得低否则总压损失直接吃掉发动机可用功。我自己做进气道设计时快速验证进气道总压恢复的经验是 1 - 0.075(M1-1)^1.35这个经验公式跟正激波解析解对得上。4.2 膨胀波与普朗特-迈耶函数超声速外流的“降压”通道超声速绕凸角流动不是激波而是通过一个连续的膨胀扇Prandtl-Meyer expansion让气流加速降压。刘老师这一节我推荐大家务必把 PM 函数表或函数关系代码存一份ν(M) √(γ1/γ-1) · arctan√(γ-1/γ1 (M²-1)) - arctan√(M²-1)实际设计超声速风洞喷管或火箭喷管这个函数是计算出口马赫数的核心。做二维超声速喷管型面时从喉部 Ma1 到出口 Ma3通过 PM 反函数得到转折角型面就按等膨胀波法画出来。不少 CFD 从业者觉得超声速很难其实第一步的特征线法在这里就已经可以脱离网格估算气流角了用来验证 CFD 结果非常方便。4.3 高超声速的牛顿冲击理论为什么分子效应最终占上风当 Ma5激波贴近壁面流场可以用“自由分子碰撞”视角简化——激波层极薄壁面压力近似 p ρ∞V∞² sin²θ牛顿正弦平方律。刘老师提到高超声速相似参数以马赫数无关原理组合出现Ma 难以单独起作用。实际做高超声速飞行器再入时激波层内温度极高引起振动激发、离解、电离比热比 γ 不再是常数这时量热完全气体假设失效NS 求解器要启用化学非平衡。所以高超声速 CFD 往往把门槛放在了热力学模型而不是网格上。5. 翼型与机翼气动特性把课程知识转成一张可用参数表学完了控制方程和流动现象最后一里路是落到翼型与机翼几何及气动数据的生成。刘老师的课从薄翼理论、有限翼展的升力线理论讲起再转到黏性修正。这一章按工程做法给出可复现的流程。5.1 用薄翼理论快速估算零升迎角与升力线斜率薄翼理论结论理想流体中小迎角线性升力升力线斜率约 2π/rad约 0.11/°升力系数的零迎角由弯度决定。对称翼型零升迎角为 0°带弯度翼型零升迎角为负值约 -2° 量级。工程上很多快速打样设计先用这个理论确认升力梯度的量级如果后续 CFD 算出的升力线斜率明显偏离 6.28/rad比如小于 5.5检查网格或湍流模型如果大于 7那可能边界条件或远场尺寸有问题。这是很好的 CFD 结果自检锚点。5.2 用升力线理论估有限翼展诱导阻力从二维到三维的桥有限翼展机翼的升力对应一个展向环量分布尾涡系诱导下洗翼型有效迎角减小产生诱导阻力。修正公式CL CL0 · (1 2/(π AR e))^(-1)近似对椭圆升力分布 e1CDi CL² / (π AR e)这里展弦比 AR 是翼展除以平均气动弦e 是奥斯瓦尔德效率因子典型 0.75~0.85。这门课教的核心是“椭圆升力分布有最小诱导阻力”——对应展向弦长椭圆变化的机翼。现代飞机不是椭圆翼但通过扭转和平面形状逼近椭圆载荷。我一般建议在概念设计时做一张展向升力分布图看是否接近椭圆若翼尖载荷过高则诱导阻力大还容易翼尖失速。5.3 把理论转成作业参数一个二维翼型算例的参数表下面是一套可复用的低雷诺数翼型比如 S1223 或 NACA4415二维 CFD 参数配合主流求解器跑一遍参数名推荐值说明来流速度20 m/sMa0.3不可压缩假设有效弦长基准0.1 m用于雷诺数定义雷诺数1.35e5对应低速无人机翼型实验攻角范围0°~15°步长 1°覆盖线性段到失速后湍流模型SST k-ω逆压梯度分离预测优于 k-ε近壁 y1需要加密附面层网格入口湍流度0.5%风洞低湍流环境远场距离20 倍弦长避免边界影响网格生成时附面层第一层高度估算用 y1 y · μ / (ρ · uτ)但 uτ 不知道可先用平板公式 uτ U∞ √(Cf/2)Cf0.664/√Re给一个初始值后迭代加密。y1 往往意味着第一层网格厚度在 1e-5 到 1e-4 倍弦长量级。5.4 从翼型数据到三维外形的映射一份低成本校验流程拿到二维极曲线后三维有限翼展直接用升力线、升力面或涡格法做初校。我习惯流程是先二维 CFD翼型雷诺数匹配得到极曲线再用涡格法如 Tornado给三维机翼加载输出 CL 对 α、CDi最后对重点设计点做一次三维 RANS 校核。这套流程的成本二维 CFD 一晚上可跑完一个极曲线具体看网格量涡格法是秒出。很多做轻型飞机、无人机的案头计算都这么来的。6. 进阶与自测用一小段脚本验证 PM 函数和附面层关系式最后一章给一个不依赖商业软件、可复现的公式验证路径。把 PM 函数和层流平板附面层公式写进 Python立刻检验你的理解是否到位也方便以后做快速设计估算。6.1 普朗特-迈耶函数与正激波关系的数值验证脚本import math def pm_function(M, gamma1.4): 普朗特-迈耶角输入马赫数输出偏转角弧度 sqrt_term math.sqrt((gamma 1) / (gamma - 1)) term1 sqrt_term * math.atan(math.sqrt((gamma - 1) / (gamma 1) * (M**2 - 1))) term2 math.atan(math.sqrt(M**2 - 1)) return term1 - term2 def normal_shock_total_pressure_ratio(M1, gamma1.4): 正激波总压恢复比 term1 ((gamma 1) * M1**2 / ((gamma - 1) * M1**2 2)) ** (gamma / (gamma - 1)) term2 ((gamma 1) / (2 * gamma * M1**2 - (gamma - 1))) ** (1 / (gamma - 1)) return term1 * term2 # 验证M2时PM角应约26.38°正激波总压恢复约0.72 for M in [1.2, 1.5, 2.0, 3.0]: pm_deg math.degrees(pm_function(M)) p_ratio normal_shock_total_pressure_ratio(M) print(fM{M:4.2f} PM角{pm_deg:6.2f}° 正激波总压比{p_ratio:.3f})这段脚本逻辑很直白pm_function 里两项反正切的差就是气流从声速加速到 M 所需偏转角normal_shock 函数里两条幂指数分别对应激波上下游参数组合。跑一下你会看到数值和课本气动表吻合。自己把 γ 改成 1.3高温空气看 PM 角下降就能理解为什么高超声速喷管设计给定总温后出口马赫数计算不能用常温空气参数。6.2 层流平板附面层数据的自检清单验证完超声速公式再用最小代码检查不可压缩层流平板数据import math def blasius_flat_plate(x, Re_x): Blasius层流平板边界层厚度、位移厚度、动量厚度、摩擦系数 delta 5.0 * x / math.sqrt(Re_x) delta_star 1.7208 * x / math.sqrt(Re_x) theta 0.664 * x / math.sqrt(Re_x) cf 0.664 / math.sqrt(Re_x) return delta, delta_star, theta, cf x 1.0 # 平板特征位置 1m U 20.0 # 流速 20 m/s nu 1.5e-5 # 空气运动粘度 Re_x U * x / nu print(fRe_x{Re_x:.2e}) delta, delta_star, theta, cf blasius_flat_plate(x, Re_x) print(f边界层厚度{delta:.4f} m) print(f位移厚度{delta_star:.4f} m) print(f动量厚度{theta:.4f} m) print(f摩擦系数Cf{cf:.4f})这个脚本每次做 CFD 遇到平板类算例我都拿来对比摩擦阻力。如果你的三维平板算例算出的平均 Cf 比理论值偏差 20% 以上先查 y——这几乎已经成为我排查气动仿真问题的第一动作。看位移厚度时记住一个工程经验位移厚度相当于把物体外形向外“推”了一个 δ*实际影响压力分布。这也是为什么风洞试验模型加工需要扣除附面层影响模型与实物的转捩模拟往往靠贴转捩带实现。6.3 自测检验你有没有真正入门我给自己和新人留三道题能完整答出来这门课就算真过了为什么对称翼型在 Ma0.2 和 Ma0.7 下的升力线斜率差很多提示可压缩修正用到普朗特-格劳厄特法则但激波出现后又失效附面层分离对翼型失速攻角的影响为什么说重点在于后缘压力恢复而不是前缘吸力峰值一个 Re1e5 的圆柱绕流不加湍流模型直接层流 NS 算出的阻力系数与实验约 1.2差多少原因是什么回答第三题时很多人会翻车在没注意亚临界圆柱绕流本身已经包含周期性涡脱落层流 NS 也能捕捉但二维计算天然无法重现三维涡结构所以阻力通常偏高 30% 以上。这类边界认识正是读刘老师 PDF 比单纯刷软件教程多获得的回报。希望帮到你。本文还有配套的精品资源点击获取 空气动力学基础PDF精讲：从附面层到CFD仿真的工程实战指南
简介《空气动力学基础》是北京航空航天大学精品课程配套讲义由刘沛清老师主讲面向航空航天工程及相关专业学生系统梳理流体力学与空气动力学的核心知识体系。内容从绪论出发依次涵盖流体的基本属性、流体静力学与运动学、不可压缩无粘流体平面势流、粘性流体动力学基础、边界层理论及其近似、可压缩高速流动基础等模块同时引入相对飞行原理、风洞实验原理并介绍了飞机绕流、边界层分离等现象兼顾理论推导与工程应用。这份PDF讲义共1个文件压缩包大小约19.65MB内容层次分明、图表丰富适合用作课程学习笔记、考研复习提纲或自学入门的系统参考资料。目前已有674人学习浏览对于希望建立空气动力学完整知识框架、理解飞行器绕流与气动力的读者而言具有较高的参考价值。1. 空气动力学基础为什么北航这门课值得你啃 PDF提到空气动力学基础从业者第一反应往往是“伯努利方程、附面层、激波”但真正吃透的人很少。不少搞飞行器设计、CFD 仿真、风电叶片、乃至赛车尾翼的工程师工作几年后发现性能瓶颈都卡在基础概念上——升力到底是靠“流速差”还是“环量”附面层分离为什么单用雷诺数判断总出错这时候回头找一门能系统讲清楚的课北航刘沛清老师的《空气动力学基础》是绕不开的经典。这份 PDF 不是简单 PPT 堆砌它能把“会背公式”升级成“会算流动”。本文按一线应用视角拆解这门课程的知识骨架、可复现的算例流程、关键参数设置以及自学过程中最容易翻车的几个坑。适合准备进入飞行器设计、CFD 仿真、风电气动设计的新手也想帮工作几年的熟手把零散经验接回理论体系。2. 从流场描述到控制方程先搞懂这门课的底层语言空气动力学的问题本质上都在回答“给定外形和来流求受力与热流”。刘沛清老师的课程前三分之一就是在建立这套描述语言——从流线、迹线、流管这些可视化工具一路推到欧拉方程与纳维-斯托克斯方程。很多人读 PDF 时着急跳到翼型升力跳过这一部分结果后面看激波边界层干扰的图像时完全看不懂。2.1 欧拉观点与拉格朗日观点CFD 后处理里最容易混淆的视角我在做 CFD 时见过太多新人在后处理软件里问“为什么这个粒子轨迹和流线不一样”根源就是没分清欧拉观点和拉格朗日观点。刘老师在课里明确说流线Streamline是某一瞬时速度场的切线族属于欧拉描述迹线Pathline是一个流体微团随时间走过的轨迹属于拉格朗日描述。只有在定常流动中两者才重合。提示后处理软件里Paraview 画 Stream Tracer 要看“当前瞬时场”而 FLUENT 的 Pathlines粒子追踪数值上就是迹线。如果计算是瞬态的两线必然分离。判断一个近似是否可用标准是 Strouhal 数St fL/U。当特征频率 f 趋近于 0定常时St→0流线与迹线才等价。这个参数在旋翼、涡轮机械的仿真里必须检查如果稳态假设下 St 过高流线图会误导你对分离区的判断。实际项目中我一般先确认旋翼转速对应的叶尖频率是否远低于来流特征频率否则就直接切到瞬态计算别省那点算力。2.2 控制方程的无量纲化三个必须记死的特征数CFD 工程师的判断力很大程度来自对无量纲数的直觉。刘老师的课把方程无量纲化的目标讲得极清楚消除单位影响让相似流动归并。实际操作中你至少要把下面三个特征数刻进概念里无量纲数定义物理意义常见量级参考雷诺数 ReρUL/μ惯性力/粘性力民用飞机机翼 1e6~1e7小型无人机 1e5马赫数 MaU/a流速/当地声速巡航客机 0.78~0.85低速风洞 0.3克努森数 Knλ/L分子自由程/特征长度标准大气低空 1e-6高空稀薄流 0.01 以上实战里我有个习惯拿到一个新几何模型先算 Re 和 Ma判断流动是层流、湍流、还是可压缩区域。Re 决定了附面层状态层流转捩在工程上通常取 Re≈5e5Ma 决定能不能用不可压缩假设。刘老师在课上强调的“连续介质假设失效”条件对应 Kn0.01在高空高超声速飞行器设计中是选 DSMC 还是 NS 求解器的门限这个在一般工程教程里很少讲透。2.3 从欧拉方程到位函数低速计算的最短路如果只做低速不可压缩流动Ma0.3连续方程简化为 ∇·V0再配合无旋假设∇×V0速度场可以写成 V∇φ整个问题塌缩成一个拉普拉斯方程 ∇²φ0。这就是刘老师课里反复强调的位函数方法——它让手动估算和早期 CFD 都成了可能。我一般会用它做二维翼型升力的快速估算在升力线理论之前先用叠加法均匀流源汇偶极子算对称翼型再用库塔条件定环量。这样算出的升力系数在小迎角5°下有不错的量级误差约 5%~10%用来预判 CFD 初场很省计算资源。把位函数方程离散成有限差分格式在网格上跑一遍高斯-赛德尔迭代不需要任何商业软件一张二维方格网就能验证“绕过圆柱有升力库塔-儒可夫斯基定理”这个反直觉结论——这个验证演示了环量产生的本质。3. 附面层与粘性流动把“边界层分离”钉死在工程直觉里附面层这一章是整个课程的枢纽也是工程里吃掉最多 debug 时间的地方。无粘位势流动算不出阻力达朗伯悖论只有引入粘性边界层才能解释压差阻力和摩擦阻力。刘老师的讲法是从普朗特 1904 年的思想出发高雷诺数流动中粘性只在壁面附近薄层内重要把外流看成有滑移的无粘流把层内看成剪切流两者通过边界层边缘速度 Ue 衔接。3.1 层流附面层动量积分用手算也能求出平板摩擦阻力工程里有时不需要立刻开大涡模拟而是先用动量积分关系式判断趋势。这是刘老师在课程中单列的重点von Kármán 动量积分方程dθ/dx (2H) θ/Ue · dUe/dx Cf/2其中 θ 是动量厚度H 为形状因子位移厚度 δ*/动量厚度 θCf 为壁面摩擦系数。对零压力梯度平板取层流自相似剖面Blasius 解得到 δ/x ≈ 5.0/√Re_xCf 0.664/√Re_x。这段我经常用来给年轻同事做“手算前理解”机翼前缘雷诺数 Re_x 对应位置如果 Re_x 从 1e5 变成 1e6层流附面层厚度减小 √10 倍摩擦阻力系数也减小 √10 倍。这个量级直觉在后处理核验 CFD 残差时很管用——如果你看到平板摩擦阻力比理论值大一个量级先怀疑网格 y 或者湍流模型而不是怀疑物理模型。3.2 转捩判断很多 CFD 算不准升阻比的元凶层流向湍流转捩的判据工程上最常用的是临界雷诺数。但直接拿 Re_cr5e5 套所有流动就是最常见的翻车现场平板自然转捩 Re_x≈5e5 到 1e6有逆压梯度的翼型可能提前到 3e5而表面粗糙或自由流湍流度高可能掉到 1e5 以下。刘老师的课给出两个稳定的工程抓手动量厚度雷诺数 Re_θ ≈ 200~400自然转捩形状因子 H 达到 2.6 附近时层流剖面趋于分离实际算翼型时我一般先看失速迎角附近的压力分布。如果层流分离泡出现Re_θ 判据会被局部逆压梯度破坏这时需要用 e^N 方法配合线性稳定性理论。在商业 CFD 里就是打开转捩模型给自由流湍流度一般取 0.5%~2%和转捩雷诺数这样翼型升力系数的攻角曲线才不会比实验提前失速。3.3 分离与失速为什么大迎角下升力突然掉附面层分离是逆压梯度造成的压力沿流向增加dP/dx0近壁流体动能被不断消耗最终在壁面速度梯度降为零并形成回流。失速的本质不是来流“冲角太大”而是分离区扩展导致环量衰减。刘老师课里用圆柱绕流的分离点位置讲了几种状态亚临界 Re约 1e5分离点约在从前驻点算起 80° 的位置超临界 Re 约在 120°分离点后移阻力系数出现“阻力危机”。实战中我用这个来理解汽车外流场车顶后缘的分离位置决定尾部低压区大小直接影响气动阻力。想延迟分离的常见手段是控制逆压梯度翼型上表面后段设计成“缓压恢复”汽车尾部加小鸭尾或增加尾部上翘角本质就是给分离区补一个压力恢复缓冲。这些设计反馈在刘老师的“分离控制”一节里都有物理依据不是单纯造型美学。4. 可压缩流动与激波高亚声速和超声速绕不开的修正民航客机巡航 Ma≈0.8军机超声速火箭更是高超声速全覆盖。刘老师的可压缩流动章节提供了从亚声速到高超声速的统一框架。工程上要抓的核心当地 Ma1 时会发生什么激波两侧参数如何突变。4.1 正激波关系式记住总压损失与熵增的非单调性超音速气流通过正激波必然减速到亚音速M21总压下降最厉害。刘老师给出了完整关系式组工程上最关心的是总压恢复比p02/p01 [ (γ1)M1² / ((γ-1)M1²2) ]^(γ/(γ-1)) · [ (γ1) / (2γM1²-(γ-1)) ]^(1/(γ-1))当 M11.5总压恢复约 0.93M12.0约 0.72M13.0约 0.33。这个衰减趋势决定了一台超声速进气道在 Ma2.5 后设计斜率要压得低否则总压损失直接吃掉发动机可用功。我自己做进气道设计时快速验证进气道总压恢复的经验是 1 - 0.075(M1-1)^1.35这个经验公式跟正激波解析解对得上。4.2 膨胀波与普朗特-迈耶函数超声速外流的“降压”通道超声速绕凸角流动不是激波而是通过一个连续的膨胀扇Prandtl-Meyer expansion让气流加速降压。刘老师这一节我推荐大家务必把 PM 函数表或函数关系代码存一份ν(M) √(γ1/γ-1) · arctan√(γ-1/γ1 (M²-1)) - arctan√(M²-1)实际设计超声速风洞喷管或火箭喷管这个函数是计算出口马赫数的核心。做二维超声速喷管型面时从喉部 Ma1 到出口 Ma3通过 PM 反函数得到转折角型面就按等膨胀波法画出来。不少 CFD 从业者觉得超声速很难其实第一步的特征线法在这里就已经可以脱离网格估算气流角了用来验证 CFD 结果非常方便。4.3 高超声速的牛顿冲击理论为什么分子效应最终占上风当 Ma5激波贴近壁面流场可以用“自由分子碰撞”视角简化——激波层极薄壁面压力近似 p ρ∞V∞² sin²θ牛顿正弦平方律。刘老师提到高超声速相似参数以马赫数无关原理组合出现Ma 难以单独起作用。实际做高超声速飞行器再入时激波层内温度极高引起振动激发、离解、电离比热比 γ 不再是常数这时量热完全气体假设失效NS 求解器要启用化学非平衡。所以高超声速 CFD 往往把门槛放在了热力学模型而不是网格上。5. 翼型与机翼气动特性把课程知识转成一张可用参数表学完了控制方程和流动现象最后一里路是落到翼型与机翼几何及气动数据的生成。刘老师的课从薄翼理论、有限翼展的升力线理论讲起再转到黏性修正。这一章按工程做法给出可复现的流程。5.1 用薄翼理论快速估算零升迎角与升力线斜率薄翼理论结论理想流体中小迎角线性升力升力线斜率约 2π/rad约 0.11/°升力系数的零迎角由弯度决定。对称翼型零升迎角为 0°带弯度翼型零升迎角为负值约 -2° 量级。工程上很多快速打样设计先用这个理论确认升力梯度的量级如果后续 CFD 算出的升力线斜率明显偏离 6.28/rad比如小于 5.5检查网格或湍流模型如果大于 7那可能边界条件或远场尺寸有问题。这是很好的 CFD 结果自检锚点。5.2 用升力线理论估有限翼展诱导阻力从二维到三维的桥有限翼展机翼的升力对应一个展向环量分布尾涡系诱导下洗翼型有效迎角减小产生诱导阻力。修正公式CL CL0 · (1 2/(π AR e))^(-1)近似对椭圆升力分布 e1CDi CL² / (π AR e)这里展弦比 AR 是翼展除以平均气动弦e 是奥斯瓦尔德效率因子典型 0.75~0.85。这门课教的核心是“椭圆升力分布有最小诱导阻力”——对应展向弦长椭圆变化的机翼。现代飞机不是椭圆翼但通过扭转和平面形状逼近椭圆载荷。我一般建议在概念设计时做一张展向升力分布图看是否接近椭圆若翼尖载荷过高则诱导阻力大还容易翼尖失速。5.3 把理论转成作业参数一个二维翼型算例的参数表下面是一套可复用的低雷诺数翼型比如 S1223 或 NACA4415二维 CFD 参数配合主流求解器跑一遍参数名推荐值说明来流速度20 m/sMa0.3不可压缩假设有效弦长基准0.1 m用于雷诺数定义雷诺数1.35e5对应低速无人机翼型实验攻角范围0°~15°步长 1°覆盖线性段到失速后湍流模型SST k-ω逆压梯度分离预测优于 k-ε近壁 y1需要加密附面层网格入口湍流度0.5%风洞低湍流环境远场距离20 倍弦长避免边界影响网格生成时附面层第一层高度估算用 y1 y · μ / (ρ · uτ)但 uτ 不知道可先用平板公式 uτ U∞ √(Cf/2)Cf0.664/√Re给一个初始值后迭代加密。y1 往往意味着第一层网格厚度在 1e-5 到 1e-4 倍弦长量级。5.4 从翼型数据到三维外形的映射一份低成本校验流程拿到二维极曲线后三维有限翼展直接用升力线、升力面或涡格法做初校。我习惯流程是先二维 CFD翼型雷诺数匹配得到极曲线再用涡格法如 Tornado给三维机翼加载输出 CL 对 α、CDi最后对重点设计点做一次三维 RANS 校核。这套流程的成本二维 CFD 一晚上可跑完一个极曲线具体看网格量涡格法是秒出。很多做轻型飞机、无人机的案头计算都这么来的。6. 进阶与自测用一小段脚本验证 PM 函数和附面层关系式最后一章给一个不依赖商业软件、可复现的公式验证路径。把 PM 函数和层流平板附面层公式写进 Python立刻检验你的理解是否到位也方便以后做快速设计估算。6.1 普朗特-迈耶函数与正激波关系的数值验证脚本import math def pm_function(M, gamma1.4): 普朗特-迈耶角输入马赫数输出偏转角弧度 sqrt_term math.sqrt((gamma 1) / (gamma - 1)) term1 sqrt_term * math.atan(math.sqrt((gamma - 1) / (gamma 1) * (M**2 - 1))) term2 math.atan(math.sqrt(M**2 - 1)) return term1 - term2 def normal_shock_total_pressure_ratio(M1, gamma1.4): 正激波总压恢复比 term1 ((gamma 1) * M1**2 / ((gamma - 1) * M1**2 2)) ** (gamma / (gamma - 1)) term2 ((gamma 1) / (2 * gamma * M1**2 - (gamma - 1))) ** (1 / (gamma - 1)) return term1 * term2 # 验证M2时PM角应约26.38°正激波总压恢复约0.72 for M in [1.2, 1.5, 2.0, 3.0]: pm_deg math.degrees(pm_function(M)) p_ratio normal_shock_total_pressure_ratio(M) print(fM{M:4.2f} PM角{pm_deg:6.2f}° 正激波总压比{p_ratio:.3f})这段脚本逻辑很直白pm_function 里两项反正切的差就是气流从声速加速到 M 所需偏转角normal_shock 函数里两条幂指数分别对应激波上下游参数组合。跑一下你会看到数值和课本气动表吻合。自己把 γ 改成 1.3高温空气看 PM 角下降就能理解为什么高超声速喷管设计给定总温后出口马赫数计算不能用常温空气参数。6.2 层流平板附面层数据的自检清单验证完超声速公式再用最小代码检查不可压缩层流平板数据import math def blasius_flat_plate(x, Re_x): Blasius层流平板边界层厚度、位移厚度、动量厚度、摩擦系数 delta 5.0 * x / math.sqrt(Re_x) delta_star 1.7208 * x / math.sqrt(Re_x) theta 0.664 * x / math.sqrt(Re_x) cf 0.664 / math.sqrt(Re_x) return delta, delta_star, theta, cf x 1.0 # 平板特征位置 1m U 20.0 # 流速 20 m/s nu 1.5e-5 # 空气运动粘度 Re_x U * x / nu print(fRe_x{Re_x:.2e}) delta, delta_star, theta, cf blasius_flat_plate(x, Re_x) print(f边界层厚度{delta:.4f} m) print(f位移厚度{delta_star:.4f} m) print(f动量厚度{theta:.4f} m) print(f摩擦系数Cf{cf:.4f})这个脚本每次做 CFD 遇到平板类算例我都拿来对比摩擦阻力。如果你的三维平板算例算出的平均 Cf 比理论值偏差 20% 以上先查 y——这几乎已经成为我排查气动仿真问题的第一动作。看位移厚度时记住一个工程经验位移厚度相当于把物体外形向外“推”了一个 δ*实际影响压力分布。这也是为什么风洞试验模型加工需要扣除附面层影响模型与实物的转捩模拟往往靠贴转捩带实现。6.3 自测检验你有没有真正入门我给自己和新人留三道题能完整答出来这门课就算真过了为什么对称翼型在 Ma0.2 和 Ma0.7 下的升力线斜率差很多提示可压缩修正用到普朗特-格劳厄特法则但激波出现后又失效附面层分离对翼型失速攻角的影响为什么说重点在于后缘压力恢复而不是前缘吸力峰值一个 Re1e5 的圆柱绕流不加湍流模型直接层流 NS 算出的阻力系数与实验约 1.2差多少原因是什么回答第三题时很多人会翻车在没注意亚临界圆柱绕流本身已经包含周期性涡脱落层流 NS 也能捕捉但二维计算天然无法重现三维涡结构所以阻力通常偏高 30% 以上。这类边界认识正是读刘老师 PDF 比单纯刷软件教程多获得的回报。希望帮到你。本文还有配套的精品资源点击获取