数値計算による中性粒子枯渇現象の密度限界・密度振動への影響の解明

通过数值计算阐明中性粒子耗尽现象对密度极限和密度振荡的影响

• 批准号: 22K14020
• 负责人: 諌山 翔伍
• 金额: $ 3万
• 依托单位: Kyushu University
• 依托单位国家: 日本
• 项目类别: Grant-in-Aid for Early-Career Scientists
• 财政年份: 2022
• 资助国家: 日本
• 起止时间: 2022-04-01 至 2025-03-31
• 项目状态: 未结题
• 来源: https://kaken.nii.ac.jp/en/grant/KAKENHI-PROJECT-22K14020/
• 关键词: ヘリコンプラズマ; 中性粒子密度枯渇; プラズマ密度限界; 3流体数値シミュレーション; 電気推進

中性粒子枯渇現象の基本的な特性を理解するため、空間2次元の数値計算モデルの構築と数値計算コードの高速化を行うことを本年度の目標とした。本流体コードは電子流体について、電子慣性を無視（ドリフト近似）した運動方程式、電子のエネルギー方程式、ポアソン方程式をセルフコンシステントに取り入れている。そのため陽解法の場合、時間ステップはプラズマ周波数と衝突周波数の比で決まる誘電緩和時間（電子の運動が電場をスクリーニングする特徴的なタイムスケール）よりも小さくとる必要がある。特に10^19 m^-3の密度を超えるヘリコンプラズマ生成の数値計算では、誘電緩和時間が光速で決まるCFL条件よりも十分に小さくなる。これまでの先行研究では、電子とイオンの両方の運動方程式をドリフト近似する手法やポアソン方程式を半解法で解く方法が提案されている。しかしながら、中性粒子の運動はイオンの運動によって決まるため、我々の問題の場合には、イオンの慣性項を無視するドリフト近似は適当でない。また、今回我々はポアソン方程式を半解法で解く方法を試したが、密度が10^19 m^-3程度になると数値不安定を引き起こす結果となった。したがって、現状本モデルでは、電子流体については完全な陰解法を用いており、これまでにMPI並列(30cpu数を使用)により約12倍の高速化を実現している。ヘリコン放電が定常状態(~100 μs)になるまでの所要時間はおおよそ30日間である。現在は本モデルを用いて中性粒子枯渇現象の解析を進めている。また、今後の機械学習とGPUを活用したコードの高速化を念頭において、ニューラルネットワークにより流体計算を行う手法も確立した。次年度については、中性粒子枯渇現象の解析を引き続き進め、更なるコードの高速化を図る。

Phenomenon of neutral particles is 渇 の な characteristics of basic を understand す る た め, 2 dimensional space の the numerical computing モ デ ル の build と the numerical calculation コ ー ド の high speed line を う こ と を の target this year と し た. This fluid コ ー ド は electronic fluid に つ い て, electron inertia を ignored (ド リ フ ト approximation) し た movement equation, electronic の エ ネ ル ギ ー equations, ポ ア ソ ン equation を セ ル フ コ ン シ ス テ ン ト に take り れ て い る. そ の た め Yang solution の occasion, time ス テ ッ プ は プ ラ ズ マ cycle for と conflict cycle for の than で definitely ま る induced electric time (electronic の movement が electric field を ス ク リ ー ニ ン グ す る of 徴 な タ イ ム ス ケ ー ル) よ り も small さ く と る necessary が あ る. Trevor に 10 19 m ^ - ^ 3 の density を super え る ヘ リ コ ン プ ラ ズ マ generated の the numerical calculation で は, induced electric time で が light definitely ま る CFL condition よ り も very small に さ く な る. こ れ ま で の leading research で は, electronic と イ オ ン の struck party の movement equation を ド リ フ ト approximate す る gimmick や ポ ア ソ ン equation を half solution く で solution methods proposed が さ れ て い る. し か し な が ら, neutral particles moving の は イ オ ン の movement に よ っ て definitely ま る た め, I 々 の problem の occasions に は, イ オ ン の inertia item を ignore す る ド リ フ ト approximate は appropriate で な い. ま た, now back to my 々 は ポ ア ソ ン equation を half solution く で solution methods try を し た が, density が 10 19 m ^ - ^ 3 degree に な る と the numerical instability を lead き up こ す results と な っ た. し た が っ て, the present situation the モ デ ル で は, electronic fluid に つ い て は completely を な Yin method using い て お り, こ れ ま で に MPI parallel (30 CPU number を use) に よ り about 12 times the high speed の を be presently し て い る. Youdaoplaceholder0 <s:1> discharge が stationary state (~100 μs)になるまで <s:1> time required おおよそ おおよそ30 days である. Now the モデ モデ を を uses the 渇 phenomenon of neutral particles in モデ て to analyze を into めて る る る. ま た, future の rote learning と GPU を use し た コ ー ド thought high speed の を に お い て, ニ ュ ー ラ ル ネ ッ ト ワ ー ク に よ り fluid line calculation を う gimmick も establish し た. Annual に つ い て は withered 渇 phenomenon, neutral particles の parsing を lead き 続 き into め, more な る コ ー ド high speed の を 図 る.

项目成果

The mechanism of the helicon plasma production from low-density to high-density regime

螺旋等离子体从低密度到高密度状态的产生机制

• 发表时间: 2022
• 作者: Yohsuke T Fukai;Kyogo Kawaguchi;Yohsuke T. Fukai;諌山翔伍;諌山翔伍
• 通讯作者: 諌山翔伍

Dynamic profile formation of neutral depletion in helicon plasma

螺旋等离子体中性耗尽的动态剖面形成

• 发表时间: 2023
• 作者: Yohsuke T Fukai;Kyogo Kawaguchi;Yohsuke T. Fukai;諌山翔伍
• 通讯作者: 諌山翔伍