用高效矢通量分裂格式求解二维叶栅的跨音流动

提出了一种高效隐式矢通量分裂格式。该格式通过在显格式上添加合适的高阶小量, 避免了近似因子分解和系数矩阵运算, 大大地减少了计算量。同时该格式在收敛性及激波捕捉精度等方面均保持了矢通量分裂格式所具有的优势。将该格式用到了二维叶栅的无粘跨音流场的数值模拟中去, 结果令人满意, 证实了无矩阵运算、无近似因子分解的隐格式构造有一定的推广意义。      

平面大头叶栅跨声速流动的一种有效解法

准确数值预估叶片头部表面的压力分布,一直是叶片气膜冷却设计中十分重要而又尚未圆满解决的问题。本文针对这个问题提出了一个高准确度、高效率的计算方法。 在文献[1]的基础上,本文从任意正交曲线坐标系中的基本方程出发,引入Von Misses变换,导出了相应的流线控制方程,并具体地在拟边界层坐标系中计算了RA叶栅流场,结果与实验符合甚好,据作者所知,本文方法是已见到的计算同类叶栅方法中准确度最高、花费机时最少的。     

计算跨声速平面叶栅流场的压力方法

以压力和速度分量为原始变量,将不可压流的ＳＩＭＰＬＥ算法应用到可压缩流动,成功求解了具有严重非正交网格及周期性边界条件的跨声速平面叶栅流场,计算结果与实验符合较好。     

三维振荡叶栅内部非定常流动数值模拟研究

进一步发展了非定常雷诺平均N-S方程求解程序,并基于影响系数法对三种不同振型下的三维涡轮振荡叶片绕流问题进行了数值模拟研究.数值模拟结果表明,所发展的程序对振荡叶栅流动模拟具有较好的精度;振型对振荡叶栅内部非定常流动以及叶片表面的非定常气动力有着较为重要的影响,其中叶片运动导致叶栅各通道中流量的再分配以及出口气流角的变化等对叶片的负荷分布会产生较为显著的影响.     

带有机匣处理结构的叶栅流场数值计算

完成了一种数值计算带有机匣处理装置叶栅流动的方法。计算采用三维、非定常、可压缩Ｎ－Ｓ方程组及ｋ－ε紊流模型方程，选择３种典型的处理结构，周向槽式机匣、轴向缝式机匣、轴向倾斜缝式机匣为算例，并引入参数周向平均的概念，以与稳态实验测量的结果做定量对比，结果发现符合较好。本文为分析机匣处理的扩稳机理提供了理论基础，也为工程设计机匣处理装置提供了新的数值工具。     

两个涡轮叶栅跨声速流场的数值模拟

采用CFD软件CFXTASCflow时两个涡轮叶栅跨声速流场进行了分析，其中一个为高压涡轮导向器环形叶栅，另一个为高压涡轮动叶中截面平面叶栅，计算结果与试验数据吻合良好。     

求解跨声速压气机叶栅粘性流动反问题的数值解

阐述了一种以有限体积的时间推进方法为基础求解跨声速压气机叶栅设计的反问题方法。应用这一技术设计叶型,规定叶片表面无量纲目标速度分布,通过比较目标速度分布和计算获得的速度分布来修改叶片压力面和吸力面的坐标,最终获得要求的叶片形状。计算程序使用分布体力方法模拟粘性,用局部时间步长和多重网格方法加速收敛。还介绍了这种方法求解的基本方程系统,并给出了采用这种方法设计的无激波超临界跨声速压气机叶栅的设计结果。     

平面压气机叶栅通道分离流动高精度直接数值模拟

采用高精度有限差分格式直接求解二维Navier-Stokes方程组,数值模拟V103平面压气机叶栅分离流动,数值结果表明:在瞬时流场中,吸力面后部发生流动分离,在分离区前端存在大尺度分离涡,分离涡下游是由二次涡和脱落涡交替形成的涡串,直至叶片尾缘,形成以脱落涡为主结构的尾迹;在时均流场中,吸力面后部存在短分离泡,分离区压力分布存在明显压力平台。与逆压力梯度下平板边界层分离流动相比,瞬时和时均流场结构相似;叶栅通道内无量纲涡脱落频率是前者的两倍。与文献计算结果对比表明:叶片表面压力分布除分离区外吻合很好;非定常计算所得分离区轴向长度比定常计算大41%。在分离区内三个二阶统计量均达到最大值,表明流场强非定常性集中在分离区。     

吸附式压气机平面叶栅数值模拟研究

本文对采用吸附式压气机概念设计的一套超声速平面叶栅采用粘性全三维数值模拟方法进行了研究。数值模拟再现了吹风试验中发现的一些现象．并进一步研究了抽吸位置、抽吸气流量对叶栅性能的影响。     

涡轮叶栅叶型损失的数值模拟方法

给出了一个计算亚、跨音涡轮叶栅叶型损失的数值计算方法。主流采用时间推进有限体积法求解积分型欧拉方程,并采用了局部网格修正技术;附面层采用全隐格式求解有限差分方程;在叶栅出口与远后方均匀流之间进行了叶片尾迹与主流的掺混损失计算。算例表明本文的数值方法可准确地预测涡轮叶栅的叶型损失。     

一种新型多维通量函数在求解平面叶栅跨音速流动的应用

采用一种改进的Ｒｏｅ格式对平面叶栅的无粘流动进行了数值模拟,并与Ｒｏｅ格式的计算结果进行了比较。采用四阶Ｒｕｎｇｅ－Ｋｕｔａ时间差分格式,空间三阶有限体积法。网格生成采用贴体网格近壁面修正的方法。计算结果表明,该格式在减少格式对网格的依赖性,抑制数值粘性方面具有明显的作用。     

透平跨音速叶栅的优化设计

将遗传算法用于透平跨音速叶栅的优化设计,提出了相应的数学模型和优化方法。优化算例表明,本文所提出的模型及算法能有效的优化叶型方案,避免了一般优化算法的局部最优或维数灾难问题,从而在新的领域内将优化计算在叶轮机械设计中的应用进行了有益的探索。      

喷管二维跨声速两相湍流流场的数值模拟

采用ＳＩＭＰＬＥ算法,结合分散颗粒轨道（ＰＳＩＣ）模型,在同位网格上实现了喷管中跨声速两相湍流流场的数值模拟。根据Ｒｈｉｅ,Ｃｈｏｗ提出的动量插值概念,以协变物理速度分量作为求解变量,完成跨声速气相流场的数值计算。湍流模式采用ｋ－ε双方程模型,计算结果与提供的实验数据相符合。     

固体火箭喷管两相粘性跨音速流场计算

本文进行了固体火箭喷管两相粘性湍流跨音速流场计算.粘性的气相控制方程用隐式近似因子分解法求解,粒子方程采用跟踪粒子轨迹的特征线法求解,粘性湍流选用代数模型,气相和凝相充分地偶合.CFL数可以取500左右,收敛速度快,使粘性两相跨音速喷管流场计算耗费的机时达到工程计算可以接受的程度.通过计算,获得了在粘性和粒子同时作用下的流场参数分布.这对固体火箭发动机流场需要同时考虑粒子和粘性作用的专题研究大有裨益.     

跨音速喷管流场波系分布的数值仿真

应用分离系数矩阵方法对跨音速喷管流场进行数值仿真, 同时计算喷管内马赫波分布。采用5151个网格点的计算结果清晰反映出喷管内膨胀波和压缩波的生成和反射。数值仿真与实验结果在宽广马赫数范围内吻合一致, 并揭示了波系分区分布, 反映出波系影响下跨音速喷管流场参数分布的规律性。      

用N-S方程计算三维叶栅流

发展了基于微元中心有限体积方法,采用雷诺平均N-S方程求解三维叶栅流动的程序,特别注意了格式中人工粘性的加法;程序采用Baldwin-Lomax模型。格式经过平板层流及湍流流动的校验,层流与Blasius解、湍流与实验经验公式及Wilcox的计算值符合很好,最后以Goldman涡轮导向器为例进行了叶栅流动计算校核,结果符合于实验以及Chima等人的计算。      

跨音速低展弦比轴流风扇转子内部流场的数值研究

采用三维雷诺平均Ｎ－Ｓ方程程序求解跨音速压气机风扇转子内部流场。这个程序采用有限体积显式时间推进方法求解控制方程,方程自变量选取在控制体顶点,采用简单的Ｈ型网格离散方程,离散后的代数方程用两步Ｒｕｎｇｅ－Ｋｕｔａ法进行时间推进求解,使用多重网格和当地时间步长方法加速计算收敛。程序中使用了Ｂａｌｄｗｉｎ－Ｌｏｍａｘ模型模拟湍流影响。用该方法计算了ＮＡＳＡ６７＃跨音速风扇转子的性能曲线,并重点计算最高效率点工况下内部流场。计算结果与实验结果的比较,证明该法可以用来模拟跨音速风扇转子内部流动。     

跨音叶栅流动的多级有限元分析

将各向异性张力有限单元和Taylor-Galerkin有限元法在多级有限元法的思想下结合起来, 构成了收敛速度和稳定性均较好的适用于跨音流动计算的各向异性多级有限元法。利用这一方法, 基于Euler方程, 研究了平面跨音速叶栅流动, 并与实验结果作了比较。结果表明, 本方法是跨音速计算的强有力手段。      

直接的矩阵分裂计算跨音速喷管流场

给出了分离系数矩阵求解欧拉方程组时系数矩阵的直接分离公式，减少和简化了应用该方法时的矩阵运算，参照Ｂｅａｍ－Ｗａｒｍｉｎｇ显式格式构造与ＳＣＭ方法一致的内点差分格式，并分析了它的数值特性，对有斜激波的跨音速内流场进行了数值仿真。计算中采用变系数的当地时间步长，加快了流场达到的收敛速度，计算无需引入任何人工参数，计算值与实验结果吻合。     

考虑顶隙的压气机单转子内三维跨音速粘性流场结构的数值模拟

借助于考虑叶片顶部间隙的压气机转子内跨音速湍流流场的三维数值解,得到了该转子流场中的复杂涡系与波系结构。给出了叶片表面及轮毂面上的壁面流型,分析并揭示了转子叶道内的马蹄涡、通道涡、角涡、尾涡等复杂涡系及其与激波的相互干涉;特别是,分析并揭示了转子叶道内顶隙射流的特点,得到了顶隙附近激波与泄漏流动的干涉形态及泄漏流动的发展,并着重给出了其与激波的干涉形态。