Abstract
In this study, Galerkin scheme and SU/PG scheme of Petrov-Galerkin family were applied to the shallow water equations and a finite element model for shallow water flow was developed. Numerical simulations were conducted in several flumes with convection-dominated flow condition. Flow simulation of channel with slender structure in the water course revealed that Galerkin and SU/PG schemes showed similar results under very low Fr number and Re number condition. However, when the Fr number increased up to 1.58, Galerkin scheme did not converge while SU/PG scheme produced stable solutions after 5 iterations by Newton-Raphson method. For the transcritical flow simulation in diverging channel, the present model predicted the hydraulic jump accurately in terms of the jump location, the depth slope, and the flow depth after jump, and the numerical results showed good agreements with the hydraulic experiments carried out by Khalifa(1980). In the oblique hydraulic jump simulation, in which convection-dominated supercritical flow (Fr=2.74) evolves, Galerkin scheme blew up just after the first iteration of the initial time step. However, SU/PG scheme captured the boundary of oblique hydraulic jump accurately without numerical oscillation. The maximum errors quantified with exact solutions were less than 0.2% in water depth and velocity calculations, and thereby SU/PG scheme predicted the oblique hydraulic jump phenomena more accurately compared with the previous studies (Levin et al., 2006; Ricchiuto et al., 2007).
본 연구에서는 천수방정식에 Galerkin법과 Petrov-Galerkin 기법의 일종인 SU/PG 기법을 적용하여 유한요소 천수흐름 해석 모형을 개발하고, 다양한 실험수로에서 이송이 지배적인 흐름을 수치 모의하였다. 수로 내부에 얇은 판 형태의 구조물이 존재하는 경우 Fr 수와 Re 수가 매우 낮은 흐름에서는 Galerkin 기법과 SU/PG 기법이 동일한 결과를 나타냈으나, Fr=1.58인 경우 Galerkin법은 발산하여 해를 얻을 수 없었다. 이 경우 SU/PG법은 Newton-Raphson법에 의한 5회의 반복에 의해 수렴된 유속결과를 구할 수 있었다. 사류와 상류가 혼재하여 천이류가 나타나는 단면확대 수로 모의에서 SU/PG 기법을 이용한 본 연구의 경우 상류단 수심조건이 잘 유지되며, 도수가 발생하는 지점 및 도수에 의한 수심 경사, 도수 후의 수심이 모두 Khalifa(1980)의 실험결과와 매우 근사하였다. 이송이 지배적인 사류(Fr=2.74)에 의한 사각도수 모의의 경우에도 Galerkin 기법은 최초 모의시간의 첫 번째 반복 이후 발산하였으나, SU/PG 기법은 도수 경계면을 수치진동 없이 잘 포착하였으며, 해석해와 비교한 수심 및 유속의 최대 오차는 0.2% 이내로 나타나 기존 연구(Levin 등, 2006; Ricchiuto 등, 2007)에 비해 더욱 정확한 결과를 도출하였다.