[논문 리뷰] Numerical validation of blow-up solutions with quasi-homogeneous compactifications

Kaname Matsue, Akitoshi Takayasu|arXiv (Cornell University)|2017. 07. 19.

Mathematical Biology Tumor Growth참고 문헌 9인용 수 7

한 줄 요약

이 논문은 무한대에서 점 渐近적 준동차성(asympotic quasi-homogeneity)을 가지는 유한차원 상미분방정식에서 폭발 해를 위한 엄밀한 수치 검증 방법을 제시한다. 새로운 준포물형(compactification)을 사용하여 발산하는 궤도를 수축된 무한대의 평형점에 도달하는 전역 궤도로 변환한다. 엄밀한 수치 적분, 리아푸노프 함수 검증, 시간에 따라 변하는 정규화(desingularization)를 결합함으로써, 오차 제어가 보장된 폭발 시간의 상한 및 하한, 해의 형태 봉투를 계산한다.

ABSTRACT

We provide a numerical validation method of blow-up solutions for finite dimensional vector fields admitting asymptotic quasi-homogeneity at infinity. Our methodology is based on quasi-homogeneous compactifications containing a new compactification, which shall be called a quasi-parabolic compactification. Divergent solutions including blow-up solutions then correspond to global trajectories of associated vector fields with appropriate time-variable transformation tending to equilibria on invariant manifolds representing infinity. We combine standard methodology of rigorous numerical integration of differential equations with Lyapunov function validations around equilibria corresponding to divergent directions, which yields rigorous upper and lower bounds of blow-up times as well as rigorous profile enclosures of blow-up solutions.

연구 동기 및 목표

• 무한대에서 점 渐近적 준동차성(asympotic quasi-homogeneity)을 가지는 ODE 시스템에서 폭발 해를 엄밀하게 검증하기 위한 수치 방법을 개발하는 것.
• 이전의 컴팩티피케이션 기반 검증 방법을 점 渐近적 동차 시스템을 초월하여 더 넓은 범주인 벡터장에까지 확장하는 것.
• 준동차 벡터장의 스케일링 구조를 유지하는 새로운 종류의 컴팩티피케이션—특히 준포물형 컴팩티피케이션—을 구성하는 것.
• 구간 산술과 리아푸노프 함수 기법을 사용하여 폭발 시간과 해의 형태를 검증 가능한 방식으로 계산하는 것.
• 폭발이 물리적·수학적으로 중요한 시스템, 예를 들어 Keller-Segel 모델과 같은 시스템에 적용 가능한 프레임워크를 제공하는 것.

제안 방법

• 방법은 단위 구조로 된 다양체에 위상공간을 매립시키기 위해 준포물형 컴팩티피케이션을 사용하여, 무한대를 경계 다원체로 매핑한다.
• 시간에 따라 변하는 변환(정규화)을 통해 시스템을 재매개변수화하여, 폭발 해가 수축된 경계의 평형점으로 수렴하는 전역 궤도가 되도록 한다.
• 유한 시간 간격 동안 해 궤도를 경계하는 데 구간 산술과 선형 산술을 사용한 엄밀한 수치 적분을 수행한다.
• 경계 평형점 주변에서 리아푸노프 함수를 검증하여 수렴 속도를 추정하고 무한대에 도달하는 데 소요되는 시간을 경계한다.
• 리아푸노프 함수의 감쇠 속도에서 유도된 오차 추정치를 통합하여 폭발 시간의 검증 가능한 상한 및 하한을 계산한다.
• 변환과 벡터장 성분의 다항식 근사치를 사용하여, 수축된 시스템이 원래 역학과 얼마나 떨어져 있는지 경계한다.

실험 결과

연구 질문

• RQ1 준동차 ODE 시스템에서 폭발 해를 컴팩티피케이션 기반 방법으로 엄밀하게 검증할 수 있는가?
• RQ2 준동차 벡터장의 점 점근적 스케일링 구조를 유지할 수 있는 적절한 컴팩티피케이션을 어떻게 구성할 수 있는가?
• RQ3 리아푸노프 함수는 검증 가능한 오차 경계를 갖는 폭발 시간 추정에 어떤 역할을 하는가?
• RQ4 이 방법은 동차 시스템을 초월하여 더 일반적인 준동차 역학으로 확장될 수 있는가?
• RQ5 검증 가능한 수치 적분과 구간 산술은 정규화와 어떻게 조합되어 폭발 프로파일과 시간을 경계할 수 있는가?

주요 결과

• 논문은 준동차 시스템에 대해 일반화된 Poincaré형 컴팩티피케이션을 초월하는 새로운 종류의 컴팩티피케이션—준포물형 컴팩티피케이션—을 도입한다.
• 폭발 해는 수축된 벡터장의 전역 궤도로서 수축된 경계의 평형점으로 수렴하는 것으로 특징지어지며, 이는 엄밀한 검증이 가능하게 한다.
• 리아푸노프 함수 감쇠 추정치의 통합을 통해 오차 경계가 보장된 폭발 시간의 검증 가능한 상한 및 하한이 계산된다.
• 경계 평형점 주변에서 해 궤도를 경계함으로써 폭발 해의 검증 가능한 프로파일 봉투를 확보한다.
• $k=1$인 경우, 폭발 시간 상한은 명시적으로 $t_{\text{max}} \leq t_N + \frac{1}{c_{\tilde{N}}} \left\{ 2\|v_1\|c_1^{-1/2}L^{1/2} + \sum_{j=2}^{\max\{2\beta_i\}} \frac{2}{j}\|v_j\|_{\infty}c_1^{j/2-1}L^{j/2} \right\}$로 계산되며, 여기서 $L$은 리아푸노프 함수 값이다.
• 이 방법은 $h(u,v) = u^2 - v$와 같이 표준 동차 컴팩티피케이션으로 다루기 어려운 비동차 스케일링을 가진 시스템에도 적용 가능하다.