[논문 리뷰] Can a traveling wave connect two unstable states? The case of the nonlocal Fisher equation

Grégoire Nadin, Benoı̂t Perthame외 1명|arXiv (Cornell University)|2010. 11. 20.

Mathematical and Theoretical Epidemiology and Ecology Models참고 문헌 6인용 수 6

한 줄 요약

이 논문은 비국소 피셔 방정식에서 이동파 해를 조사한다. 여기서 안정 상태 $ u=0 $ 와 $ u=1 $ 모두 동역학적으로 불안정하다. $ 0 $ 은 단안정성 동역학으로 인해, $ 1 $ 은 투링 불안정성으로 인해 불안정하다. 유한 도메인에서 인공 경계 조건을 사용한 새로운 수치 알고리즘을 개발하여 이동파를 계산함으로써, 저자들은 $ 1 $ 이 투링 불안정하더라도 단조 증가 이동파가 여전히 $ 0 $ 에서 $ 1 $ 으로 연결될 수 있음을 보여주며, 비국소성에 의해 유도된 불안정성 하에서 이러한 파동이 존재하지 않을 것이라는 기대에 도전한다.

ABSTRACT

This note investigates the properties of the traveling waves solutions of the nonlocal Fisher equation. The existence of such solutions has been proved recently in \cite{BNPR} but their asymptotic behavior was still unclear. We use here a new numerical approximation of these traveling waves which shows that some traveling waves connect the two homogeneous steady states $0$ and $1$, which is a striking fact since $0$ is dynamically unstable and $1$ is unstable in the sense of Turing.

연구 동기 및 목표

• 비국소 반응-확산 시스템에서 두 불안정한 안정 상태를 연결하는 이동파가 존재하는지 여부에 대한 열린 문제를 해결하기 위해.
• 안정 상태 $ u=1 $ 가 투링 불안정한 경우 비국소 피셔 방정식에서 이동파 해의 점근적 행동을 명확히 하기 위해.
• 비국소 상호작용과 불안정한 경계 조건에도 불구하고 높은 정확도로 이동파 해를 계산할 수 있는 수치 방법을 개발하고 검증하기 위해.
• 특히 $ \mu $ 가 클 경우, $ 1 $ 이 투링 불안정할 때 파동이 $ u=1 $ 에 연결되는지 아니면 안정한 주기적 상태에 연결되는지 확인하기 위해.
• 비국소 효과로 인해 이동파의 단조성 파괴가 발생하는 임계 임계값 $ \mu_c \approx 8.9 $ 를 조사하기 위해.

제안 방법

• 유한 도메인 $[x_l, x_r]$ 에서 인공 경계 조건 $ v(x_l) = 1 $, $ v(x_r) = 0 $, $ v(0) = \epsilon $ 을 갖는 이동파 방정식을 해결하기 위한 수치 알고리즘을 개발하였다.
• 계산 도메인을 두 개의 구간으로 나누었다: $ I_1 = [x_l, 0] $ 과 $ I_2 = [0, x_r] $, 각각 딜리클레 조건을 갖는 두 개의 별도 타입의 타원형 문제를 풀었다.
• 컨볼루션 항 $ \phi \star v $ 는 $ v $ 를 $ (-\infty, x_l) $ 에서 1로, $ (x_r, \infty) $ 에서 0으로 연장하여 비국소 평가를 가능하게 하였다.
• 해의 연속성을 확보하기 위해 $ x=0 $ 에서 도함수의 점프를 최소화하도록 슈팅 방법을 사용하여 파동 속도 $ \sigma^L $ 를 조정하였다.
• 수렴성은 공간 그리드 $ \Delta x $ 를 세밀하게 하고 도메인 길이 $ L $ 을 늘림으로써 평가되었으며, 오차는 $ x=0 $ 에서 $ |E(0)| $ 를 통해 측정되었다.
• 알고리즘은 균일 커널 (5), 비음성 및 컴 pact 지지 집합 커널 (6a), (6b), 가우시안 커널 (7) 등 다양한 커널에 대해 테스트하여 강건성과 커널 의존성의 정도를 평가하였다.

실험 결과

연구 질문

• RQ1 비국소 피셔 방정식에서 동역학적으로 불안정한 상태 $ u=0 $ 와 투링 불안정한 상태 $ u=1 $ 를 연결하는 이동파 해가 존재할 수 있는가?
• RQ2 상태 $ u=1 $ 에서 투링 불안정성이 존재할 경우, 단조 증가 이동파의 존재를 방해하는가, 아니면 여전히 존재할 수 있는가?
• RQ3 비국소 효과로 인해 이동파의 단조성이 상실되는 데에 해당하는 $ \mu $ 의 임계 값은 무엇인가?
• RQ4 특히 푸리에 변환값이 음수인 커널들, 예를 들어 $ \sin \xi / \xi $ 와 같은 커널의 선택이 파동의 형태와 안정성에 어떤 영향을 미치는가?
• RQ5 컴 pact 지지 집합 초기 조건을 갖는 코시 문제의 해들이 이동파에 영향을 받는가, 아니면 펄레이팅 프론트와 같은 더 복잡한 구조와 공존하는가?

주요 결과

• 비국소 피셔 방정식에서 $ u=1 $ 이 투링 불안정하더라도, $ \mu=1000 $ 과 $ \mu=2500 $ 인 커널 (5) 에 대해 수치적으로 $ u=0 $ 에서 $ u=1 $ 으로 연결되는 이동파 해가 존재함을 보여주었다.
• $ \mu=10 $ 인 경우, 파동은 거의 단조 증가이며 최댓값이 $ 1.0028 $ 으로 약간의 비단조성성을 보이며, $ \mu=10 $ 에서는 최댓값이 약간 $ 1 $ 을 초과한다.
• 비단조성성 상실의 임계 값은 수치적으로 $ \mu_c \approx 8.9 $ 로 확인되었으며, $ x=-\infty $ 근처의 선형화로부터 유도된 분석적 추정과 일치한다.
• 가장 $ \mu $ 가 증가함에 따라 파동 꼬리 부분에 감쇠 진동이 나타나며, 이러한 진동의 진폭은 도메인 길이 $ L $ 이 증가함에 따라 천천히 감소함을 보이며, $ 1 $ 에서 $ 0 $ 으로 연결되는 파동로 수렴하는 것으로 나타났다.
• 모든 테스트된 커널들—특히 푸리에 변환값이 음수인 커널들(예: $ \sin \xi / \xi $)과 가우시안 커널—에 대해 이동파 해는 여전히 $ u=1 $ 에 연결되어 있으며, 주기적 상태로 연결되지 않았다.
• 수치 알고리즘은 $ \Delta x $ 가 감소하고 $ L $ 이 증가함에 따라 오차 $ |E(0)| $ 가 감소하는 것을 성공적으로 계산하였으며, $ \mu=64 $ 에서 $ \sigma^L \to 15.367 $ 로 수렴하며 이론적 속도 $ 2\sqrt{64} = 16 $ 에 가까워졌다.