第一届航空CFD可信度研讨会总结

总结了第一届航空可信度研讨会(AeCW-1)的数值模拟结果。AeCW-1会议由中国空气动力学会计算空气动力学专委会主办,并在第一届空气动力学大会上做了专题研讨。本次研讨会吸引了14家单位参与,共提交了26组计算结果。活动的主要目的是针对典型运输机巡航构型,评估国内CFD当前技术状态,探索CFD进一步发展的研究方向,促进CFD验证和确认工作的开展,为大型客机的研制提供技术支撑。AeCW-1选择的研究模型为中国空气动力研究与发展中心(CARDC)设计的单通道运输机模型(CHN-T1),CHN-T1包括机身/机翼/平尾/立尾等部件,分别在CRADC的2.4m×2.4m风洞(FL-26)、欧洲的2.4m×2.0m风洞中(European Transonic Wind Tunnel,ETW)开展了风洞试验。AeCW-1组委会提供了系列结构网格和非结构网格,研究工况主要包括网格收敛性研究、抖振特性研究、雷诺数影响研究,风洞模型支撑装置和静气动弹性变形对数值模拟结果的影响是重要的研究内容。本文对提交组委会的数值模拟结果进行统计分析,并与FL-26及部分ETW的风洞试验结果进行了对比,对未来的研究工作提出了意见和建议。     

单通道客机气动标模CHN-T1设计

发展单通道客机标模,特别是针对国内客机研制的标模,一方面利于确认CFD软件计算的可信度,从而促进CFD软件能力提升,为客机气动外形设计提供可靠的分析工具;另一方面有利于确认风洞试验品质,校对测量仪器,改良干扰修正方法,研发先进测量技术。本文介绍了用于确认风洞试验和CFD可信度的标模CHN-T1的气动外形设计。该标模包含机身、机翼、平尾、立尾、短舱、吊挂等部件。对机翼开展了详细优化设计,并配置单通道的窄体机身、满足稳定控制需求的平尾和立尾、翼下吊挂的通气模型短舱,组成全机干净构型,具有现代单通道客机的典型几何外形特征。标模设计马赫数为0.78,设计升力系数为0.5。机翼采用高气动效率的超临界翼型,在有/无短舱/吊挂组件的影响下,机翼均展示了良好的气动性能。相关信息可为即将召开的"第一届航空CFD可信度研讨会(AeCW-1)"提供基础。     

高阶精度格式WCNS在三角翼大攻角模拟中的应用研究

采用5阶精度的加权紧致非线性格式(WCNS-E-5)数值模拟了65°三角翼的大攻角绕流流场,主要目的是考核高阶精度格式WCNS在大攻角旋涡流动方面以及跨声速流场的激波附面层干扰、涡破裂位置的模拟能力,重点研究不同网格规模和湍流模型对尖前缘三角翼涡系之间的相互作用的影响。通过求解任意坐标系下的雷诺平均N-S方程,采用5阶精度的加权紧致非线性格式(WCNS-E-5)和多块对接结构网格技术,两种湍流模型分别是一方程SA和两方程SST湍流模型,在与相应试验结果对比的基础上,详细研究了WCNS-E-5格式在跨声速大攻角旋涡流动中的表现,以及不同网格规模、两种湍流模型对主涡二次涡相互作用、涡破裂位置和表面压力分布的影响。本文的研究结果表明,高阶精度格式WCNS-E-5能成功应用于三角翼的跨声速大攻角流动,网格规模的增加进一步提高流场分辨率,SST湍流模型相对SA湍流模型在三角翼大攻角流动中具有更好的适用性。     

不同襟翼偏角梯形翼构型气动特性数值模拟

基于雷诺平均Navier-Stokes(RANS)方程和结构网格技术,采用二阶离散精度的单调迎风格式(MUSCL),结合剪切应力输运(SST)两方程湍流模型和γ-Reθ转捩模型,研究了梯形翼高升力构型襟翼偏角变化对气动特性的影响。主要目的是进一步确认Trisonic Platform(TRIP)软件模拟高升力梯形翼不同襟翼偏角引起微小气动特性变化的能力。首先,简要介绍了采用的计算方法;其次,介绍了两种襟翼偏角的梯形翼模型及风洞试验;最后,在网格收敛性研究的基础上,采用全湍流和转捩两种方式模拟了梯形翼构型不同襟翼偏角对气动特性的影响。与试验数据的对比结果表明,采用全湍流和转捩两种方式均可以较好地模拟不同襟翼偏角对气动特性的影响量,采用γ-Reθ转捩模拟方式可以提高梯形翼构型气动特性的模拟精度,失速迎角附近的气动特性模拟需要进一步研究。     

高阶精度方法下的湍流生成项对低速流动数值模拟的影响研究

基于雷诺平均的Navier-Stokes方程和结构网格技术,采用五阶空间离散精度的WCNS格式,开展了SST两方程模型不同湍流生成项对低速流动数值模拟的计算分析。主要目的是为高阶精度格式在复杂外形上的应用提供技术支撑。计算模型包含了低速NLR7301两段翼型和Trap Wing高升力构型,研究内容主要包括不同湍流生成项对收敛历程、边界层湍流粘性系数分布、边界层速度分布、压力系数分布、气动特性的影响。在与试验数据对比的基础上,计算结果表明:对于低速二维流动,不同湍流生成项对收敛历程有比较明显的影响,对附面层湍流粘性系数分布和速度型影响不明显,不同湍流生成项主要影响主翼前缘的吸力峰值,进而影响升力系数和压差阻力系统;对于低速三维流动,不同湍流生成项对低速流动的收敛特性影响不明显,对翼梢涡的模拟精度有比较明显的影响,进而影响翼梢站位的压力分布和总体气动特性。     

高性能计算支持的静气弹优化设计方法

为扩大国产超级计算机在国内航空航天工业领域的应用,面向大展弦比机翼气动结构综合优化设计的实际问题,基于具有大量工业应用和自主知识产权的CFD软件TRIP和CSD软件SiPESC,应用改进的三维类别/形状变换（CST）参数化方法、差分敏度并行计算策略和移动渐近线（MMA）等大规模优化算法,通过TRIP+SiPESC耦合气弹分析软件预测气动载荷,发展了一套机翼型架外形的优化设计方法。以CRM标模内插机翼为例,针对无弯扭的初始型架构型,以气动外形展向几何扭转分布和结构有限元模型一万个壳单元厚度为设计变量,实现1 g载荷气动优化和2.5 g载荷结构优化的综合;在不降升阻比、满足最大应力/位移的约束下,俯仰力矩系数改善近40%,针对本文假设材料和初定的尺寸参数,机翼减重近9%,验证了方法和软件的有效性。该方法直接采用产品级高保真的自主静气弹分析软件预测气动力/载荷,将计算资源耗费大的静气弹优化和设计变量/约束规模大的结构优化分治,有较好工程应用前景。     

CRM-WB风洞模型高阶精度数值模拟

通过求解雷诺平均Navier-Stokes （RANS）方程,采用5阶空间离散精度的WCNS和多块对接结构网格技术,开展了CRM翼身组合体风洞试验模型的高阶精度数值模拟,计算构型和计算状态来自AIAA第6届阻力预测研讨会（DPW Ⅵ）。主要目的是采用高阶精度方法,评估静气动弹性变形和模型支撑装置对CRM翼身组合体数值模拟结果的影响。通过与刚性外形的计算结果和NASA NTF风洞试验结果的对比,高阶精度数值模拟结果表明：迎角为3.0°时,静气动弹性变形使得机翼上表面激波位置前移并显著降低了外侧机翼上表面激波位置前的负压;迎角为3.0°时,模型支撑装置使得机翼上表面激波位置进一步前移,并导致翼梢处机翼上表面流动结构发生变化;迎角为4.0°时,计算模型中没有包含模型支撑装置是导致升力系数下降的主要原因;计算模型中包含机翼静气动弹性变形和模型支撑装置的数值模拟结果更加接近试验结果。     

无尾布局支撑干扰数值模拟

采用TRIP2.0_SOLVER软件,开展了高、低速来流条件下,无尾飞翼布局支撑干扰的数值模拟技术研究.高速来流状态下支撑干扰的数值模拟,以尾支撑作为主支撑,腹支撑作为辅助支撑,完全模拟了实验中采用的两步法;低速来流状态下的支撑干扰的数值模拟,采用了有无尾支撑的方式得到支撑干扰量.笔者介绍了支撑干扰的典型高、低速数值模拟结果,并与高速实验结果做了初步的对比.     

高阶精度方法下的湍流生成项对跨声速流动数值模拟的影响研究

利用五阶空间离散精度的WCNS格式和多块结构网格技术,通过求解雷诺平均NS方程,开展了SST两方程模型不同湍流生成项组合方式对跨声速流动数值模拟影响的计算分析。研究的主要目的是为高阶精度格式在复杂外形上的工程应用提供技术支撑。计算模型采用了RAE2822超临界翼型和DLR-F6翼身组合体构型。研究内容主要包括不同湍流生成项对残差收敛历程、边界层湍流粘性系数分布、边界层速度分布、压力系数分布以及模型整体气动力特性的影响。不同湍流生成项组合方式的流场计算结果还与风洞试验数据进行了对比。研究结果表明:对于小迎角不存在明显分离的跨声速流动,不同湍流生成项对流场的高精度计算结果的影响很小,可以不用考虑。     

低压储能的升浮一体飞行器总体参数研究

升浮一体飞行器是一种综合了太阳能无人机和平流层飞艇特点的新型临近空间飞行器,是航空航天领域研究的热点,可再生燃料电池是升浮一体飞行器最有发展前景的储能装置。针对现有燃料电池中氢气和氧气存储方式的缺陷,充分利用升浮一体飞行器巨大的机身体积,提出了采用低压气囊储气的燃料电池方案,建立了升浮一体飞行器总体参数计算模型,提出常规燃料电池和气囊储气燃料电池的理论计算方法。研究了不同储气方式对升浮一体飞行器的能源系统以及飞行器总体参数的影响。研究表明,采用低压储气方法后,氢气的质量储气密度可提高至13.0%,燃料电池存储能量为1 489.7 kWh条件下,能量密度可达到1 000 Wh/kg以上,对减小升浮一体飞行器总重和外形尺寸、提高飞行器载荷能力有显著作用。