跨音速压气机级三维紊流流场数值模拟

在对跨音速压气机级的动／静双叶排中的三维定常紊流流场数值的模拟中，利用三阶高分辨率ＮＮＤ格式和ＬＵ－ＳＧＳ隐式推进迭代法，既保证了流场中激波的模拟质量和粘性流动特征的正确预估，也实现了求解过程的高效率。对某单级跨音速压气机的计算及其与实验数据的对比证实了本文方法的有效性。将计算推广至多级轴流叶轮机械的情况将会是十分简便的     

一种分离流数值求解方法

本文以弯管内分离流场为例，根据分离区与未分离的湍流动能差别较大的特点，提出了在分离区和未分离区采用不同的湍流模型来模拟流动的湍流效应，利用Ｎ－Ｓ方程进行了数值模拟，结果表明：利用这一方法，可以加速迭代的收敛，节省机时。     

高速粘性内流的高分辨率高精度迎风型杂交格式

本文在有限体积离散和ＬＵ分解的基础上，构造出一个新的隐式迎风型杂交格式，并用于求解定常流动的稳态解。     

分层Reynolds应力模式

本文根据切变湍流结构的特征提出一种分层Reynolds应力模式。在切变湍流中存在湍能生成项达到狭峰的极大值层次(本文称它为活跃层),根据这一特征,在活跃层中保留微分Reynolds应力模式,但它可以简化为常微分方程。在两层活跃层之间以及活跃层和壁面之间的Reynolds应力分布,用满足衔接相容条件的多项式近似。时均湍流场将由以上得到的Reynolds应力代入Reynolds方程求出。用本方法计算了两个比较复杂的湍流例子,结果表明,本方法较常用的k-ε模式有明显的优越性,它的计算时间并不比k-ε模式多,而计算结果与实测结果的吻合则较k-ε模式好得多。另一方面本方法较之完全的微分Reynolds应力模式可大大节省计算时间。     

绕大攻角翼型低速非定常流动的SUPG有限元数值模拟

采用既有良好稳定性又有较高精度的ＳＵＰＧ有限元方法,通过数值求解涡量一流函数形式的不可压缩流ＮＳ方程组,模拟大攻角翼型在高Ｒｅ数下几种典型情况的非定常绕流流场,同时验证了ＳＵＰＧ有限元方法在求解对流占优的二维非定常对流的－扩散方程时的有铲性。     

用双流体—轨道模型模拟四角喷燃模型炉内三维湍流两相流动和煤粉燃烧

目前炉内两相流动和煤粉燃烧数值模拟中多半用颗粒随机轨道模型和单流体无滑移模型，这些模型都难以完整地给出三维空间内颗粒速度，浓度，湍流度分布的信息。主采用双流体－轨道模型（颗粒相连续介质－轨道模型）对一个四角喷燃模型炉内三维湍流两相流动及煤粉燃烧进行了模拟。此模型基于欧拉气相方程组、欧拉颗粒连续方程组和动量方程组以及拉氏颗粒能量方程和质量变化的方程，并使用k-ε-kp两相湍流模型，EBU－Arrhenius湍流燃烧模型，离散坐标辐射传热模型，煤粉颗粒的水分蒸发，热解挥发模型和焦炭燃烧的扩散－动力模型等。热态模拟中，为了减小为信散造成的影响，采用了扭转坐标法（将坐标扭转一定的角度使之与煤粉射流方程一致）。为了检验数值模拟，采用三维相位移普勒测速仪（PDPA）对于冷态模型炉内湍流两相流场进行了测量，得到了两相速度，湍流脉动及颗粒浓度的分布。分别对冷态模炉内两相流动和热态模型炉内三维两相流动和煤粉进行了模拟，冷态两相流动的计算与实验结果的对比表明预报的两相流场是合理的，热态模拟的结果给两相速度，气相温度、组分浓度及壁面热，显示出靠近出口处气相速度和温度分布不对称，造成一个局部高温区。     

直接加热式换热器中若干问题的探讨

目前广泛使用的直接加热式换热器,存在着热效率低、振动和噪声大的问题。本文通过对汽泡运动与热交换的理论分析和数值计算,以及对蒸汽激励与防噪声的分析,提出既可提高热效率又可大大降低振动与噪声的结构方案。实验证明了这种结构方案的有效性及其在工程中应用的可行性。     

用 Euler 方程计算某些特殊的二维分离流动

采用非定常Ｅｕｌｅｒ方程、总体时间步长、隐式时间推进，计算了几种特殊类型钝体的二维有分离绕流。一类是光滑的圆柱体，在跨音速时计算出因激波诱导分离而周期性地脱落在尾迹上的非对称涡排。另一类是带尖角的物体，如三角形柱体，计算出从尖角处脱落在尾迹上的分离旋涡，形成卡门涡街。对这两类物体，计算的平均阻力系数，分离涡脱落的Ｓｔｒｏｕｈａｌ数以及流态都与实验接近。结果表明，对于这类绕流问题，用非定常Ｅｕｌｅｒ方程计算，具有一定工程实用意义。     

跨流域三维复杂绕流问题的气体运动论并行计算

通过研究求解描述跨流域三维绕流问题的Boltzmann模型方程气体运动论耦合迭代数值格式,分析气体运动论数值计算方法的内在并行性;从变量依赖关系、数据通信与并行可扩展性三方面开展基于离散速度空间区域分解计算的研究,发展求解稀薄流到连续流跨流域三维复杂绕流问题的并行算法。通过对不同Knudsen数、不同马赫数、不同攻角三维球体及返回舱绕流的并行计算,计算结果与实验数据和理论分析吻合较好。研究表明,该并行算法负载平衡和并行可扩展性较好,对不同并行计算机系统具有很好适应性,显示直接求解分子速度分布函数的气体运动论计算方法有良好的并行计算特性。     

乘波飞行器气动力、热特性的数值模拟研究

采用多块分区网格以及并行计算技术对给定乘波构形的高超声速飞行器进行了数值模 拟,分析了飞行器前缘小半径钝化对飞行器气动性能的影响,计算了前缘钝化后飞行器表面 的热流分布状况。结果表明,前缘钝化对飞行器的升力影响不大,对阻力和升阻比的影响较 大。对于曲率半径为1cm的钝化前缘,与原尖前缘飞行器相比,其升力降低了0.78％,阻 力增加6.96％,升阻比下降7.21％。前缘钝化后,乘波飞行器仍具有较好的气动性能, 飞行器前机身可为发动机提供比较均匀的气流,飞行器整体仍基本保持了乘波的状态,热流 主要集中在飞行器前缘上。为了有效防热,需要采用主动式冷却技术。