肖茂栋¹, 张卓^1,2, 阮兵¹, 李兴华¹
1. 北京航空航天大学空间与环境学院, 北京 100191;
2. 中钢集团武汉安全环保研究院, 武汉 430081
收稿日期: 2019-10-21; 修回日期: 2020-06-10; 录用日期: 2020-06-10
基金项目: 国家重点研发计划(No.2017YFC0211404)；国家自然科学基金(No.41775121)
作者简介: 肖茂栋(1995-), 男, E-mail:xmdlarry@163.com
通讯作者（责任作者）: 李兴华(1972—), 男, 现为北京航空航天大学空间与环境学院副教授, 研究领域为大气污染物的形成、排放、影响与控制.E-mail: lixinghua@buaa.edu.cn

摘要：为研究挥发性组分对气溶胶特性的影响，设计加工了一种热扩散管(Thermal Denuder，以下简称TD)，由解析管和吸附管组成，调节温度范围为50~350℃；与文献中已有设计相比，本设计改进了测量TD管内部温度的位置，使得测温、控温更加准确；采用两段加热的方式使得温度分布更为稳定，延长有效停留时间.TD管性能测试表明：TD管内温度分布较均匀，可有效去除挥发性组分；当颗粒物的粒径大于60 nm时，颗粒物在不同温度的TD管内的传输效率变化很小，可近似看作是一个常数.
关键词：热扩散管性能测试温度分布颗粒物传输效率挥发性能
Design and performance test of a thermal denuder
XIAO Maodong¹, ZHANG Zhuo^1,2, RUAN Bing¹, LI Xinghua¹
1. School of Space & Environment, Beihang University, Beijing 100191;
2. Wuhan Institute of Safety and Environmental Protection, Sinosteel Group, Wuhan 430081
Received 21 October 2019; received in revised from 10 June 2020; accepted 10 June 2020
Abstract: In order to study the effect of volatile components on aerosol characteristics,a thermal denuder (TD) was designed,which is composed of a denuder tube and an adsorption tube. The operation temperature of this TD ranges from 50℃ to 350℃. Compared with the existing designs in the literatures,the internal temperature measurement position in our design was optimized to improve the accuracy of temperature measurement and control. Furthermore,the two-stage heating method was adopted to improve the stability of internal temperature distribution,and elongate the effective residence time of particles within the analytical tube. The middle connecting flange of the adsorption tube is insulated by heating and glass wool respectively,which can effectively reduce the temperature drop at the outlet of the denuder tube. The TD performance test showed that temperature inside the TD tube is well-distributed and could effectively reduce the volatile components. The transport efficiency of particles larger than 60 nm inside the TD tube showed litte change at differential temperatures,which can be approximately regarded as a constant.
Keywords: thermal denuderperformance testtemperature distributionparticle transport efficiencyvolatilization performance
1 引言(Introduction)挥发性是大气气溶胶重要的物理性能, 挥发性物质在大气气溶胶中占有相当一部分比例, 可同时存在于气相与颗粒相中, 影响着颗粒相-气相分配系数, 从而改变着大气气溶胶中不同组分的浓度、大气寿命等(Huffman et al., 2009).
对气溶胶挥发性进行分析和测量时, 常采用热分离(Thermal fractionation)的方法；采用该方法, 不同的化学物质会在与其蒸气压、沸点和蒸发焓相关的特征温度下迅速蒸发.随着气溶胶被加热, 挥发性较强的物质将会被优先蒸发, 而挥发性较弱的物质仍然会被保留为固相.一般情况下, 热分离方法采用加热管实现, 该加热管一般被称为热扩散管(Thermal Denuder, 简称TD管).大多数TD管采用石英管或不锈钢管, 并采用加热带包裹或置于烤炉中的形式进行升温加热.气溶胶组分在TD管不同设定温度下分离出来, 大致可分为以下几个阶段：①T=30~125 ℃, 高挥发性有机碳(Highly volatile organic carbon, HVOC)和硫酸；②T=125~175 ℃, 硫酸铵和硫酸氢铵；③T=175~300 ℃, 低挥发性有机碳(Lowly Volatile Organic Carbon, LVOC)；④T >300 ℃, 难熔气溶胶, 如元素碳(Elemental Carbon, EC)、矿物粉尘和氯化钠等(Burtscher et al., 2001).
目前, TD管已成为研究大气气溶胶挥发性的常用工具, 常与其他仪器相结合使用, 如高分辨率飞行时间气溶胶质谱仪(High Resolution Time of Flight Aerosol Mass Spectrometer, AMS)、扫描电迁移率粒径谱仪(SMPS)等(Wehner et al., 2004；Huffman et al., 2009；黄聪妮等, 2015；李园园等, 2015).此外, TD管还被广泛应用于多种形式的挥发性串联差分电迁移率分析仪(Volatility Tandem Differential Mobility Analyzer, VTDMA)中, 并应用于研究外场大气气溶胶的挥发性, 以及实验室内发生的有机气溶胶的挥发性和组分(Villani et al., 2008).
近年来, 研究人员不断进行着TD管的开发与改进(表 1).例如, Burtscher等(2001)设计了一种TD管(以下简称为Burtscher TD), 采用不锈钢制成解吸管(Desorber tube), 并在解吸管的薄壁上设计了冷却管道, 加快解吸管的降温, 可快速调节解吸管温度；解吸管后紧接有一根吸附管(Adsorber tube), 并采用活性炭填充吸附管内部, 吸附热分离出来的气体, 避免气体二次冷凝到气溶胶表面.Wehner等(2002)在Burtscher等的基础上, 对TD管进行了相应的改进(以下简称为Wehner TD), 采用填充细沙的方式提高加热的均匀性, 并在加热段的两端侧板添加了圆形陶瓷片, 减少了加热段的热扩散, 增强了保温效果.Huffman等(2008)在Wehner等和Burtscher等的设计基础上, 在外壳上增加了两个风扇, 通过吹扫空气以便快速进行温度调节；并将加热段分成了三部分加热区域, 每一区域的温度都由比例-积分-微分控制器(Proportion Integral Differential, 以下简称PID)控制；整体提高了温度的分布均匀性, 并且增加了颗粒物挥发的停留时间(以下简称为Huffman TD).
表 1(Table 1)

表 1 自制和商业化TD管各种参数 Table 1 Various parameters of self-made and commercial TD

表 1 自制和商业化TD管各种参数 Table 1 Various parameters of self-made and commercial TD TD管 流量/(L·min-1) 温度范围/℃ 解析管尺寸/nm 有效停留时间/s 内径 管长 Burtscher TD 0.3 50~400 20 500 31 Wehner TD 0.6 50~300 20 500 16 Huffman TD 0.6 54~230 22 550 21 TSI 3065 0.2~2 0~400 - - 1~1.5 Dekati 10~20 0~300 10 950 0.3 BH TD 0.1~1.0 0~350 20 500 35.5(0.3 L·min-1)

此外, TSI和Dekati两家公司还提供了两款商业化的TD管.美国TSI公司的3065型TD管气体流速为0.2~2.0 L · min^-1, 最大加热温度为400 ℃, 该TD管主要与该公司生产的扫描电迁移率粒径谱仪(SMPS)配套使用.芬兰Dekati公司的TD管(型号ELA-111、ELA-230)允许流速为10~20 L · min^-1, 最大加热温度为300 ℃, 该TD管通常与该公司生产的荷电低压撞击器(ELPI)或质量监测仪(Dekati Mass Monitor DMM-230)配套使用.上述TD的停留时间相对较小, 如TSI的TD管停留时间为1~1.5 s, 而Dekati公司的TD管停留时间为0.3 s(An et al., 2007).
为了研究大气气溶胶及燃烧产生气溶胶的挥发性, 本研究参考Burtscher TD的结构设计了一款TD管(以下简称为BH TD), BH TD针对上述自制和商业化的热扩散管类产品中存在的温度测量不够精确、有效停留时间过短的一些缺点进行了改进, 并对其性能进行测试与表征, 以期为开展气溶胶的挥发性研究提供参考.
2 TD管的设计(Design of thermal denuder)本研究设计的TD管主要包括解析管和吸附管, 其结构如图 1所示.解析管中间为一根内径20 mm、长500 mm的不锈钢管.TD管采用两处加热带缠绕加热的方式进行加热, 第一处位于解析管不锈钢管外表面, 同时包裹玻璃棉；加热温度由PID控制, 由一根插入到不锈钢管中间的k型热电偶进行温度测量, k型热电偶密封安装在定位支架上, k型热电偶与定位支架内壁存在一定间隔.解析管由室温升温至300 ℃, 所需时间小于10 min, 可控温度范围为50~350 ℃, 此处在下文统称为控温处；第二处位于连接法兰外表面, 同时包裹保温棉, 加热温度由PID控制, 控温热电偶位于连接法兰外表面, 可控温度范围为30~170 ℃, 从室温加热至可控温度范围内任意温度所需时间小于5 min, 此处在下文统称为保温处.吸附管中间是一根不锈钢纱布围成的内径20 mm、长500 mm的圆管, 吸附管外壳与不锈钢纱布圆管中间填充活性炭, 且出口安装测温热电偶, 用于监测出口气体温度.解析管和吸附管内的两根圆管通过连接法兰固定连接, 连接处总长100 mm.连接法兰内部套有高温密封圈, 法兰盘由8根螺钉锁紧.BH TD管的冷却通过泵将压缩空气吹入管内部实现.
图 1(Fig. 1)

图 1 热扩散管结构示意图 Fig. 1Schematic diagram of thermal denuder

本设计的优点在于：①加热的控温热电偶位于不锈钢管中间, 可直接监测流过解析管内部气体温度, 相较于其他将测温点布置在不锈钢外壁的设计, 测量更准确；②在解析管与吸附管连接法兰处单独缠绕加热带并且配备温度控制仪, 可以有效减少散热, 使解析管出口仍然保持着较高的温度, 挥发的气溶胶不会在解析管出口处由于存在较大的温度下降而再凝结；③吸附管出口处设有测温点, 可对出口温度进行监测.
3 TD管的性能测试(Performance test of thermal denuder)3.1 TD管温度分布为了确定解析管和吸附管内的温度分布, 将BH TD加热到设定温度, 保温处也加热到100、150、200、250、300、350 ℃, 用一根k型热电偶分别从解析管入口和吸附管出口沿着轴线水平插入管内进行测量.不同设定温度下测量得到的温度分布结果如图 2a所示, 气体温度在距离入口125 mm的时候达到了设定温度允许范围(±10%).在设定温度分别为100、150和200 ℃的情况下, 解析管出口温度依然保持在设定温度的允许范围内；在设定温度分别为300、350 ℃时, 解析管出口温度分别为260、310 ℃, 略低于设定温度允许范围.在进入吸附管之前, 温度已经基本上降到80 ℃以下, 对活性炭的吸附性影响不大.
图 2(Fig. 2)

图 2 0.3 L · min-1时不同设定温度下BH TD管内的温度分布(a)及不同体积流量下TD管内温度分布(b) Fig. 2Temperature distribution in BH TD tube at different set temperatures at 0.3 L · min-1(a) and temperature distribution in TD tube at different volume flow(b)

本研究设计的BH TD在后续研究中主要与AMS、SMPS和VTDMA联用, 考虑到这些仪器常用的流量及保证足够的停留时间, 本文选择0.1~1.0 L · min^-1的体积流量范围.本研究考察了在此流量范围内, 将温度设定为300 ℃时的BH TD管内的温度分布.从图 2b中可以看到, 体积流量越大, 距离入口处较近的地方(0~250 mm)温度上升越缓慢, 而出口处的温度分布和体积流量之间不存在显著的相关关系；在不同的流速下, BH TD管内出口温度维持在一个比较稳定的水平.由图 2b可知, 在0.1~1.0 L · min^-1的体积流量范围内, 体积流量变化对TD管内的温度分布影响不大, TD管能够满足该流量范围内的加热要求.但体积流量不同, 气溶胶在BH TD管内停留的时间也不同, 体积流量越小, 停留时间越长, 挥发越充分.
将本文设计的BH TD管与Burtscher TD管在同样条件下(加热温度280 ℃, 体积流量0.3 L · min^-1)进行了管内温度分布比较, 结果见图 3.从图 3可以看出, 两者的温度分布曲线在前半部分比较接近, 但在后半部分存在较大的差异；本文设计的TD管后半部分温度分布相对较为平缓, 而Butcher TD管出口温度下降显著；到达解析管出口时, 两者温度差距可达100 ℃.由此可知, 本文设计的BH TD管在0.3 L · min^-1的工作状况下, 相较于Butcher TD管, 出口能保持更好的温度分布, 减少了解析管出口处可能发生的再冷凝.本文设计的解析管出口处包裹的加热带由PID控制, 也可根据需要调整出口的气体温度, 使得出口的温度分布更加灵活.此外, 本文设计的TD管有效加热距离比Burtscher TD管的有效加热距离长约100 mm, 使得气溶胶的有效停留时间增长, 能够更好地解析出挥发性有机物.在室温下, 当管内气体体积流量为0.3 L · min^-1时, 气体在管内的停留时间为35.5 s, 相较于同等流速下Burtscher TD管的31 s, 停留时间有了一定的提升.
图 3(Fig. 3)

图 3 BH TD管与Burtscher TD管的温度分布比较 Fig. 3Comparison of temperature distribution between BH TD and Burtscher TD

3.2 BH TD管的颗粒物传输效率3.2.1 BH TD管的颗粒物传输效率测量颗粒物在TD管的传输过程中, 由于沉降、热泳及扩散等现象会导致颗粒物损失.已有研究表明, 由于沉降导致的细颗粒物损失可以忽略(Burtscher et al., 2001).测量TD管的颗粒物损失对于后期实验数据的校正具有重要意义, 为此本文搭建了如图 4所示的测量系统.气溶胶雾化发生器(Model 3076, TSI Inc.)产生多分散的气溶胶, 先经过硅胶扩散干燥管将相对湿度控制在10%以下, 再经过不同设定温度下的BH TD管, 然后经过中和器, 其次通过差分迁移率分析仪(Differential Mobility Analyzer, Model 3081, TSI Inc, 以下简称DMA)筛选得到不同粒径的单分散气溶胶粒子, 最后进入凝结核粒子计数器(Condensation Particle Counter, Model 3775, TSI Inc, 以下简称CPC)得到此粒径下的气溶胶数浓度N₁, 与经过旁路室温时此粒径下的颗粒物数浓度N进行对比, 由此计算出传输效率p=N₁/N, 即可得到不同粒径不同温度下的颗粒物传输效率.
图 4(Fig. 4)

图 4 TD管的颗粒物损失测量系统流程图 Fig. 4Flow chart of particulate loss measurement system for thermal denuder

由于NaCl在比较宽的粒径范围内具有热稳定性, 其挥发温度达600 ℃, 因此, 选择NaCl气溶胶进行颗粒物损失测量.实验采用超纯水和NaCl(分析纯)配制0.1 g · L^-1的NaCl溶液.
在实验测量之前, 对气溶胶雾化发生系统发生的NaCl气溶胶进行了稳定性检验(徐学哲等, 2016).气溶胶雾化发生器在2 h内连续发生NaCl气溶胶, 每隔0.5 h测5组, 一共测4次.2 h内测得的NaCl粒子谱分布几乎完全重合, 气溶胶浓度最大差值不超过5%, 稳定性较好, 表明气溶胶雾化发生系统可用于BH TD管颗粒物传输效率的测量.
在0.3 L · min^-1的体积流量下TD管的颗粒物损失测量结果如图 5a所示, 图中的点为测得的实际传输效率, 实线为拟合的传输效率曲线.从图中可以看到, 在粒径小于60 nm时, 传输效率随粒径增大呈上升趋势；当颗粒物粒径大于60 nm时, 传输效率近似于一个常数；随着TD管设定温度的不断升高, 颗粒物传输效率不断减小.25 ℃的情况下, 管内的颗粒物损失主要是扩散损失造成的.颗粒物粒径为25 nm时, 扩散损失约为20%；当粒径超过60 nm时, 扩散损失维持在10%左右.当管内温度从25 ℃逐渐升高时, 扩散损失会随之增加；此时, 气体从解析管向吸附管传输的过程当中, 气体温度高于管壁温度, 气体和管壁之间存在温度梯度, 此时气体中的颗粒物便会受到热泳力作用, 向着管壁方向移动, 造成热泳损失, 且热泳损失随着温度梯度增大而增加.通过理论计算发现, 350 ℃的情况下热泳损失在25 nm时损失最大, 可达25%左右, 而由于温度升高导致的扩散损失只增加了2%左右(Hinds, 1982), 因此, 随着温度的升高, 颗粒物损失主要源于热泳损失.BH TD管内温度为350 ℃时, 粒径为25 nm的颗粒物损失高达46%；粒径超过60 nm时, BH TD管颗粒物损失为36%左右.除去热泳损失和扩散损失之外, 还存在一小部分沉降损失, 沉降损失随着颗粒物粒径的增大而增大.本文设计的BH TD管主要传输600 nm以下的颗粒物, 在600 nm时, 理论估算颗粒物的沉降损失最多有5.6%(Hinds, 1982).
图 5(Fig. 5)

图 5 不同温度下TD管的传输效率(a)及多种设计的TD管颗粒物传输效率比较(b) Fig. 5Particulate loss for TD at different temperatures(a) and comparison of particulate loss in various TD designs(b)

3.2.2 多种TD管颗粒物传输效率比较本文还比较了BH TD管与文献中其它TD管颗粒物的传输效率(图 5b), 其中, Huffman TD、Wehner TD管采用0.6 L · min^-1的体积流量, Burtscher TD、BH TD采用0.3 L · min^-1的体积流量.在比较传输效率之前, 通过理论扩散损失计算得到了同种设定温度不同体积流量下不同粒径(25~600 nm)颗粒物传输效率的差值, 发现0.6 L · min^-1体积流量下的传输效率要比0.3 L · min^-1体积流量下的传输效率要高, 这是因为气体流速快、扩散损失小导致的；在25 nm时传输效率相差最大, 在2%左右(Hinds, 1982).从图 5b中可以看到, 粒径大于60 nm的颗粒物在不同设计的TD管中传输效率保持相对稳定, 没有明显的上升趋势.本文设计的BH TD管与Burtscher TD管相比, 在常温和100 ℃下传输效率略低于Burtscher TD管；在250 ℃和350 ℃的情况下, BH TD管的传输效率均比Burtscher TD管要高.
根据图 5a的实际测量结果可知, 当颗粒物的粒径大于60 nm时, 在一定温度下的颗粒物在TD管内的传输效率变化很小, 可近似看作是一个常数.通过建立颗粒物传输效率与温度之间的关系, 由此可获得特定温度下的颗粒物传输效率.图 6线性拟合了实际传输效率与温度之间的关系, 并与Burtscher等(2001)、Wehner等(2002)、Huffman等(2008)拟合的线性关系式进行了比较.从图 6可以看出, Burtscher TD管的线性拟合度最好, 其次是Huffman TD管, Wehner TD管的线性拟合度最差；Wehner TD管的斜率最大, Huffman TD管的截距最大；BH TD管的传输效率在250 ℃以前是最低的, 在250 ℃以后高于Burtscher TD管, 这与图 5中的结果相对应.
图 6(Fig. 6)

图 6 颗粒物粒径大于60 nm时平均传输效率与温度的关系 Fig. 6Relationship between average transport efficiency and temperature when particulate matter is larger than 60 nm

3.3 BH TD管的挥发性能测试(NH₄)₂SO₄是大气气溶胶中主要的无机组分, 在较低温度下几乎不挥发, 在中等温度范围内呈半挥发性, 在较高温度下完全挥发, 具有良好的重复性, 因此, 本文选择(NH₄)₂SO₄表征BH TD管的挥发性能.实验采用超纯水和(NH₄)₂SO₄ (分析纯)配制0.1 g · L^-1的(NH₄)₂SO₄溶液.通过气溶胶雾化发生器(Model 3076, TSI Inc.)产生稳定的(NH₄)₂SO₄气溶胶, 进入到扩散干燥管中将相对湿度控制在10%以下, 然后经过不同设定温度下的BH TD管, 最后进入SMPS进行粒径测量.以250 ℃下SMPS测量得到的(NH₄)₂SO₄质量浓度作为量化依据, 100、150、250、350 ℃下(NH₄)₂SO₄质量浓度分别剩余78%、12.5%、10%、7%.图 7为(NH₄)₂SO₄粒子挥发谱分布, 从图中可以看到, 随着温度的升高, (NH₄)₂SO₄的峰值粒径向左逐渐偏移, 且数浓度呈下降趋势；当BH TD管温度达到350 ℃时, (NH₄)₂SO₄的浓度几乎为0, 这说明在此温度下的(NH₄)₂SO₄粒子已经基本上挥发完全, 这与之前的研究结果一致(Brooks et al., 2002；Gysel et al., 2002；钱小东, 2017).(NH₄)₂SO₄分解后的产物既可能附着在活性炭管上, 也可以冷凝在气溶胶颗粒上；Saleh等(2011)利用无量纲参数Cn表征气相中的蒸气被液滴清除的速率与其被冷却段壁吸收的速率的比值, 当Cn < 0.7时, 标志着冷凝现象不明显, 即分解后的产物大部分被活性炭管吸附.本研究计算得到TD管的Cn远小于0.7, 因此, 可以判断(NH₄)₂SO₄分解后的产物被活性炭管吸附.
图 7(Fig. 7)

图 7 (NH4)2SO4粒子挥发谱分布 Fig. 7Distribution of particle volatilization spectrum of ammonium sulfate

4 结论(Conclusions)1) 本研究设计了一种TD管, 主要由解吸管和吸附管两部分组成；解吸管为加热区域, 调节温度范围为50~350 ℃；脱附管为冷却吸附区域, 装有活性炭, 可避免挥发性物质再冷凝和形成新的粒子；采用两段加热的方式, 加热处与已有设计相比, 改进了测量TD管内部温度的位置, 使得测温、控温更加准确；保温处使解析管的出口温度较高, 避免了解析管出口处的二次冷凝, 且不会影响活性炭的吸附性能.
2) 对设计的BH TD管进行了性能测试, 发现BH TD管的温度分布在不同流量下较均匀；同时, 获得了不同温度下BH TD管的传输效率, 并建立了颗粒物传输效率与温度之间的关系；气溶胶中的(NH₄)₂SO₄气溶胶的挥发性能表明, BH TD管的挥发性能良好.

参考文献

An W J, Pathak R K, Lee B H, et al. 2007. Aerosol volatility measurement using an improved thermodenuder:Application to secondary organic aerosol[J]. Journal of Aerosol Science, 38(3): 305-314. DOI:10.1016/j.jaerosci.2006.12.002

Brooks B J, Smith M H, Hill M K, et al. 2002. Size-differentiated volatility analysis of internally mixed laboratory-generated aerosol[J]. Journal of Aerosol Science, 33(4): 555-579. DOI:10.1016/S0021-8502(01)00192-6

Burtscher H, Baltensperger U, Bukowiecki N, et al. 2001. Separation of volatile and non-volatile aerosol fractions by thermodesorption:Instrumental development and applications[J]. Journal of Aerosol Science, 32(4): 427-442. DOI:10.1016/S0021-8502(00)00089-6

Chen Y.2012.Characterization of carbonaceous aerosols from biofuel combustion: emissions and climate relevant properties[D]. Urbana, IL: University of Illinois at Urbana-Champaign

Goldstein I S.2017.Organic Chemicals From Biomass[M].Boca Raton: CRC Press, Inc

Gysel M, Weingartner E, Baltensperger U. 2002. Hygroscopicity of aerosol particles at low temperatures.2.Theoretical and experimental hygroscopic properties of laboratory generated aerosols[J]. Environmental Science & Technology, 36(1): 63-68.

Hinds W C.1982.Aerosol Technology: Properties, Behavior, and Measurement of Airborne Particles[M]. New Jersey: John Wiley & Sons, Inc

Housiadas C, Drossinos Y. 2005. Thermophoretic deposition in tube flow[J]. Aerosol Science and Technology, 39(4): 304-318. DOI:10.1080/027868290931069

黄聪妮, 黄晓锋, 王少霞, 等. 2015. 华北地区夏季大气亚微米颗粒物的半挥发性特征[J]. 环境科学研究, 28(9): 1345-1352.

Huffman J A, Ziemann P, Jayne J, et al. 2008. Development and characterization of a fast-stepping/scanning thermodenuder for chemically-resolved aerosol volatility measurements[J]. Aerosol Science and Technology, 42(5): 395-407. DOI:10.1080/02786820802104981

Huffman J A, Docherty K S, Aiken A C, et al. 2009. Chemically-resolved aerosol volatility measurements from two megacity field studies[J]. Atmospheric Chemistry & Physics, 18(1): 7161-7182.

Jonsson A M, Hallquist M, Saathoff H. 2007. Volatility of secondary organic aerosols from the ozone initiated oxidation of α-pinene and limonene[J]. Journal of Aerosol Science, 38(8): 843-852. DOI:10.1016/j.jaerosci.2007.06.008

Kanakidou M, Seinfeld J H, Pandis S N, et al. 2005. Organic aerosol and global climate modelling:a review[J]. Atmospheric Chemistry and Physics, 5(4): 1053-1123. DOI:10.5194/acp-5-1053-2005

Kumar P, Fennell P, Symonds J, et al. 2008. Treatment of losses of ultrafine aerosol particles in long sampling tubes during ambient measurements[J]. Atmospheic Environmental, 42: 8819-8826. DOI:10.1016/j.atmosenv.2008.09.003

李园园, 黄晓锋, 曾立武, 等. 2015. 基于热扩散管的深圳大气气溶胶半挥发性分析[J]. 中国环境科学, 35(5): 1281-1287.

钱小东.2017.大气气溶胶VH-TDMA装置建立及吸湿和挥发特性研究[D].合肥: 中国科学技术大学

钱小东, 张启磊, 徐学哲, 等. 2017. 气溶胶吸湿和挥发特性测量的VH-TDMA装置研究[J]. 中国环境科学, 37(4): 1269-1275.

Saha P K, Khlystov A, Grieshop A P. 2015. Determining aerosol volatility parameters using a "Dual Thermodenuder" System:Application to laboratory-generated organic aerosols[J]. Aerosol Science and Technology, 49(8): 620-632. DOI:10.1080/02786826.2015.1056769

Saleh R A, Shihadeh, Khlystov A. 2011. On transport phenomena and equilibration time scales in thermodenuders[J]. Atmospheric Measurement Techniques, 4: 571-581. DOI:10.5194/amt-4-571-2011

Villani P, Picard D, Michaud V, et al. 2008. Design and validation of a volatility hygroscopic tandem differential mobility analyzer (VH-TDMA) to characterize the relationships between the thermal and hygroscopic properties of atmospheric aerosol particles[J]. Aerosol Science and Technology, 42(9): 729-741. DOI:10.1080/02786820802255668

Wang J, Flagan R C, Seinfeld J H. 2002. Diffusional losses in particle sampling systems containing bends and elbows[J]. Aerosol Science, 33: 843-857. DOI:10.1016/S0021-8502(02)00042-3

Wehner B, Philippin S, Wiedensohler A. 2002. Design and calibration of a thermodenuder with an improved heating unit to measure the size-dependent volatile fraction of aerosol particles[J]. Journal of Aerosol Science, 33(7): 1087-1093. DOI:10.1016/S0021-8502(02)00056-3

Wehner B, Philippin S, Wiedensohler A., et al. 2004. Variability of non-volatile fractions of atmospheric aerosol particles with traffic influence[J]. Atmospheric Environment, 38: 6081-6090. DOI:10.1016/j.atmosenv.2004.08.015

徐学哲, 赵卫雄, 方波, 等. 2016. 标准气溶胶发生系统的建立与性能评估[J]. 环境科学学报, 36(7): 2355-2361.