The main approach for simulating the motion of fluids under various conditions has traditionally been to discretize the Navier-Stokes equations using finite element, finite differences, or finite volume methods. The basis for this project, is to study the fluid on a so-called mesoscopic level. Fluids are represented in terms of probability densities, and the Boltzmann equation will be used to study the behaviour. These models have been shown to converge to solutions of the Navier Stokes equation. The aim is to investigate if the convergence rate depends on the spatial dimension and the lattice scheme. The plan is to implement a parallel solver for the numerical method to be able to test this convergence in three space dimensions and under different conditions such as changing lattice and boundary conditions.