[논문 리뷰] Maximum bound principles for a class of semilinear parabolic equations and exponential time differencing schemes

Qiang Du, Lili Ju외 2명|arXiv (Cornell University)|2020. 05. 23.

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

한 줄 요약

이 논문은 반선형 포아송 방정식에서 최대 bound 원칙(MBP)에 대한 이론적 프레임워크를 수립하며, 초기 조건과 경계 조건이 이러한 bound를 만족할 경우 해가 절대값에서 일관되게 유계로 유지됨을 증명한다. 이는 시간에 관계없이 MBP를 무조건적으로 유지하는 지수 시간 차분(ETD) 스킴을 개발하여 안정성과 고차 정확도를 동시에 확보한다. 또한 스칼라, 벡터 및 행렬값 시스템에 대해 엄밀한 오차 추정과 에너지 안정성 증명을 수행한다.

ABSTRACT

The ubiquity of semilinear parabolic equations has been illustrated in their numerous applications ranging from physics, biology, to materials and social sciences. In this paper, we consider a practically desirable property for a class of semilinear parabolic equations of the abstract form $u_t=\mathcal{L}u+f[u]$ with $\mathcal{L}$ being a linear dissipative operator and $f$ being a nonlinear operator in space, namely a time-invariant maximum bound principle, in the sense that the time-dependent solution $u$ preserves for all time a uniform pointwise bound in absolute value imposed by its initial and boundary conditions. We first study an analytical framework for some sufficient conditions on $\mathcal{L}$ and $f$ that lead to such a maximum bound principle for the time-continuous dynamic system of infinite or finite dimensions. Then, we utilize a suitable exponential time differencing approach with a properly chosen generator of contraction semigroup to develop first- and second-order accurate temporal discretization schemes, that satisfy the maximum bound principle unconditionally in the time-discrete setting. Error estimates of the proposed schemes are derived along with their energy stability. Extensions to vector- and matrix-valued systems are also discussed. We demonstrate that the abstract framework and analysis techniques developed here offer an effective and unified approach to study the maximum bound principle of the abstract evolution equation that cover a wide variety of well-known models and their numerical discretization schemes. Some numerical experiments are also carried out to verify the theoretical results.

연구 동기 및 목표

• 추상적인 반선형 포아송 방정식에 대해 시간에 관계없이 일정한 최대 bound 원칙(MBP)을 보장하기 위해 선형 소산 연산자 𝒮와 비선형 연산자 f에 대한 충분한 분석 조건을 수립하기 위해.
• 시간 스텝 크기와 무관하게 항상 MBP를 유지하는 1차 및 2차 지수 시간 차분(ETD) 스킴을 개발하기 위해.
• 제안된 ETD 스킴의 스칼라 및 벡터/행렬값 설정에서 엄밀한 오차 추정과 에너지 안정성 증명을 수행하기 위해.
• 디리클레, 뉴먼, 주기적 경계 조건을 포함한 시스템으로 추상적 프레임워크를 확장하고, 다양한 행렬 노름(2-노름 및 프로베니우스 노름)에 대한 MBP를 분석하기 위해.
• 적분 인자(IF) 스킴에서 시간에 관계없이 MBP를 유지할 수 있는지 탐색하고, ETD 방법과 비교하기 위해.

제안 방법

• 수축 반군 이론에 기반하여, 선형 연산자 𝒮와 비선형 연산자 f에 대한 분석 조건을 유도하여, 모든 시간에 걸쳐 해가 절대값에서 유계로 유지됨을 보장한다.
• 수축 반군의 생성자를 사용하여, 시간 이산화 설정에서 MBP를 무조건적으로 유지하는 ETD 스킴을 구성한다.
• 지수 적분기를 적용하여 시간적으로 1차 및 2차 정확도를 확보하면서도, 시간 스텝 제약 없이 MBP를 유지한다.
• 에너지 안정성 분석과 오차 추정을 통해 다양한 경계 조건 하에서 ETD 스킴의 수렴성과 강건성을 검증한다.
• 벡터 및 행렬값 시스템으로 프레임워크를 확장하고, ℓ², ℓ∞ 및 행렬 2-노름과 프로베니우스 노름에 대한 MBP를 분석한다.
• ETD와 적분 인자(IF) 스킴을 비교하여, IF 스킴은 시간 스텝 제약 하에만 MBP를 유지할 수 있지만, ETD 스킴은 무조건적으로 MBP 유지가 가능함을 보여준다.

실험 결과

연구 질문

• RQ1 선형 연산자 𝒮와 비선형 연산자 f에 대해 어떤 조건이 반선형 포아송 방정식이 시간 연속 설정에서 최대 bound 원칙(MBP)을 만족하는가?
• RQ2 시간 스텝 크기와 무관하게 항상 MBP를 유지하는 지수 시간 차분(ETD) 스킴을 구성할 수 있는가?
• RQ3 제안된 ETD 스킴의 스칼라, 벡터 및 행렬값 시스템에서의 오차 추정과 에너지 안정성 특성은 무엇인가?
• RQ4 행렬값 시스템에서 다양한 행렬 노름, 특히 2-노름과 프로베니우스 노름에 대해 MBP는 어떻게 행동하는가?
• RQ5 적분 인자(IF) 스킴은 무조건적으로 MBP를 유지할 수 있는가? 만약 그렇지 않다면 어떤 조건에서 조건부로 유지할 수 있는가?

주요 결과

• 논문은 초기 및 경계 데이터가 이러한 bound를 만족할 경우, 추상적인 반선형 포아송 방정식의 해가 모든 시간에 걸쳐 절대값에서 유계로 유지됨을 보장하는 𝒮와 f에 대한 충분한 조건을 수립한다.
• 1차 및 2차 ETD 스킴이 무조건적으로 MBP를 유지하도록 구성되었으며, 이는 시간 스텝 크기와 관계없이 항상 bound를 유지함을 의미하며, 안정성 제약 없이 적용 가능하다.
• 이론적 오차 추정을 유도하여, 해와 데이터에 적절한 정규성 조건이 만족될 경우 ETD 스킴의 최적 수렴 속도를 보여준다.
• ETD 스킴에 대해 에너지 안정성이 증명되었으며, 이는 이산 에너지가 시간이 지남에 따라 증가하지 않음을 보장하여 장기적인 수치적 신뢰성을 지원한다.
• 행렬값 시스템에서는 ETD 스킴이 2-노름에 대한 MBP를 무조건적으로 유지하지만, 비선형 항을 유계로 둘 수 있는 도전 과제로 인해 프로베니우스 노름에 대한 경우는 여전히 열린 문제로 남아 있다.
• IF 스킴은 시간 스텝 제약 하에만 MBP를 유지할 수 있지만, 적절한 안정화 요소를 도입하면 무조건적 MBP 유지가 가능할 수 있음을 제안하며, 향후 연구 방향을 제시한다.