高速粘性绕流的有限体积TVD格式

本文通过求解Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程,解决二维超声速及高速声速流动问题。通过采用坐标交换技术,可以方便地将有限体积ＴＶＤ格式应用到各种形状的流动问题中,各种典型算例的计算表明该方法数值计算稳定,捕捉激波分辩率高,能很好地模拟超声速和高速超声流动中激波反射及激波与边界层的干扰等各种复杂流动现象。湍流计算中采用ＢａｌｄｗｉｎＬｏｍａｘ湍流。数值结果与实验结果及理论解的对比令人满意。     

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

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

运动激波绕尖劈流动的研究

在激波管中进行双尖劈二维外形的激波绕流数值计算和实验研究,入射激波,Ｍ＝１．８,用激光全息双爆光技术定量测定各瞬时的密度场。计算采用欧拉方程和有限体积法进行离散并采用高精度的ＴＶＤ差分格式,计算与实验两者之间的比较,表明本文采用的计算方法对于解决尖点的绕流十分奏效,根据计算所获得的运动激波系,能够判断局部区域实验测量的密度场梯度方向,从而使密度场的定量测定获可靠度的结果。     

运动激波绕射的数值模拟

用高分辨率的ＴＶＤ格式对复杂的激波反射过程做数值模拟。激波斜壁反射计算结果与实验干涉照片、压力分布、密度分布均符合得很好。通过对运动激波绕射凸起物的数值模拟,得到非定常反射流场发展过程及复杂流场间断的清晰图谱。     

高超声速球锥粘性激波层压力分裂法数值求解

本文将压力分裂法引入细长钝体高超声速粘性激波层数值求解中,较好地解决了这种流场远后身部分的推进求解问题。针对非解析球锥体绕流数值求解的困难,作者提出了实用型初始激波形状计算方法,达到了简化数值求解步骤,提高数值计算效率的目的,数值求解结果与实验结果相符合。     

N-S方程及其简化形式在圆柱绕流分离点邻域的数学性质

证明了在圆柱绕流分离点的邻域上 ,N-S方程和高智的简化 N-S方程是正则的 ,即存在级数解 ;而 Davis和 ГКП的简化 N-S与经典的边界层方程一样存在 Goldstein奇性     

激波绕拐角运动流场的数值模拟

本文采用时间和空间方向均为二阶精度的ＴＶＤ有限差分格式来数值模拟冲击波绕三角形、半菱形障碍物的非定常复杂流场，算出了与理论分析一致的流场结构，成功地捕获了激波，膨胀波，二次弱反射激波等物理特征。     

高超声速粘性激波层高雷诺数流场的数值求解方法

本文将拉伸变换引入Hosny的高超声速粘性激波层的数值求解方法中,数值计算结果与Little所给出的实验结果相符合。从而不论在雷诺数范围上,或在后身流场的推进计算上,本方法均比Hosny原方法有明显改进。     

带尖针杆的钝体粘性绕流的数值模拟

用无波动、无自由参数耗散的差分（ＮＮＤ）格式数值求解了Ｎ－Ｓ方程,模拟了尖锥柱与球头组合体的高超音速粘性流动;清晰地给出了激波、膨胀波、剪切层和旋涡等流动图象。计算给出的定常流场与实验图象甚为一致。通过数值模拟研究了激波、旋涡的形成和发展过程以及相互作用情况。发现分离区中的旋涡存在分叉、合并及二次分离等演变历程。钝体主激波和分离激波从开始形成到相互撞击;发展到在与分离区相互干扰下向前、向后运动和变形;最后趋于定常状态。     

弹丸超声速绕流的Navier-Stokes模拟

本文轴对称Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程组出发,采用高分辨率对称ＴＶＤ格式数据模拟了包括底部回流区在内的弹丸超声速绕流流场。     

N-S方程数值模拟弹体不可压脱体涡绕流

应用在连续方程中增加压力对时间导数项的拟压缩性方法 ,数值求解三维定常不可压N- S方程。采用隐式 Beam- Warming近似因式分解格式求解。以修正的 Baldwin- Lomax代数模型模拟湍流流动。对绕正切拱形旋成柱体的大迎角脱体涡流动进行了数值模拟 ,并分析了层流粘性和湍流粘性对脱体涡及二次分离涡的影响。在头部施加微小的非对称射流扰动可得到与实验相符的非对称脱体涡     

紧致格式数值模拟超音速粘性绕流问题

通过数值求解可压Navier-Stokes方程的方法计算了球锥粘性绕流问题。方程中粘性项按通常办法逼近,为改善收敛速度利用了文[1,2]中之算子附加修正方法。方程中之流向导数利用了四阶精度的紧致差分。这是由于该方向上流动参量光滑和边界条件易于处理的缘故。在法向方向上利用了中心差分。对差分方程之隐式部分利用了近似因式分解法。由于选用了文[3]中之特殊Jacobiall系数矩阵分裂法而使计算工作量大为减少。用这种方法计算了超音速粘性绕流问题。计算结果与实验结果相符很好。     

用隐式多重网格法计算三维粘性流动

求解高雷诺湍流流动时,边界层法向网格间距较流向和展向相比非常小,因此边界层中存在高度伸缩的网络,这将大大降低多重网格的求解效率。通过谐调的处理多重网格过程的各个细节,既提高了解的精度,又克服了网格展弦比的影响,使得计算效率提高了6～7倍。     

一种求解S形进气道三维粘性流场的多重网格算法

提出了一种基于多重网格技术,可快速求解Ｓ形进气道粘性内流场的实用工程算法：从Ｂｒａｎｄｔ的多重网格理论出发,采用隐式Ｂｅａｍ－Ｗａｒｍｉｎｇ时间推进格式,结合“锯齿式”迭代技术交替在粗、细两种网格上对雷诺平均三维非定常Ｎ－Ｓ方程进行求解。用此算法数值模拟了某型号Ｓ形进气道的内流场,与试验结果对比显示：计算得到的流场符合一般Ｓ形管道流动特性,所捕获的二次流图谱清晰合理,进气道出口截面处总压恢复系数与试验值较为接近。与普通单层网格计算方法相比,计算结果精度高,收敛速度快。     

可压缩粘性流绕翼型的N─S方程自适应网格解

从Ｎ－Ｓ方程出发,采用ＬＵ－ＡＤＩ隐式分解方法求解绕翼型的可压缩粘性流动。湍流模式采用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模式和ｑ－ω二方程微分模式以研究湍流边界层的非平衡效应;为更好地捕捉激波,对网格进行了自适应调整。通过对ＮＡＣＡ００１２和ＲＡＥ２８２２两种翼型在亚音速和跨音速不同马赫数和雷诺数下进行大量计算,表明该方法对改进湍流计算、提高激波分辨率有较好的效果     

考虑粘性相互作用的跨音速机翼分析

本文将一些已有的三维全位势方程计算方法和一种三维积分边界层方法作了弱相互作用的耦合计算。计算了三种机翼并与已有的试验数据作了比较。