Improving divergence cleaning in cosmological SPMHD simulations

기술 요약: 우주론적 SPMHD 시뮬레이션에서의 발산 제거(Divergence Cleaning) 개선

문제 정의
우주론적 스무디드 입자 자기유체역학(SPMHD) 시뮬레이션에서, 솔레노이달 조건($\nabla \cdot \mathbf{B} = 0$)을 유지하는 것은 지속적인 수치적 과제이다. Gadget과 같은 코드에서 사용하는 표준 방식은 일반적으로 Powell 등(1999)의 "8-wave" 접근법에 의존한다. 이 방법은 발산 오차를 유체 흐름과 함께 이송시키지만, 이를 능동적으로 제거하지는 못한다. SPH 기반 코드에서 잔류 $\nabla \cdot \mathbf{B}$ 오차는 1에 가까운 차수에 도달할 수 있으며, 이는 비물리적인 자기장 증폭과, 자기 압력이 가스 압력을 초과할 때 운동 방정식에서 가짜 소스 항이 지배하게 되는 "인장 불안정성(tensile instability)"을 초래한다. Dedner 등(2002)이 제안한 쌍곡선/포물선형 발산 제거 방식은 해결책을 제공하지만, 그 원래의 정식화는 역동적이고 개별 타임스텝을 사용하는 환경, 특히 강한 밀도 변화가 있는 환경의 SPMHD에서 불안정한 것으로 판명되었으며, 종종 비물리적인 자기 에너지의 원천으로 작용한다.

방법론
저자들은 제약 조건이 있고 에너지를 보존하는 버전의 Dedner형 제거 알고리즘(Tricco & Price 2012; Tricco et al. 2016a)을 우주론적 SPMHD 코드인 OpenGadget3에 구현하였다. 주요 방법론적 적응 사항은 다음과 같다:

• 에너지 보존: 총 에너지 보존을 보장하도록 지배 방정식을 재정식화하였다. 여기에는 보조 스칼라 $\psi$와 제거 속도 $c_h$의 비율인 $\psi/c_h$로 제거 필드 변수를 진화시키고, 진화 방정식에 특정 소스 항 $-\frac{1}{2}\psi(\nabla \cdot \mathbf{v})$를 포함하는 과정이 포함된다. 이는 제거 필드의 포물선형 감쇠가 무한한 성장을 유발하지 않고 자기 에너지를 제거하도록 보장한다.
• 우주론적 적응: 방정식을 공변 좌표계와 개별 입자 타임스텝에 맞게 조정하였다. 제거 속도 $c_h$는 국소 빠른 자기음속($\sqrt{v_A^2 + c_s^2}$)으로 설정되어, 제거 과정이 표준 SPH 타임스텝 제약 내에서 본질적으로 Courant-안정성을 갖도록 하였다.
• 시뮬레이션 설정: 저자들은 균일한 자기장($10^{-14}$ G comoving)이 심어진 거대 은하단($M_{200c} \approx 10^{15} M_\odot$)에 대한 일련의 시뮬레이션을 수행하였다. 다섯 가지 뚜렷한 실행 케이스를 비교하였다:
  1. NIFTY-noclean: Powell 이송 방식만을 사용하는 참조 실행.
  2. NIFTY-clean: 엔트로피 역반작용을 포함한 전체 제약형 제거 방식이 적용된 기본 실행.
  3. NIFTY-clean2: 위와 동일하지만 오버클리닝 계수 $\zeta=2$ (국소 파동 속도의 두 배로 제거)를 적용한 실행.
  4. NIFTY-advect: 속도-발산 소스 항을 끈 상태의 제거 방식.
  5. NIFTY-heat: 소산된 제거 에너지가 엔트로피 방정식으로 피드백되지 않고 버려지는 제거 방식.

주요 결과

• 발산 오차 감소: 주요 결과는 Powell 방식에 비해 클러스터 부피 전체에서 상대적 발산 오차($|\nabla \cdot \mathbf{B}|/B$)가 2~3 자릿수(orders of magnitude) 감소했다는 것이다. 이 감소는 클러스터 중심부와 외곽 모두에서 균일하게 나타났다.
• 자기장 증폭: 자기장이 약하고 해상도가 낮은 클러스터 외곽 지역($r \approx 1\text{--}3 \, h^{-1}$ Mpc)에서, 능동적 제거 방식은 Powell-only 실행에 비해 자기장을 5~10배 증폭시켰다. 저자들은 이를 기존의 Powell 방식이 발산 오차로 인한 비물리적인 로런츠 힘을 제거함으로써, 이 지역의 소규모 다이너모(small-scale dynamo)를 억제했던 문제를 해결했기 때문이라고 설명한다.
• 중심부 안정성: 클러스터 중심부의 전역 밀도 구조와 자기장 기하 구조는 제거 방식과 비제거 방식 간에 거의 변하지 않았으며, 이는 제거 방식이 중심부 다이너모의 포화 상태를 인위적으로 변경하지 않음을 나타낸다.
• 열역학적 일관성: 엔트로피 역반작용(소산된 제거 에너지를 열로 다시 공급하는 것)의 포함이 열역학적 일관성을 위해 결정적임을 보여주었다. 이 역반작용이 없는 실행(NIFTY-heat)은 더 차가운 중심부의 단열 압축으로 인해 더 낮은 중심 엔트로피와 더 높은 포화 중심 자기장을 초래했다.
• 계산 비용: 이 구현은 측정 가능한 런타임 오버헤드를 발생시키지 않는다. 개선된 제거는 기존의 Courant 타임스텝 제약 내에서 사실상 "공짜"로 이루어진다.

의의 및 주장
본 논문은 정확한 발산 오차 제어가 은하단의 저밀도, 저해상도 외곽 지역에서 자기장 증폭을 포착하는 데 필수적이라고 주장한다. 저자들은 표준 Powell 방식의 통제되지 않은 발산 오차가 수치적 노이즈로 작고하여, 포화되지 않은 영역에서의 다이너모를 억제한다고 주장한다. 제약 조건이 있고 에너지를 보존하는 제거 방식을 구현함으로써, 시뮬레이션은 계산적 대가 없이 MHD 물리를 더 충실하게 구현한다.

저자들은 제약형 제거 방식이 다음을 수행한다고 결론짓는다:

1. 수치적 발산 오차를 극적으로 감소시킨다 (2~3 자릿수).
2. 해상도가 낮은 영역에서 자기장 증폭의 더 정확한 표현을 가능하게 한다.
3. 더 매끄럽고 일관된 자기장 구조를 생성한다.
4. 다양한 해상도에서 자기장 진화의 안정성을 높여, 중간 해상도에서도 물리적 충실도를 향 향상시킨다.

본 연구는 우주론적 SPMHD 응용 분야에서 자연스러운 제거 속도($\zeta=1$)를 사용하고 에너지 수지를 맞추기 위해 포물선형 에너지 역반작용을 활성화한 상태로 이 제거 방식을 기본값으로 채택할 것을 권장한다. 이 구현은 현재 OpenGadget3에서 공개적으로 사용할 수 있다.