In order to improve the accuracy of numerical simulation of magnetic fluid thermal coupling of cylindrical heat source, a numerical simulation method of magnetic fluid thermal coupling of cylindrical heat source based on LBM grid algorithm was studied. Based on the theoretical model of natural convection and radiative heat dissipation of cylindrical heat source, the boundary conditions were set, and the thermal coupling performance of flow field and temperature field in the process of natural convection and radiative heat dissipation of cylindrical heat source was simulated by LBM grid algorithm. Simulation results show that the relative error between simulation results and experimental results is less than 1.5%. The results show that this method can accurately simulate the magnetofluid-thermal coupling of cylindrical heat source. With the increase of Rayleigh number, different numerical solutions appear: steady state solution when Rayleigh number is 105, periodic oscillation solution when Rayleigh number is 106, chaotic solution when Rayleigh number is 107.