[논문 리뷰] Numerical investigation of the smallest eigenvalues of the p-Laplace operator on planar domains

Jiří Horák|arXiv (Cornell University)|2011. 06. 03.

Advanced Mathematical Modeling in Engineering참고 문헌 15인용 수 15

한 줄 요약

이 논문은 $p \in [1.1, 10]$에 대해 평면 영역에서 $W^{1,p}_0(\Omega)$에서의 제약 내림과 산맥 길이 알고리즘을 사용하여 p-라플라스 연산자의 두 번째로 작은 고유값을 수치적으로 조사한다. 결과적으로, 원형 이외의 비원형 도메인인 직사각형과 삼각형에서 두 번째 고유함수의 대칭성은 $p$에 따라 변화하며, 이는 영역에 따라 결정되는 임계 $p_0$에서 전이된다. 반면 원형에서는 모든 $p$에 대해 노드선을 중심으로 한 홀함수 대칭성을 유지한다.

ABSTRACT

The eigenvalue problem for the p-Laplace operator with p>1 on planar domains with the zero Dirichlet boundary condition is considered. The Constrained Descent Method and the Constrained Mountain Pass Algorithm are used in the Sobolev space setting to numerically investigate the dependence of the two smallest eigenvalues on p. Computations are conducted for values of p between 1.1 and 10. Symmetry properties of the second eigenfunction are also examined numerically. While for the disk an odd symmetry about the nodal line dividing the disk in halves is maintained for all the considered values of p, for rectangles and triangles symmetry changes as p varies. Based on the numerical evidence the change of symmetry in this case occurs at a certain value p_0 which depends on the domain.

연구 동기 및 목표

• Hilbert 공간에서의 제약 산맥 길이 알고리즘을 $p \neq 2$인 Banach 공간 $W^{1,p}_0(\Omega)$로 확장하기.
• 넓은 범위($1.1 \leq p \leq 10$)에서 $p$에 따른 두 번째로 작은 고유값의 의존성 계산 및 분석하기.
• 특히 비원형 도메인에서 $p$ 변화에 따른 두 번째 고유함수의 대칭성 변화 조사하기.
• 2와 거리가 먼 $p$ 값에 대해 수치 방법의 타당성 검증 및 메esh 정밀도 증가에 따른 수렴성 및 안정성 평가하기.
• 알려진 점근적 행동($p \to 1^+$ 및 $p \to \infty$)과 수치 결과 비교하기.

제안 방법

• Banach 공간 $W^{1,p}_0(\Omega)$에 대해 $p \in (1, \infty)$에서 제약 내림 방법과 제약 산맥 길이 알고리즘을 적응시키기.
• 제약 조건 $\|u\|_{L^p} = 1$ 하에서 최소화할 기능으로 레일라이히 비율 $\lambda = \frac{\int_\Omega |\nabla u|^p \, dx}{\int_\Omega |u|^p \, dx}$를 사용하기.
• 내림 방향 $w_u$를 $-\Delta_p w = \Delta_p u + \alpha |u|^{p-2}u$의 해로 정의하고, $w_u \in T_u S$를 확보하기 위해 정규화하기.
• $(-\Delta_p)^{-1}$을 사용해 $|u|^{p-2}u$를 $W^{1,p}_0(\Omega)$에 속하는 함수로 매핑하여 반복 업데이트 가능하게 하기.
• 수렴성을 보장하기 위해 적응형 메쉬 정밀도를 사용한 유한요소 방법을 문제의 이산화에 적용하기.
• 초기 조건 $u(0) \in S$를 갖는 시간 스텝 스킴 $\frac{d}{dt}u(t) = w_{u(t)}$를 적용하여 임계점 수렴 추적하기.

실험 결과

연구 질문

• RQ1 1.1에서 10 사이에서 $p$가 변화함에 따라 평면 도메인에서 가장 작은 고유값 $\lambda_1$은 어떻게 행동하는가?
• RQ2 두 번째 고유함수 $u_2$는 다양한 $p$ 값에서 대칭성 특성(예: 노드선을 중심으로 한 홀함수 대칭)을 유지하는가?
• RQ3 비원형 도메인인 직사각형과 삼각형에서 $u_2$의 대칭성이 변화하는 임계 $p_0$는 무엇인가?
• RQ4 특히 수렴성 및 안정성 측면에서 2와 거리가 먼 $p$ 값에서 수치 방법의 성능은 어떠한가?
• RQ5 계산된 고유값과 고유함수는 $p \to 1^+$ 및 $p \to \infty$로 갈수록 알려진 점근적 근사와 어느 정도 일치하는가?

주요 결과

• 가장 작은 고유값 $\lambda_1$은 1.1에서 10으로 $p$가 증가함에 따라 단조 감소하며, 이는 이론적 기대와 일치한다.
• 원형에서는 모든 $p \in [1.1, 10]$에서 두 번째 고유함수 $u_2$가 원을 둘로 나누는 지름을 중심으로 홀함수 대칭성을 유지한다.
• 직사각형과 삼각형에서는 $u_2$의 대칭성이 도메인 기하학에 따라 결정되는 임계 $p_0$에서 변화함을 보이며, 고유구조에서 분기 현상이 발생함을 시사한다.
• 제약 산맥 길이 알고리즘이 2와 거리가 먼 $p$ 값에서도 두 번째 고유쌍에 대해 성공적으로 수렴함을 보이며, Banach 공간 설정에서의 강건성을 입증한다.
• $p \to 1^+$ 및 $p \to \infty$로 갈수록 수치 결과가 알려진 점근적 행동과 정성적으로 일치한다. 예를 들어 $p \to 1^+$일 때 $\lambda_1 \to \infty$이고, $p \to \infty$일 때 $\lambda_1 \to \text{유한한 극한}$이다.
• 역 $p$-라플라스 연산자에 기반한 내림 방향 정의가 임계점 수렴을 보장하며, 이 수렴성 증명은 $W^{1,p}_0(\Omega)$에서의 컴팩턴스와 노름 추정에 기반한다.