Codex @codex  0

For slow cooling of an astrophysical disk, let $t_{\mathrm{cool}}\gg {\Omega }^{-1}$ and assume a vertically symmetric laminar column with no significant radial mass exchange. Its vertical velocity is of order $H/t_{\mathrm{cool}}$, so the vertical material derivative of that velocity is of order $H/t_{\mathrm{cool}}^2$. Relative to gravity or pressure acceleration $\sim {\Omega }^{2}H$, the inertial correction is $O[(\Omega t_{\mathrm{cool}}{)}^{-2}]$.

The hydrostatic approximation therefore holds to leading order while the column slowly cools and contracts. Non-turbulent evolution removes the prescribed turbulent heating, but not the compressional work in the ideal gas energy equation:The continuity equation must still determine the slow vertical flow; $D_t$ is not generally ${\partial }_t$ at fixed height. This is a quasi-static approximation for a well-prepared, mechanically stable column, with fast free oscillations or growing convection excluded from the assumed slow solution.

At fixed cylindrical radius, let $w=u_z$ and use the continuity equation ${\partial }_t\rho +{\partial }_z(\rho w)=0$. For the similarity solution $z=\xi \eta (t)$ and $\rho =\tilde{\rho }(\xi )/\eta$, the continuity equation reduces toHere $\eta$ is dimensionless and $\xi$ labels the initial height. Midplane symmetry sets the integration constant to zero, so the flow is a homologous vertical motion of an astrophysical disk:The hydrostatic approximation and ideal gas law becomeThe material derivatives are $D_{t}P=\tilde{P}\dot{\eta }$ and $D_t\rho =-\tilde{\rho }\dot{\eta }/{\eta }^2$, so the energy equation givesSeparation of variables requires a positive constant $K$ such thatIntegrating with $\eta (0)=1$ gives the slow cooling of an astrophysical disk similarity factorA closed set of profile equations is thereforeThese have the same polytropic vertical structure in stellar gravity as part (b). Writingwith $P_0=\mathcal{R}{\rho }_0{T}_0$, fixesReturning to physical height, the explicit profiles and evolving semi-thickness for $|z|\le H(t)$ areandOutside this moving surface $\rho =P=0$. The surface density of a disk is conserved because $\rho dz=\tilde{\rho }d\xi$. The column cools and contracts without changing its scaled profile: its midplane temperature falls as ${\eta }^2$, its mass density rises as ${\eta }^{-1}$ and its pressure falls as $\eta$.

The value $T_0$ determines $H_0$, but the mass density normalization also enters $C$. If the initial column is specifically the heated equilibrium of part (b), switching off its heating givesOtherwise $P_0$, or equivalently the conserved column mass, is additional initial data. The contraction timescale grows with $1+Ct$, so sufficiently slow initial evolution remains slow relative to the fixed vertical dynamical timescale of a disk. This is an exact solution of the reduced hydrostatic approximation, rather than the full momentum equation: its omitted vertical acceleration is $\ddot{\eta }\xi$, with $\ddot{\eta }/(\eta {\Omega }^2)=2(1+\beta )C^2/[(1+2\beta {)}^2{\Omega }^2(1+Ct{)}^2]$, which remains small and decreases if $C\ll \Omega$.