工程科学与技术2025，Vol.57Issue(2)：104-114,11.DOI:10.12454/j.jsuese.202300374
一种改进的一维河网水动力解法及其应用
An Improved Hydrodynamic Solution for One-dimensional River Networks and its Applications
摘要
Abstract
Objectives A refined approach is introduced to effectively simulate hydrodynamics in one-dimensional(1D)river networks,enhancing the com-putational efficiency of the traditional junction water stage iteration method. Methods This approach utilizes the lattice Boltzmann Bhatnagar-Gross-Krook(BGK)model,specifically designed for 1D unsteady open-channel flows.The 1D Saint-Venant equations are solved using the lattice Boltzmann BGK model,where discrete velocity vectors are represented by the D1Q3 model.Equilibrium distribution functions are established to recover the 1D Saint-Venant equations through Chapman-Enskog analysis.The method approximates the spatial derivative of the cross-sectional area in the equilibrium distribution functions using a second-order central differ-ence scheme,achieving an overall second-order spatial accuracy.The coupling of the river network is facilitated through the junction water stage iteration method,initially segmenting the network into various sub-reaches.Water level boundary conditions are imposed at the upper and lower boundaries of each sub-reach.At each time step,the coupling water level is iteratively adjusted to maintain connectivity conditions at all junc-tions.Given that the lattice Boltzmann BGK model inherently employs an explicit algorithm,it eliminates the necessity for multiple iterations to compute the coupling water level. Results Standard verification scenarios,including steady flows in a looped channel network and unsteady flows in combined looped and treelike networks,were simulated to evaluate the accuracy of the refined method.The hydrodynamics of the middle and lower reaches of the Han River during the autumn flood of 2021 were simulated to validate its utility in real-world applications.The results are as follows:1)For the steady flow scenario in a looped channel network,the simulation outcomes were benchmarked against theoretical analyses,achieving a relative error of up to 7.43%.2)In scenarios involving unsteady flows in combined looped and treelike networks,the refined method precisely simulated the flooding process,with the flood hydrograph's shape,peak discharge and flood rise and fall times aligning well with prior studies.Various metrics,includ-ing total simulation time,total iteration count,average iterations per time step,and average iteration time of the refined method,were compared to those obtained using the traditional junction water stage iteration method and the junctions group method.This comparison revealed that the total and per-step iteration counts of the refined method match those of the traditional method.However,the refined method's average iteration time was significantly reduced due to eliminating multiple global evolutions per iteration.3)the simulated results were compared to measured data and calculations performed using MIKE software in a real-world case.A remarkable agreement was achieved,with the maximum global relative error of 10-2-10-3.A comparison of inflow and outflow discharges at a specific river network junction,obtained using the proposed method,the meth-od of characteristics,the Lax diffusive scheme,and the original junction water stage iteration method,revealed that the internal boundary condi-tion of river network junctions cannot be strictly met when employing either the method of characteristics or the Lax diffusive scheme.However,inflow and outflow discharges at river network junctions can be accurately balanced using the proposed method and its original counterpart. Conclusions These findings highlighted the importance,within the context of the junction water stage iteration method,of ensuring that the meth-od utilized to solve the 1D Saint-Venant equations was consistent with that employed to iteratively solve the coupling water level at each river net-work junction.Although traditional explicit methods,such as the method of characteristics and the Lax diffusive scheme,enhanced computation-al efficiency by circumventing the global evolution within each iteration,computational accuracy cannot be guaranteed unless a uniform numeric-al method was employed for the global evolution.The proposed method improved computational efficiency without compromising computational accuracy and ensured that the connection conditions at each river network junction were rigorously maintained.关键词
计算水力学/一维河网非恒定流/格子Boltzmann方法/汊点水位迭代法/汉江流域Key words
computational hydraulics/unsteady flow in one-dimensional river networks/Lattice Boltzmann method/junction water stage itera-tion method/Han River分类
水利科学引用本文复制引用
吴家阳,徐学军,程永光..一种改进的一维河网水动力解法及其应用[J].工程科学与技术,2025,57(2):104-114,11.基金项目
国家重点研发计划项目(2022YFC3002703 （）
2021YFC3000102) （）
