|
Cococubed.com
|
| Rarefaction wave (Mader) |
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 simplest test of detonation is the one-dimensional gamma-law rarefaction wave burn, for which a slab of material is initiated on one side and a detonation propagates to the other side. For a Chapman-Jouget detonation speed of 0.8 cm/ s, it takes 6.25 us for the detonation to travel 5 cm. The rich structure of a multi-dimensional detonation is absent in the one-dimensional test problem, and a simple rarefaction wave follows the detonation front (Fickett & Davis 1979). Expansion of material in the rarefaction depends on the boundary condition where the detonation is initiated, which is usually modeled as a freely moving surface or a moving piston. For the Mader problem, a stationary piston is used. In this case, the head of the rarefaction remains at the detonation front since the flow is sonic there, and the tail of the rarefaction is halfway between the front and the piston. This code released under LA-CC-05-101.
 simple-minded detonation here is a less-simple one |
 analytical and numerical |
 convergence study |
 evolution movie |
rarefaction.f solves the Mader problem |
 evolution movie |
Density uniform mesh:
 density 200x200 cells |
 density 400x400 cells |
 density 800x800 cells |
 density error 200x200 cells |
 density error 400x400 cells |
 density error 800x800 cells |
 density asymmetry 200x200 |
 density asymmetry 400x400 |
 density asymmetry 800x800 |
Pressure uniform mesh:
 pressure 200x200 cells |
 pressure 400x400 cells |
 pressure 800x800 cells |
 pressure error 200x200 cells |
 pressure error 400x400 cells |
 pressure error 800x800 cells |
 pressure asymmetry 200x200 |
 pressure asymmetry 400x400 |
 pressure asymmetry 800x800 |
Material Speed uniform mesh:
 speed 200x200 cells |
 speed 400x400 cells |
 speed 800x800 cells |
 speed error 200x200 cells |
 speed error 400x400 cells |
 speed error 800x800 cells |
 speed asymmetry 200x200 |
 speed asymmetry 400x400 |
 speed asymmetry 800x800 |
Density adaptive mesh:
 density 200x200 cells |
 density 400x400 cells |
 density 800x800 cells |
 density error 200x200 cells |
 density error 400x400 cells |
 density error 800x800 cells |
Pressure adaptive mesh:
 pressure 200x200 cells |
 pressure 400x400 cells |
 pressure 800x800 cells |
 pressure error 200x200 cells |
 pressure error 400x400 cells |
 pressure error 800x800 cells |
Material Speed adaptive mesh:
 speed 200x200 cells |
 speed 400x400 cells |
 speed 800x800 cells |
 speed error 200x200 cells |
 speed error 400x400 cells |
 speed error 800x800 cells |