비수리 댐 붕괴 해석의 실제성: 비선형 얕은 물 방정식의 한계와 가능성

초록

본 연구는 전통적인 댐 붕괴 문제를 NSWE 검증 표준으로 삼는 관행에 의문을 제기한다. 두 유체 Navier‑Stokes 시뮬레이션과 NSWE의 해석·수치 해를 비교하여, 물 깊이 양성 가정이 깨지는 건조‑습윤 전이 구역에서 NSWE의 정확도가 급격히 저하됨을 확인한다. 결과는 NSWE가 중·대규모 홍수 예측에는 유용하지만, 초기 건조 지면 침투·파동 전파와 같은 현상은 별도 물리 모델이 필요함을 시사한다.

상세 분석

NSWE(Nonlinear Shallow Water Equations)는 수면 변위가 수심에 비해 작고, 수심이 양수이며 수직 속도분포가 거의 균일하다는 가정 하에 유도된다. 이러한 전제는 전통적인 댐 붕괴와 같은 급격한 수위 상승 상황에서도 어느 정도 타당성을 보이지만, 건조 지면을 침범하거나 물이 얕아지는 구역에서는 가정이 무너지게 된다. 논문은 두 가지 접근법을 동시에 수행한다. 첫 번째는 고전적인 NSWE를 이용해 초기 수위 차이와 중력 가속도만을 고려한 해석적 해와, 고해상도 유한체적·유한요소 기반의 2‑phase Navier‑Stokes(VOF) 시뮬레이션을 수행한 것이다. 두 번째는 동일한 초기 조건과 경계 조건을 적용해 두 모델의 시간·공간 전개를 정량적으로 비교한다.

비교 결과는 다음과 같이 요약된다. (1) 물이 충분히 깊고 흐름이 거의 수평인 구간에서는 NSWE가 Navier‑Stokes 결과와 거의 일치한다. 파동 전파 속도, 전단 응력, 에너지 보존 등 주요 거시량이 차이가 2 % 이하로 나타난다. (2) 물이 얕아지면서 수면이 건조 지면을 덮는 초기 단계에서는 NSWE가 물 깊이 음수 혹은 0이 되는 상황을 수학적으로 회피하기 위해 인위적인 ‘wetting‑drying’ 알고리즘을 도입한다. 이때 발생하는 인위적 압력 보정은 실제 물리적 압력 분포와 크게 달라져, 파동 전파 속도가 과대평가되고, 전단 응력이 과소평가된다. (3) Navier‑Stokes 시뮬레이션은 점성, 표면 장력, 공기와 물 사이의 인터페이스 전이 효과를 모두 포함한다. 특히 초기 붕괴 순간에 발생하는 급격한 가속도와 난류 발생은 NSWE가 전혀 포착하지 못한다. (4) 수치 실험에서 NSWE는 건조‑습윤 경계에서 ‘수위 음수’ 현상을 방지하기 위해 인위적인 최소 수심(예: 10⁻⁴ m)을 설정한다. 이 최소 수심이 실제 물리적 현상을 왜곡시키며, 특히 작은 규모(수십 미터 이하)의 댐 붕괴 시 시뮬레이션 결과에 큰 오차를 초래한다.

이러한 결과는 NSWE가 ‘대규모·중·장시간’ 홍수 예측에는 충분히 강력하지만, ‘초소형·단기·건조 지면 침투’와 같은 특수 상황에서는 보완 모델이 필요함을 명확히 보여준다. 논문은 또한 현재 NSWE 커뮤니티가 수치 스키마(고차 정확도, 보존형 스킴 등) 개발에 집중하고 있는 반면, 물리 모델 자체의 적용 범위 검증이 소홀히 이루어지고 있음을 지적한다.