PELOTON을 이용한 펠릿 로켓 가속 모델링 및 실험 검증

초록

본 논문은 3차원 라그랑지안 입자 코드 PELOTON에 펠릿 로켓 가속을 위한 직접 수치 모델을 구현하고, JET에서 수행된 파편화 펠릿 주입(SPI) 실험과 비교 검증하였다. 고자기장 측면(HFS)과 저자기장 측면(LFS) 사이의 전자 흐름 비대칭과 구름 전하 불균형이 압력 차이를 유발해 펠릿 파편을 LFS 방향으로 가속한다는 메커니즘을 정량화하였다. 순수 중수소와 0.5 % 네온을 함유한 복합 펠릿을 비교한 결과, 네온 함량이 로켓 가속을 억제함을 확인하였다.

상세 분석

본 연구는 펠릿이 고온 플라즈마에 주입될 때 형성되는 비대칭적인 승화 구름이 grad‑B 드리프트에 의해 HFS와 LFS 사이에 전자 입사 강도가 달라지는 현상을 정밀히 모델링한다. 저자기장 측면에서는 전자 차폐가 강해 전자 흐름이 억제되지만, 고자기장 측면에서는 차폐가 약해 전자 입사가 증가하고, 이로 인해 구름 내부 압력이 비대칭적으로 상승한다. 압력 차이 ΔP = P_HFS − P_LFS는 펠릿 표면에 작용하는 로켓 가속력 a = ΔP·r_p²/(π m_p) 로 전환된다. 논문은 이 압력 차이를 직접 계산하기 위해 전자 열 플럭스의 비균일성을 기술하는 식(9)–(18)을 도입하고, 전자 차폐 전위 Φ를 비선형 방정식(16)으로 해결한다. 또한, 저자기장 Reynolds 수가 작은 MHD 근사(δB/B ≪ 1)를 적용해 구름의 전자·이온 흐름, 압력, 온도 변화를 라그랑지안 형태의 연속·운동량·에너지 방정식(4)–(7)으로 기술한다. 여기서 전류 J는 σ⊥ u × B 로 계산되며, 전도도 σ⊥는 구름 온도와 이온화 상태에 따라 동적으로 변한다. 중요한 점은 grad‑B 드리프트 가속도 식(8)에서 구름 길이 L과 자기장 곡률 R을 통합해 전반적인 구름 이동을 정량화한다는 것이다.

플라즈마 배경 상태는 새로운 냉각 모델을 통해 구름이 차지하는 플럭스 체적에 물질을 분배하고, 이온화 손실, 옴 가열, 이온·전자 에너지 교환을 모두 포함한다. 이 모델은 JOREK와 INDEX 시뮬레이션 결과와 비교해 온도·밀도 프로파일을 재현한다.

핵심 물리적 인사이트는 다음과 같다. 첫째, 로켓 가속은 압력 차이에 의해 지배되며, 전자 열 플럭스의 비대칭이 그 원천이다. 둘째, 구름 내부의 전하 불균형(전기적 차폐)과 grad‑B 드리프트가 결합해 HFS 쪽에 더 높은 전자 밀도와 압력을 만든다. 셋째, 파편 간 구름 겹침(cloud overlap)은 압력 차이를 상쇄하거나 증폭시켜 가속도에 큰 변동을 초래한다. 넷째, 배경 플라즈마의 온도·밀도 구배가 로켓 가속에 보조적인 역할을 하며, 높은 온도·밀도 구역에서는 가속이 강화된다. 마지막으로, 소량의 네온(0.5 %)을 첨가하면 전자 차폐가 강화되어 전자 열 플럭스가 감소하고, 결과적으로 압력 차이가 감소해 로켓 가속이 억제된다.

수치 실험에서는 다양한 파편 배치(선형, 원형, 무작위)와 구름 겹침 정도를 변형해 가속도와 궤적을 조사하였다. 결과는 JET 실험에서 측정된 파편 궤적과 정량적으로 일치했으며, 특히 네온 함유 펠릿이 순수 중수소 펠릿보다 궤적 편향이 작다는 점을 재현했다. 향후 ITER 플라즈마 조건에 대한 예측과 스케일링 법칙 도출을 목표로 하고 있다.