186 | 186 | Now the EXACT solution at t=0.035 could be computed from the Exact Riemann Solver alone, and this would be by looping through the grid at different dx/dt = dx/0.035, to determine the relative position of this 'grid speed' to the shocks and contact waves generated at the boundary. If the grid speed is ahead of all shocks (far enough to the right, given all waves in this solution are traveling to the right), then the initial state is not yet disturbed, and so the EXACT solution there just is the initial condition there. The same holds true if the grid speed were far enough to the left, such that the entire system of generated waves has not interacted with the initial state. Now the tricky thing is to get the solution inside the wave structure, and honestly, it seems like something you'd ever only want to do with the help of a computer program, looping through tests and conditions to determine the new values of the disturbed initial condition. Here is the plot from Toro of the exact solution of Test 4 using the Exact Riemann Solver: |