[논문 리뷰] Energy stable and accurate coupling of finite element methods and finite difference methods

Tuan Anh Dao, Ken Mattsson외 1명|arXiv (Cornell University)|2021. 03. 12.

Advanced Numerical Methods in Computational Mathematics참고 문헌 37인용 수 13

한 줄 요약

이 논문은 연속 갈레르킨 유한요소(FE) 및 고차 유한차분(FD) 방법 간의 증명 가능하게 안정적이고 정확하며 에너지 보존 성질을 갖는 결합 방법을 제안한다. 이는 합성-통합-부분(SBP) 연산자와 동시에 approximation 항(SAT)을 사용한다. 이 기법은 최소한의 수정으로 비일치하는 다블록 결합을 가능하게 하며, 비대칭 노름 SBP 유지 보존 보간 연산자를 사용해 경계 연속성을 약하게 강제함으로써, 비선형 보존 법칙과 질량 럼프링 없이도 연속 FE 질량 행렬을 갖는 경우에도 엄격한 안정성과 보존 성질을 유지한다.

ABSTRACT

We introduce a hybrid method to couple continuous Galerkin finite element methods and high-order finite difference methods in a nonconforming multiblock fashion. The aim is to optimize computational efficiency when complex geometries are present. The proposed coupling technique requires minimal changes in the existing schemes while maintaining strict stability, accuracy, and energy conservation. Results are demonstrated on linear and nonlinear scalar conservation laws in two spatial dimensions.

연구 동기 및 목표

• 복잡한 기하구조를 위한 연속 갈레르킨 유한요소(FE) 및 유한차분(FD) 방법 간의 안정적이고 정확하며 에너지 보존 성질을 갖는 하이브리드 방법을 개발하는 것.
• 연속 갈레르킨 FE 방법에서 발생하는 비대칭 질량 행렬로 인해 FE 방법과 고차 FD 스킴을 결합하는 데 발생하는 과제를 해결하는 것.
• 기존 스킴에 대한 질량 럼프링이나 상당한 수정 없이도 엄격한 안정성과 보존 성질을 유지하는 것.
• 기존 FD 및 FE 구현에 최소한의 변경으로 비일치하는 다블록 결합을 가능하게 하는 것.
• SBP-SAT 프레임워크를 FE-FD 결합으로 확장하여 비선형 스칼라 보존 법칙을 포함하는 것.

제안 방법

• SBP-SAT 프레임워크를 사용해 결합 시스템에서 엄격한 안정성과 에너지 보존 성질을 보장한다.
• SBP 유지 보존 보간 연산자를 사용한 동시에 approximation 항(SAT)을 통해 경계 연속성을 약하게 강제한다.
• FE-FD 결합을 위한 비대칭 노름 보간 연산자를 구성하여 질량 럼프링이 필요 없도록 한다.
• 두 가지 구축 기법을 도입: 하나는 직접 행렬 구조 해법에 기반하고, 다른 하나는 기존의 비대칭 노름 보간 방법을 활용할 수 있도록 하는 대칭 노름 '브릿지' 층을 사용한다.
• SAT 수식을 FD 및 FE 측 모두에 일관되게 적용하며, 결합 스킴이나 그리드에 독립적인 일반적인 FE 수식을 제공한다.
• 구조적 격자에서는 고차 SBP 유한차분 연산자를, 비구조적 또는 비일치 격자에서는 연속 갈레르킨 유한요소를 사용한다.

실험 결과

연구 질문

• RQ1 비일치하는 다블록 설정에서 고차 유한차분과 연속 갈레르킨 유한요소 방법 간에 안정적이고 정확한 결합이 달성될 수 있는가?
• RQ2 FE 방법에서 유도되는 비대칭 질량 행렬을 사용할 경우, 경계 연속성을 약하게 강제하는 방식으로 엄격한 안정성과 에너지 보존 성질을 유지할 수 있는가?
• RQ3 질량 럼프링이나 대칭 노름 근사가 필요 없이 SBP-SAT 프레임워크를 FE-FD 결합으로 확장할 수 있는가?
• RQ4 FE 및 FD 블록 간의 비일치 경계에서 정확성과 안정성을 유지하기 위해 필요한 보간 연산자는 무엇인가?
• RQ5 제안된 결합 방법은 선형 및 비선형 보존 법칙, 특히 이동-확산 방정식과 버거스 방정식에 대해 어떻게 성능을 발휘하는가?

주요 결과

• 제안된 결합 방법은 선형 및 비선형 문제 모두에서 2차 수렴을 달성하며, 고차 FD 스킴은 정확성과 수렴률을 향상시킨다.
• 수치 결과에 따르면, 동일한 설정에서 4차 FD 스킴은 2차 FD 스킴 대비 l2 오차를 약 50% 감소시킨다.
• 비일치 경계를 통해 선형 이동-확산 방정식과 비선형 버거스 방정식 모두에서 엄격한 에너지 안정성과 보존 성질을 유지한다.
• 비대칭 FE 질량 행렬을 사용하더라도 질량 럼프링 없이도 안정성과 정확성이 손상되지 않으며, 이는 방법의 강건성을 입증한다.
• 직접 행렬 해법과 브릿지 층 방법이라는 두 가지 보간 구축 기법 모두 안정적이고 정확한 결과를 도출하며, 후자는 기존의 대칭 노름 보간 도구 재사용을 가능하게 한다.
• 이 방법은 일반적인 스칼라 보존 법칙에 적용 가능하며, 비선형 시스템, 곡선 격자, 고차 유한요소로의 확장도 가능하다.