Shock Tube Simulation by the Smooth Particle Hydrodynamic (SPH) Method.

摘要

A smooth particle hydrodynamic code (SPHC) is evaluated for performing shock wave simulations by application to 1-D shock tube problems. Results of a shock tube test case with a compression ratio of 10 are compared against a Riemann shock tube problem and theoretical predictions of shock tube behavior to validate the SPH code. A Lagrangian hydrodynamic code is validated in a similar fashion. The resolution capabilities of both codes are compared using 100, 200 and 500 particles for SPHC and 100, 400 and 800 cells for the Lagrangian code. The SPH code exhibits a sharp spike in density at the contact discontinuity for a shock tube test case run with a compression ratio of 100. This behavior is not reported in the literature and not seen in the Lagrangian code results. Run time scaling is investigated for both codes. SPHC is found to scale between N log N and N2, where N is the number of particles. The Lagrangian code scales O(N2). Computation times for the SPH code are greater than run times for the Lagrangian code by a factor of four for N 500 to achieve similar resolution.