边界层转换阶段直接数值模拟的影响矩阵法

提出一种求解时间相关、三维不可压完全Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程的影响矩阵法。考虑流动区域为上半空间的平极边界层流动,在水平方向,利用周期性条件,采用Ｆｏｕｒｉｅｒ离散,在半无穷方向,定义了一种特殊映像并利用Ｃｈｅｂｙｓｈｅｖ多项式展开。由此将Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程转化为以水平波数（ｋ１,ｋ２）为参数的一组法向Ｈｅｌｍｈｏｌｔｚ方程组,同时可更有效地求解压力。本文讨论了连续问题的影响矩阵法与原问题的等价性。影响矩阵法的离散形式,以及影响矩阵法的算法实现等。最后给出了一维Ｈｅｌｍｈｏｌｔｙ问题的数值实验结果。     

战术导弹零升阻力雷诺数效应及修正方法

叙述了雷诺数对战术导弹零升阻力的影响,并给出了适合战术导弹的零升阻力系数雷诺数效应修正方法,即变雷诺数试验外推修正方法及工程计算方法。修正结果表明,修正方法是可行的,变雷诺数试验外推法得到的修正量比工程计算得到的修正量更为合理,对于外形简单的战术导弹工程计算仍有较好的精确度。     

跨声速和超声速流中激波/边界层干扰数值模拟

对SST湍流模型中的Bradshaw常数a1进行了修正,并对跨声速和超声速流中激波/边界层干扰进行了数值模拟研究,空间离散采用二阶精度差值的低耗散通量分裂格式(LDFSS),时间离散采用对称高斯-赛德尔(SGS)算法。结果表明:在跨声速流动中,计算得到的壁面压力分布、分离区长度和速度剖面都与实验值吻合较好,而且很好地模拟了典型的λ激波结构;在超声速流动中,修正后模型的计算精度较原始模型有了较大改善,计算得到壁面压力分布和分离点的位置都和实验值吻合较好。     

电弧加热流场湍流度对尖锥边界层转捩影响的研究

电弧加热流场的热环境特性直接影响热防护系统的地面试验数据,由于电弧加热器高温气流和参数波动的原因,直接测量湍流度非常困难。为研究电弧加热流场湍流度对于边界层转捩的影响,采用红外热图热像仪,在电弧加热流场中进行了5°尖锥模型边界层转捩研究。结合数值计算,将试验结果与常规风洞的尖锥边界层转捩结果进行了比较。结果表明：马赫数影响的雷诺数转捩判别准则可以用于计算电弧加热流场的转捩雷诺数;电弧加热流场的尖锥边界层转捩雷诺数显著小于常规风洞的转捩雷诺数,表明在该试验条件下,电弧加热流场的湍流度显著大于常规风洞。     

脊状表面的加工成型技术及其减阻试验研究

为了研究脊状表面的减阻特性,采用数控加工技术在有机玻璃平板上制备脊状表面.在制备过程中,成功解决了脊状结构尺寸的确定、材料的选择、表面粗糙度的提高等问题.利用热线风速仪在低速风洞中分别对光滑表面和多个脊状表面进行了流场测试,获得了湍流边界层内流动参数的分布情况.试验结果表明,脊状表面改变了边界层内的流场结构,抑制了湍流...     

小扰动在二维超声速边界层中的弱非线性传播

采用非线性抛物化稳定性方程(PSE)研究了非平行边界层的弱非线性稳定性。通过在流向采用一阶向后差分,法向采用四阶中心差分,并用预估 校正迭代来耦合非线性项,发展出了求解非线性PSE的数值方法。研究了不同幅值扰动及其高倍谐波的演化情况,发现随着基频扰动幅值的增大,高倍谐波可能失稳,而且失稳的位置会随着基频扰动的增大而向上游移动。通过观察涡量的发展可以发现涡破碎现象。算例与Navier-Stokes方程的直接数值模拟（DNS）结果作了比较,初步检验了其正确性。     

高超声速边界层的转捩问题

取40km高空处的气体参数,对来流马赫数为8、10和12等三种不同情况的平板边界层和楔角为20°的楔型体边界层,做了转捩发生位置的预测。所用方法为改进的eN方法。结果发现,壁面条件对是第一还是第二模态波决定转捩位置有很大关系。如对于传统的eN方法,如果取N为10作为转捩判断标准,对马赫数为12的楔体边界层,等温壁条件下由第二模态波决定的转捩位置距前缘13m左右,而平板边界层中两模态波所决定的转捩位置距前缘均超过50m。而采用改进的eN方法,平板边界层中两模态波所决定的转捩位置距前缘均超过10m,而对马赫数为12的楔体边界层,两模态波所决定的转捩位置距前缘均出现小于5m的情况。     

用PIV实验研究边界层内结构函数标度律

湍流标度律的变化情况显示了湍流边界层内的流动状况。利用高分辨率、高帧率的粒子图像测速技术（PIV）对平板湍流边界层进行了实验研究。研究结果显示,当脉动速度增量结构函数的阶数为正整数且逐渐增大时,由实验所得的脉动速度增量结构函数的标度指数逐渐偏离K41和SL标度指数,出现奇异标度律;而当结构函数的阶数为正分数时,实验所得的标度指数和K41及SL标度指数较为接近;当结构函数的阶数由0趋向于-1时,实验的标度指数明显偏离K41和SL标度指数,出现奇异标度律。这也说明了用低阶结构函数的标度指数也能说明标度律的奇异性及边界层内存在间歇性。对于法向速度增量和涡量分量来说,不管结构函数的高阶或低阶的标度指数都不随壁面位置的不同而改变,具有普适性。     

多段翼型地效风洞实验的地板边界层控制研究

采用有限体积法求解雷诺平均N-S方程和SSTk-ω两方程湍流模型,数值模拟二维多段翼型在不同风洞地板边界条件下流动,并对固定地板边界层吹吸气控制技术进行了数值研究。计算结果分析了地板边界条件对多段翼型气动力、表面压力分布及流动结构的影响,证明了吹吸气控制技术可以有效减小地板边界层厚度,并对均匀吸气与切向吹气两种方法的特点进行了比较。     

MHD控制激波诱导湍流边界层分离的机理分析

为了研究磁流体动力学(Magnetohydrodynamics:MHD)加速边界层对激波-湍流边界层相互作用的影响,用高阶有限差分法求解了小磁雷诺数近似的MHD湍流方程。其中,无粘通量采用WENN格式离散、粘性通量采用Roe平均中心差分离散,时间采用半隐式推进,并采取追赶法求解。计算给出了湍流、电场、磁场和电导率等参数对边界层分离的影响,数值结果显示:在同样的逆压梯度下,湍流边界层分离能更快地趋于稳态流场,且分离区比层流小;通过施加洛仑兹力加速,边界层速度型面变得更加饱满、位移厚度减小、分离点和再附点向激波与固壁的交点靠近,分离区尺寸减小甚至最终被消除。