平面跨音叶栅中粘流-无粘流的强相互干扰计算

为了设计高效率高负荷的跨音叶轮机械,气动计算中必须计入粘性的影响。当叶片表面流动不发生分离,或分离区很小时,流道内的流动可以分为叶片表面以及尾缘下游的粘性薄层——边界层和尾迹区以及除此以外的无粘流区。因而可以分别应用欧拉方程组及边界层方程组求解这两个区内的流动。无粘流-粘流的相互影响是通过无粘流计算及边界层计算的一系列迭代计算来实现的:由给定的气动参数及叶栅几何形状进行无粘计算,得到叶片及尾迹表面的流动参数分布,作为首次边界层及尾迹计算的输入数据,求得边界层及尾迹的排挤厚度δ_1,并以此对无粘流动的几何边界进行修正,作为下一次无粘计算的流场边界。这样的迭代重复进行,直到满足全场的总体收敛条件为止,其迭代框图如图所1示。     

叶轮机械中三维跨音流场的数值模拟

由于航空发动机和大功率凝汽式蒸汽透平研制的实际需要,近年来对于叶轮机械中带激波的三维跨音流场的研究、分析与计算十分重视,并已取得了相当大的进展。这些工作都是基于时间推进法来进行数值求解。[5]由微分型的基本方程组出发,采用MacCormack格式进行计算,得到了压气机叶列中的三维跨音流场。该方法的不足之处是,要求的运算时间太长,这可能是由于使用微分方程,为保证所需精度在计算区域中安排了多达25000个网格点的缘故。[1—4]及[6、7]从积分型的基本方程组出发,     

高分辨率TVD格式的应用

在Jameson格式的基础上，利用修正通量法，构筑TVD（Total Variation Diminishing)格式。通过三个典型的一维和二维算例的计算结果表明，TVD格式对激波具有较高的分辨率，它为具有复杂激波系的喷管、叶栅等跨音速流场计算提供了一种有力的工具。     

关于跨音速叶栅物面边界条件的数值试验

本文讨论了跨音速叶栅物面边界的数值处理方法,目的是寻求与我们所用内点格式相匹配的边界上流动参数的算法,以提高计算精度。本文首先推导了适合于物面边界计算的特征关系式,其次总结了五种边界处理方法,进行了数值试验。计算结果表明,配合我们的内点计算格式,以y向抛物外推的Denton修正法计算效果较好,能够较精确地反映跨音速叶栅带激波流场的物理特点。     

网格生成软件GRGES的设计与实现

首先,本文说明了研制网格生成软件对于促进计算流体力学研究和应用的重要意义。其次,本文论述了数值网格生成的基本概念和基本步骤,提出了网格生成软件GRGES的设计思想,简单介绍了GRGES的总体结构和五种网格生成方法。最后给出了三个网格生成算例,它们均可用于计算实际流场。     

跨音涡轮中三维流动的数值试验

 本文从积分型的基本方程组出发,以我们的关于平面叶栅和任意旋成面叶栅的工作为基础,将文献（3、4）中的计算方案推广到三维流动情况。编制了相应的计算机程序,进行了一些数值试验,得到了初步的计算结果。 数值试验表明,拟三维计算结果可以导致相当大的误差,从而说明了完全三维流动计算的必要性。另一方面,如果取2000个左右的计算点,则在国产的013机上算一个例子,约需100分钟。从实际应用的观点来看也是可行的。     

隐式通量分裂线松驰迭代法求解平面叶栅无粘绕流问题

本文工作是在 Mac Cormack通量分裂格式 [2 ] 基础上发展的一种有限面积通量分裂隐式格式 ,其特点是在求解隐式离散化方程时 ,采用往返扫描一次的 Gauss- Seidel线松驰迭代方法 ,避免了对时间步长的任何限制。为提高定常解精度 ,格式的显式右端项采用二阶精度的离散。在数值求解跨音速涡轮平面叶栅问题中 ,对壁面边界作了较仔细的隐式处理。数值计算表明本方法保持了 Mac Cormack格式具有的高收敛速率 ,(约 30个时间步即可达到定常解 )而每一时间步计算量约减少一半。数值结果与实验结果符合得很好。