A Fortran90 code that computes this exact solution using a Newton-Raphson method has been written and maybe found in the repository. The velocity spectra of the exact solution is also computed (using an FFT).
While the code may be easily modified, it currently assumes that the initial velocity profile is a single mode of the form
For such an initial condition it can be shown that the shock develops at the timeAn example simple wave an .avi animation shows the density up to breaking as computed by the exact solution, while a simulation using the TCD-WENO scheme computes past breaking. The parameters for this example were
This exact solution with the is used to study the convergance of the skew-symetric form of the solver in the Simple Wave Simulations section.
-- DavidHill - 20 Dec 2004