[논문 리뷰] Simulating dust grain-radiation coupling on a moving mesh

R. McKinnon, Rahul Kannan외 4명|arXiv (Cornell University)|2019. 12. 05.

Astrophysics and Star Formation Studies참고 문헌 2인용 수 6

한 줄 요약

이 논문은 arepo-rt 복사유체역학 코드에서 실시간으로 동작하는 먼지 입자를 비정형 이동 격자에 연결하는 새로운 하이브리드 방법을 제안한다. 먼지 투과도, 복사압, 광자 흡수 및 열방출을 자성적으로 계산함으로써, 局부 열역학 평형을 가정하지 않고 실제 먼지-기체-복사 에너지 및 운동량 교환을 가능하게 한다. 이는 정확한 속도 및 온도 변화를 보장하는 해석적 해와의 검증을 통과하였다.

ABSTRACT

We present a model for the interaction between dust and radiation fields in the radiation hydrodynamic code AREPO-RT, which solves the moment-based radiative transfer equations on an unstructured moving mesh. Dust is directly treated using live simulation particles, each of which represent a population of grains that are coupled to hydrodynamic motion through a drag force. We introduce methods to calculate radiation pressure on and photon absorption by dust grains. By including a direct treatment of dust, we are able to calculate dust opacities and update radiation fields self-consistently based on the local dust distribution. This hybrid scheme coupling dust particles to an unstructured mesh for radiation is validated using several test problems with known analytic solutions, including dust driven via spherically-symmetric flux from a constant luminosity source and photon absorption from radiation incident on a thin layer of dust. Our methods are compatible with the multifrequency scheme in AREPO-RT, which treats UV and optical photons as single-scattered and IR photons as multi-scattered. At IR wavelengths, we model heating of and thermal emission from dust. Dust and gas are not assumed to be in local thermodynamic equilibrium but transfer energy through collisional exchange. We estimate dust temperatures by balancing these dust-radiation and dust-gas energy exchange rates. This framework for coupling dust and radiation can be applied in future radiation hydrodynamic simulations of galaxy formation.

연구 동기 및 목표

• 복사유체역학 시뮬레이션에서 먼지-복사 결합을 물리적으로 정확하고 자성적으로 일관된 모델을 개발하기 위해.
• 고정된 먼지-기체 비율 또는 局부 열역학 평형을 가정하는 하위격자 모델의 한계를 극복하기 위해.
• 먼지가 복사압과 저항력에 반응하는 실시간 시뮬레이션 입자로 간주되어 동적 먼지 거동을 가능하게 하기 위해.
• 다중 주파수 복사장에서 먼지의 가열, 열방출 및 광자 흡수 모델링을 향상시키기 위해.
• 구형 대칭 복사장에서 복사압에 의해 구동되는 먼지와 얇은 층에서의 광자 흡수를 포함한 해석적 테스트 문제에 대한 프레임워크의 검증을 위해.

제안 방법

• 먼지가 각기 질량, 속도 및 입자 크기에 따라 달라지는 단면적을 가진 실시간 시뮬레이션 입자로 모델링된다.
• 복사압력은 입자로부터 기체 요소로 보간된 먼지 단면적을 사용하여 계산되며, 이는 마찰을 통해 기체로 운동량 전달을 가능하게 한다.
• 광자 흡수량은 먼지 투과도를 기체 요소로 커널 스무딩 보간하여 계산되며, 복사장이 자성적으로 갱신된다.
• arepo-rt의 M1 폐쇄 체계는 다주파수 복사를 처리하며, 단일 산란된 UV/광학 및 다중 산란된 적외선 광자를 구분한다.
• 먼지와 기체 온도는 별도로 진화시키며, 충돌 과정과 복사 흡수 및 열방출을 통한 먼지 가열/냉각을 통해 에너지 교환을 수행한다.
• 이 방법은 입자-셀 하이브리드 접근법을 사용한다: 먼지 입자는 유체역학과 함께 이동하며, 복사는 비정형 이동 격자에서 해결된다.

실험 결과

연구 질문

• RQ1 구형 대칭 복사장에서 하이브리드 입자-셀 방법이 먼지 입자에 작용하는 복사압을 얼마나 정확하게 시뮬레이션할 수 있는가?
• RQ2 고정된 투과도를 가정하지 않고도 이동 격자 복사유체역학 코드에서 먼지의 광자 흡수를 자성적으로 모델링할 수 있는가?
• RQ3 비평형 먼지-기체 에너지 교환을 포함할 경우, 먼지 및 기체의 역학적 및 열적 진화에 어떤 영향을 미치는가?
• RQ4 입자 크기 분포의 변화가 복사 dominant 환경에서 복사압 및 운동량 전달에 어떤 영향을 미치는가?
• RQ5 복사압과 흡수에 의해 구동되는 먼지 가속도 및 기체와의 속도 동기화를 위한 해석적 해를 모델이 얼마나 잘 재현하는가?

주요 결과

• 고정된 빛의 강도를 가진 소스에 의한 먼지 가속도에 대해 해석적 해를 성공적으로 재현하였다. 초기에는 복사압으로 인해 먼지와 기체 속도가 분리되며, 이후 마찰을 통해 점차 동기화된다.
• 광자 흡수 없이 복사압만 작용할 경우, 먼지와 기체 속도는 무한히 증가하며, 약 108 km s⁻¹의 속도 차이를 보이며 해석적 예측과 일치한다.
• 광자 흡수를 고려할 경우, 먼지와 기체 속도는 약 44 km s⁻¹에서 동기화되며, 운동량 전달 및 에너지 흡수에 대한 해석적 해와 일치한다.
• 얇은 먼지 층을 통과하는 광자 흡수 현상을 정확히 캐릭터라이즈하였으며, 광학적 깊이 의존성이 해석적 기대와 일치한다.
• 먼지와 기체가 局부 열역학 평형에 있다고 가정하지 않으며, 별도의 온도를 진화시키며 충돌 과정을 통해 에너지 교환을 수행한다.
• 지역적 입자 농도에 기반한 자성적인 먼지 투과도 계산을 가능하게 하여 먼지-기체 비율 또는 분포에 대한 가정을 피한다.