천체 플라즈마 동역학을 위한 하이브리드‑키네틱 PIC 코드 ‘페가수스’

초록

**
페가수스는 하이브리드‑키네틱 입자‑인‑셀(PIC) 방식을 채택한 천체 플라즈마 전용 시뮬레이션 코드이다. 에너지 보존형 입자 적분기와 2차 정확도의 3단계 예측‑예측‑보정 스키마, 제한 전송(constrained transport) 기반 무발산 자기장 유지, δf 저노이즈 스킴, 그리고 회전·전단을 위한 시어링‑시트와 궤도 이동(orbital advection) 기법을 통합한다. Athena와 유사한 모듈형 구조와 수천 코어 병렬화에 최적화된 설계로, 1‑3차원 테스트를 통해 정확도와 확장성을 검증하였다.

**

상세 분석

**
페가수스는 비상대론적, 준중성, 수소성 플라즈마를 대상으로 하는 하이브리드‑키네틱 모델을 구현한다. 전자는 질량이 0인 MHD‑유사 유체로 가정하고, 이온은 완전 입자‑입자(PIC) 방식으로 다룬다. 이 접근법은 전자와 이온의 스케일이 크게 분리된 천체 플라즈마(예: 뜨거운 흡입 흐름, 은하단 내부, 태양풍 등)에서 전자 플라즈마 진동이나 빛의 속도 제한을 제거함으로써 계산 비용을 크게 낮춘다.

코드의 핵심 알고리즘은 세 가지로 요약할 수 있다. 첫째, 에너지 보존형 입자 푸시를 위해 수정된 Boris 스킴을 사용한다. 여기서는 전기·자기장 값을 시간 단계 중간값(tⁿ⁺¹/₂)에서 평가하고, 속도‑자기장 상호작용을 Crank‑Nicolson 형태로 반암시적으로 처리한다. 가속도 항 a(v,B)=v×B−2Ω ẑ×v+σvₓ ŷ를 포함해 회전·전단 효과를 정확히 반영한다. 둘째, 자기장의 무발산성을 보장하기 위해 제한 전송(constrained transport) 방식을 채택한다. 이는 ∇·B=0 조건을 기계적 정밀도로 유지하면서, 알벤·휘슬러 파동을 모든 스케일에서 안정적으로 전파한다. 셋째, δf 스킴을 도입해 초기 분포함수와의 작은 편차만을 추적함으로써 입자 잡음(noise)을 크게 감소시킨다. 이는 특히 전단 흐름에서 지속적인 자유 에너지 공급이 존재할 때 “weight growth” 문제를 완화한다.

시어링‑시트(formalism)와 궤도 이동(orbital advection) 기법은 회전·전단이 강한 원반 물리학 문제를 효율적으로 다룰 수 있게 한다. 배경 전단 흐름 u₀(x)=−σ₀ x ŷ에 대한 입자와 자기장의 advection을 별도 연산으로 분리함으로써 경계 조건을 주기적으로 적용하고, 전단에 의한 y‑방향 스트레칭을 정확히 구현한다.

아키텍처 측면에서 페가수스는 Athena 코드와 동일한 모듈식 구조를 갖는다. 물리 모듈(입자 푸시, 전류 계산, 전기장/자기장 업데이트)과 수치 모듈(보간, 입자 deposit, 경계 처리)이 명확히 구분돼 사용자가 새로운 물리 모델을 추가하거나 기존 모듈을 교체하기 쉽다. 또한 MPI 기반 도메인 분할과 OpenMP 스레드 병렬화를 결합해 수천 코어 규모의 슈퍼컴퓨터에서도 효율적인 스케일링을 보인다.

검증 테스트는 1D 전자·이온 플라즈마 파동, 2D 알벤 파동 전파, 3D 전단 불안정성(MRI) 시뮬레이션 등으로 구성된다. 모든 테스트에서 2차 정확도와 에너지 보존이 확인되었으며, δf 스킴을 사용한 경우 입자 잡음이 최소화된 것이 눈에 띈다. 특히, 전단 흐름에서 whistler 파동이 불안정해지는 기존 2‑step 스키마와 달리, 페가수스의 3‑단계 예측‑예측‑보정 스키마는 안정적으로 파동을 전파한다.

전체적으로 페가수스는 천체 플라즈마 연구에 필요한 물리적 정확도와 수치적 안정성을 동시에 제공하면서, 대규모 병렬 컴퓨팅 환경에 최적화된 실용적인 도구로 자리매김한다.

**