고속 이온이 핵심 난류와 유리면에 미치는 영향

초록

글로벌 비선형 자이로킥틱 시뮬레이션을 이용해 고속 이온이 q=1 유리면 근처의 ITG 난류에 미치는 역할을 조사하였다. 고속 이온은 자기전단 구조를 강화하고, 희석 효과와 준공명 상호작용을 동시에 일으키며, 특히 전단 구조의 불안정 임계값을 낮춤으로써 난류를 최대 90% 억제한다. 또한 n=m=1 피시본 모드가 발생하면 비트‑구동 전단 구조가 추가적으로 난류를 감소시킨다.

상세 분석

본 연구는 GENE 코드의 글로벌 비선형 모드를 활용해 고속 이온(FI)이 핵심 플라즈마의 ITG(이온 온도 구배) 난류에 미치는 복합적인 메커니즘을 정량적으로 분석한다. 첫 번째 주요 메커니즘은 FI가 자기전단(E×B 전단) 구조를 강화하는 것이다. 전단 구조는 난류 자체의 비선형 자기상호작용에 의해 생성되며, q=1과 같은 저차 유리면에서 전단 구조가 특히 강하게 발달한다. FI가 존재하면 전단 구조의 불안정 임계값이 낮아져, 전단이 더 쉽게 성장하고, 결과적으로 난류 에너지 흐름을 억제한다. 두 번째 메커니즘은 FI에 의한 열·밀도 프로파일 희석이다. FI가 플라즈마에 추가되면 전하 중성 조건을 만족하기 위해 열 이온·전자 밀도가 감소하고, 이는 ITG 모드의 드라이버인 온도·밀도 구배를 약화시켜 선형 성장률을 20% 정도 낮춘다. 논문에서는 FI를 ‘희석 종’으로만 취급한 시뮬레이션과 실제 FI를 포함한 시뮬레이션을 비교함으로써 이 효과를 명확히 구분하였다. 세 번째 메커니즘은 FI와 ITG 모드 사이의 준공명(quasi‑resonant) 상호작용이다. FI의 드리프트 주파수 ω_r가 ITG 모드의 실주파수와 일치할 때 에너지 교환이 일어나며, 이는 FI 온도 T_f가 낮을수록(특히 T_f≈2) 성장률을 증가시키는 현상으로 나타난다. 논문은 Eq.(2)를 이용해 이 공명 조건을 정량화하고, 속도공간에서의 FI 열 플럭스가 T_f≈2에서 최대가 됨을 확인하였다.

또한, n=m=1 피시본 모드가 활성화되면 비트‑구동 전단 구조가 추가로 생성된다. 이 전단 구조는 피시본 모드와 난류 사이의 비선형 상호작용을 매개하며, 피시본이 강하게 성장하면 전단이 과도하게 평탄화되어 열·밀도 프로파일이 크게 변형될 수 있다. 그러나 피시본이 적당히 억제된 경우, 비트‑구동 전단이 기존 전단과 결합해 난류 억제 효과를 더욱 강화한다.

시뮬레이션 설정은 원형 플라즈마(반경 a=1, 대반경 R=10)와 두 종류의 q 프로파일(기준 q와 이동된 ˜q)을 사용했으며, 전자·이온·FI 3종을 포함한다. 격자 해상도는 (256,48,24,24,16)이며, 경계에서는 Krook 연산자를 이용해 플럭스와 프로파일 변형을 최소화하였다. 선형 스캔 결과는 n≈25에서 최대 성장률을 보이며, FI가 포함될 경우 성장률이 약 20% 감소한다. 전자·이온·FI 각각에 대한 전기·자기 열 플럭스와 E×B 전단률을 시간 평균하여 q=1 근처에서 90% 이상의 난류 억제를 확인하였다.

결과적으로, 고속 이온은 전단 구조 강화, 열·밀도 희석, 그리고 준공명 상호작용이라는 세 가지 경로를 통해 핵심 난류를 억제한다. 특히 전단 구조의 불안정 임계값 감소가 가장 효율적인 억제 메커니즘이며, 피시본 모드가 적절히 제어될 경우 추가적인 전단 생성으로 억제 효과가 극대화된다. 이러한 발견은 향후 알파 입자와 외부 가열에 의해 생성되는 고속 이온을 활용한 토카막 플라즈마의 전송 제어 전략에 중요한 시사점을 제공한다.