摘要
针对具有干河床、急缓流态交替的河网水动力模拟,采用MUSCL-Hancock有限体积格式离散Saint-Venant方程组,运用汊点水位预测-校正法处理河网汊点连接条件,建立了断面与河道完全数值解耦的复杂河网水动力模型。运用OpenMP和OpenACC编程模式分别实现了CPU,GPU并行计算。采用经典算例对模型计算精度进行了验证,采用珠三角河网实际算例对不同并行计算方式的加速效果进行了分析。算例研究表明,该模型具有良好的稳定性和计算效率,可适用于陡坡河道水动力模拟,具有较好的推广应用价值。
A hydrodynamic model of complex river network in which the computational cross-sections and reaches are fully decoupled is constructed for river network with dry beds and transition between supercritical and subcritical flows. The MUSCL-Hancock finite-volume scheme is adopted to discrete the Saint-Venant equations, and the method of Junction-Point Water Stage Prediction and Correction(JPWSPC) is used to calculate the hydrodynamic relations at junctions. The OpenMP and OpenACC applications are used to realize CPU-based and GPU-based parallel computing, respectively. Model accuracy is validated by benchmark cases with exact solutions, and the performance of parallel computing is demonstrated by real flood simulation of the Pearl River Delta. Results show that the proposed model is stable and effective, and is applicable for hydrodynamic simulation for steep river, thus has bright application prospects.
引文
[1] 王宗志,程亮,王银堂,等.高强度人类活动作用下考虑河道下渗的河网洪水模拟[J].水利学报,2015,46(4):414-424.
[2] 王建华,肖伟华,王浩,等.变化环境下河流水量水质联合模拟与评价[J].科学通报,2013,58(12):1101-1108.
[3] 李晓,唐洪武,王玲玲,等.平原河网地区闸泵群联合调度水环境模拟[J].河海大学学报(自然科学版),2016,44(5):393-399.
[4] 方神光,张文明,徐峰俊,等.西江中游河网及梯级水库水动力整体数学模型研究[J].人民珠江,2015,(4):107-111.
[5] PREISSMANN A.Propagation des Intumescences dans les Canaux et les Rivieres[C]//Proceedings of the 1961 First Congress of the French Association for Computation,AFCAL Grenoble,France,1961:433-442.
[6] 张二俊,张东生,李挺.河网非恒定流三级联合算法[J].华东水利学院学报,1982,(1):1-13.
[7] 吴寿红.河网非恒定流四级解法[J].水利学报,1985,(8):42-50.
[8] 陈炼钢,施勇,钱新,等.闸控河网水文-水动力-水质耦合数学模型——I.理论 [J].水科学进展,2014,25(4):534-541.
[9] 吴晓玲,向小华,牛帅,等.急缓流态交替的陡坡河道水动力模型[J].水科学进展,2017,28(4):564-570.
[10] 朱德军,陈永灿,王智勇,等.复杂河网水动力数值模型[J].水科学进展,2011,22(2):203-207.
[11] 张大伟,权锦,马建明,等.应用Godunov 格式模拟复杂河网明渠水流运动[J].应用基础与工程科学学报,2015,23(6):1088-1096.
[12] 向小华,吴晓玲,牛帅,等.通量差分裂方法在一维河网模型中的应用[J].水科学进展,2013,24(6):894-900.
[13] 王船海,曾贤敏.Windows 环境下河网水流多线程并行计算[J].河海大学学报(自然科学版),2008,36(1):30-34.
[14] 刘荣华,魏加华,朱德军,等.基于JPWSPC 的河网一维水动力学模拟并行计算研究[J].应用基础与工程科学学报,2015,23(2):213-224.
[15] 张国义,房明惠,黄刘生,等.树状河网非恒定流的一则并行直接数值解法[J].系统仿真学报,2004,16(12):2673-2676,2679.
[16] 张大伟,程晓陶,黄金池,等.复杂明渠水流运动的高适用性数学模型[J].水利学报,2010,41(5):531-536.
[17] YING X Y,WANG S S Y.Improved Implementation of the HLL Approximate Riemann Solver for One-dimensional Open Channel Flows [J].Journal of Hydraulic Research,2008,46(1):21-34.
[18] SANDERSB F.High-resolution and Non-oscillatory Solution of the St.Venant Equations in Non-rectangular and Non-prismatic Channels [J].Journal of Hydraulic Research,2001,39(3):321-330.
[19] 宋利祥,周建中,王光谦,等.溃坝水流数值计算的非结构有限体积模型.水科学进展,2011,22(3):373-381.
[20] SONG L,ZHOU J,GUO J,et al.A Robust Well-balanced Finite Volume Model for Shallow Water Flows with Wetting and Drying over Irregular Terrain [J].Advances in Water Resources,2011,34(7):915-932.
[21] 何沧平.OpenACC并行编程实战[M].北京:机械工业出版社,2017.
[22] ISLAM A,RAGHUWANSHI N S,SINGH R,et al.Comparison of Gradually Varied Flow Computation Algorithms for Open-channel Network [J].ASCE Journal of Irrigation and Drainage Engineering,2005,131(5):457-465.
