This paper documents the 2D numerical study of magnetohydrodynamic unsteady natural convective heat transfer in a circular enclosure with four heating cylinders in both the horizontal and the vertical mid-plane. The fluid is an incompressible Newtonian fluid. The main transport equations based on the conservation of mass, momentum, and energy are calculated and solved using a finite element numerical solver with the following parameter ranges: dimensionless distance between cylinders S = 0.05–0.29, Rayleigh number Ra=103–106, and Hartmann number for Ha = 0–120. COMSOL Multiphysics, a numerical simulation program, was used to solve the governing equations. It was demonstrated that for lower Ra values, heat transfer through an applied magnetic field is unaffected for a specific S value because the mechanism of transport is diffusion, whereas for larger Ra, there is a complex interaction among magnetic field and physical thermal properties. The features of the heat transfer rate are determined by the interaction. The Nusselt number virtually stays constant as Ha rises at smaller Ra. However, at high Ra, the Nusselt number initially declines with Ha and thereafter essentially stays constant, and at high Ra values, the switch from conduction to convective heat transfer takes place. Additionally, Nu rises slightly with S at increasing Ra.