AthenaK 기반 조석 파괴 현상(TDE) 통합 수치 시뮬레이션 프레임워크 공개
연구진은 GPU 가속 유한체적 코드 AthenaK를 확장해, 별이 블랙홀의 조석력에 찢긴 뒤 잔해가 강착원반을 형성하기까지의 과정을 자체 중력을 포함한 단일 시뮬레이션으로 추적하는 프레임워크를 제시했다. 이동 좌표계, 적응 격자 기반 멀티그리드 푸아송 솔버, 재결합을 포함한 수소-헬륨 상태방정식, 이중 에너지 갱신, 국소 적응 시간 적분(LAT) 등이 추가됐으며 LAT만으로 계산 속도가 3배 이상 빨라졌다. 중력 솔버는 16.3 역학시간 동안 Lane-Emden 구의 반경 방향 밀도 오차를 2.0e-3로 유지했고, 이중 에너지 기법은 마하 7.75e7 엔트로피 파동의 압력 오차를 15.7에서 4.0e-12로 낮췄다. 시연으로 태양질량 1배·태양반지름 1배의 별이 10^3 태양질량 블랙홀에 beta=1로 파괴되는 경우를 계산해, 잔해가 속박·비속박 물질로 균등하게 나뉘고 후기 낙하율이 t^(-5/3)에 접근함을 확인했다. 근점 노즐에서는 고해상도 계산의 열에너지 획득량이 수직 운동에너지 손실과 5% 이내로 일치한 반면, 기준 해상도 계산은 수치 소산으로 초기 얇은 스트림을 과열시켰다.
arXiv astro-ph.HE (고에너지)
