[논문 리뷰] Error estimates and a two grid scheme for approximating transmission eigenvalues

Yidu Yang, Jiayu Han외 1명|arXiv (Cornell University)|2015. 06. 22.

Advanced Numerical Methods in Computational Mathematics참고 문헌 18인용 수 11

한 줄 요약

이 논문은 헬름홀츠 산산 문제에서의 전파 고유값을 계산하기 위해 이중 그리드 유한요소 방법을 제안한다. 비자기적이고 이차형 고유값 문제를 선형화 기법을 통해 자기를 만족하고 강력한 형태로 변환함으로써 엄밀한 오차 분석이 가능하게 한다. 이 방법은 실수 및 복소 고유값을 모두 고급 정확도로 계산하며, 최적의 수렴 속도를 달성한다. 이는 이론적 오차 추정과 다양한 영역에서의 수치 실험을 통해 검증되었다.

ABSTRACT

In this paper, using the linearization technique we write the Helmholtz transmission eigenvalue problem as an equivalent nonselfadjoint linear eigenvalue problem whose left-hand side term is a selfadjoint, continuous and coercive sesquilinear form. To solve the resulting nonselfadjoint eigenvalue problem, we give an $H^{2}$ conforming finite element discretization and establish a two grid discretization scheme. We present a complete error analysis for both discretization schemes, and theoretical analysis and numerical experiments show that the methods presented in this paper can efficiently compute real and complex transmission eigenvalues.

연구 동기 및 목표

• 헬름홀츠 산산 문제에서 전파 고유값을 계산하기 위한 강력하고 정확한 유한요소 방법을 개발하기 위해.
• 전파 고유값 문제의 비자기성과 타원성 부족 문제를 극복하기 위해.
• 일반 조건 하에서 실수 및 복소 고유값에 대해 유효한 완전한 오차 분석을 제공하기 위해.
• 계산 효율성을 높이면서도 정확도를 유지하는 이중 그리드 이산화 체계를 설계하기 위해.
• 볼록 및 비볼록 영역에서의 수치 실험을 통해 방법의 효과성을 입증하기 위해.

제안 방법

• 논문 [33]의 기법을 사용하여 비자기적 이차 고유값 문제 (5)를 등가의 비자기적 선형 고유값 문제 (11)로 선형화함으로써 좌변이 자기를 만족하고 연속적이며 강력한 형태로 보장한다.
• 선형화된 문제에 $H^2$-일치 유한요소 이산화를 적용하여 양의 정부호 헤르미트 블록 대각 행렬 형태의 강성 행렬을 도출한다.
• 바부쉬카-오즈본 스펙트럼 근사 이론을 적용하여 유한요소 근사에 대한 엄밀한 오차 추정을 유도한다.
• 비대칭 우변에 포함된 도함수를 다루기 위해 보조 문제 두 개를 도입하고, 낮은 노름 오차 분석에서 니츠체 기법을 적용한다.
• 먼저 거친 그리드 문제를 풀고, 그 해를 이용해 미세 그리드 근사를 보정하는 이중 그리드 이산화 체계를 제안하여 계산 비용을 감소시킨다.
• 삼각형 메esh에서 아르기리스 및 BFS 유한요소를 사용하여 방법을 구현하였으며, 비볼록 영역에서의 특이성을 고려해 메쉬의 적응적 세분화를 적용하였다.

실험 결과

연구 질문

• RQ1 전파 고유값의 $H^2$-일치 유한요소 근사에 대해 실수 및 복소 고유값, 다重 고유값을 포함한 완전한 오차 분석을 수립할 수 있는가?
• RQ2 전파 고유값 문제의 비자기성과 비타원성 문제는 어떻게 극복하여 신뢰할 수 있는 수치 계산을 가능하게 할 수 있는가?
• RQ3 제안된 이중 그리드 체계의 실수 및 복소 전파 고유값에 대한 수렴 거동은 어떠한가?
• RQ4 비볼록 영역에서 고유함수의 특이성은 유한요소 근사의 수렴 속도에 어떤 영향을 미치는가?
• RQ5 제안된 방법은 실수 및 복소 전파 고유값을 고정밀도로 효율적으로 계산할 수 있는가?

주요 결과

• 제안된 $H^2$-일치 유한요소 방법은 볼록 영역에서는 고유값에 대해 약 2의 최적 수렴 순서를 달성하며, 비볼록 L자형 영역에서는 특이성으로 인해 고유값의 수렴 속도가 약 1.4 수준이다.
• 이중 그리드 체계는 계산 비용을 크게 감소시키면서도 높은 정확도를 유지한다. 거친 그리드 고유값 $k_{j,H}$ 는 $h \to 0$ 일 때 미세 그리드 해 $k_{j,h}$ 로 수렴한다.
• 단위 정사각형과 L자형 도메인에서의 수치 실험은 이론적 수렴 속도를 확인하였으며, 메쉬 크기 $h$ 가 세분화될수록 상대 오차가 감소함을 보였다.
• 조각별 상수 $n(x)$ 를 가진 원형 도메인의 경우 아르기리스 요소는 [19]에서 제공한 해석 기준과 일치하는 고유값을 도출하여 방법의 정확도를 검증하였다.
• 복소 전파 고유값, 예를 들어 $n=8+x_1-x_2$ 에서 $\pm 0.566i$ 와 같은 허수부를 가진 고유값도 여러 메쉬 세분화 단계에서 고정밀도로 성공적으로 계산하였다.
• 이산화 시스템은 양의 정부호 헤르미트 블록 대각 강성 행렬을 생성하여 양호한 대수적 성질과 안정적인 반복 해법기를 보장한다.