矩形管道湍流流动与换热的直接数值模拟.pdfVIP

工程计算 363 矩形管道湍流流动与换热的直接数值模拟 马良栋 李光熙 陶文锉 李增耀 (西安交通大学能源与动力工程学院动力工程多相流国家重点实验室，陕西 西安 710049) 【摘要】本文对矩形管道内充分发展湍流流动和换热进行了数值模拟，湍流雷诺数Re(Rem=6200)为 400，普朗特数为0.71。网格数为128x128x128o湍流统计结果与已有试验及计算数据进行了比较，得 到了比较好的一致性。分析了近壁面湍流特性及雷诺应力的变化，以及分析平均温度沿壁面平分线的变化 关系和局部努谢尔特数Nu沿壁面的变化。为进一步分析和评价湍流模式理论提供了依据。 关键词:矩形管道，湍流，直接数值模拟，雷诺应力 DirectNumericalSimulationofTurbulentFlowandHeatTransfer inaSquareDuct MaLiangdong,LiGuangxi,TaoWenquan,LiZengyao (StateKeyLaboratoryofMultiphaseFlowinPowerEngineering,XianJiaotongUniversiyt,Xian710049,China) Abstract:Inthispaper,adirectnumericalsimulationofafulldevelopedturbulentflowandheattransferisstudiedin asquareductwithanimposedtemperaturedifferencebetweenthehorizontalwallsandtheperfectlyinsulatedonthe verticalwall.Thenumericalschemeemploysatime-splittingmethodtointegratethethreedimensionalincompressible Navier-Stokesequation.TheunsteadyflowfieldwassimulatedataReynoldsnumberof400(Rem6200)basedonthe meanrficationvelociytandtheductwithusing128x128x128gridpointsandPr--0.71.Theturbulentcharacteristics andtheReynoldsstressclosingthewallisanalyzed.ThemeantemperaturealongthewallbisectorandthelocalNu alongthewallisstudied,andaccordingtothesedatedifferentturbulentmodelscanbeassessedandimproved. Keywords:squareduct;turbulent;directnumericalsimulation;Reynoldsstress 1引言 [8-9]分别对Re:为300和600的矩形管道充分 发展湍流进行了直接模拟。但这些工作只对湍流 在实际工程中，如热交换器，通风及空调系 流动进行了研究，关于传热问题研究甚少，仅文 统管道，核反应堆等，矩形及非圆形管道广为采 献 【10]对湍流换热进行了研究。本文采用直接模 用。在矩形管道中由于非圆管道截面会引起二次 拟的方法对湍流雷诺数为400的充分发展流动和 流动的产生。普朗特第二类二次流是微弱的，仅 换热进行了数值计算，计算中采用基于MPI的并 为主流方向速度的2一3%11，但对主流速度有很 行环境，得到了与文献一致的结果，对切向和法 大的影响。很多研究者对矩形管道流动进行了实 向雷诺应力以及对其所产生的压力进行了详细分 验和数值研究2[[,31，文献 【不〕对其进行了详细的 析，揭示了二