翼型模型几何误差对气动性能影响的自动微分分析方法

基于自动微分原理和NS方程有限体积方法建立了一套翼型敏感性导数计算方法和程序,可以一次性获得翼型不同气动力系数、压力分布对模型几何外形误差的敏感性导数和不确定度。计算结果表明,在跨声速范围内,即使翼型的外形误差只有63μm(弦长1m),也可以给翼型压力分布带来0.312(以来流动压为参考)量级的不确定度,而激波处的流动最为敏感。这种自动微分方法对于分析数值模拟结果的分散性、风洞试验结果的分散性或不确定性具有很好的指导意义。     

基于多项式混沌法的翼型不确定性分析及梯度优化设计

耦合伴随方法和非嵌入式多项式混沌法，发展了高效、可靠的不确定性梯度优化设计方法。利用伴随方程法求解目标函数对不确定性变量的导数，发展了一种梯度增强型多项式混沌法。通过亚声速和跨声速下等多种算例可以证明该方法可以提高不确定性分析的效率和精度。同时，利用基于方差分解的全局敏感性分析方法对不确定性变量的敏感性进行了量化。建立了多项式混沌耦合伴随方程的统计矩梯度求解方法，并结合梯度增强型多项式混沌法搭建不确定性梯度优化设计系统。基于该优化设计系统对二维低亚声速和跨声速翼型开展确定性及不确定性优化设计研究。优化结果显示，相比于确定性优化设计，不确定性优化设计通过合理权衡确定性性能和不确定性性能，可提高抵抗马赫数和迎角不确定性扰动的能力，同时优化性能均值和标准差。其中阻力系数均值最大可降低17%，阻力系数标准差最大可降低80%。而确定性优化设计可能导致性能鲁棒性的降低。     

自然层流翼型高雷诺数风洞试验研究

针对所设计的自然层流翼型,在高速风洞中开展了自然层流翼型高雷诺数风洞试验研究。试验运用了测压和红外成像等技术手段,在马赫数Ma=0.74和雷诺数Re=1.5×10~7条件下获得了压力分布、升力系数、阻力系数、俯仰力矩系数和红外图像。试验结果表明,在试验条件下,该翼型上表面能够保持良好的顺压梯度;该翼型具有良好的纵向气动特性;在迎角α=0°时,上、下翼面流动均能保持约60%的层流面积区。     

翼型纵向俯仰非定常气动力研究

飞行器飞行时受到的非定常气动力是飞行器设计以及机动飞行中须考虑的一个重要因素。采用CFD方法数值模拟了NACA0012翼型在马赫数从0.2~5.0范围变化时做强迫振荡的非定常流场。通过积分法计算了在不同马赫数时翼型的俯仰动导数。计算结果表明,迟滞效应随马赫数的增大先增大后减小,在跨声速区域迟滞效应达到最大;迟滞环的绕向在跨声速区域为顺时针方向;迟滞环顺时针绕向时,动导数大于零,模型动不稳定,逆时针绕向时,动导数小于零,模型动稳定;旋转角速度越大,翼型表面非定常气动效应越明显。     

特殊大气环境下格尼襟翼的气动特性

临近空间飞行器和火星探测飞行器首要面临低雷诺数、跨声速的特殊气动问题,针对此特殊环境下低雷诺数翼型的格尼襟翼增升方案,基于CFD方法开展了数值计算研究。对不同雷诺数、马赫数条件下,后缘加装不同高度格尼襟翼的Eppler387翼型的气动特性进行了对比分析,结果表明格尼襟翼增大了翼型的前缘吸力和后缘压差,从而显著增大了翼型环量和升力。在较高的马赫数下,格尼襟翼使翼型上表面低压区的范围扩大,对激波位置略有推迟。合适高度的格尼襟翼提高了翼型的最大升阻比,在中高升力系数下能够明显增大翼型的升阻比。研究结果能够为临近空间飞行器和火星探测飞行器的设计及改进提供技术支持。     

牛顿流理论及其在谐振机翼上的应用

根据牛顿流理论,文中讨论了高马赫数、任意迎角、一般剖面在均匀来流中作沉浮和俯仰谐振的翼型和机翼。导得翼型和机翼表面的牛顿撞击压强和安定性导数公式。给出楔、平板三角翼、双抛物线剖面三角翼俯仰刚性导数(-C_(mθ)和阻尼导数(-C_(m))的数值例子。数值结果与现有的牛顿-布兹曼流理论、薄激波层理论、欧拉方程摄动理论的结果作了比较,算例表明本文的结果令人满意。     

现代自然层流翼型的设计方法

本文介绍一种设计跨声速自然层流翼型的计算流体力学(CFD)方法。本方法采用“正反迭代、余量修正”设计原理,通过将跨声速翼型设计软件NPU-TD2D中的反设计程序进行改进,并与含有层、湍流混合边界层修正的跨声速层流翼型计算程序DLRBGKWALZ耦合,实现了在跨声速粘性流动条件下直接设计层流翼型。亚、超临界的设计实例和风洞验证表明,本方法可以在几个设计迭代内设计出压力分布、转捩位置及气动参数均准确收敛于设计目标的新翼型,是一种设计现代自然层流翼型的有效而实用的CFD方法。     

类X-37B航天器气动力天地相关性数值模拟

随着高性能计算机的发展,CFD已成为飞行器设计和流场分析不可缺少的重要手段,风洞试验与飞行数据的天地相关性问题正是其中一项重要的研究内容,X-37B作为继航天飞机之后美国发展的最成功的可跨大气层在轨飞行器,从气动特性角度分析其大气层内飞行走廊的状态对中国类似航天器的研制具有重要的借鉴意义。首先,对计算类X-37B布局飞行器的网格无关性及网格修正开展了研究,在此基础上提出的网格规模影响修正方法对该类飞行器的计算结果修正经过验证是可信的;然后,分别对比分析了雷诺数的影响和试验状态支架干扰的影响,完成了基于数值模拟的高空飞行与风洞试验气动特性差异分析。结果表明,网格规模主要对亚声速来流计算状态压差产生的轴向力影响较大,对法向力系数、俯仰力矩系数和纵向压心影响较小;雷诺数对该类飞行器气动特性特别是轴向力系数、阻力系数和升阻比有较大的影响,但随着马赫数的增加,影响特性开始变的非常复杂;由于风洞试验状态支杆存在,亚跨声速来流条件对该类飞行器的底阻影响很大,需要采取一定的方法和手段对支杆影响进行修正。     

超临界自然层流机翼设计及基于TSP技术的边界层转捩风洞试验

为了实现绿色航空节能减排的目标,层流设计技术成为飞行器设计者的研究热点。对于跨声速客机而言,超临界自然层流机翼设计技术将显著减小飞行阻力,提升气动性能,减少燃油消耗和污染物排放。首先,基于高精度边界层转捩预测技术耦合翼型优化设计系统,实现超临界自然层流翼型设计;经过合理的翼型配置,形成超临界自然层流机翼。转捩数值模拟分析结果表明,超临界自然层流机翼的层流流动特性良好。然后,以比例为1：10.4的试验模型在荷兰高速低湍流度风洞进行边界层转捩风洞试验,使用温度敏感材料涂层（TSP）技术拍照获得机翼表面在不同马赫数、雷诺数和迎角工况下的层流-湍流分布。最后,通过超临界自然层流机翼边界层转捩试验结果,探讨了该类型机翼的转捩特性随来流参数的变化规律,总结了超临界自然层流机翼设计的关键因素。此外,该模型也可用来验证边界层转捩预测技术在超临界、高雷诺数工况下的预测精度。     

三维副翼铰链力矩计算

 数值模拟了带副翼偏转构型的二维、三维亚、跨、超声速来流情况的绕流流场,采用有限体积方法和AUSM+迎风格式数值求解N-S方程,四步龙格-库塔时间推进,湍流模型为Baldwin-Lomax和Spalart-Allmaras湍流模型。二维模拟压强分布与实验数据以及CFD软件模拟结果吻合较好。采用交错对接网格,数值模拟了不同副翼偏角和缝隙、来流马赫数、迎角和雷诺数的三维构型流场。通过舵面铰链力矩与实验数据的对比分析表明了该数值方法模拟此类复杂流场的可靠性和舵面铰链力矩计算的有效性。发现并研究了亚声速来流时铰链力矩随迎角的反向变化趋势,初步分析了副翼的缝隙和雷诺数效应。     

基于非定常CFD方法跨声速民用飞机标模俯仰气动导数计算研究

采用非定常计算流体力学(CFD)方法和动网格技术,计算了跨声速民用飞机全机构型TCR气动标准模型在0°~30°迎角范围内气动力系数的时均值、俯仰静导数和俯仰动导数;通过与试验结果对比,研究了非定常CFD方法在较大迎角范围内对包含鸭翼、垂尾等部件的复杂外形气动导数的计算精度。与试验结果对比表明,在16°~24°迎角范围内,受鸭翼和机翼上流动分离再附以及内外翼分离涡相互作用等复杂流动结构的影响,俯仰力矩系数的动导数与试验数据有相对较大的差异。     

2.4m跨声速风洞试验质量影响因素分析及改进措施研究

影响风洞试验质量的因素很多,如流场品质、测量系统误差、支撑干扰以及洞壁干扰等。主要对模型姿态角、马赫数、模型支撑系统等影响因素进行了改进研究。通过改进使模型迎角测量精度达到0.03°、Ma数控制精度达到0.003,并有效降低了支撑干扰影响,提高了2.4m跨声速风洞的试验质量。     

高雷诺数下跨声速自然层流短舱优化设计

为解决高雷诺数下大涵道比发动机自然层流短舱高维优化设计问题,提取短舱3个基准面实现轴对称自然层流优化设计。通过获得较大范围层流区域,从而降低短舱表面摩擦阻力。采用类别形状函数（CST）参数化、γ-Re_θt转捩模型和遗传算法建立自然层流（NLF）短舱自动优化流程。表明通过优化短舱基准面建立三维NLF短舱设计方法的可行性。进而通过CATIA二次开发构建三维非轴对称NLF短舱,解决了基准面优化后大量数据点的高效导入和曲面生成问题。针对设计的非轴对称NLF短舱,进行了设计点附近迎角、侧滑角、来流马赫数以及湍流度的转捩敏感性分析。结果表明:跨声速状态下,迎角增大层流范围减小;来流马赫数增大层流范围扩大;侧滑角和湍流度对层流范围影响很小。      

翼吊发动机机翼跨声速颤振特性研究

针对现代民用飞机翼吊发动机机翼跨声速颤振问题,建立了带超临界翼型、大展弦比后掠和翼吊一个发动机构型的机翼模型,通过偶极子格网(DLM)气动力修正方法、升力线斜率系数修正方法、ZTAIC方法和高速颤振模型风洞试验方法研究了飞机翼吊发动机跨声速颤振特性。分析了马赫数和阻尼对跨声速颤振特性的影响。结果表明:翼吊发动机机翼具有三种典型颤振,即机翼弯曲扭转耦合模态颤振、机翼面内振动模态颤振和发动机与机翼耦合模态颤振;马赫数对翼吊发动机机翼跨声速颤振会产生不利影响,主要表现在跨声速区域显著降低机翼跨声速临界颤振动压;结构阻尼对机翼-发动机耦合小阻尼颤振模态的影响显著,增大阻尼可以显著提高其颤振动压。     

跨声速高雷诺数圆球绕流的实验研究

本文采用无支撑干扰、洞壁影响很小的带加压抽空系统的弹道靶设备,进行了雷诺数在10~5相似文献   

典型超声速/高超声速动导数计算方法研究

以当地流活塞理论为基础完成超声速/高超声速动导数计算公式推导,发展了一种高效的动导数计算方法.选取1个超声速、2个高超声速动导数算例标模,计算了超声速/高超声速下飞行器动导数变化规律.研究结果表明:BFM超声速流动算例的俯仰组合动导数绝对值随马赫数的增大而减小,计算重心位置后移亦会导致俯仰组合动导数绝对值减小,纵向动稳定性降低;在0° ～20°迎角范围内,高超声速流动中的尖锥和钝头旋成体俯仰组合动导数绝对值随迎角的增大而增大,纵向动稳定性增强;所提方法计算时间约为双时间动导数计算方法的1/27.     

小迎角高超声速非旋转钝锥非线性运动分析

采用分段拟合技术对小迎角高超声速(马赫数Ma=6)非旋转钝锥双平面拍摄风洞自由飞试验结果进行了分析,获得了非旋转钝锥在小迎角高超声速下气动导数非线性的具体形式,分析了非旋转钝锥在小迎角下的非线性运动特点。研究发现在小迎角范围内钝锥的动导数系数均呈现明显的非线性,而静导数系数非线性较弱,可近似为线性。各组试验的静、动导数系数的非线性形式表明,除模型Ⅱ外的其余模型均为两方向振幅不同的极限圆锥运动,模型Ⅱ为极限平面运动。每组试验两个平面的静导数系数在小迎角范围内保持基本相等。而无论出现极限圆锥运动的四组试验还是出现极限平面运动的一组试验,模型在俯仰和偏航两个方向的动导数系数均存在较大差异。     

来流参数对防热瓦横缝旋涡结构及热环境的影响

针对高超声速飞行器表面缝隙内部流动,通过求解可压缩Navier-Stokes方程,自主研发了一套能够较好模拟缝隙流动特性的计算流体力学(CFD)软件。利用该软件研究了来流参数对防热瓦横缝旋涡结构及热环境的影响。计算结果表明:随着来流雷诺数的增加,缝内旋涡结构呈现主涡个数增多形态趋于饱满的变化趋势,缝隙壁面绝对热流和无量纲热流增加;随着来流马赫数的增加,缝内主涡个数、形态基本不变,但主涡旋转速度增加,缝隙壁面绝对热流增加,无量纲热流基本不变;随着来流迎角的增加(迎角较小时),缝内旋涡结构和热流变化规律基本与增加来流雷诺数相同。由此分析可知,涡量向下传递并形成旋涡的距离,即形成所谓"死水区"的深度,主要由来流雷诺数和来流迎角决定。     

跨声速翼型多目标优化设计方法

本文通过采用多目标分级优化方法并基于求解全速势方程的跨声速粘流翼型计算方法，研究发展了一种多参数、多约束和多目标的跨声速翼型数值优化设计方法，应用该方法可以从普通低速翼型和超临界民办型出发通过多目标优化后得到在跨声速区的多个马赫数下阻力系数最小化的翼型几何外形，设计实践表明，该方法具有收敛快，调用目标函数次数少等优点。     

用保角转绘坐标解二维和轴对称可压流Euler方程

本文推导出用保角转绘正交曲线坐标表示的二维与轴对称Euler方程,方程中的未知量以曲线坐标方向的速度分量表示,并有不含未知量导数的源项。方程组的形式与直角坐标系中的接近,因而简化了计算。本文将该方程组用于计算二维任意翼型和旋成体轴对称绕流,采用分辨率较高的Eberle特征平均显式有限体积差分算法求解,在不同区域中用不同的显式推进步数加速收敛。本方法已分别对翼型和旋成体轴对称绕流编制了计算程序,可用于自由流为亚声速、跨声速和超声速的绕流计算。大量计算结果表明,本方法精度好,激波分辨率高,编程序简单。