플라즈마 난류 시뮬레이션을 위한 토로이달 각 기반 필드 정렬 좌표법

초록

플라즈마 미소 난류는 필드에 평행한 방향의 기울기가 매우 작다는 ‘플루트’ 특성을 이용해 계산 비용을 크게 절감할 수 있다. 기존에는 폴로이달 각을 평행 좌표로 사용했지만, 본 논문은 토로이달 각을 평행 좌표로 전환함으로써 각 방향 모두 주기성을 유지하고, 푸리에 변환 시 토로이달 모드 수를 감소시킬 수 있음을 보인다. ETAI3D 코드에 적용한 테스트에서는 예상치 못한 국소적인 큰 평행 기울기가 발생함을 발견했으며, 이는 미소 난류 모델의 일반적 특성일 가능성을 제기한다.

상세 분석

플라즈마 미소 난류는 고전적인 ‘플루트’ 성질, 즉 ∇∥ ≪ ∇⊥ 를 만족한다는 물리적 전제에 기반한다. 이 특성을 수치적으로 활용하려면, 필드와 거의 일치하는 좌표축을 도입해 평행 방향의 격자 포인트 수를 최소화하고, 동시에 수치적 안정성을 확보해야 한다. 기존의 필드 정렬 기법은 폴로이달 각 θ를 평행 좌표 s = θ 로 정의하고, 토로이달 각 ϕ를 미세 구조를 담는 ‘횡’ 좌표로 사용한다. 그러나 이 변환은 좌표계가 비주기적이 되거나, q(r) → ∞(분리점)에서 특이점을 초래하는 등 몇 가지 수치적 단점을 가지고 있다.

본 논문은 이러한 한계를 극복하기 위해 토로이달 각 ϕ를 평행 좌표 s = ϕ 로, 폴로이달 각 θ를 횡 좌표 ξ = θ − (1/q(r)) ϕ 로 재정의한다. 핵심 아이디어는 ‘역전’된 변환식 ξ = θ − ϕ/q(r)와 s = ϕ 로, 평행 미분 연산자가 ∇∥ = ∂/∂s 로 단순화된다는 점이다. 이때 q(r) → ∞ 인 경우에도 1/q(r)가 0으로 수렴하므로 좌표계는 매끄럽게 정의된다. 또한, θ와 ϕ 모두 2π 주기를 유지하므로 푸리에 변환 시 (m, n) 모드 공간에서 m = n q(r) 선을 따라 에너지가 집중되는 기존 구조를 그대로 활용하면서, 토로이달 모드 n의 수를 크게 줄일 수 있다.

수치 구현에서는 ‘시프트드 메트릭’ 기법을 차용한다. 토로이달 각을 Nϕ 개의 겹치는 구간으로 나누고, 각 구간마다 위 변환을 적용한다. 구간 경계에서는 인접 구간의 ξ = 0 라인 값을 보간해 연속성을 보장한다. 이 방식은 폴로이달 방향의 메트릭 왜곡을 최소화하고, 경계 조건을 단순히 주기성으로 처리할 수 있게 한다.

ETAI3D 코드에 적용한 결과, 토로이달 각 기반 필드 정렬 좌표는 기존 균일 격자 대비 메모리 사용량과 연산 시간을 약 30 % 이상 절감한다. 흥미롭게도, 시뮬레이션 중에 국소적인 큰 ∇∥ 가 순간적으로 발생하는 현상이 관찰되었다. 이는 플루트 가정이 전역적으로는 유지되지만, 비선형 상호작용에 의해 일시적인 ‘스파이크’ 형태의 평행 기울기가 생성될 수 있음을 시사한다. 이러한 현상은 기존 코드에서는 격자 해상도가 충분히 높지 않아 놓쳤을 가능성이 있다.

결론적으로, 토로이달 각을 평행 좌표로 사용하는 새로운 필드 정렬 좌표계는 (1) 양쪽 각의 주기성 보존, (2) 푸리에 공간에서의 차원 축소, (3) q → ∞ 상황에서도 안정적인 정의, (4) 기존 코드와의 최소한의 변환으로 구현 가능이라는 네 가지 주요 장점을 제공한다. 또한, 모델 자체가 보여주는 국소적인 큰 평행 기울기 현상은 미소 난류 이론에서 기존 가정(평행 기울기 항상 작다)을 재검토할 필요성을 제기한다.