*
Cococubed.com


Chapman-Jouget Detonations

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
Given
  • the thermodynamics of the fuel (the fuel's temperature, density and composition)
  • the ashes exist in their equilibrium state (NSE in the nuclear case)
then the Chapman-Jouget (1890) detonation solution follows from conservation of momentum, energy, mass, and charge across the flame front. This gives four nonlinear, algebraic conditions whose simultaneous solution is desired. That is, a four dimensional root find is done.

The CJ solution tells you the
  • speed of the detonation
  • thermodynamics of the ashes
     
The CJ solution doesn't tell you the
  • width of the fuel-to-ash region
  • spatial variations of the variables
  • if the solution is a self-sustaining detonation
If you want what the Chapman-Jouget solution doesn't tell you, a ZND detonation will. Besides doing CJ detonations, the program below will also compute the strong and weak detonation solutions in the event that one chooses to drive the detonation at a specified mach number.


image
detonation speed 		image
why speed is non-monotonic 		image
detonation regimes


cjdet.f
plus these include files:
implno.dek
const.dek
network.dek
vector_eos.dek
cjdet.dek
plus a 13 MB eos data file
helm_table.dat plus a 6 MB nuclear reaction data file from Rob Hoffman:
BDAT


 



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.