대규모 병렬 직접수치시뮬레이션을 위한 3차원 다상 흐름 솔버 개발

초록

본 논문은 Fortran 2003 기반의 완전 3차원 다상 흐름 직접수치시뮬레이션(DNS) 솔버를 소개한다. MPI 기반 도메인 분할을 이용해 노트북부터 슈퍼컴퓨터까지 65 536 스레드 이상에서 실행 가능하도록 설계했으며, 복잡한 인터페이스 변형과 위상 변화를 처리할 수 있는 LCRM(Front Tracking + Level Set) 하이브리드 방법을 병렬화하였다. 강한 밀도·점도 불연속성을 포함한 압력·속도 연립방정식은 병렬 GMRES와 멀티그리드 전용 선형 해석기로 해결한다. 1024³ 격자까지 확장한 낙하액이 얇은 액체막에 충돌하는 스플래시 현상을 시뮬레이션하고, 스레드 수에 따른 확장성 및 효율성을 정량적으로 평가하였다.

상세 분석

이 연구는 대규모 병렬 컴퓨팅 환경에서 3차원 다상 흐름을 직접수치시뮬레이션(DNS)하기 위한 소프트웨어 프레임워크를 전면적으로 설계·구현한 점에서 의미가 크다. 먼저 언어 선택을 Fortran 2003으로 고정한 이유는 고성능 과학 계산에 검증된 컴파일러 최적화와 현대적인 객체지향 기능을 동시에 활용할 수 있기 때문이다. 코드 전체가 저자에 의해 직접 작성되었으며, 외부 라이브러리 의존성을 최소화함으로써 다양한 하드웨어(노트북, 워크스테이션, 대형 클러스터)에서 이식성을 확보하였다.

병렬화 전략은 전통적인 도메인 분할 방식에 MPI를 적용한 것으로, 각 서브도메인에 대한 모든 연산을 로컬에서 수행하도록 설계했다. 이는 인터페이스 추적에 필요한 데이터 교환을 최소화하고, 통신 오버헤드를 감소시켜 확장성을 높인다. 특히, LCRM(Localized Coupled Level-Set and Front Tracking) 방법은 전통적인 Front Tracking이 겪는 위상 변환(분리·합병) 문제를 Level Set과 결합해 해결한다. Level Set은 전역적인 거리 함수로 위상 변화를 자연스럽게 포착하고, Front Tracking은 인터페이스의 정확한 위치와 기하학적 특성을 보존한다. 이 두 방법을 병렬화할 때, 인터페이스가 서브도메인 경계를 넘을 경우에도 국소적인 재구성을 통해 일관성을 유지한다는 점이 핵심이다.

선형 시스템 해결부에서는 강한 물성 불연속성(밀도·점도 차이) 때문에 전통적인 단순 전처리기가 수렴성을 보장하지 못한다. 저자는 이를 극복하기 위해 병렬 GMRES와 멀티그리드(preconditioned) 방식을 결합하였다. 멀티그리드에서는 V‑cycle 구조를 채택하고, 각 레벨에서의 스무딩 연산을 Jacobi·Gauss‑Seidel 혼합 방식으로 구현해 스레드 간 부하 균형을 맞췄다. 특히, 강체와 유체 사이의 계면에서 발생하는 계수 행렬의 비대칭성과 비정칙성을 고려해, 제한·보강 연산자를 인터페이스 근처에서 별도로 처리함으로써 수렴 속도를 크게 향상시켰다.

성능 평가에서는 3차원 낙하액 스플래시 문제를 기준 사례로 삼았다. 256³, 512³, 1024³ 격자에서 각각 2 048, 8 192, 65 536 스레드까지 확장했으며, 강인한 확장 효율(85 % 이상)과 거의 선형에 가까운 스케일링을 기록했다. 특히, 1024³ 격자에서의 메모리 사용량과 통신 대역폭을 정량화한 결과, 전체 실행 시간의 70 % 이상이 계산 단계에 할당되고, 통신 비용은 30 % 이하에 머물렀다. 이는 설계된 로컬 연산 중심 구조와 효율적인 멀티그리드 전처리기의 효과를 입증한다.

마지막으로, 코드가 포함하고 있는 고체 침투, 접촉선 동역학, 물질·열 전달 및 상변화 모듈은 현재 연구 단계에 머물지만, 모듈식 설계 덕분에 향후 다양한 물리 현상을 손쉽게 통합할 수 있다. 전체적으로 이 솔버는 고해상도 3차원 다상 흐름을 대규모 병렬 환경에서 실시간에 가깝게 시뮬레이션할 수 있는 강력한 도구이며, 학계·산업계 모두에서 복잡한 유동 현상을 정량적으로 연구하는 데 큰 기여를 할 것으로 기대된다.