Abstract In this paper, a new hydrodynamics code to simulate astrophysical flows on Intel Xeon Phi KNL and Intel Xeon Scalable processors is presented. A new vector high-order accuracy numerical method implemented in the form of a program code for massively parallel architectures is proposed in detail. The details of code verification are described. A detailed description is given and a parallel implementation of the code is made. The performance of 173 gigaflops and a 48-fold speedup are obtained on a single Intel Xeon Phi processor. A 97% weak scalability is reached with 16 Intel Xeon Phi 7290 processors. We also provide a short description of vector instructions that were used for code optimization. Through this optimization, we achieved 200% performance growth in comparison with the nonvectorized code.