[논문 리뷰] Numerical analysis of nonlinear eigenvalue problems

Éric Cancès, Rachida Chakir외 1명|ArXiv.org|2009. 05. 11.

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

한 줄 요약

이 논문은 $\|u\|_{L^{2}}=1$ 조건을 만족하는 비선형 타원형 고유값 문제 $-\mbox{div}(A\nabla u)+Vu+f(u^{2})u=\lambda u$의 기본 상태 고유쌍에 대한 변분 근사의 사전 오차 추정을 제시한다. 에너지, 고유값, 고유벡터에서의 오차 수렴 속도를 확립하며, 일반 조건 하에서 $|\lambda_\delta - \lambda|$가 $\|u_\delta - u\|_{H^{1}}^2 + \|u_\delta - u\|_{L^{2}}$ 속도로 수렴하고, 더 강한 가정 하에서는 $\|u_\delta - u\|_{H^{1}}^2$ 속도로 수렴함을 보이며, 이는 음의 소볼레프 노름 추정을 통해 이루어진다.

ABSTRACT

We provide a priori error estimates for variational approximations of the ground state eigenvalue and eigenvector of nonlinear elliptic eigenvalue problems of the form $-{div} (A abla u) + Vu + f(u^2) u = λu$, $\|u\|_{L^2}=1$. We focus in particular on the Fourier spectral approximation (for periodic problems) and on the $¶_1$ and $¶_2$ finite-element discretizations. Denoting by $(u_δ,λ_δ)$ a variational approximation of the ground state eigenpair $(u,λ)$, we are interested in the convergence rates of $\|u_δ-u\|_{H^1}$, $\|u_δ-u\|_{L^2}$ and $|λ_δ-λ|$, when the discretization parameter $δ$ goes to zero. We prove that if $A$, $V$ and $f$ satisfy certain conditions, $|λ_δ-λ|$ goes to zero as $\|u_δ-u\|_{H^1}^2+\|u_δ-u\|_{L^2}$. We also show that under more restrictive assumptions on $A$, $V$ and $f$, $|λ_δ-λ|$ converges to zero as $\|u_δ-u\|_{H^1}^2$, thus recovering a standard result for {\em linear} elliptic eigenvalue problems. For the latter analysis, we make use of estimates of the error $u_δ-u$ in negative Sobolev norms.

연구 동기 및 목표

• 비선형 타원형 고유값 문제에서 기본 상태 에너지, 고유값, 고유벡터에 대한 변분 근사의 사전 오차 추정을 수립하는 것.
• 문제 데이터의 다양한 정규성 조건 하에서 고유벡터와 고유값의 오차 수렴 속도를 $H^1$ 및 $L^2$ 노름에서 분석하는 것.
• 음의 소볼레프 노름 추정을 활용하여 선형 고유값 문제의 표준 수렴 결과를 비선형 설정으로 확장하는 것.
• 푸리에 스펙트럼, $\mathbb{P}_1$, $\mathbb{P}_2$ 유한요소 방법과 같은 특정 이산화 방법에 대한 수렴 거동을 검증하는 것.

제안 방법

• 비선형 고유값 문제의 변분 공식화를 기반으로 하며, 기본 상태는 $L^2$-정규화 제약 조건 하에서 에너지 함수의 최소화자로 특징지어진다.
• 문제는 고유값 문제 $A_u u = \lambda u$ 및 $\|u\|_{L^2} = 1$를 만족하는 자기수반 연산자 $A_u = -\mbox{div}(A\nabla \cdot) + V + f(u^2)$로 재정의된다.
• 이산 해와 연속 해의 차이 $u_\delta - u$를 이용하여 오차 추정을 도출하며, 이는 $H^1$ 및 $L^2$ 노름에서 분석된다.
• 음의 소볼레프 노름 추정을 활용하여 고유값 오차 분석을 정밀화하며, 특히 더 강한 가정 하에서 $\|u_\delta - u\|_{H^1}^2$ 수렴 속도를 회복한다.
• 수렴 분석은 푸리에 스펙트럼 및 $\mathbb{P}_1$, $\mathbb{P}_2$ 유한요소 이산화에 적용되며, 그 결과는 $A$, $V$, $f$의 정규성에 따라 달라진다.
• 기능 $F(t) = \int_0^t f(s)\,ds$는 $[0,\infty)$에서 $C^1$이며 $(0,\infty)$에서 $C^2$이며 $F'' > 0$를 만족하며, $f = F'$는 정밀 추정을 위해 (8)형 허더 조건을 만족한다.

실험 결과

연구 질문

• RQ1 비선형 고유값 문제에서 이산화 매개변수 $\delta \to 0$일 때 고유벡터와 고유값의 오차는 어떻게 수렴하는가?
• RQ2 고유값 오차 $|\lambda_\delta - \lambda|$가 $\|u_\delta - u\|_{H^1}^2 + \|u_\delta - u\|_{L^2}$ 속도로 수렴하는 조건은 무엇인가?
• RQ3 선형 고유값 문제에서의 표준 $\|u_\delta - u\|_{H^1}^2$ 수렴 속도가 비선형 케이스에서도 회복될 수 있는가?
• RQ4 음의 소볼레프 노름 추정은 비선형 문제에서 고유값 오차 분석을 어떻게 정밀화하는가?
• RQ5 고유항의 계수 $A$, $V$, $f$의 정규성에 따라 수렴 속도는 어떻게 달라지는가?

주요 결과

• 일반적인 $A$, $V$, $f$에 대한 가정 하에서 고유값 오차는 $|\lambda_\delta - \lambda| \lesssim \|u_\delta - u\|_{H^1}^2 + \|u_\delta - u\|_{L^2}$를 만족한다.
• 더 강한 가정, 특히 $f$에 대한 (8)형 허더 조건이 만족될 경우, 고유값 오차는 $|\lambda_\delta - \lambda| \lesssim \|u_\delta - u\|_{H^1}^2$ 속도로 수렴하며, 선형 문제의 수렴 속도를 회복한다.
• 타원형 정규성과 하르나크 부등식에 의해 기본 상태 해 $u$는 유일하고, 양수이며, 어떤 $\alpha \in (0,1)$에 대해 $C^{0,\alpha}(\overline{\Omega})$에 속한다.
• 고유값 $\lambda$는 연산자 $A_u$의 단순 고유값이며, 비선형 고유값 문제의 다른 해는 반드시 더 큰 고유값을 가져야 한다.
• 에너지 함수는 $L^2$-노름이 1인 음이 아닌 함수의 집합에서 엄격히 볼록하므로, 최소화자 $\rho_0 = u^2$의 유일성을 보장한다.
• 수렴 분석은 푸리에 스펙트럼 및 $\mathbb{P}_1$, $\mathbb{P}_2$ 유한요소 방법에 적용되며, 그 수렴 속도는 문제 데이터의 정규성에 따라 달라진다.