发表VG-核融合量子科学技术研究开发机构.ppt

大型ヘリカル装置における実座標を用いた粒子軌道追跡モンテカルロコードの開発 關 良輔, 松本 裕, 鈴木康浩1), 渡邊清政1) 北海道大学大学院工学研究科, 1)核融合科学研究所 大型ヘリカル装置（LHD） 高ベータ実験（低磁場Bax = 0.5 T） ? NBIによる加熱 NBIに起因するプラズマ圧力（ビーム圧力）の同定が 高ベータプラズマでの平衡, 安定性解析において必要 ? LCFSの内側のみを解析 ?粒子の損失境界は最外殻磁気面(LCFS) はじめに 従来の高エネルギー粒子の軌道をもとにした分布関数, 圧力の解析 LHDでは高エネルギー粒子の軌道が複雑 ? 軌道追跡が特に重要 ビーム圧力の同定 ? 軌道追跡が重要 磁気座標を使用 はじめに（続き） LCFS 周辺磁場領域 Re-entering Re-entering粒子 LCFS外側の周辺磁場領域に出ても再び LCFS内部に戻ってくる粒子 目的 ○Re-entering粒子も考慮可能な, 実座標を用いた粒子軌道追跡に基づくモンテカルロコードの開発 〇高ベータプラズマ 〇低磁場中 ?LCFS外側に出る粒子が増加 Re-entering粒子が重要 ?LCFSが小さい ?周辺磁場領域が厚い ?ドリフトによるズレが大きい 実座標を用いた粒子軌道解析 ○分布関数算出の流れ 開発したモンテカルロコード 熱速度の3倍以下になった粒子が緩和 個の高エネルギー粒子の案内中心を追跡 各位相空間中のmeshにおける滞在時間 を算出 エネルギー緩和, 軌道損失まで追跡 ? 分布関数算出 真空容器壁に衝突した粒子が損失 加熱入力 初期エネルギー エネルギー緩和には, OFMCコードと同様の衝突オペレータを採用 ○δt秒後の速度 , 衝突オペレータ　(Plasma Phys. Control. Fusion 49 (2007) 1955) 平均: 分散: , drift-kinetic方程式の定常解を求めるコードを作成 ○drift-kinetic方程式 高エネルギー粒子源 衝突項 損失項 熱化による吸収項 ○温度Tbのマクスウェル分布を持つ 　b種粒子に衝突したa種粒子δt秒間の速度変化 c: pitch angle, h: gyro-phase は, 以下の平均と分散を持つ正規乱数により算出 横長断面上の磁束 の磁気面上 ○使用磁場 ○初期ピッチ角 ○初期位置 ○モンテカルロ粒子数 1,000 テスト計算1（衝突オペレータのみと開発したコードの比較） のみ与えた ○粒子の初期エネルギー 200 keV ○粒子 陽子 ○back ground(軽水素プラズマ) 温度 1 keV 密度 1020 m-3 開発したコード 速さ-分布関数 ピッチ角-分布関数 テスト計算1（衝突オペレータのみと開発したコードの比較）の結果 緩和時間 16.2 ms 16.2ms 開発したコード　　　衝突オペレータのみ　 衝突オペレータ テスト計算2 横長断面上の磁束 の磁気面上 ○使用磁場 ○初期ピッチ角 ○初期位置 ○モンテカルロ粒子数 1,300 一様に分布 ○粒子の初期エネルギー 200 keV ○粒子 陽子 ○back ground(軽水素プラズマ) 温度 1 keV 密度 1020 m-3 ピッチ角-分布関数 速さ-分布関数 磁束-分布関数 磁束(位置), 速さ, ピッチ角に対する分布関数 ?初期の粒子位置 でピーク ?プラズマ中心部ならびにLCFS近傍まで分布関数が到達 ? で分布関数がピークを持つ ? で分布関数が０になる ? の粒子を緩和と仮定 ? と において, 分布関数が減少 ? 損失領域が存在 磁束-分布関数 速さ-分布関数 ピッチ角-分布関数 中性粒子の密度 中を1個の荷電粒子が速度 で飛ぶとき Re-entering 粒子 ? 中性粒