A high-order finite-difference scheme to model the fluid-structure interaction in pneumatic seismic sources

•The fluid-structure interaction inside a pneumatic seismic source is modeled by the nonlinear 1D Euler equations.•A provably stable numerical scheme modeling the fluid-structure interaction is derived.•Convergence studies confirm the accuracy of the scheme. The obtained convergence rates agree with theoretical results.•Simulation results captures many of the main features of the experimental data from a pneumatic seismic source.