激光光梯度力消除雾霾的实验研究

2018-10-15张德峰陈晓伟王海容杜军姜久龙吴九汇

西安交通大学学报 2018年10期

张德峰,陈晓伟,王海容,杜军,姜久龙,吴九汇

(1.西安交通大学机械结构强度与振动国家重点实验室,710049,西安;2.西安交通大学机械制造系统工程国家重点实验室,710054,西安;3.空军工程大学航空航天工程学院,710038,西安)

雾霾天气是一种大气污染状态。雾霾是对大气中各种悬浮颗粒物含量超标的笼统表述,其中霾是由灰尘、硫酸、硝酸、碳氢化合物等非水成物组成的气溶胶系统;雾是近地面的微小水滴或冰晶凝结形成的气溶胶系统。雾霾使大气浑浊、能见度下降并对人体呼吸系统和心血管系统等造成严重危害。此外,雾霾污染已经阻碍了中国吸引外商投资,并对国际旅游业造成冲击,影响远远超过经济利益的损失[1]。鉴于雾霾治理的紧迫性,国务院总理李克强在2017年两会期间表示,如果有科研团队能够把雾霾的形成机理和危害性真正研究透,提出更有效的应对良策,愿意拿出总理预备费给予重奖。

当前,政府在治理雾霾等环境污染方面下定决心,但取得的成效与人们的期待相差甚远。针对雾霾这个难题,一方面,国内外许多专家做了大量研究,Chameides等研究了气溶胶光学厚度和区域霾污染对水稻和小麦产量的影响[2],Shakya等测定了尼泊尔城市居民区总悬浮颗粒物、PM10、PM2.5、黑碳和颗粒物的化学成分[3],姚青和蔡子颖等[4]研究了天津冬季雾霾天气下颗粒物质量浓度分布与光学特性,表明雾霾颗粒物的散射作用是造成大气能见度降低的主要原因;另一方面,政府也投入了大量资金及人力来研究雾霾治理。

但是,目前这些治理以及防范雾霾的方法及措施,对减轻雾霾污染的效果微乎其微,并且对雾霾粒子能够长时间悬浮的内在本质一直没有深入研究分析。因此,针对雾霾治理这一大难题,我们课题组基于雾霾粒子悬浮的内在原因,从力学角度对其受力进行了深入分析[5],率先提出了基于激光光梯度力破坏气溶胶雾霾粒子力平衡体系进而消除雾霾的新机理。

本文基于激光光梯度力消除雾霾的新机理,设计了一个小型雾霾模拟实验装置,构造了现实中的雾霾环境,即常温常压下,空气、雾和霾共存的满足湿度、颗粒物浓度等各项雾霾指标的3项体系[6]。雾霾模拟实验舱能够产生特定温湿度、雾霾浓度和粒径分布与成分的雾霾环境,同时能监测雾霾的参数,并且首次使用激光光梯度力破坏雾霾颗粒的平衡体系,从而使雾霾颗粒沉降的方法进行实验研究。

1 激光消除雾霾机理

雾霾颗粒在大气流中所受力为:曳引阻力、旋转升力、Saffman升力、Basset力、范德瓦尔斯斥力、静电力、虚假质量力、重力、浮力等。忽略其他次要力,将雾霾粒子等效为球体,大气流中雾霾粒子的受力以及与激光作用后的沉降过程如图1所示。粒子与粒子之间在范德瓦尔斯力的作用下,互相排斥形成了稳定的网状力平衡状态,悬浮微粒难以扩散而被阻滞在低空和近地面,因此消除雾霾的本质就是破坏雾霾粒子的力平衡体系。微粒在激光范围内受到一均匀时变场的作用,在场作用力下,微粒有向光场较强区域移动的倾向,这种场作用力就是指向光轴的梯度力。当激光产生的光梯度力足以打破雾霾粒子的力平衡体系时,粒子便会脱离原有的平衡状态,在光轴处聚合坠落。

图1 激光作用下雾霾粒子的受力情况

1.1 光梯度力对雾霾粒子捕获能力分析

1970年,Ashkin等人在美国贝尔实验室发现了微米级聚苯乙烯粒子在光梯度力的作用下被拉入光轴[7],从此开创了光镊领域的研究[8-9],光镊不仅可以俘获小至5 nm[10-11]的物体,还可以施加超过100 pN[12-13]的力。假设雾霾粒子的折射率n1大于周围介质的折射率n2,当光子与微粒相互作用时,光子的部分动量将转移到雾霾粒子上,这种动量转移[14-15]正是产生光梯度力的机理[16]。对光学捕获的恰当物理描述:依赖于被捕获粒子的尺寸d相对于捕获光的波长λ,一般来讲,当d≫λ时应用射线光学理论,当d≪λ时应用瑞丽散射理论。由于雾霾粒子的尺寸与激光波长相当,因此较为正确的描述应介于这两种理论之间。为了能够进一步对雾霾粒子的光学捕获原理进行一个定性的理解,可以应用射线光学理论,更加直观地解释光梯度力。

(a)高斯激光束横向梯度力 (b)高斯激光束纵向梯度力图2 射线理论中光学捕获粒子的梯度力

应用射线光学理论,可从两种介质折射率不同的角度来理解捕获力,图2定性地说明了横向梯度力和纵向梯度力的产生原因。横向梯度力可解释如下:当光线I1和I2的强度不同时,光线传递的动量量级也不同,将传递到雾霾粒子上的动量-dp1和dp2沿x轴分解,从而在雾霾粒子中产生一个反作用力,方向指向光场场强最大区域[17],即高斯光束的光轴,这个反作用力的x轴分力趋向于遏止粒子偏离光轴。轴向梯度力的产生同样是由于光折射后的动量传递,在雾霾粒子中产生一个指向焦点的回复力。在z轴方向上,粒子被牢牢地束缚在焦点附近,激光即可实现对雾霾粒子的三维捕获。

1.2 激光消除雾霾机理分析

将雾霾颗粒等效为球体,大气流中雾霾颗粒主要受到的力如下。

(1)气流曳引阻力为

(1)

式中:ρg为流场中气体的密度;Ag为颗粒物在与空气主流流速垂直面上的投影;CD为黏性阻力系数,取CD=24/Re;ug、up分别为气流和雾霾粒子的速度。

(2)颗粒自转形成的旋转升力为

(2)

将雾霾颗粒等效为球体时,需引入修正系数K

(3)

式中:dp为雾霾粒子直径;ρp为雾霾粒子密度;ω为雾霾粒子旋转速度。

(3)流场中的速度梯度导致颗粒在运动时上下部分的流速不同,因此颗粒将受到一个沿垂直方向力的作用,当颗粒上表面速度大于下表面速度时,形成一个向上的Saffman升力Fs,表达式如下

(4)

式中:μg为气体动力黏度,Pa·s。

(4)范德瓦尔斯斥力为

(5)

式中:w为黏附功;h为雾霾颗粒的间隙;z0为范德瓦尔斯力平衡间距,取z0=0.4 nm。

光梯度力[18]为

Fgrad(r,t)=[p(r,t)]E(r,t)

(6)

Fgrad(r,t)=[αE(r,t)]E(r,t)

(7)

利用(E(r,t)E2-E(E)和麦克斯韦方程E=0,式(7)可改写成

(8)

所以,激光场的光梯度力为

(9)

对大气流中雾霾颗粒主要受到的力进行计算,计算结果如图3所示,在不同粒子浓度下,功率为2.6 mW的氦氖激光器对粒径分布在0.3～2.5 μm的雾霾粒子产生的光梯度力远大于其处于悬浮状态时所受的各种力。光梯度力捕获雾霾粒子情况如下。

(1)激光光梯度力对雾霾粒子的捕获能力远大于大气对雾霾粒子的悬浮能力,雾霾粒子间原有的力平衡体系被打破,从而加速了雾霾粒子的沉降。

(2)光梯度力随着粒子半径递增,并随着雾霾浓度的增大而增大,在激光逐渐衰减后而递减。

(a)雾霾粒子浓度为200 μg/m3

(b)雾霾粒子浓度为500 μg/m3图3 雾霾粒子所受的各种力

2 激光消除雾霾实验设计及验证

实验最直接的方法是在雾霾天气时选取一片空旷的场地并铺满白纸,向大气中的雾霾辐射激光,在激光光梯度力的作用下,雾霾粒子斥力平衡体系被打破从而相互聚合、掉落,对比在有、无激光照射区域下白纸上颗粒的数目,分析激光是否能够消除雾霾。但是,考虑到风力的影响,雾霾粒子在坠落过程中会被吹离激光辐射区域,导致颗粒数目不准确。此外,只能选取雾霾天气的自然因素限制了实验进行,针对单一的雾霾天气状况也难以分析激光光梯度力消除雾霾的主要技术参数。鉴于上述诸多因素的影响,实验难以在室外进行,因此本文实验的关键在于搭建一个雾霾粒子能够长时间悬浮的雾霾模拟实验舱,并配合颗粒物发生装置、温湿度计、颗粒稀释器、激光粒子计数器等实验设备,实时监测实验区域的颗粒物浓度、粒径分布、温湿度等参数指标。

2.1 实验布局及实验原理

根据国家标准GB/T 18801—2015《空气净化器》,搭建了进行烟雾粒子悬浮实验的雾霾模拟实验舱,实验舱内尺寸为1.4 m×1.4 m×1.5 m,容积为3 m3,实验舱示意图、实物图如图4所示。实验舱主要分为控制区、进样区和观察采样区。舱内设有高效空气过滤器,用于净化实验舱内空气,在输入烟雾之前,开启循环风道的风阀,使舱内气体循环通过高效空气过滤器,将颗粒物的背景浓度降低到标准值以下。搅拌风扇位于舱内顶部,主要用于使输入的烟雾进入搅拌风扇形成的空气涡流中,保证舱内烟雾粒子分布均匀。在实验舱的底部放置硅胶加热垫,保持温度恒定在25 ℃。进样区与颗粒物发生装置相连接,选用标准颗粒物:香烟烟雾(红塔山经典150),焦油量为8 mg,将一支标准香烟放入香烟燃烧器内,香烟烟雾出口链接一根穿过实验舱壁的输送管路,空气压缩机持续提供支持燃烧的空气和动力[19],烟雾经过搅拌风扇进入空气涡流中。观察采样区与颗粒稀释器和激光粒子计数器相连接,经过颗粒稀释器稀释采集的气体,稀释采样流量为20 mL,采用的稀释比为960∶1。将稀释后的气体通入粒子计数器,记录颗粒数目,粒子计数器的示数打印间隔设为1 min,计数器的采样流量为2.83 L/min,与颗粒稀释器相匹配。

(a)实物图(b)示意图图4 雾霾模拟实验舱

激光器布置在实验舱进样区对面的工作台面上,由激光固定底座夹持,激光通过进样口从实验舱外部向实验舱内部辐射。为保证实时测量烟雾浓度,调整激光固定底座支架的角度,使激光入口、烟雾采集口与激光束保持在同一水平高度上。入射的激光光束径直辐射烟雾采集口,若光梯度力与悬浮颗粒相互作用而起到沉降粒子的效果,则一定是采集口处颗粒浓度最先发生变化,避免了舱内粒子因再次混合造成的时间延迟。使用的激光器为氦氖激光器(HN-250),功率为2.6 mW,波长为632.8 nm。实验中需要凸透镜对激光扩束,扩大激光光束的照射范围,用以提高粒子沉降效率,凸透镜固定在舱壁的激光入口处,其中凸透镜焦距为5 cm,采集口到凸透镜的距离为0.5 m,此时光斑恰好覆盖整个采集口区域。

设计的实验方案为在烟雾粒子悬浮的中途开启激光器,观测各个阶段粒子衰减速率。在实验舱内进行多次烟雾悬浮实验时,由于难以精确控制进烟过程的各个因素,例如香烟的燃烧状态不同导致粒径分布的差异,多次实验之间的对照分析不够合理。因此,实验中采用的方案为在单次实验的中途开启激光器,同时观察颗粒浓度的变化,通过对比每个阶段曲线曲率的大小,分析激光与烟雾是否产生相互作用,从而分析光梯度力是否能够打破烟雾粒子的力平衡体系,达到沉降颗粒的目的。

2.2 雾霾模拟方法

雾霾模拟技术一直是大气环境模拟的一个难点,本文采用燃烧法来模拟雾霾颗粒,选择香烟作为燃烧物来模拟雾霾环境,模拟的参数主要有:雾霾浓度、雾霾成分、粒径分布、相对湿度和温度,实验舱进烟时各个阶段如图5所示。

(1)雾霾浓度。依照国家标准GB 3095—2012《环境空气质量标准》HJ618标准,取PM10的二级日平均标准,按照一般空气、轻度雾霾(150 μg/m3)、中度雾霾(300 μg/m3)、重度雾霾(450 μg/m3)划分。通过控制香烟燃烧量来控制雾霾浓度,根据香烟颗粒不易聚合、稳定悬浮浓度高等特点,选用的测量仪器为激光粒子计数器。

(2)雾霾成分。雾霾化学成分主要有硫酸盐、硝酸盐、铵盐,有机碳氢化合物粒子以及矿物颗粒等[20]。香烟经燃烧后可生成细颗粒物以及烟焦油,其中烟焦油是有机质在缺氧条件下不完全燃烧的产物,是众多烃类及烃的氧化物、硫化物及氮化物的极其复杂的混合物。经分析得出,香烟烟雾与雾霾在化学成分上具有一致性。

(3)粒径分布。雾霾天气下PM2.5和PM10质量比大于74%,说明细污染物是造成冬季污染的主要因素[21-22]。香烟烟雾中PM2.5和PM10质量比约为84.44%,且0～1.0 μm的颗粒物约占PM10总质量的68.89%[23],与雾霾粒径分布相似。

(4)相对湿度和温度。雾霾环境的相对湿度在80%～90%,取85%为试验湿度,雾霾环境的温度取室外空气温度25 ℃。

(a)进烟前(b)进烟1 min

(c)进烟15 min(d)进烟30 min图5 模拟雾霾的悬浮情况

2.3 实验结果及分析

选用以散射法为测量原理的激光尘埃粒子计数器记录烟雾粒子的浓度。光学传感器的探测激光经尘埃粒子散射后被光敏元件接收并产生脉冲信号,微粒散射光的强度随微粒的表面积增加而增大,将测定散射光的脉冲信号放大,进行数字化处理,并与标准粒子散射信号进行比较。其中,电脉冲数量对应于微粒的个数,电脉冲的幅度对应于微粒的大小,即可推知微粒的个数和大小。采用在雾霾自然衰减过程中开启激光器的方式开始颗粒物的衰减实验,分别观测烟雾粒子在自然状态下的衰减速率以及在激光辐射状态下的衰减速率,实验步骤如下:

(1)将采样点位置布置好,避开进出风口,离玻璃壁距离为0.5 m,相对实验室地面高度为0.5 m,并与实验舱外采样器相连接;

(2)开启高效空气过滤器,净化实验舱内空气,使颗粒物粒径在0.3 μm以上的粒子背景浓度小于1 000 L-1,同时启动温度控制装置,使室内温度达到规定状态;

(3)待颗粒物背景浓度降低到合适水平,记录颗粒物背景浓度,关闭高效空气过滤器,启动搅拌风扇,烟雾输送完毕后,关闭烟雾输送管阀门,搅拌风扇工作10 min,使颗粒污染物混合均匀后关闭搅拌风扇;

(4)待搅拌风扇停止转动后,用颗粒计数器测定颗粒物的初始浓度C0;

(5)实验舱内的初始浓度(自然衰减的第1个取样点)测定后,每1 min测定1次,并记录颗粒物的浓度;

(6)在规定时间内开启激光器,观测烟雾粒子浓度在激光梯度力的影响下的衰减速率。

本次实验在观测烟雾粒子浓度时,同时对直径为0.3、0.5、1.0 μm的颗粒浓度进行计数,用以分析光梯度力对不同大小微粒的作用情况。在16～39 min时间段内开启激光器,实验结果如图6所示。

(a)直径为0.3 μm

(b)直径为0.5 μm

(c)直径为1.0 μm图6 雾霾浓度随时间的变化情况

由图6可见,开启激光器,即在16 min处,直径分别为0.3、0.5、1.0 μm的粒子浓度曲线均有非常明显的转折,粒子浓度剧烈下降,并且在激光开启过程中持续下降。3种不同大小的粒子衰减斜率均大于0～16 min的自然衰减曲线斜率,说明激光对粒子沉降起到了显著作用,其中对1.0 μm直径的粒子沉降效果最好。

激光器开启30 min后,烟雾粒子在有激光作用的情况下衰减不明显,有以下两条原因:第一,在30～39 min时间段内,若此时舱内底部的硅胶加热垫处于加热状态,受到激光作用而沉降的烟雾粒子,在自然对流的影响下被舱内流体向上托起,导致采样数偏大,因此曲线没有明显的下降趋势;第二,衰减速率与粒子浓度有关,舱内粒子浓度较低时,受到激光作用的粒子数目减少,因此曲线趋向于平缓。另外,从曲线中可以发现粒子浓度上升的现象,原因可能在于加热垫此时正在加热,从而影响了粒子浓度。因此,在后续的实验中,实验舱在保温可以采用设计外舱的形式,也可以采用设计保温层的形式,避免采样区受自然对流的影响。