面向非结构混合网格高精度阻力预测的梯度求解方法

针对非结构混合网格的特点,通过改进传统的Green-Gauss梯度求解方法,提出了一种可提高非结构混合网格黏性计算精度的节点型Green-Gauss梯度求解方法。利用改进后的方法,完成了DLR-F4翼身组合体算例的计算和对比分析。改进后的梯度求解方法残差收敛更好,下降量级更多,阻力系数和试验吻合更好,激波区域压力分布和分离区域流场细节的模拟更精确,说明改进后的梯度求解方法有效提高了程序的鲁棒性和阻力预测精度,验证了方法的有效性。采用改进后的方法对第5届AIAA阻力预测研讨会的通用研究模型（CRM）进行了详细的模拟分析,结果表明：改进的梯度求解方法更加适用于非结构混合网格的黏性计算,计算精准度达到国际同类CFD软件水平,进一步验证了改进方法的可靠性。     

高速气流冲击柱状泡沫多孔体的传热特性分析

采用流体/多孔区域一体化单区域算法,数值研究了高速绕流条件下前置于圆柱体前缘表面的柱状泡沫多孔体内部的传热特性。基于蒙特卡罗法考虑多孔域内的辐射热效应,分析了变化多孔区域长度和多孔阻力特性对模型激波阻力和前缘多孔区域气动热的影响。结果表明:在圆柱体前缘安置一定长度及带有适当阻力特性的泡沫多孔材料,可同时减小整体激波阻力并降低前缘表面的气动热效应。在模拟工况下,无量纲长度1.0、黏性阻力系数0.2×107m-2及惯性阻力系数200m-1的前缘泡沫多孔可减小激波阻力13.5%,降低约75%的前缘表面的平均气动热流密度。保持无量纲长度不变,减小泡沫多孔区域惯性阻力系数会降低激波阻力,但会略微增加前缘壁面气动热流密度。      

旋转整流罩积冰过程中水膜流动和脱离的数值研究

为了更准确地模拟旋转整流罩积冰过程,在现有三维积冰与冰层表面薄水膜流动耦合模型基础上,基于功平衡分析的方法引入了旋转部件表面水膜脱离模型,并发展了相应的计算方法,给出了水膜脱离的判定依据:当气流曳力做的功和由于离心力使水膜增加的潜能之和大于黏附功时整流罩表面的水膜会发生脱离。对旋转整流罩积冰进行数值模拟,计算结果与实验结果吻合得较好,验证了该模型的合理性和计算方法的可行性。之后分析了转速和来流速度对整流罩表面水膜脱离和积冰的影响,结果表明:转速和来流速度越大,水膜发生脱离的比例越大。在研究范围内,转速为3000r/min和6000r/min时,因水膜脱离导致积冰总量分别减少13.4%和15.8%;来流速度为40、50m/s和60m/s时,因水膜脱离导致积冰总量分别减少为12.2%、13.4%、14.2%。      

仿生减阻翼型的气动性能

利用一种基于滑移边界理论的仿生肋条结构的模化方法,对仿生减阻翼型的气动性能开展了数值模拟研究。结果显示:肋条结构同时布置在翼型压力面和吸力面湍流区域时,仿生减阻翼型的阻力可降低1.73%~3.07%,升阻比可提高2.10%~4.08%,其中黏性阻力的减小是总阻力下降的主要原因。对流场分析表明,仿生减阻翼型的速度曲线出现一定抬升,使得湍流边界层中的黏性子层增厚,进而导致了仿生减阻翼型气动性能的提升。      

基于大涡模拟的平流层浮空器气动特性分析

为研究不同球体间距对“球形囊体型”平流层浮空器(SLV)气动特性的影响,采用有限体格式、结构网格和大涡模拟(LES)计算方法求解不可压缩的Navier-Stokes(N-S)方程,对超临界雷诺数“球形囊体型”平流层浮空器绕流进行数值模拟,并对不同球体间距下的数值计算结果进行详细的分析比较.通过对比试验数据,单球体数值模拟的阻力系数时均值与Achenbach的试验数据一致,验证了计算方法分析超临界雷诺数球体绕流问题的准确性.研究不同间距的双球体阻力变化规律以及振动频谱特性,随上下游球体间距G的增加,合阻力先增大后减小,上游球体的阻力占优振动频率逐渐减小;G=1.5D(D为球体直径)和G=2D时,上下游球体有相同的占优振动频率.随间距G的增加,两球体相互作用与上游球体对下游球体的尾涡结构影响逐渐减弱.      

双后掠鸭翼气动特性的数值模拟

首先针对具有中等前缘后掠角梯形鸭翼的缺点提出双后掠鸭翼概念,然后分别对安装梯形鸭翼和双后掠鸭翼的近距耦合鸭式布局的气动性能进行数值模拟研究,分析影响双后掠鸭翼气动性能的流动机理。研究表明：在大迎角时,对于双后掠鸭翼,具有较大前缘后掠角的外翼段可以使鸭翼涡在涡核破裂后仍能形成稳定集中涡并保持较高的强度,增加鸭翼本身的失速迎角,并通过诱导作用改善机翼外翼段流场,进而提高全机大迎角性能,但在小迎角时会破坏鸭翼附着流或前缘气泡涡的发展,造成略微的升力损失。拥有较大失速迎角的双后掠鸭翼在小迎角时具有较大的可用偏度,可以增强布局的抬头控制能力。双后掠鸭翼在满足隐身约束的前提下,超声速阻力较小,具有较好的超声速性能。     

亚声速翼身融合无人机概念外形参数优化

为了兼顾翼身融合（BWB）布局无人机（UAV）的气动、隐身和结构重量要求,应用优化方法研究了某亚声速翼身融合无人机概念方案的外形设计问题。外形优化设计流程包括全机参数化几何外形模型、气动分析、机翼根部弯矩计算、雷达散射截面（RCS）分析、代理模型的建立和外形参数优化计算。选择了三种不同优化目标研究翼身融合无人机外形优化问题：①未配平状态升阻比最大;②配平状态升阻比最大;③配平状态升阻比尽量大和机翼弯矩尽量小。通过优化结果的对比分析,揭示了配平约束和机翼弯矩目标对优化设计结果的影响。研究结果表明：计入配平约束能够有效提高配平升阻比;将配平状态升阻比尽量大和机翼根部弯矩尽量小作为优化目标能够获得合理的优化外形。     

旋成体表面沟槽减阻试验研究

在旋成体模型表面粘贴铝基沟槽蒙皮,研究沟槽表面对模型阻力的影响.给出基于来流单位雷诺数的沟槽尺寸无量纲公式,用于设计沟槽齿高和齿宽,依据试验风洞的单位雷诺数范围,选定3种尺寸的沟槽进行减阻试验.结果表明,在旋成体表面湍流流动区域顺流向布置沟槽,当沟槽齿高和齿宽的无量纲值h+、s+均小于25时,在小迎角下具有减阻作用,当s+约为15时减阻效果最明显,可减小约3%~4%的阻力.     

流体横掠水滴形微针肋热沉流动和传热特性

对传统圆形微针肋进行了优化,设计了3种不同尾角的水滴形微针肋热沉,并以去离子水为工质,实验研究了各热沉流动阻力和传热特性.结果表明:3种尾角针肋中,尾角为60°时减阻效果最好.水滴形针肋的流线型结构可以改善尾部流动分布,推迟流动由层流向过度区流的转变,且尾角越小效果越明显.不同体积流量下,水滴形针肋的最优尾角有所不同.在实验中,雷诺数范围在200～1000内,尾角为60°的水滴形针肋热沉强化换热效果最好.当尾角为30°时,太长的尾部结构受到下一排针肋的影响,造成较大的流动阻力,导致其整体换热效果较差.      

含有SMA弹簧驱动器的可变倾斜角翼梢小翼研究

 针对传统翼梢小翼在非设计状态减阻效果不佳的缺点,提出一种含有形状记忆合金(SMA)弹簧驱动器的变体翼梢小翼结构,它能根据飞行状态主动调整小翼的倾斜角,实时优化飞机的阻力特性.采用力-热-应变耦合法设计了所需的SMA弹簧驱动器,并通过有限元仿真与风洞试验验证了变体翼梢小翼的变形能力,最后初步研究了变体翼梢小翼的闭环控制方法.研究结果表明,在飞机的起飞阶段(自由来流流速为26 m/s,迎角为3°),变体翼梢小翼的倾斜角能在1 min内自主完成预定变化过程,倾斜角的最大变化量为23°,控制精度的最大误差为12%,各项指标均符合设计要求.