This paper presents a new, efficient, accurate, and unconditionally stable second-order time-stepping method for the incompressible thermal micropolar Navier–Stokes equations (TMNSE) using mixed finite elements. The method linearizes the nonlinear convective terms in the momentum equation, microrotation equation, and temperature equation, requiring the solution of a linear problem at each time step. The discrete curvature of the solution is added as a stabilizing term for linear velocity u, microrotation velocity w, pressure p, and temperature T in the equations, respectively. Curvature stabilization (un+1−2un+un−1) is a new concept in computational fluid dynamics (CFD) aimed at improving the commonly used velocity stabilization (un+1−un), which only has first-order time accuracy and has adverse effects on important flow quantities such as drag coefficients. We derive a priori error estimates for the fully discrete linear extrapolation curvature stabilization method. The theoretical results and effectiveness of the new method are verified through a series of numerical experiments for θ=12, 34, 56, and 1 in 2D and 3D, respectively. In particular, this work considers the thermal cavity-driven flow experiment to validate the numerical scheme and obtains good results.