 Solution

For stable shear-driven turbulence, let $P$ be the rate of production of turbulent kinetic energy by mean shear and $B_m>0$ the rate consumed by buoyancy work. The Flux Richardson number isWith an upward coordinate and vertical turbulent buoyancy flux $\overline{w^{\prime }{b}^{\prime }}$, $B_m=-{\rho }_l\int \overline{w^{\prime }{b}^{\prime }}d\mathcal{V}$, while shear production is $P=-{\rho }_l\int \overline{u^{\prime }{w}^{\prime }}{\partial }_{z}Ud\mathcal{V}$. In a locally steady balance without net turbulent-energy transport or storage, $P=B_m+\mathcal{D}$, so $0\le {\mathrm{Ri}}_f\le 1$ and the dissipation-normalized mixing coefficient is $B_m/\mathcal{D}={\mathrm{Ri}}_f/(1-{\mathrm{Ri}}_f)$. A global ratio $C_W={\dot{E}}_P/P_I$ need not equal this local production-normalized Flux Richardson number when energy transport, storage, or bulk acceleration is important.

The interfacial Richardson number in this tank isIn model S,Thus the relative importance of interface stratification changes continually. If mixing efficiency depends on interfacial Richardson number, constant $C_W$ is not generally self-consistent over substantial deepening; it can at best describe a regime where that dependence is weak. In model K,A local efficiency depending only on this ratio can therefore remain constant, provided the fixed interface thickness and the other empirical coefficients remain appropriate.

The energetic reason to expect intermediate behavior is that entrainment must both lift and homogenize dense fluid and mobilize it. The whole upper-layer kinetic energy budget has the form $P_{\mathrm{lid}}=\dot{K}+{\dot{E}}_P+\mathcal{D}$, apart from additional boundary losses. Keeping the characteristic speed constant makes $K\propto h$ and requires a continuing power supply to accelerate new fluid, together with dissipation in an ever larger layer. Keeping $K$ constant instead forces the speed to decline as $h^{-1/2}$. A finite lid forcing plausibly produces an initial regime near the former and a deeper regime near the latter, rather than indefinite constant-speed entrainment.

One can make this expected bracket quantitative after matching the models' initial interface speed $u_0$. Suppose the actual speed lies between $u_0\sqrt{h_0/h}$ and $u_0$, and the same local coefficients are used. Set $D_0=2C_W{c}_D{u}_0^3/B$. ThenIntegrating these inequalities givesEquivalently, if $K\propto h^q$ with $0\le q\le 1$, the local power closure yields a late-time depth exponent $2/(5-3q)$ between $2/5$ and $1$. This is the intended energetic interpolation between K and S, not a universal bound following from the printed assumptions alone. Different efficiencies or unmatched initial coefficients can invalidate pointwise ordering; quantitative prediction requires the lid-to-interface power transfer and the actual energy budget.