본 연구에서는 지표수 흐름 해석을 위한 2차원 유한요소모형을 개발하였다. 개발된 수치모형에서는 2차원 흐름 해석을 위한 천수방정식을 Galerkin법과 Newton-Raphson법에 의해 이산화하였으며 지형에 따라 삼각망과 사각망을 혼용하여 사용할 수 있도록 하였다. 구성된 대수방정식의 해는 유한요소법을 푸는데 매우 효율적인 프론탈 기법에 의해 구하였다. 개발된 모형을 두 개의 만곡부를 가지는 실험실 사행수로에 적용하여 만곡부에서의 횡방향 유속 및 수심 분포를 분석한 결과 만곡부에서 편수위 현상을 잘 재현함을 알 수 있었다. 또한 내측의 유속이 외측에 비해 빠르게 나타났고 유속분포는 두 만곡부의 중심을 기준으로 대칭적이었으며, 모의결과와 수리실험에 의한 실측값이 매우 잘 일치하였다. 본 모형을 완경사, 급경사, 역경사, 급격한 하상변화가 있는 수로와 위어를 포함한 수로에 적용하여 12개의 점변화류 수면곡선을 모의한 결과 수리학적으로 적합한 수면형상이 도출되었다.
A finite element model for analyzing surface water flow was developed. Shallow water equation was discretized and solved by Galerkin and Newton-Raphson method. Triangular or rectangular elements can be mixed together to construct meshes. The algebraic equation was solved by frontal method which is very efficient in finite element problem. The developed model was applied to rectangular meandering channel with two bends and transverse velocities and water depth distributions were examined. High velocity was located near the inner bank at the apexes of the bends and velocity distribution was symmetrical about the centerline at the midsection of two bend and super elevation also occurred. Simulation results showed very good agreement with measured data. Another numerical simulation was carried out in mild, steep, adverse and abrupt bottom change slope and channels with weir. 12 water surface profiles of gradually varied flow were correct in terms of hydraulic interpretation.
