|
Cococubed.com
|
| Exact Riemann Solver |
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
... EZ stellar evolution
... FLASH code
... Mesa code
Astronomy talks
Astronomy images
Outreach
Family album
Bicycle adventures
Artwork
Contact us:
J.D. Maldonado
F.X.Timmes, my vitae
The code exact_riemann.f solves the riemann problem of discontinuity across an interface.
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 analytical solution for the at 0.4 s is shown below.
Density as a function of position |
Pressure profile |
Fluid speeds |
| exact_riemann.f |
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.