非固结多孔介质干燥的双尺度孔道网络模型与模拟：I．模型建立.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：083010 非固结多孔介质干燥的双尺度孔道网络 模型与模拟：I . 模型建立 1, 2 1, 3 2 1* 袁越锦 ，杨彬彬 ，徐英英 ，刘相东 （1.中国农业大学，北京 100083；2.陕西科技大学机电工程学院，西安 710021； 3.国家质检总局进 出口食品安全局，北京 100088） (Tel:029 Email: yyjyuan1@163.com) 摘要：以农产品颗粒物料中常见孔隙直径范围（10－4 ～10－3m ）的非固结多孔介质为研究对象，运用孔 道网络方法、多孔介质渗流物理和传递过程原理等知识，建立了考虑骨架吸湿、汽相对流、温度梯度 和孔道结构特征等因素对干燥过程影响的非固结多孔介质双尺度孔道网络干燥模型。与传统干燥模型 不同的是：该模型考虑了多孔介质中不同尺度下孔隙内热质传递机理的差别，且可直接输入多孔介质 物理结构特征参数，如孔隙率、粒径分布等，作为模型参量预测相应物料的干燥特性。 关键词：非固结多孔介质；干燥；孔道网络模型；双尺度；骨架吸湿性 0 引言 非固结多孔介质干燥涵盖的范围非常广泛，涉及到化工、食品、农业、轻工、建筑、 煤炭等众多领域。迄今为止，由于多孔介质本身的不均匀性和随机性，内部几何拓扑结 构的复杂性等原因，传统的干燥理论[1-4]大多都将这类物料直接视为虚拟的连续湿固体， 不考虑物料堆中颗粒体与孔隙的差别。而事实上，非固结多孔介质随机堆积，结构复杂， 在大部分情况下更接近于离散介质，其内包含多种不同级别的孔隙体系，如：在整体上 -3m 级别左右；另外，骨架颗粒自身 由颗粒堆积形成的孔隙体系，其孔隙半径一般在 10 -5m 以下。显然，干燥时在这些不同级别的孔隙 带有的孔隙体系，其孔隙半径一般在 10 体系中进行的热质传递将有很大差别，不能混为一谈。 孔道网络方法[5]在 20 世纪 90 年代初由 Daian[6]等人引入多孔介质干燥领域，经 Prat 和 Yortsos 等人的努力，逐渐形成了应用于干燥过程的孔道网络方法[6-14] 。这一方法认为 多孔介质的湿分迁移在孔隙中具有很大优势，采用节点（孔）和喉道描述多孔介质内部 的骨架拓扑结构，在孔道等级上研究多孔介质内的热质传递过程。该方法能很好的模拟 多孔介质干燥过程中出现的干、湿斑等现象；但在该方法的应用过程中也还存在着一些 亟待解决的问题，如：如何考虑多孔介质中不同尺度下的孔隙差别？不同尺度下的热质 传递机理需做哪些调整？等等；此外，该方法在处理吸湿性骨架问题时通常将物料简化 资助项目：国家自然科学基金（No）；陕西科技大学科研启动基金（No.BJ07-06 ） 为刚性多孔介质[8, 10 ]，而实际上自然界中大部分多孔介质骨架是吸湿的，如粮堆、土壤 等。从本质上相对多孔介质整体而言，骨架吸湿性问题实际上是一个在更小尺度上的质 量传递过程。 本文以农产品颗粒物料中常见孔隙直径范围（10－4 ～10－3m ）的非固结多孔介质为 研究对象，运用孔道网络方法、多孔介质渗流物理和传递过程原理等知识，建立考虑骨 架吸湿、汽相对流、温度梯度和孔道结构特征等因素的自然非固结多孔介质双尺度孔道 网络干燥模型，以期为后续模拟分析非固结多孔介质干燥特性的研究工作打下理论基 础。 1 二维双尺度孔道网络物理模型 1.1 物理模型构建 如前所述，非固结多孔介质内部孔隙拓朴结构非常复杂，为便于数学处理，假定其 骨架为球体，如图 1(a)所示。一般地，描述非固结多孔介质结构特征的参数有：孔隙率、 孔隙个数面密度、粒径分布、孔隙直径分布等。孔道网络模型的构建过程一般包括网格 提取、规则化网格和建立物理模型等步骤( 图 1 （b ）、（c ）、（d