The dynamical evolution of neutrino flavor in supernovae can be modeled by an all-to-all spin Hamiltonian with random couplings. Simulating such two-local Hamiltonian dynamics remains a major challenge, as methods with controllable accuracy require circuit depths that increase at least linearly with system size, exceeding the capabilities of current quantum devices. The eigenstate thermalization hypothesis predicts that these systems should thermalize, a behavior confirmed in small-scale classical simulations. In this work, we investigate flavor thermalization in much larger systems using random quantum circuits as an empirical tool to emulate the non-local dynamics, and demonstrate that the thermal behavior can be reproduced using a depth independent of the system size. By simulating dynamics of over one hundred qubits, we find that the thermalization time grows approximately as the square root of the system size, consistent with predictions from semi-classical methods. Beyond this specific result, our study illustrates that near-term quantum devices are useful tools to test and validate empirical classical methods. It also highlights a new application of random circuits in physics, providing insight into complex many-body dynamics that are classically intractable.