[논문 리뷰] FreeFem++, a tool to solve PDEs numerically

Georges Sadaka|arXiv (Cornell University)|2012. 05. 07.

Advanced Numerical Methods in Computational Mathematics참고 문헌 1인용 수 6

한 줄 요약

FreeFem++는 변분 형식을 사용하여 편미분방정식(PDE)을 수치적으로 해결하기 위한 오픈소스이자 사용자 친화적인 유한요소 소프트웨어입니다. 2차원 및 3차원 문제를 지원하며, 이방성 적응을 갖춘 자동 메쉬 생성 기능과 다양한 유한요소(불연속 및 혼합 요소 포함)를 제공하여 선형 및 비선형 타원형, 쌍곡형, 포물형 PDE를 고수준 스크립팅 언어와 내장된 솔버(예: UMFPACK 및 GMRES)를 통해 효율적으로 해결할 수 있습니다.

ABSTRACT

FreeFem++ is an open source platform to solve partial differential equations numerically, based on finite element methods. It was developed at the Laboratoire Jacques-Louis Lions, Universit\\'e Pierre et Marie Curie, Paris by Fr\\'ed\\'eric Hecht in collaboration with Olivier Pironneau, Jacques Morice, Antoine Le Hyaric and Kohji Ohtsuka. The FreeFem++ platform has been developed to facilitate teaching and basic research through prototyping. FreeFem++ has an advanced automatic mesh generator, capable of a posteriori mesh adaptation; it has a general purpose elliptic solver interfaced with fast algorithms such as the multi-frontal method UMFPACK, SuperLU . Hyperbolic and parabolic problems are solved by iterative algorithms prescribed by the user with the high level language of FreeFem++. It has several triangular finite elements, including discontinuous elements. For the moment this platform is restricted to the numerical simulations of problems which admit a variational formulation. We will give in the sequel an introduction to FreeFem++ which include the basic of this software. You may find more information throw this link http://www.freefem.org/ff++.

연구 동기 및 목표

• 과학계 계산 분야에서 교육 및 프로토타이핑을 위해 특화된 접근하기 쉬운 오픈소스 도구를 제공하기 위해.
• 저수준 프로그래밍 없이도 연구자와 학생들이 빠르게 유한요소 방법을 구현하고 테스트할 수 있도록 하기 위해.
• 비선형 및 시간에 의존하는 문제를 포함한 복잡한 PDE 시스템을 변분 형식을 통해 지원하기 위해.
• 적응 메시 생성, 질량 행렬의 럼프화, 효율적인 희소 행렬 솔버와 같은 고급 수치 기능을 통합하기 위해.
• PDE 분야에서 모델 검증, 수학적 분석 및 알고리즘 개발을 위한 수치 실험을 촉진하기 위해.

제안 방법

• 수학적 표기와 유사하게 변분 형태로 PDE를 정의할 수 있는 고수준 스크립팅 언어를 사용합니다.
• 사후 이방성 적응을 갖춘 데라운레이-보르노이 기반 자동 메쉬 생성기를 활용합니다.
• 연속 및 불연속 라그랑주 요소(P1, P2, P1dc, P2dc)와 라비아르-톰슨 요소를 포함한 다양한 유한요소 유형을 지원합니다.
• UMFPACK, SuperLU, CG, GMRES와 같은 빠른 희소 행렬 솔버를 통합하여 선형 시스템을 효율적으로 해결합니다.
• 사용자 정의 반복 스킴을 통해 시간 적분을 가능하게 하며, 안정성 분석이 수행된 θ-스키마를 통해 포물형 문제를 해결합니다.
• gnuplot 및 기타 외부 도구를 통한 메쉬 간 보간, 오차 추정, 시각화 도구를 제공합니다.

실험 결과

연구 질문

• RQ1 고수준 오픈소스 플랫폼은 연구자와 교육자들이 PDE의 수치적 해법을 어떻게 간편화할 수 있나요?
• RQ2 FreeFem++를 사용하여 열 방정식에 θ-스키마를 적용할 경우, 공간과 시간에서의 수렴 특성은 어떻게 되나요?
• RQ3 자동 메쉬 적응 및 질량 럼프화 기법은 시간에 의존하는 PDE의 안정성과 효율성 향상에 얼마나 효과적인가요?
• RQ4 변분 형식을 통해 FreeFem++는 비선형성 및 결합된 항을 포함한 복잡한 PDE 시스템을 어느 정도 지원할 수 있나요?
• RQ5 Poisson 방정식 및 비선형 타원형 방정식과 같은 벤치마크 문제를 해결할 때 FreeFem++의 정확도와 견고성은 어떠한가요?

주요 결과

• θ = 0(명시적) 또는 θ = 1(역유클리드)일 경우, θ-스키마는 공간에서 2차 수렴, 시간에서 1차 수렴을 달성하며, θ = 0.5(Crank-Nicolson)일 경우 공간과 시간 모두에서 2차 수렴을 달성합니다.
• θ ∈ [0, 0.5)일 경우, θ-스키마는 Δt ≤ CFL·(Δx)²/(4μ·(1−2θ)) 조건을 만족할 때 안정합니다. 수치 실험 결과 Δx = Δy일 때 시간 수렴률은 약 1.0, 공간 수렴률은 약 2.0에 가까운 수렴 속도를 확인했습니다.
• qf1pTlump를 통해 질량 럼프화를 적용하면 질량 행렬이 대각선이 되어 쉽게 역행렬이 가능해지며, 이는 계산 효율성을 향상시키고 명시적 스킴의 CFL 조건을 연장합니다.
• 열 방정식에 대한 수치 결과는 공간에서 약 1.99992, 시간에서 약 0.99996의 수렴률을 보였으며, 이는 이론적 기대와 일치합니다.
• θ = 1/2일 경우, 공간 수렴률은 2.00195, 시간 수렴률은 2.00195를 기록하여 양방향에서 2차 정확도를 확보함을 보여줍니다.
• 메쉬 해상도 향상과 시간 간격 감소에 따라 L² 오차가 비례적으로 감소함을 통해 소프트웨어가 정확한 해를 고도로 재현함을 입증했습니다.