[논문 리뷰] Lampray: Multi-group long characteristics ray tracing for adaptive mesh radiation hydrodynamics

Troels Frostholm, Troels Haugbølle외 1명|arXiv (Cornell University)|2018. 09. 14.

Computational Fluid Dynamics and Aerodynamics참고 문헌 1인용 수 4

한 줄 요약

Lampray는 RAMSES 코드에서 적응 메쉬 복사유체역학을 위한 다중군 장특성 빛자국 추적 방법을 도입하여 그림자, 충돌하는 광선, 비평형 화학 반응을 정확하게 처리할 수 있게 한다. H-He-C-O 광화학 반응 네트워크(36종류의 종과 240개 반응)를 통합하고 H₂와 CO의 자가차폐 효과를 고려하며, 먼지 온도 및 얼음 형성도 포함하여 빈자리 분해와 해시 테이블 검색을 통해 효율적인 병렬 처리를 달성한다.

ABSTRACT

We present Lampray: a multi-group long characteristics ray tracing method for adaptive mesh radiation hydrodynamics in the Ramses code. It avoids diffusion, captures shadows, and treats colliding beams correctly, and therefore complements existing moment-based ray tracing in Ramses. Lampray includes different options for interpolation between ray and cell domain, and use either integral, Fourier, or an implicit method for hydrogen ionization to solve the radiative transfer. The opacity can either be tabulated or computed through a coupling to the general non-equilibrium astro-chemistry framework Krome. We use an H-He-C-O network with 36 species and 240 reactions to track the photo-chemistry in the interstellar medium across 6 and 10 orders of magnitude in temperature and density. Self-shielding prescriptions for H$_2$ and CO are used together with a new model for the diffuse interstellar UV-field. We also track the dust temperature, formation of H$_2$ on grains, and H$_2$O and CO ices in detail. Lampray is tested against standard benchmarks for molecular cloud and star formation simulations, including the formation of a Strömgren sphere, the expansion of an ionization front, the photo-evaporation of a dense clump, and the H-He-C-O chemistry in a static photo-dissociation front. Efficient parallelisation is achieved with a separate domain decomposition for rays where points along a ray reside in the same memory space, and data movement from cell- to ray-domain is done with a direct hash-table lookup algorithm. Point sources are treated without splitting rays, and therefore the method currently only scales to a few point sources, while diffuse radiation has excellent scaling.

연구 동기 및 목표

• 적응 메쉬 구조에서 그림자, 충돌하는 광선, 비확산적 운반을 정확히 반영하는 고정밀 복사전달 방법을 개발한다.
• 세밀한 비평형 천체화학(36종의 종과 240개 반응을 포함하는 H-He-C-O 네트워크)을 복사전달과 결합하여 간성간 매질 조건을 정확하게 모델링한다.
• H₂와 CO의 자가차폐, 확산성 간성 우주 자외선 복사장 모델링, 먼지 온도 변화, 그리고 먼지 표면에서의 얼음 형성과 분해 과정을 포함한다.
• 빛자국 영역 분할과 직접 해시 테이블 검색을 통한 세포-빛자국 간 보간을 통해 빛자국 추적의 효율적 병렬 처리를 실현한다.
• 점원천과 확산 복사장이 존재하는 상황에서 복사유체역학 시뮬레이션에 대해 확장 가능한 솔루션을 제공하며, 특히 별 형성 및 원형성 성간 환경에 초점을 맞춘다.

제안 방법

• 연속적인 빛자국을 따라 복사 운반을 계산하기 위해 다중군 장특성 빛자국 추적을 사용하여, 유속 제한된 확산 방법의 비물리적 퍼짐을 피한다.
• 수치적 적분, 푸리에 변환 또는 암묵적 c2ray 방법을 활용해 수소 이온화를 처리하며, 빛자국과 세포 영역 간의 보간 옵션을 제공한다.
• 시간에 따라 변화하는 화학 네트워크를 6–10개 온도 및 밀도 체계에서 해결하기 위해 KROME 천체화학 프레임워크와 통합한다.
• 수정된 Lyα 기반의 제도를 사용해 H₂와 CO의 자가차폐 효과를 모델링하고, 최신 파rameter를 반영해 확산성 간성 우주 자외선 복사장을 모델링한다.
• 먼지 온도 변화를 추적하고, 먼지 표면에서 H₂ 형성 및 H₂O/CO 얼음의 형성과 분해 과정을 위한 표면 화학 반응을 포함한다.
• 빛자국 영역에 대해 별도의 도메인 분할을 통해 빛자국 추적을 병렬 처리하며, 세포에서 빛자국 영역으로의 데이터 이동은 직접 해시 테이블 검색 알고리즘을 사용한다.

실험 결과

연구 질문

• RQ1 적응 메쉬에서 장특성 빛자국 추적은 복사유체역학 시뮬레이션에서 그림자와 충돌하는 광선을 정확히 포착할 수 있는가?
• RQ2 자가차폐 및 먼지 과정을 포함한 비평형 H-He-C-O 화학 반응은 원형성 성간 환경과 원성간 성간 환경의 구조와 진화에 어떤 영향을 미치는가?
• RQ3 병렬 적응 메쉬 프레임워크에서 점원천과 확산 복사장에 대해 빛자국 추적의 성능 스케일링은 어떻게 되는가?
• RQ4 자가차폐 및 최신 자외선 복사장 모델이 간성간 매질 내 이온화 및 분자 화학 반응의 정확도에 어떤 영향을 미치는가?
• RQ5 먼지 온도 및 얼음 형성의 고려가 복사유체역학 시뮬레이션의 물리적 현실성에 얼마나 기여하는가?

주요 결과

• Lampray는 복사유체역학에서 그림자와 충돌하는 광선을 정확히 포착하여, 유속 제한된 확산 방법의 비물리적 퍼짐을 피한다.
• 스트로마르크 영역 형성, 이온화 front 확장, 고밀도 클러스터의 광학적 증발과 같은 벤치마크 테스트를 정확히 재현한다.
• 36종의 종과 240개 반응을 포함한 H-He-C-O 네트워크를 통해 6–10개 온도 및 밀도 체계에서 비평형 화학 반응을 자동으로 일관되게 모델링할 수 있다.
• H₂와 CO의 자가차폐 효과는 효과적으로 모델링되어 저광학적 두께 영역에서 정확도가 향상되고 기존 방법 대비 과도한 이온화가 감소한다.
• 먼지 온도 변화 및 얼음 형성/분해 과정이 일관되게 추적되며, 온도 및 밀도에 따라 변화하는 비율을 기반으로 고체화 및 탈착 속도를 계산한다.
• 해시 테이블 검색을 통한 빛자국 영역 분할은 효율적인 병렬 처리를 가능하게 하며, 확산 복사장에서는 뛰어난 스케일링 성능를 보이지만, 점원천 처리 시 스케일링에 제한이 있다.