Solution (source code)

= Solution

For $\beta=1/2$, the exact profiles are
$$
T=T_0(1-x^2)^{1/2},
\qquad
\rho=\rho_0(1-x^2),
\qquad
P=P_0(1-x^2)^2,
\qquad x=\frac zH.
$$
Integration through both disk faces gives
$$
\boxed{\Sigma=\rho_0H\int_{-1}^1(1-x^2)\,dx
=\frac43\rho_0H}.
$$
Since $\mu=\alpha P/\Omega$ and $\bar\nu\Sigma=\int\mu\,dz$,
$$
\boxed{\bar\nu\Sigma
=\frac{\alpha P_0H}{\Omega}
\int_{-1}^1(1-x^2)^2\,dx
=\frac{16}{15}\frac{\alpha P_0H}{\Omega}}.
$$