This thesis studies Monte Carlo simulations of QCD heavy flavor production processes (p{bar p} {yields} Q({anti Q})X) at hadron colliders. ISAJET bottom quark cross-sections are compared to the O({alpha} {sub s}{sup 3}) perturbative calculation of Nason, Dawson, and Ellis. These Monte Carlo cross-sections are computed from data samples which use different parton distribution functions and physics parameters. Distributions are presented in the heavy quark`s transverse momentum and rapidity. Correlations in rapidity and azimuthal angle are computed for the heavy flavor pair. Theory issues which arise are the behavior of the cross-section at low and high values of transverse momentum and the treatment of double counting problems in the flavor excitation samples. An important result is that ISAJET overestimates bottom quark production cross-sections and K factors. These findings are relevant for estimates of rates and backgrounds of heavy floor events.