超声速化学反应流动的LU-SGS伪时间迭代空间推进求解

介绍了求解超声速多组分有限速率化学反应流动的伪时间迭代(lower-upper symmetric-Gauss-Seidel,LU-SGS)方法.空间推进求解多组分抛物化(parabolized Navier-Stokes,PNS)方程时,在一个推进面上采用修正的LU-SGS方法迭代至收敛,把得到的结果作为初值赋给下一个推进面.沿推进面依次迭代直至求解完整个流场.采用伪时间迭代LU-SGS方法求解化学反应PNS方程,计算结果的准确性和时间迭代求解完全Navier-Stokes(N-S)方程相当,求解效率提高一个数量级.      

非定常流动计算的混合时间推进方法研究

结合ADI方法的二阶时间精度特性和高斯-赛德尔迭代法的良好收敛性特性,提出了一种ADI-高斯-赛德尔混合迭代格式,应用该格式,对NACA0012翼型的跨声速俯仰振荡进行了EULER方程求解,空间离散格式采用ROE’S FDS格式,计算结果与实验及成熟的LUSGS双时间步迭代方法计算结果进行了分析比较。计算结果表明,该格式精度能达到LUSGS双时间步迭代方法完全收敛的二阶精度,且该格式精度对内迭代次数的依赖性较低。分析表明,本文提出的方法具有较高的精度和较好的收敛性。     

一种计算非定常二维流动的无网格算法

主要目的是发展一套求解非定常流动的无网格算法。计算区域的离散方面,提出了一种按区域进行填充布点的点云自动生成方法;发展了一种点云的运动技术来实现离散点云对物面边界的随体运动;在点云离散的基础上,采用最小二乘法求解矛盾方程的方法来求取空间导数,进而获得数值通量;采用双时间方法进行时间离散推进,其中物理时间迭代采用二阶隐式格式,伪时间迭代采用四步龙格一库塔显式格式,为了加速收敛,在伪时间迭代中采用了当地时间步长和隐式残值光顺等加速收敛措施。最后,利用本文算法模拟了NACA0012翼型和NACA64A010翼型的跨音速非定常流动,并将计算结果与实验测量结果进行了对比分析,验证了上述方法的正确性和实用性。     

二维带动力吸气式高超声速飞行器绕流的PNS-NS混合求解

为了提高吸气式高超声速飞行器绕流的求解效率,采用空间推进求解抛物化PNS(parabolizedNavier-Stokes)方程和时间迭代求解Navier-Stokes(N-S)方程的混合计算流体动力学(CFD)方法来求解高超声速飞行器整机绕流.在超声速占主导的流动区域采用空间推进求解抛物化N-S方程的方法,在亚声速和分离区采用时间迭代求解N-S方程的方法.对于求解二维带化学反应的吸气式高超声速飞行器绕流,混合CFD方法和完全时间迭代方法相比,可得到同等准确的数值模拟结果,并且求解效率提高了数倍.      

双时间法求解大迎角翼型绕流非定常N-S方程

本文采用双时间推进法、顶点格式有限体积法求解非定常ＮＳ方程、应用多层网格技术计算了ＮＡＣＡ００１２翼型大迎角振荡绕流,并对伪时间推进收敛精度,物理时间步长,多层网格迭代方式等对计算结果与计算ＣＰＵ时间的影响进行了讨论,另外还讨论了采用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ湍流模型时用固定转捩或全湍流模式对计算结果的影响,并与实验结果进行了比较。     

基于笛卡尔网格的非定常气动力计算

在非结构自适应笛卡尔网格基础上构建了一种基于笛卡尔动网格的非定常气动力计算方法。为了适应非定常流场对网格品质要求较高的特点,网格生成采用了基于几何外形、模型表面曲率为基础的自适应生成及加密方法。采用显格式龙格库塔推进方法和中心格式空间离散的有限体积法对NACA0012翼型,三维机翼的跨音速刚性振动进行了数值模拟。对伪时间收敛精度控制,真实时间步长对结果的影响进行了讨论。将计算结果与试验数据进行了比较,得到了良好的效果。     

基于时间谱方法的振荡翼型和机翼非定常黏性绕流数值模拟

 传统的双时间方法在非定常计算中长时间的过渡迭代推进求解导致其计算效率相对较低,针对周期性非定常流动问题的流动特征,发展了一种基于离散傅里叶变换的高效时间谱方法,用于求解振荡翼型和机翼的非定常黏性绕流。在时空耦合的雷诺平均Navier-Stokes(RANS)方程的求解中,对流项的离散应用了Roe的通量差分格式,物理时间项的离散方法为时间谱方法,伪时间推进采用了隐式LU-SGS(Lower-Upper Symmetric Gauss-Seidel)格式。考虑到湍流的时空耦合效应,时空耦合的Spalart-Allmaras一方程湍流模型的物理时间项同样采用时间谱方法进行离散。为了进一步提高计算效率,当地时间步长和多重网格技术等加速收敛的措施均被采用。算例对俯仰振荡NACA0012翼型和Lann机翼的周期性非定常流场进行了数值计算。结果表明:对于周期性非定常流场的数值模拟,相比于传统的双时间方法,用时间谱方法近似物理时间项,不仅能够提高流场的计算精度,而且更能够大幅度提高计算效率。     

绕HERMES外形高超声速无粘流动的数值模拟

本文采用时间相关法和隐式推进迭代方法求解了绕ＨＥＲＭＥＳ外形高超声速无粘流动。将抛物型方程和代数方程相结合生成计算网格，既节省了网格生成的时间，又能得到较合理的网格分布。在时间相关法计算中，采用了显式ＮＮＤ格式；在隐式推进迭代计算中，采用了推广的隐式ＮＮＤ格式。     

用Newton子迭代方法计算前飞旋翼粘性绕流

通过求解Navier-Stokes方程数值模拟了直升机旋翼前飞非定常流场，为了同时保证计算的时间精确性和计算效率，时间推进格式采用了双时间推进方法，在该方法中，子迭代过程由十分高效的LU-SSOR方法完成，且使子迭代过程成为Newton子迭代，空间上应用中心平均的有限体积法进行离散，为了模拟前飞桨时间的相对运动，网格布局采用了运动嵌套网格方法，应用本文方法对-悬停流场进行了数值计算，计算结果与实验吻合较好，尽管缺乏实验数据的验证，对-有升力前飞状态的数值模拟结果是可信的。     

更准确、更精确、更高效——高超声速流动数值模拟研究进展

从准度、精度和效率3方面回顾了近几十年来高超声速流动数值模拟研究的进展。在物理模型方面,介绍了高超声速数值模拟中高温气体效应、稀薄气体效应以及湍流效应的建模与模拟,基于雷诺平均Navier-Stokes(RANS)方程重点对现阶段较为关注的高超声速边界层转捩的模式理论研究进行了介绍。在空间离散算法方面,主要介绍了高超声速数值模拟中常用的二阶精度迎风格式以及高阶精度格式的发展及其应用。在时间推进方面,主要回顾了隐式时间推进方法的发展及其应用。在误差和不确定度估计方面,主要介绍了其概念、来源以及常用的分析方法,同时给出了迭代误差估计、Richardson外插法以及敏感性导数方法等初步研究结果。最后,讨论了高超声速流动数值模拟中下一步需关注的问题。     

高超声速飞行器单壁膨胀喷管的自动优化设计

为了提升高超声速飞行器单壁膨胀喷管的升力和推力,对其进行了优化设计。介绍了超声速流场的分区推进求解技术、流动的时间迭代和超声速流动的空间推进求解方法。对高超声速飞行器的单壁膨胀喷管进行了参数化描述。实现了喷管二维网格的自动生成。借助成熟的单目标和多目标优化软件,对喷管的推力和升力进行了优化。优化后的喷管在推力和升力方面有了大的提高。所采用的超声速流场分区推进求解和空间推进求解技术使得优化过程在单个CPU上能快速完成。     

有限体积时间推进法在跨音速扩压器湍流流场中的应用

本文提出了一种有限体积时间推进法,将其应用于均化N-S方程的求解,模拟跨音扩压器的湍流流场,并采用空间变时间步长和多重网格技术加快收敛速度。计算结果与有关文献的实验数据吻合良好。     

基于Rao方法的二维单壁膨胀喷管优化设计

针对高超声速飞行器后体尾喷管的优化设计问题,发展了与多目标优化程序NSGA-Ⅱ相结合的自动优化流程。利用Rao喷管的近似方法建立喷管优化模型,采用区域推进求解PNS方程的流场解算器和网格自动化生成技术获得尾喷管的性能;引入目标约束设计和自适交叉技术,自动加速目标空间的搜索,大幅度提高了多目标优化设计效率;获得了给定来流和喷流条件下的一种具有良好性能的二维后体尾喷管设计方案。     

基于非结构网格空间推进的高超声速飞行器背脊线优化

针对高超声速飞行器背脊线优化进行了研究。结合自由曲面变形(FFD)方法和Hicks-Henne型函数建立背脊线参数化模型,保证了飞行器背部曲面光滑连续。为减少计算耗时,采用非结构网格空间推进方法求解高超声速无黏流场,通过钝锥算例对该方法计算效率和精度进行了参数分析,并对飞行器流场计算进行了验证,对比表明表明:非结构网格空间推进算法计算效率提高4.5倍以上,计算结果与常规时间迭代法符合较好。采用优化拉丁超立方抽样建立响应面模型,分析了各设计变量对目标函数的影响。采用非支配排序NSGA-Ⅱ(non-dominated sorting genetic algorithm-Ⅱ)算法以阻力系数和升阻比为目标对背脊线优化,约束飞行器体积变化不超过20%。结果表明,飞行器阻力系数和升阻比呈负相关,优化后阻力系数降低7.1%,升阻比提高16.9%。      

计算网格对气动力气动热数值模拟影响的研究

以二维圆柱钝头超声速流场为研究对象,应用计算流体力学方法研究了计算网格在气动力及气动热数值模拟中的影响.计算程序中对流项、扩散项分别采用AUSM格式和中心差分格式离散,并用多步Runge-Kutta显式时间推进法求解空间离散后的控制方程.计算结果表明,壁面法向第一层网格间距、网格长宽比以及扩张度等因素都会对壁面摩阻、热流的计算产生影响,尤其是壁面法向第一层网格间距影响最为突出.边界层的计算结果与实验数据进行了仔细对比,包括壁面压力分布,摩阻系数Cf分布,斯坦顿数St分布等.高马赫数滞止气流导致壁面附近产生巨大速度梯度以及温度梯度,这些梯度的准确计算需要合适的壁面法向第一层网格间距以及合理的网格长宽比、扩张度.     

使用GAO—YONG湍流模型数值研究管道凸起流动

用GAO—YONG可压缩湍流方程数值模拟了Delery管道凸起跨音流场中的激波湍流边界层干扰现象。分析了GAO-YONG可压缩湍流方程组对湍流的非平衡、多尺度、各向异性等特性的描述能力。计算中对流项、扩散项分别采用二阶ROE格式和二阶中心差分格式离散,并用多步Runge-Kutta显式时间推进法求解了空间离散后的控制方程。计算很好地模拟到了压力平台区、“入”波结构等典型激波湍流边界层干扰的流动现象,也得到了壁面压力分布、平均速度剖面以及雷诺应力分布等,并与相应的实验数据进行了对比分析,两者符合很好。     

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

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

飞机亚跨超绕流的数值模拟研究

本文以Euler方程为数学模型,采用一种高精度的TVD(Total Variation Dimishing)离散格式及一种含近似因式分解的推进迭代方法,求解亚跨超绕流’流场。通过若干算例的试算,证明方法是可行的,可以用来模拟飞机的复杂流场。