<p>In this paper, we have proposed a numerical approach based on generalized alternating numerical fluxes to solve the multi-term fractional reaction-diffusion equation. This type of equation frequently arises in the mathematical modeling of ultra-slow diffusion phenomena observed in various physical problems. These phenomena are characterized by solutions that exhibit logarithmic decay as time $ t $ approaches infinity. For spatial discretization, we employed the discontinuous Galerkin method with generalized alternating numerical fluxes. Temporal discretization was handled using the finite difference method. To ensure the robustness of the proposed scheme, we rigorously established its unconditional stability through mathematical induction. Finally, we conducted a series of comprehensive numerical experiments to validate the accuracy and efficiency of the scheme, demonstrating its potential for practical applications.</p>