We propose two-equation models in order to capture the dynamics of a turbulent plasma undergoing compression and experiencing large viscosity variations. The account for possible relaminarization phases rapid changes through closures dependent on Reynolds Froude numbers. These are determined from data-driven approach using eddy-damped quasi-normal Markovian simulations. best model is able mimic...