FEMPAR 입문: 병렬 유한요소 라이브러리 활용 가이드

초록

본 논문은 오픈소스 병렬 유한요소 프레임워크 FEMPAR v1.0.0의 사용법을 튜토리얼 중심으로 소개한다. 직렬 Poisson 문제, 적응형 옥트리 메쉬, 그리고 분산 메모리 환경에서의 BDDC 사전조건까지 세 단계로 구성된 예제를 통해 연속·비연속 Galerkin, h‑적응, 도메인 분해 솔버의 기본 흐름을 설명한다.

상세 분석

FEMPAR는 Fortran 2003/2008 표준을 기반으로 약 30만 라인의 객체지향 코드를 제공한다. 핵심 설계는 ‘triangulation’, ‘fe_space’, ‘fe_function’, ‘discrete_integration’, ‘fe_affine_operator’ 등 추상 타입을 정의하고, 이들을 다형적으로 확장할 수 있게 하는데 있다. 사용자는 weak form을 구현하는 discrete_integration 파생 클래스를 작성하고, 이를 fe_affine_operator에 연결해 전역 행렬·벡터를 자동 조립한다. 연속 Galerkin과 Interior Penalty DG 모두 동일한 인터페이스로 처리되며, 셀·면 적분은 참조 셀/면 위에서 고차 사분법을 이용해 수행된다. 메쉬는 구조화된 격자와 n‑cube, n‑simplex을 모두 지원하고, h‑적응을 위한 hanging node 처리와 p4est 기반의 병렬 옥트리 관리 기능을 제공한다. 특히, 대규모 병렬 실행 시 전역 자유도 번호 매김과 DOF 소유권 분배가 자동화되어 MPI 프로세스 간 통신을 최소화한다. 선형 시스템 해결에는 자체 구현된 BDDC(다중 레벨 균형 도메인 분해) 사전조건이 포함되어 있으며, 1.75 백만 MPI 태스크까지 확장된 사례가 보고된다. FEMPAR는 PETSc와 같은 외부 라이브러리와도 연동 가능하지만, 핵심 선형대수 모듈이 내장돼 있어 별도 의존성을 크게 줄인다. 튜토리얼 구조는 Fortran 프로그램 하나에 네 부분(모듈 임포트, 파라미터 선언, 메인 실행, 헬퍼 서브루틴)으로 일관성을 유지한다. Docker 컨테이너 기반 설치 가이드는 복잡한 의존성 문제를 회피하고, 사용자는 git 저장소에서 추가 튜토리얼을 손쉽게 내려받아 비선형, 시간 적분, 다중 물리 등 확장 기능을 학습할 수 있다. 전체적으로 FEMPAR는 고성능 병렬 유한요소 해석을 위한 완전한 스택을 제공하면서도, 초보 사용자를 위한 단계적 학습 경로를 마련한 점이 큰 강점이다.