Abstract
In this work, two numerical schemes were developed to overcome the problem of shock waves that appear in the solutions of one/two-layer shallow water models. The proposed numerical schemes were based on the method of lines and artificial viscosity concept. The robustness and efficiency of the proposed schemes are validated on many applications such as dam-break problem and the problem of interface propagation of two-layer shallow water model. The von Neumann stability of proposed schemes is studied and hence, the sufficient condition for stability is deduced. The results were presented graphically. The verification of the obtained results is achieved by comparing them with exact solutions or another numerical solutions founded in literature. The results are satisfactory and in much have a close agreement with existing results.