BBM‑유사 방정식의 여행파 해를 위에르스트라스 함수로 풀어내다

본 논문은 비선형 파동 방정식인 Benjamin‑Bona‑Mahony(BBM) 방정식의 일반화 형태인 B(m,n) : u_t+u_x+a(u^m)_x−(u^n)_{xxt}=0 (m,n>1, m≠n)를 대상으로, 여행파 해를 체계적으로 도출하는 새로운 방법론을 제시한다. 연구는 크게 네 단계로 진행된다. 1. **여행파 차원 축소** u(x,t)=φ(ξ), ξ= h x+ w t 라는 전형적인 여행파 가정을 도입한다. 이를 원 방정식에 대입하고 한 번 적분하면 φ에 대한 2차 비선형 상미분방정식 φ_{ξξ} − A φ − B φ^m + D=0 (A,B,D는 h,w,a,R에 의해 정의) 를 얻는다. 여기서 R는 적분 상수이며, A=h+ w h²/w, B=a h/w, D=R h²/w 로 명시된다. 2. **변수 변환 및 팩터화** φ^n(ξ)=W(θ), ξ=θ 로 치환하면 d²W/dθ² − A W^{1/n} − B W^{m/n} + D=0 가 된다. 이는 β=0 인 형태 d²W/dθ²+F(W)=0 와 동일하므로, (d/dθ−f₂)(d/dθ−f₁)W=0 라는 팩터화 형태를 적용한다. 일관성 조건 f₁f₂=F(W)·W와 f₂+∂(Wf₁)/∂W=0 로부터 f₁(W) 를 구하면 f₁(W)=±√{A W^{1+1/n}+B W^{1+m/n}+D W^{-1}+C} 가 된다. 여기서 C는 적분 상수이며, 최종적으로 1차 ODE dW/dθ = ±√{A W^{1+1/n}+B W^{1+m/n}+D W^{-1}+C} 를 얻게 된다. 3. **다항식 형태로 변환 및 ℘ 함수 해법** 해의 일반성을 확보하기 위해 W=φ^p (p≠0,1) 로 치환한다. 그러면 (dφ/dθ)² = α₀ + α₁ φ^{k₁}+α₂ φ^{k₂}+α₃ φ^{k₃} 형태의 다항식이 된다. 여기서 지수 k_i 는 n, m, p 에 따라 정수(0~4) 로 제한된다. 논문은 가능한 (p,m) 조합을 체계적으로 열거하고, 각각에 대해 3차 혹은 4차 다항식 P(φ) 를 도출한다. 다항식 P(φ) 의 해는 위에르스트라스 ℘ 함수의 표준 해법을 이용한다. 일반적인 형태 (dφ/dθ)² = a₀φ⁴+4a₁φ³+6a₂φ²+4a₃φ+a₄ 로부터 불변량 g₂, g₃ 를 계산하고, ℘(θ;g₂,g₃) 로 해를 표현한다. 특히, g₂=g₃=0 (레미네식), g₂<0,g₃=0 (의사‑레미네식), g₂=0,g₃>0 (등조화) 등 특수 경우를 상세히 다루어, 각 경우에 대한 명시적 해를 ℘ 함수와 그 변형으로 제시한다. 예를 들어, m=n+½, p=−2n/(1−n) 인 경우에는 (dφ/dθ)² = A + B φ 형태가 되어 ℘(θ;0,0)=1/θ² 로부터 단순한 포물선 해를 얻는다. 다른 경우에는 ℘ 함수의 주기적 특성에 따라 정규 파동, 특이 파동, 혹은 컴팩트 서포트 파동이 나타난다. 4. **라그랑지안·해밀토니안 구조와 에너지 보존** 얻어진 2차 ODE는 라그랑지안 L = ½(W_θ)² + A/(n+1)W^{1+1/n} + B/(m+1)W^{1+m/n} − D W 로부터 유도된다. 정준 모멘텀 P = ∂L/∂W_θ = W_θ 를 정의하고, 해밀토니안 H = P W_θ − L = ½P² − A/(n+1)W^{1+1/n} − B/(m+1)W^{1+m/n} + D W 를 얻는다. θ 가 명시적으로 나타나지 않으므로 H는 보존량이며, H=E 로 두 상수 I_± 로 분해한다: I_± = W_θ ∓ √{A W^{1/n}+B W^{m/n}−D} · e^{±S(θ)}. 이 구조는 팩터화 과정에서 도출된 1차 ODE 와 일치한다. **결과 및 의의** - 팩터화 기법은 기존의 직접 적분법, tanh‑method, 사인‑코사인 방법보다 더 넓은 파라미터 영역(m≠n, 정수·비정수 모두)에서 해를 제공한다. - 위에르스트라스 ℘ 함수의 특수 형태(레미네식, 의사‑레미네식, 등조화)를 이용해 해를 간결하게 표현함으로써, 파라미터에 따른 파동 형태(주기적, 특이, 컴팩트 서포트)를 명확히 구분한다. - 라그랑지안·해밀토니안 분석을 통해 에너지 보존과 팩터화된 1차 ODE 사이의 깊은 연관성을 밝히며, 비선형 시스템의 보존 구조를 이해하는 데 기여한다. - 수치 플롯(그림 1~5)으로 제시된 예시들은 이론적 해가 실제 파동 형태와 일치함을 시각적으로 확인한다. 전반적으로, 본 연구는 BBM‑유사 방정식의 m≠n 경우에 대한 해를 체계적이고 일반적인 방법으로 제공함으로써, 비선형 파동 이론 및 응용 분야(예: 해양 파동, 비선형 광학)에서 새로운 분석 도구를 제시한다.