Accelerating a barotropic ocean model using a GPU

The two-dimensional barotropic vorticity equation is one of the basic equations of ocean dynamics. It is important to have efficient numerical solution techniques to solve this equation. In this paper, we present an implementation of a numerical solution using a Graphics Processing Unit (GPU). The speedup of the calculation on the GPU with respect to that on a CPU depends on the grid size, but reaches a factor of about 50 for resolutions from 2049 × 2049 up to 4097 × 4097. It may therefore be efficient to use green GPUs in future high-resolution ocean modelling studies.

[1]  Inanc Senocak,et al.  CUDA Implementation of a Navier-Stokes Solver on Multi-GPU Desktop Platforms for Incompressible Flows , 2009 .

[2]  J. Pedlosky.  Geophysical Fluid Dynamics , 1979 .

[3]  Markus Püschel,et al.  Algebraic Signal Processing Theory: Cooley–Tukey Type Algorithms for Real DFTs , 2008, IEEE Transactions on Signal Processing.

[4]  Dinesh Manocha,et al.  LU-GPU: Efficient Algorithms for Solving Dense Linear Systems on Graphics Hardware , 2005, ACM/IEEE SC 2005 Conference (SC'05).

[5]  Zesheng Tang,et al.  Efficient magnetohydrodynamic simulations on graphics processing units with CUDA , 2009, Comput. Phys. Commun..

[6]  Manish Vachharajani,et al.  GPU acceleration of numerical weather prediction , 2008, 2008 IEEE International Symposium on Parallel and Distributed Processing.

[7]  James C. McWilliams,et al.  Large-Scale, Low-Frequency Variability in Wind-Driven Ocean Gyres , 1999 .

[8]  William H. Press,et al.  The Art of Scientific Computing Second Edition , 1998 .

[9]  Henk A. Dijkstra,et al.  Temporal variability of the wind-driven quasi-geostrophic double gyre ocean circulation: Basic bifurcation diagrams , 1997 .

[10]  L. Perelman,et al.  A finite-volume, incompressible Navier Stokes model for studies of the ocean on parallel computers , 1997 .

[11]  David Patterson,et al.  The Top 10 Innovations in the New NVIDIA Fermi Architecture, and the Top 3 Next Challenges , 2009 .

[12]  Peng Wang,et al.  Adaptive mesh fluid simulations on GPU , 2009, 0910.5547.

[13]  J. Butcher.  Numerical methods for ordinary differential equations , 2003 .

[14]  Steven G. Johnson,et al.  FFTW: an adaptive software architecture for the FFT , 1998, Proceedings of the 1998 IEEE International Conference on Acoustics, Speech and Signal Processing, ICASSP '98 (Cat. No.98CH36181).