A coupled smoothed particle hydrodynamics-finite volume approach for shock capturing in one-dimension

A novel approach coupling the finite volume method code Pyro2 and the smoothed particle hydrodynamics code PySPH is introduced and applied to one-dimensional shock problems. The finite volume mesh models the bulk of the system, while regions with discontinuous fluid values are identified and populated with PySPH particles to model the fluid around shocks. The approaches are coupled with boundary cells and ghost particles, with linear interpolation used to extract fluid properties at each timestep in the respective boundary regions. Results from three shock problems using Pyro2, PySPH, and a hybrid approach are presented. The hybrid approach preserves accuracy and computational efficiency in strong shock problems. Additionally, the coupling approach provides a potential avenue for including SPH functionality in FVM simulations, including modeling complex geometries and fluid-structure interaction. The preliminary investigation of the hybrid approach in one dimension highlights the potential efficiency gains of applying coupled FVM-SPH methods to larger blast simulations in two and three dimensions.