We present detailed numerical simulations of modal instabilities in high-power Yb-doped fiber amplifiers using a time-dependent temperature solver coupled to the optical fields and population inversion equations. The temperature is computed by solving the heat equation in polar coordinates using a 2...