[논문 리뷰] Discrete maximal regularity and the finite element method for parabolic equations

Tomoya Kemmochi, Norikazu Saito|arXiv (Cornell University)|2016. 02. 22.

Advanced Numerical Methods in Computational Mathematics참고 문헌 23인용 수 3

한 줄 요약

이 논문은 선형 및 양형 포아르 방정식에 적용된 유한요소법에 대해 이산 최대정규성의 버전을 수립한다. 순허수거듭제곱을 포함한 연산자 이론적 기법을 사용하여 질량 래핑 유한요소 스킴에 대해 $L^p(0,T;L^q(\theta))$ 노름에서 최적의 오차 추정치를 도출하며, 이는 선형 열 방정식과 국소 리프시츠 비선형성을 가진 양형 문제에 응용된다.

ABSTRACT

Maximal regularity is a fundamental concept in the theory of partial differential equations. In this paper, we establish a fully discrete version of maximal regularity for a parabolic equation. We derive various stability results in $L^p(0,T;L^q(Ω))$ norm, $p,q\in (1,\infty)$ for the finite element approximation with the mass-lumping to the linear heat equation. Our method of analysis is an operator theoretical one using pure imaginary powers of operators and might be a discrete version of G.~Dore and A.~Venni (On the closedness of the sum of two closed operators. \emph{Math.\ Z.}, 196(2):189--201, 1987). As an application, optimal order error estimates in that norm are proved. Furthermore, we study the finite element approximation for semilinear heat equations with locally Lipschitz continuous nonlinearity and offer a new method for deriving optimal order error estimates. Some interesting auxiliary results including discrete Gagliardo-Nirenberg and Sobolev inequalities are also presented.

연구 동기 및 목표

• 포아르 방정식에 대해 유한요소법을 사용하여 최대정규성 이론의 이산적 형태를 개발한다.
• 선형 및 양형 열 방정식에 대해 $L^p(0,T;L^q(\theta))$ 노름에서 안정성과 최적의 오차 추정치를 도출한다.
• 강한 정규성 조건을 요구하지 않는 국소 리프시츠 비선형성을 가진 양형 문제를 위한 새로운 오차 분석 프레임워크를 제공한다.
• 분석을 위한 보조 도구로 이산 가가르도-니레버그 및 소볼레프 부등식을 수립한다.

제안 방법

• 순허수거듭제곱을 포함한 연산자 이론적 접근법을 사용하며, Dore와 Venni(1987)의 연구에 영감을 받는다.
• 선형 열 방정식의 공간 반디스크리티제이션에 질량 래핑 유한요소 방법을 사용한다.
• 시간에 대한 가중 평균(θ-방법)을 포함한 크랭크-니콜슨 또는 후진 오일러 시간 이산화를 통해 완전 이산 스킴을 정의한다.
• 스펙트럼 이론과 $H^\infty$-함수 해석학을 적용하여 이산 연산자 성질을 분석한다.
• $p,q\in(1,\infty)$에 대해 $L^p(0,T;L^q(\theta))$에서 이산 최대정규성 추정치를 도출한다.
• 양형 문제에서 비선형성의 전역 리프시츠 연속성 조건을 요구하지 않는 새로운 방법을 통해 최적의 오차 경계를 확립한다.

실험 결과

연구 질문

• RQ1 유한요소법을 포아르 방정식에 적용할 때 완전 이산 최대정규성의 버전을 수립할 수 있는가?
• RQ2 질량 래핑을 사용할 경우 선형 및 양형 포아르 방정식 문제에 대해 $L^p(0,T;L^q(\theta))$ 노름에서 최적의 오차 추정치는 무엇인가?
• RQ3 시간 이산화 선택(예: $\theta=0, 1/2, 1$)이 이산 최대정규성 프레임워크 내 수렴 차수에 어떤 영향을 미치는가?
• RQ4 이 방법을 사용하여 오직 국소 리프시츠 연속 비선형성을 가진 양형 포아르 방정식에 대해 최적의 오차 추정치를 도출할 수 있는가?
• RQ5 이산 가가르도-니레버그 및 소볼레프 부등식은 유한요소 스킴의 오차 분석에서 어떤 역할을 하는가?

주요 결과

• $\theta = 1/2$ 또는 $\theta = 1$일 때, 적절한 시간 간격 조건 하에서 질량 래핑 유한요소법에 대해 $L^p(0,T;L^q(\Omega))$에서 $O(h^2 + \tau)$의 최적 오차 추정치가 확립된다.
• $\theta = 0$일 경우, $\tau = O(h^2)$ 조건 하에서 $O(h^2)$ 수렴을 달성하며, 이는 이론적 경계와 일치한다.
• 스킴 (25)는 이론적 예측과 일치하는 최적의 수렴 속도를 보이며, $\theta = 1/2$ 및 $\theta = 1$일 때 관측된 수렴 차수는 $O(h^2)$이다.
• 반면, 보조 스킴 (97)은 이론적으로 $O(h)$ 오차 경계만을 제공하지만, 수치 결과에서는 실제로 $O(h^2)$ 수렴을 보이며, 이는 경계가 최적적이지 않음을 시사한다.
• $\theta = 1$일 때 스킴 (97)을 사용한 수치 실험에서는 $h$가 작아질수록 수렴 속도가 둔화되며, $\alpha \in [1,2)$에 대해 $O(h^\alpha)$로 수렴하는 것으로 나타나, 최적의 수렴 차수를 상실할 가능성이 있음을 시사한다.
• 논문은 $L^p(L^q)$ 프레임워크 내 오차 분석에 필수적인 이산 가가르도-니레버그 및 소볼레프 부등식의 이산 형태를 수립한다.