干燥发展的新态势及一种新方法.docVIP

干燥发展的新态势及一种新方法可以 应用在干燥微观领域研究 战洪仁 张建伟 李冬艳 （沈阳化工学院机械工程学院，辽宁 沈阳，110142） 摘 要：分析了当前干燥学科发展的新态势及干燥微观领域研究状况，并简单介绍了格子气自动机法的发展状况及FHP模型和扩展模型LG模型原理，并用LG模型对液体的相变过程进行模拟，从模拟结果可以清晰地看到气泡的产生、生长、分解过程。在模拟相变过程的基础上，分析了LG模型应用在干燥微观领域研究的可行性。 关键词：干燥；新态势；格子气自动机； LG（liquid-gas）模型；相变；微观；数值模拟 1前言 干燥是一门庞大的学科，由于干燥涉及的面广，在国民经济中占有尤为重要的地位，所以近年来在全球范围倍受人们关注格子气模拟流场，是直接从Boltzmann碰撞模型出发作完全离散处理，FHP模型是将流场划分成间距为1的正三角形网格，时间也离散成整时间步（t=0,1,…, …） 图2 FHP模型的基本碰撞规则 依据节点上有无静止粒子，FHP模型有FHPⅠ模型（节点上无静止粒子）；FHPⅡ模型（节点上至多一个静止粒子）；FHP Ⅲ模型（节点上可允许有多个静止粒子）。 FHPⅡ模型网格节点的状态可按粒子运动方向采用二进制数的方式来表示。假设用ｓ(ｘ,ｔ)表示节点ｘ在时刻ｔ的状态，即：s(x,t)=( nα(x,t); α=0,1,2,3,4,5,6) 其中nα(x,t)为boole变量，表示x节点t时刻以速度eα运动的粒子数占有的情况，对应于nα=1或0，eα={cos(πα/3),sin(πα/3)}, α=0,…,5.静止状态的粒子速度为e6={0，0}。即FHPⅡ是具有7个速度的二维空间模型，用7位二进制数表示节点的状态如图3所示的格子中在α=0和α=2的方向上有速度，其状态可表示为：s(x,t)={1，0，1，0，0，0，0}。 图3 FHPⅡ格子气模型的节点状态编码与二进制表示 5 LG（liquid-gas）模型的原理 干燥过程不但有液态湿分的移动，还会伴随湿分的气化过程， FHP模型是以流体的一种粒子（液体粒子）为研究对象，为了模拟湿分的气化现象，在FHP模型基础上，建立了以气液两相粒子为研究对象的LG（liquid-gas）模型。LG模型相邻的粒子的碰撞规律与FHP模型完全一致，相距为r个单位长度的粒子之间有相互作用的引力和斥力，粒子在整个流体场遵循质量、动量、能量守恒定律，随着时间的发展反复地进行着碰撞、移动、相互作用，以实现相变过程[15.16]。相互作用规则如图4所示（实心箭头表示相互作用前的方向，实心点表示静止粒子，空心箭头表示相互作用后的方向）。 6 LG模型模拟的相变过程 如图5 所示，模拟区域是在一个密闭容器里，容器里充满液体，整个容器被304×512个节点所划分，图中是不同液体密度，不同粒子相互作用半径的LG（liquid-gas）模型模拟出来的液体相变过程。从图中可以看出，模拟随着时间步的进行，液体中逐渐有气泡产生，随着时间步的增加，气泡越来越多，并且小气泡逐渐结合成为大气泡。 模拟过程中发现，改变相互作用半径r，气化的难易程度也要发生变化，在液体密度不变的情况下，r越大，产生气泡的速率慢，气化越难发生，反之，r减小，产生气泡的速率快，气化容易发生。说明相互作用半径r是温度的倒函数，这与范德瓦尔方程是符合的，从而也证明了此模型能够模拟气液相变过程。 t=54 t=220 (a)无量纲液体密度为ρ=0.2粒子/节点，粒子相互作用半径为6 t=27 t=180 (b)无量纲液体密度为ρ=0.3粒子/节点，粒子相互作用半径为4 图5 不同密度、不同相互作用半径的模拟结果 7 LG模型模拟微观干燥过程的可行性分析 在FHP模型基础上，增加相互作用粒子的LG模型，可以模拟液体的相变过程，LG模型是以流体的分子运动论描述为基础的，从Boltzmann碰撞模型出发作完全离散处理，模型符合基本的物理定律，即质量、动量、能量守恒定律。由统计物理学知道，复杂系统的宏观特征是由微观粒子相互作用，而在整体上表现出来的统计规律，所以模型的宏观平均行为符合宏观的微分方程。实际上这是介于微观和宏观之间的中介模型，因为严格的说它并不是分子层次上的模型，但由于此类模型源于分子运动论描述，分子动力学的一些优点便得以保持，所以他能用来计算和模拟复杂系统中物质的分布和变化[17]。 微观的干燥过程，是典型的化学工