 Solution

Let $a=({\rho }_b-{\rho }_t)/H={\rho }_0{N}^2/g$. The initial density is $\rho (z)={\rho }_t-az$. Conservation of mass in the material homogenized over $-h<z<0$ fixes its mean density to ${\rho }_u={\rho }_t+ah/2$, while the untouched lower layer has interfacial density ${\rho }_i={\rho }_t+ah$. Hence the reduced gravity difference across the interface isThe reference contribution involving ${\rho }_t$ cancels from the difference.

Only the upper layer changes its potential energy. The potential-energy cost of homogenizing a linear stratification issoThe positive sign reflects work against stable stratification.

The rotating boundary supplies energy on scales set by $\Omega$ and $R$. Turbulent drag and dissipation remove energy, and entrainment adds initially stationary fluid that must be accelerated. If upper-layer energy grows, the increasing turbulent loss provides a restoring tendency; after a short adjustment, input and losses can approximately balance while the interface moves on a slower time scale. This gives a physical rationale for a fixed-energy turbulent mixed layer, not a consequence of mass conservation alone. Its characteristic kinetic energy is $K\sim {\rho }_0\pi R^{2}h{u}_u^2/2$. With the stipulated $N$-independent energy closure,so the characteristic speed decreases as $h^{-1/2}$ rather than remaining constant.

The stress does work at a rate proportional to $\pi R^2\tau u_u=c_D{\rho }_0\pi R^2{u}_u^3$. Combine this with the energy increase and absorb fixed drag/geometrical factors into $C>0$:Here the exponents $\alpha ,\beta$ are unrelated to the plume entrainment coefficient in Question 2. The interfacial Richardson number isSet $q=\Omega /N$, $d=\delta /R$ and $r=h/R$. Then $u_u^2/(N^2{h}^2)=d/(2r{\mathrm{Ri}}_i)$ and $r=q(2A{\mathrm{Ri}}_i/d{)}^{1/2}$. Substitution givesThe assumption that $E$ depends on the local Richardson number alone excludes separate dependences on $q,d$. The entrainment exponent from a local interfacial Richardson closure is consequentlyThis comparison treats $q$ and $d$ as independently variable external controls; fitting only a single apparatus would not by itself determine the exponents.

Finally, use the fixed-energy closure again:Integrating from a positive initial mixed-layer depth $h_0$ at time $t_0$ givesWhen the growth term dominates the initial offset, the rotating-disc mixed-layer depth law isTaking the effective time origin to be zero gives the stated scaling. Its validity is limited to the regime of the closures and $h<H$; the singular velocity predicted at $h=0$ is not a model for the initial boundary-layer formation.