N—S方程数值解用于SB—03模型的开孔洞壁干扰修正

将一个孔及其邻近区域做为开孔壁流态的基本单元，对其求解Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，找出穿透流量对附近诸如压力、速度等流动参数的依赖关系。根据这种透气规律提出的开孔壁边界条件改善了高速风洞洞壁干扰修正的准确度。用椭圆型差分方程组生成了贴体座标网格。在有限体积构架下，用隐式的ＬＵ－ＡＤＩ因子分解对角化方法、显式和隐式ＴＶＤ方法、以及显式的多步龙格－库塔时间推进方法求解完全的Ｎａｖｉｅｒ－Ｓｔｏｋｅ     

小轿车绕流场数值计算和分析

本文求解了雷诺平均Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程,湍流的模拟采用了标准的两方程模型,修正的两方程模型,亚格子模型的大涡模拟三种要用有限体积方法离散求解Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程。对流项的计算采用三阶迎风格式。为了提高求解过程的稳定性,采用显隐式混合格式来获得主对角元素区域占优。     

求解Euler／Navier—Stokes方程的有限体积龙格库塔方法

提出了一个求解定常Ｅｕｌｅｒ／Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程的计算方法,求解的主管方程是积分形式的Ｅｕｌｅｒ方程或雷诺平均的Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程,计算方法是有限体积多步龙格库塔方法,通过在格式中加入自适应的二阶加四阶人工粘性项以保证数值计算的稳定性,用空间变时间步长、焓阻尼来加速定常问题的收敛速度。     

拟压缩性方法数值模拟绕弹体的不可压缩层流脱体涡流动

本文应用在连续生方程中增加压力对时间导数项的拟压缩性方法，数值求解三维定常不可压缩Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程。数值求解过程采用隐式的Ｂｅａｍ－Ｗａｒｍｉｎｇ近似因式分解格式。计算绕尖头旋成柱体的大迎角脱体涡流动，得到与实验相符的结果。     

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

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

平板绕流的抛物化Navier－Stokes方程分析

对于绕有限长平板不可压粘性层流问题,循抛物化Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ（ＰＮＳ）方程,进行ＰＬＫ坐标摄动分析,推得与Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ（ＮＳ）方程结果完全一致的摩擦阻力规律。该规律在１５＜Ｒｅ＜１０＿４范围内成立,证实ＰＮＳ方程适应中等Ｒｅｙｎｏｌｄｓ数流动。     

考虑粘性时激波反射及壁面热流计算

本文采用显式ＴＶＤ格式，求解薄层Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，计算了激波与楔面相互干扰的流场情况，得出了马赫反射的各种等强度级及固壁上详细的物理分布，根据计算得出的流场，计算了壁面摩擦阻力及热流分布。     

再入飞行器迎风面边界层转捩的数值预估

基于求解三维Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程,首先数值预估了高超声速三维流动的边界层特性;采用已有的高超声速三维边界层转捩准则,即在转捩位置Ｒｅθ／Ｍｅ＝常数,以及通过修正由ＯＮＥＲＡ提出的代数转捩模型,建立了预估三维边界层转捩的半经验数值技术并应用于三维Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程的求解。利用ＨＡＬＩＳ模型及相关的实验数据,首次进行了在Ｍ∞＝６时不同雷诺数及迎角的验证试验。预估结果表明：在不考虑     

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

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

旋转矩形通道中道中粘性流动的数值模拟

通过有限差分方法求解时间相关的Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，对旋转矩形通道中粘性流动进行了数值研究，旋转轴方向与主流方向垂直，给出了流分发展时矩形通道中的二次流结构及主流速度分布，得到的数值解与实验吻合良好。     

飞船高超声速粘性绕流的数值模拟

本文采用稳式ＮＮＤ格式，通过求解三维薄层近似Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，数值模拟了飞船高超声速粘性流场。在计算中充分考虑了网格雷诺数和网格分布对计算结果的影响，不但得到了满意的表面压力、热流分布和流场的物理量分布，而且得到与实验观察一致的表面油流图像和分离结构。     

旋转矩形通道中粘性流动的数值模拟

通过有限差分方法（ＭＡＣ－ＡＦ１）求解时间相关的Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，对旋转矩形通道中粘性流动进行了数值研究，旋转轴方向与主流方向垂直，给出了充分发展时矩形通道中的二次流结构及主流速度分布，得到的数值解与实验吻合良好。     

高分辨率有限差分—有限元混合方法及其在气动热计算中的应用

在文献［１］结合ＮＮＤ格式思想＾［２］提出的有限元格式基础上，给出一种有限差分－有限元混合方法。通过求解完全Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，得到了高超声速情况下钝头体轴对称绕流的较满意结果，同时讨论了网格Ｒｅ数对驻点热流的影响。     

一种铁磁流体密封系统的理论研究

利用磁路方法和简化的Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，导出了锥形空间内磁场强度分布和旋转时铁磁流体稳定界面的计算方程。分析了密封间隙、密封角和转速对密封能力的影响。     

高空台扩压器的数值模拟

从雷诺平均的非定常Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程出发,利用时间相关法,采用ＭａｃＣｏｒｍａｃｋ两步显格式,结合Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型,数值模拟了高空台扩压器的内部流场,并计算了扩压器的性能。结果表明,采用求解Ｎ－Ｓ方程的方法计算高空台扩压器性能是一种有效的方法。     

双三角翼大迎角绕流的并行计算

本文在ＬＵ－ＳＧＳ方法求解可压缩Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程的串行程序基础上,发展了适合ＭＰＰ系统类型的并行计算机的并行算法程序,并针对双三角翼大迎角绕流问题,就不同的并行计算方法进行了对比研究,探讨了优化并行策略,提高计算效率和计算加速比等问题。     

几何微扰动对脱体涡非对称性影响的数值分析

利用拟压缩性方法，Ｂｅａｍ－Ｗａｒｍｉｎｇ近似因式分解格式数值求解三维定常不可压缩Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程采用Ｄｅｇａｎｉ－Ｓｃｈｉｆｆ修正的Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型，计算了绕尖头正切拱型旋成柱体的大迎角大雷诺数脱体涡流场，计算结果表明，在柱体头部施加微小的非对称几何扰动时，可以得到与实验相符的非对称脱体涡。     

边界层转换阶段直接数值模拟的影响矩阵法

提出一种求解时间相关、三维不可压完全Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程的影响矩阵法。考虑流动区域为上半空间的平极边界层流动,在水平方向,利用周期性条件,采用Ｆｏｕｒｉｅｒ离散,在半无穷方向,定义了一种特殊映像并利用Ｃｈｅｂｙｓｈｅｖ多项式展开。由此将Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程转化为以水平波数（ｋ１,ｋ２）为参数的一组法向Ｈｅｌｍｈｏｌｔｚ方程组,同时可更有效地求解压力。本文讨论了连续问题的影响矩阵法与原问题的等价性。影响矩阵法的离散形式,以及影响矩阵法的算法实现等。最后给出了一维Ｈｅｌｍｈｏｌｔｙ问题的数值实验结果。     

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

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

跨声速操纵面嗡鸣的数值研究

我国的一些高速歼击机，在飞行中曾多次出现操纵面嗡鸣问题，它严重影响了飞机的性能和飞行安全，本文采用数值计算方法，对此问题做较深入的分析和研究。气动力计算采用Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，并采用ＬＵ－ＴＶＤ混合差分格式求解由此气动力解操纵面的弹性振动方程，求出现嗡鸣的Ｍ数及振动的频率、振幅等参数。计算中采用Ｂ－Ｌ代数湍流模型，网格生成采用代数方法及保角变换方法，文中讨论了隐式算法在非定常气动力计算