Biomagnetic fluid dynamics (BFD) is an important application in medical sciences and bioengineering research. Due to this, biomagnetic fluid flow through a stenosed bifurcated artery is numerically studied. A biomagnetic fluid can be found in a living creature and its flow is influenced by the present of a magnetic field. Blood is a typical biomagnetic fluid due to the interaction of intercellular protein, cell membrane and the haemoglobin. This study considered the flow to be incompressible, laminar, two-dimensional (2D), fully developed viscous flow of a Newtonian biomagnetic fluid (blood) in a stenosed bifurcated artery under the effect of a spatially varying magnetic field. A simplified mathematical model of BFD was developed only for isothermal case. Numerical results are obtained using COMSOL Multiphysics 5.2 based on finite element method (FEM). Results concerning the different values of magnetic field intensity produce a considerable effect on the blood flow characteristics such as the velocity profiles and the streamlines patterns. It is shown that the vortex at lower wall extends vertically while vortex at upper wall becomes shrink as the magnetic field strength increases. The location of the magnetic source also can affect the velocity at the daughter artery where the velocity at the lower wall of daughter artery is lower than the upper wall.Biomagnetic fluid dynamics (BFD) is an important application in medical sciences and bioengineering research. Due to this, biomagnetic fluid flow through a stenosed bifurcated artery is numerically studied. A biomagnetic fluid can be found in a living creature and its flow is influenced by the present of a magnetic field. Blood is a typical biomagnetic fluid due to the interaction of intercellular protein, cell membrane and the haemoglobin. This study considered the flow to be incompressible, laminar, two-dimensional (2D), fully developed viscous flow of a Newtonian biomagnetic fluid (blood) in a stenosed bifurcated artery under the effect of a spatially varying magnetic field. A simplified mathematical model of BFD was developed only for isothermal case. Numerical results are obtained using COMSOL Multiphysics 5.2 based on finite element method (FEM). Results concerning the different values of magnetic field intensity produce a considerable effect on the blood flow characteristics such as the velocity pr...