The dynamic viscosity under pressure of three mixtures of pentaerythritol ester lubricants (PEs) has been measured using a rolling-ball viscometer for several temperatures with an experimental uncertainty of 3%. The first one is a multicomponent mixture of several PEs named in the present work as PEC5-C9 lubricant; the second one is a binary mixture of pentaerythritol tetra(2-ethylhexanoate), PEB8, and pentaerythritol tetraheptanoate, PEC7, with a PEB8 mole fraction of 0.6670; and the third one is another binary mixture of PEB8 and pentaerythritol tetrapentanoate, PEC5, with a PEB8 mole fraction of 0.6911. The two binary mixtures, xPEB8 + (1 - x)PEC7 and xPEB8 + (1 -x)PEC5, have been prepared with the same viscosity grade as the PEC5-C9 lubricant (VG32). A total of 1176 experimental measurements of the rolling time have been performed at pressures up to 60 MPa for the determination of 196 dynamic viscosity data points. The viscosities of these binary mixtures have been compared with the predicted values obtained by using several viscosity models (Grunberg-Nissan and Katti-Chaudhri mixing laws, self-referencing model, hard-sphere theory, and free-volume model). All methods predict dynamic viscosity values for the two binary mixtures that agree with the experimental data within an average mean deviation of 10% over the entire temperature and pressure ranges. The best predictions were found with the free-volume model, for which the average mean deviation for both mixtures is lower than 4%. Parameter values for the self-referencing model were determined from experimental viscosity data of several pure PEs. These parameters permit the estimation of viscosity values of PE lubricant of unknown composition, when a viscosity value at any temperature and pressure is available. This model predicts the viscosities of PEC5-C9 lubricant with an average deviation of 4%.