  We use three dimensional magnetohydrodynamic (MHD) simulations to model the
supernova remnant SN 1006. From our numerical results, we have carried out a
polarization study, obtaining synthetic maps of the polarized intensity, the
Stokes parameter $Q$, and the polar-referenced angle, which can be compared
with observational results. Synthetic maps were computed considering two
possible particle acceleration mechanisms: quasi-parallel and
quasi-perpendicular. The comparison of synthetic maps of the Stokes parameter
$Q$ maps with observations proves to be a valuable tool to discern
unambiguously which mechanism is taking place in the remnant of SN 1006, giving
strong support to the quasi-parallel model.
