오픈소스 고성능 몰입형 유체‑구조 상호작용 프레임워크

초록

**
본 논문은 MFEM과 FEBio 두 개의 검증된 유한요소 라이브러리를 결합한 오픈소스 몰입형 FSI(Fluid‑Structure Interaction) 프레임워크를 제시한다. 가상 도메인·분산 라그랑주 승수 방식과 변분 다중스케일(VMS) 안정화 기법을 이용해 큰 변형과 접촉을 포함한 심장 판막 시뮬레이션을 효율적으로 수행한다. GPU 가속 및 대규모 분산 메모리 병렬성을 지원한다.

**

상세 분석

**
이 연구는 기존 ALE 기반 FSI가 겪는 격자 왜곡·재메시징 문제를 회피하기 위해 가상 도메인(Fictitious Domain, FD) 접근법을 채택한다. FD 방식은 유체 방정식을 전체 배경 영역에 확장하고, 고체는 동일 영역에 삽입함으로써 인터페이스 메쉬를 별도로 관리할 필요가 없으며, 분산 라그랑주 승수(Distributed Lagrange Multiplier)로 속도와 트랙션 연속성을 강제한다. 이러한 구조는 대변형 및 접촉 현상이 빈번히 발생하는 심장 판막과 같은 생체 시스템에 적합하다.

수치적 안정성을 높이기 위해 변분 다중스케일(VMS) 안정화 기법을 적용한다. VMS는 저해상도 격자에서도 압력·속도 경계층을 정확히 포착하도록 돕고, 인공 유체 영역에서 발생할 수 있는 수치적 오차를 최소화한다. 또한, 완전 암시적(monolithic) 스키마를 사용해 유체와 고체의 비선형 연동을 하나의 비선형 시스템으로 동시에 해결한다. 이는 강하게 결합된 FSI 문제에서 수렴성을 크게 향상시키며, 시간 적분에 있어 큰 시간 단계도 안정적으로 사용할 수 있게 한다.

프레임워크의 핵심은 두 개의 오픈소스 라이브러리인 MFEM과 FEBio의 전략적 결합이다. MFEM은 고성능 GPU 및 대규모 분산 메모리 병렬화를 지원하는 최신 유한요소 툴킷으로, 특히 대규모 유체 해석에 최적화되어 있다. 반면 FEBio는 생체 조직을 위한 고급 비선형 고체 역학 모델(하이퍼엘라스틱, 비스코엘라스틱, 접촉 알고리즘 등)을 제공한다. 두 라이브러리를 C++ 레벨에서 인터페이스하고, 데이터 구조와 솔버 스케줄링을 통합함으로써, 사용자는 복잡한 고체 모델을 그대로 활용하면서도 MFEM의 스케일러블 유체 솔버를 이용할 수 있다.

구현 측면에서는 모듈형 아키텍처를 채택해 새로운 물리량(예: 전기·화학·성장 모델)이나 요소 유형을 손쉽게 추가할 수 있다. 또한, MPI 기반의 분산 실행과 CUDA/OpenCL을 통한 GPU 가속을 지원해, 3D 심장 판막과 같은 고해상도 시뮬레이션도 실용적인 시간 안에 수행 가능하도록 설계되었다.

검증 사례로는 전통적인 베눗-샤프(benchmark) 문제, 유동-구조 결합 테스트, 그리고 실제 3D 반월판(semilunar) 심장 판막 시뮬레이션이 제시된다. 특히 판막 사례에서는 고체의 복잡한 접촉·마찰 모델과 비선형 조직 강성을 그대로 적용했으며, 유동 영역에서는 혈류의 비정상적인 압력 파동과 와류를 정확히 재현했다. 성능 평가에서는 CPU 대비 GPU 가속 시 5~10배, 노드 수 증가에 따른 거의 선형 확장성을 확인하였다.

전반적으로 이 프레임워크는 고성능 컴퓨팅 환경에서 대변형·접촉을 포함한 생체 FSI 문제를 해결하기 위한 가장 포괄적인 오픈소스 솔루션 중 하나로 자리매김한다.

**