[논문 리뷰] Study of a chemo-repulsion model with quadratic production. Part I: Analysis of the continuous problem and time-discrete numerical schemes

Francisco Guillén‐González, María Ángeles Rodríguez‐Bellido외 1명|arXiv (Cornell University)|2018. 03. 06.

Mathematical Biology Tumor Growth참고 문헌 19인용 수 6

한 줄 요약

이 논문은 유계 영역에서 이차 생산을 가지는 화학적 반발력 모델을 연구하며, 1차원 및 2차원에서 약한 해의 전역 존재성과 유일성을 증명한다. 두 가지 조건 없이 에너지 안정인 일阶 시간 스킴—역행 오일러와 선형화된 결합 스킴—을 제안하여 해의 존재성, 수렴성 및 오차 추정을 입증한다. 선형 스킴는 양의 보존성이 없음에도 불구하고 수렴성을 보여준다.

ABSTRACT

We consider a chemo-repulsion model with quadratic production in a bounded domain. Firstly, we obtain global in time weak solutions, and give a regularity criterion (which is satisfied for $1D$ and $2D$ domains) to deduce uniqueness and global regularity. After, we study two cell-conservative and unconditionally energy-stable first-order time schemes: a (nonlinear and positive) Backward Euler scheme and a linearized coupled version, proving solvability, convergence towards weak solutions and error estimates. In particular, the linear scheme does not preserve positivity and the uniqueness of the nonlinear scheme is proved assuming small time step with respect to a strong norm of the discrete solution. This hypothesis is reduced to small time step in $nD$ domains ($n\le 2$) where global in time strong estimates are proved. Finally, we show the behavior of the schemes through some numerical simulations.

연구 동기 및 목표

• 유계 영역에서 이차 생산을 가지는 화학적 반발력 모델에 대한 전역 약한 해의 존재성을 확립하는 것.
• 규칙성 기준 하에 1차원 및 2차원 영역에서 해의 유일성과 전역 정칙성을 증명하는 것.
• 비선형 역행 오일러와 선형화된 결합 스킴을 포함한 두 가지 조건 없이 에너지 안정인 시간 이산화 스킴을 설계하고 분석하는 것.
• 두 스킴에 대해 해의 존재성, 약한 해로의 수렴성 및 오차 추정을 증명하는 것.
• 수치 시뮬레이션을 통해 스킴의 거동을 조사하며, 특히 에너지 안정성과 수렴성 특성에 중점을 두는 것.

제안 방법

• 이차 생산 $ f(u) = u^2 $를 가지는 화학적 반발력 시스템을 설정하고, 시스템을 분리하기 위해 보조 변수 $ \boldsymbol{\sigma} = \nabla v $ 를 도입한다.
• 유한한 영역에서의 노이만 경계 조건과 초기 자료를 포함한 $ u $, $ \boldsymbol{\sigma} $, $ v $ 를 포함하는 재구성된 시스템 (3)을 유도한다.
• 시간 이산화를 위해 비선형 역행 오일러 스킴 (104)을 제안하여 세포 보존성과 조건 없이 에너지 안정성을 확보한다.
• 계산 효율성이 뛰어난 대안으로 선형화된 결합 스킴을 도입하지만, 이는 양의 보존성을 보장하지 못한다.
• 레일레이-슈아더 고정점 정리를 적용하여, 컴actness와 연속성의 논증을 통해 시간 이산화 스킴의 해 존재성을 증명한다.
• 에너지 추정과 $ L^4 $-노름 유계를 사용하여 수렴성 및 오차 추정을 도출하며, 비선형 스킴의 유일성은 강한 노름에서 작은 시간 간격 조건이 필요하다.

실험 결과

연구 질문

• RQ1 이차 생산을 가지는 화학적 반발력 모델은 유계 영역에서 전역 약한 해를 갖는가?
• RQ2 약한 해가 유일하고 전역적으로 정칙해지기 위한 조건은 무엇인가? 특히 1차원 및 2차원에서의 경우에 대해.
• RQ3 이 모델에 대해 조건 없이 에너지 안정인 일계 시간 이산화 스킴을 구성할 수 있는가?
• RQ4 제안된 스킴들—특히 선형화된 버전에 대해—수렴성 및 오차 추정은 무엇인가?
• RQ5 스킴의 수치적 거동은 안정성, 에너지 감쇠 및 해의 프로파일 측면에서 어떻게 나타나는가?

주요 결과

• 이차 생산 $ f(u) = u^2 $를 가지는 화학적 반발력 모델은 유계 영역에서 전역 약한 해를 갖는다.
• 규칙성 기준 하에 1차원 및 2차원 영역에서 해의 유일성과 전역 정칙성이 확립된다.
• 비선형 역행 오일러 스킴은 조건 없이 에너지 안정이며, 세포 보존성도 유지되며, 약한 해로 수렴하고 오차 추정이 가능하다.
• 선형화된 결합 스킴은 조건 없이 에너지 안정이며 약한 해로 수렴하지만, 양의 보존성을 보장하지 못한다.
• 비선형 스킴의 유일성은 강한 노름에서 작은 시간 간격 조건 하에 증명되며, 이는 전역 강한 추정이 유도되는 $ nD $ 영역($ n \leq 2 $)에서는 작은 시간 간격 조건으로 축소된다.
• 수치 시뮬레이션은 두 스킴의 에너지 안정성과 수렴성 거동을 확인하여 이론 분석을 검증한다.