[논문 리뷰] Selfgravitating disks in binary systems: an SPH approach -- I. Implementation of the code and reliability tests

Luis Diego Pinto, R. Capuzzo‐Dolcetta외 1명|arXiv (Cornell University)|2019. 06. 30.

Astrophysics and Star Formation Studies참고 문헌 59인용 수 3

한 줄 요약

이 논문은 이중성계 내 자가중력 작용을 받는 원행성원반을 시뮬레이션하기 위해 설계된 새로운 스무스드 파article 유체역학(SPH) 코드인 GaSPH를 제시한다. 이 코드는 안정성 및 성능 테스트를 철저히 통과하였으며, Sedov-Taylor 폭발파 검증을 포함하여 신뢰할 수 있는 수치적 성능을 입증하였다. 이중성에 의한 교란 하에서 원반 이심률과 근일점 각도의 진화에 대해 수렴성을 보이며, 고립되지 않은 환경에서의 원반 역학을 신뢰할 수 있게 모델링할 수 있다.

ABSTRACT

The study of the stability of massive gaseous disks around a star in a non-isolated context is not a trivial issue and becomes a more complicated task for disks hosted by binary systems. The role of self-gravity is thought to be significant, whenever the ratio of the disk to the star mass is non-negligible. To tackle these issues we implemented, tested and applied our own Smoothed Particle Hydrodynamics (SPH) algorithm. The code (named GaSPH) passed various quality tests and shows good performances, so to be reliably applied to the study of disks around stars accounting for self-gravity. This work aims to introduce and describe the algorithm, making some performance and stability tests. It constitutes the first part of a series of studies in which self-gravitating disks in binary systems will be let evolve in larger environments such as Open Clusters.

연구 동기 및 목표

• 자기중력 작용을 고려한 안정적인 SPH 코드를 개발하여 비고립 환경, 특히 이중성계 내 원행성원반을 시뮬레이션할 수 있도록 하는 것.
• 근접 입자 상호작용과 장거리 힘을 정확히 처리할 수 있는 트리 기반 중력 해법을 구현하고 테스트하는 것.
• 표준 유체역학 및 중력 기준 테스트를 사용하여 코드의 수치적 성능, 안정성 및 정확도를 검증하는 것.
• 자기중력 작용을 고려한 원반-은하계-이중성 상호작용을 모델링함으로써 향후 개방된 성단 내 원반 진화 연구를 가능하게 하는 것.
• 밀집된 은하 환경에서의 별-원반 피드백을 조사하기 위한 기초를 마련하는 것.

제안 방법

• 유체의 보간을 위해 입자형태의 SPH 수식을 사용하고, 입자 간 거리의 두 번째 차수 정확도를 확보하기 위해 삼차 스퍼린 커널 함수를 적용한다.
• 원거리 입자 군집은 단극자 및 이차극자 항을 사용하여 다중극 전개 방법으로 근사한다.
• 근접한 입자 간 중력 상호작용은 직접 계산하여 소프트닝 문제를 피하고 밀도가 높은 영역에서도 정확도를 유지한다.
• 점질량(은하계 또는 행성)의 운동은 고차수의 명시적 시간 적분기법을 사용하여 통합한다.
• 가스의 자기중력과 가스-은하계 상호작용을 자성적으로 처리하며, 복사 전달은 고려하지 않는다.
• 원반의 특성, 예를 들어 이심률과 근일점 각도는 이심률 벡터와 단위 질량당 운동량을 사용하여 계산한다.

실험 결과

연구 질문

• RQ1 트리 기반 중력 해법을 갖춘 SPH 코드는 이중성계 내 자가중력 원반의 진화를 얼마나 정확하게 모의할 수 있는가?
• RQ2 입자 해상도를 변화시킬 때 원반 이심률과 근일점 각도의 수치적 수렴 수준은 어느 정도인가?
• RQ3 Sedov-Taylor 폭발파와 같은 알려진 유체역학 기준을 높은 정밀도로 재현할 수 있는가?
• RQ4 이심률 e = 0.4인 이중성 동반성의 교란 하에서 원반 궤도 매개변수는 어떻게 변화하는가?
• RQ5 2D 모델과 비교했을 때 수직 중력과 3D 원반 구조는 원반 역학에 어떤 역할을 하는가?

주요 결과

• GaSPH 코드는 높은 정확도로 Sedov-Taylor 폭발파를 재현하여 강한 충격파 역학을 다룰 수 있음을 확인하였다.
• 초기 혼란 단계 이후 원반 이심률은 약 0.075의 평균값으로 안정화되며, 이는 이전 연구와 일치한다.
• 원반 이심률은 이중성 궤도 주기(Pbin ≈ 134 년)에 의해 조절되는 작은 진동을 보이며, 동반성의 근일점 통과와 관련이 있다.
• 입자 수를 증가시킨 시뮬레이션(20,000에서 100,000)에서 원반의 근일점 각도 ωdisk가 수렴함을 확인하여 수치적 안정성을 입증하였다.
• 코드는 느린 동역학과 폭발적인 동역학 모두에서 양호한 수치 성능, 안정성 및 정확도를 보였다.
• 중력의 3D 처리 방식은 수직 운동량의 이동을 가능하게 하여 2D 모델에 존재하지 않는 새로운 자유도를 도입하였다.