[논문 리뷰] Numerical study of blow-up in solutions to generalized Korteweg-de Vries equations

Christian Klein, Ralf Peter|arXiv (Cornell University)|2013. 07. 02.

Advanced Mathematical Physics Problems참고 문헌 4인용 수 19

한 줄 요약

이 논문은 비선형성의 임계($n=4$) 및 초임계($n=5$)인 일반화된 코르티제-데브리스(gKdV) 방정식에서 폭발 현상에 대한 수치적 조사 결과를 제시한다. 동적 스케일링과 스펙트럼 방법을 사용하여, 해가 마르텔, 메를, 라파엘에 의해 예측된 것과 일치하는 자세한 폭발 프로파일을 확인하였으며, 폭발 시간 $t^*$는 분산 매개변수 $\epsilon$에 대해 지수적 의존성을 보이며, $\epsilon \to 0$의 극한에서 허프 임계 시간보다 큰 최소 폭발 시간 $t_0^* \approx 1.0536$를 가진다. 본 연구는 해의 질량이 솔리톤 질량보다 훨씬 큰 경우에도 폭발 역학의 보편성을 확인한다.

ABSTRACT

We present a detailed numerical study of solutions to general Korteweg-de Vries equations with critical and supercritical nonlinearity. We study the stability of solitons and show that they are unstable against being radiated away and blow-up. In the $L_{2}$ critical case, the blow-up mechanism by Martel, Merle and Raphaël can be numerically identified. In the limit of small dispersion, it is shown that a dispersive shock always appears before an eventual blow-up. In the latter case, always the first soliton to appear will blow up. It is shown that the same type of blow-up as for the perturbations of the soliton can be observed which indicates that the theory by Martel, Merle and Raphaël is also applicable to initial data with a mass much larger than the soliton mass. We study the scaling of the blow-up time $t^{*}$ in dependence of the small dispersion parameter $ε$ and find an exponential dependence $t^{*}(ε)$ and that there is a minimal blow-up time $t^{*}_{0}$ greater than the critical time of the corresponding Hopf solution for $ε o0$. To study the cases with blow-up in detail, we apply the first dynamic rescaling for generalized Korteweg-de Vries equations. This allows to identify the type of the singularity.

연구 동기 및 목표

• 비선형성의 $n \geq 4$인 일반화된 코르티제-데브리스(gKdV) 방정식의 해에 대한 안정성과 폭발 역학을 조사한다. 특히 $L^2$-임계($n=4$) 및 초임계($n=5$) 영역에서 조사한다.
• 마르텔, 메를, 라파엘에 의해 제안된 $L^2$-임계 gKdV 방정식에 대한 이론적 폭발 메커니즘을 수치적으로 검증한다.
• 작은 분산 매개변수 $\epsilon$에 대한 폭발 시간 $t^*$의 의존성, 특히 작은 $\epsilon$ 근처에서의 행동을 분석한다.
• 솔리톤의 섭동에서 관찰된 동일한 폭발 프로파일이 솔리톤 질량보다 훨씬 큰 질량을 가진 초기 자료에 대해서도 나타나는지 확인한다.

제안 방법

• gKdV 방정식을 폭발 특이점의 수치적 추적을 가능하게 하는 형태로 변환하기 위해 동적 스케일링을 적용한다.
• 디알리아징을 적용한 의사-스펙트럼 푸리에 방법을 사용하여 공간에서 gKdV 방정식을 해석함으로써 진동 및 폭발 역학을 고정밀도로 해석한다.
• 해가 폭발에 접근함에 따라 정확도를 유지하기 위해 적응형 시간 스텝 및 해상도 제어를 구현하며, 특히 해상도 손실을 감지하기 위해 푸리에 계수의 크기를 모니터링한다.
• 로그 노름($\ln||u_x||_2^2$, $\ln||u||_\infty$)을 $\alpha \ln(t^* - t) + \beta$에 적합하여 $t^*$를 추정하고 자세한 폭발 행동의 타당성을 검증한다.
• 최소 제곱법을 사용해 $\ln(t^*/t_c)$를 $\gamma \epsilon + \delta$에 적합함으로써 $\epsilon \to 0$에서의 $t^*$의 점근적 행동을 규명하며, 지수적 스케일링을 확인한다.
• 다양한 $\epsilon$ 값을 가진 $n=4$ 및 $n=5$에 대해 $u_0 = \text{sech}^2 x$의 초기 자료를 수치 시뮬레이션하여 산란 충격 형성과 폭발 순서를 관찰한다.

실험 결과

연구 질문

• RQ1 마르텔, 메를, 라파엘이 $L^2$-임계 gKdV 방정식($n=4$)에 대해 예측한 폭발 메커니즘이, 반드시 솔리톤에 가까운 초기 자료가 아니어도 수치적으로 해의 폭발에 나타나는가?
• RQ2 작은 분산 매개변수 $\epsilon$에 대한 폭발 시간 $t^*$는 $n=4$ 및 $n=5$의 경우에 어떻게 의존하는가? 그리고 $\epsilon \to 0$일 때 지수적 스케일링을 보이는가?
• RQ3 솔리톤 섭동에서 관찰된 동일한 자세한 폭발 프로파일이 솔리톤 질량보다 훨씬 큰 질량을 가진 초기 자료의 해에서도 관찰되는가? 이는 폭발 역학의 보편성을 시사하는가?
• RQ4 작은 $\epsilon$ 영역에서 항상 산란 충격이 폭발 이전에 형성되는가? 이는 허프 해의 임계 시간 $t_c$에 비해 폭발 시간에 어떤 영향을 미치는가?
• RQ5 $\epsilon \to 0$ 극한에서 최소 폭발 시간 $t_0^*$는 해당 허프 해의 임계 시간 $t_c$보다 엄밀히 큰가?

주요 결과

• 모든 $n=4$의 경우, 마르텔, 메를, 라파엘이 예측한 자세한 폭발 프로파일이 수치적으로 확인되었으며, $t \to t^*$일 때 $||u_x||_2 \sim ||Q'||_2 / (l_0 (t^* - t))$이며, 해는 스케일링된 솔리톤 프로파일로 수렴한다.
• 모든 $n=5$의 경우, 폭발 역학이 $n=4$의 경우와 동일한 보편적 프로파일을 따르는 것으로 수치적으로 관찰되었으며, 이는 마르텔, 메를, 라파엘의 이론이 솔리톤 질량보다 훨씬 큰 질량을 가진 초기 자료에도 적용된다는 것을 시사한다.
• 모든 $n=4$ 및 $n=5$의 경우 폭발 시간 $t^*$는 $\epsilon$에 대해 지수적 의존성을 보이며, $\ln(t^*/t_c)$는 $\gamma \epsilon + \delta$에 선형적으로 적합되며 상관계수 $r=0.9999$를 보이며, 이는 $\lim_{\epsilon \to 0} t^* = t_0^* \approx 1.0536 > t_c \approx 0.5341$임을 의미한다.
• 작은 $\epsilon$ 근처에서 항상 산란 충격이 폭발 이전에 형성되며, 충격 구조 내에서 처음 나타나는 솔리톤이 결국 폭발하는 데 기여한다.
• $\ln||u_x||_2^2$ 및 $\ln||u||_\infty$를 $\alpha \ln(t^* - t) + \beta$에 적합한 결과, 다양한 노름 간 일致한 $t^*$ 추정치를 얻었으며, 스케일링된 시간에서의 지수적 스케일링 덕분에 $n=5$의 경우 $n=4$보다 더 좋은 수렴성을 보였다.
• $\epsilon = 0.01$일 때 폭발 시간은 $||u_x||_2^2$에 대해 $t^* \approx 1.1475$, $||u||_\infty$에 대해 $t^* \approx 1.1434$로 추정되었으며, 이는 일致성과 점근적 영역의 시작을 확인한다.