求解可压缩N－S方程的无粘－粘性分数步法

用按物理过程作时间分裂的无粘－粘性分数步法求解可压缩Ｎ－Ｓ方程，无粘步用高效欧拉算法解欧拉方程，粘性步用全隐格式解抛物方程；并采取了有效措施消除由分步法产生的时间和空间不相容性，文中给出了轴对称喷管内外流场的数值模拟结果。     

粘性跨音速喷管流场计算

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

二元喷管无粘流的有限体积法分析

本介绍了一种计算收敛－扩散喷管内流特性的有限体积法，它将求解域分成许多有限体积单元，将欧拉方程在小体积单元上积分，从而得到半离散方程，然后采用Mac Cormack格式进行时间推进，求得定常解。计算结果与试验值吻合良好，是一种有效的计算具有喉道小曲率半径不连续的收敛－扩散喷管特性的方法。     

二元矢量喷管流场的N－S方程计算

对二元矢量喷管在不同工况下的流场进行了数值模拟。求解的控制方程为三维非定常Ｎ－Ｓ方程。空间方向用有限体积离散，时间方向用Ｒｕｎｇｅ－Ｋｕｔｔａ法推进。粘性项的计算上采用了一种有限体积框架下新的离散化方法。湍流模型采用了广泛使用的Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数模型。计算表明，本文方法能准确地计算二元矢量喷管、非矢量喷管在设计与大设计落压比下的流场。     

二维非结构网格上的高阶间断有限元方法研究

对二维非结构网格上的求解欧拉方程的高阶间断有限元方法进行了研究.运用间断有限元基本理论,采用多项式基函数对解进行近似描述,欧拉方程的数值通量项运用高斯积分公式和Roe格式求取,时间推进采用传统四步Runge-Kutta方法.运用发展的间断有限元方法对NACA0012翼型在亚跨音速情况下的无粘流动和圆柱的低速无粘流动进行了数值计算,结果表明:间断有限元方法具有良好的收敛特性、很小的数值耗散和很好的激波捕捉能力.     

突然起动圆柱的非定常绕流计算

本文用分区方法定量模拟高Re数情况下突然起动圆柱的初期对称流动。绕圆柱的二维不可压流场可分为无粘区、附着粘性区和分离区。无粘区有解析解，用非定常层流边界层方程解附着粘性区，用离散涡方法处理随时间变化的分离区。试用了混合的Rankine和Lamb涡模型及二次离散涡。给出了Re=300、9500两个算例，基本符合有关实验结果。     

三维带操纵面机翼欧拉方程数值模拟

本文给出了用欧拉方程求解带操纵面机翼的数值结果，采用有限体积法，四步Ｒｏｎｇｅ－Ｋｕｔｔａ时间推进，结合不连续面通量守恒传递技术，对有操纵面偏转的三维机翼流场进行求解，计算结果与实验符合较好。     

二元喷管跨音速流场的数值计算

本采用有限体积法进行空间离散，Runge-Kutta法进行时间推进的办法，通过求解Euler方程来模拟二元收敛－扩散喷管跨音速流场。运用隐式残差平均技术和当地时间步长法加速计算收敛。为了消除气流参数的波动和激波前后的振荡，在方程中添加了自适应耗散项。在100个时间步以内便可获得稳态解，计算结果与试验数据和其他数值方法吻合良好。     

三维跨声速粘性-无粘相互作用的半逆方法

本文为计算三维跨声速粘性-无粘相互作用提供一半逆方法版本。湍流边界层采用积分反方法。证明了对于跨声速流,选取外部流线角α和等价的不可压缩形参作为输入量,边界层积分方程总是双曲型的。无粘外流采用FL027程序。通过粘性-无粘迭代求得小分离区情况下的与实验符合的收敛解。     

直升机旋翼悬停流场的欧拉方程计算

描述了三维Ｅｕｌｅｒ方程数值模拟旋翼悬停流场的过程,欧拉方程的求解采用了风格中心有限体积法、五步Ｒｕｎｇｅ－Ｋｕｔｔａ时间推进格式。     

跨声速翼型有粘／无粘强干扰计算

本文介绍一种二元跨声速激波-边界层强干扰的计算方法。边界层计算采用湍流边界层积分反方法,它借助Whitfield和Swafford提供的既适合附着流,也适合分离流的速度剖面表达式。跨声速无粘流用全速势方程模拟。通过边界上排溢速度来考虑粘性的影响,用有粘/无粘迭代得到粘性流解。本方法计算的结果与其它方法以及实验的结果进行了比较,证明该方法可以在工程上推广使用。     

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

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

计算矢量喷管流场的有限元法

本文通过采用扫描有限元法及多步龙格－库塔时间推进格式求解非定常平均Ｎ－Ｓ方程来模拟二元矢量喷管的内部跨速流场,成功地计算了二维有激波的跨声速喷管内部流场及仰 仰角３０°的三维矢量喷管内部流场。     

火箭发动机塞式喷管流场的数值研究

用时间相关法和张涵信的无振荡、无自由参数ＮＮＤ显式格式,对火箭发动机塞式喷管二维轴对称定常流场进行了数值模拟。计算得到了有粘和无粘中两种情况下,喷管流场的马赫数压强、温度、流线以及速度矢量分布图,计算网络采用代数方法生成,并通过求解Ｌａｐａｃｅ方程对网格进行修正。计算结果表明：ＮＮＤ格式用于计算塞式喷管流场具有较高的效率和良好的分辨率;无粘和有粘情况得到的流场结构有较大的差别;基础喷管的出口马赫数     

不同马赫数的无粘和粘性流动高阶精度隐式计算方法

 应用隐式时间推进法对不同马赫数的无粘和粘性流动进行数值分析,给出了基于预处理方法的高阶精度隐式求解方法。方程离散采用改进的高阶精度对流迎风分裂格式。该方法通过求解曲线坐标系可压缩Euler和 Navier-Stokes方程,对低速到超音速范围内的无粘和粘性流动的典型问题进行了数值分析。计算结果与文献的数值结果或实验数据对比分析表明,该方法对低速到超音速范围内的无粘和粘性流动问题进行数值分析是可行而有效的。     

三维可压缩流欧拉方程组的隐式时间推进法

本文提出了求解三维定常无粘可压缩流场问题的一种隐式时间推进法。在计算空间以逆变速度分量为未知变量的欧拉方程组为控制方程,从而简化壁面边界的处理。在Beam-Warming的AF-方法基础上,用对角化及矢通量分裂等方法提高了计算速度及解的精度和稳定性。作为算例,对伴有激波的方形收放喷管跨声速流动进行了计算,以显示方法的有效性。     

高效欧拉方程解法求解轴对称引射喷管三股流问题

推导了三维流动的广义欧拉方程黎曼增量梯度形式的欧拉方程及其轴对称形式表达，针对引射喷管三股流问题的特点在运动网格下作了H变换。在此基础上将高效欧拉方程解法应用于求解轴对称引射喷管三股流问题，为解决次流通道狭窄带来的问题，采用了弱耦合法处理滑流边界。     

跨声速机翼流场的N—S方程计算

本文给出了三维机翼粘性流场的数值模拟过程，控制方程为时均化的三维可压缩薄层Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，湍流模型采用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ两层代数湍流模型，空间离散采用中心有限体积格式，时间采用Ｒｕｎｇｅ－Ｋｕｔｔａ多齿格式进行式积分。     

时间相关法计算喷管跨音速流动

本文用时间相关法完成了定常、无粘跨音速喷管流场的计算,内部点参数用Mac Cormack差分格式计算、壁边界点参数用简化了的特征边界条件计算,轴上点参数用反射原理计算。计算表明:计算是收敛的,计算结果与实验数据符合良好。 定常、无粘、跨音速喷管流动的控制方程是混合型的(即在亚音速区域是椭圆型,而在超音速区是双曲型的),给数值求解带来很大的困难。为了克服这一困难,广泛采用时间相关法,即认为定常流动方程的解是相应的非定常流方程的解在时间趋于无穷时的渐近解。因为非定常流的控制方程不论在亚音速区,还是在超音速区域,都是双曲型的,可以用统一的方法来求解,而且易于求解;另一方面,它可适用于形状比较复杂的喷管。 本文首先阐述了时间相关法的计算方法,然后列出算例喷管的计算结果。通过与实验数据的比较,证明计算是符合实际的。     

任意翼型有小分离气泡跨音速绕流的迭代算法

 本文给出一种带小分离气泡的任意翼型粘性跨音速绕流的计算方法,采用有粘-无粘干扰迭代的概念。无粘流的全速势方程用AF差分格式在保角变换法生成的计算网格中求解,粘流附面层方程用C-S盒式法求解,用逆算法消除分离点处的奇性。本文对Ma_∞=0.8,Re_∞=2×10~6,迎角α=3.5°和4°的NACA64A010翼型粘性绕流进行了计算,结果与实验相比较,吻合良好。