[논문 리뷰] An implementation of a locally modified finite element method for interface problems in deal.II

Stefan Frei, Thomas Richter외 1명|arXiv (Cornell University)|2018. 10. 11.

Advanced Numerical Methods in Computational Mathematics참고 문헌 31인용 수 3

한 줄 요약

이 논문은 메esh 피팅을 피하기 위해 국소적인 자유도를 고정된 사각형 패치 내에서 매개변수적으로 통합함으로써, 약한 불연속성을 갖는 인터페이스 문제를 해결하는 국소적으로 수정된 유한요소법을 deal.II에서 제시한다. 이 방법은 최적 수렴성을 달성하며, CG/PCG와 같은 반복적 해법을 통해 효율적인 해를 도출할 수 있으며, 두 가지 수치 실험과 조건수 분석을 통해 성능이 검증되었다.

ABSTRACT

This source code contains an implementation of the locally modified finite element method first proposed by S. Frei, T. Richter, SIAM Journal on Numerical Analysis 52 (5), 2315-2334. This method is a simple approach to solve interface problems, as arising in fluid-structure-interactions or multi-phase flows, or problems on complex boundaries. The implementation is contained in the files locmodfe.cc/h, together with a user program (step-modfe.cc) that is similar to many deal.II steps (www.dealii.org) and a problem-specific file (problem.h) with specific data for some numerical tests provided in the src folder. The source code is based on the deal.II.8.5.0 version.<br> In the test folder, we provide the parameter files parameters_test_case_1/2.prm, where algorithmic, numerical and physical parameters can be changed. In the current version, parameters are pre-set for the two numerical tests. For the meaning of the parameters, see the detailed comments in the parameter files.<br> <br> Moreover, we provide the output that the program writes to a terminal in the files test_case_1/2.dlog in the same folder, for a comparison of the terminal output results. In the build folder, the file CMakeLists.txt can be used to create a Makefile using cmake, see the explanations given in README.txt. Moreover, the README.txt file provides detailed installation instructions and a compact overview of all source files. Algorithmic descriptions and details of the delivered source code are provided in S. Frei, T. Richter, T. Wick; 2018, arXiv preprint arXiv:1806.00999

연구 동기 및 목표

• 불연속적인 계수를 갖는 타원형 인터페이스 문제를 위한 강력하고 효율적인 유한요소법을 개발한다.
• 메쉬 재구성 없이도 이동하거나 곡선형 인터페이스를 정확하게 해석할 수 있도록 국소적 매개변수적 강화를 제공한다.
• 표준 코드에 최소한의 수정으로 deal.II 유한요소 라이브러리 내에 이 방법을 구현한다.
• 해결된 선형 행렬에 대해 반복적 선형 해법(CG/PCG)의 효과성을 입증한다.
• 메쉬 정밀도 향상과 이방성 조건 하에서의 조건수 및 성능을 분석한다.

제안 방법

• 사각형 패치로 구성된 고정 배경 메쉬를 사용하며, 인터페이스에 의해 잘린 요소들에 대해서만 국소적으로 8개의 삼각형 또는 4개의 사각형으로 정밀화한다.
• 기본 요소와 물리적 패치 사이의 매개변수적 사상에 의해 국소 자유도를 통합함으로써 명시적 메쉬 생성을 피한다.
• 해당 방법은 이산 수준의 레벨셋 함수를 사용하여 부분도메인 소속을 결정하고, 적분 점에서 변수 확산 계수를 적용한다.
• 계수 수치의 안정성과 조건수의 경계를 확보하기 위해 계층 기반을 구성함으로써 시스템을 안정화시킨다.
• 표준 deal.II 인프라를 재사용하며, 유한요소 값과 적분 평가에만 국소적 수정을 가한다.
• 동일 차수의 유한요소를 지원하며 압력 안정화 기법을 적용할 수 있으며, 고차원 및 3차원 문제로의 확장도 가능하다.

실험 결과

연구 질문

• RQ1 메쉬 재구성 없이도 국소적으로 수정된 유한요소법이 인터페이스 문제에 대해 최적 수렴성을 달성할 수 있는가?
• RQ2 메쉬 정밀도 향상과 인터페이스 위치 변화 조건 하에서 시스템 행렬의 조건수가 어떻게 변화하는가?
• RQ3 CG 및 PCG와 같은 반복적 해법이 결과 선형 시스템에 효과적으로 적용될 수 있는가?
• RQ4 고대비 계수를 갖는 문제에서 정확도 및 수렴 속도 측면에서 이 방법의 성능은 어떠한가?
• RQ5 구현 복잡성과 안정성 측면에서 이 매개변수적 강화 전략은 기존의 강화된 FEM과 비교해 어떻게 다른가?

주요 결과

• 수치 실험을 통해, 점프 계수를 갖는 타원형 인터페이스 문제에서 최적 수렴 속도를 달성함을 확인하였다.
• 스케일링된 계층 기반을 사용할 경우 시스템 행렬의 조건수가 유한하게 유지되어, 임계 인터페이스 위치 근처에서의 악조건화를 방지한다.
• CG 및 PCG와 같은 반복적 해법은 효과적이며 확장 가능하며, 메쉬 정밀도 향상과 이방성 조건 하에서도 수렴 행동이 안정된다.
• 표준 deal.II 코드에 대한 수정이 거의 필요하며, 모든 수정 사항이 패치 수준의 유한요소 평가에 국한되어 있다.
• 이 방법은 고차원 및 3차원 문제로의 확장이 가능하며, 병렬 처리 및 경계 보정 기법에 대한 후속 연구가 진행 중이다.
• 이전의 응용 사례에서 보여졌듯이, 큰 변형을 수반하는 유체-구조 상호작용과 같은 복잡한 다물리학 문제에 적합하다.