In this paper, we present a cluster algorithm for the numerical simulations of nonadditive hard-core mixtures. This algorithm allows one to simulate and equilibrate systems with a number of particles two orders of magnitude larger than previous simulations. The phase separation for symmetric binary mixtures is studied for different nonadditivities as well as for the Widom-Rowlinson model [B. Widom and J. S. Rowlinson, J. Chem. Phys. 52, 1670 (1970)] in two and three dimensions. The critical densities are determined from finite size scaling. The critical exponents for all the nonadditivities are consistent with the Ising universality class. (C) 2005 American Institute of Physics.