ADS钨球流态固体散裂靶换热器的数值模拟研究..docx

ADS钨球流态固体散裂靶换热器的数值模拟研究刘捷*，刘欢，卢文强中国科学院大学工程科学学院，北京，10049*Tel：010 Email:nauty@ucas.ac.cn摘要：本文在优化后的换热器模型的基础上数值模拟方法计算了钨球多颗粒在竖直放置的平板换热器中与氦气的换热过程，追踪捕获了颗粒在换热器中的运动轨迹，发现颗粒与换热器的碰撞次数对其下落过程影响很大，碰撞次数越多，颗粒在换热器中运动时间越长，出口温度越低。并对比了不同速度氦气与钨球颗粒的换热效果，最终确定了换热器最低氦气速度，上述换热器数值模拟结果为ADS流态固体散裂靶换热器的可行性及实验提供了可靠的数据参考，为未来ADS固体散裂靶的发展奠定了基础。关键词：颗粒换热器；CFD-DPM；DEM；多项流；数值模拟引言加速器驱动次临界堆系统(Accelerator Driven System, ADS)是未来安全处理核废料的首选技术途径[1-3]。ADS系统由加速器、散裂靶和反应堆三大子系统构成，其中散裂靶连接加速器和反应堆，是整个系统的“心脏”[4]。由于欧洲铅铋合金气液散裂靶存在研制和运行费用高、温度-材料腐蚀效应严重并伴有次生放射性产物毒性高等缺点，中国科学院詹文龙副院长[1,3]于2013年提出了钨球流态固体散裂靶的新设计方案（如图1所示），并进行重点探索和研究。钨球流态固体散裂靶的主要功能是作为靶材料被质子束轰击产生中子，中子再轰击到布置在靶体周围的核废料，从而引起核废料的嬗变，降低其放射性。然而钨球流态固体散裂靶中的钨球颗粒在被质子束轰击之后，发生散裂反应，会产生大量的热沉，使散裂靶出口的钨球温度过高，远高于靶进口要求的温度，无法直接实现钨球的循环利用。因此需要在钨球出口和入口之间设计一个换热器，使钨球的温度降低以达到入口所要求的温度。图1 钨球流态固体散裂靶示意图[1]颗粒换热器中颗粒与流体的换热方式主要包括颗粒与流体之间的对流换热、颗粒与颗粒之间的换热及颗粒与换热器壁面之间的换热，这三种换热过程涉及导热、对流和辐射三种传热机理，并与颗粒的运动紧密关联。颗粒换热与颗粒运动的这种极强耦合性导致对其换热机理的研究非常困难，但由于其广泛存在于化工、环境、能源和生命科学等领域，并与流化床、固定床、燃烧等应用密切相关，因此一直吸引着国内外学者对其进行大量研究。在颗粒运动方面，刘马林[5]、柳朝晖[6]和郑梦娇[7]等使用格子Boltzmann方法对单颗粒、双颗粒沉降进行数值模拟研究，结果与实验符合很好。曲长福[8]等获得了多颗粒均匀分布系统中单颗粒所受到的气动力与空隙率的关系。Chi Zhu[9]等实验研究了小直径球形单颗粒在以不同速度流动的流体中的非定常运动，发现Basset force在颗粒的非定常运动中贡献很小，可忽略不计；并获得了与实验结果符合很好的阻力系数CD-Re公式。S.M. Dash[10]等基于浸入边界和格子Boltzmann法分别数值模拟气体与静止颗粒和运动颗粒的相互作用，发现随着Re的增加，颗粒后面回流区出现非对称现象，流动呈现不稳定性。在颗粒换热方面，朱立平[11]等基于离散单元法(DEM)提出了综合考虑颗粒碰撞传热、颗粒内部导热以及颗粒与气体间的对流换热的模型，并数值模拟固定床中丝状颗粒的传热过程，发现对流换热对整体传热量的贡献较大。Niels G. Deen[12]等给出了计算稠密颗粒-气体两相系统的平均换热系数及当地换热系数的公式，并用数值模拟验证了其公式的准确性。J.Sun和M.M.Chen[13]引入弹性碰撞理论，从理论推导出两颗粒在碰撞过程中的传导热量。国际著名颗粒学者、澳大利亚新南威尔士大学余艾冰教授[14-15] 考虑了非接触颗粒间的气体导热，对Sun和Chen颗粒碰撞的导热公式进行了修正，并用有限单元法进行了验证。C.S. Bu[16]建立了包含颗粒内部的传热、粗糙表面的传热、颗粒表面气膜的传热和紧密接触粗糙表面的孔隙传热四种导热机理的理论模型并用并用离散单元法进行计算，指出当颗粒温度在700K以下时，颗粒的辐射换热可以忽略。Z.Y. Zhou[17-19]等在颗粒尺度计算中采用含导热、对流和辐射三种方式的模型计算了流化床中的换热系数，并给出了上述三种换热方式计算的经验公式，与实验结果吻合较好。颗粒碰撞之间的瞬时性及其表面变形给颗粒碰撞换热研究带来了困难，已发表的成果中并无颗粒碰撞换热统一的模型及公式，这导致准确描述颗粒与气体之间的换热，尤其是稠密流动颗粒系统仍充满了挑战。本文在理论计算确定的换热器模型基础上，对气-固两相换热过程进行了数值模拟研究，通过追踪颗粒的运动，研究多颗粒系统的气固换热过程，并对比了气体速度对两相系统换热的影响，为未来ADS流态固体散裂靶颗粒换热器的实验及发展提供了参考数据。1 数理模