Simulation of Shock Waves by Smoothed Particle Hydrodynamics

TL;DR

This paper presents a simulation method for shock waves using smoothed particle hydrodynamics (SPH), comparing results with analytical solutions and providing a flexible Fortran code adaptable to various scientific fields.

Contribution

It introduces a SPH-based simulation algorithm for shock waves, with a detailed program implementation and potential for diverse applications.

Findings

SPH accurately simulates isothermal and adiabatic shocks

Results align well with analytical solutions

Flexible Fortran code provided for various applications

Abstract

Isothermal and adiabatic shocks, which are produced from fast expansion of the gas, is simulated with smoothed particle hydrodynamics (SPH). The results are compared with the analytic solutions. The algorithm of the program is explained and the package, which is written in Fortran, is presented in the appendix of this paper. It is possible to change (to complete) the program for a wide variety of applications ranging from astrophysics to fluid mechanics.