Nonthermal radiation has been seen in the winds around O-stars and Wolf-Rayet stars in the radio which is attributed to shock accelerated cosmic rays.  The shocks thought to be doing the accelerating are either internal shocks in the wind due to instabilities or shocks due to a colliding wind binary system.  Currently very few numerical Diffusive Shock Acceleration (DSA) simulations exist for these systems due to the complicated, multidimensional nature of the winds.  We present  the first 2-D Magnetohydrodynamic DSA (MHD-DSA) simulations of massive stellar winds using AstroBEAR (Astronomical Boundary Embedded Adaptive Refinement) which has recently been modified to include cosmic rays.  This work is supported at the University of Minnesota by NSF, NASA and the Minnesota Supercomputing Institute.