Navier-Stokes方程的参量样条分步解法

作者建议了应用参量样条求解N-S方程的按坐标方向进行分步计算的隐式格式。在每个坐标方向上,仅需求解一三对角形方程组。与SADI方法相比,具有更大的适应性,可节省计算机时,减少内存量。应用于带有斜激波的粘性流和高雷诺数下的平面Cavity流动问题的数值结果是令人鼓舞的。     

对输运方程“人工粘性”的认识和旋转的逐向分析格式

本文提出对涡输运方程差分格式的精度和所谓“人工粘性”的新认识,指出过去各种差分格式在高网络Re数情况下失败的原因在于没有真正迎风。本文提出了真正迎风的、旋转的逐向分析格式,并用它计算丁Re=100,400,1000,2000条件下的平面Cavity流动,用10×10和20×20的粗网格,从低Re数到高Re数均获得收敛迅速并且比较精确的数值解。     

燃气源开槽导管三维流场的数值模拟

为了解某火箭弹控制系统燃气源开槽导管内的气流流动特性,通过ＳＩＭＰＬＥ方法对其复杂流场进行了数值模拟。采用隐式时间分裂格式、幂函数方案及交错网格修正压力等方法,获得了燃烧室后部和导管内气流的速度、压力、密度和温度的三维分布,计算结果与实验现象的规律一致     

应用TVD格式分离求解不可压N-S流动

提出一种新的将连续方程与动量方程分离求解的TVD引用方法。对动量方程中的对流项采用三阶精度、迎风TVD格式离散, 粘性项和压力项采用中心差分格式离散。压力采用压力替代法求解。计算所用的网格是具有交错网格主要特点的修正非交错网格, 并采用一种压力修正法和G-S中心对称迭代法加速收敛。给出了二维管流和二维空穴流动在不同状态下的计算结果, 并与其它解法的计算结果或精确解进行比较。结果表明, TVD格式不仅能用于不可压流数值计算, 而且还能得到比非TVD差分法精度更高的结果。      

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

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

新的解跨音速Euler方程的隐式杂交方法

 本文提出一个新的杂交格式,它由Jameson-Turkel LU分解法和Harten TVD（total variation diminishing）格式的修正数值通量耦合而成。LU分解用于Euler方程隐式部分的离散化而Harten的数值通量控制所谓的稳态解残差。精心构造LU算子,可使矩阵算子L和U很容易逐点推进求逆;以Harten的数值通量计算杂交格式右边的通量,这样做在a=1/2的条件下,可保证收敛的稳态解具有二阶精度并能高分辨率地捕获激波。数值实验表明,本文格式对于稳态计算是有生命力的并能保证较高的激波分辨率。     

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

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

样条Lax-Wendroff格式和样条跳蛙格式

1.引言 用三次样条求解粘性流体力学问题的数值方法已有了相当的发展。为了探讨用该方法求解无粘流动问题的可能性,本文尝试把差分方法中一些典型的计算格式扩展到样条方法,提出了样条Lax-Wendroff格式和样条跳蛙格式。以一维模型方程为例,按照Warming等的方法导出了它们的修正方程,分析了格式的稳定条件及数值弥散和耗散效应。然后扩展到对一般非线性双曲方程求解的研究。最后,讨论了某些数值解例。     

求解双曲型守恒律的五阶松弛格式

给出了一种求解一维双曲型守恒律的五阶松弛格式.该格式以五阶WENO重构和显隐式Runge-Kutta方法为基础.本文格式保持了松弛格式简单的优点,即不用Riemann解算器和计算非线性通量函数的雅可比矩阵.用该格式对一维Euler方程进行了数值试验,并与三阶和四阶松弛格式的计算结果进行了比较,结果表明本文的格式具有更低的数值耗散和更高的分辨率.     

低耗散TVD格式及叶轮机内低马赫数流动模拟

推导了条件化预处理矩阵在守恒变量下的形式,以及三维任意曲线坐标系下的相应的系统特征值与特征矩阵。根据预处理技术修正了任意曲线坐标下的Harten-TVD格式,得到低耗散TVD格式,并讨论了抑制了压力速度失耦的方法。分别使用修正后低耗散格式与未修正的格式,对高负荷低压涡轮叶栅T106内低速流动:进口马赫数分别为0.1,0.01和0.001,进行数值模拟。计算表明,未修正格式不能得到有意义的结果,而修正后的低耗散格式,即使在极低的马赫数流动下,也工作良好。同时可以注意到:对于无粘流动,三种进口马赫数条件下,流场内等压线分布几乎是一致的。      

实际机身跨音速流场分析计算

本文应用修正的跨音速小扰动(MSD)方程和疏密网格交替迭代法,将Jameson格式在非等步长网格中的推广形式与Engquist-Oslier格式相结合,在机身侧边局部横向网格线上使用横向线松弛法,经数值实验得到了一个收敛速度比较快的计算方案。许多算例均证明本方案迭代次数少,计算结果与实验结果符合得也很好。     

两种湍流模型及其可压缩修正在高超压缩拐角中的评估(英文)

采用在高阶精度差分格式对SA湍流模型和SST湍流模型及相应的可压缩修正模型在高超拐角流中进行了评估。可压缩修正方式考虑了密度梯度、压力膨胀和湍流马赫数等方法。为了减小数值误差与模型误差之间的混淆,控制方程的对流项采用了5阶精度的加权紧致非线性格式(WCNS-E-5),粘性项采用了一种半结点/结点交错的4阶中心格式。通过对马赫数为9.22的15度拐角和34度拐角湍流的模拟,考察了原始湍流模型及其修正模型的效果。计算表明:原始SST模型对高超拐角湍流的预测比原始SA模型准确,这种准确主要体现在对分离区预测、再附点附近压力和热流峰值预测上。通过混合采用Catris和Shur等的方法对 SA模型进行可压缩修正可以大大改进模拟效果。在SST模型的可压缩修正方法中,Catris的修正方法最好;考虑压力膨胀修正后得到的分离区远远偏大,比本原始SST模型更差;考虑湍流马赫数的修正方法得到的分离区偏小。本文还给了部分低阶格式的计算结果,高阶格式与低阶格式相比,对分离区大小、再附点附近的压力和热流峰值等的预测准度有所改进。     

跨音速定常小扰动压力方程的混合差分法

本文提出直接求解跨音速定常小扰动压力方程的数值方法。对于研究某些洞壁干扰问题,与传统的速势方程相比,用压力方程作为求解跨音速流场的主管方程,边界条件为Dirichlet形式,易于处理,且待求变量为压力,可减少积累误差,提高计算精度。 本文采用混合差分法求解压力方程,通过数值试验,确定合适的差分格式及迭代线化方法。其收敛解与相应的以速势方程为主管方程求到的解相比吻合得比较好,从而证实了本文方法的可行性。 最后给出应用本文方法计算鉴定跨音速翼型风洞壁干扰以及由给定的压力分布计算翼型外形的典型算例。     

基于浸入式边界方法的运动物体声辐射数值模拟

为了在声学计算中简化网格的处理,从而发展出适用于运动复杂边界问题模拟的计算格式,基于浸入式边界方法,提出了一种适用于无黏流动、无穿透边界条件的处理方法,并对一些简单的运动物体声辐射问题进行了数值模拟.数值结果与解析解得到了令人满意的对比,初步验证了所提出的研究方案的可行性.     

复杂外形气动力计算的阻力修正方法

基于气动力计算的位流线化理论,发展了一种亚音速附面层阻力修正方法。采用联合流场求解方法获得复杂外形的表面压力分布,在此基础上,将绕复杂外形物体沿流向面元划分为片条,在片条上,把三维附面层转换为二维附面层,对二维附面层进行数值求解,得到附面层粘性阻力修正。算例表明,该方法计算量小,与实验和N-S数值解结果吻合较好。为复杂气动外形飞行器初期设计、方案选型及评估提供了较为理想的工程设计手段。     

带扰流片的翼型绕流流场的计算

文章用基于伪压缩方法的有效不可压Navier-Stokes解算器，（该解算器可推广应用于非定常问题），以时间精度的方式对带扰流片的翼型绕流流场进行了数值模拟。对大多数计算而言，空间离散多采用三阶通量分裂格式，对得到的矩阵方程用点SGS松弛格式来求解。用单区和多区网格生成方法分别对流场进行了数值计算，并对得到的数值计算结果进行了比较。该程序由流过NACA 0012翼型的定常湍流进行了进一步验证，并且计算得到的表面压力系数与物理试验和其它数值计算的结果符合得很好。而且该程序还显示出很好的矢量化和并行化特点得到进一步证明。计算用基本翼型是先进的波音运输机预研翼型。其几何尺寸在网格生成时稍微有些改动。文章还研究了扰流片偏转角固定为60。。并研究了攻角对流场的影响。扰流片端部和翼型后缘非定常涡脱落过程也得到了很好的模拟，同时还显示了在不同时间步时的流型。在湍流模型的局限内，计算得到的平均量和物理试验结果符合得很好。     

冲压发动机进气道掺混段与弹体内外流场一体化数值模拟

采用二阶迎风隐式TVD格式分区求解可压缩N-S方程,数值模拟了导弹、进气道和掺混段内外流一体化通气模型的复杂流场.本文给出了掺混段出口压力与来流压力比为P/P∞=3.5,4.2和4.9情况的数值模拟结果.并分析了流场特性及激波波系结构,总压恢复系数σ和流量系数ψ.同风洞试验结果进行比较,表明计算结果是合理的.文中还对该CFD分区计算内边界处理所采用的数据结构进行了阐述.     

二维超／高超声速进气道流场数值模拟

对超／高超声速三级压缩进气道流场进行了数值模拟,来流马赫数为4、6,进气道内流动为层流状态,根据二维Navier-Stokes方程,采用二阶精度Roe格式进行离散。按照流场特点,合理地设计网格分布及调整不同黏性范围的熵修正,防止了壁面附近过大的数值耗散,使计算结果更加合理。在进气道模型的各级压缩折转角处,获得了清晰的激波结构,在进气道内部的各种波系的相交、反射和激波诱导的边界层分离等现象都得到合理的描述。计算得到的压力分布,在各级压缩斜板上同简单波理论结果十分接近。用本文方法计算了另一个二级压缩进气道,沿上、下壁面的压力分布与试验比较符合得较好。     

二维高亚声速空腔流激振荡的数值模拟研究

采用基于P .L .Roe的近似Riemann解的修正Osher&Chakravarthy (MOC)三阶TVD有限差分格式 ,数值求解二维雷诺平均全Navier Stokes方程 ,并用Cebeci Smith代数湍流模型 (对腔内区域作修正 )来模拟湍流效应 ,时间方向的积分采用四阶Runge Kutta方法 ,通过对GAMM超音速前台阶绕流的计算 ,验证了格式及程序的有效性 ,对高亚声速来流下的空腔流动作了数值模拟研究 ,取得了较好的效果     

一种新型多维通量函数在求解平面叶栅跨音速流动的应用

采用一种改进的Ｒｏｅ格式对平面叶栅的无粘流动进行了数值模拟,并与Ｒｏｅ格式的计算结果进行了比较。采用四阶Ｒｕｎｇｅ－Ｋｕｔａ时间差分格式,空间三阶有限体积法。网格生成采用贴体网格近壁面修正的方法。计算结果表明,该格式在减少格式对网格的依赖性,抑制数值粘性方面具有明显的作用。