图1 充氮工艺示意图
图2 灌装机三维模型
其中,ρ 为密度,kg/m 3 ;u、v 和 w 为速度矢量在 x、y 和 z 三个方向上的分量,m/s;t 为时间,s。流体在 x、y、z 三个坐标轴方向上的动量守恒方程表达式为:
其中,ρ 为密度,kg/m 3 ;t 为时间,s;u 是速度矢量,m/s ;ρ 是流体微元体上的压力,N ;μ 是动力黏度,N · s/m 2 ;Su 是源项;i 为 x,y,z。
其 中,t 为 时 间,s ;V 为 体 积,m 3 ;a 为面积矢量,m 2 ;i 为组分指数;ρ 为总体密度,kg/m 3 ;v 是速度,m/s ;u t 为湍流动力粘度;σ t 为湍流施密特数;S Y i 为指定区域源项;J i 为层流扩散通量。(2)湍流模型根据设备实际运动情况,本文选 择基于湍动能k及耗散率ε 的标准k-ε湍流模型进行分析。其中,k 和 ε 通过以下方程获得:
式中,G k 为由于平均速度梯度产生的湍流动能生成项;G b 为浮力引起的湍流动能生成项;Y M 表示可压缩湍流中脉动膨胀对总耗散率的贡献;C 1ε 、C 2ε 和 C 3ε 为常数;σ k 和 σ ε 分别是 k和 ε 的湍流普朗特数;S k 和 S ε 为用户定义的源项。
图3 氮气浓度分布云图图4 展示了压塞工位瓶中氮气浓度变化曲线,由图可知,在极短时间内,瓶内氮气会有少量溢出,在氮幕的作用下,对瓶内氮气进行快速补充,氮气浓度升高,其中,靠近后充氮工位的两个西林瓶瓶内氮气浓度升高的较快,远离此工位的瓶子内部氮气浓度升高较慢。结合图3氮气浓度分布云图,2s 后,靠近后充氮工位的西林瓶瓶内氮气浓度最高,最高体积占比为 93%,远离后充氮工位的西林瓶瓶内氮气浓度最低,最低体积占比为 82%,也就是说空气占比约 7% 和 18%。考虑到氮气浓度越低药液被氧化的风险越高,以远离后充氮工位的西林瓶顶空残氧量作为分析依据,众所周知空气中存在 21% 的氧含量,在不考虑药液中溶解氧因素的影响,可知此充氮系统可满足的最大顶空残氧量约为 3.78% [5] 。很明显此设计方案距离国际先进水平(残氧量< 1%)还有很大的差距。
图4 压塞工位瓶中氮气浓度变化曲线
图5 速度分布云图b、由于加塞孔为非封闭环境,且充氮孔的倾斜角度较大,导致大量氮气在短时间内流失到加塞孔外,使得加塞孔内缺乏足够的氮气;c、氮幕机构与胶塞、瓶口之间留有的配合间隙较大,导致压塞区域的氮气流失严重。参考以上因素并结合实际结构的复杂程度,可提出以下优化方案,如图6所示。在原设计方案的基础上,适当增加充氮孔的数量,增大充氮孔的直径,并调整充氮孔的充氮角度。同时,合理控制氮幕机构与瓶口、胶塞之间的配合间隙,增加单位时间进入加塞环境的氮气含量,并改善加塞区域的氮气分布。
图6 优化方案在与原设计方案相同的仿真条件下,对优化方案进行了 CFD 流体分析。图7 和图8 分别展示了仿真得出的优化方案的氮气浓度分布云图及压塞工位瓶中氮气浓度变化曲线图。从图中可以看出,在充氮时间内,压塞区域的氮气浓度快速升高并置换瓶中的空气,且压塞工位不同位置的西林瓶中氮气含量较均匀。瓶内氮气的最大浓度为 99.6%,最小为 98.8%,即空气占比约为 0.4% 和 1.2%,最大顶空残氧量约为 0.252%。通过对比图 5 和图 9 方案优化前后的速度分布云图,可以发现优化方案在瓶口上方加塞区域、顶空的气流速度大于原设计方案,从而可以快速置换顶空及瓶口附近的空气,短时间内提高压塞区域的氮气浓度,充氮效率得到提升。
图7 优化方案氮气浓度分布云图
图8 优化方案压塞工位瓶中氮气浓度变化曲线
图9 优化方案速度分布云图
图10 实验现场图| 欢迎光临 环球过滤分离技术网(http://guolvfenlitech.com)是中国过滤.分离.净化.纯化、超洁净技术行业领先的前沿科技、产品和技术信息供求,技术解决方案,技术交流分享以及过滤、分离、纯化、超洁净产品和技术供需服务平台。 (http://www.guolvfenlitech.com/) | Powered by Discuz! X3.3 |