功能化涂层的制备在血管支架开发过程中发挥着关键作用, 其技术难点在于对微米级的支架结构实现对涂层总量的精确控制并确保涂层的稳定和均匀. 对血管支架涂层制备技术进行系统的研究, 不仅可以指导血管支架涂层制备的过程, 合理的控制涂层的质量, 还可以实现涂层形貌的优化, 从而提高支
架的生物医学性能.
2005 年, 超声雾化喷涂技术开始应用于血管支架涂层制备, 与浸涂法(Dip-Coating)
[1]、包覆法(Wrap-Coating)[2]、气相沉积法(Gas-phase Deposition)[3]等相比, 具备涂层制备高效、涂层质量可控、喷涂方式灵活等优势. 目前研究者更多关注涂层材料的生物医学性能[4~6], 较少关注血管支架涂层存在的形貌特征缺陷(图 1), 所制备的涂层质量千差万别[7~9].随着技术发展, 一些新型超声喷涂系统陆续出现, 对不同的超声雾化喷涂系统, 其雾化形式、喷雾特性及应用范围不同, 尤其在血管支架这种对涂层质量要求很高的微观场合, 其工艺制定过程较为复杂, 参数选型亟需可靠依据, 同时对超声雾化理论机制与微观过程研究较少, 理论体系尚不完善[10~12].
针对上述问题, 本文提出对制备血管支架表面涂层的超声波喷涂工艺进行系统研究. 研究推导超声波雾化粒径方程, 并基于计算流体动力学仿真(CFD)研究超声雾化的微观过程, 分析影响超声雾化粒径的关键参数; 在上述基础上, 研究构建实验检测喷雾粒径尺寸的平台, 实现喷雾微粒测量, 验证数值分析的正确性, 最后依据理论和数值研究成果, 给出制备优良血管支架超声喷涂涂层制备的工艺参数,并进行喷涂试验验证, 解决支架喷涂中存在的质量缺陷问题, 为制备涂层稳定和均匀的血管支架提供技术支撑.
1 超声波雾化的物理模型
超声波雾化是指液膜在高频振动的激励下, 形成小液滴脱离液面的现象. 超声波雾化粒径大小是表征雾化性能的重要指标. 本文采用数量平均直径和索特平均直径描述喷雾粒径大小, 其表达式如下:
代表喷雾中任意粒子直径.以整个液膜作为研究对象, 其中雾化前能量(Pf)包括: 振动能(Pvf)、内能(Pif)、表面能(Psf)以及势能(Ppf); 雾化后能量(Pl)包括: 液滴动能(Pvl)、内能(Pil)、表面能(Psl)以及势能(Ppl). 根据能量守恒原理, 雾化前后总能量不变, 则有.vP(3)将液膜看作弹簧振子, 当毛细波角频率为?, 振幅为 A 时, 假设雾化后的雾滴半径为 rx, 单位时间雾化的液膜质量为 mf, 则雾化前振动能 Pvf满足表面能方面, 雾化前液膜表面能Psf为0; 雾化后根据表面能方程, 喷雾总表面能(P
sl)计算式为超声雾化过程中, 液膜内能变化与温度有关, 课题组在进行超声雾化试验过程中, 用红外温度测试仪检测液体温度(Tl)和喷嘴温度(Tn)随时间变化差异较小, 忽略雾化前后的内能变化, 即 P
if≈ Pil. 依据Sindayihebura 等人[13]的实验结果, 雾化粒子速度很小, 忽略粒子动能, 即Pvl=0. 同时根据 Edwards 等人[14]的研究, 超声雾化毛细波波数较大, 可忽略重力影响,即 Ppf= Ppl.综上, 推导索特平均直径 d32满足其中??为表面张力, 为密度, f 为超声波频率. 早期超声雾化研究中, Lang[15]依据 Kelvin 方程[16], 在实验基础上提出了超声激励毛细波波长(c)公式:1/ 328π.cf?(7)在此, 假设液膜雾化瞬间振幅与毛细波波长成正比[17], 比值为, 则 d32的表达式为(8)可见喷雾平均粒径与液体表面张力、液体密度、超声波频率以及毛细波波长有关. 借助 Lacas 等人[18]的实验数据, 利用最小二乘法拟合式(8)曲线如图 2所示, 此时 const=5.6986, ??=0.231, 索特平均直径满足1/ 3
(9)式表明超声雾化粒子尺寸与毛细波波长成正比. 依据(9)式绘制出图 3 不同液体喷雾粒径随超声波频率的变化曲线, 其中水的表面张力与密度比最大, 曲线位于最上方, 随着超声波频率的增加, 粒径减小速度趋于缓慢. 从喷涂工艺角度, 当喷涂液体系一定时, 频率调节相对表面张力更加灵活, 在硬件条件允许的情况下, 可以通过提高超高频率获得细小雾滴.
的圆盘, 圆盘内液膜厚度为 50 ?m, 上方为高 150 ?m,直径1 mm的空气区域, 如图4所示. 工作时, 沿轴向圆盘处于高频振动状态, 驱动液膜表面形成毛细波并产生雾化液滴.对于三维轴对称流场,在 Gambit 中建立长 1 mm,宽 0.2 mm 的二维流域, 采用三角形单元进行网格划分, 定义底部和侧面为壁面, 上部为压力出口, 计算域如图 5 所示.考虑到存在空气和液体两相, 采用多相流(VOF)模型, 设置壁面对液体无滑移剪切作用. 在设置振动边界时, 采用 Fluent动网格模型, 将简谐振动的 UDF函数加载至整个流体区域. 流场底部受简谐激励, 轴向速度 v(x, t)满足其中, ??(x)为端面振幅. 为研究不同参数对超声雾化的影响, 拟定喷嘴端面振幅分别为 2, 4, 8, 10, 12 ?m;基于常用溶剂如水、四氢呋喃(THF)等物理学参数拟定差异化的表面张力和动力黏度; 选择 20~120 kHz范围制定超声频率. 仿真结束后, 对毛细波波长、雾化粒径以及临界振幅进行测量统计, 结果如附表 1 所示
2 超声波雾化的 CFD 模拟与分析
在超声喷涂工艺中, 超声雾化过程决定了粒子的初始尺寸, 也决定其是否满足支架的微尺度空间要求. 本小节利用CFD 模拟分析超声波频率、雾化功率、液体表面张力和黏度等参数对超声雾化的影响.
2.1 计算模型
以喷嘴端面为研究对象, 参考 Barreras 等人[19]的超声雾化实验装置, 设计直径为1 mm, 高度为200m
2.2.2 雾化功率
通过计算流场模型获得不同振幅毛细波形态如图 8 所示: 当振幅为 2 ?m 时, 液膜表面无毛细波出现, 随着振幅增加, 毛细波波高增加, 形成圆射流.根据附表 1 结果, 毛细波波长与超声雾化功率无关,始终维持在39 ?m左右; 为了进一步分析功率对雾化粒径的影响, 绘制雾化粒径随振幅变化的关系如图 9所示, 随着振幅的增加, 雾化粒径先小幅降低, 后保持不变. 综上所述, 振幅是决定超声雾化能否进行的关键因素, 当超过临界振幅后, 其数值大小对雾化粒径影响不大.分析不同振幅下超声波雾化过程, 当振幅为 4 ?m时, 毛细波未达到临界振幅, 有少量粒子出现. 定义低于临界振幅的雾化为亚雾化现象(Incomplete AtomizationPhenomenon, IAP), 如图 10 所示, 亚雾化粒子形成跨两个振动周期, 初期为纺锤体形状, 未出现主副粒子分离现象, 粒子直径较正常值大. 如果以沿 Y 轴正向运动为正,亚雾化粒子速度矢量为负值, 雾化后粒子停留在液面附近, 部分被后续增长的波峰吞噬. 亚雾化现象是超声波雾化不完全的体现, 其本质是毛细波动能不足, 当振幅增大后, 该现象消失. 超声雾化中出现的亚雾化现象使雾化粒子平均直径增大, 粒子数量减少, 影响喷雾均匀性, 和Tsai等人[20]的实验研究结果一致. 因此, 合理控制输入功率的大小, 保证其高于临界振幅是获良好超声喷雾的关键.2.2.3 表面张力根据 CFD 模拟结果, 不同表面张力对应毛细波形貌如图 11 所示, 随着表面张力的增大, 毛细波波
长增加. 根据测量相分布获得的不同表面张力下的临界振幅(见附表 1), 随着表面张力增加, 超声雾化
临界振幅增加. 此外, 雾化粒径随着表面张力的增加显著增大, 如图 12 所示.通过对表面张力为0.015 N/m时毛细波雾化过程进一步分析, 结果如图 13 所示, 柱形波峰表面形成对称波纹, 其端部首先分离形成主液滴, 中间部分随之碎裂成许多小液滴, 直径大约为主液滴直径的一半.决定表面波碎裂成小液滴的关键因素之一是碎裂长度, 在超声雾化中即射流雾化所需的波高, 根据Miesse[21]给出的紊流雾化区碎裂长度公式, 当射流直径为 dl时, 推导超声波雾化条件下表面张力波断裂长度(Lu)满足
2.2.4 液体黏度
根据CFD模拟结果, 不同黏度下毛细波在t=38 ?s时刻波形如图 14 所示, 随着黏度增加, 毛细波波长基本不变, 波高不断降低. 当黏度较小时, 同一时刻到达临界振幅的粒子数较多, 随着黏度的增加, 毛细波动能不足, 达到临界振幅的波峰数减少, 出现亚雾化现象. 根据相分布测量结果(附表 1), 随着黏度的增加, 临界振幅保持在 40 ?m 左右, 表明黏度对临界振幅没有明显影响.
在 CFD 仿真计算中, 不同黏度对应的雾化粒径如图 15 所示, 随着黏度增加, 雾化粒径整体上呈微弱增大的趋势, 这一点与黏度变化导致亚雾化现象增多趋势一致; 本质上讲, 粒径增大是由于雾化粒子中亚雾化粒子比例的增加导致的.
2.3 基于 CFD 仿真结果的参数分析与拟合
2.3.1 雾化粒径
根据参数影响规律, 比较理论值和仿真测量数值如表 1 所示, 其中BOE, d32_BOE分别为(7)和(9)式(Based on known equation, BOE)预测的毛细波波长和雾化粒径, ?CFD, d32_CFD分别为 CFD 仿真获得的毛细波波长和雾化粒径. 根据表 1 中参数, ?CFD≈?BOE, 而d32_BOE与 d32_CFD存在较大差异.
在 1.2 小节中, 推导出在黏度造成的内部剪切效应影响较小的情况下, 雾化粒子索特平均直径与毛
为了验证这一结论, 绘制?CFD-d32_CFD关系如图 16 所示, 表面波波长与雾化粒径呈近似线形关系, 利用 d32_CFD=k′?CFD拟合该曲线, 利用最小二乘法求得 k′=0.3, 即 CFD 模拟条件下超声
雾化粒子索特平均直径(d32_CFD)方程为1/ 332_CFD CFD28πd0.3 0.3 .f
2.3.2 临界振幅与功率
为探求临界振幅和毛细波波长可能存在的规律,绘制不同参数对应的毛细波波长和临界振幅Amcrit-?CFD关系图(图 17), 采用多项式 y=ax2+bx 进行曲线拟合, 利用最小二乘法求得 a= ?0.0033, =0.860, 则临界振幅 Amcrit满足
Amcrit =CFD CFD2 / 3 1/ 32 20.0033 0.868π 8π0.0033 0.86 .Af f
2.3.3 喷雾流量
根据超声波雾化的微观模式, 定义三个雾化区:亚雾化区、理想雾化区、射流雾化区, 如图 18 所示,
当功率过低时, 进入截至区, 超声雾化无法形成; 亚雾化区和理想雾化区以临界功率(Pi)作为分界点, 随着功率的增加, 雾滴直径逐渐减小, 标准偏差(StdDev)先逐渐增大, 在射流雾化区到达最大值, 后逐渐减小. 当超声喷嘴表面毛细波处于理想雾化区时, 粒子雾化频率与超声波频率相等. 假设毛细波沿纵向和横向均匀分布, 形成毛细波矩阵, 如图19所示, 其中圆圈代表波峰, 每两个波峰之间间隔为一个波长,单位面积波峰总数量 N 为?
3.2 喷雾粒径尺寸与分布规律
根据图像采集结果, 统计数量平均直径(d10)并计算索特平均直径(d32), 以标准偏差(Std Dev)表征粒
子分布均匀性, 结果如附表 2 所示. 可见超声雾化索特平均粒径在9.71~13.14 ?m 范围, 平均粒径为d32_exp=10.71 ?m, 根据(9)和(12)式, 求得索特平均直径 d32_LaF=8.5 ?m, d32_CFD=11.33 ?m, 从数值大小看, d32_CFD更接近于实际值.
采用 Rosin-Rammler 函数[22]拟合粒子尺寸分布曲线, 表达式为
3.1 实验方案设计
采用离线粒径检测方案, 流程如图 20 所示, 实验参数组如附表 2 所示, 采用 MediCoat 超声波喷涂系统(SONO-TEC, USA), 超声波频率为 120 kHz. 采用瞬间冷凝法收集喷雾微粒, 雾化液采用溴代十六烷(C16H33Br), 冷凝液为低温去离子水, 在烧杯外部包覆保温层作为采集器. 最后, 采用 Zoom645 体式显微镜和 CV-A1 扫描 CCD 进行图像采集, 并对粒子图像进行测量和统计。
3.2 喷雾粒径尺寸与分布规律
根据图像采集结果, 统计数量平均直径(d10)并计算索特平均直径(d32), 以标准偏差(Std Dev)表征粒子分布均匀性, 结果如附表 2 所示. 可见超声雾化索特平均粒径在9.71~13.14 ?m 范围, 平均粒径为d32_exp=10.71 ?m, 根据(9)和(12)式, 求得索特平均直径 d32_LaF=8.5 ?m, d32_CFD=11.33 ?m, 从数值大小看, d32_CFD更接近于实际值.
4 血管支架表面涂层超声喷涂工艺试验
采用本实验室开发的 S 型钴铬合金支架, 基于Medicoat 喷涂系统, 以 30%PEO-PSF 为涂层材料,以四氢呋喃作为溶剂, 基于上述实验研究参数方程, 确立包括超声波频率、液体流量、输入功率、辅助 气体压力以及喷涂液浓度等参数, 如表 2 所示.依据表 2 参数, 进行支架的喷涂实验, 并对喷涂后的支架利用 Quanta 200 环境扫描电镜(ESEM, FEI)进行观测, 如图 31 所示, 当放大倍率为 100×时, 支架表面没有出现厚边、桥接等缺陷, 宏观质量较优;当放大倍率为 200×时, 支架外表面涂层均匀, 内表面涂层与外表面涂层保持一致; 当放大倍率为 2000×时, 可见支架 S 形位置涂层光滑平整, 没有出现凹凸、橘皮、缩孔等缺陷, 微观质量非常出色. 综合上述情况表明, 优化后的支架涂层表面质量优, 喷涂率高, 可靠性强, 完全达到了优化的预期目的.
为了进一步检测支架的纳观形貌, 采用原子力显微镜(AFM, Asylum Research)扫描喷涂涂层表面,如图 32 所示, 分别为电化学抛光后的裸支架、以及超声波喷涂后涂层在 2 ?m × 2 ?m 的表面形貌, 根据测量结果, 抛光后裸支架金属表面粗糙度大约为12.4 nm, 超声喷涂获得的涂层表面粗糙度大约为1.72 nm, 涂层质量优良, 达到了较高水准, 涂层表面比喷涂前的金属表面光滑得多, 表明通过控制喷雾质量可以大大提高获得的涂层形貌质量, 实现显著提高血管支架涂层血液相容性的目标.
5 总结
本文从理论模型建立、流体力学仿真、喷雾粒径检测以及超声喷涂试验四个方面研究了制备性能优良的血管支架涂层所涉及的工艺参数. 定义了面向常温低黏度流体超声雾化过程的亚雾化、理想雾化和
射流雾化三种模式. 解决了 Pohlman 方程对水、四氢呋喃等常温低黏度流体超声雾化振幅方程预测不准的问题. 拟合了面向常温低黏度流体超声波雾化的临界振幅方程, 并在此基础上建立了理想雾化模式
下的雾化体积和功率需求方程, 从而为制备涂层稳定和均匀的血管支架提供了理论支撑。
超声喷雾粒径检测实验表明大部分情况下MediCoat 超声喷涂系统工作在理想雾化和射流雾化区间, 当超声波频率达到 120 kHz 时, 喷雾粒径尺寸可控制在 10 ?m 左右; 同等流量下, 降低喷涂液黏度和表面张力, 增大雾化功率有利于减小喷雾粒径, 以适应血管支架微小尺寸空间均匀喷涂要求; 同时Rosin-Rammler 分布函数拟合超声波喷雾均匀度指数普遍大于 7, 喷雾质量明显优于常规喷雾, 且增加流量有利于提高喷雾粒径分布均匀性. 超声波喷雾在飞行过程中存在少部分粒子碰撞结合的现象, 但总体稳定, 表明血管支架在喷雾场中的位移调节较为自由和灵活.最后, 在前面研究基础上制定并优化喷涂工艺参数, 并依据优化后的工艺参数进行血管支架涂层制备试验, 试验结果消除了血管支架表面涂层常见的各种缺陷, 获得了性能优良的血管支架表面涂层。
参考文献略
本站文章未经允许不得转载;如欲转载请注明出处,北京桑尧科技开发有限公司网址:http://www.sunspraying.com/
|