Implementation of method of lines to predict water levels due to a storm along the coastal region of Bangladesh

Abstract

In this study, the method of lines (MOL) has been applied to solve two-dimensional vertically integrated shallow water equations in Cartesian coordinates for the prediction of water levels due to a storm surge along the coast of Bangladesh. In doing so, the partial derivatives with respect to the space variables were discretized by the finite difference (central) method to obtain a system of ordinary differential equations (ODEs) with time as independent variable. The classical fourth-order Runge–Kutta method was used to solve the obtained system of the ODEs. We used a nested finite difference scheme, where a high resolution fine grid model (FGM) capable of incorporating all major islands along the coastal region of Bangladesh was nested into a coarse grid model (CGM) covering up to 15°N latitude of the Bay of Bengal. The boundaries of the coast and islands were approximated through proper stair step. Appropriate tidal condition over the model domain was generated by forcing the sea level to be oscillatory with the constituent M 2 along the southern open boundary of the CGM omitting wind stress. Along the northeast corner of the FGM, the Meghna River discharge was taken into account. The developed model was applied to estimate water levels along the coast of Bangladesh due to the interaction of tide and surge associated with the April 1991 storm. We also computed our results employing the standard finite difference method (FDM). Simulated results show the MOL performs well in comparison with the FDM with regard to CPU time and stability, and ensures conformity with observations.