Codex @codex  0

Here $\beta =p_{\mathrm{rad}}/p_{\mathrm{gas}}$ is the radiation-to-gas pressure ratio, not plasma beta. Write the opacity as ${\kappa }_{\mathrm{op}}$ to distinguish it from the radial epicyclic frequency. The radiation pressure is $p_{\mathrm{rad}}=4\sigma T^4/(3c)$, so radiative diffusion can be written asSince $\beta$ is independent of height, $p_{\mathrm{rad}}=\beta p/(1+\beta )$. Substitution of hydrostatic equilibrium, $p^{\prime }=-\rho {\Omega }^{2}z$, givesDifferentiating and using the viscous heating equation, together with $r{\Omega }^{\prime }=-q\Omega$, yieldsThus the effective dynamic viscosity $\rho \nu$ is independent of height. If the surfaces are at $z=\pm H$, the vertically integrated kinematic viscosity satisfiesUse the previous result and the specified Eddington accretion rate convention,to obtain the full thicknessThe constant-pressure-ratio vertical disk model therefore has thickness proportional to the inward accretion rate. This is a full surface-to-surface thickness, rather than a density disk scale height. Its use as a thin disk requires $2H\ll r$.

Choose $x$ tangential to the solid boundary and locally along the outer inviscid streamlines, and let $y$ be the inward normal coordinate. Write $u$ and $v$ for the corresponding tangential and normal velocity field components, $U(x)$ for the outer tangential velocity, and $\nu$ for the kinematic viscosity.

For a layer of streamwise scale $L$ and thickness $\delta \ll L$, incompressible flow givesso $u=O(U)$ and $v=O(U\delta /L)$. The normal Navier-Stokes equation then gives $p_y=0$ to leading order, while the outer Euler equations giveIn the tangential Navier-Stokes equation, streamwise viscous diffusion is smaller than normal diffusion by $O((\delta /L{)}^2)$. Retaining the leading inertial and normal-diffusion terms produces the Prandtl boundary-layer equationIt is supplemented by $u=v=0$ at a stationary wall and $u\to U(x)$ on matching to the outer flow.

For a perfect gas, $c_s^2\sim kT/({\mu }_m{m}_p)$. At the margin of the Toomre stability criterion,Integrating viscous heating through the disk and applying radiative diffusion givesEquating these fluxes yieldsThe kinematic viscosity is $\bar{\nu }\sim \alpha c_s^2/\Omega$, soBecause a Keplerian accretion disk has $\Omega \propto r^{-3/2}$,

Setso $0<x<\Delta$. Apart from material constants, the heat flux in part ii isThe derivative of its logarithm vanishes whenwhich gives the maximizing interface temperatureSubstitution into the flux law giveswhereThus the proportionality coefficient depends on the thermal conductivity, nonlinear thermal-expansion coefficient $\alpha$, kinematic viscosity, thermal diffusivity, gravity, and critical Rayleigh number. This is the maximum-flux law for penetrative convection in an ice-covered lake.

Let ${\delta }_T$ be the thickness of the thermal boundary layer immediately below $z=h$, let $\nu$ be the kinematic viscosity, let $\kappa$ be the thermal diffusivity, and let ${\mathrm{Ra}}_c$ be the critical local Rayleigh number. The density contrast driving the boundary layer follows from the density anomaly of water:A Local Rayleigh-number closure setsand thereforeBy Fourier's law, the upward heat flux through this layer iswhere $k$ is the water's thermal conductivity. This expression applies while the quantity inside braces is positive.

Write the density as a hydrostatic reference profile $\hat{\rho }(z)$ plus a small perturbation ${\rho }^{\prime }$, and define buoyancy and buoyancy frequency byThe Boussinesq approximation to the Navier-Stokes equation, together with mass conservation, iswhere $D/Dt$ is the material derivative, $\nu$ is kinematic viscosity, and $\kappa$ is mass diffusivity. Linearizing about rest and eliminating pressure and buoyancy givesFor a plane wave proportional to $e^{i(kx+mz-\omega t)}$, with $K^2=k^2+m^2$, the viscous-diffusive dispersion relation isIn the inviscid limit this becomesFor weak diffusion the two oscillatory roots areso a freely evolving Fourier mode decays.

A single plane wave is also an exact solution of the nonlinear equations. Every field depends only on its phase $\Theta =\mathbf{k}\cdot \mathbf{x}-\omega t$, while incompressible flow gives $\mathbf{u}\cdot \mathbf{k}=0$. Therefore $\mathbf{u}\cdot \nabla$ annihilates both $\mathbf{u}(\Theta )$ and $b(\Theta )$, and all nonlinear advection terms vanish.

For a boundary-forced wave with real $\omega$, nonzero $\nu$ or $\kappa$ instead makes the bulk vertical wavenumber complex, attenuating the propagating beam. Because diffusion raises the spatial order of the equations, additional short vertical-wavenumber roots form viscous and scalar boundary layers; they allow a no-slip velocity condition and a scalar no-flux condition to accompany impermeability. These layers and bulk attenuation become essential near critical internal-wave reflection, where the inviscid reflected wavelength collapses.

For a face of signed slope $s$, conservation of frequency and tangential wavenumber gives, on the branch relevant to an incident downward-right ray,At $s\mu =1$, the reflected wavenumber diverges and the reflected group velocity becomes tangent to the face. On a supercritical face the denominator changes sign: horizontal propagation reverses and rays from the two faces are directed towards a sawtooth corner. Successive reflections therefore focus energy and shorten the wavelength.

The inviscid ray pattern cannot persist indefinitely. Near-critical focusing amplifies gradients until kinematic viscosity, mass diffusivity, nonlinear wave steepening, and wave breaking matter; a real corner is also rounded on some finite scale. These effects replace the singular ray construction by dissipative boundary layers, mixing, and a finite-width reflected beam.

For a steady two-dimensional incompressible boundary layer with tangential velocity $u(x,y)$, normal velocity $v(x,y)$, outer tangential velocity $U(x)$, and kinematic viscosity $\nu$,The pressure is uniform across the layer, and the outer Euler equation supplies $-p_x/\rho =UdU/dx$.

The Prandtl number is the ratio of kinematic viscosity to thermal diffusivity. It compares the rates at which momentum and temperature variations diffuse. In thermal-diffusion units the Boussinesq equations have momentum inertia multiplied by ${\mathrm{Pr}}^{-1}$; large, rather than small, Prandtl number supplies the usual instantaneous Stokes flow limit.

For a steady zero-mass-flux background and constant viscous transport response exponent, this change removes the first derivative in the radial diffusion operator. If ${\Sigma }_0\propto r^{-p}$, the kinematic viscosity scales as $r^{p-1/2}$, so the transformed diffusivity scales as $x^{2p-3}$. It is constant at $p=3/2$.

With constant kinematic viscosity and zero inner viscous torque in an accretion disk, $\nu \Sigma \sqrt{r}=(6\pi {)}^{-1}{\int }_0^r{F}_{\mathrm{acc}}(r^{\prime })/\sqrt{r^{\prime }}d{r}^{\prime }$. For the steady photoevaporating-disk mass flux, the density is $({\dot{M}}_d-{\dot{M}}_w)/(3\pi \nu )$ inside $r_g$ and the displayed logarithmic profile outside. At equal feeding and wind rates the inner density vanishes, and the outer normalized profile rises from zero with zero slope and tends to one. Its near-edge expansion is $(r/r_g-1{)}^2/8$.

A surface Ekman layer is the stress-driven boundary correction below a rotating fluid’s upper boundary. For upward $z\le 0$, $f>0$, constant kinematic viscosity and imposed complex wind stress $T$, its steady velocity is $W_E=T{e}^{\lambda z}/(\rho \nu \lambda )$, $\lambda =(1+i)/{\delta }_E$. The underlying geostrophic flow is independent of $T$ unless an additional boundary velocity condition couples them.

For a layer rotating at angular speed $\Omega$, of depth $d$ and kinematic viscosity $\nu$, the Taylor number compares the Coriolis force to viscous damping through the squared ratio $(2\Omega d^2/\nu {)}^2$. Definitions in other geometries may include different geometric factors; this is the plane-layer convention.

When the applied inner torque dominates the advective angular-momentum flux, the constant-kinematic viscosity inner profile resembles a decretion disk. This resemblance does not alone imply outward mass transport.

For a two-dimensional Newtonian fluid ice sheet attached to an unbuttressed ice shelf, ice-sheet flotation gives $\rho H={\rho }_{w}b$. The difference between the depth-integrated ice and ocean hydrostatic pressures is $\rho g{H}^2/2-{\rho }_{w}g{b}^2/2=\rho g^{\prime }{H}^2/2$, where $g^{\prime }=g(1-\rho /{\rho }_w)$. Shelf force balance integrates $4\mu H{u}_x-\rho g^{\prime }{H}^2/2$ to a constant; absence of a buttressing force makes that constant zero. The shallow-shelf approximation therefore gives the displayed boundary condition, or $u_x=g^{\prime }H/(8\nu )$ with kinematic viscosity $\nu =\mu /\rho$.

For the constant gas-radiation ratio vertical disc structure, the integrated alpha prescription $\bar{\nu }\Sigma =\int \mu dz=(\alpha /\Omega )\int pdz$ gives the displayed mean kinematic viscosity and $\lambda /(1+\lambda )=\alpha H\Omega \kappa \Sigma /(8c)$. Uniform dynamic viscosity and nonuniform pressure are compatible because this prescription is an integrated equality. Imposing the corresponding pointwise equality would not preserve the exact vertical solution.

When kinematic viscosity and mass diffusivity are equal to $ϵ$, a monochromatic internal-wave beam with wavenumber $k$ and ray angle $\theta$ has leading stream-function amplitudeThe cubic dependence on $k$ makes slope-focused, short internal waves dissipate especially rapidly.

Momentum diffuses across a distance $L$ in the characteristic time $L^2/\nu$, where $\nu$ is the kinematic viscosity.