超临界压力下对传热的湍流模型研究.pdfVIP

A Dissertation Submitted to Shanghai Jiao Tong University for the Degree of Master STUDY ON TURBULENCE MODEL UNDER SUPERCRITICAL PRESSURE Author ：Zhang Hao Specialty: Nuclear Science and Engineering Advisor: Prof. Yang Yan-Hua School of Mechanical Engineering Shanghai Jiao Tong University Shanghai, P. R. China February, 2010 超临界压力下对流传热的湍流模型研究 摘 要 超临界压水堆是第四代堆型之一，极具发展前景。顾名思义，超 临界压水堆是在超临界压力下运行，因此，研究超临界压力下流体的 对流传热特性具有重要的意义。 在超临界压力下，物性随温度发生剧烈变化，特别是在拟临界点 温度附近。物性的剧烈变化对传热的影响很大。根据传热机理的不同， 可分为以下几类情况：大流量、小热流密度下的强迫对流，比热容的 变化引起传热强化；小流量、大热流密度下的混合对流，密度随温度 的变化导致产生强大的浮升力，对传热产生很大的影响；大流量、大 热流密度下的对流传热，由于密度的变化，产生强大的加速效应，从 而使传热能力变弱。 作者编写了 CFD 程序，通过数值模拟的方法研究了竖直通道中超 临界流体强迫对流和混合对流的传热特性。因为超临界流体可视为单 相流体，所以数值计算工作的主要内容是对湍流模型的评价以及改 进。低 Reynolds 数 k-ε模型是非常具有发展前景的湍流模型。本文评 价了多个低 Reynolds 数 k-ε模型的表现，并根据浮升力作用机理，对 湍流模型进行了一系列改进。 对于强迫对流传热，计算结果表明传统的低 Reynolds 数 k-ε模型 能很准确的模拟壁面温度。通过对数值计算结果的分析，作者认为其 传热特性是由近壁面处 Prandtl 数的变化所决定。近壁面处的 Prandtl 数变大，温度边界层的粘性底层变薄。当 Prandtl 数达到最大值时， 温度边界层的粘性底层变得最薄，传热能力达到最强，即传热强化现 象。 第 I 页 根据浮升力作用方向的不同，混合对流可分为浮升力助推流动和 浮升力阻滞流动。对于浮升力助推流，在浮升力的作用下，传热能力 依次经过弱化、回复、再弱化。传统的湍流模型虽然能模拟出壁面温 度的变化趋势，但是无法准确地计算壁面温度。DNS 数据显示，此 时k 方程中的浮升力项占据很重要的地位，而传统的浮升力项模型极 大地低估了浮升力的作用，导致其无法准确地计算壁面温度。作者研 究了浮升力项模型，并对湍流模型进行了改进，提出：k 方程中的浮 升力项采用改进的 AFM 模型，湍流模型采用 k-ε-k -ε 模型。计算结果 t t 表明作者改进后的模型极大地改进了计算结果，能较准确地计算壁面 温度。理论分析以及计算结果表明，浮升力作用决定了其传热特性。 一方面浮升力对湍流的间接作用（外部效应）使剪切应力和速度梯度 发生变化，当浮升力的作用使剪切应力和速度梯度变小，传热能力变 弱，当剪切应力的绝对值开始变大以及 M 型的速度分布出现时，传 热能力回复。另一方面浮升力对湍流的直接作用（结构效应）也影响 着传热，k 方程中的浮升力项起着很重要的作用，表明浮升力作用是 湍动能产生的主要来源之一。 对于浮升力阻滞流，浮升力的作用使传热能力强化。理论分析和 DNS 数据表明，k 方程中的浮升力项不占重要地位，因此浮