Gallium Nitride (GaN) high electron mobility transistors (HEMTs) exhibit superior electrical properties for power and radio frequency applications, but performance is compromised by localized Joule heating, increasing channel temperatures. Precise thermal analysis during design is essential for optimizing device architecture and management strategies. Traditional methods like the Fourier heat diffusion equation (HDE) and the phonon lattice Boltzmann method (PLBM) with the D2Q8 scheme inadequately model phonon ballistic transport at high Knudsen numbers. This study introduces a nongray multi-speed PLBM integrated with the drift-diffusion model to analyze electro-thermal processes in GaN HEMTs. Validated through simulations of thermal conductivities in two-dimensional GaN thin films, the approach examines internal temperature rise in GaN HEMTs under different gate voltages, comparing results with HDE and gray BTE models, and emphasizing the need for non-Fourier effects in thermal analysis. It also evaluates the impact of GaN layer thickness on temperature distribution, providing a robust solution for thermal analysis in GaN HEMTs and other field-effect transistors.