论文部分内容阅读
Rayleigh-Benard(RB)热对流是在一个封闭腔体中,下壁加热,上壁冷却,四周壁面绝热,在上下温差驱动下形成的流动系统,其边界条件简单,但传热系数努塞尔数(Nu)和运动强度雷诺数(Re)与代表驱动力的上下板温差瑞利数(Ra)、流体物性普朗特数(Pr)和宽高比(Γ)的关系十分复杂。研究该系统的对流传热机制对环境、大气、地球物理等人类社会活动有重要的科学价值。长久以来,RB热对流研究沿用的是本世纪初Grossmann和Lohse建立的描述全局物理量关系的理论(GL理论),该理论将RB系统简化为边界层(boundary layer)和中心流动(bulk region),推出全局物理量的关系式,因此无法精确刻画复杂因素对全局换热系数的影响。本文应用佘振苏教授近年来提出的结构系综理论,结合同伦分析和壁射流相关理论,对湍流RB热对流各区域的流动结构耦合机制及其对整体热流的贡献开展了定量研究并构建了各流动区域的二维自相似模型。本文首先获得了大尺度环流以及角涡的自相似多层结构模型。完成了 RB热对流三维(Ra=1 × 107~5 × 109,Pr=0.08~50)和二维(Ra=5 × 107~1 × 1010,Pr=0.01~104)直接数值模拟,基于流动物理特征将流场划分为“大尺度环流”、“角涡”、“斜射流再附区”、“逆压剪切区”、“羽流发射区”等五个区域。基于数值模拟流动的几何相似性,对中等Ra数与Pr数情况下的角涡和大尺度环流,构建了同伦变换的几何相似变量,结合流场时均流函数,提出了角涡及大尺度环流的运动相似解,基于所得参数唯一地定义了流动的特征雷诺数。基于结构系综理论,构建了大尺度结构在近壁区的多层结构函数。从而,建立了可刻画任意二维大尺度涡结构从中心到壁面的完整相似解。此外,将Castaing等人[1]提出的“混合区理论”延拓到角涡流动,在给定角涡尺度标度rcr~Ra0.085情况下,精确刻画了角涡特征雷诺数ReCr~Ra0.25与传热系数Nucr~Ra1/3关于Ra数的标度律关系。进一步,提出了大尺度环流诱发斜射流再附壁面的“壁射流”机制并获得了流动自相似解。发现了以壁射流局部动量率作为近壁流动的特征量归一化最大速度与特征高度所满足的自相似标度关系式。通过分离变量法,推导出壁射流传热系数指数衰减律Nuimp=Numaxexp(-x#)。通过热量输运的动力学平衡以及角涡的标度关系,获得了 Nu数与Ra数的标度关系,Numax~Ra0.2925。基于结构系综理论的对称性分析,给出了逆压梯度边界层的速度和温度剖面随流向变化的解析表达式和羽流发射区的温度分布。基于湍流普朗特数流向与Ra数不变性假设,导出了热流与摩阻以及动量和热量输运涡尺度比值的关系式。根据羽流发射的平衡机制,确认局部传热系数标度律Nu~Ra0.369。由区域空间占比的加权平均局部热流准确计算出全局Nu数与Ra数的定量关系,继而成功推广到其他几何工况下的对流传热过程。最后,应用结构系综理论刻画速度和温度边界层的Pr数效应。研究确认了速度边界层的应力长序函数的两层结构与Pr数无关,还发现随着Pr数减小,压力梯度的作用将逐渐减弱,特征涡尺度减小而粘性底层厚度在增厚。对于温度边界层,应力长满足三层结构,且随着Pr数增大,导热底层以及温度缓冲层厚度增加。应用结构系综理论的参数分析方法,获得了热卡门系数的流向变化规律以及Pr数效应,其结果与计算数相吻合。综上所述,本研究将结构系综理论推广到有热流存在的湍流热对流并提出了二维大尺度分离流动和传热的相似解,给出精确刻画湍流热对流局部平均速度和热流的相似理论。