단층 파동 모델 적분 이산화 연구

본 논문은 Camassa‑Holm 방정식의 단파 모델(SCHE)과 2차원 Toda 격자(2DTL) 사이의 깊은 연관성을 밝히고, 이를 토대로 연속, 반이산, 완전 이산 형태의 SCHE를 구축한다. 서론에서는 SCHE가 Harry‑Dym 계열에 속하며, κ=0 일 때는 Hunter‑Saxton 방정식과 동일함을 언급하고, κ≠0 일 때는 스케일 변환을 통해 일반화될 수 있음을 제시한다. 이어서 SCHE를 2DTL의 이중 감소 형태와 연결시키기 위해 Hirota D‑연산자를 이용한 bilinear식 (5)를 도입하고, τ‑함수를 Casorati 행렬식으로 정의한다. ψ_i는 두 개의 지수항으로 구성되며, 분산관계 ∂ψ/∂x₋₁=ψ_{n‑1}, ∂ψ/∂x₁=ψ_{n+1}을 만족한다. 두 감소 조건 p_i=−q_i 를 적용하면 τ_{n‑1}=Π p_i^{-2} τ_{n+1} 가 되고, τ₀=f, τ₁=g 로 두어 u=g/f 로 정의하면 bilinear식 (9)-(10)으로 변환된다. 여기서 ρ=u² 로 두고 hodograph 변환 x=2y−2(ln g)_s, t=s 를 수행하면 SCHE의 표준 형태 (∂_t+w∂_x)w_{xx}=2w_x(1−w_{xx}) 를 정확히 재현한다. 따라서 SCHE는 2DTL의 이중 감소에서 유도된다는 것이 핵심적인 결과이다. 다음으로 N‑cuspon 해를 제시한다. ψ_i를 (7)식의 구체적인 파라미터 a_i, p_i, q_i 로 정의하고, Casorati 행렬식 τ_n=g 로부터 w=−2(ln g)_{ss} 와 x=2y−2(ln g)_s 를 얻는다. p_i=−q_i 로 두면 N‑cuspon 파동이 형성되며, κ≠0 일 때는 x와 t에 κ 스케일을 적용해 일반화한다. 해의 정규성을 위해 τ‑함수가 양정인 조건을 강조한다. 구체적인 1‑cuspon과 2‑cuspon 예제가 그림 1에 제시되어 파라미터에 따른 파형 변화를 확인한다. 반이산화에서는 격자 인덱스 k 를 도입하고, ψ_i에 대한 차분·미분 관계 Δ_k ψ=ψ_{n+1}, ∂_s ψ=ψ_{n‑1} 를 만족하도록 설계한다. 이에 따라 반이산 2DTL의 Bäcklund 변환식 (22)를 이용해 bilinear식 (23)-(24)를 얻고, u_k=g_k/f_k 로 치환하면 (25)-(27)식이 도출된다. 여기서 discrete hodograph x_k=2ka−2(ln g_k)_s 와 δ_k=x_{k+1}−x_k 를 정의하고, ρ_k=u_k² 로 두면 연산을 정리해 (31)-(33)식, 즉 1/δ_k (w_{k+1}−w_k)−1/δ_{k‑1}(w_k−w_{k‑1}) = ½(δ_k+δ_{k‑1})−2a²(1/δ_k+1/δ_{k‑1}) dδ_k/dt = w_{k+1}−w_k 을 얻는다. 연속극한 a→0, δ_k→0 에서는 원 SCHE (4)를 정확히 회복한다. 반이산 해는 g_k와 f_k 를 각각 ψ_i의 Casorati 행렬식으로 정의하고, w_k=−2(ln g_k)_{ss} 로 표현한다. 또한 κ 파라미터를 포함한 형태 (37)와 그에 대응하는 N‑cuspon 해 (38)를 제시한다. 완전 이산화에서는 두 격자 변수 (k,l) 를 도입하고 ψ_i에 대한 차분·시간 차분 관계를 포함한 Casorati 행렬식 (39)를 구성한다. bilinear식 (40)-(41)와 두 감소 조건을 적용하면 네 개의 bilinear 방정식 (42)-(45)이 얻어진다. 이를 로그 미분하여 (46)-(49) 형태로 변환하고, w_{k,l}=−2(ln g_{k,l})_{ss}, x_{k,l}=2ka−2(ln g_{k,l})_s 로 정의한다. 격자 간 거리 δ_{k,l}=x_{k+1,l}−x_{k,l} 를 도입하고, 연산을 정리하면 (61)-(66) 두 차분식이 도출된다. 최종적으로는 Δ² w_{k,l}= (1/δ_{k,l}) M(δ_{k,l}−4a²δ_{k,l}) 와 δ_{k,l+1}−δ_{k,l}= b·