대규모 유역 규모 물 흐름을 위한 오픈소스 병렬 리차드스 방정식 솔버

초록

RichardsFOAM은 OpenFOAM 기반의 오픈소스 병렬 솔버로, 리차드스 방정식을 대규모(수 km², 수십 년~수백 년) 시공간 해상도로 풀 수 있다. 1024 프로세서에서 약 90%의 병렬 효율을 달성했으며, 강·약 스케일링 분석을 통해 성능 유지 조건을 제시한다. 토양 오염 물질 이동, 토지 침하에 따른 수자원 변화, 유역 풍화 과정 등 다양한 환경·수리학 응용에 활용 가능하다.

상세 분석

RichardsFOAM은 비포화 토양 흐름을 기술하는 비선형 편미분 방정식인 Richards 방정식을 Finite‑Volume 방식으로 이산화하고, OpenFOAM의 기존 메쉬 분할·MPI 통신 구조를 그대로 활용한다. 비선형성 해소를 위해 기본적으로 Picard 반복을 사용하지만, 필요시 Newton‑Krylov 방법과 다중 수준 전처리기(ILU, GAMG 등)를 선택할 수 있다. 시간 적분은 암시적 백워드‑Euler와 가변 시간 단계 제어를 결합해, 급격한 습윤·건조 전이에서도 안정적인 수렴을 보장한다.

병렬 성능 측면에서 저자들은 도메인 분할에 따른 부하 균형을 강조한다. 메쉬가 균일하게 분포된 경우에는 단순한 정사각형 블록 분할만으로도 90% 이상의 효율을 얻지만, 실제 유역 모델에서는 토양 물성치가 공간적으로 크게 변동한다. 이를 보완하기 위해 메쉬 파티셔닝 단계에서 METIS/ParMETIS 기반의 그래프 분할을 적용하고, 각 서브도메인에 대한 로컬 메모리 사용량을 모니터링한다. 또한, I/O 병목을 최소화하기 위해 병렬 HDF5 포맷으로 결과를 기록하고, 체크포인트 파일을 비동기식으로 저장한다.

성능 테스트에서는 강 스케일링(고정 문제 크기에 프로세서 수 증가)과 약 스케일링(문제 크기를 프로세서 수에 비례해 확대) 모두에서 1024 코어까지 85~92%의 효율을 달성했다. 특히, 선형 시스템 해결 단계에서 GAMG 전처리기와 BiCGStab 솔버를 조합했을 때 가장 높은 스케일링을 보였으며, 전처리기 파라미터(스무딩 단계, 코어 수당 레벨 수) 조정이 전체 실행 시간에 미치는 영향을 정량적으로 분석하였다.

코드 공개는 GitHub와 CPC 프로그램 라이브러리를 통해 이루어졌으며, CMake 기반 빌드 스크립트와 상세한 사용자 매뉴얼이 제공된다. OpenFOAM의 모듈식 구조 덕분에 기존의 다상 흐름, 열전달, 화학 반응 모델과의 연계가 용이하며, 향후 기후·수문 모델링과의 통합도 기대된다. 이러한 기술적 구현과 성능 검증은 대규모 유역 수준에서 물 흐름을 메커니즘 기반으로 시뮬레이션하려는 연구자와 엔지니어에게 실질적인 도구적 가치를 제공한다.