一种改进的非定常激波装配算法

基于非结构动网格技术和格心型有限体积方法,提出一种改进的非定常激波装配算法,进一步拓展了其在包含有运动激波的非定常流场的应用范围。首先,针对激波在直/曲壁面传播这类问题,分别建立了壁面间断节点的运动模型;其次,为保证激波在大范围运动时装配阵面不产生失真,基于Bézier曲线拟合方法实现了间断节点分布的自动重构;接着,通过嵌入局部网格自动重构模块,提高了算法的计算效率和自动化程度;最后,对于激波相交点的运动,设计了一种根据位移推算速度的方法进行装配。数值算例表明,所提算法能够有效地处理激波传播问题,相比激波捕捉方法可以提取更多的流场信息,同时可以获得流场间断更加直观清晰的图谱。     

面向三维激波问题的装配方法

给出了一种基于非结构动网格技术的三维激波装配方法。在该方法中,三维激波面由被标记为激波属性的网格点连接构成,标记为激波属性的网格点称为激波点。激波点具有两组参数分别代表激波的上下游,利用激波点上下游参数求解R-H关系式获得激波点运动速度。非结构动网格技术的使用允许激波大幅度运动,降低了对初始激波位置的要求。通过引入网格属性定义避免了对计算网格进行分区,增加了装配激波的灵活性。通过球柱体绕流问题验证了该三维装配方法的合理性,针对三维激波装配中比较困难的交点装配问题,通过对三维激波反射以及三维激波相交等算例进行研究找到了可用的三维激波交点运动速度的确定方法,保证了激波运动过程中交点运动与流场求解之间的相容性,获得了相应的装配结果。     

动态间断装配法模拟斜激波壁面反射

基于非结构动网格技术和边界装配思想提出了动态间断装配法,该方法能够应用于求解含有间断的流动问题。无论入射激波还是反射激波都是作为边界进行处理,激波运动速度由兰金-许贡纽(Rankine-Hugoniot)关系确定。激波作为动网格的一部分,其运动由动网格技术实现。采用该方法模拟了超声速二维流场中激波与壁面相交问题,并且与捕捉法进行比较,二者的流场结构符合良好,但是在细节上还是存在明显差异。通过对流动结构的分析,得出采用装配方法得到的流场要优于捕捉方法的结论。激波壁面反射的问题模拟,也说明了边界激波装配方法对于复杂的激波相交问题是具有处理能力的。     

基于非结构动网格的非定常激波装配法

从ALE方程出发,通过修正体积的新思路构建离散几何守恒律,基于非结构动网格的信息传递方法得到网格重构以后新网格的流场参数,采用一阶格式捕捉法得到的流场确定装配法的初始激波位置,成功地把激波装配法应用到非定常流动的模拟。模拟超声速飞行环境下整流罩分离的非定常流固耦合过程,在马赫数6条件下装配法和捕捉法的计算结果完全一致,随着马赫数增加,二阶精度格式的捕捉法需要进行调整限制器参数、时间步长等人工干预,但是装配法直到马赫数20计算过程没有出现异常。该基于非结构动网格技术的激波装配方法,与传统的激波装配方法相比,能够用于复杂外形产生的非规则形状激波和非定常流动产生的运动激波,解决二阶精度有限体积法用于高马赫数飞行条件下多体分离过程的流动模拟时遇到的稳定性问题。     

LU-AUSMPW混合格式在可压缩无粘和粘性流动数值分析中的应用

本文首先将 AUSMPW格式与三阶 MUSCL格式融合, 给出了其在任意曲线坐标下的三维形式, 并将其与 LU-SGS格式结合, 应用于可压缩 Euler和 Navier-Stokes方程的求解。其次, 构造了一种新的通量限制器。最后, 为了验证 LU -AUSMPW混合格式的性能, 将平面叶栅跨音速无粘流动以及喷管超音速粘性流动作为算例。本文计算结果与文献计算结果和实验数据相符很好, 表明采用 LU-AUSMPW混合格式数值模拟可压缩流场, 具有较高的计算精度、较快的收敛速度和良好的稳定性。      

Jameson有限体积法对非结构网格推广的改进

本文将Ｄｅｌａｕｎａｙ法用于复杂物形的非线性结构网格生成,并利用黄明格提出的布点原理避免边界穿透。绕流的计算采用Ｅｕｌｅｒ方程和Ｊａｍｅｓｏｎ的有限体积法。本文对三维问题提出新的激波敏感器,对方法做了改进,并用于ＡＧＡＲＤ Ｍ６机冀、导弹模型和飞机模型的绕流计算。     

二维非结构网格的一个TVD型有限体积方法

考虑双曲型守恒律方程,对二维非结构三角形网格给出一种 TVD型有限体积方法,主要思想是在一阶单调格式的基础上,在每一个单元上对变量作单调线性重构函数,时间离散采用二阶 Runge-Kutta方法。通过计算分析了该方法的精度,对平面激波反射和空穴流动的计算结果表明该方法是成功的。     

基于非结构网格的不可压N-S方程多矩有限体积法

提出了一种基于三角形及四面体非结构网格的有限体积法(FVM),用以鲁棒且精确地求解不可压粘性流动问题.与传统的FVM方法仅将体积分平均值(VIA)作为计算变量的做法不同,本文提出的方法将VIA及点值(PV)同时作为计算变量并在每个迭代步进行计算更新.VIA以通量形式进行计算以确保数值守恒,PV可以通过控制方程的不同形式进行求解更新,无需守恒,因此可以采用非常高效的方法进行求解.将PV作为增加的变量使得紧致网格模板得以实现更高阶精度的重构,而且由此获得的数值模型对于非结构网格变得更鲁棒.本文针对二维/三维的三角形/四面体非结构网格提出了数值格式,给出了几个基准测试算例,验证了本文提出的数值方法在采用非结构网格求解不可压粘性流动问题时的精确性和鲁棒性.     

分区和有限容积法计算长尾喷管中的湍流流动

为了给长尾喷管热防护结构设计提供热边界条件，从轴对称二维雷诺平均N－S方程出发，采用分区算法和有限容积法计算了长尾喷管中的湍流流动，分析了流场的结构特性，提出了长尾喷管中燃气温度，压强和马赫数分布的某些规律以及燃气温度沿壁面的分布，计算和分析结果提出了长尾喷管热防护结构设计应注意的事项。     

模拟飞机尾喷焰辐射传输的非结构有限体积法

有限体积法因其物理意义明确,在流动与传热的数值计算中有广泛的应用。在利用有限体积法求解辐射传输问题时,不仅能够保证在每一控制体内辐射能量守恒,还可以保证在每一个控制立体角内辐射能量守恒,因而使得模拟结果具有较高的精度;同时,非结构化网格技术是处理复杂几何形状物理问题的一种有效方法。本文编写了求解辐射传输的非结构化有限体积法程序,将计算结果与其他文献的算例进行了对比,证明了模型和程序的正确性,分析了各种因素对程序计算精度和计算时间的影响。在此基础上,将求解辐射传输的非结构化网格有限体积法程序集成到Fluent软件中,计算了飞机尾喷焰的温度场分布。     

斜网格修正的低速计算方法及应用

采用无坐标变换的一种有限体积法对RANS方程进行离散求解,离散方法考虑了对网格非正交情况的二阶修正,通过构造一个动量插值实现了同位网格下的压力修正算法,并发展了相应的分区算法.利用以上方法对轿车绕流和三峡坝区主流摆动问题进行了数值计算,计算结果与实验吻合较好.     

一种用于分离流动的网格自适应算法

对于可压缩粘性流动,提出利用流场速度的紊乱度作为指示变量进行网格自适应.Jameson中心格式有限体积法、五步Runge-Kutta时间推进法/双时间推进法求解定常/非定常N-S方程.基于雷诺平均N-S方程模拟紊流,选用SA一方程模型.在数值求解二维静态失速和动态失速问题过程中,加入网格自适应算法,提高数值模拟对流动分离特性的捕捉和分辨能力.算例结果表明在流场发生失速后,运用本文的自适应算法能够在增加少量网格单元的情况下明显提高计算精度.     

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

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

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

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

基于两相流理论滑动轴承动力特性求解

将计算流体动力学(CFD)两相流与动网格技术应用于滑动轴承动力特性数值求解,建立了基于CFD两相流滑动轴承动力特性求解模型,该模型无需设定油膜破裂边界条件且更能准确模拟滑动轴承流场特性.比较了单相流与两相流滑动轴承压力分布特性,计算分析了滑动轴承气穴分布特征及其影响因素,研究了两相流模型对滑动轴承动力特性的影响.计算结果表明:气化比例随着转速、偏心率和气化压力的增加而迅速增大,随进口压力的增加而缓慢减小.考虑两相流后,直接刚度系数增加,交叉刚度系数减小,直接与交叉阻尼系数均减小.随着偏心率的增加,单相流与两相流动力特性系数求解结果偏差增大.      

基于TLBM-FVM耦合的飞行器舱内热环境跨尺度预测方法

舱内热环境的有效预测是优化飞行器热控与防热设计、减小系统冗余并保障飞行器热安全的重要基础。在国家数值风洞（NNW）工程支持下,针对目前舱内热环境多尺度、多机制复合传热特点及其数值预测面临的精度与效率提升难题,发展了多区域协同推进的时空耦合模型及流/固界面的自适应分辨率识别算法,建立了基于热格子Boltzmann方法（TLBM）与有限体积法（FVM）相互耦合的舱内复合传热跨尺度预测方法,开展了典型飞行器仪器舱的综合热分析,验证了耦合方法的计算精度及效率。研究表明,相关方法可实现舱内热环境的局部精细化与整体大规模的协同模拟,用于开展整体自然对流与设备局部热量传递的多尺度数值模拟,掌握不同环境参数对舱内热质传递过程的影响规律,从而为飞行器热防护/热管理一体化设计提供重要技术支撑。     

分层嵌套重叠网格自适应树结构动态组装方法

采用重叠网格可以有效地进行复杂流动的大规模数值模拟,特别是包含运动部件(如旋翼、投弹)的动态模拟。本文将树结构的自适应直角网格用于重叠网格组装过程中的切割和贡献单元的搜索,大大加快重叠网格的组装速度。通过二叉树自适应直角网格对物体外形进行离散,实现切割过程的快速定位;采用八叉树自适应直角网格对流场区域进行离散,高效地搜索贡献单元。使用基于壁面距离准则的重叠区域最小化方法和分层嵌套重叠策略,能提高重叠网格组装的效率和质量。对于具有运动部件的动态重叠网格问题,采用多个二/八叉树减少组装过程中信息更新的冗余计算,从而大幅度减少重叠网格组装的时间消耗。实际算例的重叠网格组装结果说明本文发展的重叠网格组装方法具有很高的计算效率,可以满足运动边界复杂流动问题的动态计算要求。     

栅格翼导弹流场混合网格N-S方程数值计算

采用结构和非结构混合网格技术,对栅格翼导弹黏性流场进行数值模拟,预测其气动特性.计算的马赫数为0.7～2.5.用有限体积法和LU-SGS算法求解N-S方程.计算与实验数据以及栅格翼附近流场结构进行了比较.     

大型客机涡扇发动机动力特性模拟

机体/发动机干扰问题是现代大型客机设计中必须考虑和解决的核心问题之一。长期以来,商用和in-house计算流体力学（CFD）软件,都是通过在发动机进气口设置质量流量比,排气口和外涵道设置总温比/总压比,建立动力特性模型来等效模拟机体/发动机干扰流场。在现有动力特性模型基础上,借鉴特征边界思想,将外涵道指定为特征边界,建立了一种新的发动机动力特性模型。采用轴对称超高涵道比涡扇发动机模型、轴对称涡轮动力模型以及某型客机带发动机模型对两种动力特性模型进行了系统的验证和确认,结果表明：两种发动机动力特性模型均能很好地模拟涡扇发动机动力效应,且计算结果与试验数据吻合较好,说明了两种动力特性模型的正确性、可靠性以及工程适用性。此外,所建立的新动力特性模型特别适合于大涵道比涡扇发动机的动力效应模拟,且可以在外涵道总温比、总压比等参数未知的情形下,评估发动机的动力特性,在实际工程中应用更广泛。