[논문 리뷰] Asymptotic Behaviours of Solutions for Finite Difference Analogue of the Chipot-Weissler Equation

Houda Hani, Moez Khenissi|arXiv (Cornell University)|2015. 02. 10.

Differential Equations and Numerical Methods참고 문헌 4인용 수 3

한 줄 요약

이 논문은 반응 및 기울기 의존성 확산 항을 포함한 비선형 포물형 편미분방정식인 Chipot-Weissler 방정식의 유한차분 수치근사법을 분석한다. 매개변수와 초기 자료에 대한 제약 조건 하에서 수치해가 연속해로 수렴하는 것을 입증하지만, 수치적 및 해석적 폭발 집합 간의 불일치를 드러내며, 수치적으로 점점 커지는 폭발 행동을 포착하는 데에 한계가 있음을 시사한다.

ABSTRACT

This paper deals with nonlinear parabolic equation for which a local solution in time exists and then blows up in a finite time. We consider the Chipot-Weissler equation. We study the numerical approximation, we show that the numerical solution converges to the continuous one under some restriction on the initial data and the parameters of the non linearity. Moreover, we study the numerical blow up sets and we show that although the convergence of the numerical solution is guaranteed, the numerical blow up sets are sometimes different from that of the PDE

연구 동기 및 목표

• 유한차분 스킴을 사용하여 Chipot-Weissler 방정식의 해의 수치적 행동을 연구하기 위해.
• 초기 자료와 매개변수 p 및 q에 대한 특정 제약 조건 하에서 수치해가 진짜 해로 수렴하는 지를 조사하기 위해.
• 수치적 폭발 집합을 연속 PDE의 실제 폭발 집합과 비교하기 위해.
• 연속 및 이산 설정에서 해의 점근적 행동을 분석하기 위해.

제안 방법

• Semi-implicit 유한차분 스킴을 Chipot-Weissler 방정식 $ u_t = u_{xx} + u^p - |u_x|^q $ 에 적용함: $ x \in (-1,1), t > 0 $.
• 균일한 격자로 공간과 시간을 이산화하고, 공간 도함수는 중심차분으로, 시간 도함수는 후진 오일러로 근사함.
• 초기 자료와 매개변수 $ p > 1 $, $ 1 \leq q \leq \frac{2p}{p+1} $ 에 대한 제약 조건 하에서 수치해가 연속해로 수렴함을 증명함.
• 최대값을 시간과 공간에 따라 추적하여 수치적으로 폭발 집합을 분석함.
• 에너지 추정과 비교 원리를 조합한 이론적 분석을 통해 해를 유계로 유지하고 폭발 시간을 제어함.
• 이산 및 연속 해의 폭발 집합을 비교하기 위해 수치 실험을 수행함.

실험 결과

연구 질문

• RQ1 초기 자료와 매개변수 p 및 q에 대해 적절한 조건 하에서 유한차분 스킴이 Chipot-Weissler 방정식의 진짜 해로 수렴하는가?
• RQ2 수치적 폭발 집합은 연속 PDE의 실제 폭발 집합과 어떻게 비교되는가?
• RQ3 특히 유한시간 폭발을 포함한 해의 점근적 행동은 수치 스킴에 의해 정확하게 포착될 수 있는가?
• RQ4 매개변수 p와 q는 수치해의 수렴성과 폭발 역학에 어떤 역할을 하는가?

주요 결과

• 적절한 초기 자료 제약 조건과 매개변수 $ p > 1 $, $ 1 \leq q \leq \frac{2p}{p+1} $ 하에서 수치해가 연속해로 수렴함.
• 수렴은 이루어지지만, 수치적 폭발 집합는 PDE의 실제 폭발 집합과 다를 수 있으며, 이는 이산 근사에서 폭발 구조를 유지하지 못함을 시사함.
• 수치해의 폭발 시간은 진짜 폭발 시간으로 수렴하지 않을 수 있으며, 특히 기울기 항 $ |u_x|^q $ 가 강할 경우 더욱 그렇다.
• 이론적 분석은 국소 해가 존재하고 유한시간 내에 폭발함을 확인하며, PDE의 행동과 일치함.
• 수치 시뮬레이션은 해의 수렴이 노름 기준으로 이루어져도 폭발 위치가 이산화 과정에서 이동할 수 있음을 보여줌.
• 폭발 집합 간의 불일치는 반응 항과 기울기 항 간의 비선형 결합으로 인해 발생하며, 이는 유한차분 스킴에서 완전히 유지되지 않기 때문이다.