그래픽 처리 장치 가속 적응형 격자 정밀 시뮬레이션 코드

초록

GAMER는 3차원 적응형 격자(AMR) 구조와 옥트리 트리를 기반으로, 수소역학과 포아송 방정식을 GPU에서 병렬 처리하도록 설계된 코드이다. 비동기 메모리 복사와 CPU‑GPU 연산 겹침을 통해 데이터 전송 오버헤드를 최소화하고, 다중 GPU 클러스터에서도 높은 확장성을 보인다. 표준 테스트와 대규모 우주론 시뮬레이션에서 1 GPU당 최대 12배, 16 GPU 구성에서는 10배 이상의 속도 향상을 입증하였다.

상세 분석

GAMER는 기존 AMR 코드가 직면한 두 가지 핵심 병목 현상을 해결한다. 첫째, 격자 패치 단위의 연산을 GPU 코어 수천 개에 동시에 할당함으로써 패치당 연산량을 크게 늘렸다. 이를 위해 3차원 Relaxing TVD(총변량 보존) 스킴을 GPU 전용 커널로 구현했으며, 고차원 데이터 구조를 연속 메모리 블록에 재배치해 메모리 접근 패턴을 최적화하였다. 둘째, 중력 포텐셜을 구하는 다중 레벨 완화(멀티레벨 릴랙세이션) 방식은 각 레벨마다 경계 조건을 교환해야 하는데, GAMER는 이 과정을 CPU와 GPU 사이의 비동기 복사 스트림을 이용해 겹쳐 수행한다. 즉, GPU가 현재 패치를 계산하는 동안 CPU는 다음 패치의 고스트 영역을 준비하고, 완성된 결과는 즉시 다음 단계로 전파된다. 이러한 파이프라인 구조는 데이터 전송 대기 시간을 실질적으로 0에 가깝게 만든다.

또한, GAMER는 멀티 GPU 환경에서 MPI 기반 도메인 분할을 적용했으며, 각 노드가 담당하는 옥트리 서브트리를 독립적으로 관리한다. GPU 간 직접 통신은 없고, CPU가 중계 역할을 수행해 복잡한 동기화 문제를 회피한다. 실험 결과, 4096³ 격자 규모에서 1 GPU 사용 시 12.19배, 8192³ 규모에서 16 GPU 사용 시 10.47배의 가속을 달성했으며, 이는 전통적인 CPU 전용 AMR 코드 대비 1~2자리 수의 성능 향상이다.

코드 정확도 검증을 위해 Sod 충격파, Sedov 폭발, 블랙홀 주변 중력 포텐셜 등 표준 테스트를 수행했으며, 수치 오차는 기존 고정 격자 코드와 동등하거나 더 낮았다. 특히, 포아송 솔버의 다중 레벨 완화는 경계 조건 불연속성을 최소화해 중력 장의 정밀도를 유지한다.

이러한 설계는 향후 복잡한 물리 모듈(자기장, 방사선 전달, 화학 반응 등)을 GPU에 추가하기 위한 기반을 제공한다. 현재는 순수 바리온 물리만 구현돼 있지만, 모듈화된 구조와 추상화된 데이터 인터페이스 덕분에 확장이 용이하다.

요약하면, GAMER는 GPU의 대규모 병렬성을 AMR의 계층적 데이터 구조와 결합해, 메모리 전송 최소화와 연산 겹침을 통한 효율적인 파이프라인을 구현함으로써 천체 물리 시뮬레이션의 계산 비용을 획기적으로 낮추었다.