LU—AUSM^＋混合格式在三维可压缩流动数值分析中的应用

将一维ＡＵＳＭ＾＋格式扩展到三维任意曲线坐标。为了获得高阶精度,采用了三阶ＭＵＳＣＩ格式。将ＡＵＳＭ＾＋格式与ＬＵ－ＳＧＳ格式结合,对凸包通道跨音速无粘流动、平面叶栅跨音速和超音速无粘流动,以及喷管超音速粘性流动进行了数值模拟。计算结果与文献计算结果和试验数据相符很好。     

叶栅激波流场的数值模拟

给出了ＡＵＳＭ 格式在任意曲线坐标下的形式,将ＡＵＳＭ 格式和三阶精度ＭＵＳＣＬ格式融合,并结合ＬＵ－ＳＧＳ隐式时间推进法,对典型双圆弧叶栅无粘跨音速和超音速流动进行了数值模拟．本文计算结果与文献计算结果相符     

AUSM+格式在亚音、跨音和超音速流动数值模拟中的应用

在任意曲线坐标系下,采用有限体积法计算可压缩流动。利用AUSM+格式替换传统的Rhie-Chow动量插值方法,给出有限体积界面的速度和压力求解方式。为了适用于亚音速流动,通过预处理的方法给出了速度-压力耦合的表达式。算例考核结果表明,该算法在不可压流动到超音速的全范围内都有良好的适用性,而AUSM+格式的计算结果均优于Rhie-Chow动量插值方法。     

非线性双曲型Euler方程的数值分析

采用 L U- AUSML W算法 ,对非线性双曲型 Euler方程进行求解 .该算法综合了 L U- SGS方法与 AUSMPW格式的优点 .为了验证 L U- AUSML W算法的有效性 ,在曲线坐标系中 ,对亚音速、跨音速和超音速流动进行了数值模拟 .数值试验的计算结果与文献计算结果相符很好     

组合的LU分解迎风方法的多重网格收敛

为提高计算收敛速度,以新的高收敛率的LU型隐式格式和高精度、高分辨率的MUSCL TVD迎风格式为基础,对可压缩无粘和粘性流动问题计算的多重网格方法进行了探讨。采用二维凸包通道和VKI LS 59跨音速透平叶栅为算例,通过与单网格上的计算迭代性能做对比,证明这种算法在原有基础上更大幅度地提高了计算收敛速度,节省了CPU时间,而且方法本身也简单易行。     

以压力为变量解可压流动的非交错网格SIMPLE方法研究

提出一处以压力为基本求解变量的非交错网格ＳＩＭＰＬＥ方法，用来计算可压缩流动，本方法的特点在于将适用球可压缩流动计算的ＳＩＭＰＬＥ方法加以扩展，使之也适用于可压缩流动的计算，因而有较大的适用范围，对亚音速、跨音速及超音速等３个较典型的算例进行了数值计算，计算结果与实验及文献的结果符合较好。     

一种高精度隐式计算格式的构造及其应用

AUSMPW格式虽然精度很高,但计算时间稍长.AUSM+up格式的计算速度较快,然而精度不及AUSMPW格式,而且具有较强的参数局限性.因此,结合AUSM+up格式和SD-SLAU格式等计算格式的优点提出了MAUSM+up格式.通过比较隐式格式AF方法和LU-SGS方法,选用计算速度更佳的AF方法与改进的MAUSM+up格式结合,对已有的三维黏性流计算程序加以改进,并在H-O-H型网格中计算三维黏性流场.结果表明编写的计算程序是一种较为理想的涡轮叶栅流场数值模拟程序.     

数值模拟入射斜激波／平板湍流边界层干扰流动

应用GAO-YONG可压缩湍流模型数值计算了入射斜激波／平板湍流边界层相互干扰现象。计算程序中的对流项、扩散项分别采用AUSM格式和中心差分格式离散,并用多步Runge-Kutta显式时间推进法求解空间离散后的控制方程。计算中包含了无分离流动、初始分离流动以及较大分离流动等多个情况,比较了平板壁面压力、法向平均速度剖面、壁面摩阻系数Cf以及壁面斯坦顿数St等的计算结果与实验值。结果发现：GAO-YONG可压缩湍流模型能够很好地预测入射斜激波／平板湍流边界层相互干扰下的无分离以及小分离流动,对高马赫数下的大分离流动也能得到较合理的结果,但再附点之后的壁面摩阻系数以及斯坦顿数的计算值不够理想。     

基函数法在黏性可压缩流动中的研究与应用

研究基函数法在二维和三维轴对称情形下可压缩黏性流动中的应用, 并利用此法数值地计算了二维圆柱的超音速黏性绕流和三维球头的轴对称超音速黏性绕流两个算例。取一阶三角函数为基函数, 构造出导数的中心格式和迎风格式。对于可压缩N-S 方程中的对流项, 采用通量分裂法及中心格式与迎风格式相结合的技术;对于可压缩N-S方程中的黏性项, 则采用中心格式进行处理。由此, 构造出了数值求解黏性可压缩流动的一阶三角函数类型的基函数格式。两个算例的数值计算结果表明, 基函数法不仅在处理无黏可压缩流动时, 是一种高精度、高分辨率的新型计算方法, 而且在处理黏性可压缩流动问题时也一样行之有效。     

基于高阶谱差分和嵌套网格方法的三维无粘流动数值模拟

数值方法采用高阶谱差分格式,单元边界上的通量采用基于Roe格式的近似黎曼解,在嵌套网格上求解三维Euler方程.嵌套边界上使用了一种高阶紧致的插值格式.对三维亚音速定常无粘流动进行数值模拟,使用嵌套网格与单块网格得到的计算结果吻合很好.     

通风分裂算法在有限元非结构网格中的推广

将在结构网格中发展起来的通风分裂算法推广到有限体积非结构网格中,并成功地进行了定常跨声速二维欧拉方程的数值计算．此格式具有流通量矢量分裂格式的简单性及通量差分分裂格式的高阶精度．最后给出的ＮＡＣＡ００１２翼型绕流的结果与实验结果符合的较好．     

预处理算法在非定常流动计算中的应用

介绍了一种可应用于全速度范围的非定常N-S方程数值求解方法,并全面考核该方法的适用性。该方法通过引进伪时间导数项,对伪时间导数进行预处理以解决低速时方程系数矩阵刚性过大的问题;同时用双时间步推进算法保证时间精度,将该方法推广应用到非定常问题的求解。文中给出了若干算例,覆盖了从极低马赫数到超音速、无粘/层流/湍流、二维/三维、定常/非定常情况。算例结果表明,该方法在较宽的使用范围内均能得到理想的计算结果。     

跨音速下旋成体干扰流场计算

应用粘性-无粘性干扰方法计算了亚跨音速下旋成体的轴对称绕流流场,并着重研究了船尾处的后体流场.将流场划分为无粘的位流区和粘性影响集中的边界层区,利用有限差分法和积分关系式法分别求解之后,以边界层位移厚度为粘性匹配条件进行耦合计算,求出弹身表面的压力分布.该方法可直接用于计算导弹燃气喷流产生的“阻塞”效应对船尾-后体流场的影响。     

新型高分辨率格式及其在CFD的应用

 利用小波奇异分析构造了一种新型的高分辨率高效计算格式,借助小波奇异分析所找到的流场中奇异点区域和光滑区域,合理选择高分辨率、高精度WENO格式和高精度中心差分格式进行流场的离散求解;利用该格式对超声速二维前台阶流动问题、跨声速RAE 2822典型翼型的二维外流场、跨声速VKI-LS59涡轮叶栅的二维流场以及NASA Rotor 37跨声速压气机转子和NASA Rotor 67跨声速风扇转子的三维流场进行数值模拟,计算结果表明,采用这种新型格式可以用于内外流复杂流场计算。数值计算还发现,在采用相同网格的情况下,本文格式不仅具有通常WENO格式高分辨率的特点,而且计算效率较WENO格式高出3～5倍。     

用时间推进法数值求解叶轮机械跨音叶栅的粘性绕流问题

应用一种高效，可靠的数据方法来模拟叶轮机叶栅内部的实际流动，获取大量真实的流动信息，据此，在工程上可以总结出各种实用的规律，这样，可以为现代叶型设计提供一种快速，可靠，实用的途径，数值试验表明：用文中提出的二阶精度隐式矢量单步差分格式求解叶栅的跨音流场，具有精度高，稳定性好，计算量少，收敛快等优点，在流场的计算中，较好地捕获了激波，与实验比较，结果令人满意。     

TVD有限体积法在喷流计算中的应用

该文根据ＷａｎｇＪ．Ｃ．Ｔ．的二阶ＴＶＤ有限体积法求解守恒型Ｅｕｌｅｒ方程，对原文中的数值格式进行了部分修改。计算网格单元是非均匀的任意四边形。文中计算了楔形风洞中的斜激波反射问题和无粘喷流近场流动结构问题。数值计算结果表明，该文方法具有激波分辨率高、计算稳定、伪振荡小等优点。同时该文方法可以处理任意复杂边界的二维及轴对称流动问题的计算，也可以将该文的有限体积法思想推广到三维具有复杂边界的流动的计     

二维无旋超音速扩散对称流动控制方程及其特征线方程和相容性方程

通过改变止推气体轴承等间隙的流道结构,对变气膜厚度圆盘止推气体轴承中二维无旋超音速扩散对称流动控制方程进行推导,得出其特征线方程和相容性方程,对所得特征线方程和相容性方程用等间隙气膜内一维超音速流动的解析解进行验算.结果表明,当内点间距合适时,特征线法数值求解与解析解结果十分接近,推导的特征线方程和相容性方程正确.