1 500 MHz 5-cell超導(dǎo)腔的仿真優(yōu)化

楊 峻 馬震宇 王 巖 侯洪濤 李 正 羅 琛毛冬青 是 晶 劉建飛

1 500 MHz 5-cell超導(dǎo)腔的仿真優(yōu)化

楊 峻1,2,3馬震宇1,3王 巖1,3侯洪濤1,3李 正1,3羅 琛1,3毛冬青1,3是 晶1,3劉建飛1,3

1（中國科學(xué)院上海應(yīng)用物理研究所 嘉定園區(qū) 上海 201800）
2（中國科學(xué)院大學(xué) 北京 100049）
3（上海市低溫超導(dǎo)高頻腔技術(shù)重點(diǎn)實(shí)驗(yàn)室 上海 201800）

能量回收直線加速器技術(shù)是近年來一項(xiàng)重要技術(shù)。針對能量回收直線加速器的主直線加速器，本文使用多種射頻軟件優(yōu)化設(shè)計(jì)了一種新型的1 500 MHz 5-cell超導(dǎo)腔。該腔加速模TM010的R/Q值達(dá)到550?，可以有效降低低溫?fù)p耗。該腔的端cell采用擴(kuò)大束管傳輸有害高次模方案將高次模傳輸出超導(dǎo)腔，良好地抑制了高次模。新腔型保持了較低的表面峰值場強(qiáng)與平均加速梯度的比值：Epk/Eacc=2.06，Bpk/Eacc=4.22 mT/(MV·m?1)，避免了場致發(fā)射的風(fēng)險(xiǎn)，并且在中等加速梯度區(qū)域(15?20 MV·m?1)無二次電子倍增現(xiàn)象。模擬結(jié)果顯示，該超導(dǎo)腔能夠滿足連續(xù)波、準(zhǔn)連續(xù)波或高重復(fù)頻率穩(wěn)定運(yùn)行的要求。

加速器，超導(dǎo)腔，高次模，模擬計(jì)算

射頻超導(dǎo)技術(shù)是目前及將來先進(jìn)加速器的關(guān)鍵技術(shù)之一。尤其是在新一代高能加速器、連續(xù)波自由電子激光(Continuous Wave Free Electron Laser, CW-FEL)或強(qiáng)流能量回收直線加速器(Energy Recovery Linac, ERL)等設(shè)施上，射頻超導(dǎo)技術(shù)因?yàn)槠淇蓪?shí)現(xiàn)高效率粒子加速的特性更具有不可取代的優(yōu)勢。目前加速器領(lǐng)域的前沿研究機(jī)構(gòu)都相繼采用射頻超導(dǎo)技術(shù)。如美國的LCLS-II (Linac Coherent Light Source-II)，其已計(jì)劃采用1 300 MHz的超導(dǎo)腔來實(shí)現(xiàn)4 GeV、0.3 mA高重復(fù)頻率（兆赫茲以上）的束流[1]；歐洲的自由電子激光EURO-XFEL (European X-ray Free Electron Laser)，其已采用了1300 MHz超導(dǎo)腔實(shí)現(xiàn)脈沖型硬X射線FEL出光[2]；還有美國CORNELL大學(xué)的ERL光源，在多年預(yù)研的基礎(chǔ)上其研究機(jī)構(gòu)確定了適合強(qiáng)流的1 300 MHz 7-cell超導(dǎo)腔的腔型結(jié)構(gòu)[3]。

中國科學(xué)院上海應(yīng)用物理研究所已經(jīng)在深紫外自由電子激光方面取得了顯著的成果[4]，不久必將開展連續(xù)波或者高重復(fù)頻率運(yùn)行的自由電子激光技術(shù)研究。為實(shí)現(xiàn)連續(xù)波自由電子激光的高效率粒子加速，相關(guān)的超導(dǎo)腔是其唯一的選擇。結(jié)合已有的研究基礎(chǔ)，本文主要設(shè)計(jì)了一種可工作在ERL型FEL裝置上的1 500 MHz的5-cell超導(dǎo)腔SH-ERL (Cavity for Shang Hai ERL)，并通過模擬結(jié)果表明，該超導(dǎo)腔能夠滿足連續(xù)波或高重復(fù)頻率穩(wěn)定運(yùn)行的要求。

1 頻率選擇

超導(dǎo)高頻腔的工作頻率對腔體的性能參數(shù)和體積大小有巨大影響，因此在頻率的選擇上需綜合考慮加速性能與制造、運(yùn)行成本。在頻率的選擇上首先參照當(dāng)年TESLA的報(bào)告，其顯示1 500 MHz為最優(yōu)頻率[5]，但是TESLA由于其1 300 MHz的功率源技術(shù)更成熟，選擇了1 300 MHz作為工作頻率。其次，固定平均加速梯度下，單位長度的熱負(fù)載與頻率成反比，如下：

式中，Pdiss為熱損耗；Eacc為平均加速梯度；R/Q為腔特性阻抗，只與腔形狀有關(guān)。選擇高頻率可以降低單位長度的熱損耗。采用的諧振頻率越高，腔體表面積越小，腔的成本因加工所消耗鈮材減少而降低。同時更小的腔體積也使得高溫處理和相關(guān)測試的操作相對簡單。并且，中國科學(xué)院上海應(yīng)用物理研究所已經(jīng)在250 MHz頻段的VHF (Very HighFrequency)電子腔與500 MHz單cell[6]和5-cell[7]的超導(dǎo)腔上開展研究與實(shí)驗(yàn)，它們可用作1 500 MHz超導(dǎo)加速腔的電子源和注入器。從另一方面考慮，選擇相對高的頻率雖然也造成了相對高的表面BCS電阻，但是可以通過降低工作溫度或采用新的表面處理方法[8]，如氮摻雜[9]和高溫處理的方法降低表面BCS電阻?；谝陨弦蛩兀疚倪x擇1 500 MHz作為主加速器的工作頻率。

2 腔型設(shè)計(jì)的目標(biāo)與要求

腔型優(yōu)化希望達(dá)到更小的損耗、更高的加速梯度，并且避免場致發(fā)射等限制超導(dǎo)腔加速性能的因素。超導(dǎo)腔的低溫?fù)p耗可由式(2)給出：

式中，Pdiss為功率損耗；Vc為加速腔壓；Rs為表面電阻。根據(jù)式(2)，G×R/Q值越高[3]，損耗越小。

但是讓所有的高頻(Radio Frequency, RF)參數(shù)同時達(dá)到最優(yōu)值是不可能的，因此對于不同機(jī)器參數(shù)的超導(dǎo)腔，腔的優(yōu)化設(shè)計(jì)需要在保持其他高頻參數(shù)在可以接受的范圍時側(cè)重追求某一高頻參數(shù)。當(dāng)前國際上L波段腔型按照優(yōu)化高頻參數(shù)的不同側(cè)重方向分為三類：(1) 低功率損耗腔型(Low Loss Type)，如CEBAF升級改造項(xiàng)目中的腔型[10]和Sekutowicz為ILC設(shè)計(jì)的腔型[11]；(2) 最小化Epk/Eacc參數(shù)以降低場致發(fā)射可能性的腔型（超導(dǎo)腔在15?20 MV·m?1的中等加速梯度區(qū)間運(yùn)行容易發(fā)生場致發(fā)射[3]，所以運(yùn)行在中等加速梯度區(qū)間的腔型更加側(cè)重最小化Epk/Eacc參數(shù)），如TESLA[12]和CORNELL大學(xué)為ERL設(shè)計(jì)的腔型[3]等；(3) 具有最小化的Bpk/Eacc，以實(shí)現(xiàn)更高加速梯度的腔型，如CORNELL大學(xué)的重入式腔[13]等。

本文設(shè)計(jì)的1 500 MHz 5-cell超導(dǎo)腔運(yùn)行時的加速梯度選擇為15?20 MV·m?1的中等加速梯度區(qū)間，參數(shù)優(yōu)化時側(cè)重最小化Epk/Eacc，以降低場致發(fā)射的風(fēng)險(xiǎn)。綜合考慮下提出目標(biāo)和要求為：(1) Epk/Eacc≤2.1；(2) Bpk/Eacc≤4.3 mT/(MV·m?1)；(3) 盡量高的G×R/Q；(4) 滿足強(qiáng)流ERL對高次模抑制的要求；(5) 曲率半徑大于6 mm（小于6 mm難以加工）[3]。

針對SH-ERL設(shè)計(jì)的具體目標(biāo)和要求，首先使用Superfish對腔型進(jìn)行優(yōu)化設(shè)計(jì)；然后用CST (Computer Simulation Technology)微波工作室來對高次模進(jìn)行分析；最后用Multipac對超導(dǎo)腔的二次電子倍增進(jìn)行檢驗(yàn)。本文最終目的是通過一系列模擬運(yùn)算，使該超導(dǎo)腔能夠滿足連續(xù)波或高重復(fù)頻率穩(wěn)定運(yùn)行的要求。

2.1 中間腔的優(yōu)化設(shè)計(jì)

使用Superfish軟件來進(jìn)行腔型的模擬計(jì)算以初步確定其高頻參數(shù)，為SH-ERL設(shè)計(jì)的1 500 MHz中間腔示意圖如圖1所示。

圖1 中間腔腔型示意圖Fig.1 Schematic of a middle cell.

如圖1所示，Req表示赤道半徑，Rai表示束孔半徑，Ai、Bi、ai、bi分別表示腔型的赤道處橢圓以及束孔處橢圓橫縱半軸，Li表示腔型的半腔長度，α為傾斜角。通過單獨(dú)改變每個幾何參數(shù)，并調(diào)整Req使得基模頻率為1 500 MHz，得到各幾何參數(shù)的變化對高頻參數(shù)的影響如圖2所示。

由圖2(a)所示，R/Q在ai取10?11 mm處會取得一個最大值。腔型模擬結(jié)果在該區(qū)間達(dá)到最小的腔耗，因此ai取值在10?11 mm。由圖2(b)所示，Bpk/Eacc、Epk/Eacc分別在赤道橢圓橢圓率eequator為0.85和0.9附近取得最低值。eequator取值選為0.85?0.9時設(shè)計(jì)出腔型能夠有更高的加速梯度。由圖2(c)所示，束孔孔徑橢圓橢圓率eiris小于1.75時，Epk/Eacc最小，能更好地降低場致發(fā)射的風(fēng)險(xiǎn)。所以eiris取值選在1.65?1.75。由圖2(d)所示，Epk/Eacc在α大于12°時顯著減小。為降低場致發(fā)射的風(fēng)險(xiǎn)，α取值為12°?14°?；谝陨戏治?，中間腔腔型幾何尺寸參數(shù)取值如表1所示。

圖2 ai (a)、eequator (b)、α (c)、eiris (d)變化對超導(dǎo)腔高頻參數(shù)影響Fig.2 RF parameters as a function of ai (a), eequator (b), α (c), eiris (d).

表1 中間腔腔型尺寸參數(shù)Table 1 Dimensions of a middle cell.

2.2 端cell優(yōu)化設(shè)計(jì)

多cell超導(dǎo)腔中，端cell連接中間腔和束管部分。在端cell的設(shè)計(jì)上不僅需要將腔體頻率維持在1 500 MHz，而且要滿足場平坦度良好[14]的要求，還需要考慮如何更好地將高次模(Higher Order Mode, HOM)傳輸出超導(dǎo)腔。另外在具體優(yōu)化過程中，因端cell連接中間腔，Req固定不變，需通過改變Le將頻率調(diào)節(jié)到1 500 MHz。圖3為采用擴(kuò)大束管傳輸高次模方案的端cell腔型示意圖，其中Ae、Be、ae、be分別表示腔型端cell的赤道處橢圓以及束孔處橢圓橫縱半軸，Le表示腔型的半腔長度，at、bt、c、Rbp表示束管過渡段的結(jié)構(gòu)參數(shù)。擴(kuò)大束管傳輸高次模原理為：通過選擇合適的束管半徑，使得束管截止頻率低于高次模的諧振頻率，讓高次模能夠傳出腔體被高次模吸收器吸收。據(jù)式(3)：

式中，fc為束管截止頻率；c為光速；r為束管半徑；X對于TM波為vnm，對于TE波為μnm。最低頻率的高次模頻率約為1 820 MHz，將其代入式(3)計(jì)算得到束管半徑為47.67 mm。

由于不對稱的端cell設(shè)計(jì)有利于將高次模傳出超導(dǎo)腔[15]，本文采用不對稱端cell設(shè)計(jì)結(jié)構(gòu)。幾何尺寸取值如表2所示。

圖3 端cell腔型示意圖Fig.3 Schematic of an end cell.

表2 端cell尺寸參數(shù)Table 2 Dimensions of end cells and beam pipes.

2.3 腔數(shù)目的選擇與高次模分布

優(yōu)化腔型結(jié)構(gòu)時一方面希望盡可能地提高加速模式性能，另一方面需要盡可能地降低高次模對帶電粒子束的不利影響，尤其是強(qiáng)流ERL對高次模的抑制比直線加速器要求更高。隨著腔體cell數(shù)目的增加，雖然腔體有效加速長度增加，但是高次模將會更難被耦合出腔體[16]。所以腔體數(shù)目的選擇需要在更高的有效加速長度與更好的高次模的抑制之間進(jìn)行權(quán)衡考慮。

針對于ERL的特點(diǎn)，參考CORNELL大學(xué)ERL對高次模阻尼的要求：其認(rèn)為如果偶極模的外部品質(zhì)因子(Qe)小于104，四極模的外部品質(zhì)因子小于1010，那么高次模對于束流的影響可以忽略[17]。而在給定的ERL機(jī)器物理參數(shù)下束流崩潰不穩(wěn)定性的閾值電流Ith(式(4))，如需得到百毫安量級的閾值電流以滿足運(yùn)行要求，根據(jù)Cornell大學(xué)提出的判據(jù)[17]，對于偶極模，如式(5)所示，對于四極模，如式(6)：

本文使用CST微波工作室進(jìn)行模擬計(jì)算，建立的5-cell模型如圖4所示，束管端假設(shè)為加載了理想高次模吸收器。

圖4 用來計(jì)算外部品質(zhì)因子的微波工作室模型Fig.4 Microwave studio model for Qe calculation.

通過CST計(jì)算了cell數(shù)目分別為5、7、9的腔典型高次模的外部品質(zhì)因子Qe與特性阻抗R/Q，計(jì)算結(jié)果見表3和圖5。表3給出了cell數(shù)目分別為5、7、9的腔典型高次模外部品質(zhì)因子分布。由表3可以看出，對于同一模式，隨著cell數(shù)目的增大，高次模的外部品質(zhì)因子相應(yīng)升高。雖然表3中所有腔型的四極模的外部品質(zhì)因子都小于1010，但只有cell數(shù)目選擇為5時偶極模的外部品質(zhì)因子都小于104。

表3 不同cell數(shù)高次模的外部品質(zhì)因子Table 3 External quality factor of the higher order modes in cavities with different cell numbers.

腔 cell 數(shù)目選擇為5時(R/Q)×(Qe/f)的分布如圖5所示。由圖5可見，偶極模的(R/Q)×(Qe/f)均小于105，四極模的(R/Q)×(Qe/f)均小于106。所以1 500 MHz的5-cell新腔型的高次模能夠滿足束流崩潰不穩(wěn)定性的百毫安量級電流閾值要求。

圖5 SH-ERL腔偶極模與四極模的(R/Q)×(Qe/f)分布圖Fig.5 (R/Q)×(Qe/f)'s spectra of the dipole higher order modes (HOMs) and the quadrupole HOMs in SH-ERL cavity.

另外，如式(7)定義的腔敏感系數(shù)αs高低會影響腔內(nèi)場平坦度和加工公差的要求高低，其數(shù)值越低越好[18]。

式中，kcc為腔間耦合系數(shù)，由cell尺寸決定；N為cell數(shù)目；加速粒子為電子時=1β。從式(7)可見，對于確定kcc的腔型結(jié)構(gòu)N越大，敏感系數(shù)αs越大。所以綜合考慮有效加速長度、ERL對高次模阻尼的要求以及敏感系數(shù)，本設(shè)計(jì)的超導(dǎo)腔cell數(shù)目選擇為5。

當(dāng)超導(dǎo)腔腔型幾何尺寸及cell數(shù)目確定后，可采用Superfish模擬進(jìn)行計(jì)算，所得到的SH-ERL超導(dǎo)腔加速模TM010-π的電場分布如圖6(a)，相應(yīng)軸向場分布如圖6(b)。通過軸向場分布計(jì)算可得其場平坦度為98%。

圖6 Superfish (field flatness=98%) 1 500 MHz 5-cell腔型的電場分布(a)和軸向場分布(b)Fig.6 Electric field distribution (a) and on-axis field profile (b) of the 1 500 MHz 5-cell cavity calculated by Superfish (field flatness=98%).

通過一系列優(yōu)化設(shè)計(jì)，得到的1 500 MHz 5-cell超導(dǎo)腔的主要高頻性能參數(shù)如表4。將新腔型高頻性能參數(shù)與國際上主流的L波段超導(dǎo)腔進(jìn)行比較，主要高頻性能參數(shù)基本相當(dāng)。SH-ERL腔保持了較低的表面峰值場強(qiáng)與平均加速梯度的比值為Epk/Eacc=2.06，Bpk/Eacc=4.22 mT/(MV·m?1)，同時G×R/Q值較高。

表4 不同腔型高頻參數(shù)對比Table 4 RF parameters comparison of different cavities.

還用ABCI (Azimuthal Beam Cavity Interaction)計(jì)算了不同束團(tuán)長度下腔的縱向損耗因子，不同束團(tuán)長度的縱向損耗因子如圖7所示，可通過縱向損耗因子計(jì)算高次模功率如式(8)所示：

式中，PHOM為高次模功率；k//為縱向損耗因子；ω為角頻率；為加速模式的損耗因子。

圖7 ABCI中SH-ERL腔對于不同束團(tuán)長度的損耗因子Fig.7 Integrated loss factor for the SH-ERL cavity with different bunch lengths calculated by ABCI.

2.4 二次電子倍增

在超導(dǎo)腔的優(yōu)化設(shè)計(jì)過程中還需要考慮一些限制超導(dǎo)腔加速性能的因素，其中二次電子倍增會吞噬高頻能量和引起高能量部分（比如耦合器、窗口、高流強(qiáng)超導(dǎo)腔）崩潰。為檢查二次電子倍增是否發(fā)生，本文使用Multipac 2.1來模擬腔內(nèi)的二次電子倍增效應(yīng)。設(shè)置二次電子碰撞次數(shù)為30時計(jì)算得到：在峰值電場0?40 MV·m?1區(qū)間表征二次電子發(fā)射系數(shù)的參數(shù) e30/c0＜1。該結(jié)果表明在1 500 MHz 5-cell SH-ERL超導(dǎo)腔中不會發(fā)生二次電子倍增。

3 結(jié)語

本文針對SH-ERL上的主直線加速器對超導(dǎo)腔的需求，設(shè)計(jì)了一種1 500 MHz 5-cell新型超導(dǎo)腔。其加速模特性阻抗R/Q =550 ?，并且在良好抑制高次模的同時保持了較低的表面峰值場強(qiáng)Epk/Eacc=2.06，Bpk/Eacc=4.22 mT/(MV·m?1)。主要高頻參數(shù)與國際上主流的L波段超導(dǎo)腔相當(dāng)。

本文還對1 500MHz 5-cell SH-ERL超導(dǎo)高頻腔的高次模進(jìn)行了初步的分析。分析了SH-ERL腔典型偶極模(TE111，TM110)、四極模(TE211，TM210)的外部品質(zhì)因子，判斷出其能達(dá)到ERL高次模抑制要求（偶極模外部品質(zhì)因子Qe大于104，四極模外部品質(zhì)因子Qe大于1010）。并通過進(jìn)一步對高次模的(R/Q)×(Qe/f)值計(jì)算表明新腔型還能夠滿足ERL百毫安量級閾值電流的要求。

1 Galayda J. The LCLS-II project[C]. Processing of Conference TUCA01, 2014

2 Schneidmiller E A, Yurkov M V. Photon beam properties at the European XFEL[R]. Deutsches Elektronen-Synchrotron (DESY), Hamburg, Germany, 2011

3 Valles N. The main linac cavity for Cornell's energy recovery linac: cavity design through horizontal cryomodule prototype test[J]. Nuclear Instruments and Methods in Physics Research A, 2014, 734: 23?31

4 王興濤, 蘭太和, 劉波, 等. 上海深紫外自由電子激光裝置(SDUV-FEL)種子激光傳輸與注入系統(tǒng)[J]. 核技術(shù), 2013, 36(5): 050103

WANG Xingtao, LAN Taihe, LIU Bo, et al. Seed laser transmission and injection system for the Shanghai Deep-Ultraviolet Free-Electron Laser test facility[J]. Nuclear Techniques, 2013, 36(5): 050103

5 Edwards D A. TESLA test facility linac: design report[M]. Version 1.0. Hamburg: Dt. Elektronen-Synchrotron, 1995

6 Liu J F, Hou H T, Mao D Q, et al. Great progress in developing 500 MHz single cell superconducting cavity in China[J]. Science China Physics, Mechanics and Astronomy, 2011, 54(2): 169?173

7 Hou H T, Ma Z Y, Mao D Q, et al. Studies of LL-type 500MHz 5-cell superconducting cavity at SINAP[J]. Chinese Physics C, 2015, 39(4): 047001

8 Ma Z Y, Liu J F, Hou H T, et al. Surface preparation processing for superconducting cavities[J]. Nuclear Science and Techniques, 2014, 25: 060102

9 Grassellino A, Romanenko A, Sergatskov D, et al. Nitrogen and argon doping of niobium for superconducting radio frequency cavities: a pathway to highly efficient accelerating structures[J]. Superconductor Science and Technology, 2013, 26(10): 102001

10 Daly E F, Campisi I E, Henry J, et al. Improved prototype cryomodule for the CEBAF 12 GeV upgrade[C]. Particle Accelerator Conference, Proceedings of the IEEE, 2003, 2: 1377?1379

11 Sekutowicz J, Ko K, Ge L, et al. Design of a low loss SRF cavity for the ILC[C]. Particle Accelerator Conference, Proceedings of the IEEE, 2005: 3342?3344

12 Aune B, Bandelmann R, Bloess D, et al. Superconducting TESLA cavities[J]. Physical Review Special Topics-Accelerators and Beams, 2000, 3(9): 092001

13 Geng R L. Review of new shapes for higher gradients[J]. Physica C: Superconductivity, 2006, 441(1): 145?150

14 Tang Z B, Liu J F, Hou H T, et al. Frequency control and pre-turning of a large aperture 500 MHz 5-cell superconducting RF cavity[J]. Nuclear Science and Techniques, 2014, 25(3): 030102

15 Li Y M, Zhu F, Quan S W, et al. The design of a five-cell high-current superconducting cavity[J]. Chinese Physics C, 2012, 36(1): 74?79

16 Liepe M, Knobloch J. Superconducting RF for energy-recovery linacs[J]. Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, 2006, 557(1): 354?369

17 Valles N, Liepe M. Seven-cell cavity optimization for Cornell's energy recovery linac[C]. Proceedings of SRF, 2009

18 Calaga R R. Linear beam dynamics and ampere class superconducting RF cavities@ RHIC[D]. New York: Stony Brook University, 2006

CLC TL503.2

Simulation of a 1 500 MHz 5-cell superconducting cavity

YANG Jun1,2,3MA Zhenyu1,3WANG Yan1,3HOU Hongtao1,3LI Zheng1,3LUO Chen1,3MAO Dongqing1,3SHI Jing1,3LIU Jianfei1,3
1(Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Jiading Campus, Shanghai 201800, China)
2(University of Chinese Academy of Sciences, Beijing 100049, China)
3(Shanghai Key Lab of Cryogenics & Superconducting RF Technology, Shanghai 201800, China)

Background: Energy Recovery Linacs (ERLs) are being developed over the world to realize future light sources and other applications. Purpose: The aim is to design a new 1 500 MHz 5-cell superconducting cavity for main linac in the ERL. Methods: The new cavity geometry was optimized by several radio-frequency softwares. Results: The new cavity is endowed with the following advantages: with the R/Q of the fundamental mode TM010to be 550 ?, cryogenic losses could be reduced effectively; in order to reduce the risk of field emission, the parameters Epk/Eacc, Bpk/Eaccwere optimized respectively to 2.06 and 4.22; a larger beam pipe was employed to propagate the higher order modes (HOMs) out of the cavity, increasing the damping efficiently for the dangerous HOMs. Conclusion: It was found that radio frequency property of the new cavity could meet the requirements of ERL-based Free Electron Laser (FEL) very well.

Accelerator, Superconducting cavity, Higher order modes, Simulation

TL503.2

10.11889/j.0253-3219.2015.hjs.38.060103

No.11335014)資助

楊峻，男，1990年出生，2008年畢業(yè)于山東大學(xué)，現(xiàn)為碩士研究生，研究領(lǐng)域?yàn)楹思夹g(shù)應(yīng)用

劉建飛，E-mail: liujianfei@sinap.ac.cn

2015-02-11，

2015-04-13