轴对称跨声速绕流的凸角流动分析

本文以锥-柱-船尾形旋成俸跨声速绕流为背景,讨论了轴对称踌声速绕流的凸角流动。给出了该流动特征、规律、并着重与相应平面二维Prandtl-Meyer型膨胀流动作了对比,指出了二者异同之处。     

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

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

NACA4412翼型低速绕流数值计算中湍流模型对比

使用Spalart-Allmaras(S-A)、SSTk-ω、Gao-Yong 3种湍流模型对NACA4412翼型低速绕流进行了数值计算,研究了尾迹流动松弛效应。对流项采用Roe格式离散,扩散项采用二阶中心格式离散,离散后的控制方程用多步Runge-Kutta显示时间推进法求解。计算中对翼型尾缘流松弛效应进行了分析,比较了翼型表面压力系数、速度剖面、雷诺应力等的分布,3种湍流模型总体上能够较好地模拟NACA4412翼型低速绕流。SSTk-ω模型对流动细节及升力系数计算最好,Gao-Yong模型对翼型平均速度剖面及雷诺剪切应力分布计算最准确。     

航天飞机头部无粘绕流的无波动,无自由参数格式的数值计算

本文采用气动中心张涵信教授提出的无波动、无自由参数耗散格式(NND)进行了如下的计算:航天飞机头部段对称剖面内的二维无粘绕流;航天飞机头部段三维无粘绕流;轴对称凹陷外形的三维无粘绕流。在物面坐标系下,用代数方程生成计算网格,利用特征相容关系式处理激波和物面边界。为了比较,还用SCM格式进行了计算。结果表明,本文所采用的NND格式,在捕捉流场内激波时不会出现非物理振荡,其收敛速度和所获得的计算结果都是令人满意的。     

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

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

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

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

高超声速电离绕流的数值模拟

本文从三维N－S方程出发，采用稳式L－U分解算法和7组分15个反应的化学模型，数值模拟高超声速电离空气绕流。首先采用对称TVD格式、AUSMPW＋格式和Van Leer的矢通量分裂格式计算了高超声速球头绕流，并对它们的计算结果做了对比分析。然后用前两种格式，对RAM－C飞行试验模型三个再入高度（81km、71km、61km）的流场进行了数值模拟，计算的流场电子数密度值和试验测量数据符合较好。     

绕流问题的区域收缩法及其数值解

本文提出了一种适于绕流问题的数值方法,它考虑了流场中不同区域内不同的物理现象,给出了外部边界条件的良好逼近。使用这种方法可以节省大量的计算机存贮单元和计算时间,也可以推广到三维绕流问题的数值计算。     

应用于翼型绕流的线性/非线性湍流模式的研究

本文选取了四个线性湍流模式、四个非线性涡粘性湍流模式和一个显式代数应力模式对绕翼型的不可压缩分离流动进行了数值模拟.因计算鲁棒性的需要,其中部分模式在壁面附近耦合了一方程模式.通过与实验结果的比较,对翼型在大攻角情况下流动产生分离的气动特性进行了评估.计算结果表明,非线性模式能够较好地反映湍流的各向异性和曲率影响.     

汽车绕流的数值分析及计算机模拟

本文简要介绍了汽车绕流数值分析的有关理论基础和计算机模拟几个重要问题,涉及基本方程、有限体积法的应用、交错网格系统及边界条件的处理等,并给出了应用算例。对有关计算的相应技术处理也进行了探讨。     

钝体绕流流场分离涡的数值模拟

 采用近似因式分解和特殊Jacobin系数矩阵分裂法,求解可压缩Navier-Stokes方程,数值模拟矩形柱体绕流流场中分离涡的形成和发展。给出了Ma_∞=0.3,0.8;Re=10~4,10~6的计算结果。从结果中可清楚地看到柱体非对称涡形成的非定常过程。     

钝体化学非平衡粘性流数值模拟

用有限差分法求解Ｎ Ｓ方程,数值模拟了二维钝楔头部和近尾迹非平衡化学反应流场。模拟条件为Ｍａ∞＝２５,飞行高度ｈ＝７０ｋｍ,半楔角θ＝２０°,钝度比为０６０;气体组分为７个组元,共１１个化学反应式。头部流场结果与已有计算结果符合较好。在Ｍａ∞＝２５条件下,底部近尾迹较大区域内具有较高的温度和较大的电子密度,且下降趋势缓慢。     

可压缩后台阶流的数值模拟

基于可压缩Navier-Stokes方程,采用Beam-Warming矢通量分裂格式以及vanLeer的矢通量分裂方法,分别通过Johnson-King及Baldwin-Lomax湍流模型对亚音速和超音速的后台阶流进行了计算,计算结果与实验符合较好。     

绕翼型低速湍流的数值模拟

采用求解低速流动的SIMPLE算法,用标准K-ε两方程模型、非线性K-ε两方程模型和一种改进的K-ε两方程模型来对绕两种翼型（A翼型和GAW-1翼型）在接近失速攻角情况下的低速湍流流动进行了数值模拟。计算结果表明,非线性K-ε两方程模型和改进的K-ε两方程模型较好地模拟出了翼型表面上的分离流动特性,较准确地预测出了分离点位置以及翼型的气动力系数,对翼型与机翼的工程设计具有一定的参考价值。     

二维埋入式进气道的数值分析

采用并改进了适用于分离流动的 Ｊｏｈｎｓｏｎ Ｋｉｎｇ 紊流模型,结合 Ｌ ＵＳ Ｇ Ｓ近似隐式格式和基于 Ｍ Ｕ Ｓ Ｃ Ｌ非线性插值的 Ｒｏｅ 格式编写了雷诺平均的二维 Ｎ-Ｓ方程的求解程序。然后借助于该流场计算程序完成了二维埋入式进气道的数值研究,分析了几何设计参数对进气道性能的影响。     

翼型大攻角绕流的数值模拟

以求解二维N-S方程数值模拟NACA0012翼型大攻角状态的可压绕流特性;N-S方程是在贴体坐标系中给出的,以代数方法生成C型网格系统。采用LU-ADI格式计算,为提高格式的稳定性在隐式和显式部分分别添加了2阶和4阶人工粘性项。应用BaldwinLomax湍流二层代数模型模拟了大攻角时前缘分离涡的形成,旋涡对流及其非定常现象。在某些Mach数和攻角下NACA0012翼型的湍流解具有周期性。通过比较,本文数值计算结果同实验及国外相应的数值计算结果基本吻合。     

矩形风洞收缩段流场的计算和分析

本文结合西工大低速风洞设计进行了矩形三元收缩段流场的计算和分析,计算采用差分方法和贴体坐标,以AF1格式进行离散化并利用ADI方法求解。对三种常用的收缩曲线(即Witozinsky曲线,五次方曲线和双三次方曲线)的三元流场进行了具体计算和比较。计算结果正确反映了三元收缩段流动的特点,给出了流动三元性对壁面逆压梯度和出口均匀性影响的数值结果,表明本文方法是三元收缩段设计的有效分析和设计手段。     

化学非平衡边界层绕流计算

本文介绍用化学非平衡边界层方程组,数值计算高超声速钝头体绕流流场、组元浓度和电子密度分布情况。计算中考虑了电离空气的七种主要成分以及相互间的化学反应,比较了变熵、普朗特数及Lewis数的影响,讨论了完全催化壁和非催化壁两种情况,计算结果与文献及试验进行了比较,结果令人满意。     

三维湍流室内气流数值解

本文基于有限差分法,k-ε模型和p-v法提出了计算三维湍流室内气流问题的数学模型,且用它求解了某水电厂厂房的流场。同时还做了相应的模型试验,并把二者进行了比较,获得了满意的结果。     

燃气舵外围流场计算

对燃气舵沿气流方向的横截面进行其外围的二维流场计算。计算中将燃气舵外围流场划分成三个区域,I区为钝头体头部绕流脱体激波后的亚跨声速区,Ⅱ区为脱体激波后的超声速区,Ⅲ区为膨胀波后的流动区。对各个流动区,分别用条带法、特征线法和普朗特尔-迈耶公式计算,获得各区的流动参数,为燃气舵的气动力、气动加热、烧蚀、传热、应力等计算提供流场数据。