✅作者简介：热爱科研的Matlab仿真开发者，擅长数学建模、数据处理、建模仿真、程序设计、完整代码获取、论文复现及科研仿真。
🍎 往期回顾关注个人主页：Matlab科研工作室
👇 关注我领取海量matlab电子书和数学建模资料
🍊个人信条：格物致知,完整Matlab代码获取及仿真咨询内容私信。
🔥 内容介绍
跨声速工况下翼型的临界马赫数直接决定飞行器的巡航效率与激波出现阈值，是高速无人机、亚声速民机气动设计的核心性能指标。本研究以临界马赫数最大化为核心目标，采用CST（类形状函数变换）方法完成翼型几何的参数化光滑建模，结合二阶涡旋面板法实现翼面流场的快速求解，通过卡门-钱学森公式完成亚声速可压缩性修正，最终引入改进自适应粒子群算法开展全局约束寻优，构建出一套低计算成本、高可靠性的翼型优化全流程。研究以RAE2822超临界基准翼型为初始样本，在约束翼型最大相对厚度不低于12%、设计升力系数不小于0.6的前提下，优化后翼型的临界马赫数从初始的0.73提升至0.784，相对增幅达7.4%，跨声速阻力发散马赫数同步提升6.3%。该方法无需依赖大规模CFD网格生成，单步翼型流场求解耗时控制在1秒以内，可在个人计算平台上快速完成全流程优化，为翼型方案的初期快速迭代提供了高效可行的技术路径。
关键词：CST参数化；二阶涡旋面板法；临界马赫数；自适应粒子群优化；亚声速气动设计
一、引言
在现代航空装备的亚声速巡航与跨声速突防场景中，翼面局部超声速区的出现会直接引发激波阻力骤增、边界层分离、操纵性恶化等一系列问题，而临界马赫数正是翼面首次出现局部声速时对应的自由来流马赫数，是衡量翼型跨声速气动潜力的核心指标。长期以来，传统翼型设计依赖大量风洞试验与高精度RANS数值模拟，虽然计算精度高，但迭代周期长、计算成本高昂，难以支撑多变量全局优化的快速探索需求。
近年来，翼型参数化方法与快速流场求解技术的结合，为低成本高性能翼型设计提供了新的思路。其中CST参数化方法凭借少量控制变量即可实现翼型几何的高精度连续描述，有效避免了传统样条参数化中容易出现的几何不光滑、约束能力弱的问题；而面板法基于势流理论构建求解框架，无需生成三维体网格，仅通过翼面离散面元即可快速得到亚声速无黏流场的速度与压力分布，计算效率较传统CFD方法提升两个数量级以上。在此基础上引入智能优化算法，可在短时间内完成多变量设计空间的全局探索，快速定位临界马赫数最大化的翼型构型。本研究即围绕这一技术链路展开系统性研究，验证该方法在工程翼型优化中的实用性与有效性。
二、核心理论与方法体系
2.1 临界马赫数的物理判定逻辑
临界马赫数的物理本质是翼面局部流速达到当地声速时的自由来流状态。对于二维亚声速定常流动，不可压缩流条件下的翼面压力系数可通过势流理论直接求解，而亚声速可压缩性的影响可通过卡门-钱学森“切线气体近似”完成精准修正，该方法通过构建假想不可压缩流与真实可压缩流之间的速度映射关系，无需求解复杂的非线性速度势方程，即可在全亚声速区间得到误差远小于普朗特-格劳厄特修正的压力分布结果。
本研究中临界马赫数的计算采用迭代扫描法：首先固定来流迎角，从低马赫数开始逐步抬升自由来流速度，每一步通过面板法求解流场并经可压缩性修正得到全翼面压力分布，逐点计算翼面各位置的局部马赫数，当翼面任意点的局部马赫数首次达到1时，对应的自由来流马赫数即为该翼型的临界马赫数。该判定过程物理边界清晰，计算稳定性高，完全适配优化迭代中的批量计算需求。
2.2 CST翼型参数化建模原理
本研究采用6阶CST参数化方案，翼型上下表面各分配7个独立的形状函数控制系数，合计14个设计变量，仅需这组变量即可精准重构RAE2822、NACA0012等经典翼型，几何拟合误差控制在10−510−5弦长量级，同时可直接施加前缘半径、后缘厚度、最大厚度位置等工程约束，从根源上避免优化过程中生成翼面交叉、局部凹陷等非物理畸形构型。
2.3 二阶涡旋面板法求解框架
本研究采用的二阶涡旋面板法，将翼型表面离散为150段连续的二次曲线面元，在每个面元上布置线性分布的涡旋奇点，通过物面不可穿透边界条件与库塔后缘条件构建线性方程组，求解得到各面元的涡强分布后，叠加均匀来流速度即可得到全流场任意点的速度矢量。相较于传统常数面板法，二阶面板法的面元几何与奇点分布均采用二阶精度描述，在相同面元数量下，翼面速度分布的计算精度可提升约40%，有效降低了离散误差对临界马赫数计算结果的影响。
流场求解完成后，通过伯努利方程直接导出翼面不可压缩压力系数，再代入卡门-钱学森公式完成可压缩性修正，最终得到亚声速可压缩流条件下的全翼面压力分布，为后续临界马赫数的精准计算提供可靠数据支撑。
