小展弦比飞翼标模雷诺数影响数值模拟研究

采用数值模拟方法开展小展弦比飞翼标模的雷诺数影响研究。使用自行研制的多块对接结构网格亚跨超声速流场解算器程序Mbflow,计算了试验雷诺数、二倍试验雷诺数和飞行雷诺数等三种雷诺数情况下小展弦比飞翼标模的流场。通过对计算结果的分析研究,得到了不同马赫数(Ma=0.2、0.8和1.5)和攻角情况下雷诺数对小展弦比飞翼标模的气动特性和流场特征的影响规律。小攻角情况下,雷诺数主要影响摩阻的大小;而大攻角时,雷诺数对压阻也有明显影响。在亚声速时,雷诺数主要影响分离涡的起始位置和强度;跨声速和超声速时,雷诺数还会影响到激波的位置和强度。进一步研究了小展弦比飞翼标模的自准区雷诺数问题,发现试验雷诺数接近于小展弦比飞翼标模的自准区雷诺数。     

高超声速风洞多天平测力试验技术研究

在同一次试验中通过安装多台天平测量多个部件气动力,可以更加详细地了解飞行器的气动特性,并可以降低试验成本、提高试验效率。多天平测力技术在亚跨超声速风洞中应用已经非常成熟。但是,在高超声速风洞中,由于模型尺寸小、来流温度高,多天平的布局比较困难,所以目前国内在高超声速风洞中进行的多天平测力试验相对较少。为满足型号研制的需要,在 CARDC 超高速所的Φ1m 高超声速风洞上开展了多天平测力试验技术研究。选取某飞行器的3个控制舵为研究对象。针对试验模型径向尺寸小,而轴向尺寸相对较大的特点,确定采用铰链力矩天平轴线与测量舵转轴互相垂直的“纵轴式”布局方式;采用合理的小型化六分量天平结构形式,升降舵天平测量元件的长、宽、高尺寸为20mm×20mm×25mm;为更加真实地模拟天平的工作状态,设计了专用的加载头,并采用单元加载和组合加载相结合的方法进行校准,校准结果表明,天平主要分量的静校准度均在0．7％以内;在采取温度补偿措施之后,零点漂移明显减小;采用了天平-舵偏角变换装置一体化设计,有效提高了定位和安装精度。在Φ1m 高超声速风洞上开展的 Ma＝5条件下的风洞试验结果表明,采用的多天平测力试验技术是可行的,可以应用到高超声速风洞测力试验中。     

翼型尾缘低频大功率合成射流的动态载荷特性研究——B：动态气动载荷特性分析

针对NACA0015翼型,设计了适用于风洞研究的尾缘低频大功率合成射流致动器,对翼型尾缘合成射流作用下的非定常气动特性进行了风洞实验研究。研究表明：尾缘合成射流与横流的相互作用能够有效改变作用在翼型上的气动载荷;单侧喷口喷/吸时,喷冲程气动力系数响应幅值约为吸冲程幅值的3倍;双侧喷口同时工作时,升力和力矩系数的幅值并不是单侧喷口单独工作时喷气幅值的简单叠加,而是处于单侧喷气幅值和喷吸幅值和之间;升力和力矩响应的幅值与喷流动量系数的平方根之间存在近似的线性关系;在动量系数不变时,升力和力矩系数响应的幅值会随减缩频率的增加而减小;给定合成射流器行程,升力和力矩响应幅值与合成射流频率之间近似呈线性关系。     

有翼高超声速再入飞行器气动设计难点问题

有翼高超声速再入飞行器是近年来的研究热点,气动设计是飞行器设计的关键。为了更清楚地认识有翼高超声速再入飞行器气动设计的难点问题,对有翼高超声速再入飞行器的发展、优势及总体任务剖面进行了介绍,从5个方面详细介绍了该类飞行器气动设计的难点问题,包括多约束复杂面对称气动布局设计、高温真实气体效应对气动特性影响、天地差异与天地换算方法、反作用控制系统(RCS)喷流干扰对气动特性的影响以及气动数据不确定度等,简要阐明了这些难点问题对总体设计的重要性以及初步的解决思路,为有翼高超声速再入飞行器气动设计提供了一些参考。     

高超声速复杂气动问题数值方法研究进展

高超声速流场具有复杂流动特征,其中真实气体效应、磁流体干扰效应和力热结构耦合效应等对气动力分析产生了重要影响。将流体力学研究扩展到分子动力学、电磁流体力学以及流固耦合等交叉学科领域,这给数值模拟方法带来了巨大挑战。针对高超声速气动力/热分析的热点问题,重点关注高温效应与低密度流动效应、磁流体干扰效应和力热结构耦合效应等,结合算例分析了相应的数值求解技术;在气动热方面主要比较了3类求解方法(纯工程方法、纯数值方法和基于Prandtl边界层理论的方法),并给出了相应算例;对于气动力/热/结构耦合问题,从耦合模型及耦合计算方法两方面开展了分析。最后指出了高超声速复杂气动问题数值求解技术未来需重点关注的几个方面。     

高超声速飞行器表面温度分布与气动热耦合数值研究

针对高超声速飞行器热防护设计中的高温气体非平衡效应问题和气动热环境精确预测问题,基于流场的非平衡Navier-Stokes方程、表面的能量守恒方程和内部的热传导方程,考虑流场的非平衡效应、表面的热辐射效应、催化效应和烧蚀效应以及热防护层内部的热传导效应,建立了初步的表面温度分布与气动热的耦合计算方法,完善了高超声速飞行器气动物理流场计算软件(AEROPH_Flow)。在表面材料为碳-碳(C-C)条件下,对飞行高度为65km和飞行速度为8,10km/s的半球以及飞行高度为50km和飞行速度为8km/s的球锥模型,开展了表面温度分布与气动热的耦合计算,验证了计算方法和计算软件,分析了表面温度分布对气动热环境的影响。研究结果表明:表面温度分布对气动热的计算结果有较大影响,在气动热环境的预测中,不仅要考虑热化学非平衡效应和表面催化效应的影响,还要考虑表面温度分布的影响,最好是采用表面温度分布与气动热耦合计算的方法,以减小表面温度分布对气动热计算结果的影响。为此,需要发展完善非平衡流场/表面催化和烧蚀/热传导温度场(气/表/固)的计算模型、耦合求解技术和计算软件,实现对高超声速飞行器的真实飞行条件下高温气体非平衡效应和气动热环境的精确模拟。     

小型无人倾转旋翼机气动与操纵特性试验研究

由于倾转旋翼机飞行模式多,各部件气动干扰复杂且操纵面冗余,特别是倾转过渡模式,短舱带动旋翼系统倾转,结构布局发生改变,从理论上确定气动与操纵特性难度大。为了研究倾转旋翼机的气动与操纵特性,对某小型无人倾转旋翼机展开全尺寸、全模式吹风试验,其中不带动力试验主要用于研究倾转旋翼机在不同迎角、短舱倾角、前飞速度等飞行状态下的气动特性;带动力试验主要用于研究倾转旋翼机不同飞行模式带机翼与不带机翼时,旋翼/机翼/襟副翼相互干扰作用,以及总距、副翼、升降舵的操纵功效。根据试验数据推导出小型无人倾转旋翼机全包线飞行的操纵特性方法,对进一步完善倾转旋翼机设计以及试飞试验的成功提供了参考。     

柔性扑翼气动结构耦合特性数值研究

微型扑翼体积小、重量轻,其柔性变形对气动特性有显著的影响。通过求解雷诺平均N-S方程(ReynoldsAveraged Navier-Stokes,RANS)和结构动力学方程,对微型柔性扑翼飞行器的气动结构耦合特性进行了数值模拟研究。针对微型扑翼的大幅运动,发展了适用于扑翼的气动结构耦合数值计算方法,研究了微型扑翼的气动结构耦合特性。通过求解雷诺平均Navier-Stokes(RANS)方程得到微型扑翼的非定常气动特性;利用哈密顿原理(Hamilton Principle)推导了扑翼的结构动力学方程,采用结构有限元方法对该动力学方程进行离散并求解,得到扑翼的动态结构特性;采用松耦合方法进行迭代。计算结果与风洞实验结果相比吻合良好,验证了所发展方法的有效性。在此基础上研究了惯性力和关键运动参数对柔性扑翼气动及结构特性的影响规律,有助于比较详细、全面地了解微型扑翼的气动机理,为柔性扑翼的设计提供了参考依据。     

旋翼非定常气动特性CFD模拟的通用运动嵌套网格方法

针对直升机旋翼非定常气动特性CFD模拟中的网格生成难题,提出了一套高效、通用的运动嵌套网格生成方法.首先,基于Poisson方程求解和翻折法生成旋翼桨叶的正交贴体网格.其次,针对旋翼桨叶的扭转分布及变距、挥舞等复杂运动,建立了一套通用的洞单元识别的扰动衍射法;为保证洞包络面的封闭性,完善了挖洞过程中网格加密策略;在洞边界确立基础上,提出了一种高效、鲁棒的最小距离法贡献单元搜索的改进方法.在此基础上,建立了基于RANS(Reynolds-averaged Navier-Stokes)方程的旋翼非定常流场CFD模拟方法.最后,采用所建立的方法分别对悬停和前飞状态下的C-T(Caradonna-Tung)旋翼和7A(Helishape 7A)旋翼的气动特性、桨尖涡的位置进行了计算,计算结果与试验值误差小于5%,验证了该运动嵌套网格生成方法在旋翼非定常气动特性CFD模拟中的有效性.      

基于EFFD方法的自然层流短舱优化设计

采用extended free-form deformation(EFFD)方法研究了自然层流(natural laminar flow,NLF)短舱的气动外形优化设计方法.使用基于Bernstein基函数的EFFD方法完成了NLF短舱剖面的参数化,利用基于k-εSST(shear stress transport)两方程湍流模型的γ-θ转捩模型进行自然转捩预测,结合EFFD、一种混合动网格方法、Kriging代理模型和改进的粒子群算法(particle swarm optimization,PSO)建立了针对NLF短舱气动外形的优化设计框架.采用该框架分别对通气NLF短舱和带动力NLF短舱进行优化设计.单独通气NLF短舱优化结果的外表面实现48%的层流,阻力系数比初始通气NLF短舱减小了0.0003.带动力NLF短舱的优化结果外表面保持了41%的层流.这些结果表明采用相关技术建立的优化设计框架在NLF短舱设计中具有一定应用价值.