(b, s) 좌표로 표현된 천수방정식에 유한체적법을 사용하여 해상 및 해저 산사태 수치모의

Numerical Simulation of Subaerial and Submarine Landslides Using the Finite Volume Method in the Shallow Water Equations with (b, s) Coordinate

산사태의 시간에 따른 전파를 모의하기 위해서 천수방정식을 사용하여 산사태 수치모형을 개발하였다. 하천 및 해양에서의 산사태에 모두 해석이 가능하도록 (b, s) 좌표로 표현된 천수방정식을 개발하였다. 산사태에서 발생하는 수치적인 불연속성을 극복하기 위해서 HLL approximate Riemann solver와 total variation diminishing (TVD) limiter를 사용한 유한체적법을 사용하였다. 댐파괴 흐름와 토석류의 각 경우에 수치해석을 수행한 결과를 해석해와 실험자료와 비교를 하였다. 그 결과 서로 유사함을 확인되었다. 본 모형을 사용하여 해상 산사태와 해저 산사태를 성공적으로 모의하였다. 해저 산사태에 비해 해상 산사태의 전파속도가 더 빠르고, 바닥경사가 급할수록 또는 거칠기가 작을수록 산사태 전파속도가 더 빨라짐을 확인하였다.

A model of landslides is developed using the shallow water equations to simulate time-dependent performance of landslides. The shallow water equations are derived using the (b, s) coordinate system which can be applied in both river and ocean. The finite volume scheme employing the HLL approximate Riemann solver and the total variation diminishing (TVD) limiter is applied to deal with the numerical discontinuities occurring in landslides. For dam-break water flow and debris flow, numerical results are compared with analytical solutions and experimental data and good agreements are observed. The developed landslide model is successfully applied to predict subaerial and submarine landslides. It is found that the subaerial landslide propagates faster than the submarine landslide and the speed of propagation becomes faster with steeper bottom slope and less bottom roughness.