甘露醇喷雾干燥过程中液滴粒度分布变化的群体粒数衡算模拟和实验研究

doi:10.16085/j.issn.1000-6613.2018-0752

化工进展 ›› 2019, Vol. 38 ›› Issue (02): 772-778.DOI: 10.16085/j.issn.1000-6613.2018-0752

甘露醇喷雾干燥过程中液滴粒度分布变化的群体粒数衡算模拟和实验研究

吕凤^1,²(),张扬¹,马才云²,王学重^1,²()

1. 华南理工大学化学与化工学院，广东广州 510046
2. 利兹大学化学和过程工程学院，英国利兹 LS2 9JT

收稿日期:2018-04-12 修回日期:2018-06-05 出版日期:2019-02-05 发布日期:2019-02-05
通讯作者: 王学重
作者简介:<named-content content-type="corresp-name">吕凤</named-content>（1986—），女，博士研究生，研究方向为喷雾干燥过程模拟。E-mail：<email>celvfeng@mail.scut.edu.cn</email>。|王学重，教授，博士生导师，研究方向为过程控制和三维成像。E-mial：<email>xzwang@leeds.ac.uk</email>。
基金资助:
国家自然科学基金(91434126, 61633006);广东省自然科学基金(2014A030313228, 2017A030310262);广东省科技应用项目(2015B020232007);中央高校基本科研业务费(2017MS092)

Simulation and experimental study on the evolution of droplet size distribution during spray drying of mannitol

Feng LÜ^1,²(),Yang ZHANG¹,Caiyun MA²,Xuezhong WANG^1,²()

1. School of Chemistry and Chemical Engineering, South China University of Technology, Guangzhou 510640, Guangdong, China
2. School of Chemical and Process Engineering, University of Leeds, Leeds LS2 9JT, U. K

Received:2018-04-12 Revised:2018-06-05 Online:2019-02-05 Published:2019-02-05
Contact: Xuezhong WANG

摘要/Abstract

摘要：

利用群体粒数衡算（population balance，PB）计算机模拟和实验研究了甘露醇水溶液的喷雾干燥过程中液滴的粒度分布的变化规律。液滴干燥过程中的颗粒粒度的萎缩速率，在群体粒数衡算模型中描述为液滴的逆（或负）生长项，通过单个液滴反应动力学方法（reaction engineering approach，REA）获得。基于单个液滴干燥的反应工程方法模型REA和群体粒数衡算模型PB集成建立了PBREA模型。PBREA 模型的求解是通过高分辨率数值方法。本文模拟研究了不同工况下，不同粒径液滴的干燥时间、液滴平均含湿量以及液滴粒度分布随时间的变化。结果显示，液滴粒径越大，干燥时间越长，模型预测的颗粒平均粒径为实验值的1.0~1.5倍，粒度分布跨度是实验值的0.61~0.89倍。模拟误差主要来源于液滴及颗粒粒径分布统计精度、单个静止液滴与群体运动液滴干燥的差异、热导率及扩散系数是经验值3个方面。在使用Buchi 290 小型喷雾干燥仪进行的实验中，使用了图像采集和分析方法得到了液滴及颗粒的数密度分布，并和模拟结果做了对比。结果表明该模型可以有效地预测喷雾干燥过程中干燥颗粒的平均粒度及分布跨度。

关键词: 干燥, 粒度分布, 数学模拟, 群体粒数衡算模型, 反应工程方法

Abstract:

Computer simulation and experimental study were carried out on the evolution of droplet size distribution during spray drying of mannitol dissolved in water. The shrink rate of the diameter of a droplet during spray drying was treated as a negative growth rate in the population balance (PB) model which was obtained using the reaction engineering approach (REA).The integration of PB and REA yielded the PBREA model, which was solved using a high-resolution numerical method. The drying time of droplets of varied sizes, the change of the mean moisture content of droplets and the droplet size distribution were simulated under different spray drying conditions. The results showed that the drying time increased with the increase of droplet diameter. The ratio of the predicted and measured particle mean diameters was between 1.0 to 1.5, and the span was from 0.61 to 0.89. The errors were analyzed and attributed to three main factors: error due to the statistical analysis of droplet and particle sizes, the difference in drying a single static droplet and a group of droplets in motion, and the empirical thermal conductivity and diffusion coefficient values. Images were collected and analyzed to obtain droplet and particle size distributions with the Buchi 290 spray dryer. The simulated and experimental results showed that the PBREA model could effectively predict the mean diameter and span of particles during spray drying process.

Key words: drying, particle size distribution, mathematical modelling, population balance, reaction engineering approach

中图分类号:

吕凤, 张扬, 马才云, 王学重. 甘露醇喷雾干燥过程中液滴粒度分布变化的群体粒数衡算模拟和实验研究[J]. 化工进展, 2019, 38(02): 772-778.

Feng LÜ, Yang ZHANG, Caiyun MA, Xuezhong WANG. Simulation and experimental study on the evolution of droplet size distribution during spray drying of mannitol[J]. Chemical Industry and Engineering Progress, 2019, 38(02): 772-778.

图/表 10

表1 喷雾干燥操作参数及颗粒产品物性参数

工况	入口温度/℃	母液中甘露醇的质量分数/%	雾化气体流量/L·min^-1	出口温度/℃	产率/%	D ₅₀/μm	分布跨度	产品含湿量/%
A	135	10	4.1	86	71.3	9.36	1.62	0.38
B	135	10	5.95	85	84.4	7.24	1.39	1.22
C	135	15	4.1	86	72.7	11.76	1.57	0.70
D	135	15	5.95	85	80.2	8.01	1.44	0.74

图1 显微镜下的二甲基硅油中的样品A液滴

图2 样品A颗粒的扫描电镜图像

图3 液滴/颗粒粒径的累积密度分布

图4 4种工况下不同粒径液滴/颗粒的平均含湿量随时间的变化

图5 不同浓度及粒径液滴的干燥时间

表2 模拟值与实验结果对比

样品编号	实验值/μm			实验值分布跨度	PBREA模拟值/μm			PBREA模拟值分布跨度	REA模拟值 D ₅₀/μm	PBREA模拟值/ 实验值		REA D ₅₀模拟值/实验值
样品编号	D ₁₀	D ₅₀	D ₉₀	实验值分布跨度	D ₁₀	D ₅₀	D ₉₀	PBREA模拟值分布跨度	REA模拟值 D ₅₀/μm	D ₅₀	分布跨度	REA D ₅₀模拟值/实验值
A	3	6	12	1.5	5.00	9.00	17.00	1.33	11.58	1.50	0.89	1.93
B	2	5	12	1.3	5.00	7.00	11.00	0.86	8.86	1.40	0.66	1.77
C	3	9	17	1.6	6.00	10.50	18.00	1.14	13.40	1.17	0.71	1.59
D	3	7	12	2	3.50	7.00	12.00	1.21	8.30	1.00	0.61	1.19

表3 4组液滴干燥过程中的总萎缩率

初始液滴平均直径/μm	颗粒与液滴平均粒径的比值/%
初始液滴平均直径/μm	实验值	PBREA模拟值	REA模拟值
17	35.29	52.94	68.12
13	38.46	53.85	68.12
21	42.86	50.00	63.83
14	50.00	50.00	63.83

图6 4个样品的粒度分布随时间的变化

符号说明

A	——		液滴的表面积， m²
a ₁, a ₂, a ₃	——		拟合系数
b	——		REA模型中液滴粒径变化的系数
$c p$	——		液滴的定压比热容， kJ/(kg·K)
D	——		液滴或颗粒直径， m
D ₁₀ ，D ₅₀ ，D ₉₀	——		累积粒度分布分别为10%， 50%和90%时对应的颗粒直径
$D f$	——		扩散系数， $m 2 / s$
G	——		液滴或颗粒的萎缩速度，m/s
h	——		传热系数， $W / (m 2 ? K)$
$h m$	——		传质系数， m/s
k	——		时间步长
$k f$	——		热导率， $W / (m ? K)$
L	——		特征尺寸
m	——		质量， kg
Pr	——		普朗特数
Q	——	种群的密度函数
R	——	气体常数
Re	——	雷诺数
T	——	温度， K
t	——	时间， s
X	——	液滴或颗粒的含湿量， kg水/kg固体
z	——	尺寸步长
ΔH	——	水的气化潜热， $k J / k g$
ρ	——	密度， $k g / m 3$
$? E$	——	干燥过程中的活化能， J/mol
上、下角标
b	——	干燥气体
c	——	平衡状态下的参数
d	——	液滴
s	——	表面
sat	——	饱和状态下的参数
t	——	时间
v	——	蒸汽
0	——	初始状态下的数值
L	——	特征尺寸

符号说明

A	——		液滴的表面积， m²
a ₁, a ₂, a ₃	——		拟合系数
b	——		REA模型中液滴粒径变化的系数
$c p$	——		液滴的定压比热容， kJ/(kg·K)
D	——		液滴或颗粒直径， m
D ₁₀ ，D ₅₀ ，D ₉₀	——		累积粒度分布分别为10%， 50%和90%时对应的颗粒直径
$D f$	——		扩散系数， $m 2 / s$
G	——		液滴或颗粒的萎缩速度，m/s
h	——		传热系数， $W / (m 2 ? K)$
$h m$	——		传质系数， m/s
k	——		时间步长
$k f$	——		热导率， $W / (m ? K)$
L	——		特征尺寸
m	——		质量， kg
Pr	——		普朗特数
Q	——	种群的密度函数
R	——	气体常数
Re	——	雷诺数
T	——	温度， K
t	——	时间， s
X	——	液滴或颗粒的含湿量， kg水/kg固体
z	——	尺寸步长
ΔH	——	水的气化潜热， $k J / k g$
ρ	——	密度， $k g / m 3$
$? E$	——	干燥过程中的活化能， J/mol
上、下角标
b	——	干燥气体
c	——	平衡状态下的参数
d	——	液滴
s	——	表面
sat	——	饱和状态下的参数
t	——	时间
v	——	蒸汽
0	——	初始状态下的数值
L	——	特征尺寸

参考文献 25

1	赵志福, 朱宏吉, 于津津, 等 . 菊粉生产新技术研究进展[J]. 化工进展, 2008, 27(10): 1522-1532．
	ZHAO Zhifu , ZHU Hongji , YU Jinjin , et al ．New technical progress of the manufacture of inulin[J]. Chemical Industry and Engineering Progress, 2008, 27(10): 1522-1532．
2	郭义, 胡建华, 杨梓剑, 等 . 核黄素结晶母液的处理及回收利用[J]. 化工进展, 2016, 35(11): 324-327．
	GUO Yi , HU Jianhua , YANG Zijian , et al ．Treatment and recycling of riboflavin crystallization mother liquor[J]. Chemical Industry and Engineering Progress, 2016, 35(11): 324-327．
3	张丽丽 . 冲击式气流喷雾雾化机理及干燥过程数值模拟的研究[D]. 济南: 山东大学, 2008．
	ZHANG L L ．Study on atomization mechanism and numerical simulation of drying process of air-blast spray[D]. Jinan：Shandong University, 2008．
4	WAWRZYNIAK P , JASKULSKI M , ZBICIŃSKI I , et al . CFD modelling of moisture evaporation in an industrial dispersed system[J]．Advanced Powder Technology, 2017, 28(1): 167-176．
5	NANDIYANTO A B D , OKUYAMA K ．Progress in developing spray-drying methods for the production of controlled morphology particles：from the nanometer to submicrometer size ranges[J]. Advanced Powder Technology, 2011, 22(1): 1-19．
6	KEMP I C , WADLEY R , HARTWIG T , et al . Experimental study of spray drying and atomization with a two-fluid nozzle to produce inhalable particles[J]. Drying Technology, 2013, 31(8): 930-941．
7	ELVERSSON J , MILLQVIST FUREBY A , ALDERBORN G , et al ．Droplet and particle size relationship and shell thickness of inhalable lactose particles during spray drying[J]. Journal of Pharmaceutical Sciences, 2003, 92(4): 900-910．
8	ALI M, MAHMUD T , HEGGS P J , et al . A one-dimensional plug-flow model of a counter-current spray drying tower[J]．Chemical Engineering Research and Design, 2014, 92(5): 826-841．
9	LIU J J , MA C Y, HU Y D , et al ．Modelling protein crystallisation using morphological population balance models[J]．Chemical Engineering Research and Design, 2010, 88(4): 437-446．
10	FALOLA A , BORISSOVA A , WANG X Z ．Extended method of moment for general population balance models including size dependent growth rate, aggregation and breakage kernels[J]．Computers & Chemical Engineering, 2013, 56(56): 1-11．
11	LANGRISH TAG , KOCKEL T K ．The assessment of a characteristic drying curve for milk powder for use in computational fluid dynamics modelling[J]. Chemical Engineering Journal, 2001, 84(1): 69-74．
12	LIN S , CHEN X D . A model for drying of an aqueous lactose droplet using the reaction engineering approach[J]. Drying Technology, 2006, 24(11): 1329-1134．
13	WOO M W, DAUD W R W , MUJUMDAR A S , et al ．Comparative study of droplet drying models for CFD modelling[J]．Chemical Engineering Research and Design, 2008, 86(9): 1038-1048．
14	LIN S , CHEN X D ．Changes in milk droplet diameter during drying under constant drying conditions investigated using the glass-filament method[J]．Food & Bioproducts Processing, 2004, 82(3): 213-218．
15	HAR C L, FU N , CHAN E S , et al ．Unraveling the droplet drying characteristics of crystallization prone mannitol – experiments and modeling[J]．AIChE Journal, 2017, 63(6): 1839-1852．
16	WIJLHUIZEN A E , KERKHOF P J A M , BRUIN S ．Theoretical study of the inactivation of phosphatase during spray drying of skim-milk[J]．Chemical Engineering Science, 1979, 34(5): 651-660．
17	TOMINAGA T , MATSUMOTO S . Diffusion of polar and nonpolar molecules in water and ethanol[J]. Bulletin of the Chemical Society of Janpan, 2006, 63(2): 533-537．
18	FU N , WOO M W, SELOMULYA C , et al . Drying kinetics of skim milk with 50% initial solids[J]. Journal of Food Engineering, 2012, 109(4): 701-711．
19	LEVEQUE R J , MIHALAS D , DORFI E A , et al ．Computational methods for astrophysical fluid flow[M]. Berlin: Springer, 1998: 71-90．
20	WANG W-N , PURWANTO A , LENGGORO I W , et al ．Investigation on the correlations between droplet and particle size distribution in ultrasonic spray pyrolysis[J]．Industrial & Engineering Chemistry Research, 2008, 47(5): 1650-1659．
21	陈文武, 毕荣山, 刘振东, 等 . 气液喷射反应器内液滴粒径分布PLIF研究[J]. 化工进展, 2012, 31(4): 754-757．
	CHEN Wenwu , BI Rongshan , LIU Zhendong , et al . PILF study on droplet size distribution in gas-liquid jet reactor[J]. Chemical Industry and Engineering Progress, 2012, 31(4): 754-757．
22	郭金海, 谭心舜, 毕荣山, 等 . 压力旋流喷嘴雾化滴径分布的模型预测和实验[J]. 化工进展, 2012, 31(3): 528-532．
	GUO Jinhai , TAN Xinshun , BI Rongshan , et al ．Model prediction and experiment study on spray droplet size distribution of pressure swirl nozzle[J]. Chemical Industry and Engineering Progress, 2012, 31(3): 528-532．
23	KIEVIET F G . Modelling quality in spray drying[D]. Holland: Technische Universiteit Eindhoven, 1995．
24	WAWRZYNIAK P , ASKULSKI M , ZBICIŃSKI I , et al . CFD modelling of moisture evaporation in an industrial dispersed system[J]. Advanced Powder Technology, 2016, 28(1): 167-176.
25	ULLUM T , SLOTH J , BRASK A , et al ．Predicting spray dryer deposits by CFD and an empirical drying model[J]. Drying Technology, 2010, 28(5): 723-729.

[1]	王兰江, 梁瑜, 汤琼, 唐明兴, 李学宽, 刘雷, 董晋湘. 快速热解铂前体合成高分散的Pt/HY催化剂及其萘深度加氢性能[J]. 化工进展, 2023, 42(8): 4159-4166.
[2]	孙闫晨昊, 王维, 李一喆, 祝妍妮, 刘学武, 张大为. 橘皮油微胶囊制备及其产品质量评价[J]. 化工进展, 2023, 42(5): 2626-2637.
[3]	吴霞, 蒋勋涛, 张余晓, 吕慧园, 黄方, 于筱溪. 基于液滴微流控技术的蛋白质结晶[J]. 化工进展, 2023, 42(4): 2024-2030.
[4]	王光宇, 孟境辉, 张锴. 煤泥间歇微波干燥模拟及介电性质[J]. 化工进展, 2023, 42(4): 1779-1786.
[5]	刘畅, 李玉宝, 左奕, 李吉东. 羟基磷灰石负载黄芩苷的不同载药方法比较分析[J]. 化工进展, 2022, 41(9): 4995-5002.
[6]	杨磊, 宋金玲, 唐初阳, 于诗尧, 杨欣宇. 基于FUSION模型的碳基固废热解产物产率预测[J]. 化工进展, 2022, 41(7): 3966-3973.
[7]	姜华, 张子惠, 宫武旗, 常越勇. MVR分质提盐蒸发结晶系统设计及性能分析[J]. 化工进展, 2022, 41(7): 3947-3956.
[8]	李娟娟, 张天永, 李祥高. 电泳显示用高质量铁锰黑纳米分散系的构建[J]. 化工进展, 2022, 41(6): 3178-3185.
[9]	周亚兰, 闫雯, 罗路, 范毜仔, 杜官本, 赵伟刚. 酚醛基炭气凝胶的研究进展[J]. 化工进展, 2022, 41(4): 1970-1981.
[10]	武晨浩, 李昆锋, 李肖华, 费志方, 张震, 杨自春. 二氧化硅气凝胶常压干燥工艺的研究进展[J]. 化工进展, 2022, 41(2): 837-847.
[11]	韩芬, 杨娜, 孙永利, 姜斌, 肖晓明, 张吕鸿. 玻璃纤维聚结器脱除油中乳化水[J]. 化工进展, 2022, 41(12): 6723-6732.
[12]	张爽, 赵立新, 刘洋, 宋民航, 刘琳. 脱气除油旋流系统流场分布及分离特性[J]. 化工进展, 2022, 41(1): 75-85.
[13]	贺远鹏. 响应曲面法复配混凝剂处理海上生产返排液[J]. 化工进展, 2021, 40(S1): 439-445.
[14]	赵宁, 冯永新, 林廷坤, 谢志文. 热烟气环境下脱硫废水液滴蒸发特性研究进展[J]. 化工进展, 2021, 40(8): 4508-4514.
[15]	李金, 王军锋, 徐惠斌, 郑高杰, 孟新. 电晕放电喷雾荷电特性[J]. 化工进展, 2021, 40(3): 1300-1306.

甘露醇喷雾干燥过程中液滴粒度分布变化的群体粒数衡算模拟和实验研究

Simulation and experimental study on the evolution of droplet size distribution during spray drying of mannitol

RichHTML

PDF (PC)

可视化

摘要/Abstract

引用本文

使用本文

图/表 10

参考文献 25

相关文章 15

编辑推荐

Metrics

本文评价