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

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

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

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

含激波的湍流流动高精度大涡数值模拟方法

针对采用亚格子模型进行含激波的湍流流动模拟时会面临激波附近的精度损失问题,考虑从通过亚格子模型以及数值模拟方法两方面的改进来实现湍流流动大涡模拟的精度提高.大涡模拟采用了Yee及Sj(o)green (2009)提出的高阶低耗散方法.该方法采用自适应的流场探测器以控制计算中所需区域的数值耗散,并考虑对动力学模型采用在激波位置使用Sagaut和Germano(2005)提出的单边亚格子过滤器和(或)直接禁用亚格子项等方法加以改进.对于标准的马赫数1.5和3条件下的激波-湍流干扰问题,上述新方法相较于全区域采用亚格子模型的方法均表现出了相似的精度提升.同时实现的数值精度改进方案采用了Harten的亚单元分辨过程来定位和锐化激波,并在精确激波位置附近的网格点处采用了单边测试滤波.     

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

本文采用显式 TVD 格式,求解薄层 Navier-Stokes 方程,计算了激波与楔面相互干扰的流场情况,得出了马赫反射的各种等强度线及固壁面上详细的物理量分布,根据计算得出的流场,计算了壁面摩擦阻力及热流分布。     

一类逆风ENO线方法的构造及其对Euler方程的应用

本文利用ENO(essentially-non-oscillatory)插值思想,对非线性双曲型守恒律=0提出了一类逆风ENO线方法,使用这种格式其计算量比C.Shu在文[6]中提出的通量分裂方法明显减少,进而同Runge-kutta时间离散相结合,提出了时空高精度差分格式。本文还把这类方法推广到方程组并用于计算Euler方程组的黎曼问题和拟一维喷管问题,得到了比较满意的数值结果。     

HLL-HLLC格式的构造与应用研究

HLLC格式是一种单调高分辨率格式,能够精确捕捉激波、接触间断和稀疏波,在可压缩流动中具有很高的应用价值;HLL格式相较于HLLC格式,抹平了接触间断,具有较大的数值耗散。然而,数值计算表明,在低马赫数和跨声速计算中表现较好的HLLC格式,在高马赫数强激波附近出现了激波不稳定现象。本文意图通过研究HLLC和HLL格式的数学性质,构造出一种适合更大马赫数范围的HLL-HLLC格式。新格式在较低马赫数下表现出HLLC的性质,是一种低耗散的格式;在高马赫数时具有HLL格式的性质,能够克服激波不稳定现象。通过对高超声速双楔流动、超声速后台阶流动和高超声速钝头体流动数值模拟证明了本文构造格式克服激波不稳定现象的有效性和鲁棒性。     

广义Burgers方程声爆传播模型高阶格式离散

声爆预测技术是研制新一代环保型超声速民机的关键技术之一。针对当前基于广义Burgers方程的声爆远场传播模型的数值求解格式精度较低的情况,开展了高阶格式离散方法研究。通过分析该模型中各效应项的物理性质,应用合适的高阶精度格式对各项分别进行离散求解,并采用经典几何声学射线法计算声爆传播路径,实现了对地面声爆波形的精准预测。通过对美国F-5E机型的声爆飞行试验和第二届国际声爆预测研讨会的典型算例进行数值模拟,验证了本文高阶离散方法的可靠性。进一步的结果分析表明,采用高阶离散格式的预测方法相比传统二阶精度预测方法具有明显的低耗散特性,在同等网格数目上可获得更高分辨率的计算结果,且其网格收敛性远优于传统二阶精度格式,同时在计算效率上也有一定优势;在Burgers方程的各效应项中,非线性项的影响更为明显,因而采用高阶离散格式时其优势更为突出;热黏性吸收项对数值结果影响较小,采用高阶离散格式对计算精度的提升效果也不明显,实际计算中仍可以采用传统离散格式,甚至忽略该项的贡献。     

一种通量加权型紧致格式

将通量加权的思想引入到紧致格式中,构造了一个传统方法与紧致格式混合组成的加权型差分格式.利用该格式与二阶TVD格式分别计算了一维方波、组合波问题和一些黎曼问题,如Sod问题和Shu问题,以及一维定常激波问题.计算结果的比较表明加权格式无论在捕捉各种间断,还是在分辨各种复杂波系上,都具有较大的优势,并与精确解非常吻合.      

超燃进气道分离流动的凸起物控制机理研究

双模态超燃冲压发动机由于压力扰动可能发生不起动现象,造成推力严重下降,对飞行稳定性与飞行安全具有很强的破坏性。不起动初始阶段主要受到激波与边界层相互作用引起的流动分离影响,采用5阶特征型WE-NO格式与3阶TVD型Runge-Kutta格式的高精度数值方法,求解三维Navier-Stokes方程,研究了利用凸起物作为涡发生器的被动控制方法,及其对高超声速流动分离现象的控制效果。结果表明高精度数值格式能够捕捉到复杂精细的流动分离结构,总结了摩阻、压力等在分离再附位置的变化规律;发现凸起物通过诱导形成局部流向涡进行流动控制,能够改变压力分布,减弱分离强度,影响分离结构。     

可压缩多介质流动问题的高精度数值模拟方法

在可压缩流动问题的数值模拟领域,激波的高分辨率计算已取得重要进展.但是包含物质界面的可压缩多介质流动的数值模拟还存在诸多数值挑战,主要表现为界面处数值耗散过大和非物理振荡等问题.界面处流体性质的不连续性是造成可压缩多介质流动问题物理建模与数值方法困难的主要原因.为了建立一套高效的可压缩多介质流动问题的高精度数值模拟方案...