陸羿行 方文程 張俊強(qiáng) 郭玉森,5 趙振堂
1(中國(guó)科學(xué)院上海應(yīng)用物理研究所 上海 201800)
2(中國(guó)科學(xué)院大學(xué) 北京 101408)
3(中國(guó)科學(xué)院上海高等研究院 上海 201204)
4(上海市低溫超導(dǎo)高頻腔技術(shù)重點(diǎn)實(shí)驗(yàn)室 上海 201800)
5(上??萍即髮W(xué) 上海 201210)
射頻四極場(chǎng)加速器(Radio Frequency Qaudrupole,RFQ)[1-4]常用于低能質(zhì)子的加速,可以保證其束流品質(zhì),如今廣泛應(yīng)用于質(zhì)子加速器及相關(guān)裝置中。現(xiàn)有裝置上的RFQ工作頻率通常在200~400 MHz,而歐洲核子中心(European Organization for Nuclear Research,CERN)研發(fā)的用于質(zhì)子治療裝置的RFQ工作頻率為750 MHz[5],這種高頻率的RFQ具有更高的加速梯度,可以推動(dòng)質(zhì)子加速器的小型化。近十多年中,上海光源(Shanghai Synchrotron Radiation Facility,SSRF)也在致力于質(zhì)子治療裝置的研發(fā),為瑞金醫(yī)院腫瘤質(zhì)子中心設(shè)計(jì)、建造的首臺(tái)國(guó)產(chǎn)質(zhì)子治療裝置已經(jīng)通過(guò)驗(yàn)收,進(jìn)入臨床試驗(yàn)階段[6-8];同時(shí),提出了基于直線加速器的質(zhì)子閃療裝置[9-11],緊湊型RFQ作為裝置中的一部分,也已經(jīng)過(guò)了多年的研究。調(diào)諧是緊湊型RFQ的關(guān)鍵技術(shù)之一,傳統(tǒng)的調(diào)諧算法基于腔體理論計(jì)算得到的物理特性,用于緊湊型RFQ效果欠佳。CERN在高頻RFQ上利用實(shí)際場(chǎng)分布對(duì)調(diào)諧棒的響應(yīng)量,即響應(yīng)矩陣,來(lái)進(jìn)行調(diào)諧[12-13],在求解的過(guò)程中應(yīng)用了奇異值分解(Singular Value Decomposition,SVD)方法。然而,這種方法仍然存在不足,在求解過(guò)程中常常出現(xiàn)調(diào)諧量超出調(diào)諧棒可調(diào)節(jié)范圍的情況,即使對(duì)SVD方法進(jìn)行了改造解決了這個(gè)問(wèn)題,又會(huì)造成精度下降,導(dǎo)致調(diào)諧步驟多、二極場(chǎng)分量難以降低等問(wèn)題。
針對(duì)以上問(wèn)題,本文提出了一種在響應(yīng)矩陣的基礎(chǔ)上利用最小二乘法(Least Squares,LS)進(jìn)行求解的調(diào)諧算法。該算法可以限定求解的范圍,得到其中的最優(yōu)解;同時(shí),對(duì)腔體的二極場(chǎng)分量賦予權(quán)重,可以有效地對(duì)二極場(chǎng)進(jìn)行調(diào)諧,解決二極場(chǎng)誤差難以降低的問(wèn)題?;谶@種新型的RFQ調(diào)諧算法,在Matlab下進(jìn)行了算法的開發(fā),完成了模擬條件下的調(diào)諧實(shí)驗(yàn),搭建了基于微擾法的測(cè)量平臺(tái),編寫了控制和通信的Labview程序,并在RFQ樣機(jī)上完成了算法的最終驗(yàn)證。
經(jīng)過(guò)實(shí)驗(yàn)驗(yàn)證,新調(diào)諧算法可以在限定范圍內(nèi)得到最優(yōu)解,并有效地解決二極場(chǎng)誤差。該算法具有普適性,可以應(yīng)用到其他頻率的RFQ中,解決更多RFQ的調(diào)諧問(wèn)題;將來(lái)可以促進(jìn)質(zhì)子加速器的小型化,推動(dòng)醫(yī)用質(zhì)子裝置的發(fā)展,普及相關(guān)診療手段(質(zhì)子治療、質(zhì)子CT(Computed Tomography)[14-15]、硼中子俘獲療法(Boron Neutron Capture Therapy,BNCT)[16]、醫(yī)用同位素生產(chǎn)[17]等)。
四翼型RFQ的腔體是在諧振腔的腔壁上添加4個(gè)極頭形成的,其工作模式是類TE210模。在橫向上,4個(gè)極頭間可以激發(fā)出電四極場(chǎng),腔體中的電場(chǎng)也主要集中在極頭附近;在遠(yuǎn)離極頭的位置,主要為磁場(chǎng),方向沿軸向(圖1(a))。當(dāng)束流在軸線上通過(guò)時(shí),會(huì)受到四極電場(chǎng)的聚焦作用,并且該電場(chǎng)的極性會(huì)隨時(shí)間和粒子的位置變化,因此,每過(guò)半個(gè)射頻周期粒子就會(huì)受到相反的作用。在縱向上,極頭與軸線間的距離最小為a,最大為ma,m為調(diào)制參數(shù)。經(jīng)過(guò)調(diào)制后,極頭呈類三角函數(shù)的形狀,形成縱向的電場(chǎng),實(shí)現(xiàn)對(duì)束流的聚束和加速作用。極頭形狀變化的周期為βλ,相對(duì)極頭的“峰”和“谷”是對(duì)應(yīng)的,而相鄰極頭的“峰”和“谷”則是錯(cuò)開的(圖1(b))。
圖1 RFQ結(jié)構(gòu)及測(cè)量框架圖(a) RFQ橫向形狀及電磁場(chǎng)分布,(b) RFQ極頭縱向形狀,(c) RFQ測(cè)量方法框架圖Fig.1 Structure of an RFQ and the framework of the measurement(a) Transverse section of an RFQ and the electromagnetic field on it, (b) Longitudinal section of an RFQ, (c) Framework of measuring an RFQ
測(cè)量腔體中的場(chǎng)分布是RFQ調(diào)諧的第一步,可以利用微擾法實(shí)現(xiàn)。RFQ的最小束流孔徑一般都非常小,在1~3 mm,難以測(cè)量其軸線上的電場(chǎng)分布,因此,RFQ是通過(guò)測(cè)量4個(gè)象限中的磁場(chǎng)來(lái)得到場(chǎng)分布的。用步進(jìn)電機(jī)帶動(dòng)魚線,魚線又帶動(dòng)金屬小球在腔體中運(yùn)動(dòng),金屬小球會(huì)對(duì)腔體中的電磁場(chǎng)分布產(chǎn)生影響,探針測(cè)量得到的S12信號(hào)也會(huì)發(fā)生改變,由此可以推算出小球所在位置的磁場(chǎng)強(qiáng)度,存在關(guān)系[18]:
式中:B表示磁場(chǎng)強(qiáng)度;ΔP表示小球的擾動(dòng)引起的相位差。電機(jī)控制器和網(wǎng)絡(luò)分析儀可以與計(jì)算機(jī)進(jìn)行通信來(lái)傳遞數(shù)據(jù)(圖1(c))。
加速器調(diào)諧的基本邏輯是利用調(diào)諧器對(duì)腔體的擾動(dòng),彌補(bǔ)腔體內(nèi)場(chǎng)分布的誤差。通過(guò)對(duì)腔體的測(cè)量,得到場(chǎng)分布及其與理想狀態(tài)的差別,再利用腔體的物理特性,計(jì)算得到調(diào)諧器所需的調(diào)諧量。RFQ的主要調(diào)諧方法可以分為兩大類:一種基于理論計(jì)算得到的物理特性,例如文獻(xiàn)[19],可以稱為“理論法”;另外一種基于實(shí)際測(cè)量得到的物理特性,例如文獻(xiàn)[12-13],可以稱為“解析法”。理論法的缺點(diǎn)在于把理想狀態(tài)下計(jì)算得到的數(shù)據(jù)應(yīng)用于實(shí)際,忽略了實(shí)際腔體和理想模型之間的差別;對(duì)于敏感度更高的緊湊型RFQ而言,實(shí)際情況與理想狀態(tài)相差較大,因此理論法難以適用。解析法則可以消除理論和實(shí)際之間的誤差對(duì)調(diào)諧造成的影響,可以解決緊湊型RFQ的調(diào)諧問(wèn)題;解析法的缺點(diǎn)是需要額外測(cè)量響應(yīng)矩陣。下面將詳細(xì)介紹基于響應(yīng)矩陣和最小二乘法的新調(diào)諧算法。
分別用B1、B2、B3、B4來(lái)表示各象限中的磁場(chǎng)(圖1(a)),則可以計(jì)算得到腔體中的二、四極場(chǎng):
式中:Q為四極場(chǎng);Ds、Dt為兩個(gè)方向上的二極場(chǎng)。在實(shí)際測(cè)量和計(jì)算中,以上場(chǎng)量均使用歸一化值,其歸一化誤差定義為:
將二、四極場(chǎng)組合起來(lái)形成場(chǎng)矩陣:
式中:n為沿軸線取場(chǎng)的點(diǎn)數(shù),并滿足場(chǎng)矩陣中元素的數(shù)量大于等于調(diào)諧器數(shù)量,即3n≥ntuner,因?yàn)檫@樣可以提高調(diào)諧的精度。RFQ初始狀態(tài)的場(chǎng)矩陣記為Fi,接著分別移動(dòng)每一個(gè)調(diào)諧器并重復(fù)測(cè)量對(duì)應(yīng)的場(chǎng)分布,記為F1,F(xiàn)2,…,F(xiàn)n-tuner,可以得到不同場(chǎng)分量對(duì)調(diào)諧器的響應(yīng),即響應(yīng)矩陣:
式中:d為調(diào)諧器的移動(dòng)量。對(duì)于一個(gè)理想的RFQ,需要四極場(chǎng)的值都相等,且二極場(chǎng)大小為0,即目標(biāo)矩陣:
因此,調(diào)諧器對(duì)場(chǎng)分布的補(bǔ)償量應(yīng)為:
至此,調(diào)諧問(wèn)題轉(zhuǎn)化為求解方程:
式中:T即為需要求解的調(diào)諧量矩陣。
從上文可以發(fā)現(xiàn),響應(yīng)矩陣M通常為一個(gè)非正方形矩陣,因此無(wú)法直接求解。CERN使用SVD方法來(lái)求解,但該方法無(wú)法限定調(diào)諧量的范圍,直接求解得到的調(diào)諧量非常大,超出調(diào)諧棒可調(diào)諧的范圍,因此,在求解過(guò)程中將特征值向量中的部分元素代換為0。這種改造后的SVD方法雖然解決了結(jié)果的范圍問(wèn)題,但是舍去部分特征值會(huì)影響到結(jié)果的精度,并可能放大二極場(chǎng)誤差難以減小的問(wèn)題。
對(duì)SVD方法做嘗試時(shí)發(fā)現(xiàn)了上述問(wèn)題。針對(duì)結(jié)果的范圍問(wèn)題,可以利用最小二乘法對(duì)調(diào)諧量矩陣T中的每一個(gè)元素進(jìn)行單獨(dú)求解,得到這一元素在限定范圍內(nèi)使誤差最小的值,可以用式(9)表示:
式中:tn為矩陣T中的第n個(gè)元素。依次求解矩陣T中的元素,并對(duì)這個(gè)過(guò)程進(jìn)行迭代,直到所有誤差達(dá)到預(yù)期的目標(biāo)為止。
在這種方法中,還可以對(duì)二極場(chǎng)分量賦予權(quán)重,以解決二極場(chǎng)誤差難以減小的問(wèn)題,即將原方程變?yōu)榍蠼猓?/p>
式中:k為給二極場(chǎng)賦予的權(quán)重。
整個(gè)調(diào)諧算法如圖2所示。以上為調(diào)諧實(shí)驗(yàn)過(guò)程中針對(duì)一些問(wèn)題給出的新解法,下面從模擬和實(shí)驗(yàn)兩方面闡述對(duì)新調(diào)諧算法的驗(yàn)證。
圖2 基于響應(yīng)矩陣和最小二乘法的調(diào)諧算法流程圖Fig.2 Flowchart of the tuning algorithm based on response matrix and least squares method
為驗(yàn)證新調(diào)諧算法,首先利用三維電磁場(chǎng)模擬軟件進(jìn)行了大量模擬調(diào)諧實(shí)驗(yàn)。RFQ的三維模型基于上海光源設(shè)計(jì)的714 MHz RFQ,在正式設(shè)計(jì)中,該RFQ的長(zhǎng)度為2 m,可以將質(zhì)子從50 keV加速至4 MeV,而用于模擬實(shí)驗(yàn)的模型長(zhǎng)為1 m,橫截面形狀與正式設(shè)計(jì)相似。RFQ的實(shí)驗(yàn)?zāi)P腿鐖D3(a)所示,為金屬模型,在該腔體上共有16根調(diào)諧棒,如圖中所指部分,其余為耦合器、真空口及探測(cè)器預(yù)留位置。為實(shí)驗(yàn)?zāi)P椭匦略O(shè)計(jì)了底切,在所有調(diào)諧棒均插入5 mm深的狀態(tài)下,將其頻率調(diào)至714 MHz,且將場(chǎng)分布調(diào)至理想狀態(tài)。讀取腔體4個(gè)象限中的磁場(chǎng)分布,在不同象限中讀取的位置是對(duì)稱的,其中一個(gè)位置如圖3(b)所示。對(duì)場(chǎng)分布處理后,得到腔體的二、四極場(chǎng)分布如圖3(c)所示,即為理想狀態(tài)。圖中四極場(chǎng)上的“x”標(biāo)記為后續(xù)調(diào)諧實(shí)驗(yàn)過(guò)程中,選取場(chǎng)值用于調(diào)諧算法的位置,即n=6。
圖3 模擬條件下RFQ的模型及場(chǎng)分布(a) RFQ模型,(b) RFQ截面圖及選取場(chǎng)分布位置示意,(c) RFQ二、四極場(chǎng)分布(Qo為四極場(chǎng),Ds、Dt分別為水平和豎直方向的二極場(chǎng))Fig.3 3-D model of the RFQ and its field distribution in simulation(a) Model of the RFQ, (b) Transverse section of the RFQ and the place where the field was picked, (c) Dipole and quadrupole components of the RFQ
根據(jù)算法在Matlab中編寫了相應(yīng)的程序。在理想狀態(tài)下,將每根調(diào)諧棒的插入深度分別減少和增加3 mm,重復(fù)獲取模型中的場(chǎng)分布,得到每根調(diào)諧棒在6 mm范圍內(nèi)對(duì)場(chǎng)的影響,計(jì)算出的響應(yīng)矩陣為18×16的矩陣。
首先進(jìn)行單調(diào)諧棒擾動(dòng)模擬。在理想狀態(tài)下,將第一根調(diào)諧棒t1的深度增加2 mm,對(duì)場(chǎng)產(chǎn)生擾動(dòng),破壞理想狀態(tài),通過(guò)基于響應(yīng)矩陣和LS方法來(lái)對(duì)場(chǎng)進(jìn)行模擬調(diào)諧實(shí)驗(yàn)。
在單根調(diào)諧棒擾動(dòng)的情況下,雖然在算法中加入了對(duì)調(diào)諧量范圍的限制,但是在實(shí)際計(jì)算的過(guò)程中并未出現(xiàn)解出的調(diào)諧量大于調(diào)諧棒可調(diào)諧范圍的情況。在給二極場(chǎng)誤差加上不同的權(quán)重后,對(duì)比得到的結(jié)果可以發(fā)現(xiàn)(圖4),隨著權(quán)重的增加,調(diào)諧棒t1的解越接近-2 mm;而不論權(quán)重如何變化,其他調(diào)諧棒解出的調(diào)諧量絕對(duì)值均非常小,都在0.1 mm以內(nèi)。由于在數(shù)值上該結(jié)果就已經(jīng)十分接近擾動(dòng)量,因此沒(méi)有再將該結(jié)果應(yīng)用到模型中進(jìn)行模擬。基本可以推論,僅通過(guò)一次調(diào)諧,這些誤差就已經(jīng)大幅降低。
圖4 單調(diào)諧棒擾動(dòng)下不同權(quán)重得到的調(diào)諧量Fig.4 Solution with the different weight for dipole components under single tuners' perturbation
通過(guò)單調(diào)諧棒擾動(dòng)下的模擬調(diào)諧,初步驗(yàn)證了調(diào)諧算法的可行性,并且實(shí)現(xiàn)了較為精準(zhǔn)的調(diào)諧效果。
單調(diào)諧棒擾動(dòng)實(shí)驗(yàn)初步驗(yàn)證了算法的準(zhǔn)確性,但這種誤差情況與實(shí)際不符,因此進(jìn)一步進(jìn)行了多調(diào)諧棒擾動(dòng)模擬。通過(guò)在[-3 mm,3 mm]范圍內(nèi)生成隨機(jī)數(shù),為每一根調(diào)諧棒賦予了“擾動(dòng)”,作為調(diào)諧的初始狀態(tài)。調(diào)諧過(guò)程中的歸一化誤差用下式定義:
初始的二、四極場(chǎng)分布如圖5(b、c)所示,初始的二極場(chǎng)誤差較大,接近30%。
圖5 模擬調(diào)諧實(shí)驗(yàn)結(jié)果及對(duì)比(a) 不同權(quán)重下得到的二、四極場(chǎng)誤差比較,(b) 調(diào)諧過(guò)程中四極場(chǎng)分布比較,(c) 調(diào)諧過(guò)程中二極場(chǎng)分布比較Fig.5 Result comparisons of the tuning experiment in simulation(a) Quadrupole and dipole component errors with different weights, (b) Quadrupole field variation during the tuning,(c) Dipole field variation during the tuning
同樣利用LS調(diào)諧算法來(lái)進(jìn)行調(diào)諧。首先比較了二極場(chǎng)誤差不同權(quán)重下的結(jié)果(圖5(a)),圖5中q為四極場(chǎng)誤差,ds、dt分別為水平和豎直方向的二極場(chǎng)誤差。與單調(diào)諧棒擾動(dòng)的情況不同,當(dāng)二極場(chǎng)誤差權(quán)重達(dá)到5之后,四極場(chǎng)誤差明顯增長(zhǎng),迅速增長(zhǎng)到2%以上,當(dāng)權(quán)重增加到10以后,誤差達(dá)到了6%。為了詳細(xì)了解權(quán)重對(duì)結(jié)果的影響,補(bǔ)充了權(quán)重為2、3的結(jié)果,可以發(fā)現(xiàn)這兩種情況下,四極場(chǎng)誤差會(huì)有所增長(zhǎng),同時(shí)二極場(chǎng)誤差會(huì)降低,并且在這個(gè)范圍內(nèi),四極場(chǎng)誤差均不超過(guò)2%,二極場(chǎng)誤差也均在2%以內(nèi)。幾種權(quán)重相比而言,2倍的權(quán)重可以在保證較小的二極場(chǎng)誤差的同時(shí),讓四極場(chǎng)誤差在1%左右,是比較好的結(jié)果。以上模擬實(shí)驗(yàn)表明,多調(diào)諧棒擾動(dòng)的情況和單調(diào)諧棒是有明顯區(qū)別的,多調(diào)諧棒擾動(dòng)下進(jìn)行模擬調(diào)諧會(huì)更加接近實(shí)際情況,同時(shí)也說(shuō)明為二極場(chǎng)誤差賦予權(quán)重的方法仍然是有效的。
進(jìn)一步地,在權(quán)重為2的調(diào)諧結(jié)果的基礎(chǔ)上進(jìn)行了第二次調(diào)諧,在第二次調(diào)諧中,賦予二極場(chǎng)誤差的權(quán)重依舊設(shè)為2。兩次調(diào)諧結(jié)果,以及初始場(chǎng)、目標(biāo)場(chǎng)的比較見(jiàn)圖5(b、c)。圖5中,q為四極場(chǎng)誤差,ds、dt分別為水平和豎直方向的二極場(chǎng)誤差,前綴o表示目標(biāo)場(chǎng),i表示初始場(chǎng),l1表示第一次調(diào)諧后結(jié)果,l2表示第二次調(diào)諧后結(jié)果。從圖5中可以發(fā)現(xiàn),經(jīng)過(guò)一次調(diào)諧后,四極場(chǎng)分布已經(jīng)和目標(biāo)十分接近,同時(shí)二極場(chǎng)誤差大大減小。經(jīng)過(guò)兩次調(diào)諧后,二、四極場(chǎng)誤差均又在縮小,最終四極場(chǎng)誤差為0.8%,二極場(chǎng)誤差為0.4%,是非常理想的誤差水平。
經(jīng)過(guò)多調(diào)諧棒擾動(dòng)模擬調(diào)諧實(shí)驗(yàn),完全驗(yàn)證了調(diào)諧算法的可行性,驗(yàn)證了二極場(chǎng)權(quán)重對(duì)消除二極場(chǎng)的作用,預(yù)示著調(diào)諧方法在實(shí)際調(diào)諧中的優(yōu)秀效果。
對(duì)基于SVD解法的調(diào)諧算法也進(jìn)行了簡(jiǎn)單驗(yàn)證。不可避免地,通過(guò)SVD解法得到的調(diào)諧范圍很大,遠(yuǎn)遠(yuǎn)超出了調(diào)諧棒的可移動(dòng)范圍,必須利用改造后的SVD方法。根據(jù)算法,一共可以得到16組結(jié)果,其中前幾組結(jié)果存在同樣的問(wèn)題,而后幾組結(jié)果又存在丟失精度過(guò)多的問(wèn)題,因此挑選了其中第9組來(lái)進(jìn)行調(diào)諧。經(jīng)過(guò)調(diào)諧后場(chǎng)分布誤差有明顯改善,四極場(chǎng)誤差在1%以內(nèi),最大二極場(chǎng)誤差約2%;第二次調(diào)諧依舊選用第9組結(jié)果,經(jīng)過(guò)第二次調(diào)諧后,二極場(chǎng)誤差減?。?.8%),而四極場(chǎng)誤差增大(2%)。通過(guò)SVD方法兩次調(diào)諧后四極場(chǎng)誤差偏大,必須經(jīng)過(guò)更多的調(diào)諧次數(shù)來(lái)減小誤差;而應(yīng)用二極場(chǎng)誤差權(quán)重加LS方法,如上一節(jié)所述,兩次調(diào)諧后實(shí)現(xiàn)了較為理想的場(chǎng)分布。
最終SVD方法和LS方法得到的二、四極場(chǎng)誤差比較如圖6所示,通過(guò)比較結(jié)果,驗(yàn)證了LS方法可限定求解范圍、高效處理二極場(chǎng)誤差、縮短調(diào)諧過(guò)程等優(yōu)點(diǎn)。
圖6 不同調(diào)諧算法下二、四極場(chǎng)誤差對(duì)比Fig.6 Comparison of the quadrupole and dipole component errors with different tuning methods
以上完成了調(diào)諧算法的初步驗(yàn)證,根據(jù)模擬調(diào)諧實(shí)驗(yàn)的RFQ模型,加工了一段鋁模樣段(圖7(a)),利用RFQ樣機(jī)進(jìn)一步驗(yàn)證調(diào)諧方法。
圖7 RFQ樣段調(diào)諧平臺(tái)及實(shí)驗(yàn)結(jié)果(a) 714 MHz RFQ樣機(jī),(b) RFQ測(cè)量平臺(tái),(c) 電機(jī)控制和微波測(cè)量程序面板,(d) RFQ樣段初始場(chǎng),(e) 調(diào)諧過(guò)程中二、四極場(chǎng)誤差變化Fig.7 Testbed of the RFQ prototype and the experimental results(a) Prototype of the 714-MHz RFQ, (b) Testbed of the prototype, (c) User interface of the motor control and RF measurement,(d) Initial field of the RFQ prototype, (e) Variation of quadrupole and dipole component errors during the tuning
RFQ的測(cè)量方法如§1所述,具體平臺(tái)搭建見(jiàn)圖7(b)。在RFQ樣機(jī)的端板上開有4個(gè)測(cè)量孔,分別在RFQ的4個(gè)象限,魚線帶著金屬小球在4個(gè)象限中來(lái)回穿過(guò),并繞過(guò)固定在端板上的滑輪,利用步進(jìn)電機(jī)帶動(dòng)小球在腔體中勻速運(yùn)動(dòng)。通過(guò)腔體上預(yù)留的耦合孔、真空抽口以及探測(cè)器孔將探針伸入腔體內(nèi),探針的另一頭連接到矢量網(wǎng)絡(luò)分析儀,獲取腔體中的微波參數(shù)。步進(jìn)電機(jī)控制及網(wǎng)絡(luò)分析儀的數(shù)據(jù)傳輸均通過(guò)在LabView環(huán)境下編寫的程序來(lái)實(shí)現(xiàn)(圖7(c))。測(cè)量得到的數(shù)據(jù)先進(jìn)行預(yù)處理,再利用調(diào)諧算法完成調(diào)諧量的計(jì)算,數(shù)據(jù)處理及計(jì)算均通過(guò)Matlab程序?qū)崿F(xiàn)。RFQ樣機(jī)上有16個(gè)調(diào)諧器,調(diào)諧器中螺桿帶動(dòng)調(diào)諧頭移動(dòng),可以實(shí)現(xiàn)較為精準(zhǔn)的調(diào)諧量。
對(duì)RFQ樣段的初始場(chǎng)進(jìn)行測(cè)量。先將每個(gè)調(diào)諧器的插入深度都調(diào)節(jié)為5 mm,使用兩根探針測(cè)量腔體內(nèi)的微波參數(shù),處理后得到RFQ樣段的初始場(chǎng)分布,其中,二極場(chǎng)誤差在30%以上(圖7(d))。主要原因在于初始狀態(tài)調(diào)諧器的深度不一致,誤差主要來(lái)源于兩方面:一是調(diào)諧器與腔體間沒(méi)有焊接,導(dǎo)致調(diào)諧器與腔體的間隙不一樣,存在誤差;二是調(diào)諧器的精度未達(dá)到預(yù)想的程度,難以判斷各調(diào)諧器的深度是否一致。以上誤差只會(huì)影響初始誤差,不會(huì)影響調(diào)諧過(guò)程,因此在此樣段上進(jìn)行調(diào)諧實(shí)驗(yàn)來(lái)驗(yàn)證算法及調(diào)諧平臺(tái)的可靠性仍是可行的。
利用上述調(diào)諧算法對(duì)RFQ模型腔進(jìn)行調(diào)諧。初始狀態(tài)下二、四極場(chǎng)誤差分別為24.09%和1.57%,經(jīng)過(guò)5次調(diào)諧后兩者誤差分別降為2.33%和1.39%。此時(shí)誤差已經(jīng)達(dá)到理想的水平,對(duì)腔體進(jìn)行了一次頻率調(diào)諧,即將所有調(diào)諧器的深度增加0.6 mm,二、四極場(chǎng)誤差產(chǎn)生輕微變化,為2.52%和1.34%(圖7(e)劃線處)。之后又進(jìn)行了2次額外的調(diào)諧以驗(yàn)證調(diào)諧算法的收斂性,結(jié)果顯示誤差在小范圍振蕩,可能是誤差較小的結(jié)果已接近調(diào)諧算法的極限精度。調(diào)諧各步的誤差結(jié)果見(jiàn)圖7(e)。
針對(duì)高頻率、緊湊型RFQ在調(diào)諧過(guò)程中遇到的問(wèn)題,提出了一種基于響應(yīng)矩陣和最小二乘法的新調(diào)諧算法,可以在調(diào)諧棒可調(diào)節(jié)范圍內(nèi)求出最優(yōu)解,且為二極場(chǎng)分量增加權(quán)重可以高效地減小二極場(chǎng)誤差。在模擬環(huán)境下進(jìn)行了調(diào)諧實(shí)驗(yàn),初步驗(yàn)證了新調(diào)諧算法的可行性,與SVD解法的比較中體現(xiàn)出了新調(diào)諧算法的優(yōu)勢(shì)。根據(jù)調(diào)諧算法編寫了Matlab程序,基于Labview開發(fā)了RFQ的測(cè)量和調(diào)諧平臺(tái),通過(guò)RFQ樣機(jī)的調(diào)諧實(shí)驗(yàn)進(jìn)一步驗(yàn)證了調(diào)諧算法。這種調(diào)諧算法是針對(duì)高頻率、緊湊型RFQ開發(fā)的,同時(shí)也具備很強(qiáng)的普適性,適用于其他頻率的RFQ,將來(lái)可以促使緊湊型RFQ應(yīng)用在更多質(zhì)子加速器裝置上,推動(dòng)醫(yī)用質(zhì)子裝置的普及。
作者貢獻(xiàn)聲明陸羿行負(fù)責(zé)研究的提出及設(shè)計(jì)、數(shù)據(jù)的收集和整理、文章的起草和最終版的修訂;方文程負(fù)責(zé)項(xiàng)目的監(jiān)督和管理、文章的修訂;張俊強(qiáng)負(fù)責(zé)實(shí)驗(yàn)數(shù)據(jù)的收集;郭玉森負(fù)責(zé)實(shí)驗(yàn)數(shù)據(jù)的收集;趙振堂負(fù)責(zé)項(xiàng)目的監(jiān)督和管理、文章的修訂。