非结构网格在积冰数值模拟中的应用

基于非结构网格对翼型表面的积冰进行了数值模拟研究.通过求解气体欧拉方程得到翼型的绕流流场,使用穿透型壁面边界模拟水滴对翼型表面的撞击,采用迭代法对水滴控制方程进行数值求解,得到翼型表面的水滴收集特性;采用不同的积冰冻结模型模拟霜冰和明冰的冻结过程;利用冰层时间推进法模拟积冰过程,基于积冰法向生长假设生成积冰外形.对NA...     

复杂积冰翼型气动性能分析

为了克服复杂外形难以划分结构化网格的问题, 采用结构/非结构混合多块网格对复杂积冰翼型的气动性能-进行了分析计算.在靠近积冰边界的内层采用了非结构网格, 在非结构网格外采用结构化网格以节省存储空间和计算时间.不同类型网格之间采用了点点对接的方式以简化通量重构的处理方式.整个流动区域采用SST湍流模型, 求解了雷诺平均N-S方程.数值计算结果表明结构/非结构混合多块网格在复杂积冰翼型气动性能计算上的可行性, 积冰对翼型的气动性能造成了严重的影响.      

基于非结构动态网格的翼型积冰过程数值模拟

讨论了积冰气象条件下翼型积冰过程的非结构动态网格数值模拟方法.外部无黏绕流流场采用基于非结构网格的分布式并行计算技术,空间离散采用VanLeer迎风有限体积格式,时间推进为5步显式龙格-库塔方法.采用四阶龙格-库塔方法求解过冷水滴运动轨迹方程,使用了一种合理、简单的物面水积冰量计算方法,提高了计算效率.使用非结构动态网格技术对积冰后的翼型网格进行快速调整.在对NACA0012翼型积冰形状进行校核验证后,对某型飞机在积冰气象条件下飞行时翼型积冰过程进行了数值模拟,得到了与参考文献结果相符的冰形特性.     

应用混合网格分析复杂积冰翼型气动性能

通过发展一种点-点对接的结构/非结构混合网格生成方法,避免了复杂积冰外形难以划分结构化网格的问题,并对复杂积冰翼型的气动性能进行了分析计算.在靠近积冰边界的内层采用了非结构网格以拟合复杂的边界,在非结构网格外采用结构化网格以节省存储空间和计算时间.为了检验网格对选用湍流模型的影响,整个流动区域分别采用SST和SA湍流模型,求解了雷诺平均N-S方程.数值计算结果表明SST模型更适于模拟复杂分离流动,积冰对翼型的气动性能造成了严重的影响.     

基于非结构动态网格的翼型积冰过程数值模拟

讨论了积冰气象条件下翼型积冰过程的非结构动态网格数值模拟方法。外部无黏绕流流场采用基于非结构网格的分布式并行计算技术,空间离散采用VanLeer迎风有限体积格式,时间推进为5步显式龙格 库塔方法。采用四阶龙格 库塔方法求解过冷水滴运动轨迹方程,使用了一种合理、简单的物面水积冰量计算方法,提高了计算效率。使用非结构动态网格技术对积冰后的翼型网格进行快速调整。在对NACA0012翼型积冰形状进行校核验证后,对某型飞机在积冰气象条件下飞行时翼型积冰过程进行了数值模拟,得到了与参考文献结果相符的冰形特性。     

基于多重网格方法的跨声速颤振数值模拟研究

使用多重网格的方法对二维翼型和三维机翼的跨声速颤振进行了数值求解.流场控制方程为N-S方程,湍流模型采用SA模型网格.计算网格采用结构网格,空间离散采用中心格式,使用双时间法进行时间推进.流场与结构之间的数据通过径向基函数(RBF)插值方法进行交换.通过耦合求解流场控制方程和结构运动方程得到Isogai(case A)机翼模型和Agard445.6机翼的跨声速颤振边界.通过与使用单重网格时得到的结果相比较得出,多重网格方法能够节省计算时间,提高颤振分析计算效率.     

一种强耦合Spalart-Allmaras湍流模型的RANS方程的高效数值计算方法

在工程实际中,一方程湍流模型或两方程湍流模型的求解通常和雷诺平均Navier-Stockes (RANS)方程的求解是解耦的,也称之为松耦合求解.在松耦合求解过程中,RANS方程和湍流模型方程通常采用不同的数值方法异步求解.这种求解方式很容易产生因两者计算精度不一致而引起的额外数值耗散.为了消除这种耗散,将RANS方程与Spalart-Allmaras模型方程耦合成一个系统方程——强耦合RANS方程,并发展了一种用于求解该系统方程的高效强耦合算法,其中对流项离散采用了Roe格式,时间项的离散采用了隐式LU-SGS(Lower-Upper Symmetric Gauss-Seidel)格式,为了提高计算效率,采用了三层V循环多重网格方法.通过翼型/机翼和振荡翼型/机翼等算例验证了本文发展的强耦合算法不仅具有较好的收敛性,而且计算精度明显优于松耦合算法,特别对于阻力的预测,强耦合算法更加准确.     

共轭梯度方法在求解不可压N-S方程中的应用

采用非结构网格上的SIMPLE(Semi-implicit Method for Pressure-linked Equations)算法、k-ε湍流模型求解了三维不可压Navier-stokes方程,对低速不可压粘性流场进行数值模拟,在求解压强方程时采用共轭梯度方法,求解其它变量方程时,采用预处理的BICG(biconjugate gradients)算法.对NACA0012翼型绕流流场和飞艇绕流流场进行数值模拟并对结果进行分析,取得较好的结果.     

非结构聚合多重网格法流场数值模拟研究

 在非结构网格上应用聚合多重网格技术来加速 Euler方程的收敛过程。设计了一种高效率的网格聚合方法。采用四重网格,在每一层网格上采用有限体积法进行计算。通过对多段翼型绕流的数值求解验证了多重网格加速收敛的高效性。     

基于伴随方程方法的非结构网格自适应技术及应用

对基于伴随方程的网格自适应技术的基本原理进行了研究,构造并实现了基于流场伴随方程的非结构网格自适应探测器,利用局部加密网格的办法,建立了基于Euler方程的无粘伴随自适应能力。对于粘性流动问题,提出首先应用无粘计算和伴随自适应技术从基准网格自动生成具有较好网格分辨率的非结构自适应网格,然后采用层推进法自动生成粘性计算网格的求解策略。采用RAE2822二维翼型和ONERA-M6三维机翼对所建立的无粘伴随自适应方法进行了验证,并应用所提出的方法和策略,对第二届国际涡流试验项目(VFE-2)中的65°尖前缘三角翼大攻角涡流场进行了数值模拟。对比计算结果表明,所建立的非结构网格伴随自适应方法是可行的和有效的。     

跨声速风扇叶片变形对气动性能的影响

为了研究跨声速风扇叶片变形对气动性能的影响,发展了流固耦合方法计算风扇叶片在气动力作用下变形的数值程序,耦合求解了三维可压缩雷诺平均的Navier-Stokes (N-S)方程和描述叶片振动的结构动力学方程,并模拟了跨声速风扇Rotor 67转子叶片在气动载荷下的变形,重点分析了由于叶片变形带来的气动性能的变化.      

一种高效的CFD/CSD耦合飞行器多学科优化设计方法

建立了一种基于CFD/CSD (computational fluid dynamics/computational structure dynamics)耦合的多学科优化设计方法.该方法采用结构有限元方程进行结构特性分析,同时采用数值求解N-S方程方法获得气动载荷,保证了设计结果的可信性,并采用基于代理模型的优化设计方法解决了采用高精度模型后计算量太大的问题.采用RBF(径向基函数)插值方法进行流固耦合界面数据插值,保证了传递过程中的能量守恒.利用建立的优化设计方法对一典型的翼身组合体进行了气动和结构优化设计,算例表明在满足约束条件的前提下优化后阻力和质量分别减小了5.0%和4.2%,证明了该方法的有效性.      

混合网格方法在栅格翼数值模拟中的应用研究

栅格翼是采用一种蜂房结构设计的升力面和控制面,和传统舵面相比具有铰链力矩小、能在较大攻角下保持升力等特点,并得到了越来越广泛的应用。由于栅格尾翼的特殊结构,网格生成具有很大的难度。本文采用了结构、非结构混合网格的办法,在弹体和栅格翼的物面附近区域,生成适合粘性计算的大长宽比的结构网格,弹体网格和栅格翼网格之间采用非结构网格进行填充,满足非结构网格和结构网格交接面的完全对接。本文基于结构/非结构网格体系采用有限体积方法求解NS方程,对不同舵偏角下的栅格翼构形进行了数值模拟,并通过实验结果对数值方法进行了验证。     

应用Delaunay图映射与FFD技术的层流翼型气动优化设计

采用非均匀有理B样条(NURBS)基函数属性建立了任意空间的自由式变形(FFD)翼型参数化方法,进一步结合基于Delaunay图映射技术建立了结构对接网格变形模式,通过粒子群优化(PSO)算法进行参数化方法、网格变形模式以及计算流体力学(CFD)数值模拟技术之间的整合,研究、构建了气动优化设计系统,并对某型层流理念设计的高空长航时(HALE)飞机基本翼型进行气动优化设计。气动特性目标函数评估方法中,边界层转捩数值模拟技术采用γ-Re_θt转捩模型耦合剪切应力输运(SST)模式湍流模型。优化设计后翼型气动特性表明:采用相关技术建立的层流翼型气动优化设计系统对于层流理念设计的HALE飞机翼型的设计具备较高的优化效率。     

基于CFD/CSD耦合方法的旋翼气动弹性载荷计算分析

为提高直升机前飞状态下旋翼非定常气动弹性载荷的预估精度,在旋翼气动弹性综合分析方法中引入旋翼CFD模块,建立了一套基于CFD/CSD松耦合分析的计算方法和程序。为高效解决流固耦合方法中由于桨叶挥舞、扭转等弹性变形带来的旋翼贴体网格变形问题,采用基于代数变换方法的网格变形技术,桨叶运动变形量和旋翼气动力信息通过流固交接面传递。旋翼流场分析方法中,主控方程采用耦合S-A湍流模型的Navier-Stokes方程,围绕旋翼流场的网格采用结构嵌套网格方法生成,无黏通量计算采用Roe格式,时间推进采用双时间法。旋翼结构分析中,考虑旋翼配平,基于Hamilton变分原理和20自由度Timoshenko梁模型求解弹性旋翼非线性运动方程。分别对CSD和CFD方法进行验证,在此基础上,计算了SA349/2旋翼桨叶在前飞状态下的非定常气动力、挥舞弯矩、摆振弯矩和扭转力矩,并与飞行测试数据进行了对比。计算表明:CFD/CSD耦合方法可以显著提高旋翼非定常气动弹性载荷的分析精度,精确捕捉桨叶表面压强峰值、激波位置等,表明本文发展的旋翼CFD/CSD耦合方法可以有效地运用到旋翼气动弹性载荷的预测分析中。     

双旋翼微型飞行器的非定常气动特性

发展了一种基于计算流体力学(CFD)的非定常气动特性预测方法,计算方法包括了动量源模型、预处理方法、非结构嵌套网格和Spalart-Allmaras(S-A)湍流模型等技术.通过计算悬停Caradonna-Tung算例和俯仰振荡NACA0012算例,验证了计算方法模拟双旋翼微型飞行器动态流场的有效性.数值模拟了双旋翼微型飞行器动态流场,给出了非定常气动系数的迟滞曲线,分析了缩减频率、前飞速度和螺旋桨转速对非定常气动特性的影响.计算结果表明:力矩系数迟滞效应随缩减频率、前飞速度和螺旋桨转速增大而增大,升力系数迟滞效应随缩减频率和前飞速度增大而增大,但基本不随螺旋桨转速变化而变化.      

基于 N-S 方程的飞机风挡雨刷气动力计算

风挡雨刷在飞行过程中易出现被气流吹离复位位置、飘起现象,阻碍飞行员视野,严重影响飞行安全。选取某型机风挡雨刷作为优化对象,采用了结构、非结构混合网格的办法,在机身壁面附近生成结构化网格,在雨刷机构表面生成非结构化网格,用 Interface 面将两部分结合生成雨刷最终计算模型。利用 Fluent 软件,基于 N-S 方程和 S-A 湍流模型,采用时间平均法对不可压缩流的参数进行时均化,计算了不同飞行状态下,风挡雨刷在复位位置以不同角度停放及增加扰流板时所受气动力,对风挡雨刷的外形及停放角度进行优化,比较分析了优化前后雨刷所受气动力。仿真计算结果表明,雨刷合理的停放角度、在刷臂上加装扰流板能有效引导风挡雨刷迎风面气流的方向,从而减小雨刷迎风面所受气动力。通过试飞验证了数值计算的可靠性及优化方案的可行性,为风挡雨刷预紧力的设置及复位位置的确定提供依据。     

倾转旋翼机多部件对机翼气动干扰的分析及优化

为研究倾转旋翼机各部件对机翼的气动干扰效应，建立了一套适用于倾转旋翼机流场CFD求解方法，提出并生成了一套由三棱柱/四面体组成的适用于倾转旋翼机多种飞行状态的非结构混合网格系统.以三维Navier-Stokes方程为主控方程，采用动量源方法进行倾转旋翼的模拟，并引入了高效的OpenMP并行加速技术等，提出了一套新型的动量源项添加及搜索方法.计算分析了倾转旋翼机旋翼/机身/短舱对于机翼气动特性的干扰影响，得出了一些对设计有指导意义的结论.采用基于径向基函数（RBF）的代理模型方法，根据机翼不同展向位置的干扰程度对各段翼型配置进行气动优化.优化结果表明:考虑气动干扰作用，在巡航速度下优化后的全机升阻比增长达到了36.78%.      

弹射座椅救生性能研究

针对座椅在铅垂面内弹射的简化模型,建立了系统的非线性运动方程,然后用推进阵面法建立了人－椅系统平面的非结构网络,并用Ｅｕｌｅｒ方程作为主控方程对某型座椅的气动阻 力系数,升力系数和气动力矩系数进行了验证性计算,最后对非线性运动方程在时域内进行积分,得到人－椅系统运动轨迹,运动速度和旋转角速度等性能参数,并与座椅弹射测试结果进行了比较分析,结果表明,用非结构网络计算的座椅气动特性与座椅吹风数据吻合,计算了弹射救生性能曲线与试验数据的趋势基本一致。     

基于CFD/CSD方法的蜻蜓柔性翼气动特性分析

给出了一种基于计算流体力学/计算结构力学（CFD/CSD）的双向流固耦合方法.通过交替数字二叉树（ADT）搜索技术识别流固网格之间的宿主-受体关系.采用局部插值算法完成两套网格系统之间的数据交换,并使用Delaunay图映射方法来完成气动网格的移动.将自编的非线性结构有限元程序、接口程序与南京航空航天大学（NUAA）微型飞行器中心的流体计算程序3D2MUFS相连接,应用于蜻蜓柔性翼拍动飞行的气动计算中.计算结果表明:柔性变形使得蜻蜓翼的时均举力系数从0.31提高到0.53,时均推力系数从0.07提高到0.13,证实了柔性变形能改善扑翼的气动性能.