We describe the details of numerical implementation of the multidimensional kinetic solver using kinetic and total energy domains. This solver provides a substantially lower computational cost than the Monte Carlo method, and resolves all points in phase space with equal accuracy. The effect of quantum corrections on spatial distribution of carriers is analyzed. Simulation results for an ultra-small MOSFET will be presented.