A deterministic Boltzmann model equation solver called dugksFoam has been developed in the framework of the open source CFD toolbox OpenFOAM. The solver adopts the discrete unified gas kinetic scheme (Guo et al., 2015) with the Shakhov collision model. It has been validated by simulating several test cases covering different flow regimes including the one dimensional shock tube problem, a two dimensional thermal induced flow and the three dimensional lid-driven cavity flow. The solver features a parallel computing ability based on the velocity space decomposition, which is different from the physical space decomposition based approach provided by the OpenFOAM framework. The two decomposition approaches have been compared in both two and three dimensional cases. The parallel performance improves significantly using the newly implemented approach. A speed up by two orders of magnitudes has been observed using 256 cores on a small cluster. Program summary Program Title: dugksFoam Program Files doi: http://dx.doi.org/10.17632/zwn7t9cf5w.1 Licensing provisions: The MIT License Programming language: C++ External routines/libraries: OpenFOAM (http://www.openfoam.org) Nature of problem: Solving the Boltzmann equation with Shakhov model explicitly. Solution method: Discrete unified gas kinetic scheme (DUGKS) Restrictions: Symmetric boundary condition can only be applied at walls parallel to axis directions. (C) 2016 Elsevier B.V. All rights reserved.