In an effort to discover why the dual energy formalism is producing such unsatisfactory results on the high mach number shock tube, now I will look at a simpler test. In this post, I test an advected square wave (a density enhancement in pressure equilibrium moving across a periodic grid).
In the following tests, the pressure is everywhere, , and density is , where in the middle section of the grid, and elsewhere. ranges from 0.105 in the high density regions to 0.129 in the high density regions.
The domain goes from . For the low mach number test, , so the mach number is . For the high mach number test, , so the mach number is . The tests are run until , such that the low velocity test crosses the grid once, and the high velocity test crosses the grid 100 times.
Without Dual Energy
Figure 1: PPMP, exact solver, no dual energy,
Figure 2: Same as Figure 1 but with
PPMP fares quite well on this test, even at high mach number (not surprising, given that’s what it was designed to do). But what happens if we turn on the Enzo dual energy switch?
With Dual Energy
At low mach number, the switch doesn’t activate dual energy, and the results are identical to Figure 1. At high mach number , the switch kicks in, and the results are a mess.
Figure 3: Same as Figure 2 but with dual energy
Figure 4: Same as Figure 4 but with PPMC
The Internal Energy Update
The equation that is used to update the internal energy has two parts: a term that represents the advected flux of internal energy, and a term that represents the gain or loss of internal energy due to pressure and velocity gradients. In this case, the velocity is equal across all cells, so the second term should be 0. For the first time step, this criterion holds, but by the second, there are already a few cells which show a very small velocity gradient.
Is this because the first term isn’t acting properly? Or because it’s simply not accurate enough?