定常轴对称超音速喷管和超音速射流的计算

本文叙述轴对称超音速喷管和超音速射流的数值计算。内点用MacCormack有限差分格式,分步加粘性;边界用从特征关系列出的二阶精度格式;喷管口用Prandtl-Meyer流公式;射流边界用一个迭代过程确定。对有激波喷管和无激波喷管做了大量的数值实验,并作了比较和讨论,给出了三个喷管的计算结果。     

斜切反喷管性能计算

本文提出了斜切反喷管流场计算方法,内部点用MacCormack二步显格式数值解控制方程,入口边界和固壁边界的参数用物理边界条件和可应用的沿双特征线的相容性方程计算.固壁边界计算中采用了激波装配技术和开关格式,有效地捕获了流场内的激波系列,计算中得到的激波情况和物理分析一致.本文在流场计算的基础上,进行了反喷管的性能计算.本文提出的方法可应用于反喷管的性能预估.     

气动计算中一个简单有效的差分格式

本文采用算子附加修正方法,提出了一个气动计算中简单而有效的差分格式。此差分格式的特点是保持原显式格式的精度,计算过程有显式格式的简单,又有隐式格式的特点,即对于任意步长此格式都是稳定的。本文利用此格式求解了N-S方程,给出了二维激波与边界层干扰的分离流计算结果,并与MacCormack的隐式格式进行了比较。计算结果表明,本文提出的格式与MacCormack隐式格式相比,大大简化了计算过程,节省了计算机时。     

球锥-柱-裙组合体超音速无粘流场的数值计算

本文叙述组合体超音速无粘流场的一种数值计算方法,跨越流场内部的间断面采用MacCormack二步显式差分格式(下称内间断“捕捉法”),在物面边界上混合使用两种计算方法。文中还论述了MacCormack二步格式的“异步”程序设计,对于低存贮量的计算机,可成倍提高计算速度。本文方法适用于钝锥、钝多锥,幂次体及钝锥-柱-裙组合体等轴对称物体,也可推广于身部非对称体的计算。大量计算表明,该方法简便、省机时、结果可靠、适用性强。     

一种新的用于MHD方程逆风格式的Jacobian矩阵分裂方法

针对理想MHD方程,提出了一种新的基于MacCormack算法的雅可比矩阵分裂方法,克服了原有方法稳定性差的问题,并成功地应用于理想MHD方程的求解.控制方程在非结构混合网格上进行空间离散,其中对流项采用本文发展的逆风向量分裂格式,并引入了双曲型磁场散度清除技术,时间推进为显式5步龙格-库塔方法.对MHD激波管流动和带均匀磁场干扰的二维高超声速钝头体绕流流场进行了数值模拟,得到了与参考文献相吻合的数值结果,表明本文发展的数值分裂方法可以有效地捕捉MHD流场的流动特征,并且具有比MacCormack方法更高的稳定性和计算精度.     

向量化多重网格的Mac Cormack显格式解轴对称底部干扰流场

本文将多重网格法与MacCormack显格式的算法相结合,计算了轴对称底部流场中斜向喷流与超声速外流粘性干扰流场,提高了计算效率。流场的分析表明分辨率较好,计算合理。在YH-1上将程序向量化,大大节省了计算时间。     

介绍一种适合于计算跨音速流动的混合差分格式

本文介绍MacCormack格式和迎风格式的一种混合差分格式。它适用于跨音速流计算。本文用它计算了绕圆弧翼型的跨音速流场,并详述了处理方法和结果。     

用时间相关法计算喷管跨声速流场

本文用时间相关法完成了定常、无粘、跨声速喷管流动的计算。内部点参数用MacCormack差分格式计算,壁边界点参数用简化了的特征边界条件计算,轴上点参数用反射原理计算。计算表明:计算是收敛的,计算结果与实验数据符合良好。     

关于跨音速叶栅物面边界条件的数值试验

本文讨论了跨音速叶栅物面边界的数值处理方法,目的是寻求与我们所用内点格式相匹配的边界上流动参数的算法,以提高计算精度。本文首先推导了适合于物面边界计算的特征关系式,其次总结了五种边界处理方法,进行了数值试验。计算结果表明,配合我们的内点计算格式,以y向抛物外推的Denton修正法计算效果较好,能够较精确地反映跨音速叶栅带激波流场的物理特点。     

压气机叶片排中三维流场计算的多网格法

以三维有限差分形式的MacCormack格式为基础, 对相对圆柱坐标系下弱守恒型控制方程, 导出了残差在粗、细网格上的传播与回插公式, 形成了一种多网格算法。计算表明, 这种算法大幅度地提高了计算收敛速度, 而且方法本身简便易行。      

超声速翼—身组合体Euler方程数值解

本文采用MacCormack二步显式差分格式,用空向推进法给出了超声速翼-身组合体Euler方程的数值解。应用薄翼假设生成简易的多区网格。采用Kentzer方案处理边界。计算中通过对横流速度型的修正改善了背风面压力分布特性,解决了推进中断问题,同时模拟出背风涡。     

轴对称喷管流场分析的有限体积法

本应用有限体积法对轴对称喷管的跨音速流场进行了分析，采用显式预校两步MacCormack格式进行计算。运用隐式残差平均技术及当地时间步长来加速计算的收敛，提高格式的稳定性。对三个不同几何参数的轴对称喷管的数值分析表明，计算结果与试验数据吻合良好。     

采用二阶迎风格式的 PISO 算法研究

在ＰＩＳＯ算法基础上,发展了一种新的二阶迎风格式以提高计算精度又保证算法的稳定性。新格式应用于液体火箭发动机燃烧流动过程计算,其数值稳定性很好。     

一个边界差分格式及其在定常轴对称喷管流动计算中的应用

本文工作进一步研究了作者在[1]中提出的一维迎风格式。此守恒格式在边界上仅需把其落在边界外部部分置零,所得的方程近似于正确的特征关系式和“无反射边界条件”。对定常问题可先用格式计算出边界上的值,然后用给定的边界条件调整,这尤其对隐式格式是很方便的。§1讨论一维情况。§2将边界格式推广到二维定常问题的计算,这时特征关系对应这样的特征面,其法向量在二维空间的投影垂直于边界。为了把方法试用于定常轴对称喷管流,在§3给出控制体积离散方程和与它有关的矩阵性质,内点采用了全隐中心ADI格式;由于此内点格式的稳定性问题和超音区可能产生的斜激波问题,每层加了平滑。对若干喷管使用了正交和非正交格网,试算结果在§3给出,计算是收敛的并得到的近似定常解是合理的,即使在较粗的网格上也是如此。     

多重网格法求二元跨声速流在一种半自适应网格上的欧拉解

本文采用MacCormack两步显式格式,用有限体积法求解了二元跨声速欧拉流。推导了物面边界条件,采用了特征远场边界条件及远场环量修正。利用保角变换方法生成O型贴体网格,并修改得到了一种在激波处局部加密的半自适应贴体网格。采用多重网格及焓阻尼加速收敛技术计算了NACA0012翼型的跨声速气动特性,得到了十分满意的收敛过程和计算结果。     

振荡叶栅三维非定常流动数值分析

以全三维、非线性、非定常欧拉方程计算作用于振荡叶栅的非定常气动载荷。欧拉方程的数值离散采用了具有时、空二阶精度的MacCormack显式分裂的两步格式,用时间推进方法得到其数值解。应用动网格技术捕获叶片的运动,保证叶片振动过程中网格边界与叶片表面保持一致。以相移周期条件处理具有不同相位的叶片振动问题,同时发展了伪边界形状修正方法以简化周期边界条件处理。     

绕弯曲头锥超音速流动差分法数值计算

本文给出了绕弯曲头锥理想气体超音速流动差分数值计算方法。这是一个三维绕流问题,由于弯曲头锥的外形特点,头锥部分的数值计算采用了计算轴移动法。数值计算采用MacCormack二阶精度二步格式,物面上用一组由特征相容条件导得的方程。本方法已经用BCY算法语言编制成计算程序,在109丙电子计算机上进行运算。本文最后给出计算实例。     

时域推进与边界元混合瞬态热传导问题

对时间变量在小区间内应用多项式插值,对空间变量采用边界元求解,提出一种新的分析瞬态热传导问题的计算格式。与现行的利用含时间变量的基本解所得到的边界元计算格式相比较,具有基本解简单易得、系数矩阵容易计算以及无重复积分等优点;与向前差分消去时间因子的边界元计算格式相比较,则在算法稳定性、计算精度以及时间步长的选取等方面具有明显的优势,并且很容易推广到求解非线性问题。算例表明,采用时域推进与边界元混合方法分析瞬态热传导问题有效可行。     

MacCormack格式在计算非定常激波时的特点

前言 MacCormack格式(以下简称为MC格式)近年来在国内外用得较多。在非线性气体力学问题的计算中,它与著名的L-W格式或L-W二步格式相比有许多不同。其主要差别是MC格式的数值解中激波后的过头(Overshoot)较小。对于多维问题,MC格式对时间变量分步(Time Splitting)是十分有效的。 为了更好地了解MC格式,本文讨论MC格式的一些性质,并计算了两个典型问题。     

固体火箭发动机一维非定常两相喷管流场计算

本文采用MacCormack预校二步差分格式对一维非定常两相流动喷管进行了数值求解,得到了固体火箭发动机压力建立过程中气相和凝相颗粒参数的变化规律.