다중이온 기가스키넥틱 시뮬레이션으로 밝힌 DTHe 연소 플라즈마의 입자와 에너지 균형

초록

본 연구는 ITER와 유사한 D‑T‑He 플라즈마에서 이온 온도 구배(ITG)와 트랩 전자 모드(TEM) 구동 난류를 다중종 기가스키넥틱 Vlasov 시뮬레이션으로 조사한다. D와 T 연료 이온의 입자 플럭스가 D‑T 비율 및 He‑ash 축적에 따라 크게 불균형함을 확인하고, 이를 바탕으로 Reiter의 “steady burning” 조건(He‑ash 보존시간 < α · 에너지 보존시간)을 만족하는 최적 프로파일 영역을 최초로 제시한다. 또한, 제논 흐름(zonal flow)과 비열적 He‑ash가 최적 영역에 미치는 영향을 분석한다.

상세 분석

이 논문은 전통적인 0‑차원 전력 균형 모델이 간과하는 연료와 He‑ash의 입자 수송을, 실제 질량을 갖는 전자와 상호종 충돌을 포함한 다중종 전자기 기가스키넥틱 Vlasov 방정식(GKV)으로 직접 계산한다. 선형 분석에서는 D‑T‑He 혼합 플라즈마가 ITG와 TEM 두 모드가 동시에 불안정해짐을 보여주며, He‑ash 10 % 농도에서도 성장률과 주파수는 크게 변하지 않는다. 비선형 시뮬레이션에서는 전형적인 gyro‑Bohm 스케일의 에너지 플럭스(Q)와 입자 플럭스(Γ)가 정착되지만, D와 T의 입자 플럭스는 He‑ash 존재 시 방향이 뒤바뀌는 비대칭성을 보인다. 특히, Γ_D와 Γ_T는 각각 내부(π ≈ 3.14)와 외부(0) 영역에서 반대 부호를 가지며, 단일 유효 질량(A_eff) 모델에서는 재현되지 않는다. 이는 다중종 효과가 입자 보존식 ∑ Γ_i = Γ_e 를 만족하도록 개별 이온 플럭스를 자유롭게 조정하게 하기 때문이다.

He‑ash 농도를 0 %, 5 %, 10 %로 변화시키면서 D‑T 비율을 스캔한 결과, 에너지 플럭스(Q_D, Q_T)는 D‑T 비율이 0.5(동등 혼합)일 때 거의 동일하게 수렴하지만, 입자 플럭스(Γ_D, Γ_T)는 He‑ash이 존재하면 0.5에서 크게 불균형해진다. 즉, 연료 이온이 내부로 피닝(pinch)되는 반면, He‑ash은 외부로 배출되는 흐름이 강화된다. 이러한 불균형은 Reiter 조건을 만족시키는 핵심 파라미터인 τ_He/τ_E = ρ_τ에 직접적인 영향을 미친다.

Reiter 조건을 τ_He < α · τ_E (α≈7–15) 형태로 재정의하고, τ_He와 τ_E를 각각 He‑ash 확산계수 D_eff와 이온 열전도계수 χ_eff 로 표현한다. 시뮬레이션 결과는 Γ_He > 0, Γ_D, Γ_T < 0 (연료 내부 피닝)와 동시에 Q_i/α > η_i T_i (n_i/n_He) Γ_He 를 만족하는 프로파일이 존재함을 보여준다. 특히, D‑T 비율이 0.4–0.6 사이이고 He‑ash 농도가 5 % 이하일 때, η_i와 온도 구배가 적절히 조정되면 α ≈ 10 ~ 12 범위 내에서 조건을 만족한다.

제논 흐름의 억제와 비열적 He‑ash(고에너지 알파 입자)의 삽입을 추가 실험했을 때, 제논 흐름이 강해지면 난류 전파가 감소해 입자 피닝이 약화되고, 비열적 He‑ash는 전자 온도 구배를 완화시켜 TEM을 억제함으로써 전체 난류 강도가 낮아진다. 따라서 최적 프로파일을 설계할 때는 제논 흐름을 적절히 활용하고, 비열적 He‑ash 배출을 조절하는 것이 중요하다.

전반적으로, 다중종 기가스키넥틱 시뮬레이션은 연료 피닝과 He‑ash 배출 사이의 미세한 균형을 정량적으로 제시하며, 기존 단일종 모델이 놓친 입자 플럭스 비대칭성을 명확히 밝힌다. 이는 ITER와 차세대 장치에서 연료 주입 전략 및 He‑ash 배출 시스템 설계에 직접적인 가이드라인을 제공한다.