전체 볼륨 전자기 자이로키네틱 시뮬레이션을 위한 고차 필드 정렬 삼각형 유한요소법

초록

본 논문은 토카막 플라즈마의 핵심·경계·개방 필드 라인을 모두 포함하는 전자기 자이로키네틱 PIC 시뮬레이션을 위해, 원통 좌표계에서 C¹ 연속 삼각형 메쉬와 고차 필드 정렬 유한요소 기반의 피스와이즈(Field‑aligned) 방식을 제안한다. δf·full‑f 모두에 적용 가능하며, 혼합 변수‑풀백 스킴과 효율적인 Ampère 방정식 반복 해법을 도입해 수치 소멸을 억제한다. TRIMEG‑C1 코드에 구현한 뒤 TCV‑X21 구성을 대상으로 ITG·KBM 모드의 선형·비선형 결과를 검증, 핵심‑경계 영역을 정확히 포착함을 보였다.

상세 분석

이 연구는 전자기 효과가 β가 유의미한 토카막 플라즈마에서 미세불안정성과 난류를 정확히 기술하기 위한 첫 번째 전역 전자기 자이로키네틱 PIC 프레임워크를 제시한다. 기존의 전역 필드 정렬 좌표계는 토카막의 분리면(세퍼레이터) 근처에서 격자 왜곡과 특이점을 야기해 전체 볼륨 시뮬레이션에 한계를 보였다. 저자들은 이를 극복하기 위해 ‘피스와이즈 필드 정렬 유한요소(PFA‑FEM)’를 도입했는데, 이는 각 토로이달 서브도메인(ϕk 중심)마다 로컬 좌표 ηk를 정의하고, 그 좌표에 맞춰 고차 C¹ 삼각형 기반 함수를 배치한다. 이렇게 하면 각 서브도메인 내부에서는 필드 라인에 정확히 정렬된 격자를 유지하면서도, 전체 도메인에서는 전통적인 비정렬 메쉬(원통 좌표 R‑Z)와 동일한 정점 구조를 유지한다. 결과적으로 격자 왜곡을 최소화하고, 분리면을 포함한 전 영역을 연속적으로 계산할 수 있다.

수치적 안정성을 위해 혼합 변수(pullback) 스킴을 확장하였다. 전자기 포텐셜 δA∥를 symplectic(δAs∥)와 Hamiltonian(δAh∥) 부분으로 분해하고, u∥와 δf를 각각 변환함으로써 전자기 항에 의한 수치 소멸을 억제한다. 특히, 두 가지 pullback 옵션(이상 오옴 법을 만족하는 Scheme I와 ∂tδAs∥=0인 Scheme II)을 구현해 전자기·전기적 모드에 맞는 최적 선택을 가능하게 했다. Ampère 방정식은 skin‑depth 연산자를 포함한 연산자 (∇⊥²−∑s1/ds²)⁻¹ 로 구성된 반복 스킴으로 풀며, δAh∥의 초기값을 δAs∥와 결합해 빠른 수렴을 달성한다.

TRIMEG‑C1에 구현된 이 방법은 고차 B‑스플라인(ϕ 방향)과 18개의 C¹ 삼각형 기저함수를 결합해 3차원 변수 Y(R,ϕ,Z)를 전역적으로 표현한다. 피스와이즈 정렬 좌표 RFAk, ZFAk는 각 격자점에서 ϕK→ϕk까지의 필드 라인 적분으로 사전 계산되어 3차원 배열에 저장된다. 이렇게 사전 계산된 좌표는 입자 푸시와 필드 보간에 직접 사용돼 연산 비용을 크게 낮춘다.

검증 사례로 TCV‑X21 플라즈마를 선택했으며, 선형 전자기 ITG와 KBM 모드의 성장률·구조를 기존 전역 필드 정렬 코드와 비교해 일치함을 확인했다. 비선형 시뮬레이션에서는 전자기 플럭스 튜브와 전류 시트가 분리면 근처에서 자연스럽게 형성되고, 전자기 압력에 의한 β‑드리븐 비선형 포화 메커니즘이 재현되었다. 계산 효율 측면에서는 동일한 물리 해상도를 유지하면서도 전역 필드 정렬 좌표 대비 약 2‑3배 적은 메모리와 CPU 시간을 요구한다는 결과가 보고되었다. 이와 같이 PFA‑FEM은 복잡한 토카막 기하학(분리면·배리어·오픈 필드 라인)을 포함한 전자기 자이로키네틱 시뮬레이션에 필수적인 정확도와 효율성을 동시에 제공한다.