反舰导弹三维流场数值模拟

三维流场的数值模拟是国外军事技术发达国家在飞行器设计过程中所采用的重要方法.本文建立了导弹的真实三维几何外形,然后结合实验数据,采用有限体积法,对国产某典型反舰导弹在高亚音速、不同攻角时的三维流场进行了数值模拟.计算出了导弹在不同工况下的阻力系数、升力系数和力矩系数.把计算所得结果与实验值进行了比较,用理论公式对计算结果进行了对比分析,并简单地分析了偏差的原因和解决的可行措施.文中所用的残差光顺技术和网格自适应技术能有效地提高计算的收敛速度和精度.     

飞机三维结冰模型及其数值求解方法

 对飞机结冰外形进行纯三维数值模拟,是目前飞机结冰预测中的一大难点。为了建立三维结冰数值模拟方法,基于Messinger的二维结冰模型,提出了一种考虑液态水溢流效应的三维结冰计算模型,并针对该模型建立了表面单元内溢流水流动的分配方案、发展了相应的迭代求解方法。采用本文方法对MS-317后掠翼结冰进行了计算,并与实验和Lewice3D的计算结果进行了对比。研究结果显示：本文方法具有较好的收敛性,计算的霜冰和明冰外形均与Lewice3D计算的冰形一致;对于霜冰和结冰时间较短的明冰,本文计算的冰形与实验吻合较好;对于结冰时间较长的明冰,本文计算的冰形与实验对比还有一定差异,但冰生长的总体趋势和大致体积与实验一致。同时,对比了MS-317无后掠翼与有后掠翼结明冰的外形,发现机翼后掠导致的三维溢流效应对结冰外形有明显影响,因此,对于三维结冰分析,如果用二维截面的结果来代替三维结果,其合理性还需验证。     

三维可压和不可压缩流动数值模拟的一种预处理多重网格方法研究(英文)

将矩阵预处理方法与多重网格方法结合,发展了适用于机翼及复杂外形气动数值优化设计的高效数值模拟方法和程序。在不改变定常流动解的前提下,对 Navier-Stokes 方程的时间导数项实施 Choi-Merkle 矩阵预处理,从而改善可压缩控制方程在低速情况下的系统刚性。Choi-Merkle 预处理最初用于二维低马赫数粘性流动数值模拟,本文将其推广应用到计算三维流动,并针对提高跨音速流动计算收敛性进行了相应的改进。预处理后控制方程采用中心格式有限体积进行空间离散,Runge-Kutta 混合多步方法进行显式时间推进求解,并应用 FAS 多重网格方法加速收敛。采用所发展的方法和程序数值模拟了绕 M6 机翼和 M100 机翼的低速和跨音速流动(马赫数从0.01 到 0.839)以及绕 F4 翼身组合体流动。计算结果与实验值吻合良好,验证了所发展的方法的正确性。算例研究表明所发展方法较大幅度地提高了三维流动数值模拟的效率,且同时适用于可压和不可压流动计算,在机翼及复杂外形气动优化设计方法具有应用前景     

A Preconditioned Multigrid Method for Efficient Simulation of Three-dimensional Compressible and Incompressible Flows

将矩阵预处理方法与多重网格方法结合,发展了适用于机翼及复杂外形气动数值优化设计的高效数值模拟方法和程序。在不改变定常流动解的前提下,对 Navier-Stokes 方程的时间导数项实施 Choi-Merkle 矩阵预处理,从而改善可压缩控制方程在低速情况下的系统刚性。Choi-Merkle 预处理最初用于二维低马赫数粘性流动数值模拟,本文将其推广应用到计算三维流动,并针对提高跨音速流动计算收敛性进行了相应的改进。预处理后控制方程采用中心格式有限体积进行空间离散,Runge-Kutta 混合多步方法进行显式时间推进求解,并应用 FAS 多重网格方法加速收敛。采用所发展的方法和程序数值模拟了绕 M6 机翼和 M100 机翼的低速和跨音速流动(马赫数从0.01 到 0.839)以及绕 F4 翼身组合体流动。计算结果与实验值吻合良好,验证了所发展的方法的正确性。算例研究表明所发展方法较大幅度地提高了三维流动数值模拟的效率,且同时适用于可压和不可压流动计算,在机翼及复杂外形气动优化设计方法具有应用前景     

基于SST湍流模型的二维操纵面空化流场研究

操纵面是安装在超空泡航行器尾端的关键部件,它的基本功能是提供稳定力矩和操纵力矩。文章基于k-ω模型的SST(剪切应力输运)湍流模型,对典型二维操纵面的空化流动进行数值模拟,根据空化流仿真计算的结果,得到了超空泡外形、阻力特性以及升力特性的变化规律,并与水洞实验进行对比。研究结果表明,仿真结果能够很好地与水洞实验保持一致,说明所采用的湍流模型和数值方法能够准确预报空化流问题,为进一步研究空化流问题提供了一定的参考依据。     

襟翼表面粗糙度对多段翼型升力的影响

在襟翼上表面填充多孔吸声材料,虽然缝道的噪声有所降低,但翼型升力系数也有较大损失。为了确定是否因多孔材料造成翼面粗糙度增加而降低了升力系数,也为了进一步确定数值模拟中粗糙表面简化模型的适用性,通过风洞试验和数值模拟进行研究。数值模拟采用沿流向三角形凸起模拟吸声材料的粗糙度,采用ANSYS-Fluent的二维模块、非结构网格和k-ωSST湍流模型进行计算。风洞试验在某NF-3低速翼型风洞中完成,采用等弦长的二维测压模型,光滑模型表面是通过在吸声材料上覆盖光滑薄膜实现的。结果表明:可以采用沿流向三角形凸起模拟吸声材料的粗糙度对升力特性的影响;吸声材料的粗糙度影响翼型的升力系数。因此,若通过吸声材料降低缝道流动噪声,需要进一步研究吸声性能好且表面比较光滑的材料/结构。     

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

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

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

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

壁面粗糙度湍流扩展模型及流动数值模拟

采用有限体积法数值求解雷诺平均Navier-Stokes流动控制方程,数值模拟表面粗糙度对气动特性和热力学特性的影响。通过修正Knopp的基于切应力输运k-ω湍流模型的粗糙度湍流扩展模型中边界条件ω的值,提高了壁面附近的涡黏度,增加了总壁面切应力,从而引入壁面粗糙度的影响。通过对不同粗糙度的平板流动和NACA 65_2215翼型绕流进行数值模拟,分析粗糙度对表面摩阻系数、斯坦顿数和升力系数的影响。计算结果表明:粗糙度使得平板表面摩阻系数和斯坦顿数增大,翼型流动分离提前,升力系数减小。流经粗糙壁面的流动处于过渡状态时,改进的粗糙度湍流扩展模型改善了原模型速度偏移量与经验公式相比偏小的缺点,同时,表面摩阻系数、斯坦顿数和升力系数与实验结果吻合得更好。该改进模型不需要对壁面附近网格进行加密,能够快速、准确地模拟流经粗糙壁面的流动。     

基于eN转捩预测方法的增升装置失速特性数值模拟研究

基于三维可压缩Navier-Stokes方程耦合e~N转捩预测方法,对30P30N二维多段翼型和DLR-F11三维增升装置进行了数值模拟研究。计算结果表明,e~N方法计算得到的转捩位置与试验值符合较好,增升装置翼面特别是主翼和襟翼下翼面存在大范围层流流动,为提高增升装置绕流计算准度,计算中考虑转捩的影响是非常有必要的。相比全湍计算,e~N转捩计算升力系数与试验值符合更好,提高了增升装置大迎角失速特性计算准度。     

高升力标模确认计算研究

针对高升力标模NASA Trapwing全展襟翼构型展开确认计算研究。采用N-S方法,预测高升力模型的流动状态,分析高升力的流动现象及机理,评估自主研发的WiseManPlus软件及WoF90软件预测高升力流动的能力。通过对网格生成、计算策略及数据分析等各个环节进行分析,制定出一套适用于高升力模型工程应用的CFD计算方法,得到与实验值吻合较好的计算结果。结果表明,目前所采用的结构网格求解器WiseManPlus软件及混合网格求解器WoF90软件计算精度相当,适用于模拟高升力流动。     

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

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

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翼身组合体局部修型引起的微小气动特性变化和翼身结合部流动特性的差异。     

基于N-S方程的跨声速翼型多目标多约束优化设计

本文将粘性流场分析和数值优化方法耦合起来，发展了一种跨声速翼型设计方法，用以提高翼型在一个或多个设计点、在多种约束条件下的气动性能。由粘性流场分析程序计算得到的升力、阻力等气动参数构成目标函数，数值优化程序对其进行最小化。粘性流场分析采用了雷数平均Ｎ－Ｓ方程，这比过去翼型设计中使用的全速势方程或Ｅｕｌｅｒ方程更能模拟流动的本质，因而设计结果的可靠性大大提高了。优化算法采用传统的拟牛顿法（Ｑｕａｓｉ     

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

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

垂直/短矩起降飞机机翼内埋式风扇布局气动特性分析

以垂直/短距起降飞机过渡飞行状态为背景,针对机翼内埋式风扇布局的自由来流/风扇喷流混合型流动,基于结构/非结构混合网格使用CFD方法进行了非定常数值模拟和分析.首先使用滑移网格技术对NASA涵道螺旋桨进行算例验证,其时均计算结果与实验值的误差为5.3%,证明了计算方法的可靠性和准确性,然后数值模拟了机翼内埋式风扇布局在不同迎角下的气动性能.结果表明:风扇喷流在机翼上产生了特有的“抽吸”和“堵塞”效应,引起了机翼总升阻力的显著增加,升力最大增量达到干净机翼升力的2.6倍,阻力最大增量为干净机翼阻力的3.2倍,混合流场在机翼后缘引起了升力损失并卷起对涡.      

钝化前缘乘波布局及其一体化构型气动特性

以最大升阻比为优化目标,在锥型流场中优化设计出乘波布局,并考虑高超声速飞行器的防热需求,对乘波布局进行钝化设计,利用数值模拟和风洞实验两种手段,研究钝化前缘乘波布局的气动特性.结果表明:在一定钝化半径内,随着钝化半径的增加,乘波构型的升力特性变化仅为2%,但阻力特性增加近3倍,升阻比降低了将近50%.尽管如此,为了钝化乘波布局,仍维持了较高的升阻比,升阻比为3左右.同时,以二维顶压式进气道为基础,在多级楔锥组合体流场中,设计出满足超燃发动机进气要求的乘波前体/进气道一体化构型,并进行前缘钝化设计.针对一体化构型进行了数值验证,结果表明:此类一体化构型升阻比大于2.6,同时发动机总压恢系复数保持在40%左右,满足进气道的要求.      

螺旋桨滑流对增升装置气动特性影响研究

采用准定常N—S方法,通过单独螺旋桨模型计算结果与非定常、试验结果的对比,验证了计算螺旋桨滑流的可行性。重点对带两台螺旋桨发动机的增升装置流场特性进行了计算分析,研究了螺旋桨滑流对飞机流场和气动性能的影响。研究结果表明：螺旋桨滑流明显改变了机翼表面的压力分布,使飞机升、阻力系数增大。为螺旋桨飞机设计提供了一种快速、有效的数值模拟手段,在工程领域具有一定的应用前景。     

翼型舵面偏转非定常流动数值模拟

针对NACA0012翼型舵面偏转问题,数值模拟了不同参数对翼型气动特性的影响。基于非结构动网格技术,采用ALE有限体积描述下的二维可压缩非定常N-S方程,计算通量采用Vanleer格式、时空二阶格式,利用Venkatakrishnan限制器抑制数值振荡。非定常计算结果表明,NACA0012翼型绕1/4弦点作周期性俯仰振动的升力系数和俯仰力矩系数结果与实验数据吻合良好,验证了数值方法的准确性;在翼型舵面表面有分离区产生,升力系数和俯仰力矩系数形成滞回环,在亚声速情况下,滞回环幅值较小,进入超声速阶段以后,幅值增大,随着翼型间缝隙宽度逐渐增加,翼型升力系数和俯仰力矩系数与无缝翼型相比逐渐降低。