[논문 리뷰] MultiMesh Finite Element Methods I: Solving PDEs on Multiple Intersecting Meshes

August Johansson, Benjamin Kehlet외 2명|arXiv (Cornell University)|2018. 03. 23.

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

한 줄 요약

이 논문은 다중 교차 메시 위에서 PDE를 해결하기 위한 새로운 프레임워크인 MultiMesh 유한요소법(MultiM-FEM)을 소개한다. 이는 엔진 부품이나 다중물리계 시스템과 같은 별개의 영역을 독립적으로 메시화할 수 있도록 하여, 큰 상대 운동을 겪는 구성 요소들에 대한 해를 구하는 데 기여한다. 최적 수렴성, 얇은 교차 영역 및 반올림 오차에 대한 강건성, 그리고 큰 변형을 겪는 문제에의 적용 가능성을 입증한다.

ABSTRACT

We present a new framework for expressing finite element methods on multiple intersecting meshes: multimesh finite element methods. The framework enables the use of separate meshes to discretize parts of a computational domain that are naturally separate; such as the components of an engine, the domains of a multiphysics problem, or solid bodies interacting under the influence of forces from surrounding fluids or other physical fields. Such multimesh finite element methods are particularly well suited to problems in which the computational domain undergoes large deformations as a result of the relative motion of the separate components of a multi-body system. In the present paper, we formulate the multimesh finite element method for the Poisson equation. Numerical examples demonstrate the optimal order convergence, the numerical robustness of the formulation and implementation in the face of thin intersections and rounding errors, as well as the applicability of the methodology. In the accompanying paper~\cite{mmfem-2}, we analyze the proposed method and prove optimal order convergence and stability.

연구 동기 및 목표

• 서로 다른 메시로 독립적으로 메시화된 구성 요소들로 이루어진 계산 영역에서 PDE를 해결하는 데 도전하는 문제를 해결한다.
• 유한요소법이 자연스럽게 분리된 영역, 예를 들어 상호작용하는 고체, 유체 또는 다중물리계 구성 요소를 처리할 수 있도록 한다.
• 매우 얇거나 왜곡된 메시 교차 영역 및 수치 반올림 오차가 존재하는 상황에서도 안정적이고 정확한 수치 프레임워크를 개발한다.
• 각 영역에 대해 별도의 메시를 사용하는 시간에 따라 변화하는 문제, 특히 복잡한 기하학적 상호작용이 있는 문제를 해결하기 위한 기반을 제공한다.

제안 방법

• 각 메시가 독립적으로 정밀화되고 적응 가능한 다중 겹치는 메시를 사용하여 페르소의 약한 형태를 수립한다.
• 메시 교차 영역에서의 연속성 또는 경계 조건을 약하게 강제하기 위해 니츠체의 방법을 적용한다.
• 배경 메시와 전경 메시의 교차 영역을 통합하기 위한 적분 규칙을 사용하여 겹치는 요소에서 정확한 적분을 보장한다.
• 메시의 기하학적 복잡성을 처리하기 위해 몸체에 맞는 메시가 필요 없이 컷 유한요소법 접근법을 적용한다.
• 메시 간의 일치하지 않는 경계에서 해를 결합하기 위해 보존적인 모르타르 유사 접근법을 구현한다.
• 작은 또는 왜곡된 교차 영역이 존재하는 상황에서도 안정성을 확보하기 위해 안정화된 변분 형태를 사용한다.

실험 결과

연구 질문

• RQ1 복잡한 방식으로 교차하는 다수의 독립적으로 생성된 메시 위에서 PDE를 해결할 수 있는 유한요소법을 제안할 수 있는가?
• RQ2 특히 교차 영역이 매우 얇거나 왜곡된 경우, 메시 간의 연속성과 안정성을 어떻게 확보할 수 있는가?
• RQ3 기하학적 및 적분 문제에 직면해도, 페르소 방정식에 대해 제안된 방법이 최적 수렴률을 달성하는가?
• RQ4 수치 반올림 오차 및 작은 교차 영역 부피에 대해 이 방법은 얼마나 강건한가?
• RQ5 이 프레임워크는 구성 요소 간에 큰 상대 운동을 겪는 시간에 따라 변화하는 문제로 확장 가능한가?

주요 결과

• 제안된 MultiMesh FEM는 수치 실험을 통해 페르소 방정식에 대해 최적 순서 수렴성을 확보한다.
• 메시 간 교차 영역이 매우 얇거나 잘 맞지 않는 경우에도 수치적으로 강건성을 유지한다.
• 반올림 오차가 존재하는 상황에서도 안정성과 정확성이 유지되어 실용적 신뢰성을 입증한다.
• 별개의 영역에 대해 독립적인 메시화를 지원하여, 복잡한 다체계를 효율적으로 처리할 수 있다.
• 메시 재구성이나 복잡한 경계 추적 기법이 필요 없기 때문에, 큰 변형을 겪는 문제에 매우 적합하다.
• 이론적 분석을 통해 후속 논문에서 제안된 바와 같이, 이 구현은 다른 PDE 및 다중물리 문제로도 확장 가능하고 스케일링이 가능하다.