*
Cococubed.com


Piecewise Parabolic Method (PPM)

Home

Teaching materials
Astronomy research
Astronomy codes
... Stellar equation of states
... EOS with ionization
... EOS for supernovae
... Chemical potentials
... Stellar atmospheres
... Voigt Function
... Polytropic stars
... Cold white dwarfs
... Hotter white dwarfs
... Cold neutron stars
... Stellar opacities
... Neutrino energy loss rates
... Ephemeris routines
... Fermi-Dirac functions
... Galactic chemical evolution

... Nuclear reaction networks
... Nuclear statistical equilibrium
... Laminar deflagrations
... CJ detonations
... ZND detonations
... Fitting to conic sections
... Unusual linear algebra
... Derivatives on uneven grids
... Pentadiagonal solver
... Quadratics, Cubics, Quartics
... Supernova light curves
... Exact Riemann solutions
... 1D PPM Hydrodynamics
... Verification problems
... Plane - Cube Intersection

... FLASHcode
... MESA code
... GR1D code
Astronomy talks
Astronomy images
Outreach
Bicycle adventures
Artwork

Contact:
F.X.Timmes, my vitae
The code ppm.f, and its include file ppm.dek, implements a 1D PPM hydrodynamics solver.

While the Sod problem has become a standard test case for hydrodynamics codes it isn't a very discriminating test for modern codes. The strong shock tube problem below is a more demanding test problem, because of the stronger discontinuities across the shock interface and the narrow density peak that forms behind the shock. The code sets the initial density to 10 on the left, and 1 on the right. The initial pressure is set to 100 on the left, and 1 on the right. The PPM solution and the exact solution to the exact Riemann solution at 0.4 s is compared below.


image
Pressure as a function of position 		image
Density profile 		image
Fluid speeds


ppm_1d.f and its include file ppm.dek


 



If you want a modern multi-physics, multi-scale, multi-dimensional hydrodynamics code, check out The FLASH code

Please cite the relevant references if you publish a piece of work that use these codes, pieces of these codes, or modified versions of them. If you're nice, offer co-authorship of the publication. At best, you'll love these programs so much that you'll send great wads of cash to me.