[논문 리뷰] Gaseous Dynamical Friction: a Numerical Study of Extended Perturbers

B. R. Morton, Sadegh Khochfar외 1명|arXiv (Cornell University)|2021. 03. 29.

Fluid Dynamics Simulations and Interactions참고 문헌 5인용 수 5

한 줄 요약

이 연구는 GIZMO 코드를 사용한 고해상도 라그랑주 유체역학 시뮬레이션을 통해 기체의 천체적 마찰력을 조사하며, 현대 유체역학적 방법이 마하 수 0.7에서 2 사이에서 마찰력을 10%에서 50% 정도 과소평가하고 있음을 밝혀냈다. 특히 M=1 근처에서 그러한 경향이 두드러진다. 이 과소평가의 원인은 마하 원뿔의 앞선 가장자리에서 날카로운 밀도 프로파일이 수치적으로 퍼짐으로써 발생하며, 이는 위험대 구조를 왜곡시키고 운동량 이행을 감소시킨다. 이는 이러한 해법에 의존하는 천체론적 시뮬레이션의 정확성을 도전한다.

ABSTRACT

The process of momentum and energy transfer between a massive body and a background medium it is moving through is known as dynamical friction (DF). It is key to our understanding of many astrophysical systems. We present a series of high-resolution simulations of gaseous DF using Lagrangian meshless finite mass hydrodynamics solver, the moving-mesh MUSCL scheme, and the piecewise parabolic method (PPM) solver. We use a set of simulations of massive bodies, modelled as Plummer spheres, moving with Mach $0.2 \leq \mathcal{M} \leq 3$. We investigate at which radial distances from the perturber these solvers recover the linear point mass solution for gaseous DF. We analyse the drag force and the structure and time evolution of the wake. The different solvers agree closely. Numerical convergence is reached when the initial spatial resolution is $0.2r_s$, where $r_s$ is the softening scale of the Plummer sphere. We find that the wake structure and drag force are recovered, at the $5\%$ level, when compared beyond $4r_\mathrm{s}$. Our results predict that models using the standard linear point mass DF solution will overestimate the drag force on extended perturbers by as much as 25\%, for Mach$\sim$1. Finally, we consider DF in the context of galaxy clusters, where dark matter subhaloes move through circumgalactic media. We show that DF is typically in the linear regime for most subhaloes in hosting haloes $<10^{11}$ M$_{\odot}$ but non-linear in more massive host haloes.

연구 동기 및 목표

• 다양한 마하 수에서 현대 라그랑주 유체역학 해법이 기체의 천체적 마찰력을 얼마나 정확하게 시뮬레이션하는지 평가하기.
• 기체 매질을 통과하는 거대한 외란자에 대한 수치적 시뮬레이션과 해석적 예측 간의 마찰력 차이의 근본 원인을 규명하기.
• 수치적 마찰력 과소평가가 해상도 부족, 위험대 구조의 부정확성, 또는 유체역학 해법의 본질적 한계 때문인지를 평가하기.
• 이 마찰력 과소평가가 천체론적 시뮬레이션, 특히 암흑물질 하위구조 융합률에 미치는 영향을 규명하기.
• 해석적 해가 잘 정의되어 있고 물리적으로 관련성이 높기 때문에, 천체적 마찰력을 유체역학 코드의 표준 기준 테스트로 도입할 것을 주장하기.

제안 방법

• GIZMO 코드의 라그랑주 유체역학 해법(GSPH 및 G4PSPH)을 사용해 균일한 기체 매질을 통과하는 거대한 외란자를 모델링하기 위한 고해상도 시뮬레이션.
• 마하 수 범위를 0.2에서 3까지로 설정하고, 다양한 중력 연화 길이와 초기 조건을 사용해 결과의 탄력성 테스트.
• 천체적 마찰력의 해석적 예측은 선형 섭동 이론에서 유도되며, 특히 Ostriker(1999) 및 Just & Kegel(1990)의 해법을 기반으로 한다.
• 특히 마하 원뿔의 앞선 부분에서의 밀도 프로파일과 초과밀도 특징을 측정하여 위험대의 반경 구조 분석.
• 마찰력과 위험대 형태의 시간적 변화를 추적하여 수렴성과 장기적 거동 평가.
• 다양한 반경과 마하 수에서 수치적 마찰력과 해석적 예측을 비교하여 마찰력 과소평가의 크기를 정량화함.

실험 결과

연구 질문

• RQ1 현대 라그랑주 유체역학 해법이 기체 매질을 통과하는 거대한 외란자에 대해 해석적 마찰력 예측을 얼마나 잘 재현하는가?
• RQ2 왜 수치적 마찰력은 마하 수 범위 0.7 < M < 2에서 체계적으로 해석적 예측보다 낮게 나타나는가?
• RQ3 마하 원뿔의 밀도 프로파일이 수치적으로 퍼짐으로써 전체 운동량 이행과 마찰력 계산에 어떤 영향을 미치는가?
• RQ4 전체 마찰력의 몇 퍼센트가 원형 외부 위험대(r ≥ 4rs)에서 기인하는가? 이 비율은 시간이 지남에 따라 어떻게 변화하는가?
• RQ5 마찰력 과소평가가 천체론적 시뮬레이션에서 암흑물질 하위구조 융합률에 미치는 영향은 어느 정도인가?

주요 결과

• 마하 수 0.7에서 2 사이에서는 시뮬레이션 마찰력이 해석적 예측보다 10%에서 50%까지 체계적으로 낮게 나타나며, 특히 M=1 근처에서 가장 큰 과소평가가 발생한다.
• 마찰력 과소평가의 주요 원인은 마하 원뿔의 앞선 가장자리에서 날카로운 초과밀도 프로파일이 수치적으로 퍼짐으로써 발생하며, 이는 최대 밀도를 감소시키고 밀도를 앞으로 퍼뜨린다.
• 내부 영역인 4rs(A=0.1일 경우) 또는 8rs(A=1일 경우)에서 가장 심각한 구조적 불일치가 발생하며, 여기서 위험대의 날카로운 전이가 부드럽게 퍼져나간다.
• 원형 외부 위험대(r ≥ 4rs)는 수치적으로 잘 재현되며, 총 마찰력에 기여하는 비율은 시간이 지남에 따라 증가하지만, 여전히 내부 영역의 과소평가를 상쇄하기에는 부족하다.
• 마찰력 과소평가는 최소 150개의 교차 시간 동안 지속되며, GSPH 및 G4PSPH 해법, 다양한 초기 조건에 대해 모두 탄력적이다.
• 결과적으로, 이러한 해법을 사용하는 천체론적 시뮬레이션은 암흑물질 하위구조가 경험하는 천체적 마찰력을 과소평가할 가능성이 있으며, 이는 융합률을 과소평가하고 허브 질량 함수를 정확히 반영하지 못하게 한다.