드리프트키네틱 시뮬레이션을 위한 보존 PSM 스키마 수치 제한기 테스트

초록

본 논문은 ITER 플라즈마 시뮬레이션에 사용되는 4차원 드리프트키네틱 모델을 대상으로, 보존형 반라그랑주(P​SM) 스키마에 적용 가능한 다양한 수치 제한기(ENT, UMEDA, OSL, SLS)를 제안·검증한다. Hermite 형식을 이용해 P​SM과 LAG 스키마를 통합하고, 제한기를 통해 발생하는 비물리적 진동을 억제하면서 질량·에너지 보존성을 유지한다. 저해상도와 고해상도 실험을 통해 각 제한기의 정확도·안정성을 비교 분석한다.

상세 분석

본 연구는 6차원 Vlasov‑Poisson 모델을 5차원 gyrokinetic, 최종적으로 4차원 drift‑kinetic 형태로 차원 축소한 뒤, 보존형 반라그랑주(Conservative Semi‑Lagrangian) 스키마인 PSM(Parabolic Splines Method)과 LAG(Lagrange polynomial) 방식을 Hermite 형식으로 재구성하였다. Hermite 형식은 원시 함수 G(x)의 누적값을 셀 경계에서 정확히 계산하고, 이를 고차 보간(큐빅 스플라인 또는 3차 라그랑주 다항식)으로 연장해 셀 내부 값을 재구성한다. 이 과정에서 고차 스키마는 무조건적인 4차 정확도를 제공하지만, 급격한 기울기나 비선형 구조에서 Gibbs 현상에 의한 비물리적 진동을 야기한다. 이를 완화하기 위해 네 가지 제한기—ENT(EnTrropic flux limiter), UMEDA, OSL(Oscillation Limiter), SLS(Slope Limited Splines)—를 도입하였다. 각 제한기는 (1) 고차 스키마를 저차 확산형 스키마로 전환하는 방법, (2) 진동이 예상되는 구간을 탐지하는 기준을 제공한다. 특히 OSL은 Hermite 계수를 이용해 국소 진동 지표를 계산하고, SLS는 전통적인 기울기 제한 기법을 스플라인 계수에 적용해 부드러운 전환을 보장한다. 수치 실험에서는 1‑D 상수 속도 전파 테스트와 4‑D 드리프트키네틱 모델을 사용해 제한기의 효과를 정량화하였다. 저해상도 시뮬레이션에서는 제한기가 없을 경우 질량 보존 오차가 급격히 증가하고, 전기 포텐셜에 비정상적 스파이크가 발생한다. 반면 OSL과 SLS는 이러한 오차를 12 차수 수준으로 감소시키며, 고해상도에서는 전체 에너지 보존과 전자기장 구조 재현에 큰 차이를 보이지 않는다. 또한, 제한기 적용 시 계산 비용은 약 1015% 증가했으나, 병렬화 효율은 크게 저하되지 않아 대규모 ITER 시뮬레이션에 실용적임을 확인하였다.