[논문 리뷰] Faranoff-Riley type I jet deceleration at density discontinuities "Relativistic hydrodynamics with realistic equation of state"

Z. Méliani, Rony Keppens외 1명|arXiv (Cornell University)|2008. 08. 18.

Astrophysics and Cosmic Phenomena참고 문헌 86인용 수 30

한 줄 요약

이 논문은 활성은하핵의 상대론적 제트가 외부 밀도 불연속성에 도달할 때 갑작스럽게 감속됨을 제안하며, 실제 상태방정식을 사용한 고해상도 상대론적 유체역학 시뮬레이션을 통해 이를 검증한다. 주요 발견은 강한 밀도 대비가 반사 충격을 통해 상당한 에너지 손실을 유도하고, 제트의 감속과 팬아로프-라이얼러 유형 I 형태로의 전이를 초래함을 보여주며, 하이브리드 형태의 전파원천에서 관측되는 비대칭 제트 행동을 설명한다.

ABSTRACT

The deceleration mechanisms for relativistic jets in active galactic nuclei remain an open question, and in this paper we propose a model which could explain sudden jet deceleration, invoking density discontinuities. This is particularly motivated by recent indications from HYMORS. Exploiting high resolution, numerical simulations, we demonstrate that for both high and low energy jets, always at high Lorentz factor, a transition to a higher density environment can cause a significant fraction of the directed jet energy to be lost on reflection. This can explain how one-sided jet deceleration and a transition to FR I type can occur in HYMORS, which start as FR II (and remain so on the other side). For that purpose, we implemented in the relativistic hydrodynamic grid-adaptive AMRVAC code, the Synge-type equation of state introduced in the general polytropic case by Meliani et al. (2004). We present results for 10 model computations, varying the inlet Lorentz factor from 10 to 20, including uniform or decreasing density profiles, and allowing for cylindrical versus conical jet models. As long as the jet propagates through uniform media, we find that the density contrast sets most of the propagation characteristics, fully consistent with previous modeling efforts. When the jet runs into a denser medium, we find a clear distinction in the decelaration of high energy jets depending on the encountered density jump. For fairly high density contrast, the jet becomes destabilised and compressed, decelerates strongly (up to subrelativistic speeds) and can form knots. We point out differences that are found between cylindrical and conical jet models, together with dynamical details like the Richtmyer-Meshkov instabilities developing at the original contact interface.

연구 동기 및 목표

• 하이브리드 형태의 전파원천에서 한쪽면이 FR I 유사 제트를 가지며 다른 쪽은 FR II 유사 제트를 가지는 상황에서 활성은하핵의 상대론적 제트 감속 메커니즘을 조사한다.
• 특히 한쪽 제트만 FR I 특성을 보일 때 상대론적 제트가 FR II에서 FR I 형태로 전이되는 과정에 대한 열린 질문을 해결한다.
• 제트가 높은 초기 로렌츠 인자에서라도 외부 밀도 불연속성에 도달할 경우 급격한 감속이 발생함을 보여준다.
• 에너지 전환과 충격 역학을 정확히 모의하기 위해 상대론적 유체역학 시뮬레이션에서 실제 다항 상태방정식의 사용을 검증한다.
• 변동하는 외부 밀도 프로파일을 가진 제트-매질 상호작용의 수치 모델링을 통해 하이브리드 원천에서 관측되는 비대칭 제트 형태의 물리적 기원을 설명한다.

제안 방법

• 메리아니 등(2004)의 싱게 기체 법칙 기반 일반 다항 상태방정식을 적응 메시 리피네먼트(AMR) 방식의 상대론적 유체역학 코드 AMRVAC에 구현하였다.
• 고해상도, 메시 적응 시뮬레이션을 사용하여 수백 개의 제트 빔 반지름에 걸쳐 제트의 전파를 모의하였으며, 서브그리드 리피네먼트를 통해 충격파와 접촉 불연속면을 해상도를 높여 해석하였다.
• 극도로 높은 로렌츠 인자(γ ≈ 100)에서의 리만 문제를 정밀하게 시험하였고, 수치 결과를 정확한 해석적 해와 비교하였다.
• 다양한 입구 로렌츠 인자(10–20)를 가진 고에너지 및 저에너지 제트 모델을 포함한 원통형 및 콘형 제트 기하구조를 시뮬레이션하였다.
• 특히 반사 충격 및 전방 충격 역학에 중점을 두어 상대론적 충격에서 정렬된 운동 에너지가 열에너지로 전환되는 과정을 추적하였다.
• 충격 효율성과 압축률을 결정하는 데 중요한 역할을 하는 효과적 다항 지수(Γ_eff)의 역할을 분석하였으며, 특히 상대론적(Γ=4/3)에서 고전적(Γ=5/3) 행동으로 전이되는 영역에서의 특성에 초점을 맞췄다.

실험 결과

연구 질문

• RQ1 외부 매질에서의 밀도 불연속성이 높은 초기 로렌츠 인자에서 출발한 상대론적 제트의 갑작스러운 감속을 유도할 수 있는가?
• RQ2 효과적 다항 지수는 제트-매질 상호작용 중 충격 효율성과 에너지 소산을 결정하는 데 어떤 역할을 하는가?
• RQ3 밀도 점프를 만났을 때 원통형과 콘형 제트 기하구조 간의 제트 형태와 감속 특성은 어떻게 다를까?
• RQ4 대칭적인 제트 방출에도 불구하고 일부 하이브리드 형태의 전파원천에서 한쪽면은 FR I 유사 제트를, 다른 쪽은 FR II 유사 제트를 보이는 이유는 무엇인가?
• RQ5 밀도 대비의 강도가 제트가 비상대론적 속도로 전이되고 응집된 덩어리(노드)를 형성하는 데 얼마나 기여하는가?

주요 결과

• 외부 밀도 대비가 강할 경우(n_ext / n_jet ≫ 1), 반사 충격에서의 에너지 반사로 인해 상당한 에너지 손실이 발생하며, 이로 인해 제트는 상대론적 속도에서 비상대론적 속도로 감속된다.
• 높은 로렌츠 인자 제트(γ = 10–20)의 경우, 임계 임계치를 초과하는 밀도 점프는 강한 압축, 불안정성, 리히트마이어-메시코프 유사 불안정성에 의해 발생하는 노드 형성으로 이어진다.
• 밀도 대비가 너무 약할 경우, 고에너지 제트는 최소한의 감속으로 인해 FR II 특성을 유지하며, FR I 형태로의 전이에 대한 임계값 행동을 보임을 시사한다.
• 압축 영역에서의 효과적 다항 지수 Γ_eff는 상대론적 영역(Γ=4/3)에서부터 약간 상대론적인 영역(Γ=1.48–1.5)으로 변화하며, 이는 충격 압축과 에너지 전환 효율성에 영향을 미친다.
• 원통형과 콘형 제트 모델 간의 동적 반응은 뚜렷한 차이를 보이며, 원통형 제트는 접촉면에서 더 강한 비대칭성과 이른 불안정성 성장이 관측된다.
• 극한의 경우(예: γ ≈ 100), 반사 충격과 접촉 불연속면 사이의 압축 영역이 ΔX₃(t) ≈ 0.00315t 정도로 축소되어 수치적 도전을 야기하며, 이는 적응 메시 리피네먼트로 해결되었다.