任意回转流面上跨音速压气机叶栅含激波流场的计算

本文运用时间分裂有限体积法,进行了任意回转流面上跨音速压气机叶栅含激波流场的计算。分析了MacCormack显式两步差分格式中的几种不同运算形式,提出了对于叶栅绕流计算问题的合适的运算方案。对几种边界条件的处理方法进行了数值试验和分析比较。流场计算中采用了一种周向非等距网格,这将有助于缩短计算时间。     

粘性跨音速叶栅绕流的数值模拟及前驻点处理

一、前言 对于粘性跨音速叶栅绕流,直接用N-S方程数值求解仍有很多困难,因此利用边界层理论考虑流体的粘性效应,进行有粘—无粘迭代计算是粘性流体叶栅绕流的很好模拟。采用多重网格有限体积法进行跨音速叶栅的无粘流计算,具有易于适应复杂几何形状,保证差分格式的守恒性、能较准确地计算激波。加速收敛等良好特性。为了提高精度和简化计算,在求解边界层微分方程时,采用了Illingworth变换方法,用驻点方程的相似性解给定边值速度剖面,使用全湍流Prandtl混合长度粘性模式并且其混合长度选用Plctcher     

叶栅激波流场的数值模拟与高分辨率差分格式

一、前言 近年来,激波流场数值模拟技术的一个重要进展,是各种高分辨率总变差减小(TVD)格式的建立和发展。这类格式与常用的中心差分方法相比,其突出的优点在于,不需要人为地在差分方程中附加耗散项就能给出激波上、下游无波动的数值解,保证了数值流场具有很高的空间分辨率。目前这种方法已引起人们的高度重视,并逐渐发展成为应用计算中的流行格式。本文的目的,是探讨TVD格式在叶栅激波流场数值分析     

跨音速叶栅流场计算的多网格法

本文是对多网格法应用于求解Euler方程的探讨,文中以MacCormack显式两步差分格式为基础,提出了一种简单的求解积分型Euler方程的多网格算法,阐述了这种方法的基本思路,研究了稳定条件和一种新的人工粘性形式.文中还进行了几种跨音速涡轮叶栅流场的计算,用以说明多网格法在缩短计算时间上的效果。     

用显格式时间推进法求解跨音叶栅的粘性流动

本文提出了一种基本格式加修正环节的显式差分格式,用来求解跨音叶栅绕流的完全的Navier-Stokes方程,可得到近似于二阶精度的解.试算表明本文差分格式和给出的边界条件是成功的,它提供了透平机械叶栅中粘性流动的一种数值模拟方法.在给定了合适的湍流模型后,本文的方法能方便地推广应用到湍流的计算中去.     

应用Navier-Stokes方程的叶栅粘性流动数值分析

本文以时间相关显式格式求解完全的Navier-Stokes方程,计算了二维叶栅中的粘性流动。一种Euler方程多网格算法,被成功地应用于粘流问题计算,显著地加速了收敛过程。用含分离流动的压气机叶栅和跨音速涡轮叶栅等算例,说明了解的准确程度和多网格法的效果。     

双相各向同性介质BISQ模型弹性波方程数值模拟方法

本文从Taylor级数展开式出发,推导出了交错网格一阶空间导数的任意偶数阶精度展开式和相应差分系数计算式;从本构方程和运动方程推导出了BISQ模型双相介质一阶双曲型应力-速度弹性波方程交错网格任意偶数阶精度差分格式以及推导出二维双相各向同性介质完全吸收层边界条件公式和相应的高阶交错网格差分格式。通过数值模拟研究表明,该方法边界吸收效果好,稳定性好,能够高精度模拟双相介质中地震弹性波场,且计算效率也高。     

两种TVD格式在跨音透平叶栅流场计算中的应用

两种ＴＶＤ格式在跨音透平叶栅流场计算中的应用黄伟光，刘建军（中国科学院工程热物理研究所北京１０００８０）关键词跨音速透平，Ｅｕｌｅｒ方程，ＴＶＤ格式１引言现代高负荷透平大都工作在跨音范围，研究发展能准确预测跨音速透平叶栅内部流场的激波位置与强度以及叶...     

任意迴转面叶栅超声速进口流场特征线法计算机程序和唯一进气角的确定

根据吴仲华提出的基本方程,采用特征线法,对任意迴转面叶栅超声速进口流场的计算编制了计算机程序.计算中特征线方程中的系数是采用特征线两端的平均值,因而需要进行迭代.为了提高计算精度,并提出了一种最近三点二元插值法.叶栅前缘脱体激波的形状和位置是采用Moeckel提出的近似方法确定的.为了满足叶栅流动周期性条件,在计算中对叶栅各个通道进行逐个计算,直到相邻通道中流动情况相同为止.根据守恒条件,确定了无限远来流状态.最后讨论了根据这些计算确定唯一进气角的方法,并给出了计算实例.     

使用非正交曲线坐标与速度分量的S_1流面跨声速流函数解

在非正交曲线坐标系上采用守恒型流函数方程,用人工可压缩性方法解任意(辶回)转面上叶栅跨音流动,考虑了通过激波间断的熵增。用速度梯度方程积分的方法解决了密度双值问题。并试验比较了几种数值计算方法。     

基于Voronoi网格的扩散方程差分格式

在Voronoi网格上利用一种基于回路积分法的有限体积法构造扩散方程的的差分格式.在这种特殊的网格上离散扩散方程比通常在四边形网格上离散的格式要简单,不会引进角点未知量,提高了对网格边上的流的离散精度,及差分格式整体精度.这种Voronoi网格上的扩散计算也可以与单元中心流体力学计算耦合.数值算例表明这种格式比四边形网格上的格式精度高,且能更好的应对网格扭曲情形.     

跨声速叶栅流的激波捕获-分区计算法

本文基于文献[1]提出的建议,即先用流函数方程或势函数方程计算压气机叶栅的跨声速流场,得到大致的通道激波位置后,再对激波的上、下游区分别进行计算;最后通过对激波位置的调整以满足Rankine-Hugoniot条件,得出确切和明晰的激波形状及气流参量通过激波的突跃变化.文中对具有实验数据的一个双圆弧叶栅分别用势函数方法和流函数方法捕获通道激波并将二者所得激波的平均位置作为分区计算时进行通道激波调整的初始波形.在计算结果同实验值的比较中,还考虑了平面跨声速叶栅实验时实际存在的轴向速度密度比和沿流线熵增对计算结果的影响,所得计算结果是接近实验值的。     

气冷涡轮级叶栅非定常流场数值模拟

采用具有三阶精度TVD性质的有限差分格式、自由型曲面网格技术、分区算法以及双时间步长的方法,对某型涡轮级叶栅流场进行了非定常NS方程数值求解,考察了在有、无冷气喷射条件下涡轮级气动性能的非定常变化。结果表明,上游静叶栅是否喷射冷气对下游动叶栅超音速区域的影响具有明显的区别,有冷气喷射时动叶栅前缘气动负荷降低,级效率下降约1%,但是不同动叶通道内气动性能随时间周期性变化的幅度明显减小了。     

叶栅内定常及非定常粘性流动数值模拟

本文给出了一个模拟叶栅内准三维定常和非定常粘性流动的数值方法。对于定常流动,采用ＴＶＤ Ｌａｘ－Ｗｅｎｄｒｏｆｆ格式和代数湍流模型求解雷诺平均Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程,使用当地时间步长和多网格技术使计算加速收敛到定常状态;对于非定常流动,使用双时间步长和全隐式离散,采用与求解定常流动相似的多网格方法求解隐式离散方程。文中给出了ＶＫＩ透平叶栅内的定常流结果和１．５级透平叶栅内的非定常数值结果。     

环形燃烧室火焰简及旋流器流场计算

本文采用偏微分方程法生成贴体网格,在任意曲线坐标系下数值研究两种先进燃烧室火焰筒及其旋流器三维紊流流场。由于旋流器的形状复杂,本文采用型线定点法确定网格的边界。在非交错网格系下采用SIMPLE算法和混合差分格式对离散方程进行求解。计算结果表明计算方法合理,这计算程序进一步扩展,可用来预估环形燃烧室反应流流场。     

一维非线性对流占优扩散方程的变网格特征差分方法

针对一维非线性对流占优扩散方程,提出了一类变网格特征差分格式,该格式能够根据解的梯度变化及时对计算网格进行调整.与均匀网格格式相比,给出的变网格特征差分格式对于对流占优扩散问题有着更好的计算效果.     

解三维扩散方程的积分内插法格式

研究三维扩散方程的数值模拟.在非正规六面体网格上,使用积分内插法建立扩散方程差分格式,涉及到27个相邻网格,适用于大变形网格上带间断系数的拟线性扩散方程的计算.叙述差分格式的建立,推导通量流和网格顶点温度的计算公式,给出了数值试验结果.     

涡轮叶栅双排孔气膜冷却数值模拟

采用具有三阶精度TVD性质的有限差分格式、自由型曲面技术以及分区网格算法,对某型具有冷气孔形状的涡轮叶栅进行了全三维N-S方程数值求解,描述了相邻两排冷气射流在叶栅吸力面形成的冷却气膜以及壁面附近冷却射流运动的特点,分析了不同吹风比和喷射角度情况下冷却绝热效率的分布规律。结果表明,在较大的吹风比和喷射角下,交错排列的两排冷却射流运动规律非常复杂,在两排孔之间的区域与冷气孔下游区域冷却气膜的形成规律具有明显的区别。     

气相爆轰高阶中心差分-WENO组合格式自适应网格方法

研究一种高阶中心差分-WENO组合格式,并采用自适应网格方法进行二维和三维气相爆轰波的数值模拟.采用ZND爆轰模型的控制方程为包含化学反应源项的Euler方程组.组合格式在大梯度区采用WENO格式捕捉间断,在光滑区采用高阶中心差分格式提高计算效率.采用一种基于流场结构特征的自适应网格.计算结果,表明这种方法同时具有高精度、高分辨率和高效率的特点.     

平面叶栅气膜冷却流动的数值模拟

为了能够准确地对透平叶栅气膜冷却效率进行数值预测,本文采用了FNM形式的结构化网格,对一个平面叶栅中的气膜冷却流场进行了数值模拟。计算中采用了包括LU-SGS-GE隐式格式和改良型高精度、高分辨率的MUSCL TVD格式的时间推进算法求解三维RANS方程以及低Reynolds数q-ω双方程湍流模型。计算结果表明本文采用的模型及方法在低吹风比的条件下可以较准确地对气膜冷却效率进行数值预测。