In this study, the problem of an electrically conducting and steady incompressible micropolar fluid with a natural convection boundary layer flow in a sphere in the presence of a strong magnetic field has been considered. This type of flow has many engineering applications, such as drying processes, electronic cooling, ceramic processing, heat exchanger devices and many more. The effect of radiation, chemical reaction, magnetic source; permeability and heat source/sink on the velocity, temperature and concentration profiles of the fluid flow in the medium are investigated. The governing equations describing the problem are transformed into non-dimensional partial differential equations. The equations are solved with the recently introduced multi-domain spectral quasilinearization method. The method involves the use of quasilinearization in linearizing nonlinear functions, a domain splitting approach to split the time interval into several sub-intervals and the Chebyshev spectral collocation method in solving the system. The results obtained suggest that the various parameters tested have a significant impact on various flow profiles. While the velocity and temperature of the fluid increases with an increase in the heat source parameter, the concentration of the fluid are observed to increase with an increase in the radiation and permeability parameters, respectively. Increase in the chemical reaction parameter leads to a decrease in the fluid velocity around the boundary layer. The analysis performed suggests that the multidomain spectral quasilinearization method is efficient in solving fluid flow problems of similar dynamics.