GPU accelerated CESE method for 1D shock tube problems

Wei Ran, Wan Cheng, Fenghua Qin, Xisheng Luo*

*Corresponding author for this work

Research output: Contribution to journalArticlepeer-review

11 Scopus citations


In the present study, a GPU accelerated 1D space-time CESE method is developed and applied to shock tube problems with and without condensation. We have demonstrated how to implement the CESE algorithm to solve 1D shock tube problems using an older generation GPU (the NVIDIA 9800 GT) with relatively limited memory. To optimize the code performance, we used Shared Memory and solved the inter-Block boundary problem in two ways, namely the branch scheme and the overlapping scheme. The implementations of these schemes are discussed in detail and their performances are compared for the Sod shock tube problems. For the Sod problem without condensation, the speedup over an Intel CPU E7300 is 23 for the branch scheme and 41 for the overlapping scheme, respectively. While for problems with condensation, both schemes achieve higher acceleration ratios, 53 and 71, respectively. The higher speedup of the condensation case can be ascribed to the source term calculation which has a local dependence on the mesh point and the SOURCE kernel has a higher acceleration ratio.

Original languageEnglish (US)
Pages (from-to)8797-8812
Number of pages16
JournalJournal of Computational Physics
Issue number24
StatePublished - Jan 1 2011


  • CUDA GPU computing
  • Condensation
  • Shock tube
  • The space-time CESE method

ASJC Scopus subject areas

  • Numerical Analysis
  • Modeling and Simulation
  • Physics and Astronomy (miscellaneous)
  • Physics and Astronomy(all)
  • Computer Science Applications
  • Computational Mathematics
  • Applied Mathematics


Dive into the research topics of 'GPU accelerated CESE method for 1D shock tube problems'. Together they form a unique fingerprint.

Cite this