http://xyst.i2p/kozcg_example_sod.html
Dirichlet boundary conditions were applied on all flow variables sampling their initial values at the end-faces of the domain and symmetry (free-slip) conditions were set along the surface of the tube. The numerical solution was advanced using the CFL number 0.5 using four different meshes whose properties are list in the following table: Mesh Points Tetrahedra h 0 7,294 34,181 0.008393 1 51,794 273,448 0.004138 2 389,139 2,187,584 0.002060 3 3,014,277 17,500,672 0.001028...