  We propose two-equations models in order to capture the dynamics of a
turbulent plasma undergoing compression and experiencing large viscosity
variations.
  The models account for possible relaminarization phases and rapid viscosity
changes through closures dependent on the turbulent Reynolds and on the
viscosity Froude numbers. These closures are determined from a data-driven
approach using eddy-damped quasi normal markovian simulations. The best model
is able to mimic the various self-similar regimes identified in
\citet{Viciconte2018} and to recover the rapid transition limits identified by
\citet{Coleman1991}.
