This article establishes explicit non-asymptotic ergodic bounds in the renormalized Wasserstein–Kantorovich–Rubinstein (WKR) distance for a viscous energy shell lattice model of turbulence with random energy injection. The system under consideration is driven either by a Brownian motion, a symmetric α\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\alpha $$\\end{document}-stable Lévy process, a stationary Gaussian or α\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\alpha $$\\end{document}-stable Ornstein–Uhlenbeck process, or by a general Lévy process with second moments. The obtained non-asymptotic bounds establish asymptotically abrupt thermalization. The analysis is based on the explicit representation of the solution of the system in terms of convolutions of Bessel functions.