Abstract
This paper presents a computational simulation method for river problem. For the actual flow problem, it is needed to compute flow velocity, water elevation and water region at the same time. For the basic equation, unsteady shallow water equations are used. As the numerical approach, implicit FEM is proposed by bubble function. To control numerical stability and accuracy, stabilized bubble function method is used to solve finite element equations. A numerical test are shown to inspect of this method.