This computation was performed at 100 mesh points per ZND half-reaction length. Total number of mesh points was 1000x500. The dimensionless parameters were overdrive 1.2, activation energy 10.0, heat release 50.0, gamma 1.2.
Each frame consists of plots of (from top to bottom) pressure gradient, temperature gradient, vorticity, and gradient of fuel mass fraction. Only 650 x 500 cells of the computational domain are shown on each plot. Time is increasing from left to right and from top to bottom, with 250 time steps between each frame. Make your browser about 800 pixels wide to see all the frames in a single row.
The main transverse wave interacts with the smaller vortices along the shear layer, where a Kelvin-Helmholtz instability is evident. It finally begins to hit the large vortex at the bottom center of frame ten. The interaction is most evident on each pressure gradient frames. The position of the vortex can be more clearly seen in the temperature, vorticity and mass fraction plots. The transverse wave does seem to get deformed as it traverses across the vortex.
Click on the smaller images to see the larger ones (260 KB - gifs).
