[논문 리뷰] A Discontinuous Galerkin Method for General Relativistic Hydrodynamics in thornado

Samuel Dunham, Eirik Endeve외 3명|arXiv (Cornell University)|2020. 09. 27.

Computational Fluid Dynamics and Aerodynamics참고 문헌 55인용 수 3

한 줄 요약

이 논문은 3+1 conformally-flat approximation를 사용한 thornado 프레임워크에서 일반 상대론적 유체역학(GRHD)을 위한 고차수 불연속 갈레르킨(DG) 해법을 제시한다. 이 방법은 충격 동역학에 대해 정확하고 견고한 시뮬레이션을 가능하게 하며, 2차원 GRHD에서 처음으로 정지 수축 충격 불안정성(SASI)을 성공적으로 재현한다. 결과는 이전 연구들과 정성적으로 일치하며, 상대론적 조건 하에서도 뛰어난 성능을 보여준다.

ABSTRACT

Discontinuous Galerkin (DG) methods provide a means to obtain high-order accurate solutions in regions of smooth fluid flow while, with the aid of limiters, still resolving strong shocks. These and other properties make DG methods attractive for solving problems involving hydrodynamics; e.g., the core-collapse supernova problem. With that in mind we are developing a DG solver for the general relativistic, ideal hydrodynamics equations under a 3+1 decomposition of spacetime, assuming a conformally-flat approximation to general relativity. With the aid of limiters we verify the accuracy and robustness of our code with several difficult test-problems: a special relativistic Kelvin--Helmholtz instability problem, a two-dimensional special relativistic Riemann problem, and a one- and two-dimensional general relativistic standing accretion shock (SAS) problem. We find good agreement with published results, where available. We also establish sufficient resolution for the 1D SAS problem and find encouraging results regarding the standing accretion shock instability (SASI) in 2D.

연구 동기 및 목표

• thornado 계산 프레임워크에서 일반 상대론적 유체역학(GRHD)을 위한 고차수 불연속 갈레르킨(DG) 해법을 개발한다.
• 특히 충격 형성과 진화를 포함한 핵붕괴 초신성(CCSN) 동역학을 일반 상대론적 조건 하에서 정확하고 안정적으로 시뮬레이션할 수 있도록 한다.
• 충격 불안정성과 강한 중력 영역과 같은 도전적인 상대론적 유체역학 문제에 대해 해법의 정확성과 견고성을 시험한다.
• GRHD, 중성미자 수송, 적응 메쉬 다듬힘(AMR)을 포함한 향후 3차원 CCSN 시뮬레이션을 위한 기반을 마련한다.
• 고차수 DG 방법을 사용하여 일반 상대론이 정지 수축 충격 불안정성(SASI)에 미치는 영향을 탐구한다.

제안 방법

• 해법은 시공간의 3+1 분해와 conformally-flat approximation(CFA)을 사용하여 GRHD를 위한 아인슈타인 장 방정식을 단순화한다.
• GRHD 방정식은 밸렌시아 형태로 표현되며, 고차수 정확도를 가지는 불연속 갈레르킨(DG) 유한요소 방법을 사용하여 해결한다. 이는 컴 pact 스텐실에서 수행된다.
• 보존 변수(D, Sj, τ)는 유량과 소스 항의 약한 형태를 사용하여 진화시키며, 충격 해소를 위해 Rusanov 유량과 기울기 제한기를 사용한 재구성 기법을 적용한다.
• 해법은 hp-적응성과 AMR를 위한 AMReX 프레임워크를 통해 MPI를 이용한 병렬 처리를 구현한다.
• 원시 변수는 보존 변수에서 뉴턴-라프슨 근 찾기 방법을 사용하여 재구성되며, 충격을 넘어선 안정성을 유지하기 위해 제한기 조치를 적용한다.
• 테스트 문제를 사용하여 코드를 검증하였으며, 이는 특수 상대론적 켈빈-헬름홀츠, 리만, 일반 상대론적 정지 수축 충격(SAS) 문제를 포함한다.

실험 결과

연구 질문

• RQ1 고차수 불연속 갈레르킨 방법은 일반 상대론적 영역에서 강한 충격과 복잡한 유체역학을 얼마나 정확하게 해석할 수 있는가?
• RQ2 일반 상대론은 2차원 시뮬레이션에서 정지 수축 충격 불안정성(SASI)의 발달과 형태에 어느 정도 영향을 미치는가?
• RQ3 퍼티어션에 의해扰동된 상황에서도 thornado DG 해법이 고차수 정확도를 유지하면서 정적 충격 프로파일을 유지할 수 있는가?
• RQ4 시공간 곡률과 유체 루프터츠 인자와 같은 상대론적 효과가 뉴턴 근사와 비교해 충격 동역학에 어떤 영향을 미치는가?
• RQ5 MPI와 AMR를 사용한 다차원 GRHD 시뮬레이션에서 DG 해법의 성능과 확장성은 어떠한가?

주요 결과

• DG 해법은 과밀도 쉘에 의한 교란 조건 하에서도 정적 충격 프로파일을 유지하면서 1차원 정지 수축 충격(SAS) 문제를 높은 정확도로 재현하였다.
• 이 코드는 일반 상대론적 유체역학에서 처음으로 2차원 정지 수축 충격 불안정성(SASI)을 정확하게 포착하였다. 결과는 이전에 발표된 연구와 정성적으로 일치하였다.
• 시뮬레이션 결과, 내부 경계에서 랩스 함수 α가 1에서 10% 정도의 편차를 보이며, 이는 충격 동역학과 유체 속도에 상당한 영향을 미친다는 것을 확인하였다.
• 루프터츠 인자 W와 특별한 엔탈피 h는 모두 1에서 크게 벗어나 있어, 시공간 곡률과 유체 운동 모두 강한 상대론적 영역에서 작동하고 있음을 확인하였다.
• 포lytropic 상수(엔트로피의 대리 척도)는 2차원에서 동적으로 변화하며, 시간이 지남에 따라 충격이 북극과 남극 쪽으로 흔들리는 현상이 명확히 관측되었으며, 이는 SASI 발달과 일치하였다.
• MPI와 AMReX를 사용한 구현은 견고성과 확장성을 입증하였으며, 향후 3차원 CCSN 시뮬레이션의 가속화를 위한 GPU 이식에 유망한 결과를 보였다.