激波干扰及超声速湍流边界层的数值模拟研究

【摘要】：激波和超声速湍流边界层的研究在流动物理和工程应用中都具有重要意义。本文采用直接数值模拟方法研究了激波-激波干扰和超声速湍流边界层的流动特性。文中着重探讨了激波-激波干扰中的非定常流动现象及反馈机制、超声速湍流边界层内流动拓扑的演化和壁面温度对超声速湍流边界层统计特性的影响。主要工作和研究成果如下： (1)研究了高超声速IV型激波-激波干扰的非定常流动特性。基于对涡量和胀压输运方程源项的分析,揭示了流动结构与剪切、压缩及热力过程之间的内在联系。详细地分析了流场中存在的非定常流动现象,包括超声速“喷流”冲击壁面、“喷流”激波的非定常运动、激波-边界层相互作用和旋涡-边界层相互作用等。研究表明：“喷流”激波的振荡和变形会产生“热斑”,导致物体表面热流峰值出现瞬时极大值；“喷流”激波和超声速“喷流”所诱导的剪切层相互作用,形成向上游传播的压缩波；压缩波与上游的弓形激波作用产生向物体传播的熵波和涡波。通过对压力场的本征正交分解发现,流场中的压缩波还与激波-边界层相互作用和旋涡-边界层相互作用相关联。基于对流场演化过程的分析,发现流场非定常运动的反馈机制与超声速“喷流”相关,由向上游传播的压缩波和向下游演化的相干涡结构构成反馈回路。基于自激振荡机制,定量地分析了流动的非定常特征,包括非定常运动的频率以及入射激波角度对运动频率的影响。此外,随着雷诺数的减小,流场由非定常状态变为定常状态。 (2)研究了超声速湍流边界层内流体质点的拓扑结构演化。系统地分析了速度梯度不变量空间中不变量变化率的瞬时和平均统计特性。研究表明：不变量发生剧烈变化的可能性呈现明显的-3幂次律衰减特征；在归一化的不变量空间内,边界层外层的条件平均迹向内盘旋至原点,缓冲层的条件平均迹呈现极限环行为。基于不变量动力学演化方程统计分析了相互作用项、压力-Hessian相关项和粘性项对拓扑结构演化的贡献,揭示了压缩性对流动拓扑演化的影响。统计结果表明：拓扑结构的平均演化受压缩性的影响显著；不变量间的相互作用项决定了拓扑结构在不变量Q-R平面按顺时针方向演化；在局部压缩区和局部膨胀区,压力项的贡献起主导作用,但两个区域内的演化趋势是相反的；粘性项主要表现为耗散作用,导致速度梯度衰减,在粘性底层以外可以采用线性衰减模型来模化。 (3)研究了壁面温度对超声速湍流边界层的湍流统计特性的影响。系统地分析了平均量和湍流脉动量的统计特征、拟序涡结构的空间结构以及应变与涡量的相互作用。结果验证了Morkovin假设的合理性,即壁面温度对湍流平均特征的影响主要通过平均密度和粘性的变化所体现。通过分析联合概率密度分布,本文建立了密度脉动和温度脉动之间、流向速度脉动和温度脉动之间的合理关系式。在边界层内层,拟序涡结构与壁面的夹角随着壁面温度的增加而增大；在边界层外层,拟序涡结构与壁面的夹角基本不受壁面温度影响。通过涡拉伸项定量分析了应变与涡量的相互作用；结果表明,在边界层内层第二主应变率相关项占主导,主要来源于不稳定应变占优区域(UN/S/S)和径向收缩、轴向拉伸的涡结构(SF/S),应变占优的UN/S/S的平均涡拉伸强度大于旋转占优的SF/S;在边界层外层第一主应变率相关项占主导,主要来源于SF/S,并且SF/S的平均涡拉伸强度大于UN/S/S。