고차원 적응형 메쉬 ADER WENO 유한체 스킴

초록

본 논문은 다중 차원에서 고차 정확도를 갖는 ADER‑WENO 일단계 유한체 스킴에 셀‑단위 적응형 메쉬 정제(AMR)를 결합한다. 공간 차원에서는 WENO 재구성을, 시간 차원에서는 로컬 시공간 DG 예측자를 이용해 일단계 고차 시간 적분을 구현한다. 트리 기반 셀‑별 AMR과 MPI 병렬화를 적용해 지역 시간 보폭을 유지하면서 효율적인 계산을 가능하게 하였으며, 압축성 유체역학 및 MHD 방정식 등 비선형 보존법칙에 대한 광범위한 테스트를 통해 수렴성, 속도 향상 및 고차 정확도를 입증하였다.

상세 분석

이 연구는 기존의 2차 정확도 AMR 기법이 갖는 시간·공간 비효율성을 극복하기 위해, ADER(Arbitrary high‑order DERivative)와 WENO(Weighted Essentially Non‑Oscillatory)를 결합한 일단계 고차 스킴을 제안한다. WENO 재구성은 셀 평균값을 기반으로 고차 다항식 근사를 수행하며, 급격한 불연속 근처에서도 비진동성을 유지한다. 시간 적분은 전통적인 Runge‑Kutta와 달리 로컬 시공간 DG 예측자를 사용해 각 셀마다 독립적인 고차 다항식 해를 구한다. 이 예측자는 초기 데이터와 경계 조건을 받아, 공간‑시간 부피에서 강인한 해를 제공함으로써 전체 스킴을 일단계로 진행할 수 있게 한다.

AMR 구현은 ‘cell‑by‑cell’ 방식의 2:1 균형 트리를 채택했으며, 각 레벨마다 고유의 시간 단계 Δt_l을 부여한다. 일단계 특성 덕분에 서로 다른 레벨 간의 인터페이스에서 시간 동기화가 복잡하지 않으며, ‘시간‑정확한 로컬 타임 스텝핑(local time stepping)’이 자연스럽게 구현된다. 또한, MPI 기반 도메인 분할을 통해 다중 코어·다중 노드 환경에서 확장성을 확보했으며, 통신 오버헤드를 최소화하기 위해 경계 셀만 교환하도록 설계하였다.

수치 실험에서는 2D·3D 유체역학 및 MHD 문제를 대상으로 정밀한 수렴 테스트와 복잡한 파동·충격 파동 전파 사례를 수행했다. 정규 격자와 비교했을 때, 동일한 정확도를 유지하면서 평균 8배 이상의 연산 속도 향상을 기록했으며, 특히 고차 재구성으로 인해 얇은 경계층이나 강한 전단 영역에서도 스키밍 현상이 현저히 감소하였다. 이러한 결과는 고차 ADER‑WENO와 AMR의 결합이 고성능 컴퓨팅 환경에서 비선형 보존법칙을 효율적으로 해결할 수 있음을 강력히 시사한다.