Numerical simulations of electromagnetic instabilities, driven by a cool tenuous ion beam propagating along an ambient magnetic field, have been conducted in one and two spatial dimensions. The calculations employ particle ions, fluid electrons, and a predictor corrector scheme for solving the electromagnetic field in two dimensions that is described in some detail. While the principle features of the one-dimensional calculations (which reproduce previous work) are retained, the two-dimensional simulations show some reduction of the overall level of the magnetic field fluctuations. Enhancement of the heating of the beam ions at the expense of the core ions also occurs in the case where the beam density is sufficiently large that the right-hand nonresonant instability dominates. Implications of the results for modeling of the ion foreshock and quasi-parallel shocks are discussed.