계단형 바닥과 외력 하 얕은 물 흐름을 위한 준이중층 유한체적 스키마

** 본 논문은 비균일한 바닥과 외력이 동시에 작용하는 얕은 물 흐름을 정확히 모사하기 위한 새로운 수치 스키마를 제시한다. 서론에서는 얕은 물 방정식이 3차원 유체 방정식의 수직 평균에 의해 도출되지만, 복잡한 바닥 형상과 외력(예: 코리올리 힘) 하에서는 수직 비균일성이 무시될 수 없으며, 기존의 유한체적·Godunov‑type 방법들은 부분 침수·건조 영역을 처리하는 데 한계가 있음을 지적한다. 문헌 검토에서는 Roe‑type, 하이드로스태틱 재구성, LeVeque‑법, 표면 기울기법 등 다양한 웰밸런스 스키마와 외력 처리 방법을 소개하고, 특히 단계형 바닥에 대한 정확한 리만 문제 해법이 부족함을 강조한다. 2절에서는 얕은 물 방정식에 외력 E를 포함한 일반 형태를 제시하고, 외력을 ‘가상의 시간‑의존성 바닥(k(x,t))’으로 변환하는 수학적 변환을 전개한다. 이 변환을 통해 외력 항을 ∂k/∂x 형태로 재작성함으로써, 기존의 소스 항과 동일한 구조를 유지하면서도 바닥이 수행하는 기계적 일을 명시적으로 나타낼 수 있다. 3절에서는 ‘준이중층’ 모델의 핵심 아이디어를 상세히 설명한다. 단계형 바닥 근처에서 흐름을 상층(실제 물)과 하층(가상 물)으로 분리하고, 하층은 리만 문제 설정을 위한 보조 변수로만 사용한다. 하층 두께는 단계 높이와 상층 흐름 속도에 따라 동적으로 결정되며, 건조 구역이 발생하면 하층 두께가 0이 되어 물이 수직 벽을 적시지 못하는 ‘폭포 현상’을 자연스럽게 재현한다. 유한체적 스키마는 셀 중심값을 보존 변수(깊이 h, 모멘텀 hu, hv)로 두고, 경계면에서 flux를 계산한다. 경계에 단계가 존재하면, ‘준이중층 Riemann 솔버’를 적용해 상층과 하층의 상태를 각각 고려한 특성값을 구한다. 이때 추가된 소스 항은 바닥 기울기와 가상의 하층 두께에 의해 결정되며, 이는 바닥이 흐름에 수행하는 일(∫g k dx)으로 해석된다. 특히 부분 침수·건조 셀을 처리하기 위해 셀을 ‘완전 침수’, ‘부분 침수’, ‘완전 건조’ 세 종류로 구분하고, 각 경우에 맞는 Riemann 해와 소스 항을 적용한다. 이를 통해 전통적인 스키마가 겪는 수치 발산이나 비물리적 음수 깊이 문제를 회피한다. 4절에서는 기존의 웰밸런스 스키마와 제안된 스키마를 비교한다. 기존 방법은 바닥 기울기와 외력 항을 단순히 평균하거나 보정하는 방식으로, 건조 구역에서 에너지 손실을 과소평가한다. 반면 제안된 방법은 단계별 기계적 일을 정확히 계산해 에너지 보존성을 향상시킨다. 5절에서는 다섯 가지 시험 사례를 통해 스키마의 성능을 검증한다. (1) 경사면 위 직사각형 물기둥 붕괴 실험: 단계별 흐름 전이와 물량 보존이 정확히 재현된다. (2) 코리올리 힘을 포함한 대규모 회전 흐름: 잠재 와도 보존이 향상되고, 회전 대류 구조가 정확히 포착된다. (3) 2‑차원 경사면 댐‑브레이크 실험: 수위·속도 프로파일이 실험 데이터와 높은 일치도를 보인다. (4) 복합 삼각함수형 바닥 위 대규모 흐름: 복잡한 지형에서도 안정적인 계산이 가능하다. (5) 해안에 장애물을 둔 쓰나미 충돌 시뮬레이션: 건조 구역과 부분 침수 영역이 동적으로 변화하며, 파고 감쇠와 흐름 분산이 실험적 관측과 일치한다. 각 사례에서 제안된 스키마는 기존 웰밸런스 방법보다 평균 절대 오차가 30 % 이상 감소했으며, 특히 건조 구역 전이와 폭포 현상에서 큰 차이를 보였다. 6절에서는 연구 결과를 요약하고, 준이중층 모델이 비균일 바닥과 외력 하에서 얕은 물 흐름을 물리적으로 충실하게 재현할 수 있음을 강조한다. 또한, 향후 연구 방향으로 3차원 확장, 고차 정확도 재구성, 실시간 홍수 예측 시스템 적용 등을 제시한다. **