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

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

NURBS曲线及节点插入算法在透平叶片优化设计中的应用

将NURBS曲线及节点插入算法应用到透平叶栅的二维参数化造型,可以兼顾造型参数的实际物理意义和型线调节的灵活性;将控制点可调的NURBS曲线应用到三维叶片积叠线造型,能够构造任意形状的弯、扭、掠叶片。在此基础上搭建的透平叶片三维气动优化平台,将参数化造型、网格自动划分、CFD并行求解集成在一起,应用实验设计和响应面分析方法进行多目标寻优。对某单级实验透平分别进行单排叶片和整级优化后,单排叶片总压损失降低了0.57%;单级绝热效率提高了0.8%。     

吸力面附面层抽吸在三维高负荷扩压叶栅中的作用机制

对某矩形高负荷扩压叶栅在不同弦向位置开设全叶高抽吸槽的5组方案进行了数值研究,分析了抽吸槽弦向位置等参数对抽吸量分布规律的影响;通过叶栅实验探究了局部展向抽吸方案的效果.数值仿真的计算域包含吸附叶片内部的真空腔,边界条件按照实验条件设置.研究发现:全叶高抽吸方案的抽吸量沿展向大致呈C型分布;叶高中部和端部的主要抽吸效果都体现在叶高中部流场,端部的抽吸量对叶栅角区的回流有一定的抑制效果.抽吸量沿展向的分布规律受叶栅流道和叶片内腔流场的共同作用,因此应根据三维高负荷扩压叶栅流场的具体特性对吸力面抽吸槽/孔进行细化设计.      

环形弯扭静叶栅流场的计算与实验研究

本文应用S2流面正问题计算方法和三维Euler方程的时间推进计算方法对某涡轮弯扭静叶栅流场分别进行了分析计算, 并在环形涡轮叶栅实验台上测试了该叶栅出口截面上气动参数的分布。通过分析、比较计算结果和实验数据, 既考核了算法, 也研究了弯扭静叶栅中的流动机理, 由此得出了对弯扭叶片的气动设计计算具有重要意义的几点结论。      

压气机叶栅多点气动优化设计

采用Pareto排名策略和共享小生境技术,实现遗传算法的多点优化设计。压气机叶栅的多点气动优化设计是针对某初始叶栅同时优化不同来流攻角状态点的气动性能,使其达到攻角特性最优。压气机叶栅网格由NU-MECA/Igg软件生成,流场由NUMECA/Fine软件求解。多点气动优化设计结果显示出本文发展的优化设计平台能够较好的实现多点组合气动优化设计。     

非均匀有理B样条曲线及节点插入算法在透平叶片优化设计中的应用

将非均匀有理B样条(NURBS)曲线及节点插入算法应用到透平叶栅的二维参数化造型,可以兼顾造型参数的实际物理意义和型线调节的灵活性;将控制点可调的NURBS曲线应用到三维叶片积叠线造型,能够构造任意形状的弯、扭、掠叶片.在此基础上搭建的透平叶片三维气动优化平台,将参数化造型、网格自动划分、流体计算并行求解集成在一起,应用实验设计和响应面分析方法进行多目标寻优.对某单级实验透平分别进行单排叶片和整级优化后,单排叶片总压损失降低了0.57%;单级绝热效率提高了0.8%.      

多块协调变形的网格变形技术及其应用

基于多块网格协调变形的理论,提出了适用于叶轮机内流多物边界的更具一般性的网格变形技术,先获得块边界处的变形量,然后插值得到块内网格各点的变形量。通过二维振荡叶栅和带叶尖间隙的三维跨音风扇叶片的非定常绕流算例,证明网格变形技术在叶轮机气弹稳定性数值模拟中的应用是有效的,能够在保证网格高精度的条件下,提高网格变形速度。     

叶栅三维流场的实验研究

本文通过流动显示和详细测量,研究了叶栅中的三维流动。在平面叶栅风洞中,对非平直端壁的大弯度压气机叶栅所完成的测量,包括沿叶高的叶型表面上气流参数分布和叶栅通道中的整个流场。为了进行对比和详细研究端壁效应,还对二维叶栅完成了同样的测量内容。对所观察到的复杂流动现象进行了分析。实验结果还与叶栅全三维无粘流计算进行了对比。本文所获得的完整的三维流场实验数据,为建立叶栅流动模型,研究小展弦比叶栅中气流展向掺混的机理,以及为检验全三维流动计算方法提供了可靠的依据。      

亚音速复合倾斜透平静叶栅的理论与实验研究

采用三维全势函数多重网格数值方法求解亚音速透平静叶栅流场。以 30 0 MW汽轮机调节级静叶栅的设计为目标 ,采用子午面上外环壁收缩和叶片周向复合倾斜两项措施 ,以求大幅度地减少叶栅二次流损失。实验研究表明新设计叶栅的损失系数比直叶栅降低 42 % ,测出的叶片表面静压分布与理论计算结果十分吻合     

弯-掠叶片对压气机叶栅端壁流动的控制作用

利用数值和实验方法对气动弯掠在压气机叶栅中的作用机理进行了研究。结果表明,相对于网格数目和网格类型,进、出口边界条件和湍流模型对计算和实验结果吻合的影响要大。叶片前掠改变了叶栅的三维压力场,使得低能流体重新分配。正弯与前掠的结合进一步加强了这种控制作用,为叶片设计提供了一个自由度。     

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

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

通道中有稳定器的湍流流场分区计算方法

采用任意曲线坐标系对含有Ｖ型稳定器的通道内湍流流场进行数值计算。对含有Ｖ型槽的复杂流场中网格分区方法作了详细的分析,提出一个新的网格划分方法。数值计算采用曲线坐标系下非交错网格的ＳＩＭＰＬＥ方法,计算时对整个流场进行分区迭代,直到得到收敛结果。文中对一含有Ｖ型稳定器的直通道湍流流场进行了数值计算,结果合理。     

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

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

基于非结构动网格的非定常气动力计算

对用于非结构动网格生成的弹簧近似方法进行了研究。通过采用顶点弹簧方法,分析研究了弹簧倔强系数的取值,同时通过引入挤压倔强系数和边界修正,对标准弹簧近似方法进行了改进。改进后的方法可以大大提高网格变形能力和网格质量。应用本文发展的非结构动网格生成方法耦合求解基于Arbitrary Lagrangian-Eulerian(ALE)描述的欧拉方程,计算了谐和振动NACA0012翼型及矩形机翼的跨声速非定常气动力,计算结果与参考文献提供的结果及实验结果吻合良好。     

非结构网格生成方法的改进及气膜冷却三维静子流场的求解

讨论五面体网格下,涡轮静子叶片含多排冷却孔时三维粘性流场雷诺平均Navier-Stokes方程组的数值求解问题。考虑到叶片含多排冷却孔的工程背景,本文发展了生成五面体非结构网格的快速有效方法即仅在一个拟S₁面上生成非结构网格,而沿叶高方向仍采用结构网格的生成办法。文中还采用了新的网格优化措施以改善生成网格的质量。流场的计算采用了非结构网格下的高分辨率迎风格式。几个典型算例表明:本文发展的网格生成方法快速有效,而且网格生成质量较高;发展的计算格式能够获得较为满意的三维流场数值解,并具有较高的激波分辨率。      

三维自适应非结构网格的Euler方程解

 将 Ausm+ 迎风格式应用于三维非结构网格中求解 Euler方程。对单元变量进行重构以获得空间高阶精度,对时间域采用多步龙格库塔法推进,并采用了当地时间步长和隐式残差光顺技术来加速收敛。采用多点择优推进阵面法生成复杂曲面的三角形网格,利用推进阵面法生成四面体网格。采用网格自适应技术对网格进行局部加密,以减少总体网格数目,从而提高计算效率。最后给出了绕 ONERA M6机翼的跨音速流动及绕麻雀 导弹的超音速流动算例,结果表明了本方法的有效性。     

复杂外形飞行器多块结构贴体网格生成

通过求解椭圆型方程生成三维贴体坐标网格,采用分块技术生成多块结构网格体系。每个子区域转换成计算空间内的立体域,具有自身的曲线坐标系,从而所形成的整个网格系统不会过分扭曲和不规则,这将有利于对于复杂边界的真实飞行器绕流流场的求解。文中讨论了控制函数对网格生成的影响。     

二维结构/非结构混合网格的生成及其流场模拟

采用分区求解方法在二维结构/非结构混合网格中对粘性绕流作了数值模拟。利用层面推进法生成绕物体的O形结构网格,利用推进阵面法在结构网格外生成非结构网格。在结构网格区内,采用B-L湍流模型,求解雷诺平均N-S方程。在非结构网格区内求解Euler方程。两区的解通过在交界面上满足通量守恒条件进行耦合。数值算例表明了方法的有效性和实用性。     

结冰计算中翼面网格重构研究

对结冰数值计算中的网格生成进行了深入研究,针对结冰后几何边界的复杂性及边界依赖于时间步长变化对网格生成的影响展开了分析.网格生成是结冰计算的基础,文中基于求解椭圆形偏微分方程(PDEs)生成了未结冰时翼面周围的背景网格,采用Hilgenstock源项构造方法,合理控制了结冰边界周围的网格间距和正交性.模型表面结冰后,为...     

A 3D hybrid grid generation technique and a multigrid/parallel algorithm based on anisotropic agglomeration approach

A hybrid grid generation technique and a multigrid/parallel algorithm are presented in this paper for turbulence flow simulations over three-dimensional (3D) complex geometries. The hybrid grid generation technique is based on an agglomeration method of anisotropic tetrahedrons. Firstly, the complex computational domain is covered by pure tetrahedral grids, in which anisotropic tetrahedrons are adopted to discrete the boundary layer and isotropic tetrahedrons in the outer field. Then, the anisotropic tetrahedrons in the boundary layer are agglomerated to generate prismatic grids. The agglomeration method can improve the grid quality in boundary layer and reduce the grid quantity to enhance the numerical accuracy and efficiency. In order to accelerate the convergence history, a multigrid/parallel algorithm is developed also based on anisotropic agglomeration approach. The numerical results demonstrate the excellent accelerating capability of this multigrid method.