有限分析方法求解粘性可压缩紊流

本文将有限分析方法进一步推广到粘性可压缩紊流流场的数值研究, 在贴体座标网格上数值分析了单列和串列叶栅内粘性可压缩紊流流场。计算表明, 有限分析方法具有数值精度高, 稳定性好和收敛速度快等优点, 适于粘性可压缩紊流流场数值计算。      

数值模拟二维跨音速叶栅可压粘性流动

用有限体积法数值模拟了任意回转流面跨音速压气机叶栅中的粘性流动，四步龙格－库塔法用于N－S方程的时间推进。给出了一种抑制叶栅计算中“数值失速”的方法，并采用隐式残差光顺技术加速收敛。用所发展的方法计算了三种叶栅在宽广来流Ma数条件下的流动，计算结果与实验结果吻合良好，流场激波具有较高的分辨率。     

弦向缝隙叶栅对边界层分离的控制

本文提出了一种在Ｓ１流面上进行势流与粘性流迭代求解弦向缝隙叶栅流场的计算方法，将Ｓ１流面上势流方程组的解与经Ｉ１１ｉｎｇｗｏｒｔｈ－Ｓｔｅｗａｒｔｓｏｎ变换后的边界层方程结合起来联立求解弦向缝隙叶栅流场。计算值同实验结果相当吻合，在正大攻角时，向缝隙叶栅能效地控制边界层分离。     

LU-AUSMPW混合格式在可压缩无粘和粘性流动数值分析中的应用

本文首先将 AUSMPW格式与三阶 MUSCL格式融合, 给出了其在任意曲线坐标下的三维形式, 并将其与 LU-SGS格式结合, 应用于可压缩 Euler和 Navier-Stokes方程的求解。其次, 构造了一种新的通量限制器。最后, 为了验证 LU -AUSMPW混合格式的性能, 将平面叶栅跨音速无粘流动以及喷管超音速粘性流动作为算例。本文计算结果与文献计算结果和实验数据相符很好, 表明采用 LU-AUSMPW混合格式数值模拟可压缩流场, 具有较高的计算精度、较快的收敛速度和良好的稳定性。      

梯度加权局部分析解差分格式及其应用

提出了一种新的针对Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程符合含有梯度项的对流扩散方程的梯度加权局部分析解的有限差分格式，且给出了通用的系数表达式。ＧＷ－ＬＡＤ格式是基于积分型的有限体积法并且通过附加源项的方法来实现。能够适用于粘性和无粘性及不可压缩和可压缩的压力驱动流动。通过对叶栅流道中二维湍流流的计算表明ＧＷＬＡＤ格式与ＥＤ格式相比，具有较高的精度和很快的收敛速度并且得到更好的压力场和速度场。     

用空间推进法来解钝锥三维粘性绕流

本文将文献［１］中提出的方法推广应用于三维可压缩Ｎ－Ｓ方程，文中给出了钝锥超声速三维粘性绕流的结果。得到了背风面分离的流动图画，计算结果与实验相吻合。     

用空间推进法求解钝锥三维粘性绕流

本文将文献〔１〕中提出的方法推广应用于三维可压缩Ｎ－Ｓ方程，方中给出了钝锥超声速三维粘性绕流的结果。得到了背风面分离的流动图画，计算结果与实验相吻合。     

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

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

叶轮机械失速颤振的全三维数值分析方法

给出了建立在全三维N-S方程求解基础上的叶轮机械失速颤振分析的数值方法。包括用定常N-S方程求解三维粘性失速流和用三维N-S振幅方程求解三维粘性失速振荡流。为了验证方法的准确性, 用压气机环形叶栅的失速流实验和振荡叶片振荡流场的实验作了对比, 同时应用本方法分析了某30万千瓦汽轮机末级叶片的失速颤振。      

湍流附面层方程的有限差分解

本文用三点隐式有限差分法数值求解了轴对称流动的层流、转捩和湍流可压缩附面层方程,湍流附面层用两层概念处理,而每一层用适当的涡粘性模型代替Reynolds应力项,用常值湍流Prandtl数将湍流热流项与Reynolds应力关联起来。支配方程用迭代法求解。计算结果表明,采用本文改进的涡粘性模型和涡粘性分布的平滑办法,对文中所计算的几种情况,都给出了满意的结果。     

扫描有限元法求解二元喷管的三维可压粘性流场

发展了一种二元喷管流场计算的三维Ｎ－Ｓ方程扫描有限元方法,给出了Ｎ－Ｓ方程中二阶导数项新的离散方法,以及更完善、简洁的扫描有限元方程的表达式。开发了相应的计算程序,用所发展的方法对二维平面收扩喷管、两种大宽高比的二元喷管及其圆变方转接段的三维粘性流场进行了计算,计算结果表明本方法具有较宽广的适用范围,计算结果良好     

钝体前激波的数值模拟

本文通过直接求解N-S方程,研究激波问题,对是否可以和如何用直接求解N-S方程来处理超跨声速流进行了探索,在对粘性项的处理建立可压缩流体的压力校正方程和计算方法等方面均提出并实施了新的方法,最终得到了基本合理的数值计算结果。     

求解叶轮机粘性流场的再生解析数值方法

应用再生核方法求解叶轮机粘性流场是一种新颖的解析数值方法, 它可以解决计算中的存贮、计算速度和精度问题。本文给出了用再生核方法求解叶轮机粘性流场的思路和方法, 以及如何在曲线坐标系中应用再生核方法, 并给出了流函数积分方程的求解过程。最后对平面扩压管道进行了数值试验, 结果与有关文献是吻合的。      

可压缩流动SIMPLE算法压力修正方程边界条件处理探析

可压缩粘性流体Nivier-Stokes方程SIMPLE算法求解中压力修正方程起着关键的作用.提出从连续性方程本质出发确定压力修正方程边界条件的方法.通过对双圆弧(凸包)通道内流场的计算验证了方法的正确性.     

一种模拟喷管二维流场的新方法

采用ＳＩＭＰＬＥ方法对固体火箭发动机喷管进行了纯气相流场的数值模拟。推导了建立在同位网格之上利用协变物理速度分量为计算变量的离散方程，采用协变物理速度分量推导压力修正方程，对密度采用一阶迎风格式插值，将ＳＩＭＰＬＥ方法扩大到计算可压流动。计算结果表明：本方法计算精度高，收敛速度快，程序编写简单。     

燃气排导箱三维场的数值计算

对燃气排导箱内的典型三维稳态流动过程进行了粘性不可压数值计算,采用较为成熟的ＳＩＭＰＬＥ方法。所发展的计算程序可求解质量、动量和能量守恒方程,ｋ－ε双方程湍流模型方程,以及热辐射通量方程,并进行了程序的算例验证。三维燃气排导箱内的计算结果也得到了一些实验测量的数据的验证,得出一些有工程实用价值的结果和结论。     

粘性跨音速喷管流场计算

本文继用隐式近似因子分解法成功地计算无粘跨音速喷管流场之后,用同样的方法,结合任意曲线坐标系,通过非定常方程在相当长时间后的定常解,针对矩形截面喷管和轴对称喷管粘性跨音速流场,求解了可压缩层流薄层N-S方程,获得了核心流,特别是边界层中的流动参数.对于矩形截面和轴对称两种喷管进行了计算,其结果和实验数据相当一致.将本文的方法应用到粘性两相流动计算,可望较多的节省机时.     

基于黏性伴随方法的跨声速机翼气动优化设计

 进行了基于黏性伴随方法和Navier Stokes方程的跨声速机翼气动优化设计研究。分别推导了适用于三维跨声速机翼气动反设计和减阻设计的黏性伴随方程、边界条件和梯度求解表达式,并研究了伴随方程的数值求解方法。通过将网格生成、流场计算、黏性伴随方程数值求解、梯度求解和拟牛顿优化算法等几方面的有效结合,发展了一种跨声速机翼气动优化设计方法。为了提高计算效率,将多重网格方法应用到方程的数值求解中来加速收敛。跨声速机翼反设计和减阻设计算例验证了本文所发展的方法的正确性。采用本文的方法进行优化设计,一般通过20~30次迭代就能得到满意的结果。     

复杂区域内粘性流动的数值计算

本文应用原始变量法(SIMPLEC)计算了复杂区域内的不可压和可压缩流动问题。计算时,仍使用交错网格,但在求解动量方程时采用协变分量做为原始变量,有效地消除了在任意非正交坐标系下压力和笛卡尔速度分量不耦合而产生的压力波动。在压力修正方程中考虑了密度的影响并构造出一种密度的迎风格式,保证了计算的稳定与收敛,使之预示马赫数变化范围较大的流动成为可能。对由NACA65(12A10)10构成的平面叶栅、由NGTE10C4/30C50构成的平面叶栅的不可压缩流动和平面超声速缩放喷管进行了计算。结果表明,本方法有较好的收敛性和准确性,能够适用于从亚音速到超音速的流动场合,有一定实用价值。      

低速叶型气动反问题设计方法

 低马赫数不可压流动中声速与流速大小差别巨大,采用基于可压缩流动控制方程的计算格式求解流场时,由于数值黏性的污染,解的精度低且收敛性差,通常可使用时间预处理技术来解决这一问题。在基于控制理论的优化方法中,共轭方程的Jacobian矩阵和流动方程的系数矩阵相似,因而在低流动马赫数下,求解共轭方程存在着与求解流动方程相同的数值污染和数值刚性问题。首先推导了带有预处理的Roe格式,然后发展了适合全速度流动的共轭方程求解方法,最后选取翼型和叶栅两个典型算例进行了验证。计算结果表明所发展的方法可很好地用于低马赫数时的气动反问题设计。