基于Boussinesq方程的波浪数学模型主要用来模拟水波在近岸浅水区的非线性变形,被广泛应用于海岸工程中。近几十年来,其数学模型本身和数值求解方法都得到了很大的发展。但是,这也使得求解此模型的计算量增大。单个处理器由于计算速度、内存以及存储空间的限制,已经不能满足工程实际计算的需要。本文基于这一实际情况,提出Boussinesq方程的一个区域分解并行算法。本文首先回顾了Boussinesq方程在色散性、非线性和浅化性方面的改进,以及并行算法在CFD领域的发展及应用,并对并行算法的分类、并行程序性能评价方法及MPI作了介绍。在提出并行算法之前,本文先介绍了所选用的Boussinesq方程模型、方程的差分格式、数值造波方法、边界处理方法以及串行算法求解过程。本文的并行算法采用了区域分解法,x方向计算时采用有重叠的区域划分,y方向计算时采用无重叠的区域划分,并通过Schwarz迭代法实现子区域之间的耦合。由于两个方向计算时都沿同一方向划分,避免了计算时不断进行方向的转换。对三个典型物理模型试验的模拟结果表明并行计算结果与串行计算结果吻合得很好,和试验结果一致,由此验证并行算法的正确性。本文还分析了该并行算法的加速比和效率,并讨论了重叠区域大小对计算时间的影响。