張晶晶，李 明，姜 杰，喬曉林，張冬冬

(1. 哈爾濱工業(yè)大學(威海)，山東 威海 264209；2. 中國計量科學研究院，北京 100013)

基于有限元法的離子遷移譜電場數(shù)值模擬研究

張晶晶1，李 明2，姜 杰1，喬曉林1，張冬冬1

(1. 哈爾濱工業(yè)大學(威海)，山東 威海 264209；2. 中國計量科學研究院，北京 100013)

傳統(tǒng)線性離子遷移譜中漂移管電場的均勻性是離子遷移譜性能的主要影響因素之一。本研究利用COMSOL Multiphysics仿真平臺，基于有限元法對線性漂移管的內(nèi)電場進行數(shù)學物理建模及數(shù)值模擬分析，給出建模方法和均勻電場的標準，研究了電極厚度、絕緣環(huán)厚度、漂移管內(nèi)徑以及漂移管外部電場等因素對漂移管內(nèi)電場的影響，旨在通過模擬結果優(yōu)化離子遷移譜儀器設計，提高離子遷移譜的性能，為設計者提供參考。

離子遷移譜；電場；有限元；COMSOL Multiphysics

離子遷移譜法(ion mobility spectrometry, IMS)是一種痕量化學物質分析方法，主要通過氣態(tài)離子的遷移率來鑒別化學物質，達到對物質分析的目的[1]。IMS方法具有可在大氣壓條件下工作、檢測靈敏度高、分析時間短、體積小、重量輕和功耗低等優(yōu)點，在化學戰(zhàn)劑、毒品和爆炸物檢測等領域得到了廣泛應用，在環(huán)境、醫(yī)學、生物等領域的應用也有了長足的發(fā)展[2]。

漂移管是IMS的核心部分，起到分離離子的功能，它的結構設計對IMS的性能有直接的影響。線性漂移管一般由一系列電極環(huán)和絕緣環(huán)交替排列而成，通過等值電阻連接起來，形成漂移管電場，均勻的電場分布是離子遷移譜儀獲得高分辨、高靈敏度信號的保障。因此，通過電場模擬的結果對儀器設計進行指導變得尤為重要。

近年來，通過計算機模擬計算來指導漂移管結構的優(yōu)化已有報道。Soppart等[3]利用有限差分模擬程序研究了漂移管內(nèi)電場的分布，指出電極環(huán)的厚度與內(nèi)徑之比越小越好，對漂移管的設計有了初步的理論指導。時迎國等[4]利用電場模擬的結果優(yōu)化設計了IMS儀器的結構。Liu等[5]通過模擬電極環(huán)、絕緣環(huán)的厚度，外徑及形狀對電場均勻性的影響優(yōu)化了儀器結構。楊杰等[6]和Han等[7]利用SIMION軟件對傳統(tǒng)線性漂移管進行了多物理場模擬。以上研究大多針對自制儀器進行了電場模擬和參數(shù)優(yōu)化，并沒有給出全面的幾何設計對電場均勻性影響的一般性規(guī)律；此外上述工作大多基于有限差分法進行電場求解，雖然該方法相對比較直觀，計算速度快，但它難以處理復雜的、不規(guī)則的求解域，而有限元法更易實現(xiàn)對偏微分方程的離散，適合處理復雜區(qū)域，并能保證計算精度。

COMSOL Multiphysics是一款基于有限元法的高級數(shù)值仿真軟件，它能高精度的求解偏微分方程(單物理場時)或偏微分方程組(多物理場的情況)，在COMSOL Multiphysics中可根據(jù)需要將邊界區(qū)域的網(wǎng)格剖分加密，因此邊界處的物理量計算更精確，進而實現(xiàn)真實物理現(xiàn)象的模擬。目前，COMSOL Multiphysics已在聲學[8]、光學[9]、化學反應[10]、地球科學[11]、電磁場[12]等領域有著廣泛的應用。Barth等[13-14]和Cumeras等[15-16]利用COMSOL Multiphysics軟件對非對稱電場漂移管進行了多物理場模擬，而基于COMSOL Multiphysics對線性電場漂移管的模擬研究未見報道。本研究利用COMSOL Multiphysics軟件基于有限元法對傳統(tǒng)線性電場漂移管進行系統(tǒng)模擬，并總結關鍵參數(shù)對電場均勻性的影響規(guī)律，為IMS的設計提供參考。

1 電場數(shù)值模擬的基本理論及方法

IMS內(nèi)電場問題屬于靜電場問題，可以通過求解泊松(poisson)方程來實現(xiàn)。泊松方程的表達式為：

▽2·φ=-ρ/ε

(1)

式中φ為電勢，ρ為自由電荷密度，ε=ε0εr，εr為介質的相對介電常數(shù)，ε0為真空介電常數(shù)(ε0=8.545·10-12F/m)。

若求解域內(nèi)沒有自由電荷，則泊松方程就簡化為拉普拉斯方程：

▽2·φ=0

(2)

為了求解上述方程，還需給定求解區(qū)域邊界上的物理情況，即需要指定邊界條件。本研究利用狄利克雷(dirichlet)邊界條件，即給定邊界面上各點的電勢φ=φ0，在計算機的輔助下求得漂移管內(nèi)電場的數(shù)值解。

COMSOL Multiphysics擁有大量預定義的應用模式，本研究選用靜電應用模式。定義相關參數(shù)后，建立的幾何模型截面示意圖示于圖1。

圖1 漂移管幾何模型示意圖Fig.1 Schematic diagram of drift tube geometry model

設定電極環(huán)、絕緣環(huán)的材料分別為不銹鋼和聚四氟乙烯，漂移管內(nèi)外介質為空氣，指定各電極環(huán)的電勢φn(n=1,…,N，N為電極環(huán)的個數(shù))。此外，為了觀測漂移管電場分布情況，在漂移管外添置一個矩形作為求解域，矩形的4個邊界條件設置為“零電荷/對稱”。網(wǎng)格剖分方式選用默認的“自由剖分三角形網(wǎng)格”，在此基礎上手動控制最大單元和最小單元尺寸。利用其他軟件模擬電場時，常因無法精密刻畫電極邊界形狀而使模擬的精度受到影響，而在COMSOL Multiphysics中可以根據(jù)需要將邊界區(qū)域的網(wǎng)格剖分加密，示于圖2。因此，邊界處的物理量計算更精確，這也是基于有限元方法計算出的電場分布較有限差分法精確的原因之一。

由于此模型求解的是單物理場穩(wěn)態(tài)問題，選用最基本的穩(wěn)態(tài)線性求解器。

圖2 網(wǎng)格剖分局部放大圖Fig.2 Partial enlargement drawing of meshing sequence

2 模擬結果

本研究基于COMSOL Multiphysics數(shù)值仿真平臺分別模擬電極厚度、絕緣環(huán)厚度、漂移管內(nèi)徑以及漂移管外部電場等因素對漂移管內(nèi)部電場的影響。因為反應區(qū)的電場分布規(guī)律與漂移區(qū)的一致，為了簡化說明，本工作重點研究漂移區(qū)的電場分布情況。如未做特殊說明，下面用到的幾何參數(shù)列于表1。

2.1 電極厚度對電場分布的影響

模擬的幾何模型的坐標設置示于圖1，漂移區(qū)截面圖的水平方向定為X軸，縱向為Y軸，單位均為cm。根據(jù)IMS的工作原理，在設計漂移管時盡量使電場在X方向的值趨于一個定值，而在Y方向盡量為零。

表1 漂移管幾何模型參數(shù)列表

為了直觀的表現(xiàn)漂移管的幾何設計參數(shù)對電場均勻性的影響，在本研究中，人為定義一個衡量電場均勻性的標準，即某區(qū)域內(nèi)的電場強度(以250 V/cm為例)變化的相對偏差在4%內(nèi)，認為電場在該區(qū)域內(nèi)的分布是均勻的。那么通過數(shù)值模擬可以得到漂移管所對應的均勻電場邊界值y0，使得在0≤|y|≤y0范圍內(nèi)電場滿足下式條件：

|Ey|≤10 V/cm，

240 V/cm≤Ex≤260 V/cm

(3)

其中，Ey為電場強度的縱向分量；Ex為電場強度的水平分量。

電極厚度分別為0.1、1、3 mm時，Ey分布對比圖示于圖3，易知電極厚度越小，y0值越大，即電場均勻性的范圍越大。

圖3 不同電極厚度時電場y分量的對比圖Fig.3 Comparison chart of the electric field y-component for different thicknesses of electrode

3種條件下，Ex在距軸心線0.83 cm平行線上的分布對比示于圖4。由圖4可以看出，電極越薄，圍繞理想值振蕩幅度越小,電場均勻性越好。

圖4 不同電極厚度時電場x分量分布的對比圖Fig.4 Comparison chart of the electric field x-component for different thicknesses of electrode

為了量化模擬結果，定義比值η1和η2：

η1=y0/r，η2=S0/S=η12

(4)

其中，r為電極環(huán)半徑，單位cm，本節(jié)r為1 cm；S0為均勻電場的面積，S0=πy02；S為漂移管內(nèi)腔的橫截面積，S=πr2。η1和η2越大，均勻電場所占比例越大。通過數(shù)值模擬求出絕緣環(huán)厚度為3 mm，電極內(nèi)徑為20 mm時，對應圖3所示的3種電極厚度3、1、0.1 mm，η2所占的百分比分別為44.89%、68.89%、92.16%。

數(shù)值模擬結果說明，僅從電極厚度對電場分布影響角度考慮，設計漂移管時，在允許的前提下，應盡量減小電極厚度。

2.2 絕緣環(huán)厚度對電場分布的影響

絕緣環(huán)的厚度也是一個關鍵參數(shù)，本節(jié)研究不同絕緣環(huán)厚度漂移管內(nèi)電場的均勻狀況。電極厚度為1 mm，電極內(nèi)徑為20 mm，當絕緣環(huán)厚度分別為1、3、6 mm時，所對應的η2值分別為77.44%、68.89%、60.84%，可知絕緣環(huán)厚度越小，均勻電場的比例越大。

上述3種絕緣環(huán)厚度所對應的電場x分量在距軸心線0.83 cm平行線上的分布對比示于圖5。由圖5可知，隨著絕緣環(huán)厚度的增加，電場在x方向上的振蕩幅度變大，均勻電場的范圍變小。所以，在保證有效絕緣效果的前提下，絕緣環(huán)的厚度越小越好。

圖5 不同絕緣環(huán)厚度時電場x分量分布對比圖Fig.5 Comparison chart of the electric field x-component for different thicknesses insulating ring

2.3 漂移管內(nèi)徑對電場分布的影響

漂移管內(nèi)徑對內(nèi)部電場分布有重要影響，在設計時需結合功能、體積綜合考慮。通過模擬結果可見：漂移管內(nèi)徑越大，η2值越大，均勻電場的范圍越廣，結果列于表2。表2中的數(shù)據(jù)是在固定漂移管壁厚(電極外徑與內(nèi)徑差的一半)的前提下得到的。假設壁厚10 mm，此壁厚可保證表2中的結果是在無外部電場干擾下得到的。

如果需要在漂移管內(nèi)徑較小的情況下獲得更大的電場均勻區(qū)域，可以通過減小電極和絕緣環(huán)厚度來實現(xiàn)。漂移管內(nèi)徑1 cm，電極厚度0.1 mm，絕緣環(huán)厚度3 mm條件下，η2值較漂移管內(nèi)徑3 cm，電極厚度1 mm，絕緣環(huán)厚度3 mm條件下的還大?？梢?，雖然漂移管內(nèi)徑較小，由于使用了薄電極，同樣可以獲得較大范圍的均勻電場。由表2還可知，當電極厚度、絕緣環(huán)厚度和電極內(nèi)徑之比等于1∶1∶10時，計算出的η2值均為60.84%，這說明同比例放大或縮小漂移管，不改變電場分布的均勻性。另外結合2.2結果可以看出，當電極厚度、絕緣環(huán)厚度和電極內(nèi)徑之比等于1∶6∶20時，η2值也為60.84%，說明可以通過減小絕緣環(huán)厚度來彌補漂移管內(nèi)徑較小時電場分布均勻性變差的現(xiàn)象。

2.4 漂移管外部電場對內(nèi)部電場的影響

一般離子遷移譜的漂移管外部為接地的金屬外殼，需要考慮漂移管壁厚(電極環(huán)內(nèi)外徑之差的一半)、絕緣環(huán)厚度以及外殼到漂移管的距離等參數(shù)，使外殼不會影響到漂移管內(nèi)部的電場。本研究在數(shù)值模擬時將金屬外殼與漂移管電極的間距固定在3.0 mm，電極與外殼之間的最大電勢差為2 000 V。在無外殼或絕緣環(huán)厚度較小時，內(nèi)電場的均勻性沒有受到干擾，示于圖6a和6c；而當絕緣環(huán)厚度較大，漂移管壁厚較小時，外殼會對漂移管的內(nèi)電場有較大影響，電場的均勻性變差，示于圖6b。漂移管壁厚固定，絕緣環(huán)厚度不同時，Ex在距軸心線0.83 cm平行線上的分布對比示于圖7。由圖7可見,漂移管內(nèi)電場抗外界干擾能力與漂移管幾何設計參數(shù)有關。

表2 電極厚度、絕緣環(huán)厚度及電極內(nèi)徑之比不同時，η2值的結果對比

注：a. 無外殼；b. 外殼接地；c. 外殼接地(優(yōu)化幾何參數(shù))圖6 等勢線分布對比圖Fig.6 Comparison chart of equipotential line distribution

圖7 漂移管壁厚度固定為1 mm，絕緣環(huán)厚度不同時Ex分布對比圖Fig.7 Comparison chart of the electric field x-component for different thicknesses of insulating ring

當漂移管電極電勢低于2 000 V時，漂移管壁厚與絕緣環(huán)厚度相等時，即可保證內(nèi)電場不受外部電場的干擾；否則需要增大漂移管壁厚與絕緣環(huán)厚度的比值，才能保證內(nèi)電場不受干擾。

3 結論

本研究將IMS中的電場問題抽象成數(shù)學模型，使用偏微分方程對其進行描述，借助COMSOL Multiphysics軟件基于有限元方法求解得到了IMS漂移管中的電場分布，進而分別模擬研究了電極厚度、絕緣環(huán)厚度、漂移管內(nèi)徑和漂移管外部電場等因素對漂移管內(nèi)部電場均勻性的影響。采用量化的方式定義了穩(wěn)定區(qū)域，并且以此為依據(jù)總結出了以下規(guī)律：

1) 電極越薄，內(nèi)徑越大，絕緣環(huán)厚度越小，電場的均勻性越好；

2) 電極環(huán)內(nèi)徑較小時，可采用薄電極短間距排列的方式；

3) 同比例放大或縮小漂移管幾何尺寸不改變電場分布；

4) 漂移管的接地外殼會對內(nèi)電場產(chǎn)生影響，為了避免外部電場的干擾，應選擇合適的漂移管壁厚和絕緣環(huán)厚度。

[1] EICEMAN G A, KARPAS Z. Ion mobility spectrometry[M]. Boca Raton: CRC Press, 2005.

[2] 張東風, 孔德義, 梅 濤, 等. 離子遷移譜儀微型化的現(xiàn)狀與進展[J]. 儀器儀表學報, 2006, 27(2): 199-204. ZHANG Dongfeng, KONG Deyi, MEI Tao, et al. Present status and prospect of the miniaturized ion mobility spectrometry [J]. Chinese Journal of Scientific Instrument, 2006, 27(2):199-204 (in Chinese).

[3] SOPPART O, BAUMBACH J I. Comparison of electric fields within drift tubes for ion mobility spectrometry [J]. Measurement Science and Technology, 2000, 11(10): 473-479.

[4] 時迎國,杜明娟. 離子遷移譜儀遷移管中電場均勻性的研究[J]. 科學技術與工程, 2010,10(19): 4 752-4 754. SHI Yingguo, DU Mingjuan. Research of the uniformity of the electric field in ion mobility spectrometer drift tube [J]. Science Technology and Engineering, 2010, 10(19): 4 752-4 754 (in Chinese).

[5] LIU X, LI S L, LI M S. Optimization design of drift tube for ion mobility spectrometer based on simulation of drift electric field [J]. International Journal of Ion Mobility Spectrometry, 2012,15:231-237.

[6] 楊 杰, 曹樹亞, 郭成海, 等. 離子遷移(IMS)漂移管的多物理場模擬方法研究[J]. 現(xiàn)代科學儀器, 2011, 6: 41-46. YANG Jie, CAO Shuya, GUO Chenghai, et al. A method of multi-physics simulation for drift tube in ion mobility spectrometer (IMS) [J]. Modern Scientific Instruments, 2011, 6: 41-46 (in Chinese).

[7] HAN F, DU Y, CHENG S, et al. Computational fluid dynamics-monte carlo method for calculation of the ion trajectories and applications in ion mobility spectrometry[J]. International Journal of Mass Spectrometry, 2012, 309: 13-21.

[8] STEPANENKO D A, MINCHENYA V T. Deve- lopment and study of novel non-contact ultrasonic motor based on principle of structural asymmetry [J]. Ultrasonics, 2012, 52(7): 866-872.

[9] XIANGRU W, CAIDONG X, XIA W, et al. New side-pump scheme for large mode area fiber lasers by wrapping polished fibers conjugate caustic removed[J]. Optics Communications, 2010, 283(1): 44-48.

[10] NAVARRETE A, MATO R B, COCERO M J, et al. A predictive approach in modeling and simulation of heat and mass transfer during microwave heating. Application to SFME of essential oil of lavandin super [J]. Chemical Engineering Science, 2012， 68(1): 192-201.

[11] BUTLER S L, SINHA G. Forward modeling of applied geophysics methods using comsol and comparison with analytical and laboratory analog models[J]. Computers & Geosciences, 2012, 42: 168-176.

[12] JALAAL M, SOLEIMANI S, DOMAIRRY G, et al. Numerical simulation of electric field in complex geometries for different electrode arrangements using meshless local MQ-DQ method[J]. Journal of Electrostatics, 2011, 69(3): 168-175.

[13] BARTH S, BAETHER W, ZIMMERMANN S. System design and optimization of a miniaturized ion mobility spectrometer using finite-element analysis [J]. IEEE Sensors Journal, 2009, 9(4): 377-382.

[14] BARTH S, BAETHER W, ZIMMERMANN S. Concept of a miniaturized ion mobility spectrometer and a numerical model for fast system design and optimization[C]. Proceedings of IEEE Sensors， 2007:79-82.

[15] CUMERAS R, GRACIA I, FIGUERAS E, et al. Finite-element analysis of a miniaturized ion mobility spectrometer for security applications [J]. Sensors and Actuators B: Chemical, 2012, 170: 13-20.

[16] CUMERAS R, GR'ACIA I, IVANOV P, et al. Comsol simulation of acetone ions in planar ion mobility spectrometer[C]. Proceedings of the 2009 Spanish Conference on Electron Devices. Santiago de Compostela, Spain, 2009:323-326.

Numerical Simulation of Electric Field for Ion Mobility Spectrometer Based on Finite Element Method

ZHANG Jing-jing1, LI Ming2, JIANG Jie1, QIAO Xiao-lin1, ZHANG Dong-dong1

(1.SchoolofInformationandElectricalEngineering,HarbinInstituteofTechnologyatWeihai,Weihai264209,China;2.ChemicalMetrology&AnalyticalScienceDivision,NationalInstituteofMetrology,Beijing100013,China)

The homogeneity of electric field in the traditional drift tube is very critical to the performance of ion mobility spectrometry. In this study, numerical simulation of electric field in traditional drift tube was conducted by means of COMSOL Multiphysics software with finite element method. Additionally, modeling method and homogeneity standard for electric field were also presented. The influences of some important factors such as thickness of electrode and insulating ring, inner diameter of drift tube as well as metal enclosure of drift tube, on the homogeneity of electric field were researched. On the basis of the research, some useful laws were concluded. This paper is aimed at optimizing instrument design through the numerical simulation results, which can help designers to improve the performance of ion mobility spectrometry. Compared with the traditional research methods through experimental structure design, the simulation results can help shorten the instrument development cycles and save development costs.

ion mobility spectrometry; electrical field; COMSOL multiphysics； finite element method

2013-05-08；

2013-06-21

國家自然科學基金(21205023)；山東省自然科學基金(ZR2010BQ011)資助

張晶晶(1983～)，女(漢族)，河北獻縣人，博士研究生，信息與通信工程專業(yè)。E-mail: summer2002198320@163.com

姜 杰(1978～)，男(漢族)，吉林四平人，副教授，從事分析儀器的開發(fā)與應用研究。E-mail: hitjiangjie@gmail.com

時間：2014-01-23； 網(wǎng)絡出版地址：http://www.cnki.net/kcms/doi/10.7538/zpxb.youxian.2014.0002.html

O 657.63

A

1004-2997(2014)03-0232-06

10.7538/zpxb.youxian.2014.0002