悬停旋翼无粘流场数值模拟中的多重网格方法

发展了一种基于不同空间离散格式的多重网格算法,并应用于悬停旋翼无粘绕流的Euler方程数值模拟。由于悬停旋翼流场中存在不可压区域,同时旋翼尾涡系统的发展需要较长的时间,使得旋翼流场的收敛速度远低于固定翼流场,因此研究旋翼流场的多重网格算法具有重要意义。空间离散采用了Roe s FDS格式和Jameson中心有限体积格式,时间推进应用了五步Runge-Kutta方法。采用多重网格V循环方式,对一跨声速悬停旋翼无粘流场进行了数值计算。计算结果表明:多重网格算法可以显著加速悬停旋翼无粘流场的数值计算收敛速度;无论在激波分辨率还是在计算精度上,Roe s FDS格式都优于JST格式。     

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

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

钝体头部边界层特征谱分析和转捩途径预测

利用线性稳定性理论,分别对钝头体高超声速绕流边界层中的亚声速区域、跨声速区域和超声速区域的特征谱进行了系统分析.结果表明,根据线性稳定性理论,所分析的各速度范围流场,对所有单频波的扰动都是稳定的.但随速度的增加,稳定波的衰减率越来越低;随激波到物面距离的增加,s族扰动波相速度的实部向无量纲速度1靠拢,说明头部主激波的作用在一定意义上与平面槽道流动的上壁面相似;研究工作还表明,流场中存在着低衰减率扰动波,且无论那个速度范围内的流场,其最低衰减率扰动波的相速度实部都在0.7左右.根据钝头体特征谱分析,本文可能发生的"bypass"型转捩机理进行了理论分析,提出了几条可能的"bypass"型转捩途径.     

径向基函数在三维Euler方程数值计算中的应用

基于三维Euler方程,对非结构四面体网格给出了一种基于紧支径向基函数重构的ENO型有限体积格式,方法的主要思想是先对每一个四面体单元构造插值径向基函数,而在计算交界面的流通量采用高斯积分公式以保证格式的整体精度,时间离散采用三阶TVD Runge-Kutta方法。最后用该格式对一些典型算例进行了数值模拟,结果表明该方法计算速度快,对间断有很好的分辨能力。     

燃气自由射流的正格式数值模拟

采用二阶正格式方法对复杂燃气自由射流进行了数值模拟。将二维守恒方程的正格式方法发展到轴对称Euler方程组的求解,并对不同欠膨胀压力比下的燃气射流进行了数值计算。计算结果表明,该方法能够较好捕捉到波胞、滑移线、射流边界以及三波点等复杂射流流场波系结构,与实验照片反映的流动特征及高精度、高分辨率的三阶间断有限元方法的计算结果吻合较好,马赫盘位置与理论预估值和实验测量结果误差较小。证明该方法对激波具有较强的捕捉能力,为此类流动的数值研究提供了一种新手段。     

AUSM+方法在黏性内流模拟中的应用

简要阐述了三维雷诺平均Navier-Stokes(RANS)方程求解程序中使用的对流迎风分裂格式改进形式(AUSM+)以及Spalart-Allmaras(S-A)一方程湍流模型,然后对Sajben收缩扩散通道内跨声速流动和高速透平静叶排内部流动进行了数值模拟,并将数值结果与实验结果进行了比较,结果表明自主开发的三维R...     

可压缩混合层中的涡结构和激波

采用高精度高分辨率格式进行直接数值模拟,既得到了对流马赫数为1.2的混合层中线性发展之后的λ涡、发卡涡和涡的破碎,还得到了破碎后发展的多模态主导下的分层涡结构,以及对应动量厚度饱和下的涡结构。其流动结构的发展显示了混合层在自由空间充分发展能力的有限性。此外捕捉和显示了混合层中自然发展起来的外缘激波、对旋激波和单涡激波,激波主要是由于涡结构的高速旋转造成的相对超声速、低压区和逆压梯度区而产生的。激波的存在抑制了混合层从主流吸取质量和动量,也抑制了混合层内部流质混合的均匀性。     

全隐式TVD格式求解二维Euler方程的对角化算法

本文采用有限体积法,Chakravarthy提出的基于近似Riemann解的全隐式TVD差分格式求解二维非定常Euler方程。采用牛顿法对半离散化后的全隐式差分方程进行线性化,然后对线性化后的方程采用近似因式分解法处理,得到了两组三对角块矩阵方程组。在求解该方程组时,采用了矩阵对角化思想使原来的4×4三对角块矩阵方程组转化为4组分离的三对角代数方程组,因而大大节省了求解方程组所耗费的计算时间,提高了计算效率。通过对翼型的跨声速绕流问题及圆柱超声速流动问题的计算,证实了本文算法不但简便可靠,而且还具有较强的通用性。采用TVD格式捕捉到的激波分辨率高,上、下游无任何波动。     

块结构的MUSCL算法求解三维进气道流场

本文将ＭＵＳＣＬ格式扩展为三维算法,并用于湍流流场的数值模拟。通过求解超声双斜楔压缩面流动,得到了含有激波、嵌入激波、剪切流和湍流边界层等流动特征的复杂流动结构,数值解与实验数据吻合得很好。文中将此算法与块结构网格技术相结合,使之能求解具有三维复杂几何形状的流场。并应用该算法模拟了Ｆ－１６战斗机进气道在５种典型飞行状态下的三维湍流流场,计算结果清楚地揭示了这些情况下的三维流动现象,得出了进气道的性     

一种同位网格上的迎风有限元方法

通过Ｇｒｅｅｎ定理将对流项变量从微分算子中分离出来,从插值函数入手引入迎风格式,是对强对流问题有限元计算中对流项变量的一种新的处理方法。按这种方法采用局部斜迎风格式及速－压同位网格公式,构成了一种对高Ｒｅｙｎｏｌｄｓ数流体流动数值模拟比较有效的有限元方法。