Improvement of Fluid Simulation Runtime of Smoothed Particle Hydrodynamics by Using Graphics Processing Unit (GPU)

Wahyu Srigutomo, Ruddy Kurnia, Suprijadi Suprijadi


This study concerns an implementation of smoothed particle hydrodynamics (SPH) fluid simulation on a graphics processing unit (GPU) using the Compute Unified Device Architecture’s (CUDA) parallel programming. A bookkeeping method for the neighbor search algorithm was incorporated to accelerate calculations. Based on sequence code profiling of the SPH method, particle interaction computation – which comprises the calculation of the continuity equation and the momentum conservation equation – consumes 95.2% of the calculation time. In this paper, an improvement of the calculation is proposed by calculating the particle interaction part on the GPU and by using a bookkeeping algorithm to restrict the calculation only to contributed particles. Three aspects are addressed in this paper: firstly, speed-up of the CUDA parallel programming computation as a function of the number of particles used in the simulation; secondly, the influence of double precision and single precision schemes on the computational acceleration; and thirdly, calculation accuracy with respect to the number of particles. Scott Russell’s wave generator was implemented for a 2D case and a 3D dam-break. The results show that the proposed method was succesfull in accelerating the SPH simulation on the GPU.


Compute Unified Device Architecture (CUDA); fluid simulation runtime; graphic processing unit (GPU); Scott Russell’s wave generator; smoothed particle hydrodynamics (SPH).

Full Text:



Kurnia, R., Omata, S. & Kazama, M., The Stochastic Smoothed Particle Hydrodynamics to Overcome Energy Loss. Gakuto International series, Mathematical sciences and applications, 34, pp. 157-174, 2011.

Monaghan, J.J. & Gingold, R.A., Shock Simulation by the Particle Method SPH, J. Comput. Phys., 52, pp. 374-389, 1983.

Goswami, P., Schlegel, P., Solenthales B. & Pajarola, R., Interactive SPH Simulation and Rendering on the GPU, ACM SIGGRAPH Symposium on Computer Animation, pp. 1-10, 2010.

Amada, T., Imura, M., Yasumoro, Y., Wanabe, Y. & Chihara, K, Particle-based Fluid Simulation on GPU, ACM Workshop on General-Purpose Computing on Graphic Processors, 2004.

Monaghan, J.J., Smoothed Particle Hydrodynamics, Annu. Rev. Astron. Astrophys., 30, pp. 543-574, 1992.

Colagrossi, A. & Landrini, M., Numerical Simulation of Interfacial Flows by Smoothed Particle Hydrodynamics, J. Comput. Phys., 191, pp. 448- 475, 2003.

Monaghan, J.J., Simulating Free Surface Flow with SPH, J. Comput. Phys., 110, pp. 399-406, 1994.

Monaghan, J.J. & Kos, A., Solitary Waves on a Cretan Beach, Journal of waterway, port, coastal, and ocean engineering, 125, pp. 145-154, 1999.

Monaghan, J.J., Particle Methods for hydrodynamics, Computer Physics Report, 3, pp. 71-124, 1985.

Ashtiani, B.A. & Rezaei, A.M., Modification of Weakly Compressible Smoothed Particle Hydrodynamics for Preservation of Angular Momentum in Simulation of Impulse Wave Problems, Coastal Engineering Journal, 51(4), pp. 363-386, 2009.

Monaghan, J.J. & Kos, A., Scott Rusell’s Wave Generator, Physics of Fluids, 12, pp. 622-630, 2000.

Lo, E.Y.M. & Shao, S., Simulation of Near-shore Solitary Wave Mechanis by an Incompressible SPH Method, Applied Ocean Research, 24, pp. 275-286, 2002.

Marrone, S., Colagrossi, A., Le Touze, D. & Graziani, G., Fast Free-surface Detection and Level Set Function Definition in SPH Solvers, J. Comput. Phys., 229, pp. 3652-3663, 2010.

Doring, M., Développement d’une Méthode SPH Pour Les Applications à Surface Libre en Hydrodynamique, PhD thesis, École Centrale Nantes, 2005. (Text in French)



  • There are currently no refbacks.

Contact Information:

ITB Journal Publisher, LPPM – ITB, 

Center for Research and Community Services (CRCS) Building Floor 7th, 
Jl. Ganesha No. 10 Bandung 40132, Indonesia,

Tel. +62-22-86010080,

Fax.: +62-22-86010051;