Initial Conditions: Shock Heated Mixture

(CJ Speed)

- demo_cvCJ.m

Matlab program used to generate the sample output described below

- Numerical Solution Methods for Control Volume Explosions and ZND Detonation Structure

Documentation discussing algorithm used for this demo

Input Data | |
---|---|

Mixture | Stoichiometric H2/O2 'H2:2 O2:1 N2:3.76' |

Mechanism | h2air_highT.cti |

Initial Pressure | 1 bar = 100000 Pa |

Initial Temperature | 300 K |

Output Data | |

CJ speed | 1967.8451 m/s |

gas object at equilibrium | |

Output file |

The CJ velocity is the minimum of the plot | Curve fit to find the minimum of the data |

Constant Volume Temperature & Pressure Profiles

& Induction Time (tau) calculated as the time to the maximum temperature gradient.