적응형 라그랑지안 메쉬리스 MHD 코드 Phurbas 알고리즘 소개

초록

**
Phurbas는 이동 최소제곱(Moving Least Squares) 보간을 이용해 입자 기반으로 이상적인 자기유체역학(MHD) 방정식을 풀며, 라그랑지안 입자 이동, 자동 입자 생성·소멸, 2차 예측‑보정 시간 적분, 인공 점성·확산을 통해 안정성을 확보한다.

**

상세 분석

**
본 논문은 기존 고정격자·SPH 방식이 갖는 시간‑스텝 제한, 비갈릴리안 수치 확산, 입자 질량·해상도 비균일성 문제를 극복하기 위해 완전한 라그랑지안 메쉬리스 프레임워크를 설계하였다. 핵심은 주변 입자들의 값으로 3차 다항식을 최소제곱 피팅해 현 위치에서 함수값과 1·2차 공간미분을 얻는 Moving Least Squares(MLS) 보간이다. MLS는 SPH와 달리 영점 일관성(zero‑order consistency)을 보장하므로 상수함수와 그 미분을 정확히 재현한다. 이를 통해 압력·자기장·속도 등 모든 물리량을 고차 정확도로 평가할 수 있다.

시간 적분은 2차 예측‑보정(predictor‑corrector) 스키마를 사용한다. 라그랑지안 형태의 미분 연산을 그대로 적용하므로 입자 자체가 유체를 따라 이동하고, CFL 조건이 신호속도에만 의존하게 된다. 이는 고속 전이 흐름이나 원반과 같이 큰 전이속도가 존재하는 문제에서 큰 시간‑스텝 이득을 제공한다.

적응성은 두 단계로 구현된다. 첫째, 목표 해상도(공간별 입자 밀도)를 사전에 정의하고, 현재 입자 배치가 이를 만족하는지 검사한다. 둘째, ‘void‑clump’ 탐지를 통해 입자 간격이 과도히 크면 새로운 입자를 삽입하고, 과도히 촘촘하면 가장 불필요한 입자를 삭제한다. 삽입 시에는 주변 입자들의 MLS 보간값을 이용해 초기 물리량을 부여하고, 삭제 시에는 보존량을 최소화하도록 설계된 절차를 따른다.

수치적 안정성을 위해 동적 인공 점성·질량·열 확산 항을 도입한다. 점성 계수는 지역 흐름의 발산·전단에 따라 자동 조정되며, 충격파와 같은 급격한 변화를 억제한다. 또한, ∇·B=0 조건을 유지하기 위해 Dedner식 하이퍼볼릭‑패러볼릭 클리닝을 적용하고, 필요 시 추가적인 ‘divergence‑cleaning’ 절차를 수행한다.

알고리즘 전체 흐름은 (1) 입자 위치·물리량 초기화, (2) MLS 보간을 통한 공간미분 계산, (3) 인공 점성·확산 포함한 라그랑지안 MHD 방정식의 2차 시간미분 계산, (4) 예측‑보정 단계로 상태 업데이트, (5) 입자 적응(삽입·삭제) 및 시간‑스텝 재조정, (6) ∇·B 정리 로 구성된다.

이 설계는 고차 정확도와 라그랑지안 특성을 동시에 확보하면서도, 메쉬가 필요 없는 자유로운 입자 분포를 허용한다. 따라서 복잡한 기하학, 대규모 전이, 다중 스케일 현상을 다루는 천체물리·플라즈마 시뮬레이션에 적합한 기반을 제공한다.

**