[논문 리뷰] On the numerical approximation of $p$-Biharmonic and $\infty$-Biharmonic functions

Nikos Katzourakis, Tristan Pryer|arXiv (Cornell University)|2017. 01. 25.

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

한 줄 요약

이 논문은 p-비라멘드 및 ∞-비라멘드 함수의 수치적 근사를 위해 C⁰ 혼합 유한요소 방법을 개발하며, 이는 p-비라멘드 방정식의 약한 해로의 수렴을 증명하고, p→∞로의 엄밀한 극한을 취해 ∞-비라멘드 해를 포착한다. 이 방법은 ∞-비라멘드 함수가 단순한 경계 조건 하에서도 복잡한 인터페이스 패턴을 가진 조각별 이차 함수임을 드러낸다.

ABSTRACT

In [KP16] (arXiv:1605.07880) the authors introduced a second-order variational problem in $L^{\infty}$. The associated equation, coined the $\infty$-Bilaplacian, is a \emph{third order} fully nonlinear PDE given by $Δ^2_\infty u\, := (Δu)^3 | D (Δu) |^2 = 0.$ In this work we build a numerical method aimed at quantifying the nature of solutions to this problem which we call $\infty$-Biharmonic functions. For fixed $p$ we design a mixed finite element scheme for the pre-limiting equation, the $p$-Bilaplacian $Δ^2_p u\, := Δ(| Δu |^{p-2} Δu) = 0.$ We prove convergence of the numerical solution to the weak solution of $Δ^2_p u = 0$ and show that we are able to pass to the limit $p o\infty$. We perform various tests aimed at understanding the nature of solutions of $Δ^2_\infty u$ and in 1-$d$ we prove convergence of our discretisation to an appropriate weak solution concept of this problem, that of $\mathcal D$-solutions.

연구 동기 및 목표

• 최소한의 정규성 가정 하에 네 번째 차수의 완전 비선형 편미분방정식인 p-비라멘드 방정식을 위한 안정적이고 수렴하는 수치적 스킴을 개발하기 위해.
• 혼합 유한요소 방법이 p-비라멘드 문제의 약한 해로 수렴함을 확립하고, 추가 정규성 조건 하에서 최적의 수렴 속도를 확보하기 위해.
• p-비라멘드 방정식의 p→∞ 극한이 산술적으로도 분석적으로도 ∞-비라멘드 방정식으로 간주될 수 있음을 정당화하기 위해. 이는 발산형이 아닌 세 번째 차수의 완전 비선형 편미분방정식이다.
• 특히 조각별 정규성과 인터페이스 복잡성 측면에서 ∞-비라멘드 함수의 구조적 성질을 수치적으로 조사하기 위해.
• 고전적 해가 존재하지 않기 때문에 ∞-비라멘드 방정식에 대해 약한 해 개념(D-해)을 정의하고 수치적으로 검증하기 위해.

제안 방법

• 헤시안과 라플라스 연산자에 대한 보조 변수를 도입함으로써 C¹ 유한요소가 아닌 C⁰ 유한요소를 사용할 수 있도록 혼합 유한요소 설정을 수립한다.
• p-비라멘드 에너지 기능을 최소화하는 변분 설정을 사용하며, 이를 네 번째 차수의 편미분방정식을 분리할 수 있도록 혼합 형태로 재구성한다.
• 이전 연구들에서 영감을 얻은 inf-sup 조건이 안정성을 보장하고, 더 높은 정규성 조건 하에서 최적의 수렴 속도 분석이 가능하게 한다.
• 해리지 방식을 사용하여 다항식 근사로 이루어진 이산 시스템을 풀며, 컴actness 및 얀 측도 기법을 통해 수렴을 증명한다.
• p→∞ 극한은 L∞에서의 변분법 이론을 활용하여 분석되며, 이는 p-비라멘드 편미분방정식의 형식적 극한으로 ∞-비라멘드 방정식을 도출한다.
• 경계 조건에서 m의 다양한 값에 대해 2차원 영역에서 수치 실험을 수행하여, p가 증가함에 따라 Δu가 조각별 상수 함수로 수렴하는 경향을 관찰한다.

실험 결과

연구 질문

• RQ1 C¹ 유한요소가 필요 없이 p-비라멘드 방정식에 대해 안정적이고 수렴하는 유한요소 방법을 어떻게 설계할 수 있는가?
• RQ2 수치 해의 수렴 행동은 어떠한가? 정규성 가정 하에서 최적의 수렴 속도를 확보할 수 있는가?
• RQ3 p-비라멘드 방정식의 p→∞ 극한이 수치적·분석적으로 ∞-비라멘드 방정식으로 간주될 수 있는가?
• RQ4 ∞-비라멘드 함수는 조각별 정규성과 인터페이스 복잡성 측면에서 어떤 구조적 성질을 보이는가?
• RQ5 고전적 해가 존재하지 않기 때문에 ∞-비라멘드 방정식에 대해 약한 해 개념(D-해)을 어떻게 정의하고 수치적으로 검증할 수 있는가?

주요 결과

• 최소한의 정규성 가정 하에서 혼합 유한요소 방법이 p-비라멘드 방정식의 약한 해로 수렴하며, 컴팩트니스 및 얀 측도 기법을 통해 수렴이 증명된다.
• 추가 정규성 조건 하에서, 수렴 속도는 p의 값에 따라 달라지며, p가 클수록 안정성과 정확도가 향상된다.
• p가 증가함에 따라 p-비라멘드 방정식의 수치 해는 ∞-비라멘드 극한에 수렴하며, Δu는 조각별 상수 함수로 수렴함을 보여주어, 이는 해가 조각별 상수 우변을 가진 라플라스 방정식을 만족함을 시사한다.
• 경계 조건에서 m=1,2,3에 대한 수치 시뮬레이션은 ∞-비라멘드 함수가 영역 전체에 걸쳐 조각별 이차 함수임을 드러내며, 매우 복잡하고 비자명한 인터페이스 패턴을 보임을 확인한다.
• 이 방법은 ∞-비라멘드 방정식의 D-해 개념을 성공적으로 포착하여, 고전적 해가 존재하지 않음에도 약한 공식화의 타당성을 검증한다.
• 수치 실험에서 관측된 수렴 속도는 이론적 예측을 초과하며, 저차수 문제에서 알려진 바와 같이 준-노름 안정화를 통한 향상 가능성이 있음을 시사한다.