[논문 리뷰] Computation of Transmission Eigenvalues for Elastic Waves

Xia Ji, Peijun Li외 1명|arXiv (Cornell University)|2018. 02. 11.

Numerical methods in inverse problems참고 문헌 18인용 수 5

한 줄 요약

이 논문은 탄성파의 가장 작은 실수 전파 고유값을 계산하기 위한 새로운 수치적 방법을 제시한다. 이는 비선형적이며 자기수반적이지 않은 네째차 미분방정식 문제로서 매우 도전적인 문제이다. 고유값 문제를 전파 고유값에 해당하는 근을 가진 비선형 함수로 재구성하고, 자기수반 부분문제에 대해 H²-일치 유한요소를 사용함으로써, 이 방법은 루트를 구하기 위해 단순 반복법을 적용한다. 수렴성은 증명되었고, 혼합 유한요소 방법을 통해 검증되었다. 본 논문은 탄성 전파 고유값에 대한 첫 번째 종합적인 수치 연구이다.

ABSTRACT

The goal of this paper is to develop numerical methods computing a few smallest elasticity transmission eigenvalues, which are of practical importance in inverse scattering theory. The problem is challenging since it is nonlinear, non-self-adjoint, and of fourth order. We construct a nonlinear function whose values are generalized eigenvalues of a series of self-adjoint fourth order problems. The roots of the function are the transmission eigenvalues. Using an $H^2$-conforming finite element for the self-adjoint fourth order eigenvalue problems, we employ a secant method to compute the roots of the nonlinear function. The convergence of the proposed method is proved. In addition, a mixed finite element method is developed for the purpose of verification. Numerical examples are presented to verify the theory and demonstrate the effectiveness of the two methods.

연구 동기 및 목표

• 탄성파의 가장 작은 실수 전파 고유값을 효과적으로 계산하기 위한 수치적 방법을 개발하는 것. 이는 재료 물성 추정을 위한 역산산란에서 핵심적인 역할을 한다.
• 탄성파 전파에서 발생하는 비선형적이며 비자기수반인 네째차 고유값 문제를 해결하는 데 도전하는 것.
• 근이 정확히 전파 고유값에 해당하는 비선형 함수 F(ω)를 구성하는 것. 이를 통해 표준 반복 해법을 이용한 근 찾기 가능하게 한다.
• 결과의 수치적 신뢰성을 확보하기 위해 혼합 유한요소 방법을 사용하여 주요 방법의 결과를 검증하는 것.
• 탄성 전파 고유값 문제에 대해 수렴성과 강건성을 보장하는 컴퓨팅 프레임워크를 제공하는 것. 이는 이전에 수치적으로 탐구되지 않은 주제이다.

제안 방법

• 탄성 전파 고유값 문제를 네째차 미분방정식 형태의 이차 고유값 문제로 재구성하는 것.
• 근이 전파 고유값에 정확히 해당하는 비선형 함수 F(ω)를 구성하는 것. F(ω)는 자기수반 네째차 문제의 일반화된 고유값으로 정의된다.
• 자기수반 네째차 고유값 문제의 이산화를 위해 H²-일치 Argyris 유한요소를 사용함으로써 안정성과 정확성을 확보하는 것.
• F(ω)의 근을 계산하기 위해 단순 반복 기반 반복 방법을 적용하는 것. 이는 F(ω)가 강한 자기수반 문제의 일반화된 고유값임을 이용한다.
• Lagrange 요소를 사용한 혼합 유한요소 방법을 개발하여 주요 방법의 결과를 교차 검증하는 검증 도구로 활용하는 것.
• 경계 조건을 처리할 때 접선 및 법선 도함수를 신중하게 적용하여 경계에서 응력 자유 조건과 변위가 0인 조건을 일관되게 유지하는 것.

실험 결과

연구 질문

• RQ1 탄성 전파 고유값 문제의 비선형적이고 비자기수반적인 성질을 수치 계산에서 어떻게 극복할 수 있는가?
• RQ2 탄성파 산란에서 가장 작은 실수 전파 고유값을 효과적이고 수렴성 보장된 유한요소 방법으로 어떻게 계산할 수 있는가?
• RQ3 근이 정확히 전파 고유값에 해당하는 비선형 함수를 구성할 수 있는가? 이를 통해 표준 반복 해법을 이용한 근 찾기가 가능해지는가?
• RQ4 다른 유한요소 설정을 통해 계산된 고유값의 정확성과 신뢰성을 어떻게 검증할 수 있는가?
• RQ5 재구성된 일반화된 고유값 문제에 대해 제안된 단순 반복 기반 방법의 수렴 특성은 어떠한가?

주요 결과

• 제안된 방법은 자기수반 네째차 고유값 문제로부터 유도된 비선형 함수에 대해 단순 반복 기반 근 찾기 방법을 적용하여 탄성파의 가장 작은 실수 전파 고유값을 성공적으로 계산한다.
• 비선형 함수 F(ω)가 근 주변에서 연속적이며 엄격히 증가함을 가정할 경우, 단순 반복법의 수렴성이 엄밀히 증명된다.
• H²-일치 Argyris 유한요소 방법은 자기수반 네째차 고유값 문제에 대해 안정적이고 정확한 이산화를 제공하여 일반화된 고유값의 신뢰성 있는 계산을 가능하게 한다.
• 수치 실험은 이론적 수렴 속도를 확인하고, 다양한 기하형상과 재료 매개변수에서 방법의 강건성과 정확성을 입증한다.
• Lagrange 요소를 사용한 혼합 유한요소 방법은 신뢰할 수 있는 검증 경로를 제공하며, 주요 방법의 결과와 매우 유사한 고유값을 도출한다.
• 복소평면 상의 |Z₀(ω)|의 등고선도는 원형 곡선의 중심에 전파 고유값이 존재함을 확인하여, 반경 대칭 케이스에서의 해석적 해와 수치 결과가 일치함을 검증한다.