시간 의존 문제를 위한 적응형 도메인 분할 기법과 유체 역학 적용

초록

본 논문은 시간 의존 편미분 방정식의 수치 해석에 공간‑시간 불연속 갈루아(SDG) 방법을 적용하고, 적응형 메쉬와 뉴턴‑유사 반복을 결합한 뒤, GMRES와 두 단계 슈워츠(Schwarz) 프리컨디셔너를 이용한 도메인 분할(preconditioner) 전략을 제시한다. 비용 모델을 기반으로 서브도메인 수와 거친 격자 요소 수를 자동 선택하는 적응형 도메인 분할 기법을 설계하고, 압축성 유체 흐름 두 가지 벤치마크에서 효율성을 검증한다.

상세 분석

이 연구는 시간‑공간 불연속 갈루아(STDG) 방법의 장점인 고차 정확도와 h‑p 적응성을 유지하면서, 비선형 시스템을 뉴턴‑유사 방법으로 선형화하고, 각 선형 단계에서 GMRES와 두‑레벨 슈워츠 프리컨디셔너를 적용한다는 점에서 독창적이다. 특히, 서브도메인 수 M과 거친 격자 요소 수 n₀(=s·M)을 조절하는 두 개의 도메인 분할 파라미터를 도입해, M이 증가하면 병렬 효율은 높아지지만 수렴 속도가 저하되고, s가 커지면 전역 정보 전달이 강화돼 수렴이 빨라지지만 연산량이 늘어나는 상쇄 관계를 정량화한다. 이를 위해 저자들은 부동소수점 연산 수(FLOPs), 코어당 연산 속도, 통신 지연시간을 포함한 ‘비용 모델’을 제안한다. 모델은 각 서브도메인과 거친 격자에 대한 행렬 규모, 스파스 구조, 그리고 MPI 기반 통신 비용을 추정해, 전체 시뮬레이션 단계별 예상 비용을 사전에 계산한다.

비용 모델을 기반으로 최적의 (M, s) 조합을 선택하는 적응형 알고리즘은 다음과 같이 동작한다. ① 현재 메쉬와 시간 단계에서 예상 DoF와 연산량을 평가한다. ② 사전 실험으로 얻은 GMRES·슈워츠 수렴 특성(예: 조건수, 반복 횟수)을 파라미터 M, s에 대한 함수 형태로 보간한다. ③ FLOPs와 통신 비용을 합산한 총 비용 함수를 최소화하는 (M, s)를 선택한다. ④ 선택된 파라미터에 따라 METIS를 이용해 서브도메인을 재분할하고, 거친 격자를 재구성한다. 이 과정은 메쉬 적응이 일어날 때마다 반복된다.

수치 실험에서는 2차원 압축성 Navier‑Stokes 방정식과 비수평 대기 모델을 대상으로, (i) 고속 충격 파동 전파와 (ii) 대규모 대기 대류 현상을 시뮬레이션하였다. 두 사례 모두 적응형 메쉬가 급격한 해석학적 변화(충격 전선, 온도 구름)를 효과적으로 포착했으며, 제안된 적응형 도메인 분할이 고정 파라미터 대비 평균 30%~45%의 실행 시간 절감을 달성했다. 특히, 통신 비용이 전체 시간의 20% 이하로 억제된 점은 대규모 병렬 환경에서의 확장성을 시사한다.

이 논문의 주요 기여는 (1) STDG와 비중첩 두‑레벨 슈워츠 프리컨디셔너를 결합한 효율적인 선형 솔버 프레임워크, (2) FLOPs·통신 기반 비용 모델을 통한 사전 비용 예측, (3) 비용 모델을 활용한 자동 파라미터 선택 적응형 도메인 분할 전략, (4) 압축성 유체와 대기 역학 두 분야에 걸친 실증적 검증이다. 또한, 제시된 비용 모델은 다른 비선형 시공간 적응법에도 일반화 가능하므로, 차후 고차원·고복잡도 문제에 대한 병렬 효율 최적화 연구에 중요한 토대를 제공한다.