金银花提取液喷雾干燥过程的数值模拟
- 时间：1周前
- 作者：产品展示
与其他干燥过程相比，喷雾干燥过程中粉末与热空气之间的接触时间更短，可以有很大效果预防物料因加热时间过长而变性。尤其对于中药制药过程，粉体是一种很重要的物料形式，如中药配方颗粒的所需要的中药浸膏粉往往通过喷雾干燥获得。但是由于中药种类及来源的复杂性，目前业界对中药喷雾干燥中产品质量形成与传递的规律认识不够清楚，导致质量控制不精准，缺乏基于模型预测的前馈控制与优化方法，喷雾干燥过程的工艺调整主要依赖经验。研究各工艺变量对粉末特性的影响有必要进行实验室和中试试验，这需要消耗大量的时间和人力成本
使用计算流体力学模拟对液滴与干燥空气之间的复杂相互作用、液滴与壁面相互作用以及液滴和干燥空气之间的传热传质过程进行研究，能大大的提升我们对干燥塔中喷雾干燥过程的理解[11-12]。利用数值模拟能获得实验中难以直接测量的数据，并预测实验过程中难以直接观察到的现象，这将有利于干燥塔的更有效的设计和操作，并提升产品质量[13]。目前，还没有关于中药喷雾干燥过程机理的研究，因为中药与化学药或生物药相比，不同之处在于中药材是天然来源的，受到种植条件及采收时间的影响，导致同一中药材的批次间差异性较大，从而造成生产的全部过程中提取液的物性波动，不能通过简单的干燥动力学模型进行研究[14-16]。本实验采用欧拉-拉格朗日框架对连续相和离散相进行模拟计算，建立了金银花提取液喷雾干燥计算流体力学模型。通过分布式光纤测温系统对干燥塔内的温度场模拟进行了验证，并分析了不同初始粒径液滴的运动和干燥特性，为解释中药提取液喷雾干燥机制和工艺参数优化提供理论参考。
H-spray 5S中试型喷雾干燥设备，天津市现代创新中药科技有限公司；Mastersizer 3000型激光粒度仪，英国马尔文仪器有限公司；NKT Photonics型分布式光纤测温系统，德国LIOS Technology GmbH公司。
金银花药材，河北省安国市腾跃药业有限公司，批号211010，经天津中医药大学中药学院张丽娟教授鉴定为忍冬科忍冬属植物忍冬Lonicera japonicaThunb.的干燥花蕾或带初开的花。
2.1.1金银花提取液的制备金银花药材通过冷凝回流方式来进行提取，料液比1∶10（金银花-水），提取2次。第1次提取时间为2 h，第2次提取时间为1.5 h。2次提取结合使用旋转蒸发器进行浓缩，即得。
2.1.2模型描述金银花提取液喷雾干燥工艺流程图如图1所示。使用中试型喷雾干燥机进行实验，干燥空气被加热后由干燥塔顶部分风器进入塔内。金银花提取液通过蠕动泵进入双流体雾化喷嘴，被压缩气体分散小液滴喷入干燥塔。干燥热空气与液滴从塔顶并流向下流动，进行传热传质。干燥完成后，干燥空气和金银花颗粒在旋风分离器中完成气固分离。较大的颗粒将沉积在旋风分离器下的接收器中，极细的粉尘和空气被排出干燥塔。
2.2.1理论模型在本研究中，使用欧拉-拉格朗日模型定义了相关的守恒方程[17]。热空气和干燥的颗粒分别被视为连续相和分散相。
（1）连续相控制方程：在欧拉-拉格朗日框架中对连续相进行建模时，连续相的动量、能量和质量的平衡方程如下。
ρC为连续相密度，t为时间，UC为连续相速度矢量，Γm,CD为热空气与液滴相互作用产生的源项
ui、uj为平均速度分量，p为压力，g为重力加速度，FCD为液滴与热空气的相互作用力，τ是黏性应力张量
有效黏度（μeff）包括运动黏度（μ）和湍流黏度（μt），定义为μeff＝μ＋μt。
kt为热导率，TC为连续相温度，QC为连续相传热，Γe,CD为热空气与液滴的能量源项，hC为传热系数，该系数是从Ranz-Marshall相关性获得的。
dD为液滴直径，k∞为连续相的导热系数，ReD基于粒子直径和相对速度的雷诺数，Pr为连续相的普朗特数
k为湍流动能，xi,j为x、y方向上的单位向量，ω为比耗散率，β为热线胀系数，Gk代表产生的湍流动能，Gω代表ω的生成
（2）离散相控制方程：在离散相模型中，粒子的运动行为由力平衡方程计算得出。基于牛顿第2定律，单粒子运动方程构建为
UD为颗粒/液滴的速度，ρD为颗粒/液滴的密度，FCDdrag被定义为单粒子曳力函数，并用于确定方程（2）中牛顿第3定律上离散相和连续相的相互作用
其中a1、a2和a3是依赖于粒子雷诺数的常数，适用于粒子雷诺数，范围0～50 000
（3）液滴干燥动力学：每个液滴/颗粒内的内部传输现象由Methericher等[18-19]先前开发和验证的2级干燥动力学模型描述。含固体液滴的干燥过程分为2个干燥阶段：在干燥的第1阶段与纯液体蒸发类似。在第2个干燥阶段，干燥速率由水分从颗粒湿核通过外壳孔向颗粒外表面扩散的速率控制。这样的一个过程中湿颗粒的核心收缩，颗粒外壳区域的厚度增加。由于对干燥塔内的液滴做测量难度较大，本研究使用激光粒度仪测定，通过单液滴干燥实验反推得到液滴的初始粒径。
本实验采用特征干燥率曲线（characteristicdrying rate curve，CDC）模型计算液滴的蒸发速率[20]。CDC方法基于2个不同干燥期的假设，即恒定速率周期和下降速率周期。当外壳出现在液滴表面上时，恒定速率周期结束。此时，临界水分（Xcr）阻碍了下降速率周期中的水分去除过程。经验无量纲函数f根据颗粒水分含量X定义，根据课题组前期研究基础[13]，参数n的推导值为1.81。在计算模型中，Xcr设置为1.42。实验表明，所有悬浮液滴的最终湿度接近于零，因此，假设方程（1）中的平衡湿度为Xeq=0。
根据方程（6）和方程（13）中，传质方程用于模拟热空气条件下的单颗粒蒸发速率。
通过将Ranz-Marshall相关性乘以经验函数f来考虑饲料配方的干燥动力学[21]，即
为了确定两相之间传热所改变的能量，可以用方程（16）计算单个粒子的内部能量方程。
2.2.2模型设置在温度为448 K和绝对湿度为2.93 g/kg的干燥空气的情况下，通过干燥塔顶部的中央圆形入口以相对于垂直轴15°的角度向干燥室提供预热的空气。进风速度为6 m/s，湍流动能为0.043 2 m2/s2，比耗散速率等于180.69 s−1。液滴的喷射是通过在位于塔顶的压力喷嘴中雾化液体进料获得的。料液的组成为85%的水和15%的可溶性固形物。通过对喷雾过程进行拍照，确定喷雾锥角为30°。喷嘴出口处的液滴速度被指定为65 m/s，进料温度设置为293 K。相应的喷雾流量等于2.145 L/h。通过热平衡实验确定，从干燥塔到环境的热损失参数为1.3 W/(m2K)。干燥室出口管中的表压设置为−100 Pa。
最初模拟空气以使用耦合求解器实现稳态，随后注入粒子以达到稳定状态。仿真结果采取了压力-速度耦合算法求解。对动量、体积分数、湍流动力学和湍流耗散方程应用了二阶迎风格式。压力插值由PRESTO!求解（Fluent 2021，ANSYS，美国），通过有限体积法来离散化流场，然后使用离散连续平衡方程来计算交错控制单元的交错压力。
将喷雾干燥的实测值与数值模型的计算结果作比较，从而验证模型的准确性。实验中发现从旋风分离器逃逸的颗粒几乎能忽略不记，计算旋风分离器下端收集到的颗粒与进料量之比得到得粉率，平均粒径通过激光粒度仪获得，数值模型通过统计出口的颗粒质量及粒径从而计算得到得粉率和平均粒径的模拟值。
如表1所示，出口温度、得粉率和颗粒平均直径均在5%以内，证明模拟结果与实际的喷雾干燥现象基本一致。本研究采用分布式光纤测温系统测定了干燥塔内的温度分布，得到干燥塔3个高度处的温度分布，如图2所示。分布式光纤测温系统的测量精度为0.15 m，数值模型的计算精度0.01 m，因此测量节点之间的温度剧烈变化系统没有办法获得，但是实验与模拟大部分区域的温度分布（在0.18、0.48、0.78 m高度处）趋势均较为接近。根据实验与模拟对比，可以在某些特定的程度上反映喷雾干燥塔内部液滴的传热传质过程。
图3为YZ平面不同高度处的空气速度分布。干燥塔中心的空气速度明显高于周围，中心处越远离入口空气速度越低。
如图4、5所示，在实验和模拟中，动量、质量和能量传递方程之间有相互作用关系。当热能被液滴蒸发带走时，空气速度能由周围的热空气流场补充。随着蒸发速率的变化，空气温度和速度有明显的不规则分布。本课题组已经根据分布式光纤测量结果验证了数值模型的温度场，然后可以合理地认为模拟计算的速度场是可靠的[22-23]。从图4和图5中能够准确的看出，干燥塔的大部分区域具有几乎恒定的温度和速度。如图4所示，在塔的上部中心区域附近有一个低温区，周围有一个高温区，并且中心外部区域中的气流速度很低。低温区为推断喷雾干燥过程中液滴的运动特性和相变位置提供了证据，这对于工艺优化和设备设计至关重要。
图6是随粒子停滞时间变化的各尺寸粒子（初始粒径）的归一化速度分布图。所有注入液滴的初始速度相同。与较大的颗粒相比，较小的颗粒失去初始速度更快，因为较小的颗粒具有较小的质量，因此惯性较小，阻力的影响更大，因此较小的颗粒向下移动时减速更快。将颗粒的最小速度认定为最终下降速度，当粒子达到其最终下落速度时，它们继续以该最终速度下落。与较大的颗粒相比，较小的颗粒具有较大的最终下落速度。考虑到粒子与壁面的碰撞和沉积作用，较小的粒子有更大的可能性在塔内壁沉积。
图7是不同颗粒尺寸随粒子停滞时间的温度分布图。当液滴从雾化器喷射时，浆液液滴的初始温度达到甚至略高于湿球温度。在此阶段，较小的液滴表面发生的蒸发大于向液滴的热传递。因此较小尺寸的液滴，温度开始下降。较大尺寸的液滴初始阶段升温较慢，最终滴达到一个恒定湿球温度，随着粒子直径的增大，恒定湿球温度阶段也越长。在某一点之后，颗粒的温度开始快速上升。这是由于液滴表明产生了结壳。所有液滴尺寸都观察到这种趋势。与较大液滴尺寸相比，最小液滴的结壳形成时间要短得多。一旦结壳形成，液滴的水分损失率就会降低，导致与水分蒸发释放的热量相比，向液滴的热量传递更大。最小粒径（3.384 µm）迅速达到平衡温度，随后基本保存不变，较大的粒子达到平衡温度需要更长的时间，最终粒子保持不变证明粒子的中的水分已完全干燥。
图8是不一样的尺寸粒子的传热系数图。较小尺寸的粒子（＜20 µm）的传热系数在顶部最小，并随停滞时间增大而减小。较大尺寸的粒子的传热系数先减小后增大，最后保持恒定。对于最小的液滴，传热系数的初始值最大，因为它与直径成反比。由于速度的迅速下降，较小的粒子传质系数迅速增大，并最终随着粒子速度达到最终下降速度而变得恒定。与所有较大尺寸相比，较小粒子最终传质系数仍然是最高的，最大尺寸（63.540 µm）粒子的传热系数最小。
图9是不一样的尺寸的液滴/颗粒的归一化含水率分布图。较小的颗粒比较大的颗粒损失水分更快，这是由于较小尺寸的较大的比面积。颗粒经历水分蒸发曲线个不同斜率，代表2个干燥阶段。这些在较大尺寸颗粒（＞20 µm）的含水率变化中更明显。在第一阶段，含水率变化曲线的斜率几乎是线个阶段的恒速干燥阶段。含水率下降的速度随后变慢，这对应于液滴的第2个阶段的限速干燥阶段。最后当水分蒸发率随着颗粒水分含量接近零和温度接近气体温度而降低，因为该阶段颗粒和气体之间的温差作为蒸发的驱动力变小了。最终所有直径的颗粒在离开塔之前会失去所有水分。由于粒子表面有水分形成的水膜，水分快速蒸发，导致粒子直径在前期出现了明显降低，而干燥后期水分降低较慢，粒子直径变化较不明显，这与干燥的2个阶段的过程是类似的（图10）。
温度、传热系数、含水率和直径的分析，发现干燥塔内的温度分布与液滴的干燥和运动特性紧密关联。干燥塔中测得的温度与模拟一致，可以很好地解释干燥过程中液滴的变化。干燥塔中有3个不同的气温变化区，第1个是雾化区，由于蒸发去除水分，液滴尺寸也强烈收缩。水相变吸收大量热量，并且该区域的液滴温度上升得更慢。第2个区域处于液滴的第2个干燥干燥阶段，干燥速率由于液滴的表面结膜而降低，伴随着液滴表面温度的快速升高。最后当液滴完全干燥时，观察到温度和湿度水平稳定在干燥塔出口上方。粒子的初始粒径不同，其运动和干燥速度也不同，但是基本遵循相似的变化规律。4讨论
拉格朗日模型，以计算干燥塔中的瞬态气流、传热和传质以及颗粒运动。模拟所得结果与实验数据基本一致，最大相对误差为4.34%，证明了数值模型的准确性。在实验中使用分布式光纤测温系统测量了干燥塔的温度场，干燥塔内温度分布的模拟和实验证明了干燥塔上部存在一个低温区，粒子在此区域发生剧烈的传质传热反应。通过对不同初始粒径的速度、温度、传热系数、含水率和直径的分析，分析了液滴的蒸发和运动特性。本实验建立的数值计算模型模拟了喷雾干燥塔内的复杂气流模式和中药提取物喷雾的干燥历史，与普通建模方法相比，数值计算模型可以基于产品质量和投料物性快速优化生产的基本工艺，减低试错成本，未来能够适用于中药生产设备和工艺的优化。利益冲突
来 源：崔彭帝，薛启隆，程宇康，刘新龙，于 洋，李 正．金银花提取液喷雾干燥过程的数值模拟 [J]. 中草药, 2023, 54(15):4832-4838.
