The role of normal stresses in causing particle migration and macroscopic spatial variation of the particle volume fraction phi in a mixture of rigid neutrally buoyant spherical particles suspended in Newtonian fluid is examined for curvilinear shear flows. The problem is studied for monodisperse noncolloidal Stokes-flow suspensions, i.e., for conditions of low-Reynolds-number flow and infinite Peclet number, Pe = O(eta(gamma) over dot a(3)/k T ), where eta is the suspending fluid viscosity, (gamma) over dot is the shear rate, a is the particle radius, and k T is the thermal energy. Wide-gap Couette, parallel-plate torsional, and cone-and-plate torsional flows are studied. The entire phi dependence of the compressive shear-induced normal stresses is captured by a "normal stress viscosity" eta(n)(phi), which vanishes (as phi(2)) at phi = 0 and diverges at maximum packing in the same fashion as does the shear viscosity eta(s)(phi). Anisotropy of the normal stresses arising from the presence of the particles is modeled as independent of phi, so that ratios of any two particle contributions to the bulk normal stress components are constants, Sigma(22)(p)/Sigma(11)(p) = lambda(2) and Sigma(33)(p)/Sigma(11)(p) = lambda(3); the standard convention of (1,2,3) denoting the (flow, gradient, vorticity) directions is used so that, for example, Sigma(11)(p) the normal component of the particle stress Sigma(p) in the flow direction. Predictions for the steady and unsteady flows are presented to demonstrate the influence of variation of the normal stress anisotropy parameters lambda(2) and lambda(3), the rheological functions eta(s), and eta(n), and the sedimentation hindrance function used to represent the resistance to relative motions of the phases during migration. Comparison with available experimental data shows that a single set of parameters for the rheological model is able to describe all qualitative features of the observed migrations in the flows considered.