mPPM, Viewed as a Co-Design Effort

The Piecewise Parabolic Method (PPM) was designed as a means of exploring compressible gas dynam-ics problems of interest in astrophysics, including super-sonic jets, compressible turbulence, stellar convection, and turbulent mixing and burning of gases in stellar interiors. Over time, the capabilities encapsulated in PPM have co-evolved with the availability of a series of high performance computing platforms. Implementation of the algorithm has adapted to and advanced with the architectural capabilities and characteristics of these machines. This adaptability of our PPM codes has enabled targeted astrophysical applica-tions of PPM to exploit these scarce resources to explore complex physical phenomena. Here we describe the means by which this was accomplished, and set a path forward, with a new miniapp, mPPM, for continuing this process in a diverse and dynamic architecture design environment. Adaptations in mPPM for the latest high performance machines are discussed that address the important issue of limited bandwidth from locally attached main memory to the microprocessor chip.

[1]  S. Zalesak.  Introduction to “Flux-Corrected Transport. I. SHASTA, A Fluid Transport Algorithm That Works” , 1997 .

[2]  R. Hornung,et al.  HYDRODYNAMICS CHALLENGE PROBLEM , 2011 .

[3]  P. Woodward,et al.  The Piecewise Parabolic Method (PPM) for Gas Dynamical Simulations , 1984 .

[4]  Paul R. Woodward,et al.  First experience of compressible gas dynamics simulation on the Los Alamos roadrunner machine , 2009, Concurr. Comput. Pract. Exp..

[5]  B. Vanleer,et al.  Towards the ultimate conservative difference scheme. IV. A new approach to numerical convection , 1977 .

[6]  P. Woodward.  Shock-driven implosion of interstellar gas clouds and star formation. , 1976 .

[7]  I. Bohachevsky,et al.  Finite difference method for numerical computation of discontinuous solutions of the equations of fluid dynamics , 1959 .

[8]  Paul R. Woodward,et al.  CFD Builder: A Library Builder for Computational Fluid Dynamics , 2014, 2014 IEEE International Parallel & Distributed Processing Symposium Workshops.

[9]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme V. A second-order sequel to Godunov's method , 1979 .

[10]  L. Margolin,et al.  Implicit Large Eddy Simulation: Computing Turbulent Fluid Dynamics , 2011 .

[11]  Paul R. Woodward.  Numerics for iles: The ppm compressible gas dynamics scheme , 2007 .

[12]  David H. Bailey,et al.  The Nas Parallel Benchmarks , 1991, Int. J. High Perform. Comput. Appl..

[13]  P. Woodward,et al.  HYDRODYNAMIC SIMULATIONS OF H ENTRAINMENT AT THE TOP OF He-SHELL FLASH CONVECTION , 2013, 1307.3821.

[14]  Paul R. Woodward,et al.  Portable petaFLOP/s programming: applying distributed computing methodology to the grid within a single machine room , 1999, Proceedings. The Eighth International Symposium on High Performance Distributed Computing (Cat. No.99TH8469).

[15]  Phillip Colella,et al.  High resolution difference schemes for compressible gas dynamics , 1981 .

[16]  P. Woodward.  Piecewise-parabolic methods for astrophysical fluid dynamics , 1986 .

[17]  Jack J. Dongarra,et al.  The LINPACK Benchmark: past, present and future , 2003, Concurr. Comput. Pract. Exp..

[18]  P. Woodward,et al.  The numerical simulation of two-dimensional fluid flow with strong shocks , 1984 .

[19]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme. IV. A new approach to numerical convection , 1977 .

[20]  B. C. Curtis,et al.  Very High Resolution Simulation of Compressible Turbulence on the IBM-SP System , 1999, ACM/IEEE SC 1999 Conference (SC'99).

[21]  Paul R. Woodward,et al.  The Piecewise-Parabolic Boltzmann Advection Scheme (PPB) Applied to Multifluid Hydrodynamics , 2010 .