降膜结晶分离过程的数值分析.pdfVIP

中国工程热物理学会 传热传质学学术会议 降膜结晶分离过程的数值分析。 张建文 张政 (北京化工大学化学工程学院，北京，100029) 摘要本文在混合物理论和实验观察所获得的结晶形态及形貌特，怔的基础E建立了 描述降膜结晶过程液相区、两相区及固相区的动量、热量和质量传递的分区统一数学 模型。本模型的一个重要特点是认为在降膜结晶过程中于熔体液相区和固相区之问存 在一薄厚不匀且不断变化的二相区。以二元有机熔体(萘一硫茚)为删，针对实验管 式降膜结晶器的实际降膜结晶过程动量、热量和质量传递过程进行了数值模拟．所得 结果在未引入任何经验常数的条件下．与实验结果符合较好。 关键词：降膜结晶 传热传质 液固二相 多孔介质 数值模拟 1引言 降膜结晶分离过程(FFC process)独具特点，在工业上得到极为J“泛的应用。它 以高的分离效率和高生产能力而受到人们的日益重视。国内外业已投入大量人力、物 力进行降膜结晶过程的研究开发【“】。本文作者在实验基础上基于双流体理论和混合物 理论建立了降膜结晶过程数值模拟的液固二相分区统一模型并应用该模型对降膜结晶 过程的复杂相变耦合传热传质问题进行数值模拟，模拟结果由我们的实验结果进行验 证。 2降膜结晶过程基本假设及模型方程 2．1降膜结晶过程描述 在实验观察的基础上，降膜结晶过程可被看成是一个双组分，固液二相的连续统 内进行的流动、传热传质同时存在并具有熔体固化结晶相变的复杂物理过程。如图I所 向液相区扩大的过程，可阻认为结晶方向沿r方向。而在两相区内．晶体自固体区侧 向液体区横向(r方向)不断扩大直至最终完全固化成固相区。鉴于啦个枝晶横向(x方 向)尺寸仅有10。3～10。mpJ，与结晶区x方向以及=降制客积尺寸相比很小，在实际计算模 型中，可以不考虑枝晶的细微形貌结构上的不同，而将之视为多孔介质。 2．2基本假设： 除上述描述外，为建立基本数学模型，进一步作如下假设： (1)混台物性质可以看成是各组分性质简单按所占份额的数学加和： (2)固相不运动，即Vs=O，忽略固体内的扩敞传质(Schiel假定)[7-10]； (3)下降液膜在流动过程中仅受到重力的驱动，压力分布均匀，等于环境外压： (4)固、液相有关物性参数均为常数，不随结晶过程温度、浓度的变化而变化： 固、液相密度相等，P。=pI=p： ‘国家自然科学基金资助项目，No59876(104 Ⅲ一25 (5)忽略结晶过程初始过冷度的影响，认为结晶是平衡结晶状态。考虑到枝晶及间 隙尺if-(IO’5～10一}m)远小于可能取的微观控制体，可忽略其实际几何形态微结构的影 响，根据显微镜热台熔体定向结晶结果可近似将结晶前锋看成是锥柱：}l}，如图l所示； (6)将阿相区域晶体之间多孔介质 内渗流运动的流体视为牛顿流体； (7)任一点液圜相达到局部热平 衡，具有相同的温度五 (8)微观晶体生长为局部平衡生长 过程，在不同生长时刻，锥柱状晶体 沿生长方向向前推移，如图l阴影线所 示区域。各新生层次的局部结晶固相 浓度为相应时刻当地温度下的平衡固 相浓度，即CAs=(■。(T)。 。，麓簿鬻o_ 2．3降膜结晶数值模拟模型方程 参考图l，根据双流体理论和混合 物连续介质理论m。1…．在上述假设基础 上可得到固相不运动下降膜结晶二元 同液相结晶过程的守恒方程如下： 2．3．1质量守恒方程 1 Geometricalstructure region Figure ofmultiphage (1) illustrationof V·(g，V)=0 and c