Performance Comparison of the Three Numerical Methods to Discretize the Local Inertial Equation for Stable Shallow Water Computation Tomohiro Tanaka 1 , Hidekazu Yoshioka 2(&) , Sokly Siev 3 , Hideto Fujii 4 , Ly Sarann 5 , and Chihiro Yoshimura 3 1 Kyoto University, Kyotodaigaku-katsura, Nishikyoku, Kyoto, Japan tanaka.tomohiro.7c@kyoto-u.ac.jp 2 Shimane University, Nishikawatsu-cho 1060, Matsue, Japan yoshih@life.shimane-u.ac.jp 3 Tokyo Institute of Technology, Ookayama, Meguro-ku, Tokyo, Japan {siev.s.aa,yoshimura.c.aa}@m.titech.ac.jp 4 Yamagata University, 1-23, Wakaba-machi, Tsuruoka, Yamagata, Japan fhideto@tds1.tr.yamagata-u.ac.jp 5 Institute of Technology of Cambodia, P.O. Box 86, Phnom Penh, Cambodia ly.sarann@itc.edu.kh Abstract. The local inertial equation (LIE) is a simple shallow water model for simulating surface water dynamics. Recently, the model has been widely applied to flood simulation worldwide. Keys in numerical implementation of the LIE are the staggered spatio-temporal discretization and the stable treatment of the friction slope terms. The latter is critical for stable and ef ficient computation. Currently, several discretization methods (semi-implicit, fully-implicit, and exponential methods) for the friction slope terms with comparable computa- tional ef ficiency are available. However, their performance evaluation has been carried out only independently. We thus compare the performance of the three methods through their application to test and realistic cases. In this paper, firstly, theoretical linear stability analysis results are reviewed, indicating the highest stability of the implicit method. It is also consistent in a certain sense. Appli- cation of these methods to a 1-D test case with an advancing wet and dry interface implies that all the methods work well where the fully-implicit method has the least error. Their application to 2-D flood simulation in Tonle Sap Lake and its floodplains in South-East Asia demonstrates that the exponential method gives slightly more oscillatory results than the others. Dependence of the sim- ulated surface water dynamics on the spatial resolution is investigated as well to give a criterion of the resolution for numerical simulation with satisfactory accuracy. Keywords: Local inertial equation  Finite difference scheme Friction slope term  Tonle Sap Lake © Springer Nature Singapore Pte Ltd. 2018 L. Li et al. (Eds.): AsiaSim 2018, CCIS 946, pp. 451–465, 2018. https://doi.org/10.1007/978-981-13-2853-4_35