Abstract
This paper describes a numerical method to predict scouring process in downstream reaches of bed protection works with hollow squared blocks. Modeling the distributions of the bed shear stress measured in scour holes in a experiment, variations of scour hole profiles with time are simulated by solving numerically the continuity equation of the bed material in non-equilibrium states. In addition, calculations of water surface profiles at every stage of the scour process are carried out using the two step Lax-Wendroff scheme. The validity of the present method is conformed by comparisons of the numerical results with those of experiments under static and dynamic conditions of sediment motion.