A Boundary-Implicit Constraint Reconstruction Method for Solving the Shallow Water Equations

To improve the accuracy of second-order cell-centered finite volume method in near-boundary regions for solving the two-dimensional shallow water equations, a numerical scheme with globally second-order accuracy was proposed. Having the primary objective to overcome the challenge of accuracy degradation in near-boundary regions and to develop a robust numerical framework combining high-order accuracy with strict conservation, the key research objectives had been as follows: Firstly, a physical variable reconstruction method combining a vertex-based nonlinear weighted reconstruction scheme and a monotonic upwind total variation diminishing scheme for conservation laws was proposed. While the overall computational efficiency was maintained, linear-exact reconstruction in near-boundary regions was achieved. The variable reconstruction in interior regions was integrated to achieve global second-order accuracy. Subsequently, a flux boundary condition treatment method based on uniform flow was proposed. Conservative allocation of hydraulic parameters was achieved, and flow stability in inflow regions was enhanced. Finally, a series of numerical test cases were provided to validate the performance of the proposed method in solving the shallow water equations in terms of high-order accuracy, exact conservation properties, and shock-capturing capabilities. The superiority of the method was further demonstrated under high-speed flow conditions. The high-precision numerical model developed in this study holds significant value for enhancing the predictive capability of simulations for natural disasters such as flood propagation and tsunami warning. Its robust boundary treatment methods also provide a reliable tool for simulating free-surface flows in complex environments, offering broad prospects for engineering applications.