MPS‐Based Model to Solve One‐Dimensional Shallow Water Equations

Abstract

AbstractHere, a moving particle simulation method is presented to spatially integrate the cross‐sectional average shallow water equations using a prediction‐correction procedure for the time discretization. A density‐ratio equation is derived for the water depth computation according to the particle number density concept. The newly derived equation does not miscalculate the water depth in case an incorrect searching radius parameter is adopted, unlike the typical volume‐summation formula in meshless shallow water flows. A new one‐dimensional Spiky kernel function is developed to satisfy the unity condition employed in the Newton‐Raphson iteration to calculate the water depth. Dynamic stabilization is adopted to capture shockwave problems, a case‐independent technique based on the inelastic collision with unequal masses. The convective flux term is eliminated under the Lagrangian framework, and so the momentum is adequately conserved without the need for any special treatment. The proposed scheme maintains the exact C‐property, meaning that the water depth gradient and bed slope are hydrostatically well balanced within a discretized solution domain. Compared to analytical solutions and experimental data, the results of this study reveal that the present model is a robust numerical solver without unphysical oscillations. It can capture various shock problems, including steep gradient shock‐front, discontinuous wet bed, transcritical flow regimes, and irregular bed topography. The present model can also simulate the dry‐wet flow transition influenced by friction without experiencing any divisions by zero, negative values, or unphysical perturbations. This advantage is basically due to the water particles' absence and presence for the dry and wet regions, respectively.