基于γ-Reθt转捩模型的翼型气动特性计算精度改进研究

将结合SSTk-ω湍流模型的γ-Reθt转捩模型与二维雷诺平均NS（RANS）方程数值求解相耦合,发展了考虑转捩影响的翼型气动特性计算方法。得到的计算结果与没有考虑转捩影响的计算结果的比较表明,与实验更加吻合,计算精度得到了明显提高。     

NS方程中耦合转捩判断的风力机翼型气动力数值模拟

在Reynolds-Averaged Navier-Stokes(RANS)方程计算中耦合Michel经验判据和Chen-Thy-son转捩模型对风力机翼型S809进行计算。由RANS方程求得翼型表面压力分布作为层流边界层方程求解的输入参数,然后使用Michel经验判据分析判断层流边界层的解得到转捩点的位置,使用Chen-Thyson转捩模型得到从层流到湍流的转捩过渡区,这样随着流场的迭代求解,求解器自动更新判断出转捩点位置。对S809翼型进行全湍流和耦合转捩判断的RANS方程计算,可以看到考虑转捩判断后得到的升阻力系数与实验值吻合较好,验证了方法的可行性。     

基于放大因子与Spalart-Allmaras湍流模型的转捩预测

为了验证放大因子输运方程与Spalart-Allmaras(S-A)湍流模型耦合对转捩现象的模拟精度,选取Schubauer and Klebanoff(S-K)平板、S809低速翼型、30p30n多段翼型以及复杂的三维HiLiftPW-1构型进行自由转捩计算,并将计算结果与实验进行比较分析,其中针对S809算例,还与Langtry-Menter(L-M)转捩模型进行了比较.算例结果表明:放大因子输运方程与S-A湍流模型的耦合能够较好的捕捉转捩位置以及转捩发展过程,对分离泡诱导的转捩的模拟相比L-M转捩模型更精确,转捩位置的捕捉精度提升了10%;对比实验,多段翼转捩位置的捕捉误差最大为6.5%;针对三维高升力增升构型,以实验作为参考,全湍流计算与考虑边界层转捩的对比显示考虑边界层转捩能够更加精确的模拟气动力系数,升力和表面摩擦阻力系数的模拟精度精度提升1%.      

转捩对风力机翼型和叶片失速特性影响的数值模拟

采用基于k-ωSST湍流模型的Gamma-Theta转捩模型对S809翼型和NREL Phase VI叶片进行了考虑转捩的气动力数值模拟,研究了转捩对其气动特性特别是失速特性的影响。首先对S809翼型在迎角0°～30°范围内开展了数值模拟,比较了转捩模拟和全湍流模拟获得的翼型升力系数和流场特征,发现前缘层流分离泡的存在显著影响了翼型的失速特性;然后对NREL Phase VI叶片开展了类似的数值模拟,结果显示转捩对叶片失速特性和翼型失速特性产生影响的作用方式是相似的。经过对以上数值模拟结果加以分析认为,转捩对翼型和叶片失速特性的影响主要通过前缘层流分离泡的作用体现出来,前缘层流分离泡的存在使翼型和叶片更早地进入深失速。     

NS方程计算中耦合转捩自动判断的阻力精确计算方法初探

在Reynolds-Averaged Navier-Stokes(RANS)方程计算中耦合了流动转捩的自动判断以提高现有求解器预测翼型阻力的准确性.由RANS方程求得翼型表面压力分布作为层流边界层方程求解的输入参数,然后使用简化的eN-数据库转捩判断方法分析层流边界层的解得到转捩点的位置,这样随着流场的迭代求解求解器自动判断转捩点的位置.在对NLF0416翼型的气动性能计算中考虑流动转捩的因素后得到的翼型升阻力特性和实验吻合较好,验证了本文方法的正确性.     

翼型几何模型对边界层转捩预测精度的影响

对飞行器表面的边界层转捩现象进行准确地预测,是设计高性能层流翼型的关键。基于边界层数值模拟和现代CFD技术的特点,选取与SSTk-ω湍流模型相耦合的γ-Reθt转捩模型,对不同类型的翼型边界层转捩进行预测。结果表明:在转捩计算过程中出现的摩阻曲线震荡,是由翼型表面计算网格不光顺所引起的;提高翼型几何模型光顺程度能够在不同程度上提高翼型气动特性的计算精度,并能得到更加准确的转捩预测结果。     

基于遗传算法及转捩模型的层流翼型优化设计研究

讨论了层流翼型对飞机机翼减阻的重要意义,概述了遗传算法的优化原理及其特点。在基准翼型的基础上,以Hicks—Henne型函数的解析函数线性叠加法来描述翼型。采用遗传算法完成了层流翼型优化设计。计算中采用了转捩模型耦合SST两方程湍流模式,通过求解雷诺平均的N-S方程组模拟了翼型的转捩流动。对计算结果进行了分析,分析表明优化翼型在升力特性、阻力特性方面比基准翼型有很大的提升,达到了优化升阻比的目标。     

使用GAO-YONG湍流方程组对翼型绕流的计算

将GAO-YONG湍流方程组应用到翼型复杂绕流场的计算中,研究了机械能方程中Du/Dt项和各向异性湍流黏性系数分量中的实度系数Cs的影响,给出了控制Du/Dt项作用区域和选择实度系数Cs值的判别函数,对NACA4412翼型最大升力迎角下的低速绕流算例、AS240翼型8°迎角下的转捩算例、AS240翼型19°迎角下的分离流算例进行了模拟,与实验结果比较,得到了比较好的结果.结果表明GAO-YONG湍流方程组能够准确预测翼型前缘层流到湍流的转捩问题,翼型尾迹流动松弛问题和翼型大迎角分离流.     

转捩模型在翼型气动性能预测 中的应用研究

边界层转捩位置的准确预测对于提高飞行器气动性能的预测精度具有重要意义.选取与k-ω SST湍流模型相耦合的γ-Reθt模型,以零压力梯度平板为研究对象,通过求解基于有限体积法的雷诺平均N-S方程验证该模型自动捕捉流动转捩的准确性;将该模型应用于传统有压力梯度的NACA 0012翼型的流场特性和气动性能的研究中,并与原始k-ω SST模型的计算结果及全湍流试验数据进行比较.结果表明:远场边界距离对翼型阻力系数有较大的影响;与无转捩模型相比,γ-Reθt转捩模型对翼型阻力系数的预测精度有一定程度的提高;对于二维模型,γ-Reθt转捩模型难以捕捉翼型表面的三维效应和非定常分离特性.     

涡轮叶片气热耦合数值模拟方法

研究了气热耦合计算模型,建立了采用有限差分方法进行涡轮气热耦合计算的数值仿真平台,并利用该平台采用考虑转捩和不考虑转捩两类湍流模型对具有径向冷却通道的MARK II型叶片进行数值模拟,同时将结果与实验值进行比较,对比表明:考虑转捩的湍流模型能够准确的反映出边界层内的流态和热传导情况,因此壁面处温度和换热系数更接近于实验值,与实验值吻合较好。     

Love波飞机结冰传感器研究

Love波结冰传感器有望应用于飞机结冰预测与监测,其通过Love波在导波层中的传播变化,来敏感导波层表面冰的厚度以及状态变化。在Zimmermann理论基础上,根据Love波结冰传感器实际结构,针对导波层为冰层及压电材料为ST90°X石英的实际实验情况进行了仿真计算,计算得到了冰层厚度与相速之间的特性关系,以及冰层厚度对灵敏度的影响,并与N.Bari等的灵敏度理论的计算结果进行了比较,结果表明Love波结冰传感器的灵敏度与相速的理论计算是正确的。通过初步的实验对Love波结冰传感器进行原理性的验证研究,实验结果表明Love波结冰传感器可以明显区分空气、水和冰的状态,实验数据符合理论计算的结果。     

某型涡扇发动机部件老化对性能影响的分析与计算

详细分析了部件老化对发动机性能影响的机理.利用所能得到的统计数据, 基于部件匹配技术, 通过将考虑部件老化影响的部件特性嵌入到发动机稳态性能计算模型中, 建立了考虑部件老化的发动机性能计算模型, 采用Newton-Raphson迭代收敛技术求解非线性方程组, 获得考虑部件老化影响的发动机性能.以某型涡扇发动机为例, 定量地计算分析了部件老化对发动机性能的影响, 计算结果和发动机实际试车数据对比具有较好的一致性, 表明本文发展的方法是合理的.      

喷管流场和羽流场的N－S数值模拟

利用时间推进的方法，通过求解可压缩的、雷诺平均的Ｎ－Ｓ方程，并利用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ两层代数湍流模型和Ｐｒａｎｄｔｌ混合长度模型，数值模拟了飞机和火箭发动机喷管的内流和羽流，获得了喷管内流场和羽流场的参数分布，计算结果与试验结果及国外研究成果一致。网格划分采用了自适应网格技术。     

喷管流场和羽流场的N－S数值模拟

利用时间推进的方法，通过求解可压缩的、雷诺平均的Ｎ－Ｓ方程，并利用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ两层代数滴流模型和Ｐｒａｎｄｔｌ混合长度模型，数值模拟了飞机和火箭发动机喷管的内流和羽流，获得了喷管内流场和羽流场的参数分布，计算结果与试验结果及国外研究成果一致。网格划分采和了自适应网格技术。     

2m×2m超声速风洞引射器气动设计

采用一维气体引射器流动特性方程计算了该等面积混合引射器在不同运行状态下的性能参数,设计增压比为2.1~2.7,引射系数为0.392~0.651,综合结构设计要求确定引射喷嘴数目为24个,径向布置两层,面积比为0.235.引射器1∶13缩比模型的试验结果与计算值具有良好的一致性,5个设计状态下的实测增压比均略高于设计值.2m×2m超声速风洞调试时因整体参数调整导致引射器在某些工作状态偏离了设计点,但增压比与引射系数的对应关系仍与设计结果基本一致,引射器总体性能指标满足风洞试验需求.     

高速弹射试验机操稳性能的计算与试飞

本文介绍了为研究高速弹射试验机的操稳性能和在弹射力扰动下飞机的动态响应,把弹射力以瞬态脉冲载荷的形式作用于飞机时所建立的数学模型和采用的计算方法.把计算结果与试飞结果和原型机试飞结果作了比较,计算结果与试飞结果的一致性是令人满意的。这为高速弹射试验机的安全飞行和指导飞行员实施弹射试验提供了可靠的依据。     

粘性流体对流换热问题用直角坐标网格的解法

研究了非结构化直角坐标网格全贴体求解粘性不可压流体自然对流换热问题。该方法在于利用直角坐标网格但通过在边界附近保留不规则控制体,使得算法是完全贴体的。利用此数值方法对二维倾斜空腔内的自然对流及上固体边界移动时的对流换热问题进行计算,并与STAR-CD的计算结果进行比较。对比结果表明,该方法计算的结果与STAR-CD的结果很接近,可以满足一定精度要求。      

某型弹用涡喷发动机安装性能分析

在分析S形进气道的总压损失以及进气道，喷管／后体外部阻力计算方法的基础上，建立了针对某型号导弹涡喷发动机的安装性能的计算程序，并结合实际飞行条件，对用户附加引气，功率输出，进气道，喷管安装和环境温度对该发动机性能的影响进行了计算，结果表明，引气和提取功率影响不大，而进气道总压损失和进、排气系统的外部阻力有重要影响，环境温度有一定影响；而且由于弹用涡喷发动机的推力比较小，必须考虑安装性能计算，以便准确判断每一种安装因素所造成的性能损失，为设计和使用方提供理论指导依据。     

H—6飞机在低空风切变中的动态响应特性和下滑航迹的研究

本文用修正的哈明预测校验计算方法,计算了H-6飞机在水平风的垂直切变和三向风切变的情况下的动态响应特性和下滑航迹的变化,并根据需要给出了改善飞机对风切变适应能力的自动驾驶仪结构,用快速数字仿真方法计算出了改善结果.     

X-cor夹层结构压缩强度模型改进与实验验证

详细阐述了X-cor夹层结构固化工艺过程中残余应力的产生过程,分析结果为碳纤维Z-pin受残余拉应力.考虑工艺过程残余应力的影响,将Z-pin视为弹性基础上梁,提出一种X-cor夹层结构压缩强度改进计算模型,并计算其压缩强度.制作了不同Z-pin密度、直径和角度的X-cor夹层结构实验件用于压缩强度实验.实验值和模型计算值的对比验证了该压缩强度计算模型的正确性.工艺过程中Z-pin所受残余拉应力随Z-pin直径、植入密度增加而增加,随Z-pin植入角度增加而减小.