[논문 리뷰] An analytic technique for the solutions of nonlinear oscillators with damping using the Abel Equation

A. Ghose‐Choudhury, Partha Guha|arXiv (Cornell University)|2016. 08. 08.

Fractional Differential Equations Solutions인용 수 3

한 줄 요약

이 논문은 일반화된 리카티 방정식을 일阶 아벨 방정식으로 감소시키고 치엘리니 적분 가능성 조건을 적용하여 비선형 감쇠 진동자 클래스를 해석적으로 해결하는 방법을 제시한다. α ≥ 2√2(n+1) 인 시스템  ddot{x} + α x^{2n+1} dot{x} + x^{4n+3} = 0 에 대해 명시적인 닫힌 형태의 해가 유도되며, 이는 이전의 수치적 추측인 주기적에서 비주기적 운동으로의 전이 및 같은 시스템에서의 이소크로노니를 해석적으로 확인한다.

ABSTRACT

Using the Chiellini condition for integrability we derive explicit solutions for a generalized system of Riccati equations $\ddot{x}+αx^{2n+1}\dot{x}+x^{4n+3}=0$ by reduction to the first-order Abel equation assuming the parameter $α\ge 2\sqrt{2(n+1)}$. The technique, which was proposed by Harko extit{et al}, involves use of an auxiliary system of first-order differential equations sharing a common solution with the Abel equation. In the process analytical proofs of some of the conjectures made earlier on the basis of numerical investigations in \cite{SJKB} is provided.

연구 동기 및 목표

• 비선형 감쇠 진동자에서 주기적 운동에서 비주기적 운동으로의 전이에 대한 수치적 추측을 해석적으로 증명하는 것.
• 치엘리니 적분 가능성 조건을 일반화된 리카티 형 방정식을 포함한 고차 비선형 시스템으로의 적용 범위를 확장하는 것.
• 비선형 감쇠와 강한 비선형성의 복원력이 있는 비선형 진동자 계열에 대해 명시적인 닫힌 형태의 해를 유도하는 것.
• Lié nard 형 시스템  ddot{x} + (2n+3)x^{2n+1} dot{x} + x^{4n+3} + w_0^2x = 0 에서의 이소크로노스 진동의 추측을 해석적 방법으로 검증하는 것.

제안 방법

• 변수 치환  dot{x} = 1/v 를 통해 이阶 비선형 진동자 방정식을 첫째 순서 아벨 방정식의 첫 종류로 감소시키는 것.
• 적분 인자  exp(−∫ g_2(x)dx) 를 사용하여 아벨 방정식을 표준형  du/dx = A(x)u^2 + B(x)u^3 으로 변환하는 것.
• 고정된 상수 s 에 대해  d/dx(B/A) = sA(x) 인 치엘리니 적분 가능성 조건을 적용하여 적분에 의한 정확한 해법 가능성을 보장하는 것.
• α ≥ 2√2(n+1) 조건 하에서 감소된 아벨 방정식의 적분을 통한 명시적 해 유도.
• 아벨 방정식과 공통 해를 가지는 보조 1차 상미분방정식계를 사용하여 적분을 용이하게 하는 것.
• 변수 분리 및 변환 기법을 통해 치엘리니 조건을  ddot{x} + f(x) dot{x}^{n+1} + g(x) dot{x}^n = 0 형식의 고차 리에나르 방정식으로 일반화하는 것.

실험 결과

연구 질문

• RQ1 비선형 진동자  ddot{x} + α x^{2n+1} dot{x} + x^{4n+3} = 0 에서 주기적 운동에서 비주기적 운동으로의 전이의 해석적 메커니즘은 무엇인가?
• RQ2 Lié nard 형 시스템  ddot{x} + (2n+3)x^{2n+1} dot{x} + x^{4n+3} + w_0^2x = 0 에서의 이소크로노스 진동에 대한 수치적 추측을 엄밀하게 증명할 수 있는가?
• RQ3 치엘리니 적분 가능성 조건은 리에나르 형 방정식의 고차 비선형 감쇠 항으로 어떻게 확장되는가?
• RQ4 임계 감쇠 임계값  α_c = 2√2(n+1) 이 시스템의 정성적 행동을 결정하는 데 어떤 역할을 하는가?

주요 결과

• α ≥ 2√2(n+1) 인 경우, 시스템  ddot{x} + α x^{2n+1} dot{x} + x^{4n+3} = 0 은 아벨 방정식으로의 감소와 치엘리니 조건 적용을 통해 명시적인 닫힌 형태의 해를 가진다.
• 임계 값  α_c = 2√2(n+1) 는 아벨 방정식의 삼차항의 판별식이 0이 되는 임계점이며, 주기적에서 비주기적 운동으로의 전이를 알리는 신호이다.
• α = 2n+3 이고 w_0 ≠ 0 이면, 해는 주기적이고 이소크로노스이며, 리에나르 시스템에서의 이소크로노니에 대한 사바티니 기준을 만족함으로써 확인된다.
• α = 2n+3 이고 w_0 = 0 이면, 해는 특이점이며 유한 시간 내에 폭발하며, 명시적 형태로  x ∝ [(t_0 ∓ t)^{-1}]^{1/(2n+1)} 으로 표현된다.
• α = 2n+3, w_0 = 1 인 경우의 해는 기저 함수  {}_2F_1 로 표현되며, 이는 해석적 접근 가능성의 확인을 제공한다.
• 일반화된 치엘리니 조건  ((g/f)') = K (g^n / f^{n-1}) 는  dot{x}^{n+1} 과  dot{x}^n 감쇠 항을 포함하는 고차 리에나르 방정식의 정확한 적분을 가능하게 한다.