An implicit method with selective total variation diminishing (TVD) term inclusion was developed and used to find solutions to the shallow water equations (SWEs) for free-surface flows in natural and engineered channels. The conservative form of the SWEs was employed, and the method incorporated an algorithm for selectively underrelaxing the iterative process to maintain stability and accuracy in the presence of shock interfaces. The value of the Courant number and the frequency at which the TVD term was incorporated were constantly updated during the computation to achieve optimal speed of execution while maintaining stability. The method was tested against published results from physical experiments and from computations employing alternative algorithms, and the results obtained demonstrate both the economy and accuracy of the proposed algorithm.

In this paper, we present a GPU implementation of a two-dimensional shallow water model. Water simulations are useful for modeling floods, river/reservoir behavior, and dam break scenarios. Our GPU implementation shows vast performance improvements over the original Fortran implementation. By taking advantage of the GPU, researchers and engineers will be able to study water systems more efficiently and in greater detail.