基于多GPU并行框架的河流泥沙输移过程数值模拟方法技术

本发明专利技术的基于多GPU并行框架的河流泥沙输移过程数值模拟方法，具体为：将计算区域划分为与GPU数量相同的三角形网格计算区域；对变量进行初始化；创建2个CUDA流控制2个GPU并行计算，将每个河道的数据存储到对应的GPU中；直接进行GPU和GPU间数据通信；在每个GPU中，计算每个三角形网格的源项和通量得到水深值，再由水深计算每个网格的水位值；根据得到的水位值，将每个GPU中计算的水位结果数据存储到CPU中完成水位数据合并；计算得到每个时刻河流的每个河段和出口断面的水位值。本发明专利技术采用基于三角形网格的多GPU并行计算的河流泥沙输移过程模拟方法，能高效模拟河流区域泥沙输移过程，且能快速预测断面的水位值。且能快速预测断面的水位值。且能快速预测断面的水位值。

【技术实现步骤摘要】
基于多GPU并行框架的河流泥沙输移过程数值模拟方法


[0001]本专利技术属于河流泥沙输移过程数值模拟
，涉及一种基于多GPU并行框架的河流泥沙输移过程数值模拟方法。

技术介绍

[0002]对于河流泥沙输移过程，一方面，河流泥沙输移过程涉及到悬移质和床沙输移过程，运动过程包括了泥沙的冲刷和淤积变化过程；另一方面，河流泥沙输移运动会改变河道和河床的地形，而河道和河床的复杂地形变化会引起河流某些河段水流发生强烈的起伏震荡动演变运动，需要重点研究这些区域的泥沙输移过程。因此，通过研究河流泥沙输移过程，可对河流的泥沙输移运动规律和洪水风险进行准确评估。特别是在黄河流域，泥沙的冲刷和淤积变化过程会改变河床的地形情况，继而影响黄河下游河流水位的变化，可能会导致河流下游洪水灾害的高风险性；同时可能会威胁到整个下游区域市民的生活居住和房屋财产安全，也会破坏河流下游城市水体生态环境。综上，对河流泥沙输移过程进行研究很重要，不仅能对河流泥沙输移过程的机理进行详细的分析，也能对河流洪涝规划设计和洪水应对管理预案提供可靠的依据和科学的指导。
[0003]对于河流泥沙输移过程研究的方法，国内外学者主要是采用数学模型对河流泥沙输移过程进行模拟仿真，而以水动力模型为基础构建河流泥沙输移过程数值模型已被众多学者所应用，能准确实现了泥沙输移过程的动态模拟，通过建立水沙模型分析并考虑河流泥沙运动的洪水演进变化过程，探讨河流泥沙输移运动的时空分布情况。在模拟功能方面，现有的数值模型主要是模拟单个的河流床沙泥沙输移过程或者悬移质泥沙输移过程；在应用方面，现有的数值模型较多的模拟河流泥沙输移过程，没有考虑结合河流泥沙输移过程和洪水风险的模拟评价；在计算效率方面，尽管有部分学者采用了单个GPU并行计算技术应用在河流泥沙输移过程，但是在大规模河流区域模拟计算性能较低。在计算网格方面，现有的多GPU数值计算方法，主要是基于结构网格计算的多GPU计算方法，结构网格的计算单元不能更准确反映关键精细地形的特征，因此现有方法尚未涉及到非结构网格的多GPU计算方法。综上，河流泥沙输移过程的快速模拟对泥沙动态运动过程分析、洪水风险评估、河道河床地形环境变化等很重要。因此，需要建立一种高精度、高性能、多过程的耦合河流泥沙输移过程的高性能数值模拟方法。

技术实现思路

[0004]本专利技术的目的是提供一种基于多GPU并行框架的河流泥沙输移过程数值模拟方法，解决了现有技术中存在的河流泥沙输移过程模拟方法的模拟过程功能单一和模拟计算性能不足的问题。
[0005]本专利技术所采用的技术方案是，基于多GPU并行框架的河流泥沙输移过程数值模拟方法，具体按照以下步骤实施：
[0006]步骤1，将整个河流划分为相同长度的两段河道，再将两段河道所在两部分计算区
域细分计算区域一和计算区域二，两计算区域中间相邻的两列计算网格区域为重叠区域，计算区域一和计算区域二均由多个三角形的非结构网格组成；
[0007]步骤2，将重叠区域中间边界左右两侧的三角形网格一一对应共用一条边；
[0008]步骤3，对水流及泥沙耦合过程中的变量进行初始化，并设定计算总历时、边界条件，设定每一次迭代计算的时间步长值为Δt；
[0009]步骤4，创建两个GPU，分别将计算区域一和计算区域二存储到对应的GPU中。
[0010]步骤5，创建两个CUDA流，为每个GPU分配内存，将初始化后的变量从CPU中分别复制到对应的GPU内存中；
[0011]步骤6，在多GPU中进行河流泥沙输移过程的并行计算，根据设定的时间步长Δt计算，完成第一次时间步长计算，得到每个河段下游出口断面的水位值；
[0012]步骤7，根据设定的计算总历时，在每个GPU中更新时间步长，重复步骤6，进行循环计算，直到所有时间步长的循环计算任务完成后，将每个GPU计算的结果数据存储到CPU中进行合并后输出。
[0013]本专利技术的特点还在于，
[0014]步骤2具体按以下步骤实施：将重叠区域两侧的网格单元均划分为整齐的一列网格单元，重叠区域中间边界左右两侧的三角形网格一一对应共用一条边，保证重叠区域的共用边上的界面通量能计算正确。
[0015]步骤3中，变量包括：降雨、入流流量、地形高程、曼宁系数、床沙泥沙浓度、悬移质泥沙浓度；边界条件包括：时间步长值Δt、总历时值、河道上游入口边界、下游出口边界、河道两侧的边界。
[0016]步骤6具体为：在多GPU中进行并行计算，并将重叠区域中间边界左右两侧的网格计算结果进行GPU和GPU间数据交换通信，在两个GPU上分别同步计算每个三角形网格界面的源项和通量得到水深值，再在两个GPU上同步计算水深，得到每个河段下游出口断面的水位值。
[0017]并行计算方式为：建立了水流及泥沙过程耦合输移方程，采用有限体积法对输移方程进行离散求解；
[0018]输移方程具体表示为：
[0019][0020][0021]式中：
[0022]q——为q
x
、q
y
、h的变量矢量，其中q
x
、q
y
分别为x和y方向的单宽流量，h为水深；
[0023]f、g——为x和y方向的通量矢量；
[0024]S——为源项矢量；
[0025]η——为水位值；
[0026]z
b
——为河道的地形高程；
[0027]u、v——为x和y方向的流速，其中为q
x
＝uh和q
y
＝vh；
[0028]g——为重力系数，g＝9.81m/s2；
[0029]C1、C2——为悬移质泥沙和床沙泥沙的浓度；
[0030]β——为淤积物和水流的速度差值；
[0031]S
x1
、S
x2
——为x方向悬移质泥沙和床沙泥沙的底坡源项；
[0032]S
y1
、S
y2
——为y方向悬移质泥沙和床沙泥沙的摩阻力源项；
[0033]E1、D1——为悬移质泥沙沉积率和携沙率；
[0034]E2、D2——为床沙泥沙沉积率和携沙率。
[0035]步骤7中，在2个GPU中对时间步长进行更新计算，并将重叠区域中间边界左右两侧的网格计算结果进行GPU和GPU间数据交换通信，在2个GPU上同步并行计算水流及泥沙过程耦合的输移过程，直到所有时间步长的循环计算任务完成后，将2个GPU的计算结果数据复制存储到CPU中，结果数据包括计算区域每个网格的坐标数据、地形高程数据、水位数据，再将这些数据进行数据合并，最后提取每个时刻的水位值，计算得到多个历时的每个河段下游出口断面的水位值。
[0036]本专利技术的有益效果是：本专利技术采用基于三角形网格的多GPU并行计算的河流泥沙输移过程模拟方法，能适用于具有复杂地形的河流泥沙输移过程高精度模拟，能有效进行更精细地形河流的加速模拟，在计算精度和加速性能方面具有较强本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.基于多GPU并行框架的河流泥沙输移过程数值模拟方法，其特征在于，具体按照以下步骤实施：步骤1，将整个河流划分为相同长度的两段河道，再将两段河道所在两部分计算区域细分计算区域一和计算区域二，两计算区域中间相邻的两列计算网格区域为重叠区域，所述计算区域一和所述计算区域二均由多个三角形的非结构网格组成；步骤2，将重叠区域中间边界左右两侧的三角形网格一一对应共用一条边；步骤3，对水流及泥沙耦合过程中的变量进行初始化，并设定计算总历时、边界条件，设定每一次迭代计算的时间步长值为Δt；步骤4，创建两个GPU，分别将计算区域一和计算区域二存储到对应的GPU中。步骤5，创建两个CUDA流，为每个GPU分配内存，将初始化后的变量从CPU中分别复制到对应的GPU内存中；步骤6，在多GPU中进行河流泥沙输移过程的并行计算，根据设定的时间步长Δt计算，完成第一次时间步长计算，得到每个河段下游出口断面的水位值；步骤7，根据设定的计算总历时，在每个GPU中更新时间步长，重复步骤6，进行循环计算，直到所有时间步长的循环计算任务完成后，将每个GPU计算的结果数据存储到CPU中进行合并后输出。2.根据权利要求1所述的基于多GPU并行框架的河流泥沙输移过程数值模拟方法，其特征在于，所述步骤2具体按以下步骤实施：将重叠区域两侧的网格单元均划分为整齐的一列网格单元，重叠区域中间边界左右两侧的三角形网格一一对应共用一条边，保证重叠区域的共用边上的界面通量能计算正确。3.根据权利要求1所述的基于多GPU并行框架的河流泥沙输移过程数值模拟方法，其特征在于，步骤3中，所述变量包括：降雨、入流流量、地形高程、曼宁系数、床沙泥沙浓度、悬移质泥沙浓度；所述边界条件包括：时间步长值Δt、总历时值、河道上游入口边界、下游出口边界、河道两侧的边界。4.根据权利要求1所述的基于多GPU并行框架的河流泥沙输移过程数值模拟方法，其特征在于，步骤6具体为：在多GPU中进行并行计算，并将重叠区域中间边界左右两侧的网格计算结果进行GPU和GPU间数据交换通信，在两个GPU...

【专利技术属性】
技术研发人员：韩浩，金钊，侯精明，蔺铭益，
申请(专利权)人：西安理工大学，
类型：发明
国别省市：