Recently several works have been published on numerical simulation of an external gear pump (EGP). Such kinds of pumps are simple and relatively inexpensive, and are frequently used in fluid power applications, such as fluid power in aeronautical, mechanical, and civil engineering. Nevertheless, considerable effort is being undertaken to improve efficiency and reduce noise and vibration produced by the flow and pressure pulsations. Numerical simulation of an EGP is not straightforward principally for two main reasons. First, the gearing mechanism between gears makes it difficult to handle a dynamic mesh without a considerable deterioration of mesh quality. Second, the dynamic metal–metal contact simulation is important when high pressure outflow has to be reproduced. The numerical studies published so far are based on a two-dimensional (2D) approximation. The aim of the present work is to contribute to the understanding of the fluid flow inside an EGP by means of a complete three-dimensional (3D) parallel simulation on a cluster. The 3D flow is simulated in a linux cluster with a solver developed with the openfoam Toolbox. The hexahedral mesh quality is maintained by periodically replacing the mesh and interpolating the physical magnitudes fields. The meshing contact point is simulated with the viscous wall approach, using a viscosity model based on wall proximity. The results for the flow rate ripples show a similar behavior to that obtained with 2D simulations. However, the flow presents important differences inside the suction and the discharge chambers, principally in the regions of the pipes' connection. Moreover, the decompression slot below the gearing zone, which can not be simulated with a 2D approximation, enables a more realistic simulation of a contact ratio greater than 1. The results are compared with experimental measurements recently published.