30P-30N多段翼型复杂流场数值模拟技术研究

采用"亚跨超CFD软件平台"(TRIP2.0)数值模拟了30P-30N多段翼型的复杂流场,主要目的是考核湍流模型、转捩位置对多段翼型压力分布和典型站位速度型的影响。本文通过求解任意坐标系下的雷诺平均的N-S方程,采用多块对接结构网格技术,在与相应试验结果对比的基础上,详细研究了SA一方程湍流模型、SST两方程湍流模型、不同的转捩位置对该翼型压力分布和典型站位速度型的影响。本文的研究结果表明,采用全湍流模拟方式可以较好地模拟该多段翼型的压力分布,但对速度型的模拟精度较差;模拟试验的转捩位置可以改善主翼附面层与前缘缝翼边界层尾迹区的模拟精度;采用微吸气技术推迟前缘缝翼的转捩位置,可以进一步提高缝翼尾迹区的数值模拟精度。     

考虑转捩影响的风力机翼型气动特性计算研究

基于二维雷诺平均Navier-Stokes(RANS)方程,耦合eN数据库方法进行流动转捩判断,采用SA湍流模型计算了风力机专用翼型DU91-W2-250的气动特性,得到了雷诺数Re=1.0×106时,攻角在-11°~12°范围内该翼型的升力和阻力特性曲线及压力分布。通过对全湍流计算结果、耦合转捩判断的计算结果及实验数据的对比,表明考虑转捩因素后的计算结果与实验结果吻合更好,说明要准确模拟风力机翼型的绕流,必须考虑转捩影响。     

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

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

高超声速飞行器机体/推进一体化设计的启示

机体/推进一体化设计是吸气式高超声速飞行器的关键技术。飞行器的前体和后体既是主要的气动型面,又是发动机进气道的外压缩型面和尾喷管的膨胀型面,一体化设计直接影响飞行器的气动与发动机性能。本文阐述了吸气式高超声速飞行器的主要特点,梳理了飞行器的推阻匹配、升阻比特性、操稳匹配等主要气动设计问题。通过对国外典型高超声速飞行器机体/推进一体化设计技术的综合分析,总结了前体/进气道、后体/尾喷管、边界层强制转捩装置等关键部件的气动设计方法,获得了有意义的启示,可为后续吸气式高超声速技术研究提供重要参考。     

吸气式高超声速飞行器机体推进一体化技术研究进展

吸气式高超声速一体化飞行器最显著的特点是子系统之间的耦合较其他类型飞行器更加强烈,这使得其设计具有挑战性。所有的子系统之间部件相互干涉,包括:气动、推进、控制、结构、装载和热防护等,特别是机体与超燃冲压发动机之间的耦合最为突出。飞行器的前体和后体下壁面既是主要的气动型面,又是超燃冲压发动机进气道外压缩型面和尾喷管的膨胀型面,在产生推力的同时也产生升力和俯仰力矩。机体与发动机的强耦合作用对飞行器的推力、升力、阻力、俯仰力矩、气动加热、机身冷却、稳定性和控制特性有直接的影响。本文介绍了国内外机体推进一体化技术的研究进展,重点介绍了中国空气动力研究与发展中心(CARDC)的相关研究工作,包括:密切曲锥曲面乘波进气道和基于双激波轴对称基准流场内转式进气道设计方法、独创的大尺度脉冲式燃烧加热风洞一体化飞行器带动力试验技术和高超声速内外流耦合数值模拟技术等。对高速飞行中激波边界层相互干扰、流动分离机理、可压缩湍流转捩及其控制、超燃冲压发动机燃烧流动机理等相关基础问题也进行了研究,强调了对高效高精度计算方法的迫切需求。     

超声速平板边界层旁路转捩直接数值模拟

采用具有七阶精度的有限差分格式和基于Euler方程组的特征边界条件,通过直接数值模拟研究了由固定相位周期性吹吸扰动诱导的、Ma∞=2.25、Re∞=2.5×107/m(635 000/inch)空间发展平板边界层的旁路转捩过程,并建立了相应的直接数值模拟数据库。数值结果表明:扰动相位的随机性并不是转捩发生的必要条件,固定相位周期性吹吸扰动即可诱发边界层的转捩;当采用固定相位时,随着吹吸频率的增加,转捩位置向下游移动;当扰动频率足够大时,周期性吹吸无法诱导出边界层的转捩。随后,在所得数据库中提取了速度及热力学变量的平均值,与已有的数据库和理论分析吻合良好;速度脉动均方根与公认的不可压缩平板边界层实验基本吻合;热力学变量脉动均方根和雷诺应力与已有数据库的统计结果及实验数据具有相同的分布特征。对比发现,不同转捩的诱导因素对平均转捩位置和热力学变量均方根、雷诺应力等高阶统计量的影响较大。     

基于壁面压力分布的非对称喷管反设计

针对目前非对称喷管的设计方法上的缺陷,给出了通过指定壁面压力分布规律来反设计其膨胀面型线的方法,获得了膨胀面型线反设计程序,并结合优化算法寻找综合性能较好的喷管壁面压力分布.将采用该方法设计得到的喷管模型与最大推力喷管进行了对比研究.结果表明:在设计点,该喷管的推力系数比最大推力喷管只降低0.102%,而升力和俯仰力矩分别提升2.295%和15.774%.验证了设计思想的正确性,为非对称喷管的设计提供了一种高效的设计方法.      

非对称结构损伤飞机鲁棒容损控制器设计

当飞机发生非对称结构损伤时,飞机的质量、重心位置和气动特性都会发生突变,飞机机体的对称性遭到破坏,致使飞机的横纵向间运动产生强烈的耦合。针对飞机发生非对称结构损伤时导致的飞行控制问题,建立了非对称结构损伤飞机的损伤模型,并基于一种新型鲁棒容损控制策略,采用非线性扩张状态观测器和非线性动态逆相结合的方法,对飞机的姿态控制器进行了设计,兼顾了飞机系统的性能和对损伤的鲁棒性。最后,基于NASA的通用运输机模型,对所设计控制器的控制效果进行了仿真验证。仿真结果表明,设计的姿态控制器有效地抑制了非对称结构损伤给飞机控制系统带来的不确定性和扰动,具有较好的控制性能。     

AZ31B镁合金大应变循环变形行为研究

为了研究AZ31B镁合金在大应变幅条件下的变形机制,开展了该合金在7.嬲应变幅条件下的循环行为研究.结果表明:在拉伸阶段的最大应力值随着循环周次的增加而减小,而在压缩过程中的最大应力值随着循环周次的增加而增大,在整个循环过程中材料呈现循环应变硬化特性,拉应力是导致循环应变硬化的主要原因;随着循环周次增加,滞回曲线的不对称性基本不变.真应力-真应变滞回曲线在卸载和反向拉伸阶段出现3个拐点.在压缩过程中发生{10-12}孪生,反向拉伸过程发生去孪生行为,包申格效应对去孪生行为具有较大影响.研究表明:孪生-去孪生是大应变幅循环变形的主要变形机制;对拉伸、反向压缩过程的变形特征及机制的分析,可为低周疲劳行为的研究提供参考.     

基于边界层相似性解的放大因子输运模型

为了使二维线性稳定性理论能够适应现代CFD求解技术,通过求解Falkner-Skan边界层相似性方程获得各个形状因子下的相似速度型。基于线性稳定性理论,对各个速度型进行稳定性分析从而获得对应不同速度型的扰动放大因子包络线。最后使用标量输运方程的形式实现包络近似方法中放大因子的当地化求解,并结合原始γ-Re_θt转捩模型中的间歇因子输运方程,实现了自然转捩和分离泡转捩的建模。使用该输运模型对S&K平板、S809翼型、NLR7301翼型和DLR-F5机翼进行转捩预测,结果均与试验结果吻合较好,验证了该模型构建的合理性和可行性。