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

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

CRM-WBPN风洞试验模型数值模拟

基于雷诺平均Navier-Stokes（RANS）方程和多块对接结构网格技术,采用亚跨超CFD软件平台（TRIP 3.0）,研究了共同研究模型（CRM）风洞模型支撑装置对CRM翼/身/架/舱（CRM-WBPN）模型压力分布及气动特性的影响。构造了一个包含小、粗、中、细网格的4套网格,并进行了网格收敛性研究,在升力系数为0.50条件下,开展了固定升力系数下模型支撑对CRM-WBPN模型和CRM-WB模型压力分布及挂架短舱阻力增量的影响,以及不同来流迎角下模型支撑对CRM-WBPN模型气动特性的影响。通过与NASA国家跨声速风洞设备（NTF）试验结果和CRM-WB模型数值模拟结果的对比分析,表明模型支撑对固定升力系数下挂架短舱阻力增量影响较小;在计算构型中考虑模型支撑的影响,引起升力系数和阻力系数下降、俯仰力矩系数增加,数值模拟结果更加接近相应的试验结果。     

CRM翼身组合体模型高阶精度数值模拟

基于五阶空间离散精度的WCNS格式,开展了CRM翼身组合体模型的高阶精度数值模拟,以评估WCNS格式对复杂外形的模拟能力以及典型运输机巡航构型阻力预测的精度。首先依照DPW组委会提出的网格生成指导原则,利用ICEM软件生成了粗、中、细、极细四套网格,网格规模从"粗网格"的2 578 687个网格点逐渐扩展到"极细网格"的65 464 511个网格点。研究了设计升力系数下,网格规模对气动特性、压力分布和翼根后缘局部分离区的影响,采用"中等网格"开展了抖振特性的数值模拟研究。通过与二阶精度的计算结果、DPW V统计结果和部分试验结果的对比分析,高阶精度数值模拟结果表明,阻力系数计算结果与DPW V统计平均结果吻合较好;网格密度对机翼上表面的激波位置和翼身结合部后缘局部分离区略有影响;迎角为4°时,升力系数下降的主要原因是机翼上表面激波诱导分离区和翼身结合部后缘局部分离区的增加。     

DPWⅢ机翼和翼身组合体构型数值模拟

采用亚跨超CFD软件平台(TRIP)数值模拟了DPW Ⅲ提供的DPW_W1/W2两种机翼构型和DPW-F6/F6_FX2B两种翼身组合体构型,主要目的是通过两种机翼构型和两种翼身组合体构型的数值模拟,研究网格密度对运输机构型气动特性计算结果的影响。数值模拟采用的多块对接网格来自AIAA CFD Drag Prediction Workshop Ⅲ(DPW Ⅲ),采用National Transonic Facility(NTF)的试验结果和CFL3D的计算结果作对比。详细研究了网格密度对两种机翼构型和翼身组合体的总体气动特性和压力分布的影响。采用SST两方程模型计算两种构型均得到了网格收敛结果,网格密度主要影响压差阻力而对摩擦阻力影响较小,计算结果较好地预测了机翼和翼身组合体外形优化前后总体气动特性的变化量。     

考虑动力影响的民用飞机静气动弹性分析方法

采用了一种基于多块网格的N-S方程和结构柔度影响系数法,考虑气动、结构非线性的基于RBF插值和RBFDelaunay动网格变形技术的静气动弹性分析方法对喷流对弹性机翼的气动力影响进行了研究。利用DLR F6翼身组合体构型对静气动弹性方法进行验证,保证了计算的可信性。采用该方法对比分析了某民用飞机无喷流/有喷流构型的静气动弹性特性,表明发动机喷流会给机翼带来一个正的扭转效应,抵消一部分机翼后掠效应的影响,使机翼前后缘挠度均会有所增大,弹性变形引起的多数剖面的附加扭转角有所减小。研究表明:喷流影响会使刚性机翼表面的压力分布发生变化,升力系数有所损失;考虑喷流的机翼静气动弹性变形是一个耦合效应,发动机喷流区主要受喷流影响,外翼段主要受弹性变形影响。数值模拟结果表明:无喷流影响时机翼的弹性变形使升力系数下降约16%,升阻比下降8.4%,考虑喷流影响时,升力系数下降达到18%,升阻比下降36%。因此,对于大展弦比机翼,考虑喷流影响的静气动弹性分析十分必要。     

采用TRIP3.0模拟CHN-T1模型气动特性

采用TRIP3.0软件平台(TRIsonic Platform version 3.0),对第一届航空CFD(Computational Fluid Dynamics)可信度研讨会组委会提供的运输机标模构型(CHiNa-Transport,CHN-T1)进行了流动数值模拟。本文采用粗、中、细三套网格及百亿极细网格进行了网格收敛性研究;采用相应网格分别从气动特性、压力系数分布曲线、表面流态三个方面分析,研究了支撑装置、机翼静弹性变形以及雷诺数效应对CHN-T1构型气动特性的影响。数值模拟结果与试验结果有良好的一致性。数值模拟得到了网格收敛结果;支撑装置对力矩特性影响较大;机翼静弹性变形对气动特性影响较小;雷诺数效应对气动力特性影响较大。     

运输机标模CHN-T1的验证确认计算研究

采用自研的非结构混合网格解算器CARIA-OVERSET对我国第一届航空CFD可信度研讨会提供的运输机标模CHN-T1进行了验证确认计算。重点研究了网格形式、网格尺度、湍流模型等对运输机标模CHN-T1的总体气动特性、局部流场特征、压力分布以及计算收敛性等的影响,并把计算结果与风洞试验数据进行比较分析。分析结果表明,非结构混合网格解算器的计算结果与风洞试验数据吻合度较好。在固定升力条件下,非结构网格的网格尺度误差大于结构网格,湍流模型对迎角和俯仰力矩的影响较大,翼身结合处的分离特性模拟差异较小。带支撑构型的全机升力和低头力矩增加、阻力减小,机翼静气弹变形后升力、阻力和低头力矩均减小。全湍和自由转捩模型的计算结果存在一定差异,差异随着雷诺数的增大而变小。     

翼/短舱/吊架气动特性及其动力影响数值研究

应用计算流体动力学(CFD)数值模拟方法对DLR-F6模型分别在通流及动力条件下翼/短舱/吊架(W-N-P)组合体部分的气动干扰特性进行研究.结果表明:组合体相互干扰引起的干扰通道附增激波及短舱后部逆压导致的涡流增加了短舱外罩阻力,降低了机翼升力;动力条件较通流条件,F6短舱外罩阻力降低,机翼升阻力特性在零度攻角时有所提高.      

CRM翼/身/平尾组合体模型高阶精度数值模拟

基于五阶空间离散精度的WCNS格式,开展了CRM翼/身/平尾组合体模型的高阶精度数值模拟,以评估WCNS格式对复杂外形的模拟能力以及典型运输机巡航构型阻力预测方面的精度。首先依照DPW组委会提出的网格生成指导原则,利用ICEM软件生成了粗、中、细3套网格,对应的y⁺分别为1、2/3和4/9。通过对CRM模型的计算和分析,研究了网格规模对气动特性、压力分布和翼根后缘局部分离区的影响。通过与DPW4统计结果和部分实验结果的对比,高阶精度数值模拟结果表明,本文的气动特性计算结果与统计平均结果吻合较好;网格密度对机翼内侧的展向站位压力分布影响较小,对机翼外侧展向压力分布影响较大;网格密度对翼根后缘局部分离区略有影响。     

后缘弧形襟翼的气动特性研究

针对格尼襟翼低头力矩较大的问题,提出了后缘弧形襟翼概念作为格尼襟翼的改进构型,该构型通过襟翼背风面形成的局部高压区域来提供附加抬头力矩。基于CFD数值计算模型分析了弧形襟翼的气动力特性和襟翼附近的流场特征。计算结果表明:弧形襟翼的低头力矩明显低于相同高度的格尼襟翼;相比格尼襟翼,升力系数和阻力系数虽略有减小,但升阻比增大;尖部向流场中延伸能够进一步降低襟翼的低头力矩。相比格尼襟翼,弧形襟翼更有利于在飞行器设计中实际应用。     

带升力风扇飞翼布局开口机翼气动特性研究

为获得带升力风扇飞翼布局飞机因开口对机翼气动特性的影响规律,对平飞状态下开口机翼的气动特性进行三维气动仿真,分析升力系数、阻力系数和力矩系数随来流速度和迎角变化的特性.结果表明:随着来流速度增大,阻力和力矩呈上升趋势;来流速度一定时,随着迎角加大,升力系数增大,阻力系数先减小后增大;随着迎角增大,力矩系数先减小后增大再减小,且一直产生低头力矩.     

基于CFD和CSM耦合的通用静气弹分析方法

提出了一种适用于有限元精细化建模的流固耦合插值点选择方法,通过RBF(径向基函数)方法实现流固耦合面的数据交换,实现了基于CFD/CSM(computational fluid dynamics/computational structural mechanics)耦合的通用非线性静气弹分析方法。以HIRENASD(high Reynolds number aero-structural dynamics)风洞试验模型为验证对象,数值结果很好地与风洞试验结构变形、气动压力分布吻合,验证了所发展非线性CFD/CSM耦合静气弹求解器的精度。详细研究了HIRENASD模型在大迎角(AOA)流动下的静气动弹性特性,以及该模型弹性变形对机翼气动特性影响规律。研究表明:HIRENASD弹性模型变形后其升力小于刚性模型;在小迎角范围内刚性、弹性模型升力差随迎角增大呈线性增长;当迎角大于4°后,升力差先减小后基本保持不变,呈非线性关系。      

鸭式飞机矢量喷流对大迎角气动特性的影响

针对单发鸭式布局飞机，通过低速风洞试验，研究了矢量喷流对飞机大迎角气动力的影响特性。研究结果表明：发动机喷口直径变大使得飞机大迎角升力和阻力系数增加，并产生低头力矩系数。喷流使得飞机大迎角升力和阻力系数明显增加，并产生低头力矩系数；大喷口状态喷流影响比小喷口状态高50%左右。发动机喷管上/下偏转时，矢量喷流对飞机上下表面气流诱导不对称，喷管上偏减小升力和阻力系数、产生抬头力矩系数，喷管下偏增加升力和阻力系数、产生低头力矩系数，且喷管下偏影响明显比上偏大。在此基础上，基于数值模拟结果对喷流与飞机主流的相互作用机理进行了分析。     

基于非结构混合网格的CHN-T1标模气动特性预测

本文使用自行研制的基于非结构混合网格的亚跨超声速流场解算器MFlow,对AeCW-1提供的客机标模CHN-T1进行了数值模拟研究。介绍了非结构混合网格的生成情况,重点分析了网格收敛特性、压力分布、气动特性曲线、流动分离等。计算得到了近似线性的网格收敛特性。随着网格加密,对激波和分离气泡的模拟更精细。尾支撑对气动特性的影响非常明显,特别是对平尾气动特性有很大影响。机翼静气动弹性变形的影响主要是使升力系数和阻力系数减小。湍流模型的QCR修正对大迎角计算结果有较大的影响。计算结果表明,MFlow程序能够准确地预测客机标模的气动特性。     

动力效应对民机起飞构型气动特性影响的数值研究

采用在点对点多块结构化网格系统上求解三维可压缩雷诺平均N-S（Navier-Stokes）方程的方法,研究了发动机动力效应对民机起飞构型气动特性的影响.首先,通过涡轮动力模拟器风洞试验模型及民机高升力构型标准模型,对研究方法进行了验证,计算和试验结果的良好吻合说明研究方法用于模拟发动机动力效应及预测民机高升力构型气动特性是可行的.其次,针对分别安装通气短舱和动力短舱的某翼吊涡扇发动机民机起飞构型,研究了发动机动力效应对其气动特性的影响规律.结果表明:在通气短舱基础上预测的起飞最大升力系数、失速迎角以及阻力和力矩,在考虑动力效应后会产生明显的不同.动力效应可使外形最大升力系数增大并延迟机翼失速,但同时也带来了低头力矩增大和升阻比降低的不利影响.建议在对机翼及其高升力装置设计结果的数值校核中使用考虑进排气的动力短舱代替目前常用的通气短舱,这能为设计改进提供更可靠的依据.      

HiLiftPW-3高升力构型数值模拟

基于雷诺平均Navier-Stokes方程和拼接网格技术,采用MUSCL-Roe格式和Spalart-Almaras一方程湍流模型,对第3届高升力构型性能预测会议提供的两组高升力标模进行了数值模拟,主要目的是确认本文计算方法模拟复杂高升力构型的能力。研究内容主要包括高升力构型网格生成技术、网格收敛性研究及气动特性数值模拟。通过与JAXA （Japan Aerospace eXploration Agency）提供的测压、测力试验结果的对比分析,表明,在失速迎角之前,数值模拟得到的气动力系数和压力分布均与试验结果吻合;较好地模拟了局部外形变化引起的气动特性差量。本文建立的数值模拟方法对典型运输机三段翼布局的低速问题具有良好的适用性,可以为大飞机低速构型的气动设计及评估提供技术支撑。     

DLR-F6/FX2B翼身组合体构型高阶精度数值模拟

基于雷诺平均Navier-Stokes(RANS)方程和结构网格技术,采用五阶空间离散精度的加权紧致非线性格式(WCNS)和剪切应力输运(SST)两方程湍流模型,开展了DLR-F6和DLR-F6_FX2B 2种翼身组合体构型的高阶精度数值模拟,计算外形来自AIAA第三届阻力预测研讨会。主要目的是确认WCNS模拟跨声速典型运输机构型和预测局部构型变化引起的气动特性变化量的能力。在固定升力系数条件下,采用粗、中、细3套网格开展了网格收敛性研究,从气动力系数、压力系数分布、表面流态等方面研究了网格规模对DLR-F6和DLR-F6_FX2B翼身组合体数值模拟结果的影响;采用中等网格开展了来流迎角对2种翼身组合体气动特性的影响研究。通过与National Transonic Facility(NTF)的试验结果和CFL3D的计算结果对比,表明采用高阶精度计算方法得到了网格收敛的数值模拟结果,较好地模拟了DLR-F6翼身组合体局部修型引起的微小气动特性变化和翼身结合部流动特性的差异。     

融合式翼梢小翼减阻效应研究

有效减小诱导阻力对于飞机降低油耗、提高航程具有重要意义。针对某飞机翼身组合体构型,采用CFD数值模拟方法分析融合式翼梢小翼对飞机气动力特性的影响,尤其是其减阻效应;并给出翼梢小翼附近的空间流场。结果表明:带翼梢小翼后翼尖涡强度减弱,飞机阻力系数明显下降;固定升力系数0.5时,弯矩增加3.2%,阻力系数减小4.2%。     

翼梢装置对某机翼气动弹性行为影响研究

以大客某方案机翼为基本翼,通过数值模拟的方法研究了翼梢装置对机翼气动弹性特性影响,包括静气动弹性及颤振特性。其中通过CFD/CSD弱耦合求解的方法研究其静气动弹性响应,气动力计算采用面元法,结构响应计算采用结构有限元法,通过插值实现翼面气动力与有限元节点力之间的传递,以及有限元模型与气动网格之间的变形传递。对基本翼及带翼梢装置机翼静力学有限元模型局部修改得到动力学模型,应用MSC NASTRAN进行颤振特性分析。研究发现翼梢装置使得机翼的气动弹性特性不同程度均有降低,而不同翼梢装置对其影响又有所不同,可见,翼梢装置的设计在追求气动特性改善的同时必须关注其带来的结构特性的损失。     

舰载机折叠翼气动特性分析

折叠翼技术是舰载机与航母相匹配的关键,研究折叠翼的气动特性对舰载机飞行动态及安全性具有重要意义。基于离散化的思想对折叠翼进行建模,并利用计算流体力学CAE软件Fluent对其进行数值计算,从升力系数和阻力系数变化及压力分布分析舰载机折叠翼的气动特性。结果表明:当折叠角为75°时,外翼的升力系数及侧向力系数将达到最大,阻力系数最小;一旦折叠机构失效,外翼的折叠将导致机翼整体升力系数及阻力系数大幅度下降,气动特性变差;与折叠机构失效前相比,折叠后的机翼失速迎角变小,且在失速后升力系数有缓慢上升趋势。