심장혈관 CFD 시뮬레이션을 위한 오픈소스 솔버 lifex‑cfd 소개

초록

lifex‑cfd는 현대 C++와 deal.II 기반의 고성능 유한요소 라이브러리 lifex 위에 구축된 오픈소스 CFD 솔버로, 이동하는 심장판막·도메인·전이‑난류 흐름을 포함한 심혈관 혈류를 비압축성 Navier‑Stokes 방정식으로 모델링한다. ALE 프레임워크와 RIIS 판막 모델, 다양한 경계조건, SUPG/PSPG·VMS‑LES 안정화, MPI 기반 분산 메모리 병렬화를 제공하며, 수렴 검증과 벨라미 흐름 벤치마크에서 거의 이상적인 스케일링을 보인다. LGPLv3 라이선스로 공개되어 재현성·협업을 촉진한다.

상세 분석

lifex‑cfd는 심혈관 CFD에 특화된 기능을 통합한 종합 솔버로, 기존 오픈소스 패키지와 비교했을 때 몇 가지 핵심 차별점을 가진다. 첫째, ALE(Arbitrary Lagrangian‑Eulerian) 프레임워크를 통해 심실 수축·확장과 같은 대규모 도메인 변형을 정확히 추적한다. 도메인 변위는 조화 확장 또는 가상의 선형 탄성 모델을 이용해 풀며, 사용자는 강성 계수 α를 영역별로 조정해 메쉬 품질을 유지할 수 있다. 둘째, 심장판막은 RIIS(Resistive Immersed Implicit Surface) 방법으로 구현된다. 판막 표면을 거리 함수 φ로 정의하고, 얇은 층 ε 안에서 저항 R을 적용해 유체와 판막 사이의 무슬립 조건을 약하게 강제한다. 이는 복잡한 판막 개폐 동작을 별도의 구조해석 없이도 효율적으로 시뮬레이션하게 한다. 셋째, 수치 해석 측면에서 선형 및 비선형 시스템을 고차원 유한요소(FE) 공간에 Galerkin 방식으로 이산화하고, SUPG/PSPG, VMS‑LES 등 다양한 안정화·난류 모델을 선택 가능하게 한다. 특히 VMS‑LES는 저해상도 격자에서도 전이‑난류 현상을 포착한다. 선형 솔버는 Trilinos, PETSc, deal.II 자체 백엔드를 지원하며, GMRES, FGMRES, CG 등 Krylov 방법과 AMG, BoomerAMG, Additive Schwarz 등 사전조건자를 조합해 대규모 문제에서도 빠른 수렴을 보인다. 넷째, 병렬 구현은 MPI 기반으로, 메모리 분산과 매트릭스‑프리 연산을 활용해 수천 코어까지 거의 선형 속도 향상을 달성한다. 논문에서는 Beltrami 흐름 벤치마크와 다양한 이상·환자‑특이 모델을 통해 수렴 차수와 병렬 효율을 정량적으로 검증하였다. 마지막으로, lifex‑cfd는 파라미터 파일 기반의 사용자 친화적 인터페이스와 체크포인트·재시작 기능을 제공해 비전문가도 손쉽게 시뮬레이션을 설정·운영할 수 있다. LGPLv3 라이선스로 배포되어 코드 접근성·투명성을 보장하고, lifex 생태계 내 전기생리·기계·FSI·전기‑기계‑유체 상호작용 모듈과 원활히 연동될 수 있다. 이러한 설계 철학은 심혈관 연구·임상 전이 단계에서 재현 가능하고 확장 가능한 시뮬레이션 파이프라인을 구축하는 데 큰 장점을 제공한다.