[논문 리뷰] Robust numerical methods for local and nonlocal equations of porous medium type. Part I: Theory

Félix del Teso, Jørgen Endal외 1명|arXiv (Cornell University)|2018. 01. 22.

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

한 줄 요약

이 논문은 국소적 및 비국소적 레비 유형의 연산자를 포함하는 일반적인 비선형 확산 방정식에 적용된 단조 유한차분 스킴에 대한 통합적이고 강건한 이론을 개발한다. 최소한의 가정 하에 잘 정의됨, 안정성, 수렴성, 콤팩트성 등을 입증하며, 특히 매우 불규칙한 해에 대해서도 $L^1 \cap L^\infty$에서 분포해로의 수렴을 보인다.

ABSTRACT

We develop the theory for monotone schemes of finite difference type for general nonlinear diffusion equations, $$ \partial_t u-\mathfrak{L}^{\sigma,\mu}[\varphi(u)]=f(x,t) \qquad ext{in}\qquad \mathbb{R}^N imes(0,T), $$ where $\mathfrak{L}^{\sigma,\mu}$ is a general symmetric diffusion operator of L\'evy type and $\varphi$ is merely continuous and non-decreasing. This theory covers strongly degenerate Stefan problems and the full range of porous medium and fast diffusion equations. Examples of diffusion operators $\mathfrak{L}^{\sigma,\mu}$ are the Laplacian $\Delta$, the fractional Laplacian $-(-\Delta)^{\frac\alpha2}$, $\alpha\in(0,2)$, discrete operators, and combinations of these. The observation that finite difference type operators are nonlocal L\'evy operators, allows us to give a unified, compact, and efficient nonlocal theory for both local and nonlocal, linear and nonlinear diffusion equations. The theory includes well-posedness, stability, equicontinuity, compactness, and convergence of the methods under minimal assumptions including assumptions that lead to very irregular solutions. In other words, the schemes we introduce are robust in the sense that they converge under very unfavorable conditions. Such general results are possible because we are able to work in the largest class of weak solutions where uniqueness is known: distributional solutions in $L^1\cap L^\infty$. As a byproduct of our numerical analysis, we obtain a new and general existence result for the original PDE, a result that was previously announced. We also present some numerical tests, but extensive numerical testing is deferred to a companion paper along with a more detailed discussion of the numerical methods included in our theory.

연구 동기 및 목표

• 국소적 및 비국소적 확산 방정식에 모두 적용 가능한 단조 유한차분 스킴에 대한 일반적이고 통합적인 이론을 개발하는 것.
• 해의 정(regularity)이 낮은 경우를 포함하여 최소한의 가정 하에 수치 스킴의 수렴성과 안정성을 확립하는 것.
• 강하게 탈선하는 스텐프 문제와 다공성 매체 및 빠른 확산 방정식의 전 범위에 이론을 확장하는 것.
• 수치 해석의 부산물로 원래 PDE에 대한 새로운 일반적 존재 결과를 도출하는 것.
• 단순히 연속적이고 비감소적인 비선형성을 갖는 방정식에 적용 가능한 강건한 수치 방법의 이론적 기반을 마련하는 것.

제안 방법

• 대칭 레비 유형 연산자인 $ \mathfrak{L}^{\sigma,\mu} $를 갖는 일반 확산 방정식을 $ \partial_t u - \mathfrak{L}^{\sigma,\mu}[\varphi(u)] = f(x,t) $ 형태로 수식화한다.
• 기초 PDE와 일致하고 단조인 유한차분 형식의 스킴을 사용하여 안정성과 수렴성을 보장한다.
• 유일성에 대해 알려진 가장 큰 클래스인 $ L^1 \cap L^\infty $에서의 분포해의 공간 내에서 작업한다.
• 비국소 분석 기법을 적용하여 국소적(예: 라플라스 연산자) 및 비국소적(예: 분수 라플라스 연산자) 연산자를 통합적으로 다룬다.
• 수치 해의 등연속성과 콤팩트성을 확립하여 선형화 방법 또는 부분수열 추론을 통해 수렴성을 증명한다.
• 레비 연산자의 구조를 활용하여 이산, 연속, 하이브리드 확산 연산자를 통합적으로 다루는 데 기여한다.

실험 결과

연구 질문

• RQ1 최소한의 정규성 가정 하에 국소적 및 비국소적 비선형 확산 방정식에 대해 단일이고 강건한 수치 프레임워크를 개발할 수 있는가?
• RQ2 불연속적이거나 탈선적인 비선형성을 갖는 방정식에 대해 단조 유한차분 스킴이 약한 해로 수렴하는 조건은 무엇인가?
• RQ3 레비 유형 연산자의 이론은 국소적 및 비국소적 확산 과정의 분석을 어떻게 통합할 수 있는가?
• RQ4 수렴성과 안정성을 보장하기 위해 $ \varphi(u) $와 연산자 $ \mathfrak{L}^{\sigma,\mu} $에 대해 필요한 최소한의 가정은 무엇인가?
• RQ5 수치 해석 프레임워크가 원래 PDE 시스템에 대한 새로운 존재 결과를 도출할 수 있는가?

주요 결과

• 제안된 유한차분 스킴은 해가 매우 불규칙한 경우에도 $ L^1 \cap L^\infty $에서 분포해로 수렴한다.
• 이론은 탈선성과 비연속 비선형성 등 최소한의 가정 하에 잘 정의됨, 안정성, 등연속성, 콤팩트성 등을 확보한다.
• 레비 유형 연산자를 활용하여 국소적(예: 라플라스 연산자) 및 비국소적(예: 분수 라플라스 연산자) 확산 연산자의 분석을 통합한다.
• 수치 해석의 부산물로 원래 PDE에 대한 새로운 일반적 존재 결과가 도출되었으며, 이는 이전 결과를 확장한다.
• 이 방법은 강건하다: 고전적 정규성 또는 유일성 이론이 실패하는 경우에도 수렴성이 입증된다.
• 이론은 다공성 매체 및 빠른 확산 방정식의 전 범위에 적용되며, 강하게 탈선한 스텐프 문제를 포함한다.