王舰航,尹招琴,凃程旭,包福兵

(中国计量大学 计量测试工程学院,杭州 310018)

0 引 言

空化是指液体内局部压强低于临界值时气泡的形成、发展和溃灭的动态过程[1]。空化现象最早发现于驱逐舰的螺旋桨上,空化气泡溃灭产生的冲击波和微射流作用在螺旋桨叶片上,对叶片产生冲击作用出现空蚀,损坏螺旋桨叶片表面,进而出现振动和噪声,影响推进效果[1]。近年来由于超声清洗技术[2]、空化污水处理技术[3]、激光空化强化技术[4]以及空化技术在军事上的广泛应用[5],国际国内兴起了一股空化应用的研究热潮[6]。

空化技术在生物医学上也有广泛用途,随着基因工程和分子生物学的发展,基因治疗[7]被认为是最有前景的癌症治疗方法之一。基因治疗的关键是找到一种安全高效的基因转染方法,而超声靶向微泡破坏技术(UTMD)[8]是其中最有发展潜力的一种。超声靶向微泡破坏技术利用微气泡在超声波场中的动力学响应,诱发细胞膜声孔效应来增加膜的通透性,进而有效促进大分子物质的胞内传输。Taniyama等[9]通过电镜观测到超声作用后,细胞膜上产生可恢复的孔洞,验证了超声靶向微破坏技术的可行性及其临床应用价值。但到目前为止超声作用下气泡溃灭如何增加细胞膜通透性的机理尚未有定论。

以往的近壁气泡空化溃灭研究往往关注气泡在刚性壁面附近溃灭时产生的射流及其对壁面的冲击作用,分析气泡尺寸和近壁距离等参数的定量影响[10]。近年来的研究发现,气泡在弹性壁面附近溃灭时存在与刚性壁面不一样的行为。Chen等[11]研究了超声作用下血管内的微气泡溃灭过程,发现溃灭射流并不是射向血管壁,而是射向血管中心,这与刚性壁面的结论正好相反,他们认为是壁面材料性质造成了射流反向。Brujan等[12]通过实验研究了激光空化气泡在弹性壁面附近的溃灭现象,发现弹性壁面的弹性模量较小时,气泡溃灭射流远离壁面;在弹性模量较大时,气泡溃灭射流朝向壁面。Ma等[13]发现在空化泡生长和溃灭过程中弹性壁面的被动变形使得附近微气泡的溃灭行为有别于刚性壁面。然而细胞与弹性壁面在弹性模量等性质上存在较大区别,而且由于细胞与气泡尺寸相当,当附近空化气泡溃灭时具有更丰富的动力学响应行为。

超吸水聚合物(super absorbent polymer,SAP)具有与癌细胞较为接近的密度和弹性模量[14],SAP弹性小球在空化气泡溃灭流场中具有较好的变形行为[14],能够较好地近似UTMD 技术中细胞在附近微气泡溃灭时的动态响应过程。因此,本文拟采用激光诱导方法来生成单个可控空化气泡,研究SAP 弹性小球在空化气泡膨胀与溃灭时的形变过程。

1 实验装置与方法

图1 实验系统示意图Fig.1 Schematic diagram of experimental system

实验中,SAP 弹性小球通过细铜丝悬挂于水箱的去离子水中,水箱放置于精密位移平台上,通过调整位移平台来控制激光空化气泡中心位于小球正下方,并且可以调整气泡中心与小球的距离。实验中的SAP弹性小球为高吸收性聚合物材料,可以吸收并保留大量液体,据文献中的实验测量数据,其泊松比v＝0.49[15],杨氏模量E ＝21 k Pa[14],密度ρ＝1117 kg/m3。

实验中,固体激光器产生的脉冲激光经过扩束-聚焦镜片组后被聚焦至水箱的去离子水中,当聚焦点的激光功率密度超过去离子水的击穿阈值时,去离子水被电离击穿,产生等离子体并迅速膨胀形成空泡[4],类似微爆炸。激光聚焦点即为空化气泡中心,气泡尺寸可通过固体激光器能量耗散器调节。实验中气泡中心到水箱壁面的距离远大于气泡尺寸,因此认为壁面的影响可忽略不计[4]。实验中定义球泡间的无量纲距离为γ[16]:

式中,Rmax为空化气泡膨胀到最大时的半径,r 为SAP弹性小球的半径,L 为气泡中心到弹性小球表面的最近距离,实验参数如图2所示。

图2 实验参数示意图Fig.2 Schematic diagram of experiment parameters

实验过程中,通过高速相机捕捉记录下完整的空化气泡膨胀与溃灭过程,使用无频闪常亮光源作补光。为保证图像清晰度,实验中高速相机的采样频率为50 000帧/秒,对应两张照片之间的时间间隔为20 μs,此时图像的分辨率为320 pixel×232 pixel。后期通过图像处理软件ImageJ来获得空化气泡的最大半径以及弹性小球变形参数。在实验中,由于高速相机帧率的限制,相机拍摄到的气泡最大半径很大概率不是空化气泡的实际最大半径。而且由于激光空化时的偶然性,在激光能量一致时也不能完全保证每次空化气泡大小的一致性。

2 空化气泡溃灭演化过程

图3给出了SAP弹性小球附近激光空化气泡的溃灭演化以及弹性小球的变形过程。实验中采用的弹性小球半径r＝1.74 mm,激光空化气泡中心到小球中心距离L＝0.4 mm,实验中生成的空化气泡的最大半径Rmax＝1.59 mm,所以根据公式(1)计算可知,球泡无量纲距离γ＝1.35。实验中,根据延时触发器设置,定义相机拍摄第一张照片为t＝0 s。

图3 气泡在SAP弹性小球附近的膨胀与溃灭过程Fig.3 Process of bubble expansion and collapse near a SAP elastic ball

从图3可以看出,SAP 弹性小球附近的激光空化气泡存在多次膨胀与溃灭过程。每次可以分为两个阶段:膨胀阶段和溃灭阶段。以第一次膨胀与溃灭过程为例,图3(a～h)是气泡的一次膨胀阶段,持续140μs左右。在气泡膨胀过程中,随着气泡的发展变大,弹性小球开始被挤压变形,在气泡达到最大体积时(图3(h))弹性小球挤压变形也最大。另一方面,由于弹性小球的存在,气泡膨胀时形状也受到影响,与自由溃灭时的球形不一样,气泡呈现非球形。

当空化气泡达到最大体积之后,在内部气液相变和外部流体压力的共同作用下,气泡体积迅速变小,发生溃灭。图3(h～n)是气泡的一次溃灭阶段,持续时间120μs左右。气泡溃灭最后时刻,由于周围液体的瞬间冲击,会在气泡中心出现往四周传播的冲击波,但是在本实验中,由于相机时间分辨率和空间分辨率的限制,没有捕捉到冲击波现象。

进一步分析发现：在慢代谢组中，服用氯吡格雷的患者心血管缺血事件发生率显著高于替格瑞洛（χ2=5.506，P=0.019）；慢代谢组服用氯吡格雷的患者心血管缺血事件发生率显著高于快代谢组和中代谢组（χ2=4.192，P=0.041）。

气泡一次膨胀与溃灭总持续时间在260μs左右,当气泡一次溃灭结束后,气泡体积又重新膨胀,进入二次膨胀(图3(n～q))和二次溃灭(图3(q～u))阶段,总持续时间在140μs左右,二次膨胀最大体积小于一次膨胀时的最大气泡体积。当二次溃灭结束后,还存在三次膨胀与溃灭(图3(u～x)),总持续时间在60μs左右,而且气泡体积更小,溃灭强度更弱。实验中,由于相机时间分辨率和空间分辨率限制,气泡溃灭持续时间和气泡最大尺寸都存在一定的误差。当高速相机时空分辨率进一步提高时应该还会出现四次以及更高次的膨胀与溃灭阶段。

从图3中还可以发现,与刚性壁面[17]和刚性小球[16]不一样的是,气泡在弹性小球附近溃灭时,并没有出现射向弹性小球的微射流。在刚性壁面和刚性小球附近,当气泡溃灭阶段周围液体往气泡中心汇聚时,由于壁面的存在,气泡靠近壁面侧的液体补充不及时,而远离壁面侧的液体往气泡中心运动时会穿过近壁面侧的气泡界面,形成了指向刚性壁面的微射流[16-17]。细胞也具有较大的可变形性,在超声靶向微泡破坏技术中,微气泡在细胞附近溃灭时,并不一定会出现射向细胞的微射流,溃灭微射流可能并不是造成气泡膜通透性增加的原因。

3 弹性小球变形分析及模拟

3.1 弹性小球变形分析

在空化气泡的膨胀与溃灭过程中,附近弹性小球出现明显变形。从图3中可以发现,在一次膨胀过程中,随着气泡迅速变大,弹性小球被挤压变形,在气泡达到最大体积时(图3(h))弹性小球变形也最大。气泡一次溃灭时,周围液体迅速向气泡中心汇聚,速度较大,流体中存在较大的剪切,造成弹性小球被拉扯变形。在一次溃灭气泡体积最小时刻,弹性小球下部被拉扯出一个尖端,如图3(n)所示,此时,弹性小球表面受到最大应力。

在气泡一次膨胀与溃灭过程的弹性小球变形中,溃灭阶段时间更短,而且变形存在奇点,因此有理由相信,在超声靶向微泡破坏技术中,超声作用下10微米量级的气泡在相同尺寸细胞附近溃灭时,细胞也会出现类似变形,从而增加了细胞膜的通透性,进而实现了达到基因转染的目的。

从图中还可以看出,二次溃灭和三次溃灭对气泡的挤压和拉扯作用较弱,在这两个阶段,SAP弹性小球在其自身弹性的作用,逐渐恢复其原来的球状。

气泡的膨胀与溃灭过程伴随着周围液体的运动,气泡溃灭时由于周围液体往气泡中心急速运动,会伴随着压力波的生成[1,4]。在各向同性的均质物体中,压力波沿纵向直线传播,压力波的传播速度up由下式给出[18]:

式中K 是体积模量,G 是剪切模量,ρ 是波传播介质密度。均质各向同性线弹性材料的弹性特性可由其中的任何两个模量唯一确定,给定任意两个可以计算其他弹性模量[18]:

式中,E 为弹性小球杨氏模量,v 为弹性小球泊松比。

根据SAP弹性小球的泊松比[15]和杨氏模量[14],可以计算得知压力波在直径为3 mm 的弹性小球内部传播的时间尺度为9.3μs,而弹性小球受到气泡一次膨胀与溃灭时冲击作用时间约为260μs,压力波在弹性小球内的传播速度远高于气泡溃灭引起弹性小球形变的速度,因此排除压力波对本实验的影响,这也解释了弹性小球形变主要在靠近气泡的下半球部分,而上半球部分几乎无形变的原因。

3.2 弹性小球变形模拟

空化气泡膨胀与溃灭时和弹性小球的相互作用过程涉及到复杂的能量转换,为了分析弹性小球在气泡膨胀阶段的受力,模拟了弹性小球的受力变形过程。采用ANSYS Workbench里多体动力学分析中的瞬态动力学模块对弹性小球形变过程进行模拟,可以得到不同半径SAP弹性小球在不同压力荷载作用下的形变量。从图3可以发现,弹性小球上半部分在空化气泡一次膨胀溃灭过程中几乎没有变化,因此采用如图4所示的半球模型,通过在弹性半球上方放置平板并施加压力,观察弹性小球的形状变化情况。弹性球物性参数按照Goh等[14]的实验测量数据设置,弹性半球的半径根据该实验设置,采用5种小球(r＝1.45、1.58、1.74、2.13、2.68 mm)。

图4 弹性半球模型Fig.4 Elastic hemisphere model

弹性半球的形变量Δh,定义为弹性球达到最大形变时高度方向的形变量。图5给出了5种SAP弹性半球在不同外加载荷作用下的变形情况。随着外加压力的增加,弹性小球的形变量逐渐增大,但是变化趋势逐渐变慢。作用力比较小时,弹性小球半径越小,越容易变形。但是当形变量达到一定程度之后,相同作用力下半径越小的弹性小球形变量也越小。基于图5中气泡膨胀阶段弹性小球的最大形变量与作用的压力荷载之间的F－Δh 曲线,根据实验中测得的SAP弹性小球形变量,就可以获得作用在弹性小球上的作用力。图6给出了5种SAP弹性小球在不同球泡距离γ 时空化气泡膨胀阶段受到的最大挤压作用力。

图5 SAP弹性半球在不同压力荷载下的变形情况Fig.5 Deformation of SAP elastic hemisphere under different pressure loads

因为实验中激光空化气泡略小于弹性小球尺寸,按照方程(1)中无量纲距离γ 的定义,气泡中心离弹性小球越近,无量纲距离γ 越接近1;而气泡中心离弹性小球越远,无量纲距离γ 也越大。从图中可以看出,对于给定尺寸的弹性小球,气泡中心离弹性小球越近,气泡膨胀阶段对弹性小球的挤压作用越大。当γ＞2时,空化气泡一次膨胀时对气泡的挤压作用力小于0.01 N;而当γ＜1.5时,气泡对弹性小球的挤压作用力可达到0.05 N,气泡中心越靠近弹性小球,弹性小球受到的挤压作用越大。Zhang等[16]在研究刚性小球附近的气泡溃灭时,他们发现当无量纲距离γ＞2.5时,附近气泡溃灭与壁面几乎相互无影响,而气泡在弹性小球附近溃灭时,此时气泡溃灭还会影响弹性小球变形。

图6 弹性小球受到的冲击压力随无量纲距离γ 的变化Fig.6 Variation of impact pressure on the elastic sphere with the dimensionless distanceγ

4 结 论

本文研究了激光空化气泡膨胀和溃灭时对附近SAP弹性小球形变的影响,获得以下结论:

1)通过高速相机获取的空化气泡与弹性小球的动态演化过程发现空化气泡存在多次膨胀与溃灭过程,一次膨胀与溃灭最剧烈,持续时间最长,二次和三次溃灭逐渐减弱。

2)在气泡作用下,SAP弹性小球出现明显变形,在气泡膨胀阶段,弹性小球下半部分受迫挤压,在空化气泡达到最大体积时弹性小球变形最大。而当气泡一次溃灭时,由于受到流体强剪切的作用,弹性小球被剧烈拉扯,出现明显变形。气泡的二次和三次膨胀与溃灭过程对弹性小球变形影响较小。

3)通过弹性小球形变过程的模拟,还获得了不同球泡无量纲距离γ 时气泡溃灭时弹性小球受到的挤压作用。随着球泡无量纲距离变小,膨胀气泡对弹性小球的挤压作用增大。

4)与刚性壁面和刚性小球不同的是,气泡在弹性小球附近溃灭时没有出现射流现象。

5)根据SAP 弹性小球在附近气泡膨胀与溃灭时的形状变化过程,我们认为超声靶向微泡破坏技术中,是空化气泡溃灭时强流体剪切引起的挤压和拉扯作用,而不是溃灭射流增加了细胞膜的通透性。由于细胞结构的复杂性,SAP 弹性小球在附近气泡膨胀与溃灭时的动力学响应并不能完全揭示细胞膜通透性增加机理,细胞在超声靶向微泡破坏技术中的响应还有待进一步研究。