层流下功能流体在不同管径中冷却效果的比较.doc

层流下功能流体在不同管径中冷却效果的比较 靳 健 刘沛清 （北京航空航天大学航空科学与工程学院，北京 100083） 摘要 本文利用两相流模型研究了功能流体在不同管径的管道中流动和传热情况。数值模拟表明，在壁面附近功能流体存在一个低浓度区，而且在直径较小的管中，壁面附近浓度降低的更加明显。研究表明，当热流密度较低时，较大管径中功能流体的冷却性能较好，但当热流密度进一步增加时，较大管径中功能流体的冷却性能明显下降。本文详细分析了这种现象的原因 关键字 微胶囊相变材料； 两相流模型； 有效热导率； 热阻； 1. 引言 功能流体自从上个世纪七十年代末被提出后，受到了人们的广泛关注[1-3]。由于受实验设备复杂性和费用的限制，目前对功能流体的实验研究通常限于单管径的情况。对于模拟研究，虽然将功能流体处理成单相流的方法较为简单，也能够在一定程度上把握相变悬浊液强化传热的特性。但颗粒和液体由于物性的不同，流动状况不会完全相同，在某些工况下差别可能还会较明显，这时利用单相流模型就无法准确的把握相变悬浊液的流动本质。因此本文建立了双流体模型，并针对体积浓度为10%的功能流体，利用该模型计算分析了层流下管径变化对功能流体传热和冷却效果的影响。 2. 理论模型 本文中选用了两种管径的不锈钢直管内径，；，，两种管道的长径比均为100。管壁采用均匀热流密度加热 。在本文中，液体和相变颗粒被分别定义为液相和固相（分别用下标l和s表示），并假定两相是互相贯穿的。以文献[4]中所提出的两相流模型为基础，对两相分别建立连续方程、动量方程和能量方程，并进行改进。 连续方程： 其中和分别为液相和固相的体积分数。 动量方程： 和是液相和固相间动量交换量，是动量交换系数，是颗粒的虚拟质量力。 对于液体，动力粘性和导热系数随温度的变化都较明显，这种物性变化并定会对液体的流动和传热产生影响，本文中考虑了液体粘性和导热系数随温度的变化。 能量方程： 其中的热导率采用文献[4]中的近似模型。 对于固相的焓值： 其中L是相变颗粒的潜热总量，和是颗粒的相变起始和终止温度。 3. 数值计算 在本文的数值模拟中，离散方法采用有限体积法，控制方程中的对流项采用一阶迎风格式，压力和速度耦合采用SIMPLE算法。对于管2，径向网格数为30，轴向网格数为1000，共计。 对于固相和液相，中轴线都被定义为轴对称边界条件，计算在一个二维轴对称柱坐标中进行。针对动量方程，在入口处，固相和液相以相同的均匀速度（）计入流场。出口选用出流边界条件。对于固相和液相，管壁处均定义无滑移边界条件。针对能量方程，在管壁处，对固相和液相均定义常热流边界条件，在入口处固相和液相以相同的相变起始温度（）进入流场。在出口处，固相和液相均定义出流边界条件。 4. 计算结果 本文中在评估功能流体的冷却效果时，均是与同一流速下的纯水进行比较。而两种管道之间功能流体冷却性能的比较则在同一Re下进行。此外，在工程应用中，利用不同管径的管道对同一设备进行冷却时，很明显，加在管壁上的热流密度都是一样的，因此在本文中比较两种管道冷却效果时，均指在同一热流密度下。 4.1 固相浓度分布 由图1可以看出，MEPCM在径向的大部分区域分布均匀，但在壁面附近存在一个低浓度区，而且D = 2mm的管中，壁面附近MEPCM 浓度降低的更加明显。类似的现象已被其他研究者发现[4]，这说明颗粒存在径向迁移现象。颗粒的径向迁移的机理是较为复杂的，受剪切力，颗粒间的相互碰撞等多个因素的影响。单项流模型模拟功能流体的流动无法发现这一现象。 图1 速度完全发展段径向MEPCM浓度分布 颗粒在管壁附近存在低浓度区对功能流体的传热有两方面的影响，首先浓度的降低减少了潜热量，不利于冷却；另一方面，由于MEPCM的导热率明显低于水，固相浓度的降低意味着悬浮液导热率的增加，有利于传热。这两种影响的综合效果还需进一步研究。 4.2 壁面温度分布 为了比较的方便，定义有效雷诺数，其中为悬浮液的混合密度，为入口处速度和悬浮液粘性。定义温差，其中为功能流体对应的壁面温度，为纯水对应的壁面温度。 图2 时，两种管径下浓度10%的功能流体的冷却效果 如图2所示，两种管径下功能流体的冷却性能存