A numerical simulation scheme based on a three-dimensional boundary integral method is used for the purpose of estimating the macroscopic flow from the microscale drop deformation for emulsions of high viscosity ratio under steady, oscillatory shear, and extensional flows. The effects of drop deformation and shear rate are examined for moderate and high viscosity ratio. The accuracy and convergence of the numerical simulations is illustrated by comparing drop deformation and rheology with experimental observations and an asymptotic theory based on a second-order small deformation theory also developed in this work. We also examine an emulsion undergoing a pure extensional flow and an extensional viscosity of the emulsion is computed by the theory and the numerical simulations. In addition, we also explore regimes of linear and nonlinear viscoelasticity when the emulsion is subjected to oscillatory shear for different viscosity ratio. Emulsion Fourier modes are examined for two different viscosity ratio. For a more concentrate emulsion with an unit viscosity ratio, a static shear elastic modulus is computed as a function of drop volume fraction and compared with experimental data. The numerical simulation results are validated against theory and experimental observation of drops in shear and extensional flows. A very good agreement is observed for moderate and high viscosity ratio. The regime of emulsions with relatively high viscosity ratios has not been much explored in the current theoretical and experimental literature. The boundary integral numerical simulations developed in this work has the ability to predict what we expect on the behavior of emulsion with moderate and high viscosity ratio.