다중유체 비이상성 MHD 시뮬레이션을 위한 PIERNIK 코드 확장

초록

PIERNIK 코드는 Relaxing‑TVD 스키마를 기반으로 한 다중유체 격자 MHD 코드이며, 먼지와 확산성 우주선(CR) 가스를 독립적인 유체로 추가하고, 전단 상자 경계조건, 자체 중력, 오믹 저항성 등을 구현하였다. 본 논문에서는 특히 자기장에 의해 확산되는 CR 전파를 다루기 위한 비이상성 확산 모듈을 소개한다.

상세 분석

PIERNIK 코드는 Jin & Xin(1995)의 Relaxing TVD(RTV) 스킴을 기반으로, Trac & Pen(2003)과 Pen et al.(2003)의 구현을 확장한 형태이다. 기존 RTV는 고차 정확도와 충격파 포착 능력이 뛰어나지만, 전통적인 MHD에서는 단일 유체만을 다루었다. PIERNIK은 이 한계를 넘어, 동적으로 독립적이면서도 상호작용하는 다중 유체를 도입한다. 구체적으로는 (1) 먼지 유체, (2) 확산성 우주선(CR) 가스 유체를 추가했으며, 두 유체 모두 질량·운동량·에너지 방정식을 별도로 풀면서 전자기장과 중력에 의해 상호 결합한다. 특히 CR 유체는 유체 근사(fluid approximation)를 사용해 압력 텐서와 확산 계수를 정의하고, 비등방성 확산(자기장에 평행한 방향으로는 빠르고, 수직 방향으로는 느린 확산)을 구현한다. 이를 위해 전도 텐서를 자기장 방향에 따라 가중치 부여하는 방식으로 구성했으며, 시간 적분은 반명시적(implicit‑explicit) 스키마를 채택해 강한 확산 제한에도 안정성을 유지한다.

코드 구조는 모듈형 설계로, 새로운 유체를 추가하려면 인터페이스 정의와 초기화 루틴만 구현하면 된다. 이는 연구자들이 별, 원시성운, 은하 디스크 등 다양한 천체 물리 현상에 맞춰 맞춤형 유체를 손쉽게 삽입할 수 있게 한다. 또한, 전단 상자(shearing‑box) 경계조건을 지원해 원반 물리학 시뮬레이션에 필수적인 코리올리와 전단 흐름을 정확히 재현한다. 자체 중력 모듈은 FFT 기반 Poisson 솔버를 이용해 전역 중력 포텐셜을 계산하고, 오믹 저항성 모듈은 전기 전도도와 저항성 텐서를 통해 비이상성 마그네토수소 전도 현상을 모사한다.

병렬화 측면에서는 MPI를 이용한 도메인 분할을 적용했으며, 각 프로세스는 로컬 격자 블록을 담당한다. 통신은 주로 경계면 데이터 교환과 전역 중력 포텐셜 계산 시 FFT 변환을 위한 All‑to‑All 방식으로 이루어진다. 이러한 설계는 3차원 대규모 시뮬레이션에서도 높은 스케일링 효율을 보장한다. 특히 CR 확산 모듈은 비등방성 확산 계수를 공간적으로 가변시킬 수 있어, 관측적으로 추정되는 은하계 ISM의 복잡한 전도 구조를 재현한다. 전체적으로 PIERNIK은 기존 MHD 코드 대비 물리적 자유도가 크게 확대되었으며, 비이상성 현상을 포함한 다중 유체 상호작용을 정밀하게 탐구할 수 있는 강력한 도구가 된다.