A Box with Reflecting Walls

In this example we consider a 2-D Riemann problem inside a square box (2x2 units) with reflecting walls. Initially the density is unity everywhere and the velocities zero. The pressure is equal 1.000 inside a square (0.5x0.5 units) in the center of the box and equal 10 elsewhere.

The solution is computed using front tracking and dimensional splitting on a 400x400 grid. To reach time t=0.3 we use 1200 time steps.

Click on the link(s) below to see snapshots of temperature, density, kinetic energy, and pressure. The images are quite spectacular!

The simulation gets more polluted by numerical diffusion with time. Therefore we will not claim anything about validity of the phenomena depicted.

START the slide show.