气流式喷雾干燥过程的计算流体力学模拟.docxVIP

doi：10.13732/j.issn.1008-5548.2014.05.004气流式喷雾干燥过程的计算流体力学模拟马川川， 张丽丽（齐鲁工业大学 机械与汽车工程学院 ， 山 东 济 南 250 doi：10.13732/j.issn.1008-5548.2014.05.004 气流式喷雾干燥过程的计算流体力学模拟 马川川， 张丽丽 （齐鲁工业大学 机械与汽车工程学院 ， 山 东 济 南 250300） 流 态 化 研 究 摘 要 ： 基于计算流体力学技术和气-固 两 相 流 理 论 ，对气流式喷雾干燥 过程进行非稳态数值模拟 ，得到干燥塔内热空气的温度 、速 度 场 分 布 云 图 ，以及颗粒的运动轨迹 、停留时间和粒径分布 。 结 果 表 明 ，改 变 热 风 温 度 、适当延长干燥时间都可使料液干燥完全 ，有助于提高粉体 的 干 燥 效 率 ，改 善 产 品 质 量 。 关 键 词 ：喷 雾 干 燥 ；非 稳 态 ；计 算 流 体 力 学 ；数 值 模 拟 渗 透 ，计 算 流 体 力 学 （CFD）为解决该问题提供了一种 新思路 。 本文中利用计算流体力学技术对气流式喷雾 干燥过程进行数值模拟 ， 得到干燥塔内 热 空 气 的 温 度 、速度场分布云图 ，以 及 颗粒的运动轨迹 、停 留 时 间 和 粒 径 分 布 ，为进一步研究喷雾干燥的干燥 、雾 化 机 理，制备高质量的干燥产品提供参考 。 中 图 分 类 号 ：TQ021 文 献 标 志 码 ：A 文 章 编 号 ：1008-5548（2014）05-0015-04 1 喷雾干燥的计算流体力学模型 1.1 连续相控制方程 热空气的连续相 湍流控制方程包括质量 、 动 量 、 能量、本构方程及水气输运方程 ，其中输运方程为[1] Computational Fluid Dynamics Simulation on Spray Drying Process of Pneumatic Atomization （ρφ） +div(ρuφ)=div（Γφgradφ）+Sφ， （1） t MA Chuanchuan， ZHANG Lili （School of Mechanical and Automotive Engineering， Qilu University of Technology， Jinan 250300， China） 式 中 ： ρ 为 流 体 密 度 ，kg/m3； φ 为通用求解变量 ； t 为 时 间 ，s； u 为流体的速度矢量 ，m/s； Γφ 为 广 义 扩 散 系 数 ； Sφ 为广义源项 。 离散相控制方程 对 于 气 体 、颗粒两相流动系统 ，考 虑 作 用 在 颗 粒 上 的 各 种 力 ，根据牛顿第二定律 ，可得拉格朗日坐标 下喷雾干燥过程中颗粒的轴向速度方程 ： 1.2 Abstract: The spray drying process of pneumatic atomization was simulated using the method of computational fluid dynamics and gas-solid two-phase theory. Contours of hot air velocity and temperature as well as tracks， residence time and size distribution of particles in the spray dryer were presented and analyzed. The results show that the products with high quality and efficiently can be produced if the temperature of hot air and drying time are changed and prolonged properly. Key words: spray drying； unsteady； computational fluid dynamics； numerical simulation dvx，p =Γp（vx-vx，p）+gx（ρp-ρ）+Fx， （2） dt 式 中 ：vx， p 为 颗 粒 速 度 ，m/s； vx 为 流 体 相 速 度 ，m/s； ρp 为 颗 粒 密 度 ，kg/m3； gx 为重力加速度轴向分量 ，m/s2； Fx 为雾滴受到的其他作用力 ，N； 颗 粒 驰 豫 时 间 Γp（单 位为 s）为 气