[논문 리뷰] Cosmological Hydrodynamics with Adaptive Mesh Refinement: a new high resolution code called RAMSES

Romain Teyssier|arXiv (Cornell University)|2001. 11. 19.

Galaxies: Formation, Evolution, Phenomena참고 문헌 50인용 수 1,070

한 줄 요약

이 논문은 적응 격자 구조를 사용하는 새로운 천체역학 유체역학 코드인 RAMSES를 제시한다. 이 코드는 동적 격자 정밀도 향상을 위해 트리 기반 데이터 구조를 사용하며, 2차 Godunov 수치 해법과 TREE 기반 N-체 해법을 결합하여, 낮은 밀도 ΛCDM 우주론에서 구조 형성의 고해상도 시뮬레이션을 가능하게 한다. 이는 공식적인 해상도 8192³에 도달하며, 해결한 해상도 한계까지 수렴성을 보이며, 어둠센 물질 및 기체의 파wer 스펙트럼에 대해 헬로 모델 예측과 양호한 일치를 보인다.

ABSTRACT

A new N-body and hydrodynamical code, called RAMSES, is presented. It has been designed to study structure formation in the universe with high spatial resolution. The code is based on Adaptive Mesh Refinement (AMR) technique, with a tree based data structure allowing recursive grid refinements on a cell-by-cell basis. The N-body solver is very similar to the one developed for the ART code (Kravtsov et al. 97), with minor differences in the exact implementation. The hydrodynamical solver is based on a second-order Godunov method, a modern shock-capturing scheme known to compute accurately the thermal history of the fluid component. The accuracy of the code is carefully estimated using various test cases, from pure gas dynamical tests to cosmological ones. The specific refinement strategy used in cosmological simulations is described, and potential spurious effects associated to shock waves propagation in the resulting AMR grid are discussed and found to be negligible. Results obtained in a large N-body and hydrodynamical simulation of structure formation in a low density LCDM universe are finally reported, with 256^3 particles and 4.1 10^7 cells in the AMR grid, reaching a formal resolution of 8192^3. A convergence analysis of different quantities, such as dark matter density power spectrum, gas pressure power spectrum and individual haloes temperature profiles, shows that numerical results are converging down to the actual resolution limit of the code, and are well reproduced by recent analytical predictions in the framework of the halo model.

연구 동기 및 목표

• ΛCDM 우주론에서 대규모 구조와 소규모 복사 물리 현상을 모두 해상도 해석이 가능한 고해상도 천체역학 시뮬레이션 코드를 개발하는 것.
• 기존의 N-체 및 에일러리안 유체역학 방법의 한계를 극복하여 고역동 범위와 정확한 충격 포착 능력을 확보하는 것.
• 셀당 입자 수를 일정하게 유지하는 적응 격자 정밀도 전략을 구현하여 이중 입자 상호작용과 포isson 노이즈를 최소화하는 것.
• 표준 유체역학 및 천체역학 테스트 케이스를 사용하여 코드의 정확성과 수렴성을 검증하는 것.
• 어둠센 물질 및 기체의 파워 스펙트럼에 대해 헬로 모델의 분석 예측을 재현할 수 있는 능력을 입증하는 것.

제안 방법

• 코드는 트리 기반 데이터 구조를 사용하여 AMR 격자의 재귀적이고 셀 단위의 정밀도 향상을 가능하게 하여, 밀도가 높은 영역에서 고공간 해상도를 확보한다.
• N-체 해법은 장거리 힘을 위해 다중극 전개를 사용하고, 가장 정밀한 수준에서는 직접 입자-입자 상호작용을 수행하는 TREE 알고리즘에 기반하여 고역동 범위를 확보한다.
• 유체역학 해법은 Riemann 해법을 사용하는 2차 Godunov 방법을 적용하여 충격 포착의 정확성과 총 에너지 보존(중력 제외)을 보장한다.
• 준라그랑주 정밀도 전략은 정밀도 수준 간에 셀당 입자 수를 근사적으로 일정하게 유지하여 수치 노이즈와 이중 입자 상호작용을 감소시킨다.
• 계층적 격자 구조와 중첩된 격자를 사용하여, 적응 격자에서의 힘과 유체역학 계산을 효율적으로 수행할 수 있다.
• 시뮬레이션은 고정된 박스 크기(100 h⁻¹ Mpc)로 시작하며, 초기 격자 해상도를 다양하게 설정(64³에서 256³)하여 최고 해상도에서 공식적인 해상도 8192³를 달성한다.

실험 결과

연구 질문

• RQ1 트리 기반 데이터 구조를 갖는 고해상도 AMR 코드는 천체역학적 구조 형성 시뮬레이션에서 정확하고 수렴 가능한 결과를 도출할 수 있는가?
• RQ2 준라그랑주 정밀도 전략은 N-체 시뮬레이션에서 수치 노이즈와 이중 입자 상호작용에 어떤 영향을 미치는가?
• RQ3 고해상도에서 어둠센 물질 및 기체의 파워 스펙트럼에 대한 수치 결과가 어느 정도 수렴하며, 분석적 헬로 모델 예측과 일치하는가?
• RQ4 AMR 격자 경계를 가로질러 충격파가 전파될 때 발생하는 허구적 효과는 무엇이며, 이는 무시할 수 있는가?
• RQ5 헬로 모델을 확장하여 시뮬레이션된 은하단의 기체 밀도 및 온도 프로파일을 정확하게 기술할 수 있는가?

주요 결과

• RAMSES 코드는 100 h⁻¹ Mpc 박스에서 공식적인 공간 해상도 8192³를 달성하며, 이는 12 h⁻¹ kpc의 공진동 해상도에 해당한다.
• 어둠센 물질 밀도 파wer 스펙트럼에 대한 수치 결과는 공식적인 해상도 한계까지 수렴하며, 분석적 헬로 모델 예측과 양호한 일치를 보인다.
• 기체 압력 파워 스펙트럼의 경우, 최고 해상도 시뮬레이션에서 50 h⁻¹ kpc 이상의 스케일에서 수렴이 관찰된다.
• 헬로 모델은 모든 스케일에서 해결한 해상도 한계까지 어둠센 물질 파워 스펙트럼을 몇 퍼센트 이내로 재현한다.
• 50 h⁻¹ kpc 이상의 스케일에서 기체 파워 스펙트럼에 대해 헬로 모델은 약 2배 이내로 근사하며, 기체 밀도 프로파일에 β-모델을 도입한 확장된 모델이 일치도를 향상시킨다.
• 정적 평형과 β-모델 기반 기체 밀도 프로파일에서 유도된 분석적 온도 프로파일은 개별 은하단의 시뮬레이션 온도 프로파일과 정확히 일치한다.