用TVD方法模拟喷管内的横流流场

采用显式的ＴＶＤ ＭａｃＣｏｒｍａｃｋ格式，并结合Ｂａｌｄｗｉｎ－Ｌｏｍａｘ的代数湍流模型和低雷诺数两方程ｋ－ε模型分别求解了喷管内横流流场的二维Ｎ－Ｓ方程，得到了与理论分析结果相近似的流场结构图谱，并对代数湍流模型和两方程模型得到的结果进行了比较，表明运用两方程模型所得到的结果稍好于代数模型。     

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

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

跨音速机翼结构动力响应的计算研究

通过交替求解流体运动方程和结构动力方程来计算时域内二维机翼的结构动力响应,并研究了结构参数对系统稳定性的影响,使用控制中的相平面图对计算结果进行分析,来判断系统的结构稳定性。计算中流体运动方程采用二维非定常Ｅｕｌｅｒ方程,并用Ｊａｍｅｓｏｎ提出的有限体积方案,四步Ｒｕｎｇｅ－Ｋｕｔｔａ时间推进求解;结构模型采用二自由度ＰＡＰ模型。经分析,计算结果合理。     

气固两相流k-ε-kp-εp双流体模型及应用

以流体ｋ－ε双方程模型为基础，提出了一种包括颗粒湍动能输运方程和耗散率方程的ｋ－ε－ｋｐ－εｐ两相湍流模型。并对两个垂直上升管道内气固两相湍流进行了数值计算，计算结果和已有实验结果吻合得很好。     

用预处理方法数值模拟街区峡谷内的污染物扩散

通过耦合求解二维可压缩Navier-Stokes(NS)方程和k-ε湍流模型方程及平流扩散方程，数值模拟了街区谷（Canyon)内的低速湍流流场及大气污染物在其中的扩散。为了加快计算的收敛，对控制方程进行了预处理。预处理后的NS方程的求解采用LU－SGS隐式方法和JST中心差分格式，数据模拟结果表明预处理方法求解低速流场可加快敛，同时就街峡欲内的污染物扩散状况得出了有意义的结论。     

低速大迎角三角翼旋涡流动的有粘怀无粘数值模拟

通过人工压缩性方法求解了三维定常不可压Ｅｕｌｅｒ方程和雷诺平均Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，应用Ｂｅａｍ－Ｗａｒｍｉｎｇ近似因子分解格式及其对角化形式，采用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型。计算模型为７０°三角翼，分析了有粘与无粘流动旋涡特性的影响，尤其是对Ｅｕｌｅｒ方程模拟大迎角分离涡、涡破裂的能力及其局限性提供了数值验证，计算结果与实验结果相一致。     

对 k-ω 湍流模型的改进

为了更适合于较复杂几何外形的数值计算，本文提出仅引入１个参数来改进ｋ－ω湍流模型并采用了一种双网格系统的有限体积格式来求解Ｎ－Ｓ方程与改进的的ｋ－ω模型方程。通过平板附面层的计算结果表明了用本文提出的对ｋ－ω湍流模型的改进对控制转捩是有效的。     

二维超音速湍流混合特性研究

着重讨论顺流喷射情况下的超音速混合的流场特征及影响因素。简要介绍了顺流喷射混合的研究现状。同时，探讨了不同的来流参数对流场特征的影响及混合机理。计算采用了一种二阶ＴＶＤ格式有限体积法及Ｂ－Ｌ湍流模型求解二维湍流Ｎ－Ｓ方程。结果表明：所采用的计算方法能成功地抑制非物理振荡，适用于二维超音速内流的湍流混合特性研究。     

二维及三维贴体网格的一种生成方法

本文讨论了二维及三维情况下应用Ｌａｐｌａｃｅ方程生成贴体网格的方法，并应用ＡＦ１格式求解方程以大大加速迭代收敛的速度，对于二维及三维情况，分别讨论了相应用Ｄｉｒｉｃｈｌｅｔ边界问题解法。对于二维情况，还讨论了翼型的Ｏ－、Ｃ－和Ｈ－型三种典型的网格生成方法。着重讨论了有关割缝上的网格点的处理方法，以保证网格的光滑性。三维情况下，以简化航天飞机外形为例，讨论了网格的生成方法，尤其着重讨论了边界点的处理     

基于非交错有限体积格式和 k-ω 湍流模型模拟叶栅湍流分离流动

发展了一种求解任意坐标系下二维Ｎ－Ｓ方程的非交错有限体积法和一种高阶ＴＶＤ对流格式。采用一种ｋ－ω湍流模型和几种不同的ｋ－ω模型对叶栅湍流分离流动进行了详细的分析。数值结果与得到的实验结果比较表明,此方法不仅精度高,而且具有较好的数值稳定性。     

大迎角粘性复杂流场的数值模拟

本文使用于存在较大分离区的Ｊｏｈｎｓｏｎ－Ｋｉｎｇ湍流模型对多段翼型流场及某型歼击机外形全机流场的雷诺平均ＮＳ方程进行数值计算，研究了多段翼型及有三角翼的全机外形的最大升力和三角翼翼尖涡，二次涡等粘性特性。     

基于非交错有限体积式和k—w湍流模型模拟叶栅湍流分离流动

发展了一种求解任意坐标系下二维Ｎ－Ｓ方程的非交错有限体积法和一种高阶ＴＶＤ对流格式。采用一种ｋ－ω湍流模型和几种不同的ｋ－ω模型对叶栅湍流分离流动进行了详细的分析。数值结果与得到的实验结果比较表明,此方法不仅精度高,而且具有较好的数值稳定性。     

超燃冲压发动机燃烧室三维化学反应流场的数值模拟

采用ＬＵ隐式方法强耦合求解有限体积法离散的完全Ｎ－Ｓ方程及组分方程，数值模拟超燃冲压发动机机燃烧室机的化学反应流流场。计算中和了ＡＵＳＭ通量分裂格式，两方程的化学反应模型和Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代化学反应流场。     

二维超燃冲压发动机尾喷管优化设计

为得到二维超声速燃烧冲压发动机尾喷管最优型面，建立了二维尾喷管基于N－S方程的优化设计模型，运用复合形优化方法对该模型进行求解，得到了尾喷管型面的优化设计，结果表明，尾喷管性能参数有明显的提高，推力系数和升力系数升高，俯仰力矩系数降低。同时，随着飞行马赫数的降低，推力系数和俯仰力矩系数都会下降，而升力系数却上升，并且初步研究了尾喷管下壁长度对尾喷管性能的影响，为超燃冲压发动机尾喷管设计提供一定参考。     

驻点壁面催化速率常数确定的研究

以平衡流动作为热环境估算的依据，提出了用数值求解非平衡Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程和实验测量热流值确定模型表面材料催化速率常数的方法。用５组分１７个化学反应Ｄｕｎｎ－Ｋａｎｇ空气化学模型和轴对称热化学非平衡Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，对激波管中球头和平头圆柱模型绕流流场进行了数值模拟，给出了驻点热流随催化速率常数变化的分布，并根据激波管实验测量的热流值确定了表面材料Ｐｔ、ＳｉＯ２、Ｎｉ和某种     

氢气超音速燃烧非常过程数值模拟

用有限差分法求二维Ｎ－Ｓ方程和组分守恒方程，使用代数涡粘性湍流模型，两步化学反应模型和ＭａｃＣｏｒｍａｃｋ时间分裂法，模拟超音速气汉中横喷氢气自动着火及火焰传播过程，清晰地描述了喷氢、火焰传播到主流、稳定燃烧过程中，温度、压力、速度和油气当量比的变化。计算结果表明，高温附面层和喷嘴和回流区能促使氢气自动着火，进口气流的静温是火焰向主流传播的重要影响因素。计算还证明了喷嘴后氢与空气两股流呈现平行流动     

跨声速三维非线性洞壁干扰的数值计算

分别以固壁条件和洞壁附近的压力分布模拟各类实壁和透气壁试验段的洞壁边界条件，利用Ｅｕｌｅｒ方程和Ｎ－Ｓ方程数值求解模型在风洞中的绕流场，得出洞壁干扰对跨声速模型绕流和气动力的影响，初步的研究结果表明，该方法能较有效地模拟模型在跨声速风洞中的绕流场，经洞壁干扰修正后的ＧＢＭ－０４Ａ模型在０．６ｍ风洞中的试验结果与无干扰参考结果吻合较好     

翼—身—尾组合体绕流的Euler方程计算

本文采用重迭网格技术和Ｅｕｌｅｒ方程计算翼－身－尾组合体绕流。对翼－身与尾－身部分采用各自的Ｈ－Ｏ型网格。Ｅｕｌｅｒ方程求解采用Ｊａｍｅｓｏｎ有有限体积法，即中心差分挖和显式Ｒｕｎｇｅ－Ｋｕｔｔａ时间推进。采用前后区交替迭代使前后两区通过重迭区交换信息。本文用ＮＡＣＡＴＮ４０４１翼－身－尾模型为例，计算的空气动力特性与实验符合较好。     

跨声速开槽壁超临界翼型洞壁干扰的N－S方程模拟

本文采用二维可压缩非定常Ｎ－Ｓ方程模拟了带有上下开槽壁的跨声速二维管风洞中超临界翼型绕流的流场。网格生成中内层采用解双曲型偏微分方程，外层采用代数生成法。开槽壁采用只有压力梯度没有顶板运动的Ｃｏｕｅｔｔｅ流动的解析解近似。推导出槽间流动与开闭比的三次方关系。上下壁边界层可由无滑移条件直接模拟，也可用平板边界层近似。本文对驻室压力作了几种近似处理，结果表明，所给Ｎ－Ｓ方程的边界条件与风洞实验多接近一步，其结果就符合得更好，这为下一步用Ｎ－Ｓ方程进行洞壁干扰修正打下了基础。     

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

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