Kinetic approaches are routinely employed to simulate the dynamics of systems that are too rarified to be described by the Navier-Stokes equations. However, generally they are far too computationally expensive to be applied for systems that are governed by continuum hydrodynamics. In this paper, we introduce a massively parallelized test-particle based kinetic Monte Carlo code that is capable of modeling the phase space evolution of an arbitrarily sized system that is free to move in and out of the continuum limit. Using particle mean free paths which are small with respect to the characteristic length scale of the simulated system, we retrieve continuum behavior, while non-equilibrium effects are observed when the mean free path is increased. To demonstrate the ability of our code to reproduce hydrodynamic solutions, we apply a test-suite of classic hydrodynamic shock problems. Simulations using tens of millions of test-particles are found to reproduce the analytical solutions well.

We report on the development of a test-particle based kinetic Monte Carlo code for large systems and its application to simulate matter in the continuum regime. Our code combines advantages of the Direct Simulation Monte Carlo and the Point-of-Closest-Approach methods to solve the collision integral of the Boltzmann equation. With that, we achieve a high spatial accuracy in simulations while maintaining computational feasibility when applying a large number of test-particles. The hybrid setup of our approach allows us to study systems which move in and out of the hydrodynamic regime, with low and high particle densities. To demonstrate our code's ability to reproduce hydrodynamic behavior we perform shock wave simulations and focus here on the Sedov blast wave test. The blast wave problem describes the evolution of a spherical expanding shock front and is an important verification problem for codes which are applied in astrophysical simulation, especially for approaches which aim to study core-collapse supernovae.