This work presents a model for the simulation of transverse mode instability (TMI) in rare earth doped optical fiber amplifiers. The model evaluates the internal temperature of a fiber using a superposition of a finite number of thermal eigenmodes. This simplification greatly enhances the speed of calculation with negligible impact on calculation accuracy. This new method is described and quantitatively compared to an older model that uses standard, spatially resolved FDTD to integrate the heat diffusion equation. When tested over a range of spatial and temporal resolutions, this model reduces runtime by a factor of ∼13.9 on average relative to identical simulations using the spatially resolved model.