一种可用于机翼优化设计中的气动噪声预测模型

发展了一种可用于翼型/机翼外形设计中的气动噪声快速预测方法。相较于传统的半经验噪声预测方法,该方法以两方程非线性k-ε湍流模型模化雷诺应力的雷诺平均方程为背景,考虑了升力系数、三维流动效应以及机翼几何参数等因素对后缘噪声的影响。而相对于直接数值模拟或声类比拟方法,该方法虽不能准确预测噪声强度,但其计算量小,能给出不同翼型/机翼的相对总声压级,以及总声压级随升力系数的变化情况,易于应用于翼型/机翼气动外形优化设计中。通过计算分析二维NACA0012翼型几何参数或来流状态的改变所带来的气动噪声差异,与ANOPP软件及Brooks等计算结果进行对比,验证了该模型的可靠性。最后,计算分析二维、三维翼型/机翼气动噪声,凸显该方法在翼型/机翼气动外形优化设计中的应用价值。     

跨音速翼型上的激波／边界层干扰自适应控制计算

在激波区使用自适应壁对跨音速翼型的激波／边界层的相互作用（干扰）进行控制，可改变机翼的气动性能，这种被动控制可通过在翼型的激波区开一凹腔，其上覆盖一弹性橡胶膜柔壁来，本文给出用Ｎａｖｉｅｒ－Ｓｔｏｋｅｒ方程数值模拟这一自适应控制翼型的跨音速粘性绕流，提出了一个适应于本特殊情况（物面边界局部地区在求解过程中有变化）的处理办法。并探讨了自适应柔壁对当代跨音速翼绕流的影响。     

变可信度模型管理优化设计方法在气动外形设计中的应用

在信赖域最优化方法的基础上,将梯度搜索方法与近似模型管理结构相结合,结合N-S方程与Euler方程,建立了一种用于气动外形设计的梯度AMF优化设计方法.该方法在气动优化中,通过对采用不同的流体控制方程的计算结果进行组织和管理,从而可以在优化迭代过程中将主要搜索过程集中在低可信度模型的计算上,并利用高可信度模型监控优化过程,使得最终的优化解收敛到高可信度模型上.本文采用梯度AMF气动优化设计方法进行翼型设计,取得了良好的设计效果,实际结果表明本文提出的方法具有可行性和适用性.     

用拟压缩性方法求解不可压Euler方程

用拟压缩性方法和Ｊａｍｅｓｏｎ的有限体积算法求解了二维和三维定常可可压Ｅｕｌｅｒ方程。分别采用显、隐式时间离散推进求解;分析了人工粘性的阶数对定常解收敛性的影响,应用该方法计算了单个翼型和翼身组合体的低速绕流,结果与实验吻合较好。     

无网格算法在多段翼型流动计算中的应用

研究了一种求解欧拉方程的无网格算法，发展出了一套布点及点云自动生成的方法；在点云离散的基础上，采用最小二乘法求解矛盾方程的方法来求取空间导数，进而获得数值通量；采用四步龙格-库塔方法进行时间推进，并引入当地时间步长和残值光顺等加速收敛措施。通过对NA-CA0012翼型的跨音速流动和多段翼型复杂绕流的数值模拟，验证了上述无网格算法的正确性和实用性。     

褶皱幅度对仿生翼型滑翔气动特性的影响

蜻蜓翅膀具有独特的褶皱状形貌．研究者们致力于利用仿生学原理,设计在低雷诺数条件下具有更优气动性能的褶皱翼型．本文采用计算流体力学方法,求解二维不可压Navier-Stokes方程组,探讨了四种翼型（平板翼型、流线翼型、小幅度褶皱翼型和大幅度褶皱翼型）的气动表现．在低雷诺数条件下得到以下结果：(1) 较小幅度的褶皱结构有利于增加升力和减小阻力．(2) 雷诺数变化时褶皱翼型的升力系数呈非线性变化;在特定雷诺数区间,幅度相近的褶皱翼型会发生相对气动优势的转变．(3) 褶皱结构内的回流区通过减小粘性阻力,使得翼型总阻力下降．(4) 翼型前缘的极小区域会产生脉冲高升力,对升力表现产生较大影响．这些结果表明,调整褶皱幅度是实现褶皱翼型气动优化的有效方案．     

基于近似技术的涡轮叶片气动优化设计

根据五次多项式方法进行三维涡轮叶片的参数化建模,采用N-S方程和湍流模型进行三维流场分析计算,以K-S函数法作为优化方法,利用近似技术加速循环优化速度,建立了一种基于近似技术的涡轮叶片的气动优化方法。将气动效率和总压比作为目标函数,对涡轮叶片进行多目标气动优化、形状优化。算例表明本文提出的涡轮叶片优化设计方法是有效的。     

欧拉微分式的一种简化方法

欧拉微分式的一种简化方法冯贤桂（重庆大学，重庆６３００４４）在材料力学和弹性力学中平面极坐标问题的求解及轴对称问题的计算，常常遇到欧拉微分方程的求解．例如，轴对称厚壁圆筒的位移法求解的基本方程为￣［１］圆形薄板轴对称弯曲位移法基本方程￣［２］方程（１...     

排翼布局飞行器气动性能的实验研究

通过低速低湍流度风洞实验,研究了利用排翼布局改善充气飞机采用大厚度翼型机翼带来的气动效率偏低问题。首先比较了采用不同厚度翼型的单翼与排式双翼布局的气动特性。在此基础上,为了优化排翼布局的气动特性,研究了给后翼安装偏转角对排翼布局气动特性的影响。同时,基于NACA0030翼型,设计了波纹型外形的充气机翼,比较了此外形下单翼和排翼布局气动性能的差异。实验结果表明,采用排翼布局能够改善采用厚翼型单翼布局的气动性能,而给后翼安装一定偏转角可以进一步提高排翼布局的升力和升阻比。采用波纹外形和光滑外形机翼模型的对比结果表明,波纹外形能够在大迎角时改善充气机翼的失速性能。分析认为,造成这一现象的流动机理是由于波纹型机翼在实验条件下提前由层流转捩为湍流,使失速推迟,流动分离现象有所减弱。     

翼型、机翼非定常运动模式下动态数值模拟

基于RANS方程,通过刚性动网格技术实现对翼型和机翼典型运动模式的描述,采用双时间推进方法和Roe空间离散格式对流场求解,构建了一个非定常气动计算平台;以NACA0012翼型为算倒进行了动态数值模拟可信度验证。数值模拟结果与试验数据吻合较好,升力和俯仰力矩的最大计算误差分别为3％和10％,表明了该平台的可靠性。另外,还数值模拟了M6机翼的动态非定常流场,并分析了两种湍流模型对非定常流场激波的捕捉能力。结果表明非定常流动中S-A湍流模型对激波的捕捉较B-L模型更敏感。文中开发的非定常计算平台对进一步解决三维复杂流场的流动问题有很高的工程应用价值。     

AERODYNAMIC OPTIMIZATION OF HIGH-SPEED TRAIN BASED ON RBF MESH DEFORMATION

基于局部型函数三维参数化方法、改进的蚁群算法和改进的克里金（Kriging）代理模型,开展了列车头型的三维气动减阻优化设计研究。为了避免复杂几何外形大变形情况下千万量级网格的重复生成,提高高速列车头型优化设计的效率,引入了缩减控制点的径向基函数网格变形技术。优化结果表明：径向基函数网格变形技术在不降低网格质量的情况下可以有效缩短网格变形的时间消耗,能够用于复杂几何外形的气动优化设计;在给定的设计空间内,控制鼻锥外形的6个关键设计参数对列车气动阻力的影响呈单调递增关系;优化后,在满足约束条件的情况下,简化外形列车的整车气动阻力减小5.41%,头尾车减阻效果明显,中间车气动阻力基本不变。     

基于自适应多可信度多项式混沌-Kriging模型的高效气动优化方法

多可信度代理模型已经成为提高基于代理模型的优化算法效率和可信度水平最有效的手段之一.然而目前流行的co-Kriging和分层Kriging (HK)等多可信度代理模型泛化能力不足,缺乏对高阶/高非线性建模问题的适应性,难以广泛应用.文章基于发展的自适应多可信度多项式混沌-Kriging (MF-PCK)代理模型,在提高建模效率和对高阶/高非线性问题近似准确率的同时,建立了基于该自适应MF-PCK模型的高效全局气动优化方法.在发展的方法中,提出了基于MF-PCK模型的新型变可信度期望改进加点方法,使代理优化算法效率进一步提高.为了验证发展方法的全面表现,将其应用在经典的数值函数算例以及多个跨音速气动外形的确定性优化和稳健优化设计中,并与基于Kriging和HK模型的代理优化算法进行了全面比较.结果表明,发展的新型多可信度全局气动优化方法其优化效率相对于基于Kriging和HK模型的优化效率显著提高,结果更好也更加可靠,并且稳健优化设计效率和结果也更符合工程应用需求,证明了其相对于基于Kriging和HK模型的代理优化算法的显著优势.     

FAST OPTIMIZATION METHOD OF REENTRY TRAJECTORY CONSIDERING AERODYNAMIC PARAMETER PERTURBATION

针对传统再入轨迹优化方法收敛速度慢、对初值敏感程度高等的局限性,提出了一种基于序列凸优化的再入轨迹快速求解方法.该方法以倾侧角的变化率作为控制量,改进了现有凸化策略,考虑到抑制数值优化过程中由于数值离散方式带来的锯齿化现象,采用 B 样条曲线离散控制量,同时为避免算法在初始猜想值附近出现伪不可行的问题,增加额外虚拟控制量,通过一种"回溯直线"搜索的方法,提高算法的稳定性、快速性和寻优结果的光滑性.为研究飞行器再入过程中的气动参数扰动问题,采用采样点少、易于实现,计算效率高的广义混沌多项式理论研究方法,建立了基于广义混沌多项式和凸优化相结合的再入轨迹鲁棒优化模型,该模型在优化过程中考虑气动参数扰动对寻优结果的影响作用,避免了传统轨迹与制导律的复杂迭代设计环节,可有效降低优化轨迹对气动参数扰动的敏感程度,在气动参数不确定条件的干扰下,依然可以保证飞行器顺利安全的完成飞行任务.最后,以美国某可重复使用飞行器的再入任务为例,验证了基于序列凸优化的再入轨迹优化方法的快速性以及鲁棒优化模型对气动参数扰动的抗干扰性能力,表明了该方法具有一定的工程应用性.     

考虑气动参数扰动的再入轨迹快速优化方法

针对传统再入轨迹优化方法收敛速度慢、对初值敏感程度高等的局限性,提出了一种基于序列凸优化的再入轨迹快速求解方法.该方法以倾侧角的变化率作为控制量,改进了现有凸化策略,考虑到抑制数值优化过程中由于数值离散方式带来的锯齿化现象,采用 B 样条曲线离散控制量,同时为避免算法在初始猜想值附近出现伪不可行的问题,增加额外虚拟控制量,通过一种"回溯直线"搜索的方法,提高算法的稳定性、快速性和寻优结果的光滑性.为研究飞行器再入过程中的气动参数扰动问题,采用采样点少、易于实现,计算效率高的广义混沌多项式理论研究方法,建立了基于广义混沌多项式和凸优化相结合的再入轨迹鲁棒优化模型,该模型在优化过程中考虑气动参数扰动对寻优结果的影响作用,避免了传统轨迹与制导律的复杂迭代设计环节,可有效降低优化轨迹对气动参数扰动的敏感程度,在气动参数不确定条件的干扰下,依然可以保证飞行器顺利安全的完成飞行任务.最后,以美国某可重复使用飞行器的再入任务为例,验证了基于序列凸优化的再入轨迹优化方法的快速性以及鲁棒优化模型对气动参数扰动的抗干扰性能力,表明了该方法具有一定的工程应用性.      

融入压力分布信息的气动力建模方法

气动外形优化设计与飞行器性能分析中, 直接运用数值模拟或风洞实验获取气动力的成本高, 构建代理模型是提高外形优化和性能分析效率的重要途径. 然而, 构建模型的过程中, 研究者只关注积分后的气动力和力矩信息. 本文通过充分利用采样过程中所产生的压力分布信息, 来提高建模的精度和泛化性, 进而降低样本获取的成本. 提出了一种小样本框架下融入压力分布信息的气动力建模方法, 首先通过数值模拟或风洞试验获得不同流动参数状态下翼型表面的压力分布信息和气动系数, 其次通过本征正交分解技术对压力分布信息进行特征提取, 获取不同输入参数状态下压力分布信息对应的POD系数, 之后结合输入参数通过Kriging算法对压力分布信息进行建模, 将压力分布信息积分得到低精度气动系数的预测模型, 最后低精度气动系数结合输入参数通过Kriging算法构造高精度的气动系数预测模型. 通过同状态变翼型算例以及CAS350翼型变状态算例进行验证, 该方法相比于传统的克里金模型直接预测气动力, 有效提高了气动力的预测精度和模型的鲁棒性, 同时缩小了学习样本的数据量.      

基于PCE方法的翼型不确定性分析及稳健设计简

由于能够获得一个既经济又对参数变化不敏感的设计结果,稳健型设计在工程设计中备受关注. 不确定性分析是稳健型设计的关键. 因此研究了基于混沌多项式的不确定性分析方法,并将其与CFD 方法结合,对计算空气动力学设计中的不确定性影响进行了量化分析. 首先以RAE2822 翼型为算例,对其跨音速马赫数不确定影响进行了分析,研究了多项式阶次对计算的影响,分析了平均流场和方差. 接着结合超临界翼型的马赫数稳健型设计验证了混沌多项式方法在稳健型设计中的有效性. 优化结果表明,稳健型优化后的翼型阻力系数明显降低,同时对于马赫数的敏感性显著减小. 通过分析表明混沌多项式方法能够大幅提高稳健型优化设计效率,能很好地应用于气动稳定性设计.     

基于PCE方法的翼型不确定性分析及稳健设计

由于能够获得一个既经济又对参数变化不敏感的设计结果,稳健型设计在工程设计中备受关注. 不确定性分析是稳健型设计的关键. 因此研究了基于混沌多项式的不确定性分析方法,并将其与CFD 方法结合,对计算空气动力学设计中的不确定性影响进行了量化分析. 首先以RAE2822 翼型为算例,对其跨音速马赫数不确定影响进行了分析,研究了多项式阶次对计算的影响,分析了平均流场和方差. 接着结合超临界翼型的马赫数稳健型设计验证了混沌多项式方法在稳健型设计中的有效性. 优化结果表明,稳健型优化后的翼型阻力系数明显降低,同时对于马赫数的敏感性显著减小. 通过分析表明混沌多项式方法能够大幅提高稳健型优化设计效率,能很好地应用于气动稳定性设计.      

Unsteady aerodynamic reduced-order modeling of an aeroelastic wing using arbitrary mode shapes

Computational fluid dynamics (CFD) based unsteady aerodynamic reduced-order model (ROM) can offer significant improvements to the efficiency of transonic aeroelastic analysis. To construct a ROM based on mode shapes, one run of CFD solver is needed to compute aerodynamic responses corresponding to mode excitations. When mode shapes change with structure, another run of the CFD solver is required to construct the new ROM. The typically large computational cost associated with repeated runs of the CFD solver impedes the application of existing unsteady aerodynamic reduced-order modeling methods to transonic aeroelastic design optimization and aeroelastic uncertainty analysis. This paper demonstrates a method that can replace the CFD solver used in the process of existing unsteady aerodynamic reduced-order modeling. It can produce aerodynamic responses corresponding to mode excitations for arbitrary mode shapes within a few seconds. Computational cost can be reduced by two orders of magnitude using the mode excitations and the corresponding aerodynamic responses computed by the method to construct the ROMs used for flutter analyses in aeroelastic design optimization or aeroelastic uncertainty analysis in transonic regime compared with the existing unsteady aerodynamic reduced-order modeling methods. Results show that the method can accurately produce the aerodynamic responses corresponding to the mode excitations and predict the flutter characteristics of AGARD 445.6 wings root-attached in three different ways.     

基于响应面法的新型排翼式飞艇的气动优化设计

采用均匀试验设计的方法, 以后掠排翼式飞艇为初始气动外形, 后翼的轴、法向坐标为设计变量, 升力、阻力系数和升阻比为目标函数, 通过CFD结合响应面近似理论, 建立了气动参数的二阶响应面优化模型, 并针对排翼式飞艇的升阻特性进行了优化设计. 研究结果表明: 响应面方法是解决气动优化问题的一种行之有效的方法;优化后排翼式飞艇的气动性能得到明显改善, NACA0020, NACA0030, NACA0040排翼式飞艇的升阻比分别提高了9.61%, 6.08%, 13.08%; 飞翼的轴向和法向位置对各气动参数的影响不尽相同, 升力系数C_rm l和升阻比k对轴向位置更为敏感, 而轴向和法向位置对阻力系数C_d的影响则基本相当.