凹槽通道中流动和换热问题的分区迭代并行计算.pdf

中国工程热物理学会 传热传质学 学术会议论文 编号：093470 凹槽通道中流动和换热问题的 分区迭代并行计算 李峥，杨茉 （上海理工大学，上海 200093） (Tel: 021 Email:talent2045@163.com) 摘要：对凹槽通道中流动与换热问题进行了数值模拟，选取与文献中实验结果相同的几何参数与边界 条件。分析随着雷诺数增加凹槽通道中的流动与换热情况，并与实验解进行比较分析得出在低雷诺数 情况下，凹槽区域出现流动分区情况，即主流区、自由剪切层、涡流区，高雷诺数情况下出现文献中 介绍的自维持震荡情况。为解决凹槽问题计算时间长的问题，本文采用基于 MPI 技术的并行程序设计 方式，并对并行效率进行一定程度的分析。 关键词：数值模拟；流动与换热；凹槽通道；并行计算 0 前言 平行板通道及其下壁面具有凹槽的平行板通道在工程中有广泛的应用，一些高效换 热器、生物医学设备、电子器件等中的流动和换热问题都能抽象为这种研究模型。而通 道一般较为狭长，进行数值模拟时长度方向需要分配较多网格[1]，计算量非常大，所需 运算时间长。 与此同时，计算机系统正面临深刻的变革，传统的 von Neuman 格局已经被突破， 采用并行化结构的并行计算机[2]正日益普及，并且在科学与工程计算中发挥越来越重要 的作用。计算机系统结构的并行化蕴含着提高运算速度和增加信息存储量的巨大潜力。 计算机发展的并行化趋势，必然会促使算法设计的并行化。而FORTRAN+MPI 正是可将原 有串行程序扩展为并行程序的较成熟的语言，适应了计算机的发展趋势，可以达到“空 间换时间”的目的。 本文通过 MPICH 在一小型局域网搭建了并行计算平台，并基于区域分解原理进行 了交错网格下的并行 SIMPLE[3] [4]算法的研究，给出了该算法实施的具体步骤以及实现 数据交换的措施。最后，选取与Herman 和Kang[5]采用激光全息干涉方法对凹槽通道中 的流动和换热进行了实验研究，相同的几何条件与边界条件，讨论随着雷诺数的增加通 道内流动与换热情况，并与Herman 和Kang 所得的实验解进行比较，验证了该平行算法 的可靠性，考察了算法的并行性能。 1 物理模型与数学描述 1.1 物理模型 凹槽通道内流体对流换热问题的物理模型如图3-1 所示，其中凹槽通道是由10 个矩 形突出块所构成的。假设流动和换热满足以下条件 （1） 问题为二维不可压流体； （2 ） 通道上壁面绝热，下壁面分布10个矩形凸出块，构成凹槽通道，每个凸出块 底面为横热流q加热，凹槽底部绝热； （3 ） 出口取局部单向化条件； (4) 常物性 图1 凹槽通道内的流动和换热模型示意图 Fig.1 Flow and heat transfer model in grooved channel 1.2 数学描述 动量方程 2 2 U ∂U ∂U ∂P 1 ∂ U ∂ U +V − + ( 2 + 2 ) X Y X Re ∂ ∂ ∂ X Y ∂ ∂ 1 2 2 ∂V ∂V ∂P ∂ V ∂ V U +V − + ( 2 + 2 ) ∂X ∂Y ∂Y Re ∂X ∂Y 连续性方程（质量守恒） ∂U ∂V + 0 ∂X ∂Y 能量方程 λ