A grad‐div stabilized method using the Jacobi iteration for the thermally coupled incompressible magnetohydrodynamic system