強(qiáng)流離子束傳輸和加速束流動(dòng)力學(xué)

1、973課題一例會(huì),趙亞亮2014.05.07,目錄,參與任務(wù)情況任務(wù)完成情況代碼開(kāi)發(fā)平臺(tái)的建立已有代碼功能分析新開(kāi)發(fā)程序的物理要求程序開(kāi)發(fā)進(jìn)展問(wèn)題下一步工作計(jì)劃,參與的總體任務(wù)總體，包括物理要求、軟件集成（王志軍、趙亞亮）代碼開(kāi)發(fā)平臺(tái)建立和維護(hù)（趙亞亮）CPU模塊 （趙亞亮+兩個(gè)研究生）本次例會(huì)前應(yīng)完成的任務(wù)已有軟件的構(gòu)架和功能分析；（趙亞亮、齊新、張小虎、王志軍）提出物理要求，包含待開(kāi)發(fā)的功能模塊和發(fā)展計(jì)劃

2、；提出基本的軟件構(gòu)架。,參與任務(wù)情況,任務(wù)完成情況,代碼開(kāi)發(fā)平臺(tái)的建立完成本地的svn服務(wù)器的搭建，以及帳號(hào)、密碼、權(quán)限的設(shè)置目前已開(kāi)發(fā)的代碼已加入到svn版本庫(kù)，擁有賬號(hào)、密碼的用戶可以通過(guò)svn checkout https://192.168.36.121/svn/Track_t/ Track_t  --username user檢索到本地，也可以通過(guò)https://zhaoyl/svn/Track_t/ 從網(wǎng)頁(yè)

3、上打開(kāi)代碼更改后上傳之前，要先update到最新版本再上傳，發(fā)生沖突時(shí)先將沖突引導(dǎo)本地，確認(rèn)最終版本后再上傳建議：每次代碼更改，在更改位置標(biāo)明作者、時(shí)間，并簡(jiǎn)單描述意圖；版本上傳時(shí)，也簡(jiǎn)單描述更改內(nèi)容,,,所外登陸：vpn http://www.ihep.cas.cn/jgsz/kyxt/div7/fuwu/suonei/vpn/,,任務(wù)完成情況,已有代碼功能分析PADSCMatrix跟蹤（quad、drift、gap）多粒

4、子跟蹤（field map）匹配（根據(jù)lattice周期結(jié)構(gòu)，得到與之匹配的Twiss）均溫優(yōu)化（尚未用多粒子跟蹤校驗(yàn)）其他輔助功能（計(jì)算任意兩點(diǎn)之間的相移等）LinacOptT-code PIC算法（GPU程序）BeamPath,,,任務(wù)完成情況,已有代碼功能分析T-code PIC算法（GPU程序）GPU程序常用加速元件的多粒子跟蹤（四極鐵、drift、gap、RFQ）數(shù)據(jù)后處理,,已有代碼功能分析BeamPa

5、thend-to-end 模擬 （ 輸出方式：z-code , t-code） 接受度計(jì)算 變步長(zhǎng)計(jì)算 變頻率計(jì)算束流反向演化組合誤差分析本征發(fā)射度計(jì)算Particle-in-Cell 空間電荷效應(yīng)極化束模擬 RFQ for beam funneling,任務(wù)完成情況,新開(kāi)發(fā)程序的物理要求(參考TraceWi n)圖形界面矩陣跟蹤多粒子跟蹤誤差分析匹配軌道矯正（可屬于匹配）圖形化數(shù)據(jù)后處理O

6、thers。。。,束流跟蹤元件：四極鐵、螺線管；加速間隙、腔；Circular or rectangular aperture；束測(cè)元件；漂移段；場(chǎng)分布（靜電場(chǎng)、靜磁場(chǎng)、RF電磁場(chǎng)及其疊加）；Funneling gap；RFQ ；Space charge compensation,,束流匹配目標(biāo)：算法Owner（？）單純形法（線性規(guī)劃算法）依據(jù)矩陣跟蹤結(jié)果多粒子（partran）跟蹤結(jié)果,,,誤差分

7、析Input beam errorsbend errorscavity errorsRFQ errorsQuadrupole errors,數(shù)據(jù)后處理：在某一位置或者元件出口的粒子分布單位長(zhǎng)度相移、周期相移束流密度分布 Hofmann圖 磁鐵梯度或者強(qiáng)度Field_map viewer矯正子強(qiáng)度、周期長(zhǎng)度、腔功率、能量及其增益變化、加速場(chǎng)、腔壓、接受度、field map factor、流強(qiáng)、包絡(luò)、發(fā)射度、束損、

8、束暈、剝離損失率、tune depression,,程序開(kāi)發(fā)進(jìn)展Matrix跟蹤，元件包括drift、gap、quad、solenoid，元件矩陣參考Trace3d和TraceWin遺傳算法庫(kù)（nsga）的引入和用于匹配時(shí)的驗(yàn)證依賴的庫(kù)：CLHEP（2.1.3），nsga（2），F(xiàn)FTW（多粒子跟蹤）,,已有程序架構(gòu),Probe類,用于記錄束流狀態(tài)，以及Twiss參數(shù)，線性空間電荷力,粒子類,,,beam類,記錄束流粒子狀態(tài)，用于

9、多粒子跟蹤,用于記錄單粒子信息,,Twiss參數(shù),,,,,,已有程序架構(gòu),spacecharge類,FFT算法類,其他,橢圓邊界,開(kāi)放邊界,,,,包含的類輸入文件的格式Element文件用于遺傳算法參數(shù)的輸入文件,,,,用遺傳算法(nsga)進(jìn)行束流匹配研究全局、多目標(biāo)、多變量算法原理：參考遺傳學(xué)理論，首先產(chǎn)生一個(gè)種群，種群中個(gè)體按照一定的優(yōu)劣標(biāo)準(zhǔn)和算法進(jìn)行優(yōu)劣排序，找到這一代個(gè)體中最優(yōu)的，通過(guò)交叉、變異產(chǎn)生下一代，依

10、次類推，找到最優(yōu)的個(gè)體。 變異率越高，收斂越慢；變異率太低容易收斂到局部最優(yōu)而不是全局最優(yōu),,用遺傳算法進(jìn)行束流匹配研究可行性驗(yàn)證：跟蹤一段加速結(jié)構(gòu)（D+G+Q），得到出口的Twiss參數(shù)（ProbeF），將其中一個(gè)（或多個(gè)）四極鐵的梯度G變?yōu)?，以ProbeF的Twiss參數(shù)作為的目標(biāo)，以加速段中的quad作為變量，用nsga進(jìn)行匹配，可以很容易找到G。目標(biāo)函數(shù)：2個(gè)：橫向失配參數(shù)取得最小值限制條件：2個(gè)：橫向失配參數(shù)小于0

11、.1變量個(gè)數(shù)：4個(gè)：4塊四極鐵梯度缺點(diǎn)：循環(huán)迭代結(jié)束條件比較單一,,將Q4設(shè)置為0，初始種群大小是100，經(jīng)過(guò)100代遺傳，可以找到失配參數(shù)小于0.1的解,,,Q2、Q3、Q4=0,Q1、Q2、Q3、Q4=0,Q4=0,Q3、Q4=0,問(wèn)題,T-code 多粒子跟蹤，以元件為單位是否合理？,下一步工作計(jì)劃,將該程序遺傳算法用于ADS加速段之間的匹配針對(duì)大家對(duì)已有結(jié)構(gòu)提出的意見(jiàn)和建議，盡快完善矩陣計(jì)算部分開(kāi)始進(jìn)行多粒子跟蹤部分,,