[論文レビュー] Toutatis, the Cea-Saclay RFQ code

R. Duperrier, R. Ferdinand+2 more|arXiv (Cornell University)|Aug 17, 2000

Particle accelerators and beam dynamics参考文献 3被引用数 16

ひとこと要約

TOUTATIS は、CEA-Saclay で開発された高精度でマルチグリッドを用いた RFQ 電磁界ダイナミクスコードであり、PARMTEQM の結果を検証し、連続波（CW）高出力線形加速器における電場精度を向上させることを目的としている。適応メッシュ細分化と正確なヤコビアン保存を実現する時間刻みみ込みアルゴリズムを用い、空間電荷効果、イメージ効果、および複雑なバンプ幾何形状（含む結合ギャップ）を、CERN の実験結果と5%未満の差異、IPHI 設計試験では97%の透過率で正確にシミュレートしている。

ABSTRACT

A CW high power linear accelerator can only work with very low particles losses and structure activation. At low energy, the RFQ is a very sensitive element to losses. To design the RFQ, a good understanding of the beam dynamics is requested. Generally, the reference code PARMTEQM is enough to design the accelerator. TOUTATIS has been written with goals of cross-checking results and obtaining a more reliable dynamics. This paper relates the different numerical methods used in the code. It is time-based, using multigrids methods and adaptive mesh for a fine description of the forces without being time consuming. The field is accurately calculated through a Poisson solver and the vanes are fully described, allowing to properly simulate the coupling gaps and RFQs extremities. Differences with PARMTEQM and LIDOS.RFQ are shown.

研究の動機と目的

• RFQ ビームダイナミクスシミュレーションにおいて、PARMTEQM の信頼性の高い代替手段を提供すること。
• 結合ギャップや端効果を含むバンプ幾何形状を完全にモデル化することで、RFQ の電場精度を向上させること。
• マルチグリッド法を用いて計算時間を数日から数時間に短縮しながら、高解像度を維持すること。
• 非一様ギャップを有するセグメント化バンプ構造を含む複雑な RFQ 幾何形状の正確なシミュレーションを可能とすること。
• 正確な電場モデルにより粒子損失と構造的活性化を最小限に抑えることで、高出力 CW 加速器を支援すること。

提案手法

• 時間領域における粒子-格子（PIC）アルゴリズムを採用し、速度依存のループフロッグスキームによりヤコビアンを正確に保存することで、人工的なエミッタンス増大を回避する。
• 3次元有限差分法に重み付き係数を適用してポアソン方程式を離散化し、ステアケージアーティファクトが生じないバンプ境界の正確な表現を可能にする。
• ガウス＝ザイデル平滑化、縮小、拡大を用いたマルチグリッド V サイクルの緩和法を採用し、ポアソンソルバの収束を高速化する。
• 二次的な高解像度グリッドを用いた適応メッシュ細分化（AMR）を実装し、rms バンチサイズに合わせて解像度を動的に調整する。
• ポアソンソルバに、結合ギャップや電場歪みを正確にモデル化できるように、完全なバンプ幾何形状をディリクレ境界条件として埋め込む。
• PARMTEQM の入力ファイルを直接サポートし、シームレスな移行と検証を可能にする。

実験結果

リサーチクエスチョン

• RQ1 RFQ におけるビームダイナミクスシミュレーションを、PARMTEQM よりもさらに正確かつ信頼性の高いものにするにはどうすればよいか？
• RQ2 複雑な RFQ 幾何形状における高解像度電場分解能を維持しつつ、計算時間を著しく短縮する数値的手法は何か？
• RQ3 結合ギャップはビーム透過率とエミッタンスにどのように影響するのか？ そして、正確にモデル化するにはどうすればよいか？
• RQ4 適応メッシュ細分化は、計算コストが著しく増大するリスクを伴わず、空間電荷効果の解像度を向上させることができるか？
• RQ5 TOUTATIS は、既存のコードと比較して、実験測定結果をどの程度正確に再現できるか？

主な発見

• 65³ および 33³ グリッド解像度において、理論的電場と比較して TOUTATIS は 0.7% 未満の差異を示し、ポアソンソルバの正確性が検証された。
• IPHI RFQ 設計において、ギャップがヤングの最適位置にある場合、TOUTATIS シミュレーションでは 97% の透過率と 12% のエミッタンス増大が得られ、ずれがある場合には 28% のエミッタンス増大が発生する。
• マルチグリッド法を用いることで、計算時間をピボット 450 MHz の環境で1週間から5時間に短縮し、全 RFQ の実用的シミュレーションを可能にした。
• TOUTATIS と CERN 測定値との差異は5%未満であり、PARMULT と比較して15%の差異を示す PARMULT よりも著しく優れている。
• 結合ギャップをシミュレーションに組み込むことは不可欠である。省略すると、エミッタンスと透過率の予測が楽観的になりすぎる。
• 適応メッシュ細分化により、性能コストを最小限に抑えつつ、高解像度の電荷分布モデル化が可能となり、バンチ加速中でも正確性を維持できる。