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

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

大迎角振动三角翼跨声速粘性绕流数值分析

本文通过数值求解三维可压缩雷诺平均ＮＳ方程，研究发析了振动三角翼的三维复杂非定常流场，其中对流场项计算格式的空间离散在平行物物理面方向与物面法线方向的空间导数分别采用二阶ＮＮＤ格式和由ＮＮＤ格式高阶插值而得到的三阶迎风格式，时间推进为隐式ＬＵ分解方法，并使用了修正的代数ＢＬ湍流模型模拟高雷诺数的湍流流动。     

NND格式的推广及在粘流计算中的应用

本文把ＮＮＤ格式推广到了特征变量，守恒变量和原始变量形式。通过引入通量差分分裂方法，由守恒变量形式的ＮＮＤ格式得到了Ｒｏｅ的近拟Ｒｉｅｍａｎｎ解法，并讨论了限制器的作用。在隐式ＮＮＤ格式中引入了一种分解算法，提高了计算效率。本文的数值计算表明，ＮＮＤ格式对复杂波系干扰，边界层分离与再附等流动现象均有较高的分辨率。     

适于黏性计算的自适应笛卡儿网格生成及其应用

复杂外形的黏性网格生成已经成为计算流体力学(CFD)在工程应用中的主要瓶颈,高效而鲁棒的网格生成技术研究显得尤为迫切。使用基于物面外形自适应加密的笛卡儿网格方法生成计算网格,投影方法拟合壁面边界,将投影得到的柱形单元沿法向分层获得适于黏性计算的网格。应用特征恢复技术恢复凹面特征并改善凸面处的网格质量。流场解算中应用最小二乘法重构获得2阶精度,Roe格式计算无黏通量,Holmes-Connell方法计算黏性通量。加入Spalart-Allmaras一方程湍流模型以获得湍流解。上下对称高斯 赛德尔迭代(LU-SGS)格式隐式时间推进,对网格重排序保证格式的稳定性和平衡性。算例显示网格生成方法是快速的和鲁棒的,流场计算精度满足工程需求,因此该方法非常适合应用于实际工程计算。     

适于黏性计算的自适应笛卡儿网格生成及其应用

复杂外形的黏性网格生成已经成为计算流体力学(CFD)在工程应用中的主要瓶颈,高效而鲁棒的网格生成技术研究显得尤为迫切.使用基于物面外形自适应加密的笛卡儿网格方法生成计算网格,投影方法拟合壁面边界,将投影得到的柱形单元沿法向分层获得适于黏性计算的网格.应用特征恢复技术恢复凹面特征并改善凸面处的网格质量.流场解算中应用最小二乘法重构获得2阶精度,Roe格式计算无黏通量,Holmes-Connell方法计算黏性通量.加入Spalart-Allmaras一方程湍流模型以获得湍流解.上下对称高斯-赛德尔迭代(LU-SGS)格式隐式时间推进,对网格重排序保证格式的稳定性和平衡性.算例显示网格生成方法是快速的和鲁棒的,流场计算精度满足工程需求,因此该方法非常适合应用于实际工程计算.     

小迎角轴对称体跨声速绕流的一种Euler解法

应用小迎角展开的方法将小迎角下旋成体绕流Ｅｕｌｅｒ方程进行摄动展开，对零介和一阶方程分别采用隐式ＮＮＤ格式进行了计算，与实验进行比较得到了较为满意的结果。该方法采用二维的计算得到了三维的结果，大大提高了计算的效率。     

带配平翼钝体高超声速粘性绕流的数值模拟

采用ＮＮＤ差分格式，通过求解Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，数值模拟了带配平翼钝体的高超声速粘性绕流。文中首先以球头驻点压强和热流为参数，讨论了在差分方程左端采用ＬＵ－ＳＧＳ、ＬＵ－ＡＤＩ和对角化ＡＤＩ三种不同隐式算法收敛效率的异同。然后以带配平翼钝体的高超声速粘性绕流为模型，对所研制的程序进行了计算验证。在计算中采用了代数方法和求解椭圆型方程方法相结合的网格生成技术，针对配平翼外形给出了贴体性、     

三维自适应非结构网格的Euler方程解

 将 Ausm+ 迎风格式应用于三维非结构网格中求解 Euler方程。对单元变量进行重构以获得空间高阶精度,对时间域采用多步龙格库塔法推进,并采用了当地时间步长和隐式残差光顺技术来加速收敛。采用多点择优推进阵面法生成复杂曲面的三角形网格,利用推进阵面法生成四面体网格。采用网格自适应技术对网格进行局部加密,以减少总体网格数目,从而提高计算效率。最后给出了绕 ONERA M6机翼的跨音速流动及绕麻雀 导弹的超音速流动算例,结果表明了本方法的有效性。     

11组元轴对称热化学非平衡流场的数值研究

本文从轴对称热化学非平衡Ｎ－Ｓ方程出发,采用双温度、十一组元反应气体模型,利用隐式ＮＮＤ有限差分格式的时间预处理技术数值求解了ＦＩＲＥＩＩ飞船热化学非平衡流场。文中同时给出了七组元模型的计算结果,并进行了结果比较。     

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

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

高超声速滚转阻尼导数数值模拟

采用非定常Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程描述物体简谐振动流场,并在Ｅｔｋｉｎ理论下给出绕定轴转动时滚转阻尼导数的计算公式,定常流场的计算采用空间二阶精度的交替方向隐式分解的ＮＮＤ格式,非定常流场的计算采用时、空二阶精度的Ｒｕｎｇｅ－Ｋｕｔｔａ多步格式,采用代数和方法生成物体静、动网格。最后给出高超声速下钝体外形滚转阻尼导数的计算结果,以及滚转力矩系数随瞬时振幅的变化曲线。     

高超声速再入体表面热流数值模拟研究

本文采用计算效率较高的标量对角化隐式NND格式，通过求解Navier-Stokes方程对影响再入体表面热流计算准确的诸因素进行了综合分析，研究了Steger-Warming通矢量分裂、Van Leer通矢量分裂和通量差分分裂方法及相应熵修正方法对热流的分辨能力，并阐明了在物面边界上采用二阶中心格式、二阶中心和二阶迎风混合格式、以及一阶迎风格式等不同边界格式对热流计算的影响。在此基础上，采用通量差分形式NND格式对钝锥和钝双锥高超声速粘性绕流进行了数值模拟，计算给出了与试验结果相吻合的热流分布。     

声速再入体表面热流数值模拟研究

本文采用计算效率较高的标量对角化隐式NND格式，通过求解Navier-Stokes方程对影响再入体表面热流计算准确的诸因素进行了综合分析。研究了Steger-Warming通矢量分裂、VanLeer通矢量分裂和通量差分分裂方法及相应熵修正方法对热流的分辨能力，并阐明了在物面边界上采用二阶中心格式、二阶中心和二阶迎风混合格式、以及一阶迎风格式等不同边界格式对热流计算的影响。在此基础上，采用通量差分形式NND格式对钝锥和钝双锥高超声速粘性绕流进行了数值模拟，计算给出了与试验结果相吻合的热流分布。     

矩形/非结构混合网格技术及在二维/三维复杂无粘流场数值模拟中的应用

建立了一套模拟复杂无粘流场的矩形／非结构混合网格技术,其中非结构网格仅限于物面附近,发挥非结构网格的几何灵活性,以少量的网格模拟复杂外形;同时在以外的区域采用矩形结构网格,发挥矩形网格计算简单快速的优势,有效地克服全非结构网格计算方法需要较大内存量和较长ＣＰＵ时间的不足。混合网格系统由修正的四分树（２Ｄ）／八分树（３Ｄ）方法生成。本文将ＮＮＤ有限差分与ＮＮＤ有限体积格式有机地融合于混合网格计算,消除了全矩形网格模拟曲面边界的台阶效应,同时保证了网格间的通量守恒。该算法被成功地运用于ＮＡＣＡ００１２翼型的跨、低超声速绕流,二维弯管内运动激波的非定常流,以及三维捆绑火箭的超声速绕流的数值模拟。计算结果表明本方法在模拟复杂无粘流场方面的灵活性和高效性。     

跨音速压气机级三维紊流流场数值模拟

在对跨音速压气机级的动／静双叶排中的三维定常紊流流场数值的模拟中，利用三阶高分辨率ＮＮＤ格式和ＬＵ－ＳＧＳ隐式推进迭代法，既保证了流场中激波的模拟质量和粘性流动特征的正确预估，也实现了求解过程的高效率。对某单级跨音速压气机的计算及其与实验数据的对比证实了本文方法的有效性。将计算推广至多级轴流叶轮机械的情况将会是十分简便的     

高超声速化学非平衡欧拉方程数值模拟

本文给出了求解高超声速化学非平衡无粘流场的数值模拟方法。计算中采用守恒变量型的ＮＮＤ差分格式和全流场捕捉技术;建立了（数学）奇性轴上的控制方程,克服了通常在其上插值方法所引起的扰动,改善了收敛性和计算精度;化学反应模型由五组元五反应组成;温度和压力的计算是非迭代的,有别于国内外惯用的迭代算法;此外,还发展了背风过低压区的处理方法,从而提高了欧拉方程的应用范围,使无粘非平衡绕流的计算攻角达到３０°     

类升力体外形俯仰阻尼特性数值研究

采用有限差分方法求解薄层近似的非定常Navier-Stokes方程，定常流场采用交替方向隐式分解的NND格式，非定常流场采用四步Runge-Kutta方法，在保证时、空二阶精度前提下引入变系数残值光顺技术提高非定常流场计算效率，复杂带翼外形的空间网格通过求解抛物化的椭圆型方程生成，最后在Etkin理论下给出球锥及类升力体外形的俯仰阻尼导数计算结果。     

再人体复杂外形无粘流场数值计算研究

本文在文[1,2,3]的基础上研究了非定常Euler方程的推进迭代方法,并将无波动、无自由参数的耗散差分格式(下称NND格式)发展为隐式、迎风格式,用以计算有差曲控制翼的再入复杂飞行器的育攻角、有侧滑角的超声速无粘绕流流场。其结果是准确的。本文的方法具有以下特点: 1.可以在同一计算程序中实现亚跨超流场的计算:在超声速区域,采用推进技术;在亚跨声速区,则自动选择推进迭代求解技术。 2.NND格式不仅具有公式简单,无自由参数的特点,而且过激波时满足熵条件,不会产生非物理解,也不发生振荡。 3.隐式差分格式的无条件稳定和Gauss-Seidel迭代的快速收敛,计算省时。 本文的方法曾用于烧蚀头部凹陷外形和航天飞机简化外形的计算,能自动选择迭代和推进区域。可以预言,如计算机条件许可,本方法可以用于大攻角的无粘流场计算。     

关于建立高阶差分格式的问题

为了能在不太密的网格上捕捉到流场的细致结构。通常采用高阶精度的差分 式进行数值模拟。为能抑制计算中和在激波附近产生的虚假波动,本文从物理构思出发,提出了建立高阶格式的两个基本原则,作者称之为抑制波动的原则和稳定性原则。