Stellar wind bow shocks are structures created when stellar winds with supersonic relative velocities interact with the local interstellar medium (ISM). They can be studied to understand the properties of stars as well as the ISM. Since bow shocks are asymmetric, light becomes polarized by scattering in the regions of enhanced density they create. We use a Monte Carlo radiative transfer code calle SLIP to simulate the polarization signatures produced by both resolved and unresolved bow shocks with analytically derived shapes and density structures. When electron scattering is the polarizing mechanism, we find that optical depth plays an important role in the polarization signatures. While results for low optical depths reproduce theoretical predictions, higher optical depths produce higher polarization and position angle rotations at specific viewing angles. This is due to the geometrical properties of the bow shock along with multiple scattering effects. For dust scattering, we find that the polarization signature is strongly affected by wavelength, dust size, dust composition, and viewing angle. Depending on the viewing angle, the polarization magnitude may increase or decrease as a function of wavelength. We will present results from these simulations and preliminary comparisons with observational data.