We have investigated how the addition of ions at rest (mass loading) affects the temperature of a flowing plasma in a magnetohydrodynamic (MHD) approximation, using analytic theory and time dependent, three‐dimensional MHD simulations of plasma flow past Io. The MHD equations show that the temperature can increase or decrease relative to the background, depending on the local sonic Mach number (MS) of the flow. For flows with MS > √9/5 ( when γ = 5/3), mass loading increases the plasma temperature. However, the simulations show a non‐linear response to the addition of mass. If the mass loading rate is large enough, the temperature increase may be smaller than expected, or the temperature may actually decrease, because a large mass loading rate slows the flow and decreases the thermal energy of the newly created plasma.