绕钝体分离流动的非定常数值模拟研究

利用非定常RNGk-ε模型对三种典型的钝体绕流问题进行了研究,这三种流动分别为:方柱绕流、圆柱绕流和绕立方柱流动。研究表明了采用非定常雷诺平均方法可以研究钝体绕流的非定常流动,但所获得的结果各不相同。对于分离角固定的方柱绕流和绕立方柱流动,数值模拟结果和实验以及大涡模拟结果吻合较好,对于分离角不固定的圆柱绕流,随着雷诺数的增加,数值模拟结果和实验结果相差越大,分析认为这是由于壁面模型无法准确预测流动中不固定的分离角所致。改进了绕立方柱流动的模拟方法,采用非定常模拟方法后获得的结果要比定常模拟方法得到的结果有了根本改善。     

激波绕流与激波碰撞的数值模拟

本文应用基于Poe平均构造的最小耗散的扩散与反扩散守恒差分格式，数值模拟单激波绕流三我形障碍物，及等强度激波碰撞和不等强度激波碰撞的二维非定常复杂流场。计算结果表明，该格式具有很高的单波精度及较好的激波与接触间断分辨度，可直接用于复杂流场的数值模拟。     

用N—S方程计算翼型非定常粘性大攻角绕流

本文给出了用ＮＳ方程对翼型非定常粘性大攻角绕流的数值模拟。控制方程为二维时均可压缩完全ＮＳ方程；湍流模型采用双层代数涡粘性模型。使用近似因式分解ＡＤＩ差分格式离散求解，网格是用保角变换方法生成的相对翼型固定的Ｃ型网格。本文给出两类典型非定常绕流数值模拟结果：翼型过失速常攻角周期流动和大攻角强迫俯仰谐振非定常绕流。并与国外的实验和计算结果进行了比较，表明了本方法准确、高效的特点。     

用动态非结构重叠网格法模拟三维多体相对运动绕流

 发展了动态非结构重叠网格技术。本文提出了基于阵面推进的相邻单元搜索算法和适合于高阶空间离散格式的网格间插值关系定义方法，前者提高了非结构重叠网格在模拟多体间具有大幅相对运动的非定常问题时网格处理的效率，而后者保证了网格间边界附近流场求解精度；本文非定常Euler方程的求解采用有限体积法，通量计算采用Osher格式，时间积分采用双时间步长推进。为验证本文提出的方法，分别对定常M6机翼绕流和机翼下外挂物投放过程的非定常绕流进行了模拟,结果表明本文提出的方法提高了非结构重叠网格在模拟多体间具有大幅相对运动问题时的效率，保证了网格边界处的精度，具有较高的工程应用价值。     

电磁力控制翼型体绕流的数值研究

为考察包覆电磁激活板的翼型体绕流特性,采用数值模拟的方法研究了低雷诺数下电磁力对翼型体流场结构和绕流过程的影响,探讨了不同迎角和不同电磁力大小对翼型体尾流涡街形态和流体边界层结构的影响。计算结果表明,作用于翼型体表面的电磁力可以明显改变流体边界层的结构,抑制边界层分离,并有效提高翼型体的升力。     

基于黏性涡域法的二维钝体绕流模拟

针对传统计算流体力学方法需要划分网格且计算耗时长等问题，一种确定性涡方法(黏性涡域法)逐渐发展并日臻成熟。本文利用该方法计算模拟了二维钝体的流场、气动系数及斯特劳哈尔数，并与其他成熟CFD软件计算结果对比来验证模拟精度。黏性涡域法的模拟过程为生成涡量、涡量演变及气动计算。文章首先推导了该方法各过程的基本方程，然后基于Lua语言编制了圆柱和方柱绕流计算程序，最后利用Gnuplot实现流场可视化。利用黏性涡域法模拟不同雷诺数下圆柱绕流、不同无量纲频率和无量纲振幅的振动圆柱气动力时程、不同风向角下的方柱绕流，结果表明：雷诺数为2×10²～1×10⁵时，基于黏性涡域法的圆柱绕流气动计算结果与其他方法的计算结果基本一致；基于该方法的不同风向角下方柱绕流气动计算结果与其他文献结果趋势相同。     

吸气条件对圆柱非定常分离流影响的数值研究

本文根据二维可压非定常Ｎ－Ｓ方程，基于Ｂｅａｍ ａｎｄ Ｗａｒｍｉｇｎ及Ｓｔｅｇｅｒ的方法，在圆柱绕流出现周期分离后，施加吸气条件，研究气对于圆柱非定常绕流的影响，计算结果表明，在恰当地区域施加恰当的吸气条件，可以希望通过减小压差阻力来减小总阻力，吸气还有抑制尾迹中非对称流动的作用，从而使圆柱所受的垂直于来流方向的升力的变化幅度减小，变化频率减小，但是吸气会增加磨擦阻力，过份的吸气会导致摩阻增大太     

真实飞行器无粘绕流直角切割网格数值模拟

本文采用直角切割网格，利用MGAERO无粘绕流解算器对三角翼真实飞行器带或不带有外挂物气动特性进行数值计算。计算马赫数为0．6～1．2。研究了网格密度对飞行器气动力特性的影响，同时与风洞实验结果进行了比较，符合程度令人满意。     

绕流问题的区域收缩法及其数值解

本文提出了一种适于绕流问题的数值方法,它考虑了流场中不同区域内不同的物理现象,给出了外部边界条件的良好逼近。使用这种方法可以节省大量的计算机存贮单元和计算时间,也可以推广到三维绕流问题的数值计算。     

高超声速验证飞行器助推分离段流场数值研究

为了模拟有相对运动的多体非定常绕流，发展了基于弹簧近似的非结构动网格方法及耦合动网格的Euler解算器，作为方法验证，模拟了俯爷振动的NACA0012翼型绕流，计算结果与实验及文献结果非常接近。采用二维近似外形，对高超声速验证习行器助推分离段流场进行了数值研究，得到了不同时刻由于多体相对运动形成的干扰流场结构以及分离过程的气动力参数。     

基于非结构重叠网格的二维外挂物投放模拟

本文发展了基于非结构重叠网格的动态网格技术,并提出了新的相邻单元搜索算法。非结构重叠网格集非结构网格和重叠网格的优点于一身,适合处理复杂外形和多体相对运动问题。新的相邻单元搜索算法在碰到物面时能继续搜索下去,搜索效率大大提高。非定常Euler方程的求解,采用有限体积法,时间上采用双时间步长推进。为验证算法的有效性,先对绕振荡NACA0012翼形的流动进行了模拟;然后通过耦合求解刚体三自由度运动方程,模拟了二维外挂物投放过程,结果表明非结构重叠网格在处理多体相对问题时具有优越性。     

基于预估校正和嵌套网格的虚拟飞行数值模拟

 针对导弹虚拟飞行数值模拟问题,发展了空气动力学/飞行力学数值计算方法和软件。控制方程为非定常雷诺时均Navier-Stoker(RANS)方程和刚体六自由度运动方程;流场求解器为有限体积法结构网格求解器,时间推进采用双时间步法,湍流模型为Spalart-Allmaras一方程模型;采用Adams预估校正法实现飞行力学方程与流场控制方程的耦合计算;使用嵌套网格方法模拟多体运动。首先模拟了美国国家航空航天局(NASA)窄条翼导弹模型纵向虚拟飞行,研究耦合方式和时间步长的影响。仿真结果表明,双时间步三阶Adams耦合方法,同等精度下可以显著增大时间步长,缩短仿真时间。最后,采用该方法模拟了导弹自由摇滚特性和纵向虚拟飞行,模拟结果与试验值吻合较好。     

基于动态非结构重叠网格法的直升机前飞非定常流场数值模拟研究

 发展了非结构重叠网格技术，提出适合于直升机旋翼运动的动态非结构重叠网格算法，该算法使得在采用非结构重叠网格模拟直升机流场时实现程序自动处理网格，不需任何人工干预。流动控制方程采用非定常可压缩Euler方程，空间离散采用Jameson有限体积法，时间推进采用双时间步长法，为加速收敛，在双时间步长内循环迭代中采用预处理技术。为验证本发展的预处理解算器以及动态非结构重叠网格的算法的可信性和有效性，开展一系列的数值实验，模拟了不同速度范围的NACA0012翼型绕流,Caradonna旋翼悬停和前飞流场以及Robin直升机机身-旋翼组合的前飞非定常流场，结果表明采用本文算法得到的结果是可信的，网格处理是自动高效的，可以应用于直升机全机流场模拟，为直升机设计提供可信指导数据。     

基于非结构运动对接网格的旋翼前飞流场数值模拟

发展了一种基于线性弹簧理论的非结构运动对接网格方法,通过求解非定常Euler方程模拟了旋翼的前飞流场。利用运动区网格的转动模拟旋翼的旋转运动,并运用线性弹簧动态网格技术处理旋翼的挥舞和变距运动,方程求解中采用一种精度较高的插值和重构方法。算例表明,本文的方法能够有效地模拟旋翼的前飞流场。     

基于非结构动网格的非定常流数值模拟方法

给出了基于非结构动网格的非定常流动问题的数值模拟方法,重点说明两种动态非结构网格方法的原理以及在非定常流模拟中的算法.一种是基于非结构网格的重叠网格方法,该方法结合了非结构网格方法和重叠网格方法的优点,以物面距为准则参数对网格节点进行分类,进而建立重叠网格的网格问插值关系,通过采用动态非结构重叠网格算法,可以处理多体间...     

基于非结构动网格的跨声速非定常流场求解

对用于非结构动网格生成的弹簧近似方法进行了研究.通过采用顶点弹簧方法,分析研究了弹簧倔强系数的取值,同时通过引入防扭转、防挤压倔强系数和边界修正,对标准弹簧近似方法进行了改进.采用改进后的方法并通过耦合求解了基于ALE(Arbitrary Lagrangian-Eulerian)描述的Euler方程,模拟了做谐和振动的刚性无后掠机翼及某翼身组合体的跨声速绕流.结果表明,改进后的方法可以大大提高网格变形能力和网格质量,且计算结果与实验结果吻合良好.     

共轴式双旋翼悬停流场和气动力的CFD计算

考虑到共轴式双旋翼流场特征高度复杂、桨叶承受非定常气动载荷的特点,为更好地预测共轴式双旋翼的气动特性,把非结构嵌套网格方法和网格的自适应技术相结合,发展了一套适合于共轴式双旋翼流场数值计算的求解器.在该求解器中,采用非结构嵌套网格方案来模拟桨叶之间存在的相对运动,自适应网格技术用来捕捉尾迹对流场和气动特性的影响,求解惯性坐标系下的非定常N-S主控方程来模拟流场的非定常特性.应用该求解器,首先计算了有试验结果可供对比的一副试验旋翼的诱导速度场分布,在此基础上,计算了共轴式双旋翼的桨尖涡轨迹和拉力分布特性,并与单旋翼的计算结果进行了对比和分析,得到了一些有意义的结论.     

基于非结构嵌套网格方法的旋翼地面效应数值模拟

 建立了一个基于非结构嵌套网格的流场求解器，用来精确模拟复杂的旋翼近地流场，为更好地分析地面效应(IGE)对旋翼气动特性的影响提供一套计算方法。在该求解器中，控制方程采用惯性坐标系下的非定常N-S方程，空间方向上采用二阶迎风格式，并用背景网格的一个面模拟地面作用，以方便地面边界条件的处理。应用所建立的模型，首先针对有实验结果可供对比的旋翼无地面效应（OGE）流场进行了数值模拟，以验证计算方法。然后着重计算了IGE下流场中的桨尖涡空间位置和旋翼拉力增益，并对旋翼在小速度前飞状态下的近地流场及地面涡形成过程进行了模拟和分析。在此基础上，得出了一些有意义的结论。     

Numerical Simulation of Unsteady Flow Around Forward Flight Helicopter with Coaxial Rotors

Three-dimensional unsteady Euler equations are numerically solved to simulate the unsteady flows around forward flight helicopter with coaxial rotors based on unstructured dynamic overset grids. The performances of the two coaxial rotors both become worse because of the aerodynamic interaction between them, and the influence of the top rotor on the bottom rotor is greater than that of the bottom rotor on the top rotor. The downwash velocity at the bottom rotor plane is much larger than that at the top rotor plane, and the downwash velocity at the top rotor plane is a little larger than that at an individual rotor plane. The downwash velocity and thrust coefficient both become larger when the collective angle of blades is added. When the spacing between the two coaxial rotors increases, the thrust coefficient of the top rotor increases, but the total thrust coefficient reduces a little, because the decrease of the bottom rotor thrust coefficient is larger than the increase of the top rotor thrust coefficient.     

基于直接力控制的阵风响应及阵风减缓研究

考察某无人机模型的阵风气动响应特性和基于直接力操纵的阵风减缓效果.在求解非定常Euler方程时引入“网格速度”方法模拟阵风边界条件,通过动态嵌套网格技术实现舵面运动.首先对NACA0006翼型迎角阶跃型阵风的气动力响应进行了验证计算,计算结果与理论结果、文献计算结果吻合良好.在此基础上对某无人机模型在迎角阶跃型、One-minus-cosine型阵风作用下的气动力响应过程进行了数值模拟,并分析评估了不同舵面运动方式对飞机气动性能的影响.最终研究比较了不同舵面运动方式下阵风气动响应的减缓效果.结果表明:通过合理设计舵面运动,可以有效达到抑制阵风引起的非定常气动干扰的目的.