This work is concerned with theoretical study of hydrodynamic journal bearings lubricated with ferrofluids exhibiting non-Newtonian behavior. Based on the momentum and continuity equations for ferrofluid under an applied magnetic field, a modified Reynolds equation has been obtained. Assuming linear behavior for the magnetic material of the ferrofluid, the magnetic force was calculated. The Reynolds equation has been derived to be able to apply to any magnetic field distribution model. Using different magnetic field models, the equation has been solved numerically by the finite-difference technique with appropriate iterative technique and pressure distributions have been obtained. The boundary shape of the load-carrying active regions (positive-pressure regions) and cavitation regions (zero-pressure regions) could be then determined. The solution renders the bearing performance characteristics, namely: load-carrying capacity, attitude angle of the journal center, frictional force at the journal surface, friction coefficient and bearing side leakage. The results indicated that the flow-behavior index has a large effect on the bearing performance. When the bearing operates at high eccentricity ratios, the increase of flow-behavior index gives higher load capacity, lower attitude angle, higher frictional force, lower friction coefficient and higher side leakage. At low eccentricity ratios where the magnetic effects are significant, the effect of the flow-behavior index depends mainly on the magnetic field distribution model used.