변동 수심에서 층류 유체의 유한 진폭 파동 전파와 솔리톤 붕괴 메커니즘

초록

본 논문은 가변 계수를 갖는 Korteweg‑de Vries(KdV) 방정식을 이용해 두 층으로 구성된 변동 수심 해양에서 내부 솔리톤의 전파와 변형을 분석한다. 솔리톤은 두 개의 임계 지점에서 붕괴되며, 하나는 상하 층 두께가 동일해지는 지점, 다른 하나는 두 층 흐름이 단일 층 흐름으로 전환되는 지점이다. 하층 두께에 따른 솔리톤 진폭의 의존 관계를 도출하였다.

상세 분석

이 연구는 기존의 KdV‑계열 모델을 확장하여 수평·수직 변화를 동시에 고려한 일반화 Gardner 방정식(가변 계수형)을 도입한다. 먼저, 수직 밀도 분포를 Sturm‑Liouville 고유값 문제로부터 얻은 모드 함수 Φ(y)와 비선형 보정 함수 T(y)를 이용해 1+1 차원 비선형 파동 방정식(식 1)을 유도한다. 여기서 비선형 계수 α, β와 분산 계수 c는 각각 밀도 구배와 층 두께에 의존한다(식 5‑8). 특히 두 층 모델에서는 α₁(2차 비선형 계수)의 부호가 층 경계 위치에 따라 바뀌어 솔리톤의 극성(뾰족·오목)도 달라진다.

가변 수심을 포함한 경우, 수평 변화를 완만하게 가정하고 반사 효과를 무시하면 Gardner 방정식에 공간 의존 계수 ξ(x), Q(x) 등을 도입한 식(9‑12)를 얻는다. 이 방정식은 질량·에너지 보존법칙(식 13‑14)을 만족하므로, 솔리톤의 에너지 보존을 이용해 진폭 a(x)의 변화를 직접 계산할 수 있다.

솔리톤 해는 일반화된 형태(식 20‑21)로 표현되며, 파라미터 γ는 솔리톤 폭의 역수, B는 비선형 계수 비율을 나타낸다. α₁<0인 두 층 경우 B∈(0,1)이며, B→1일 때는 KdV 솔리톤(식 23)으로 수렴하고, B→0이면 “두꺼운” 솔리톤으로 폭이 무한히 커진다. 솔리톤 질량과 에너지는 각각 식(24‑25)와(27‑28)으로 정의되며, 진폭이 커질수록 단조적으로 증가한다.

주요 결과는 다음과 같다. (1) 솔리톤 진폭 a(x)는 하층 두께 h₂(x)와 초기 두께 h₂₀, 상층 두께 h₁의 비율 q=h₁/h₂₀에 따라 식(29)로 주어지며, h₂가 감소함에 따라 진폭이 감소하거나, q가 충분히 크면 진폭이 먼저 증가 후 감소하는 비단조적 거동을 보인다. (2) 진폭이 최대가 되는 깊이 h는 식(30)으로 구해지며, 최대 증폭 비율 a/a₀는 식(31)으로 표현된다. q>0.5이면 진폭은 지속적으로 감소하고, q<0.5이면 초기 증폭 후 급격히 소멸한다. (3) 두 임계점에서 비선형 계수 α₁이 0이 되거나 c가 사라지는 상황이 발생한다. 첫 번째 임계점(층 두께 동일)에서는 솔리톤이 소멸하고 반대 극성의 새로운 솔리톤이 생성된다(문헌 12, 13 참고). 두 번째 임계점에서는 두 층 흐름이 단일 층으로 전이되어 내부 파동 자체가 존재하지 않는다.

수치 예시(그림 2‑4)에서는 q값에 따른 최대 증폭 곡선과 하층 두께 변화에 따른 진폭 변화를 시각화하였다. 특히 h₁/h₂₀=0.1(얇은 상층)과 0.5(두꺼운 상층) 경우를 비교해, 얇은 상층에서는 진폭이 크게 증가 후 급격히 소멸함을, 두꺼운 상층에서는 진폭이 지속적으로 감소함을 확인한다. 또한, h₂→0일 때 진폭이 h₂^{-2/3} 비율로 발산하는 식(32)도 제시되어, 실제 해양에서 급경사 해저에 닿을 때 발생할 수 있는 “파동 충격” 현상을 예측한다.

이 연구는 실제 해양 수문 데이터(수직·수평 밀도·수심 변동)를 직접 입력해 계수를 자동 계산할 수 있는 소프트웨어(튜긴 등, 2011)와 연계될 수 있음을 강조한다. 따라서 현장 관측에 기반한 내부 파동 전파 예측, 해저 지형에 의한 파동 붕괴 위험 평가 등에 바로 적용 가능하다.