전방위 예측이 가능한 전자기학 전자기장 전산 시뮬레이션: TCV 삼각형 효과 검증

초록

이 논문은 Gkeyll 코드의 적응형 소싱 알고리즘을 이용해 자기장 형상, 가열 전력, 입자 총량만을 입력으로 하여 토카막 가장자리와 스크랩‑오프 레이어의 전자기학 전자기장 난류를 전자기학 전자기장 전산적으로 풀‑f 시뮬레이션한다. TCV 방전 #65125와 #65130에 대해 시뮬레이션 결과를 실험 데이터와 비교했으며, 특히 음의 삼각형(negative triangularity, NT) 형상이 E×B 전단을 약 20% 증가시켜 난류 감소와 열 방출 분포 변화를 일으키는 메커니즘을 제시한다.

상세 분석

본 연구는 전통적인 경계조건을 필요로 하는 전자기학 전자기장 모델과 달리, 입자와 에너지 소스를 자동으로 조절하는 적응형 소싱 기법을 도입함으로써 “입자 총량 보존”이라는 물리적 제약만을 만족한다. Gkeyll의 보존형 불연속 갈레킨(DG) 스킴은 격자 해상도와 무관하게 입자와 에너지 보존을 보장하고, 다중 GPU 병렬화로 수백 GPU·시간 내에 수밀리초 규모의 시뮬레이션을 수행한다.

핵심 물리 모델은 장파장(long‑wavelength) 전자기학 전자기장 방정식이며, k⊥ρi≪1 가정 하에 FLR 효과를 전위 방정식의 이온 편극항에만 포함한다. 충돌 연산자는 다종 Lenard‑Bernstein‑Dougherty(LBD) 형태를 사용해 공간·시간에 따라 변하는 충돌 주파수를 계산한다. 전위는 전자기학 전자기장 준중성식(3)을 통해 구해지며, Boussinesq 근사를 적용해 전위 라플라시안 항에 이온 편극을 포함한다.

기하학은 Miller 파라미터화를 기반으로 하며, 삼각형(δ)과 연장(κ) 파라미터를 자유롭게 조절한다. 특히 NT와 PT(positive triangularity) 두 경우를 동일한 전력·입자 조건에서 비교함으로써 형상 변화가 전자와 이온의 트랩 입자 전선(precession drift)에 미치는 영향을 정량화한다. 결과는 NT가 트랩 전자 모드(TEM)를 억제하고, 전자 열전달을 감소시키며, 동시에 E×B 전단을 증가시켜 난류 억제에 기여한다는 점을 보여준다.

시뮬레이션은 전자와 이온을 모두 동역학적으로 취급하고, 폐쇄 및 개방 필드 라인 영역을 하나의 도메인에 포함한다. 폐쇄 영역에서는 twist‑and‑shift 경계조건을, 개방 영역(리미터)에서는 전도성 셰일 경계조건을 적용해 전위 연속성을 유지한다. 적응형 소싱은 입자 손실(방출)과 에너지 손실을 실시간으로 모니터링하고, 소스 강도를 조절해 목표 입자 총량을 유지한다. 이는 실험에서 가스 주입·벽 조건을 조절하는 것과 동일한 물리적 효과를 수치적으로 구현한 것이다.

검증 결과는 Thomson scattering과 Langmuir probe 측정과 비교했을 때, 전자 온도·밀도 프로파일이 10~15% 이내의 오차로 재현되었으며, 블롭(Blob) 전송, 전위 구조, E×B 전단 프로파일 등 주요 현상이 실험과 일치한다. NT 경우 시뮬레이션은 전단이 약 20% 증가하고, SOL에서의 열 플럭스가 더 넓은 영역에 분산되는 경향을 보였다. 이는 NT가 플라즈마‑벽 상호작용을 완화하고, 향후 리액터 설계에서 열 부하 관리에 유리할 수 있음을 시사한다.

한계점으로는 전자기학 전자기장 전자기장(전기적) 근사, 단순화된 중성 입자 재순환 모델, 그리고 FLR 효과의 제한적 포함이 있다. 향후 연구에서는 전자기학 전자기장(전기·자기) 모드, 보다 정교한 중성 입자 재순환 및 부착 모델, 그리고 고전압 셰일 전위 조건을 추가함으로써 정확도를 높일 수 있다.

전반적으로 이 논문은 “입력 파라미터 최소화 → 전자기학 전자기장 전산 → 실험 검증”이라는 워크플로우를 성공적으로 구현했으며, 향후 차세대 토카막 및 DEMO 설계 단계에서 고신뢰성 예측 도구로 활용될 가능성을 제시한다.