[논문 리뷰] Solitary waves and $N$-particle algorithms for a class of Euler-Poincaré equations

Roberto Camassa, Dongyang Kuang외 1명|arXiv (Cornell University)|2014. 04. 18.

Nonlinear Waves and Solitons참고 문헌 18인용 수 4

한 줄 요약

이 논문은 비연속적이고 원뿔형 모양의 커널을 가진 일련의 옐러-파이오네퍼 방정식에 대해 N-입자 입자 알고리즘을 개발하며, 특정 초깃값 조건 하에서 이러한 시스템이 완전히 통합 가능한 두 입자 동역학을 유도할 수 있음을 보여준다. 이 연구는 결과적으로 두 차원 파동 동역학이 일차원 완전히 통합 가능한 얕은 물방정식으로 감소함을 밝히며, 단일 채널 해에서 여전히 전체 2차원 공간 의존성을 유지한다. 또한 고정밀 수치 비교를 통해 입자 방법이 허위 스펙트럼 오일러 방법과 높은 정확도로 일치함을 검증한다.

ABSTRACT

We study a class of partial differential equations (PDEs) in the family of the so-called Euler-Poincaré differential systems, with the aim of developing a foundation for numerical algorithms of their solutions. This requires particular attention to the mathematical properties of this system when the associated class of elliptic operators possesses non-smooth kernels. By casting the system in its Lagrangian (or characteristics) form, we first formulate a particles system algorithm in free space with homogeneous Dirichlet boundary conditions for the evolving fields. We next examine the deformation of the system when non-homogeneous "constant stream" boundary conditions are assumed. We show how this simple change at the boundary deeply affects the nature of the evolution, from hyperbolic-like to dispersive with a non-trivial dispersion relation, and examine the potentially regularizing properties of singular kernels offered by this deformation. From the particle algorithm viewpoint, kernel singularities affect the existence and uniqueness of solutions to the corresponding ordinary differential equations systems. We illustrate this with the case when the operator kernel assumes a conical shape over the spatial variables, and examine in detail two-particle dynamics under the resulting lack of Lipschitz-continuity. Curiously, we find that for the conically-shaped kernels the motion of the related two-dimensional waves can become completely integrable under appropriate initial data. This reduction projects the two-dimensional system to the one-dimensional completely integrable Shallow-Water equation [Camassa, R. and Holm, D. D., Phys. Rev. Lett., 71, 1961-1964, 1993], while retaining the full dependence on two spatial dimensions for the single channel solutions.

연구 동기 및 목표

• 비연속 타원형 연산자를 가진 일련의 옐러-파이오네퍼 방정식을 해결하기 위한 기초 수치 알고리즘을 수립하기 위해.
• 비균일한 '일정 유량' 경계 조건이 시스템의 분산성 및 쌍곡성 특성에 미치는 영향을 조사하기 위해.
• 커널이 리프시츠 연속성 없이 존재할 경우 N-입자 상미분계에서 해의 존재성과 유일성 분석하기 위해.
• 원뿔형 모양의 커널 하에서 두 입자 동역학의 통합 가능성과 그에 따른 파동 구조에 대한 영향을 탐색하기 위해.
• 비연속 커널에 대해 입자 기반 알고리즘과 허위 스펙트럼 오일러 방법 간의 성능 비교하기 위해.

제안 방법

• PDE에서 유도된 N-입자 유한차원 상미분계를 도출하기 위해 시스템을 라그랑주 형태(특성 기반)로 재정의한다.
• 자유 공간에서 동일한 딜레흐 경계 조건을 가진 입자 알고리즘을 개발한 후, 비균일한 '일정 유량' 경계 조건으로 확장한다.
• 커널은 원뿔형으로 가정되어 비리프시츠 연속성을 가지며, 이는 상미분계에서 해의 유일성에 영향을 미친다.
• 두 입자 동역학은 해밀토니안 역학을 사용하여 분석되며, 해밀토니안 H = c₁² + c₂²를 통해 정확한 해를 도출한다.
• 안정성과 두 입자 시스템에서 안정 다변량의 구조를 분석하기 위해 리아푸노프 함수를 구성한다.
• 6차 룬게-쿠타 방법을 사용하여 수치 검증을 수행하며, 상대 거리와 운동량에 대해 입자 알고리즘 해와 정확한 해석적 해를 비교한다.

실험 결과

연구 질문

• RQ1 비연속적이고 원뿔형 모양의 커널은 옐러-파이오네퍼 방정식의 N-입자 상미분계에서 해의 존재성과 유일성에 어떤 영향을 미치는가?
• RQ2 이러한 특이 커널에 의해 지배되는 두 입자 동역학은 완전히 통합 가능할 수 있으며, 만약 그렇다면 어떤 초깃값 조건에서 그러한 성질을 가지는가?
• RQ3 비균일한 '일정 유량' 경계 조건은 시스템의 분산성 대비 쌍곡성 특성에 어떤 영향을 미치는가?
• RQ4 입자 알고리즘이 비연속 커널을 가진 PDE 해를 허위 스펙트럼 오일러 방법 대비 얼마나 정확하게 시뮬레이션할 수 있는가?
• RQ5 두 차원 파동 동역학이 일차원 얕은 물방정식으로 감소하는 과정에서 단일 채널 해에서 전체 2차원 공간 의존성이 유지되는가?

주요 결과

• 원뿔형 커널의 경우, 적절한 초깃값 조건 하에서 두 입자 시스템은 완전히 통합 가능해지며, 정확한 해가 쌍곡함수로 표현될 수 있다.
• γ = 4c²인 정면 충돌의 정확한 해는 q(t) = -2 log[sech(ct)] 및 p(t) = ±2c / tanh(ct)로 주어지며, 솔리톤 유사 행동을 보인다.
• 6차 룬게-쿠타 방법을 사용한 수치 해는 정확한 해와 비교해 위치에 대해 2-norm 오차 2.9506e-05, 운동량에 대해 8.7397e-07을 기록한다.
• 두 입자 시스템의 안정 다변량은 z = √[(G_{b-1}(0) - G_{b-1}(q)) / (4H)]로 주어지며, b = 3/2일 때 얕은 물방정식의 알려진 형태로 감소한다.
• 입자 알고리즘은 높은 정확도로 작동하며, 표준 상미분 해 존재 정리가 실패하는 비연속 커널에 대해서도 이 방법이 타당함을 보여준다.
• 비균일한 경계 조건 하에서 시스템은 쌍곡성 유사 행동에서 분산성 행동으로 전이되며, 커널 특이성으로 인해 비트리비얼한 분산 관계가 도출된다.