[論文レビュー] Nonlinear excitation of low-n harmonics in reduced MHD simulations of edge-localized modes

I. Krebs, M. Hoelzl+2 more|arXiv (Cornell University)|May 16, 2013

Magnetic confinement fusion research参考文献 2被引用数 3

ひとこと要約

本研究では、JOREKコードを用いた低次元MHDシミュレーションを用いて、エッジ局在モード（ELMs）中の低n調波の非線形励起を調査する。非線形励起のメカニズムとして、2次非線形相互作用に基づく2次的結合モデルを提案し、主要な調波から部分的な低nモード（特にn = 1）へエネルギーが移行することを示している。このメカニズムにより、ELMsの初期非線形段階において、プラズマ境界で顕著に増幅し、空間的に局在化する低nモードの成長が説明できる。

ABSTRACT

Nonlinear simulations of the early ELM phase based on a typical type-I ELMy ASDEX Upgrade discharge have been carried out using the reduced MHD code JOREK. The analysis is focused on the evolution of the toroidal Fourier spectrum. It is found that during the nonlinear evolution, linearly subdominant low-n Fourier components, in particular the n = 1, grow to energies comparable with linearly dominant harmonics. A simple model is developed, based on the idea that energy is transferred among the toroidal harmonics via second order nonlinear interaction. The simple model reproduces and explains very well the early nonlinear evolution of the toroidal spectrum in the JOREK simulations. Furthermore, it is shown for the n = 1 harmonic, that its spatial structure changes significantly during the transition from linear to nonlinearly driven growth. The rigidly growing structure of the linearly barely unstable n = 1 reaches far into the plasma core. In contrast, the nonlinearly driven n = 1 has a rigidly growing structure localized at the plasma edge, where the dominant toroidal harmonics driving the n = 1 are maximal and in phase. The presented quadratic coupling model might explain the recent experimental observation of strong low-n components in magnetic measurements [Wenninger et al., Non-linear magnetic perturbations during edge localized modes in TCV dominated by low n mode components, submitted to Nuclear Fusion].

研究の動機と目的

• トカマクプラズマにおけるタイプI ELMsの初期段階における低nトロイダル調波の非線形励起を理解すること。
• 線形理論では予測できないが、実験的に観測された磁気測定における強い低n成分（特にn = 1）を説明すること。
• 線形的にわずかに不安定な状態から、非線形的に駆動されるn = 1モードの成長への遷移を調査すること。
• ELMs中のトロイダル調波間のエネルギー移行の本質的物理を捉える最小限のモデルを構築すること。

提案手法

• 典型的なタイプI ELMy ASDEX Upgrade放電を対象として、JOREKコードを用いた非線形低次元MHDシミュレーションを実施する。
• トロイダルフーリエスペクトルの時間発展を解析し、異なる調波成分間でのエネルギー分布を追跡する。
• トロイダル調波間の2次非線形相互作用に基づく2次的結合モデルを構築する。
• モデルの予測をシミュレーション結果と比較し、エネルギー移行メカニズムの妥当性を検証する。
• n = 1モードの空間的構造の、線形成長段階から非線形成長段階への遷移を検証する。
• n = 1モードの非線形励起に最も寄与する結合の最大値を示す位置および位相が一致する駆動調波を同定する。

実験結果

リサーチクエスチョン

• RQ1 部分的な低n調波（特にn = 1）は、線形的に部分的であるにもかかわらず、ELMsの初期非線形段階においてどのように増幅するのか？
• RQ2 主要な調波から低n成分（n = 1など）へ顕著なエネルギー移行を可能にする物理的メカニズムは何か？
• RQ3 n = 1モードの空間的構造は、線形領域から非線形領域に移行する際、どのように変化するのか？
• RQ4 単純な2次的結合モデルは、観測されたトロイダルフーリエスペクトルの非線形的時間発展を再現できるか？
• RQ5 プラズマ内のどこでn = 1モードの非線形励起が最も効果的に起こり、その局在化を決定づける要因は何か？

主な発見

• n = 1調波は、線形的に支配的でないにもかかわらず、非線形段階において、線形的に支配的な調波と同等のエネルギーにまで増幅する。
• 2次的結合モデルは、トロイダルスペクトルの初期非線形的発展をうまく再現し、2次非線形相互作用によるエネルギー移行を確認した。
• n = 1モードの空間的構造は、線形成長段階ではプラズマコアにまで広がる剛体的増幅プロファイルであったのに対し、非線形成長段階ではプラズマ境界に局在化する構造に変化する。
• 非線形的に駆動されるn = 1モードは、支配的駆動調波が最大かつ位相が一致するプラズマ境界に局在化しており、これが効率的なエネルギー移行を可能にしている。
• 本モデルは、最近の実験的観測で確認されたELMs中の磁気測定における強い低n成分の物理的説明を提供している。
• 空間的構造の変化は、非線形結合ダイナミクスがELMsの最終的なモード構造を規定する上で極めて重要であることを示している。