[发明专利]海啸数值计算模型基于GPU并行的快速执行方法

权利要求书

1.一种海啸数值计算模型基于GPU并行的快速执行方法，其特征在于，包括：

CPU接收海啸数值计算模型的启动指令后，获取与所述海啸数值计算模型相关的参数信息，并获取初始海表面形变信息；

所述CPU读取到数据传输指令时，所述CPU将所述参数信息、初始海表面形变信息通过PCI-E3.0接口一次性传输至GPU的显存；

所述GPU接收所述参数信息、初始海表面形变信息后，执行所述海啸数值计算模型的海啸从深海传播至近岸计算过程，其中，所述海啸数值计算模型采用有限差分方法求解球坐标系下二维线性方程和非线性方程，获取海啸的计算结果；所述非线性方程考虑底摩擦效应；

所述GPU读取到数据回传指令时，所述GPU将所述海啸的计算结果传输至所述CPU；

所述CPU根据输出指令，将所述海啸的计算结果进行格式化处理并输出；

其中，用于海啸数值计算的GPU作为CPU的协处理器对海啸传播计算模块进行加速；所述GPU接收所述参数信息、初始海表面形变信息后，执行所述海啸数值计算模型的海啸传播计算过程，获取海啸的计算结果的步骤，包括：

所述CPU在确定需要获取至少一个时间节点的海啸计算结果时，所述CPU向所述GPU发送至少一个时间节点的海啸中间计算结果获取指令；

所述GPU根据所述至少一个时间节点的海啸中间计算结果获取指令，将计算过程中符合所述时间节点的海啸中间计算结果传输至所述CPU；

具体地，第一、根据地震发生位置和潜在海啸影响区域确定海啸数值计算模型计算区域，设定海啸计算空间范围、空间分辨率Δx、时间分辨率Δt以及计算时长；根据空间范围和空间分辨率采用公式(10)计算生成球坐标系下正交曲线网格以及相应的网格参数，即读取原始地形水深文件，插值得到地形水深网格数据；计算Courant系数；

其中，g表示重力加速度，h_max为计算范围内的最大水深；

即利用Fortran语言编程，在CPU端声明所有变量并赋初值，编写子程序读入一个固定格式的模型参数文件，根据计算区域和空间分辨率计算正交曲面网格每一点的经纬度并插值水深；编写子程序对Courant系数是否满足稳定性条件进行判断，若不满足，对时间步长进行调整；

第二、获取海底断层破裂引起的海底形变信息，即初始海表面形变量；即利用Fortran语言编程，根据Okada断层模型，结合读入的地震和断层参数计算初始海表面形变；或者直接读入文件格式的初始海表面形变量和动量通量信息；

第三、将海啸传播计算模块需要的变量和参数通过PCI-3.0接口传递至GPU显存；即利用CUDA_C语言的cudaMalloc命令，在GPU端声明变量、分配显存，然后再用cudaMemcpy(*，cudaMemcpyHostToDevice)命令，实现从CPU向GPU传递所有参与海啸传播计算的变量和参数，包括海平面初始垂向位移、沿经度和纬度动量通量；

第四、海啸波在大洋深水区的传播计算，即数值求解线性浅水方程，包括质量连续方程、经向和纬向动量方程和定解边界条件；

即利用CUDA_C语言编写在设备端执行的内核函数__global__void**_kernel(…)，包括离散化的质量连续方程、经向和纬向动量方程和海啸波吸收边界算法；采用交错式显性蛙越方案求解方程，网格中心点是海啸波高和水深，上下及相邻四边为沿经向和纬向的动量通量；

时间上，GPU执行逐时间步长迭代计算；空间上，GPU执行沿经度和纬度两个方向求解；当执行至网格边界时，采用吸收边界算法滤波；同一时间步长，GPU必须遍历所有网格点，然后才能进入下一个时间步长计算，循环往复直至设定的计算时长；方程的解包括海啸波高和两个方向的海啸动量通量；

第五、海啸波在浅水区域的传播计算，即数值求解非线性浅水方程，包括质量连续方程、经向和纬向动量方程和吸收边界算发；动量方程中增加了底摩擦项和非线性对流项；

第六、输出保存海啸计算结果，采用NetCDF格式对结果进行输出保存。

2.根据权利要求1所述的方法，其特征在于，还包括：

所述CPU将所有的海啸中间计算结果进行格式化处理并输出。