论文部分内容阅读
对液态金属低普朗特数的流体,使用RANS湍流模型很难精确模拟湍流热扩散项。在本文研究中,使用先进的大涡模拟方法(LES),精细的网格,周期性的边界条件,数值模拟液态金属在三角形排布的堆芯棒束内(又称为子通道)充分发展状态的湍流换热。本文中,首先模拟Reτ=400,600,1400,P/D=1.2以及Reτ=600,P/D=1.05的算例,研究Pr=0.026(铅铋)的液态金属的湍流传热,接着研究不同的分子普朗特数流体(Pr=0.01,0.026,0.1,0.71)在相同流场(Reτ=600,P/D=1.2)下的湍流传热特性。该LES方法是在开源代码OpenFOAM2.3.0的基础上实现。文中分析了时间平均统计量,比如平均温度、温度脉动均方根(RMS)以及湍流热通量等,结果表明,随着雷诺数的增加,温度脉动RMS增大,而且其在子通道内的分布变得复杂。在近壁面很大的区域内,平均温度具有线性规律,同时温度脉动RMS恒定。随着Pr的增加,平均温度的线性区变小,温度脉动RMS曲线峰值的位置向壁面移动,温度脉动强度增强,这意味着温度子层的厚度变薄。液态金属在窄间隙的子通道(P/D=1.05)和宽间隙的子通道(P/D=1.2)的湍流传热特性存在很大的不同,其中,在P/D=1.05时,温度脉动大幅度增加,壁面温度不均匀性增加,增大了壁面材料的热疲劳效应,因此,液态金属的窄间隙堆芯棒束设计,需要特别注意温度的不均匀现象。湍流热通量随着Pr的改变会发生明显变化。通过比较LES计算的Nu与经验公式的计算结果,推荐使用Ushakov公式。本文的计算结果,是作者对堆芯子通道液态金属湍流传热特性研究的开始,在未来将会做更多有关液态金属湍流瞬态场以及湍流相干结构的研究。此外还数值模拟了聚变堆包层第一壁矩形子通道内的湍流对流换热,重点分析了非均匀条件下的湍流换热系数。