[논문 리뷰] Continuous Interior Penalty Finite Element Methods for the Helmholtz Equation with Large Wave Number

Haijun Wu|arXiv (Cornell University)|2011. 06. 21.

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

한 줄 요약

이 논문은 큰 파장 수를 가진 헬름홀츠 방정식에 대해 복소 수치의 양의 허수부를 가진 복소 벌점 매개수를 사용하는 새로운 연속 내부 벌점 유한요소법(CIP-FEM)을 제안한다. 이 방법은 조건부 안정성과 최적의 사전점근 오차 추정치를 달성하며, 표준 FEM에 비해 오염 오차를 크게 감소시킨다. 오차 한계는 $ C_1kh + C_2k^3h^2 $이며, 이는 $ k^3h^2 \leq C_0 $일 때이고, 그 외의 경우에는 $ C_1kh + \frac{C_2}{\gamma} $이다.

ABSTRACT

This paper develops and analyzes some continuous interior penalty finite element methods (CIP-FEMs) using piecewise linear polynomials for the Helmholtz equation with the first order absorbing boundary condition in two and three dimensions. The novelty of the proposed methods is to use complex penalty parameters with positive imaginary parts. It is proved that, if the penalty parameter is a pure imaginary number $ı\ga$ with $0 C_0$ and $kh$ is bounded, where $k$ is the wave number, $h$ is the mesh size, and the $C$'s are positive constants independent of $k$, $h$, and $\ga$. Optimal order $L^2$ error estimates are also derived. The analysis is also applied if the penalty parameter is a complex number with positive imaginary part. By taking $\ga o 0+$, the above estimates are extended to the linear finite element method under the condition $k^3h^2\le C_0$. Numerical results are provided to verify the theoretical findings. It is shown that the penalty parameters may be tuned to greatly reduce the pollution errors.

연구 동기 및 목표

• 큰 파장 수를 가진 헬름홀츠 문제에 대한 유한요소법에서 지속적인 오염 오차 문제를 해결한다.
• 메쉬 제약 없이 고파장 수 문제에 대해 안정적이고 잘 정의된 유한요소법을 개발한다.
• 복소 벌점 매개수를 도입함으로써 2차원 및 3차원 선형 유한요소법의 사전점근 오차 추정치를 향상시킨다.
• 벌점 매개수의 조정이 오염 오차를 크게 감소시켜 수치 정확도를 향상시킴을 입증한다.
• 벌점 매개수의 허수부를 0으로 취한 극한에서 표준 FEM의 이론적 결과를 회복한다.

제안 방법

• 첫 번째 흡수 경계 조건을 가진 헬름홀츠 방정식에 대해 일차 다항식 요소를 사용하는 연속 내부 벌점 유한요소법(CIP-FEM)을 도입한다.
• 모든 $ \gamma > 0 $에 대해 $ \mathbf{i}\gamma $ 형태의 복소 벌점 매개수를 사용하여 무조건적 안정성과 잘 정의됨을 보장한다.
• 일반적인 $ H^1 $-노름 오차 추정치를 유도한다: $ k^3h^2 \leq C_0 $일 때는 $ C_1kh + C_2k^3h^2 $, $ k^3h^2 > C_0 $일 때는 $ C_1kh + \frac{C_2}{\gamma} $이며, 상수들은 $ k $, $ h $, $ \gamma $에 독립적이다.
• 동일한 조건 하에서 최적의 $ L^2 $-노름 오차 추정치를 확립한다.
• 일반적인 복소 벌점 매개수를 적용하여 오차 감소를 위한 매개수 조정이 가능하다.
• 극한 $ \gamma \to 0^+ $에서 표준 FEM의 오차 추정치를 회복하며, 이는 $ k^3h^2 \leq C_0 $ 조건 하에서 유효하다.

실험 결과

연구 질문

• RQ1 복소 벌점 매개수를 가진 연속 내부 벌점 FEM이 큰 파장 수를 가진 헬름홀츠 방정식에 대해 무조건적 안정성을 달성할 수 있는가?
• RQ2 CIP-FEM의 $ H^1 $-노름에 대한 사전점근 오차 추정치는 무엇이며, $ k $, $ h $, $ \gamma $에 어떻게 의존하는가?
• RQ3 복소 벌점 매개수를 사용함으로써 고파장 수 헬름홀츠 문제에 대한 선형 FEM의 오염 오차를 크게 감소시킬 수 있는가?
• RQ4 벌점 매개수의 선택이 실질적인 수렴 행동과 상대 오차에 어떤 영향을 미치는가?
• RQ5 CIP-FEM는 더 비싼 방법들(예: IPDG)과 비교해 더 적은 자유도로 유사한 정확도를 달성할 수 있는가?

주요 결과

• 복소 벌점 매개수 $ \mathbf{i}\gamma $를 사용하는 CIP-FEM는 모든 $ k > 0 $에 대해 무조건적 안정성과 잘 정의됨을 보이며, 메쉬 제약 없이 적용 가능하다.
• $ H^1 $-노름 오차 추정치는 $ k^3h^2 \leq C_0 $일 때 $ C_1kh + C_2k^3h^2 $이며, $ k^3h^2 > C_0 $일 때는 $ C_1kh + \frac{C_2}{\gamma} $이다. 상수들은 $ k $, $ h $, $ \gamma $에 독립적이다.
• 최적의 $ L^2 $-노름 오차 추정치가 유도되어 방법의 강건성을 확인한다.
• 수치 결과는 벌점 매개수 $ \gamma $의 조정이 오염 오차를 크게 감소시킴을 보이며, 낮은 $ kh $ 값에서도 표준 FEM보다 CIP-FEM가 뛰어난 성능을 보인다.
• $ k=100 $일 때, CIP-FEM는 오직 35,971개의 자유도로 $ H^1 $-반노름에서 30%의 상대 오차를 달성하지만, 표준 FEM는 229,357개의 자유도가 필요하다.
• CIP-FEM의 경우 $ h(k,0.5) \geq 0.001 $를 만족하는 최대 파장 수는 $ k_{\text{max}} = 622 $이며, 표준 FEM보다 두 배 이상 높아 오염 효과의 심각성 감소를 입증한다.