飞船外形高超音速流场N-S方程数值解

用二阶迎风TVD格式求解N-S访方程,对简化的飞船外形高超音速粘性绕流进行了数值模拟。给出了物面压力分布及流场结构,与实验及其他计算结果符合较好。计算表明,肩部圆弧半径对后体流场影响较小,但严重影响肩部附近热流分布。当肩部圆弧半径是头部半径的20％时,肩部热流低干驻点处热流;半径为5％时,肩部热流超过驻点处热流30％。     

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

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

用解析离散法构造 WENO-FCT 格式

本文利用解析离散法的导数计算公式构造高精度ＷＥＮＯ格式,并引入ＦＣＴ通量修正技术以消除非物理解振荡。通过激波－边界层干扰,激波－旋涡相干结构的数值模拟证明方法既能够有效地捕捉激波间断,又对边界层分离及旋涡结构有很好的分辨率,是一类有效的高精度计算格式。     

考虑分子转动自由度的离散统一气体动理学格式

离散统一气体动理学格式（DUGKS）是一种适用于全流域模拟的有限体积方法。之前的研究考虑了分子平动自由度,验证了DUGKS在多尺度问题中的准确性及稳定性。本文基于Rykov模型方程构造了离散统一气体动理学格式,并采用Landau-Teller-Jeans转动能量松弛模型,可用于双原子气体从连续流动到稀薄流动的多尺度问题计算。测试了激波结构、超声速平板绕流以及超声速圆柱绕流等非平衡流动问题,计算结果显示出双原子气体分子中存在平动自由度与转动自由度对应的能量交换过程,并与统一气体动理学格式（UGKS）、直接蒙特卡罗（DSMC）方法的解以及实验值吻合较好。     

运动激波绕射过程中的粘性效应

应用 TVD格式求解二维非定常全 N- S方程 ,以研究运动激波绕射现象。考虑粘性后 ,绕凹角流动时 ,压缩拐角附近表面压力分布与实验结果更为符合 ;绕凸起物时 ,出现反射激波因边界层干扰根部加宽 ,前、后角点附近有局部分离与再附 ,并引起局部压力上升 (与无粘流场模拟相比 )等现象。结果表明 ,全 N- S方程可较好地模拟实际粘性流中复杂激波障碍物干扰流场特性     

滚转运动对乘波飞行器气动特性的影响

 乘波飞行器运动过程中的非定常气动特性是高超声速飞行中的重要物理问题之一。采用数值模拟方法模拟了乘波飞行器在固定迎角下绕其对称轴强迫滚转运动这一过程。比较了在不同频率和滚转角下乘波飞行器的气动特性。计算格式采用AUSM类格式中最新的AUSM+up格式。计算结果表明：AUSM+up能很好地模拟飞行器滚转运动这一非定常过程;滚转运动时,所设计的乘波飞行器能使高压气体很好地附着在乘波飞行器下表面从而使其具有较好的气动特性;当频率较大时,乘波飞行器由于角速度的诱导作用会导致升力出现迟滞现象;做滚转运动时,滚转力矩小于零,产生正阻尼,乘波飞行器不会产生“摇滚”运动。     

基于非结构嵌套网格和逆风格式的旋翼悬停流场数值模拟

为更好地计及桨叶间相互影响,同时捕捉桨尖附近的流场细节,发展了一套非结构运动嵌套网格方法,并针对关键的贡献单元搜寻问题提出了一种高效的搜索方法.为了可方便地扩展到前飞状态的计算,同时为提高旋翼流场涡量的捕捉精度和效率,采用双时间法求解惯性坐标系下的非定常N-S主控方程,空间方向上使用二阶迎风格式,在此基础上,将不同精度格式的计算结果与试验值进行对比,并讨论了桨尖形状变化对旋翼气动特性的影响,结果表明本文建立的CFD方法对旋翼悬停流场的模拟是很有效的.     

航空记录器的过去和未来

航空记录器自出现以来的50年间,为事故/事故征候分析、故障诊断、视情维修、飞行品质监控、试飞监控等提供了丰富的信息。近年来,随着飞机系统不断改进,飞行数据记录器(FlightDataRecorder,FDR)记录的参数己由过去少量的模拟参数发展到现在大量以数字格式记录的参数,驾驶舱话音记录器(CockpitVoiceRecorder,CVR)的记录时间也由原来0.5小时延长到2小时,而且记录器的使用功能也在发生着变化。正确了解、认识和维护航空记录器,对最大限度地利用记录器记录的信息解决飞机使用和维护中出现的技术问题是十分重要的。     

收敛-扩张喷管中运用激波矢量控制法的计算研究

采用二阶TVD格式的有限体积法耦合RNG k-ε湍流模型和非平衡壁面函数求解二维守恒型雷诺平均N-S方程，数值模拟了运用激波矢量控制法的二维收敛一扩张喷管中的流动现象，很好的模拟出了由于次流喷射所产生的激波及上下壁面的压力分布，并将计算结果与国外的实验结果进行了对比。     

ALE方法求解圆柱的涡致振动

数值求解基于拉格朗日-欧拉(ALE)描述的不可压缩流体的N-S方程,计算了较小雷诺数下圆柱的涡致振动现象.N-S方程的对流项和扩散项分别采用三阶迎风紧致格式和四阶中心紧致格式离散,计算网格采用H-O非交错网格系统,并结合分块耦合方法.柱体运动简化为弹簧-阻尼-质量系统,柱体运动方程采用经典龙格-库塔方法求解.通过模拟柱体和流体之间的非线性耦合作用,成功地捕捉到了"锁定"、"拍"和"相位开关"等现象,并与试验数据相吻合.另外,本文详细分析了Re＝200时圆柱尾涡形态、升力、阻力以及圆柱位移等随圆柱自振频率变化的过程,捕捉到尾涡结构变化的频率转折点.