불완전 압축성 Navier Stokes 방정식을 위한 단일 및 블록 겹침 Schwarz 프리컨디셔너

초록

본 논문은 불완전 압축성 Navier‑Stokes 방정식의 선형화 시스템에 적용되는 단일(monolithic) 프리컨디셔너와 전통적인 블록 프리컨디셔너(PCD, SIMPLEC, LSC)를 비교한다. 두 레벨 가산 겹침 Schwarz 방법(OSM)과 GDSW 계열의 coarse space를 활용해 고성능 병렬 구현을 수행하고, 새로운 GDSW* coarse space와 압력‑속도 전용 coarse space 조합이 높은 Reynolds와 CFL 수에서도 우수한 수렴성과 스케일러빌리티를 보임을 실험적으로 입증한다.

상세 분석

이 연구는 불완전 압축성 Navier‑Stokes 방정식의 saddle‑point 형태를 직접 다루는 단일 프리컨디셔너와, Schur complement 근사를 기반으로 하는 블록 프리컨디셔너를 체계적으로 비교한다. 단일 접근법은 전체 시스템 행렬 F 에 대해 두 레벨 가산 겹침 Schwarz(O​SM) 프레임워크를 적용한다. 첫 번째 레벨은 각 서브도메인에 대한 로컬 사전해를 수행하고, 두 번째 레벨에서는 GDSW, RGDSW, 그리고 새롭게 제안된 GDSW* coarse space를 이용해 전역적인 저주파 성분을 보정한다. 특히 GDSW*는 기존 GDSW의 자유도 선택을 재구성하여 속도 변수에 대해 더 풍부한 coarse basis를 제공하면서도 압력 변수에 대해서는 RGDSW와 결합해 효율적인 차원 축소를 달성한다.

블록 프리컨디셔너는 PCD, SIMPLEC, LSC 세 가지 대표적인 스키마를 채택한다. 각 방법은 Schur complement S ≈ −B F⁻¹ Bᵀ 를 근사하는 방식이 다르다. PCD는 압력에 대한 convection‑diffusion 연산자를 별도로 구성하고, SIMPLEC은 대각선 근사와 압력‑속도 연계 항을 순차적으로 처리한다. LSC는 최소제곱(commutator) 접근을 통해 B F⁻¹ Bᵀ를 근사한다. 이때 블록 내부의 역연산 F⁻¹ 역시 두 레벨 OSM으로 근사한다.

실험 설정은 P2‑P1, Q2‑Q1, Q2‑P1‑disc. 등 다양한 혼합 유한요소 쌍과, 구조·비구조 격자, 그리고 후방 스텝과 동맥 모델이라는 현실적인 3차원 흐름 문제를 포함한다. Reynolds 수는 200‒3200, CFL 수는 0.5‒30 범위로 변동시켜 프리컨디셔너의 강인성을 평가한다. 결과는 단일 프리컨디셔너가 특히 고 Reynolds·고 CFL 상황에서 블록 방식보다 GMRES 반복 횟수가 현저히 적고, 전체 실행 시간이 짧으며, 스케일러빌리티가 우수함을 보여준다. 특히 GDSW* + RGDSW 조합은 압력‑속도 각각에 최적화된 coarse space를 제공함으로써, 전통적인 단일 GDSW만을 사용할 때보다 20‒30% 정도의 반복 감소 효과를 얻는다.

또한, 프리컨디셔너 설정 비용(예: coarse matrix 조립, 심볼릭 팩토리 재사용)과 실제 솔버 실행 비용을 정량적으로 비교한다. 단일 프리컨디셔너는 초기 설정 비용이 블록 방식보다 약 1.5배 정도 크지만, 반복 횟수 감소와 높은 병렬 효율 덕분에 전체 시간에서는 오히려 10‒40% 정도의 절감 효과를 나타낸다. 특히 코어 수가 1024‒4096 수준으로 확대될 때, 약 0.85‒0.9 수준의 병렬 효율을 유지한다.

이 논문은 두 레벨 OSM과 GDSW 계열 coarse space가 단일 및 블록 프리컨디셔너 모두에 적용될 수 있음을 입증하고, 특히 압력‑속도 변수에 서로 다른 coarse space를 적용하는 새로운 설계가 고성능 병렬 유동 시뮬레이션에 실질적인 이점을 제공한다는 점에서 중요한 기여를 한다.