  We present the first numerical radiative transfer simulation of multiple
light scattering in dust configurations containing aligned non-spherical
(spheroidal) dust grains. Such models are especially important if one wants to
explain the circular polarization of light, observed in a variety of
astronomical objects. The radiative transfer problem is solved on the basis of
the Monte Carlo method. Test simulations, confirming the correct numerical
implementation of the scattering mechanism, are presented. As a first
application, we investigate the linear and circular polarization of light
coming from a spherical circumstellar shell. This shell contains perfectly
aligned prolate or oblate spheroidal grains. The most remarkable features of
the simulated linear polarization maps are so-called polarization null points
where the reversal of polarization occurs. They appear in the case when the
grain alignment axis is perpendicular to the line of sight. The maps of
circular polarization have a sector-like structure with maxima at the ends of
lines inclined to the grain alignment axis by $\pm 45\degr$.
