Water hammer simulation by explicit central finite difference methods in staggered grids

Four explicit finite difference schemes, including Lax-Friedrichs, Nessyahu-Tadmor, Lax-Wendroff and Lax-Wendroff with a nonlinear filter are applied to solve water hammer equations. The schemes solve the equations in a reservoir-pipe-valve with an instantaneous and gradual closure of the valve boundary. The computational results are compared with those of the method of characteristics (MOC), and with the results of Godunov's scheme to verify the proposed numerical solution. The computations reveal that the proposed Lax-Friedrichs and Nessyahu-Tadmor schemes can predict the discontinuities in fluid pressure with an acceptable order of accuracy in cases of instantaneous and gradual closure. However, Lax-Wendroff and Lax-Wendroff with nonlinear filter schemes fail to predict the pressure discontinuities in instantaneous closure. The independency of time and space steps in these schemes are allowed to set different spatial grid size with a unique time step, thus increasing the efficiency with respect to the conventional MOC. In these schemes, no Riemann problems are solved; hence field-by-field decompositions are avoided. As provided in the results, this leads to reduced run times compared to the Godunov scheme.