Proper functioning of white blood cells is not possible without their ability to adhere to vascular endothelium, which may occur only if they are close enough to vessel walls. To facilitate the adhesion, white blood cells migrate toward the vessel walls in blood flow through a process called margination. The margination of white cells depends on a number of conditions including local hematocrit, flow rate, red blood cell aggregation, and the deformability of both red and white cells. To better understand the margination process of white blood cells, we employ mesoscopic hydrodynamic simulations of a three-dimensional model of blood flow, which has been previously shown to capture quantitatively realistic blood flow properties and rheology. The margination properties of white blood cells are studied for a wide range of hematocrit values and flow conditions. Efficient white blood cell margination is found in an intermediate range of hematocrit values of Ht ≈ 0.2-0.4 and at relatively low flow rates, characteristic of the venular part of microcirculation. In addition, aggregation interactions between red blood cells lead to enhanced white-blood-cell margination. This simulation study provides a quantitative description of the margination of white blood cells, and is also highly relevant for the margination of particles or cells of similar size such as circulating tumor cells.