烧蚀外形方程差分计算方法研究

分析了烧蚀外形方程的性质,研究了人工粘性项对烧蚀外形计算结果的影响,发现人工粘性在抹平非物理波动的同时,会严重影响计算精度.借鉴于CFD中NND格式的思想,构造了求解烧蚀外形方程的无波动、无自由参数的耗散差分格式,得到了满意的计算结果.     

适于黏性计算的自适应笛卡儿网格生成及其应用

复杂外形的黏性网格生成已经成为计算流体力学(CFD)在工程应用中的主要瓶颈,高效而鲁棒的网格生成技术研究显得尤为迫切。使用基于物面外形自适应加密的笛卡儿网格方法生成计算网格,投影方法拟合壁面边界,将投影得到的柱形单元沿法向分层获得适于黏性计算的网格。应用特征恢复技术恢复凹面特征并改善凸面处的网格质量。流场解算中应用最小二乘法重构获得2阶精度,Roe格式计算无黏通量,Holmes-Connell方法计算黏性通量。加入Spalart-Allmaras一方程湍流模型以获得湍流解。上下对称高斯 赛德尔迭代(LU-SGS)格式隐式时间推进,对网格重排序保证格式的稳定性和平衡性。算例显示网格生成方法是快速的和鲁棒的,流场计算精度满足工程需求,因此该方法非常适合应用于实际工程计算。     

超声速空气介质自由扩散混合流涡的演化

在不同谐波扰动促发方式下,观察二维超声速空气/空气自由扩散混合流中涡的形成及演化.考虑含空气组分扩散效应的可压缩Navier-Stokes方程,对流项采用3阶迎风紧致格式离散,输运项采用6阶对称紧致格式离散,非定常时间推进采用3阶紧致存储显式Runge-Kutta方法.在无相差情形下,获得了大尺度基频涡结构的饱和、一次对并、二次对并、三涡对并等现象.在Ma_c=0.3低对流Mach数下,基频涡较饱满,但流向尺度较小.受空气真实气体特性影响,在Ma_c=0.8高对流Mach数下未发现涡/小激波结构.      

"回收/调节"方法在混合LES/RANS模拟方法中的应用

采用一种混合大涡模拟/雷诺平均Navier-Stokes(LES/RANS)模拟方法结合三阶加权基本无振荡(WENO)格式对马赫数为2.88的压缩斜坡流动进行了数值模拟,并采用"回收/调节"方法生成入口湍流脉动边界条件.当采用固定入口边界条件时,湍流边界层会缺乏合理的湍流能量,使之抵抗逆压梯度的能力减弱,会严重高估分离...     

机动飞行器中等攻角高超声速无粘绕流数值模拟

本文采用ＮＮＤ差分格式，利用推进－迭代方法，数值模拟了机动飞行器的中等攻角高超声速无粘绕流流场，发展了一种适用于处理中等攻下背风面流场的技术，大量计算表明，结果准确、可靠。气动力力系数和壁面压力分布与实验数据一致，流场波系结构模拟正确。     

二维Euler方程SCB格式的并行实现

本文讨论二维双曲型守恒律方程组的高分辨ＳＣＢ格式，并应用这些格式求解二维Ｅｕｌｅｒ方程组。计算过程中在江南－Ⅱ型并行计算机上进行的，计算过程运用了许多并行计算技术，计算结果令人满意。     

高分辨率格式与钝舵绕流数值模拟

 采用Beam-Warming/TVD格式，Roe二阶、三阶以及四阶通量差分分裂格式数值计算雷诺平均Navier-Stokes方程，湍流附加粘性采用Baldwin-Lomax模型计算，针对超声速绕钝圆舵流动进行了数值模拟比较研究，计算给出了流场的详细结构，物面压力分布与实验值符合良好，尤其是Roe三阶格式给出了较好的计算结果，两个超声速回流区以及湍流边界层分离也得到了很好的数值模拟。     

梯度加权局部分析解差分格式及其应用

提出了一种新的针对Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程符合含有梯度项的对流扩散方程的梯度加权局部分析解的有限差分格式，且给出了通用的系数表达式。ＧＷ－ＬＡＤ格式是基于积分型的有限体积法并且通过附加源项的方法来实现。能够适用于粘性和无粘性及不可压缩和可压缩的压力驱动流动。通过对叶栅流道中二维湍流流的计算表明ＧＷＬＡＤ格式与ＥＤ格式相比，具有较高的精度和很快的收敛速度并且得到更好的压力场和速度场。     

TVD格式数值模拟斜腔内旋涡运动

 采用ＴＶＤ有限体积格式,通过求解非定常可压缩Ｎ－Ｓ方程,对恒定出入口压力比驱动下的斜腔内部流动进行了数值模拟。计算后得到清晰的流场涡系结构及其发展过程；显示了空腔流动中的自维持振荡现象；并分析了剪切层不稳定对周期性振荡的作用机制。     

粉尘气体中平面激波在固壁反射的高分辨率数值模拟

本文在稀颗粒群双流体简化模型下，采用算子分裂技术，综合应用高分辨率ＧＲＰ格式，全二阶精度的ＦＣＴ格式及预估校正格式，取三种不同粒相固壁边界条件，数值模拟了粉尘气体激波在刚性固壁上的正反射，确证了由理想平衡混合气体分析所预示的三类激波反射结构，给出了非平衡反射时波后流场的松弛特性，并评价了不同的粒相固壁边界条件假定的适用性。