Ferrofluid flow is fascinating since their fluid properties can conveniently be manipulated by external magnetic fields. Novel applications in micro-and nanofluidics as well as in biomedicine have renewed the interest in the flow of colloidal magnetic nanoparticles with a focus on small-scale behavior. Traditional flow simulations of ferrofluids, however, often use simplified constitutive models and do not include fluctuations that are relevant at small scales. Here, we address these challenges by proposing a hybrid scheme that combines the multi-particle collision dynamics method for modelling hydrodynamics with Brownian Dynamics simulations of a reliable kinetic model describing the microstructure, magnetization dynamics and resulting stresses. Since both, multi-particle collision dynamics and Brownian Dynamics are mesoscopic methods that naturally include fluctuations, this hybrid scheme presents a promising alternative to more traditional approaches, also because of the flexibility to model different geometries and modifying the constitutive model. The scheme was tested in several ways. Poiseuille flow was simulated for various model parameters and effective viscosities were determined from the resulting flow profiles. The effective, field-dependent viscosities are found to be in very good agreement with theoretical predictions. We also study Stokes' second flow problem for ferrofluids. For weak amplitudes and low frequencies of the oscillating plate, we find that the velocity profiles are well described by the result for Newtonian fluids at the corresponding, field-dependent viscosity. Furthermore, the time-dependent profiles of the nonequilibrum magnetization component are well approximated by their steady-state values in stationary shear, when evaluated with the instantaneous local shear rate. Finally, we also apply the new scheme to simulate ferrofluid shear flow over rough surface. We find characteristic differences in the nonequilibrium magnetization components when the external field is oriented in flow and in gradient direction.