[논문 리뷰] A New Public Release of the GIZMO Code

Philip F. Hopkins|arXiv (Cornell University)|2017. 12. 04.

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

한 줄 요약

이 논문은 다중 물리 현상 시뮬레이션 코드인 GIZMO 코드의 대규모 공개를 발표한다. 이 코드는 적응형 메시 없는 유한체적 유체역학, SPH, MHD, 이방성 열전도/ viscoity, 항성 및 블랙홀 형성, 먼지-기체 동역학, 탄성/플라스틱 고체, 임의의 상태방정식을 지원하는 대규모 병렬 처리가 가능한 천체물리 시뮬레이션 코드이다. 이 공개는 100만 개 이상의 CPU 스레드까지 거의 이상적인 강한 및 약한 스케일링을 가능하게 하며, 광범위한 문서화, 테스트 문제, GADGET 도구와의 호환성을 제공한다.

ABSTRACT

We describe a major update to the public GIZMO code. GIZMO has been used in simulations of cosmology; galaxy and star formation and evolution; black hole accretion and feedback; proto-stellar disk dynamics and planet formation; fluid dynamics and plasma physics; dust-gas dynamics; giant impacts and solid-body interactions; collisionless gravitational dynamics; and more. This release of the public code supports: hydrodynamics (using various mesh-free finite-volume Godunov methods or SPH), ideal and non-ideal MHD, anisotropic conduction and viscosity, radiative cooling and chemistry, star and black hole formation and feedback, sink particles, dust-gas (aero)-dynamics (with or without magnetic fields), elastic/plastic dynamics, arbitrary (gas, stellar, degenerate, solid/liquid material) equations of state, passive scalar/turbulent diffusion, large-eddy and shearing boxes, self-gravity with fully-adaptive force softenings, arbitrary cosmological expansion, and on-the-fly group-finding. It is massively-parallel with hybrid MPI+OpenMP scaling verified up to >1 million threads. The code is extensively documented, with test problems and tutorials provided for these different physics modules.

연구 동기 및 목표

• 우주론, 은하 형성, 항성 및 행성 형성, 플라즈마 물리학 등 다양한 천체물리 문제를 위한 종합적이고 공개된, 그리고 매우 확장 가능한 시뮬레이션 코드를 제공하기 위해.
• 자기장, 피드백, 다상유체역학 등 고급 물리 모듈을 사용하여 대규모 시뮬레이션을 수행할 수 있도록 하는 탄력적이고 모듈식 프레임워크를 제공하기 위해.
• 광범위한 문서화, 튜토리얼, 테스트 문제, 기존 GADGET 워크플로우 및 데이터 포맷과의 호환성을 통해 보편적인 접근성과 재현 가능성을 확보하기 위해.
• 엑사스케일 수준의 슈퍼컴퓨터에서 거의 이상적인 강한 및 약한 스케일링을 달성하여 10억 개 이상의 해상도 요소를 포함한 시뮬레이션을 지원하기 위해.
• 공개 및 비공개 개발 브랜치를 분리하여 코드의 확장성과 미래 지향성을 유지하고, 기밀 또는 불안정한 물리 모듈의 책임 있는 공유를 보장하기 위해.

제안 방법

• GIZMO는 대규모 슈퍼컴퓨터에서 효율적인 확장성을 확보하기 위해 적응형 도메인 분할과 계층적 타임스텝 전략을 결합한 하이브리드 MPI+OpenMP 병렬 처리 기법을 사용한다.
• 코드는 메시 없는 유한체적 Godunov 방법, SPH, 고정된 직각좌표 격자 등 다양한 수치 방법을 지원하며, 사용자가 정의한 격자 운동(예: 비틀림, 팽창, 또는 회전)을 허용한다.
• 완전히 적응형 중력 소프트닝을 지원하는 임의의 라그랑주-오일러 방법을 구현하였으며, 고체에 대한 Tillotson 형상방정식과 기초 물질에 대한 Helmholtz 형상방정식을 포함한 임의의 상태방정식을 지원한다.
• 핵심 물리 모듈로는 실시간 군집 탐색, 싱크 입자 형성, 자기 일관성 있는 항성 및 블랙홀 피드백, 이방성 열전도 및 점성도, Smagorinski 모델을 통한 난류 확산이 포함되어 있다.
• 먼지-기체 동역학은 선택 가능한 마찰 법칙(에프스타인, 스토크스, 쿨롱)과 루프르 힘을 모델링하며, 기체의 운동량 방정식에 반작용 효과를 포함한다.
• 우주론적 시뮬레이션은 비표준 팽창 역사, 동적 암흑 에너지, 임의의 경계 조건(주기적, 개방형, 반사형)을 지원하며, 실시간 파wer 스펙트럼 및 군집 탐색 기능을 제공한다.

실험 결과

연구 질문

• RQ1 어떻게 한 개의 시뮬레이션 코드가 우주론에서부터 행성 형성에 이르기까지 광범위한 천체물리 문제를 효율적으로 지원할 수 있는가?
• RQ2 유체역학, 중력, 자기장, 피드백을 포함한 다중 물리 시뮬레이션에서 어떤 정도의 병렬 스케일링 성능를 달성할 수 있는가?
• RQ3 어떻게 적응형 메시 없는 방법과 라그랑주 입자 추적을 조합하여 복잡한 유체 및 고체 동역학을 모델링할 수 있는가?
• RQ4 GIZMO는 강한 자기장 난류, 이방성 열전도, 다상간 우주간 매질을 시뮬레이션할 때 성능와 정확도가 어느 정도인가?
• RQ5 어떻게 모듈식이고 확장 가능한 물리 모듈을 하나의 확장 가능하고 잘 문서화된 시뮬레이션 프레임워크에 통합할 수 있는가?

주요 결과

• GIZMO는 우주론적 시뮬레이션에서 100만 개 이상의 CPU 스레드까지 거의 이상적인 약한 스케일링을 달성하였으며, 2에서 10,000 Mpc³의 다양한 부피에서 성능이 안정되어 있다.
• 밀려는 질량과 소형 은하의 줌인 시뮬레이션 모두에서 1,500만 개의 입자를 사용하여 강한 스케일링이 거의 이상적이었으며, 25만 개의 CPU와 16,000개의 노드를 통해 효율성이 유지되었다.
• 코드는 오르사그-탕 풍압, 자기열 불안정성, 켈빈-헬름홀츠 불안정성과 같은 복잡한 현상을 고정밀도 및 고해상도로 정확하게 시뮬레이션하였다.
• GIZMO는 난류가 있는 자기장이 있는 구름에서 먼지-기체 동역학을 정확하게 모델링하였으며, 에프스타인 및 루프르 힘에 따른 입자 운동과 기체에 대한 반작용 효과를 포함하였다.
• 중력에 의한 블랙홀 수축과 피드백 유도 유출을 통해 블랙홀의 고해상도 적층을 가능하게 하였으며, ISM에 정확한 에너지 및 운동량 투입을 수행하였다.
• Tillotson 형상방정식을 사용한 탄성 및 플라스틱 동역학을 포함함으로써 거대 충돌과 행성체 충돌에서의 물질 혼합 및 항복 조건을 현실적으로 시뮬레이션할 수 있었다.