A suspension flow model based on the "suspension balance" approach has been developed. This work modifies the model to allow the solution of suspension flows under general flow conditions. This requires the development of a frame-invariant constitutive model for the particle stress which can take into account the spatially-varying local kinematic conditions. The mass and momentum balances for the bulk suspension and particle phase are solved numerically using a finite volume method. The particle stress is based upon the computed rate of strain and the local kinematic conditions. A nonlocal stress contribution corrects the continuum approximation of the particle phase for finite particle size effects. Local kinematic conditions are accounted through the local ratio of rotation to extension in the flow field. The coordinates for the stress definition are the local principal axes of the rate of strain field. The developed model is applied to a range of problems. (i) Axially-developing conduit flows are computed using both the full two-dimensional solution and the more computationally efficient "marching" method. The model predictions are compared to experimental results for cross-stream particle concentration profiles and axial development lengths. (ii) Model predictions are compared to experiments for wide-gap circular Couette flow of a concentrated suspension in a shear-thinning liquid. With minor modification, the suspension flow model predicts the major trends and results observed in this flow. (iii) Comparisons are made to experiments for an axisymmetric contraction-expansion. Model predictions for a two-dimensional planar contraction flow test the influence of model formulation. The variation of the magnitude of an isotropic particle normal stress with local kinematic conditions and anisotropy in the in-plane normal stresses are both explored. The formulation of the particle phase stress is found to have significant effects on the solid fraction and velocity. (iv) Finally, for a rectangular piston-driven flow and an obstructed channel flow, a "computational suspension dynamics" study explores the effect of particle migration on the bulk flow field, system pressure drop and particle phase composition.