浮力对竖直细微管内超临界CO2流动与传热的影响论文.pdfVIP

中国工程热物理学会 传熟传质学 第十一届学术会议 编号：053386 浮力对竖直细微管内超临界C02流动与传热 ‘ 的影响 饶政华 廖胜明 (中南大学制’辟空调研究所，长沙410075) 联系电话：(0731)2655054E—mail：smliao@csu．edu．cn 摘要本文数值研究了超临界c吼在竖直细微圆管内流动时，浮升力对流动与传热的影响，得到了管内 层流及湍流分别在加热或冷却条件下．主流向上、向下和零重力流动时的对流换热系数^和努谢尔数胁，井 计算了D佛，’值，研究结果显示：浮升力对竖直细微圆管内超临界c如流动与传热的影响，与流动的状 浮升力的影响可忽略，其上边界值2．o×1旷比Jad岱on提出的106要高。 关键词超临界c02；浮升力；细微管：数值解 1前言 有良好的热物理性质，如导热系数高，运动粘度低，单位容积的制冷量高等，有优良的流 动和传热特性。Co：跨I临界循环中，高压倒换热设备(习惯上称为气体冷却器)一般采用 细微通道管的形式，其管径达到0．79mml2】，其内部的∞2保持为超临界态。研究超临界 c02流动与传热的特性以及优化气体冷却器的设计，对提高系统循环效率有很重要的意义。 et “aoa1．【3】数值研究了竖直细微管内超临界c02层流流动换热的特性，还对水平细微管内 eta1．【5q运用4种2_D湍流模型给出了管内超l临界 的换热系数进行了实验测定【4】；D∞g c02冷却和加热流的数值解，并与他们实验测定的结果进行了比较，由此推荐了符合度较 et 好的湍流模型；Jacksona1．【7】从理论上分析了揣流条件下，大管径管内超临界流体被加 热时浮升力对传热的影响。但关于超临界c02在竖直细微管内流动与传热时，浮升力影响 的研究则很少。 强迫对流换热中，浮升力对流动与传热特性的影响也不可忽视。由于超临界c02强烈的 变物性特点(尤其在虚拟临界点附近)，浮升力的影响可能还非常大。竖直管内的浮升力 会改变管内速度、温度分布，增大或减小对流换热系数，缩短或延长达到热充分发展的时 间等，其影响的结果和程度取决于流动状态、管径大小、浮力与主流的方向等因素。本文 采用数值方法模拟了超临界C02在竖直细微圆管内层流和湍流流动与传热的情况，计算采用 REFPRoP9J。 的物性参数来自NIsT 2物理模型与数学描述 本文考虑超临界c02在管内稳态流动与传热的情况，模型及其坐标系统如图1所示。流 体自上而下、自下而上及零重力地流过冷却(或加热)的竖直圆管，管长1m。本文所指的 采用柱坐标对该物理模型进行数学描述，其通用控制方程为： 1295 旦署+；学=去c‘》+；言∽挚+。 ∽ 层流时，连续性、动量、能量方程中的当量扩散系数碍，源项母具体形式参见文献【3】。对 于湍流，本文利用标准的肛P模型进行计算，模型中的连续性、动量、能量、奴e方程也有 3数值计算方法 本文采用sIMPⅡi算法对上述控制方程组进行数值求解，帕坦卡对该方法做了详细的 描述【“。它是以控制容积为基础的，采用有限容积法将控制方程组离散成差分代数方程组， 再利用Ⅱ}MA迭代法进行求解。为了保证足够的精度，本文采用了2D0x100(轴向×径向) 的网格点，方程迭代收敛的标准为： 障：：：．：二生二J‘10～，≯=“，v，P． (2) 妒“ 4计算结果分析与讨论 层流与湍流的流动传热机理及其数学描述均不同，浮升力影响的结果和程度也不同， 因此本文分别对两种流态的流动传热情况进行了研究，得到