[논문 리뷰] A Numerical Study on the Weak Galerkin Method for the Helmholtz Equation with Large Wave Numbers

Lin Mu, Junping Wang외 2명|arXiv (Cornell University)|2011. 11. 02.

Advanced Numerical Methods in Computational Mathematics참고 문헌 40인용 수 15

한 줄 요약

이 논문은 큰 파장 수를 가진 헬름홀츠 방정식을 해결하기 위해 약한 갈레르킨(WG) 유한요소 방법을 조사하며, 이 방법의 강건성, 정확도 및 구현 용이성을 입증한다. 이 방법은 연속성이 요구되지 않는 다항식 근사와 약한 형태로 정의된 기울기를 사용하여 고파장 수에서도 안정적인 수렴을 달성한다. 상대적 $ H^1 $ 오차는 $ kh = 0.25 $일 때 증가가 최소한이 되어 오염 오차 통제가 효과적임을 나타낸다.

ABSTRACT

Weak Galerkin (WG) refers to general finite element methods for partial differential equations in which differential operators are approximated by weak forms through the usual integration by parts. In particular, WG methods allow the use of discontinuous finite element functions in the algorithm design. One of such examples was recently introduced by Wang and Ye for solving second order elliptic problems. The goal of this paper is to apply the WG method of Wang and Ye to the Helmholtz equation with high wave numbers. Several test scenarios are designed for a numerical investigation on the accuracy, convergence, and robustness of the WG method in both inhomogeneous and homogeneous media over convex and non-convex domains. Our numerical experiments indicate that weak Galerkin is a finite element technique that is easy to implement, and provides very accurate and robust numerical solutions for the Helmholtz problem with high wave numbers.

연구 동기 및 목표

• 큰 파장 수를 가진 헬름홀츠 방정식을 해결하기 위해 약한 갈레르킨(WG) 유한요소 방법의 성능을 평가하는 것.
• 여러 구형 및 비구형 영역에서 균질성 및 비균질성 매질에서 WG 방법의 정확도, 수렴성 및 강건성을 평가하는 것.
• 특히 매개변수 민감도 및 오염 오차 측면에서 기존의 불연속 갈레르킨 방법과의 비교를 통해 WG 성능을 평가하는 것.
• 메쉬 세분화 및 다항식 차수의 영향을 고려하여 고파장 수 문제에서 해의 정확도에 미치는 영향을 조사하는 것.

제안 방법

• WG 방법은 해를 위한 불연속 조각다항식 함수를 사용하며, 부분적 통합을 통해 기울기 연산자의 약한 형태를 적용한다.
• 각 요소에 대해 이산적 약한 기울기를 구성함으로써 요소 경계를 넘어서도 연속성이 요구되지 않는 국소적 미분 연산자 근사를 가능하게 한다.
• 변분 형태는 각 요소에서 잔여 항의 의미로 헬름홀츠 방정식을 강제하며, 소머펠드 방사 조건을 근사하기 위해 로빈형 흡수 경계 조건을 사용한다.
• 표준 불연속 갈레르킨 방법과 달리 이 방법은 페널티 매개변수를 필요로 하지 않아 구현이 단순해지고 매개변수 조정이 감소한다.
• 수치 실험은 조각별 상수 및 조각별 선형 요소를 사용하여 다양한 파장 수와 메쉬 크기에서 수렴성을 시험한다.
• 상대적 $ H^1 $ 오차는 수렴 속도와 오염 효과 평가를 위해 계산되며, 파장 수의 강건성을 평가하기 위해 $ kh $를 일정하게 유지한다.

실험 결과

연구 질문

• RQ1 약한 갈레르킨 방법은 큰 파장 수에서 헬름홀츠 방정식에 대해 정확도와 수렴성을 유지하는가?
• RQ2 페널티가 포함된 불연속 갈레르킨 방법과 비교할 때 WG 방법의 오차 및 매개변수 의존성은 어떠한가?
• RQ3 WG 방법은 고파장 수 문제에서 오염 효과를 어느 정도 완화하는가?
• RQ4 메쉬 세분화 및 다항식 차수는 고주기적 헬름홀츠 해의 정확도에 어떤 영향을 미치는가?

주요 결과

• WG 방법은 파장 수가 $ k = 100 $에 이르는 범위까지 안정적인 수렴을 달성하며, $ kh = 0.25 $일 때 상대적 $ H^1 $ 오차가 낮게 유지되어 오염 효과가 최소임을 나타낸다.
• $ kh = 0.25 $일 때 수렴선은 거의 수평이 되며 약간의 기울기만 존재하여 파장 수 증가에 대해 강건함을 시사한다.
• $ k = 100 $ 및 $ h = 1/200 $에서 정확한 해와의 일치가 뛰어나며, 메쉬 세분화가 진행될수록 시각적 수렴성이 향상된다.
• $ k = 100 $ 및 $ h = 1/60 $에서 조각별 선형 요소를 사용할 경우 WG 방법은 해의 진폭과 고주기적 진동을 정확히 포착한다.
• 이 방법은 페널티 매개변수를 필요로 하지 않아 표준 DG 방법보다 구현이 단순하면서도 유사한 정확도를 달성한다.
• 심지어 $ k = 240 $일 때도 $ kh $를 일정하게 유지하면 오차는 $ k $ 증가에 따라 단조적으로 증가하지만 기울기는 작아 오염 오차가 통제됨을 나타낸다.