微通道中水蒸汽凝结过程流型变化分子动力学研究.pdfVIP

中国工程热物理学会 工程热力学与能源利用 学术会议论文 编号：081083 微通道中水蒸汽凝结过程流型变化的 分子动力学研究 刘彬武 刘朝 刘娟芳 （重庆大学动力工程学院 重 庆400044） Tel: 023 E-mail:liuchao@cta.cq.cn 摘要 本文首先建立了在外力作用下流体流动和凝结的分子动力学模型，模型能够很好地模拟平行微通 道间流体流动和凝结问题。以模型为基础对不同条件下水蒸汽凝结过程的流型变化进行了研究，模拟 结果和文献中实验流型吻合较好。同时本文对流型间转换机理进行了探索，本研究是一种两相流流型 的全新研究方法，对以后开展这方面的研究具有一定的指导意义。 关键词： 分子动力学 微通道 流型 1 引言 微通道内汽液两相流的流型与常规尺度下的流型有很大的区别，主要表现为微通道 内汽液两相流中分层流不会发生，而出现弹状–环状流。从目前的研究来看，不同的研 究者观测到的流型不太一致，对机理的研究很少[1、2]，而且多是通过试验的方法进行在 线观测。本文试图通过分子动力学计算的方法开展微通道内汽液两相流流型的研究，并 试图对其中流型转换的机理作一些解释。 2 分子动力学研究方法 本文采用分子动力学方法来模拟微通道中流型的变化，这是一种全新的研究方法， 有着其他方法不可比拟的优越性，但同时也受到计算机计算能力的限制。所以不可能像 试验装置那样将通道设计得太长，因为通道越长，要模拟的分子数就越多，这是 PC 机 所不能承受的。我们采用变通的处理方法，不像试验那样将通道在流动方向上分成若干 个温度不同的区域（因为分区太多，会导致通道太长），而是采用恒定壁温。实验研究 方法能在通道不同的位置观测到不同的流型，而本文采用了恒壁温的模拟方法，能够在 凝结过程中不同的时刻观测到不同的流型，达到与试验同等的效果。 补充区域 Z X Y 图 1 分子动力学模拟模型 资助项目：国家自然科学基金资助项目（NoNo 图 1 是采用分子动力学方法对凝结过程中流型的变化进行模拟的示意图。模拟盒在 X 方向的长度为 Lx=29.023σ(8.8nm),Y 方向长度 Ly=Lx/2,Lz 分别为 Lx/8 ，Lx/6 ，Lx/4 ， Lx/2 ，2Lx/3 ，流体粒子初始个数为 864 个，壁面粒子为 2048 个。过热蒸汽在外力驱动 下在通道中从左至右沿 x 方向流动，接着对壁面降温，（壁面温度低于汽相压力下的饱 和温度）。处理方法是根据给定温度调节壁面上反射粒子的速度，相当于在壁面上固定 一个恒温器，由于温差作用，蒸汽分子在固壁上凝结。模拟中若始终保持流体分子数恒 定，随着凝结不断进行，蒸汽分子将不断减少，冷凝过程必将很快趋于平衡，流型将不 会再发生变化。为了保证凝结过程持续稳定地进行，需要对流体分子进行补充。在通道 入口处设置一定长度的补充区域，并保证该区域内流体始终处于 300℃过热态，凝结过 程中每 1000 步计算一次汽相空间的密度，如果密度小于汽相初始密度，说明有部分粒 子凝结在了壁面上，此时就需要对汽相空间进行粒子补充，补充粒子随机从补充区域读 入，并随机安放在通道入口处。值得说明的是，在凝结过程中，由于汽相空间在不断减 小，需要根据凝结状况对汽相空间大小作实时变化。 通过上面一系列的处理，使得在不同通道宽度的情况下汽相密度基本恒定。本文首 先对一个较高压力下的凝结过程进行了研究。 在研究通道宽度对流型的影响时，统一设定外力作用力为无因次 2.0