论文部分内容阅读
本论文利用一般正交曲线坐标系下的抛物化稳定性方程(PSE)和高精度激波装配法,开展了高超声速钝锥边界层稳定性特征的研究工作,主要研究成果及创新如下:
(1)对高精度激波装配法提出了一种改进的并行算法,提高了该方法的并行求解效率,扩大了应用范围,同时构造了七阶紧致迎风有限差分格式用于逼近对流项,以提高方法的分辨率.通过对钝锥绕流问题的验证计算表明,本文的计算结果与他人结果[87,108]吻合的很好,因此,当前的数值方法有效、可靠,能够用于模拟来流中10-5-10-4量级的小扰动演化问题.
(2)发展了基于一般正交曲线坐标系的抛物化稳定性方程及其计算软件,该软件包括线性PSE和非线性PSE两个求解器。与直接数值模拟(DNS)相比,该方法不仅便于对边界层内扰动波的演化特征进行分析,而且降低了至少一个数量级的计算时间。算例表明,该方法的结果与DNS以及线性稳定性理论(LST)的结果吻合甚好,验证了其可靠性.
(3)细致地分析了壁温对球头半径为3.81×10-3米,半锥角为7°的高超声速钝锥边界层稳定性的影响,得到如下有意义的结果:
i)DNS结果表明,壁温变化显著地影响高超声速钝锥边界层内广义速度拐点沿壁面法向的数量和位置;LST分析表明,降低壁温加速第二模态不稳定波的增长,同时抑制了第一模态扰动的增长.向下游推进中,当地最不稳定的第二模态扰动频率逐渐降低,且增长率持续增大;当地最不稳定的第一模态频率先升高,然后又有所降低,且增长率也持续增大.壁温越低中性曲线越靠后,且第一模态推后的比第二模态更严重.由此可见,壁温对高超声速边界层稳定性有显著的影响.
ii)当来流中加入10-4的小扰动时,DNS表明下游290头半径左右的壁面扰动幅值与壁温之间没有单调变化关系,对于冷却壁存一个临界温度T∈(280,375K),当壁温低于Tc后,冷却壁时的扰动幅值大于绝热壁的情况,反之,结论相反。这一结论与Stetson[79]的实验基本一致。
iii)在中性曲线下分支处,壁面压力扰动的感受性系数的分布表明,壁温越低,感受性系数也越小.低频扰动波的感受性系数大约在(0.1,0.2)之间,而高频扰动波的情况大约在(0.7,0.8)之间.
iv)用线性PSE分析了不同壁温情况下小扰动的演化特征,结果表明,壁温变化显著影响流向波数、扰动增长率和扰动幅值的变化特征.在壁面特定位置以前,壁温越低波数越小,当大于这个值以后,结论相反;降低壁温使第二模态增长率增大;且壁温越低最终扰动幅值越大.
v)用非线性PSE模拟边界层中有限振幅扰动演化表明,边界层内扰动演化足够长时间后高次谐波得以快速增长,且降低壁温加速了高次谐波的增长,特别是三维高次谐波具有更大的增长率,高次谐波促进了边界层不稳定的发展.