To address the severe urban flooding and the inability of urban road drainage systems to effectively resolve hydrological cycle issues, four types of permeable pavement were designed, combining the advantages of the good infiltration performance and anti-slip performance of permeable asphalt pavement. Based on the SWMM (Storm Water Management Model), road modeling and hydrological quality simulations were conducted, analyzing the surface runoff reduction rate, maximum inlet flow at the convergence node, pollutant reduction rate, and water quality purification of the four different structural permeable asphalt pavements. The results showed that the surface runoff reduction rate of the four pavements ranged from 44% to 100%, the maximum inlet flow reduction rate at the confluence node ranged from 37% to 78%, the reduction rate of the main pipe flow load ranged from 36% to 100%, the reduction rate of the hydraulic load in the pipeline ranged from 25% to 64%, the maximum water storage depth ranged from 90 mm to 177 mm, and the pollutant reduction rate ranged from 28% to 81%. This study provides the optimal combination of permeable asphalt pavements for the selection of urban pavement structures.

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.