非结构网格下含冷却孔的涡轮转子三维流场计算

在转动坐标系里,讨论非结构网格下,含多排冷却孔的涡轮转子三维流场雷诺平均Navier-Stokes方程组的求解。计算表明:含冷却孔的转子流场计算要比静子流场复杂的多。为了得到较为满意的流场数值解,本文对叶片压力面与吸力面采用了不同的湍流模式、并且还对流场迭代中初始场选取方式以及转子叶尖间隙计算站的设置等问题进行了一系列的改进。几个典型流场算例表明:初始场选取、湍流模式的合理选用对转动坐标系下获得一个收敛的三维流场数值解至关重要;另外,转子叶尖间隙站的设置直接影响着间隙附近的流场,尤其对现代高负荷涡轮转子,间隙流动的研究就更为重要。      

三级涡流器环形燃烧室化学反应流场的数值研究

在三维曲线坐标系下采用多区域耦合法对包括三级涡流器在内的环形燃烧室三维两相反应整体流场进行了数值模拟。计算中采用区域法和偏微分方程法生成三维贴体网格，采用标准k双方程紊流模型，EBU-Ar-rhenius紊流燃烧模型和六通量热辐射模型。在非交错网格体系下，气相用SIMPLE法求解，液相采用颗粒群轨道模型，并用PSIC算法对其进行数值求解。计算结果表明，程序编制可靠，建立的三维网格生成和流场计算程序可为燃烧室优化设计和研制提供有用的数据。     

基于非结构重叠网格的二维N—S方程求解与应用研究

本文发展了非结构重叠网格技术.通过构造任意单元类型网格重叠的网格间边界确定算法、宿主单元搜索算法以及网格间插值算法,把非结构重叠网格推广应用于求解NS方程,并结合湍流模型模拟二维湍流流场.N-S方程的求解采用有限体积法,通量计算采用Osher格式,湍流模型采用SA模型.为验证本文方法的正确性和有效性,模拟了RAE2822翼型和某三段翼型绕流湍流流场,结果表明采用非结构重叠网格模拟复杂外形湍流流场不但网格生成效率高,而且容易推广应用于具有相对运动的非定常问题,具有较高的工程应用价值.     

雷诺应力模型在三维湍流流场计算中的应用

从雷诺应力模型出发,通过求解雷诺平均N-S方程组获得三维湍流流场的数值解。计算中比较了多种湍流模式,并进行了相应的流场计算。本文完成了两个典型算例。从算例1的计算结果与实验值比较中发现:采用雷诺应力模型(RSM)计算的三维流场比采用k-ε模型更贴近实验值;算例2采用了RSM模型及三维非结构网格,对一典型内流问题进行了三维流场计算。算例的数值实践表明:采用雷诺应力模型可以有效的计算各向异性的湍流流场;另外,发展非结构网格有助于模拟壁面附近的流动,并节省计算机内存。      

三维多块结构椭圆型贴体网格的生成

通过求解椭圆型方程生成三维贴体网格，采用分块技术生成多块结构网格体系。通过实验验证了方法对任意复杂网络区域生成三维多块结构网格能满足流场计算的要求，其质量是相当满意的。     

三维Euler方程的自适应非结构网络计算

介绍一种三维自适应非结构网格的Euler方程求解技术。采用有限体积中心差分的格点格式对Euler方程进行空间离散,多步龙格-库塔时间推进,并采用当地时间步长、焓阻尼修正等加速收敛。采用Delaunay三角剖分技术生成三维流场的四面体网格。网格自适应技术中的网格局部加密也采用Delaunay三角剖分技术。最后用ONERAM6机翼的跨音速绕流计算说明本文方法的有效性。     

带隔热屏加力燃烧室热态流场计算

在任意曲线坐标系下对带有横向波纹隔热屏、外冷却通道和尾喷口的涡喷加力燃烧室的三维热态流场进行数值模拟.利用椭圆型微分方程和区域法生成三维贴体网格, k-ε模型预估紊流特性, EBU-Arre-henius紊流燃烧模型估算化学反应速率, 六通量模型预估辐射通量以考虑其对壁面和气流温度的影响.在非交错网格系统下应用混合差分格式离散控制方程, SIMPLE算法求解离散方程, 数值分析两种进口气流参数分布对加力室热态流场和壁温分布的影响.所得计算值与试验数据基本相符, 表明计算方法合理.      

涡轮叶栅冷气掺混数值模拟方法

采用具有三阶精度TVD性质的有限差分格式、自由型曲面网格生成技术以及分区网格算法,对某型涡轮叶栅进行了任意数目冷气孔的薄膜冷却全三维N-S方程数值求解,较为详细的阐述了冷气掺混流场中卵型涡的发展规律。结果表明,本文所开发的求解系统能够较好的实现任意分区的流场求解,对涡轮叶栅中冷气喷射流场的描述是合理的。      

基于Van Leer＋AUSM混合格式求解三维N-S方程

基于多块粘性结构网格,开展了三维N-S方程数值算法的研究。控制方程的空间离散采用有限体积法,在前人工作的基础上,发展了Van Leer＋AUSM混合格式并应用到对流通量的离散中,粘性项采用中心格式离散并利用格林定理计算粘性通量中的导数项,时间推进采用五步R-K法,湍流模型为S-A一方程模型。最后,以M6机翼和某超声速弹丸的粘性流场作为数值算例,计算表明：发展的数值算法对跨声速、超声速流场均具有较高的分辨率,适用于跨、超声速流场的数值模拟。     

基于非结构网格的多段翼型流场数值模拟

将非结构网格生成技术应用于多段翼型流场数值模拟，针对多连通域可以快捷、迅速地生成非结构网格。通过控制背景网格节点的网格生成参数，进行网格的局部加密，提高多段翼型流场分析的精度，应用有限体积法求解欧拉方程，应用各种加速收敛措施提高流场求解速度，并将数值计算结果与风洞试验结果进行了对比分析，证明该方法具有型号设计应用价值。     

矩形/非结构混合网格技术及在二维/三维复杂无粘流场数值模拟中的应用

建立了一套模拟复杂无粘流场的矩形／非结构混合网格技术,其中非结构网格仅限于物面附近,发挥非结构网格的几何灵活性,以少量的网格模拟复杂外形;同时在以外的区域采用矩形结构网格,发挥矩形网格计算简单快速的优势,有效地克服全非结构网格计算方法需要较大内存量和较长ＣＰＵ时间的不足。混合网格系统由修正的四分树（２Ｄ）／八分树（３Ｄ）方法生成。本文将ＮＮＤ有限差分与ＮＮＤ有限体积格式有机地融合于混合网格计算,消除了全矩形网格模拟曲面边界的台阶效应,同时保证了网格间的通量守恒。该算法被成功地运用于ＮＡＣＡ００１２翼型的跨、低超声速绕流,二维弯管内运动激波的非定常流,以及三维捆绑火箭的超声速绕流的数值模拟。计算结果表明本方法在模拟复杂无粘流场方面的灵活性和高效性。     

基于非结构网格有限差分法的扎染算法

应用有限差分法时遵循贴体坐标系下离散等价方程的离散准则,格式中仅用到当地网格点的坐标变换系数,据此可以构建新的非结构网格有限差分法,数值算例表明在空间二维情况下用连接离散点的3条网格线构造出来的一阶迎风格式可以稳定收敛。把这种非结构网格有限差分法推广到三维,在4条网格线基础上进行离散计算,然后提出一种局部区域多重网格、重复计算的新算法。首先在包含物体的整个区域采用直角坐标系下的均匀网格进行计算,然后消除物体内部点对物体外部流场区域的影响,最后在物面和直角网格之间很小的局部区域填补非结构网格进行计算,计算过程类似于中国传统的扎染工艺。处理流场内存在曲面物体边界问题时,相较于常规结构网格差分算法,扎染算法尽管在计算中包含需要特殊处理的无用网格点,也存在重复计算的局部区域,但是以这些计算能力的浪费为代价,换取了编程简便、内存量小、网格生成快捷、易于扩展网格规模等优势。定性分析表明,这种扎染算法非常适合研制大规模并行计算,建立的计算流体力学应用软件可以充分发挥数十万至百万处理器核心的超级计算机的效能。     

一种改进的流场计算加速收敛方法

探讨了二维非结构网格无粘流场计算的加速方法。在对残差平均和多重网格这两种主要加速方法进行分析后,提出了一种改进的兼具多重网格思想的残差平均方法。经对绕NACA0012翼型无粘流动计算的验证,收效明显。在几乎不增加计算量的情况下,比原有的残差平均方法可减少迭代步25%～30%.该思想可推广用于三维、有粘流动计算和结构化网格中。     

基于伴随方程的网格自适应及误差修正

基于流场方程的离散伴随优化理论和三维非结构网格,建立了网格自适应技术和目标函数误差修正方法。详细研究了用流动伴随变量进行目标函数的误差估计和修正技术,构造了适用于格心格式有限体积法的流场变量插值技术和网格单元剖分判据,初步实现了网格物面投影和空间单元优化,发展了适用于有限体积法的整套网格自适应方法。对NACA0012翼型和ONERA-M6机翼绕流进行了自适应数值模拟,并对升、阻力等目标函数进行了误差修正。数值结果表明,本文自适应方法能正确地捕捉到影响目标函数计算精度的敏感区域,网格自适应和误差修正两项技术显著提高了升、阻力等目标函数的计算精度。     

一种用于捕捉间断多维任意单元的有限元法

本文根据差分格式构造过程中对三阶、四阶人工粘性及熵增条件的分析，构造出适合于非结构单元的有限元法。计算实例表明，这种方法能较好地捕捉间断，激波前后的伪振荡基本消除。     

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

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

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

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

基于非结构网格的气动弹性数值方法研究

基于非结构网格和网格变形技术,研究了耦合求解流体和弹性结构体互相作用问题的数值方法.气动弹性模型采用了耦合求解可压缩雷诺平均的Navier-Stokes方程和线性结构动力学方程的方法.为了能处理在复杂几何域内的问题,采用了基于混合单元的非结构网格有限体积方法.开发的方法在二维翼型简谐振动进行了验证并应用于AGARD 445.6翼型颤振边界的数值模拟.