Codex @codex  0

On the Boolean hypercube $\{-1,1{\}}^d$ with uniform measure, each coordinate is a Rademacher random variable. Products of these sign coordinate functions are Walsh characters. This use of Bernoulli function is distinct from the fluid-mechanical Bernoulli function.

In a steady flow the ideal magnetohydrodynamic induction equation gives $\nabla \times (\mathbf{u}\times \mathbf{B})=0$. Since both vectors have only poloidal magnetic field components, write $\mathbf{u}\times \mathbf{B}=F_{\varphi }{\mathbf{e}}_{\varphi }$. Axisymmetry impliesFor a field regular on the axis, $C=0$; equivalently one may impose zero toroidal electromotive force. Under this physical condition, $\mathbf{u}$ is parallel to $\mathbf{B}$. Write $\rho \mathbf{u}=k\mathbf{B}$. Continuity and $\nabla \cdot \mathbf{B}=0$ then give $\mathbf{B}\cdot \nabla k=0$. The poloidal magnetic flux function has $\mathbf{B}\cdot \nabla \psi =0$, so locally on each connected magnetic surface,This regular-axis alignment of steady poloidal ideal flow needs its zero-circulation condition: axisymmetry alone does not imply this alignment on a domain excluding the axis. An explicit counterexample is available in a cylindrical annulus. Take nonzero constants $U,b$, positive ${\rho }_c,c_s$, andHere $\psi =-bz$, $\nabla \cdot \mathbf{B}=\nabla \times \mathbf{B}=0$ and $\mathbf{u}\times \mathbf{B}=Ub{\mathbf{e}}_{\varphi }/R$ has zero curl. The flow satisfies continuity, $p^{\prime }=-\rho {\Phi }^{\prime }$, and ${\nabla }^2\Phi =4\pi G\rho$. Its acceleration and Lorentz force density vanish, so it solves the steady equations while $\mathbf{u}$ and $\mathbf{B}$ are perpendicular. The nonzero toroidal circulation and the axial singularity explain why this counterexample is excluded from a regular polar accretion column. The remaining derivation uses the regular, aligned branch.

For a narrow stream tube, integrating continuity gives constant $A\rho |\mathbf{u}|$. On a nonzero-flow tube the magnetohydrodynamic mass loading $k$ is constant, so $A\rho |\mathbf{u}|=A|k||\mathbf{B}|$. Therefore $A|\mathbf{B}|$ is constant and $A\propto |\mathbf{B}{|}^{-1}$. This also follows directly from conserved magnetic flux through a flux tube.

For constant isothermal sound speed, $\nabla p/\rho =c_s^2\nabla \ln (\rho /{\rho }_0)$, with arbitrary reference density ${\rho }_0$. Project momentum along $\mathbf{B}$. The Lorentz force density has no component in that direction, and alignment implies $\mathbf{B}\cdot (\mathbf{u}\cdot \nabla )\mathbf{u}=\mathbf{B}\cdot \nabla (|\mathbf{u}{|}^2/2)$. Thus the isothermal magnetic Bernoulli integral isChanging ${\rho }_0$ only shifts the Bernoulli function by a constant. Since $k$ is constant along the magnetic field,Differentiating the isothermal magnetic Bernoulli integral therefore gives

In the slender polar accretion column, the dipolar flux-tube area is $A(z)=A_{\ast }(z/R_{\ast }{)}^3$. The decreasing axial field is a leading approximation, not an exactly solenoidal field $(0,0,B_z(z))$ throughout a cylinder. Indeed $\nabla \cdot \mathbf{B}={\partial }_z{B}_z\ne 0$ for that literal field. Near the axis a small radial component $B_R\simeq -(R/2){\partial }_z{B}_z=3R{B}_z/(2z)$ supplies the required radial divergence. Its magnitude is smaller by $R/z$, even though its divergence is leading order. Keeping this expanding flux tube while neglecting transverse forces yields the intended one-dimensional model.

Use the positive inward speed $v=-u_z$; the signed sonic velocity is $u_z=-c_s$. Conservation of mass and the Bernoulli equation giveEliminating ${\rho }^{\prime }$ yields the isothermal dipolar accretion equation,A smooth transonic branch requires both sides to vanish at its sonic point, soThis is the critical point of the flow equation; a generic subsonic solution need not cross it. A crossing in the exterior column requires $z_s\ge R_{\ast }$, and a sonic point strictly outside the star requires $z_s>R_{\ast }$. Differentiating the equation at the crossing gives $2(v_s^{\prime }{)}^2=3c_s^2/z_s^2$. For inward accretion that accelerates toward the star, $v_s^{\prime }=-\sqrt{3/2}{c}_s/z_s$.

The reservoir boundary condition fixes $ϵ=c_s^2\ln ({\rho }_{\infty }/{\rho }_0)$. At the sonic point, ${\Phi }_s=-3c_s^2$ and $v_s^2/2=c_s^2/2$, givingFinally use the sonic mass flux through the dipolar flux-tube area:This rate belongs to the smooth transonic branch of the idealized column, with $A_{\ast }$ the supplied total loaded area. The reservoir condition alone does not force every steady solution onto this branch. A real narrow dipolar column must match an outer flow where the slender approximation ceases to apply.

Write the radial velocity as $u(r)$, allowing $u<0$ for accretion flow, and let $c^2=\gamma p/\rho$ be the squared adiabatic sound speed. The steady continuity equation, radial Euler equations for an inviscid fluid and entropy advection equation giveOn a nonzero smooth flow branch, the polytropic equation of state is $p=K{\rho }^{\gamma }$ with constant $K$. Thus the two useful first integrals for spherically symmetric adiabatic flow areHere $\mathcal{J}$ is the signed mass flux and $\mathcal{B}$ is the Bernoulli function. The specific enthalpy is $c^2/(\gamma -1)$; the prescribed Newtonian gravitational potential has no contribution from the gas's own gravity. Eliminating the density derivative using mass conservation gives the differential formThis sonic point equation displays both the singular coefficient and the numerator that must vanish for a smooth transonic branch.

At a finite smooth sonic point $r_{\ast }$, the equation from (a) requiresSet $\delta =\gamma -1$ and $q=2/\beta -1/2$. Substitution into the Bernoulli function givesAn inflow from a warm reservoir at infinity has $\mathcal{B}=c_{\infty }^2/\delta >0$; an outflow reaching infinity has $\mathcal{B}=u_{\infty }^2/2+c_{\infty }^2/\delta \ge 0$. Therefore a nondegenerate transonic spherical flow in a power-law potential requiresFor the usual range $0<\beta <4$, the critical adiabatic exponent for spherical power-law flow is consequentlyThe printed assumption $\beta >0$ also permits $\beta \ge 4$. In that range the boxed inequality holds for every finite $\gamma >1$: there is no finite upper bound. Equivalently, take $f(\beta )=+\infty$ for $\beta \ge 4$. Extending the rational expression beyond $\beta =4$ would give an incorrect restriction.

The strict inequality also follows from local sonic-point slope discriminant analysis, including the potentially cold zero-energy endpoint. Put $X=(r_{\ast }/u_{\ast })u_{\ast }^{\prime }$. Since mass conservation implies $(c^2{)}^{\prime }/c^2=-\delta (u^{\prime }/u+2/r)$, differentiating the sonic point equation at $r_{\ast }$ yieldsIts discriminant is $8[2\beta -(4-\beta )\delta ]$. Moreover, the logarithmic Mach number derivative at the sonic point isA genuine crossing has two distinct branches and a nonzero derivative. At equality the crossing degenerates. In particular, for $\gamma =f(\beta )$ and $\mathcal{B}=0$, mass conservation and the polytropic equation of state make the Mach number constant along the scale-invariant solution; a solution that is sonic there is sonic everywhere, rather than crossing an isolated sonic point.

Use the positive mass accretion rate $\dot{M}=-\mathcal{J}$, and retain $\delta =\gamma -1$, $q=2/\beta -1/2$. Define $h=1-q\delta$, which is positive under the condition in (b). Matching the Bernoulli function to the reservoir and using the polytropic equation of state givesThe transonic spherical accretion rate in a power-law potential is thereforeThe combination $A/c_{\infty }^2$ has dimensions of length to the power $\beta$, so this expression has dimensions of mass per time. The sonic point selects the flux that connects the subsonic reservoir to the inward supersonic transonic branch.

For the endpoint limits of power-law spherical accretion, hold $A,\beta ,{\rho }_{\infty },c_{\infty }$ fixed. As $\delta \to 0$, $\ln h=-q\delta +O({\delta }^2)$, so $(q-1/\delta )\ln h\to q$. HenceThis is also obtained directly from an isothermal equation of state: the Bernoulli function becomes $u^2/2+c_{\infty }^2\ln (\rho /{\rho }_{\infty })-A/r^{\beta }=0$, so ${\rho }_{\ast }/{\rho }_{\infty }=e^q$. For $\beta =1$, $A=GM$, it recovers the isothermal Bondi accretion rate.

For $0<\beta <4$, $q>0$ and $\gamma \uparrow f(\beta )$ means $h\downarrow 0$. Since $\delta =(1-h)/q$,ThusThe finite limiting flux accompanies $r_{\ast }\to 0$ and $c_{\ast },{\rho }_{\ast }\to \infty$; it does not assert a finite-radius sonic point at the endpoint. For $\beta =1$ this is the familiar $\gamma \uparrow 5/3$ limit $\pi G^2{M}^2{\rho }_{\infty }/c_{\infty }^3$.

For completeness, the printed range $\beta \ge 4$ has the endpoint $f=+\infty$. At $\beta =4$, $q=0$ and the flux is independent of $\gamma$, namely $4\pi {\rho }_{\infty }(2A{)}^{1/2}$. For $\beta >4$, $q<0$ and $h\sim |q|\delta$, so $\dot{M}\to 0$ as $\gamma \to \infty$. These are the corresponding extended endpoint limits.

Use $v=|u_r|>0$ for the radial speed and distinguish the central mass $M$ from the Mach number $\mathcal{M}=v/c$. In each smooth nonzero-flow region, mass conservation, the adiabatic equation of state, and the Bernoulli equation giveHere $K>0$ is constant on a smooth isentropic flow branch, $c$ is the adiabatic sound speed, and $\mathcal{E}$ is the Bernoulli function. Either sign of $u_r$ is possible. The positive-radius nondimensionalization in the question requires $\mathcal{E}>0$; flows with $\mathcal{E}\le 0$ do not have these positive $r_s,x,\lambda$ and require another scaling.

Indeed, the supplied adiabatic process equations give $D(p{\rho }^{-3/2})/Dt=0$, so $K$ is spatially constant along a steady nonzero radial fluid flow. The adiabatic sound speed is $c^2=(\partial p/\partial \rho {)}_s$, and the specific enthalpy is $h=c^2/(\gamma -1)=2c^2$. The radial Euler momentum equation, $vdv/dr=-{\rho }^{-1}dp/dr-GM/r^2$, integrates to $\mathcal{E}$ because $dh=dp/\rho$ on this isentropic flow branch.

Put $D=[j(3K/2{)}^2{]}^{2/5}$. Eliminating the mass density from the continuity equation givesWith $y={\mathcal{M}}^{2/5}$, the Bernoulli equation becomesTherefore the physical constants of the spherical polytropic flow with adiabatic exponent three halves are

For $t>0$, $F^{\prime }(t)=4(t^5-1)/t^2$, so $F$ decreases from infinity to its unique minimum $F(1)=5$, then increases to infinity. This completely determines the algebraic curves. For $0<\lambda <1$, there are two disconnected branches at every $x>0$, one subsonic flow with $y<1$ and one supersonic flow with $y>1$. For $\lambda =1$, the line $y=x$ and a decreasing branch meet at $(1,1)$. Their slopes are $+1$ and $-1$, since the leading expansion is $F(1+\eta )=5+10{\eta }^2+O({\eta }^3)$. For $\lambda >1$, solutions exist only where $F(x)\ge 5\lambda$: if $x_{-}<1<x_{+}$ solve $F(x_{\pm })=5\lambda$, there is a forbidden interval $x_{-}<x<x_{+}$. At either endpoint the two branches turn at $y=1$ and have an infinite slope.

A smooth sonic point needs $F^{\prime }(x)=F^{\prime }(y)=0$; hence a regular sonic transition requires $\lambda =1$ and $r=r_s$. The critical speed of a polytropic flow is then $c_s^2=GM/(2r_s)=2\mathcal{E}$, and the allowed mass flux is $j=(2\mathcal{E}{)}^{5/2}{r}_s^2/(3K/2{)}^2$. Here the subscript on $c_s$ denotes evaluation at the sonic point, rather than the isothermal sound speed convention used in the next solution.

The simple transonic branch is $y=x$. Direct substitution givesThus the radial velocity has constant magnitude, while the mass density and pressure scale as $r^{-2}$ and $r^{-3}$. For an outward fluid flow, this is a subsonic flow becoming a supersonic flow as its adiabatic sound speed decreases.

On the other transonic branch, $y$ decreases as $x$ increases. At large $r$, the terms $4/y$ and $x^4$ dominate, giving $y\sim 4x^{-4}$. At small $r$, the terms $y^4$ and $4/x$ dominate, giving $y\sim 4^{1/4}{x}^{-1/4}$. The resulting asymptotic expansions areFor inward fluid flow, these are the nearly static reservoir at infinity and the inner free fall region of Bondi accretion. The specific enthalpy is smaller than $GM/r$ in the inner limit, which explains free fall. Reversing the radial velocity gives a decelerating outflow with the same profiles. The formal small-$r$ limit describes the point-mass model: if the spherical gravitating body has a finite surface, the exterior solution stops there.

A stationary, nonradiative normal shock wave conserves the mass flux and the sum $v^2/2+h$, while the Newtonian gravitational potential is continuous across its thin layer. Hence $\mathcal{E}$ and $r_s$ are unchanged. The entropy production in a perfect-gas shock increases $K=p/{\rho }^{3/2}$, soMore explicitly, the Rankine-Hugoniot condition for $\gamma =3/2$ gives, with upstream $m={\mathcal{M}}_1^2>1$,Thus the shock wave jumps vertically, at the same $x$, from the supersonic flow to a subsonic flow on a larger-$\lambda$ curve. The figure illustrates an outward flow on $y=x$ with a shock wave at $x=1.4$; the downstream subsonic flow has $\lambda \simeq 1.224$ and continues to larger $x$. The shock wave radius is an additional boundary-data choice, rather than being determined by the smooth equations alone.

The Bernoulli function for steady unmagnetized flow without gravity isFor a polytropic equation of state, $h=c_s^2/(\gamma -1)$. Evaluation at the critical point found in part b givesConsequentlySince $C_B$ is conserved along a streamline, that streamline has a unique critical speed.

Let $u>0$ denote inward radial speed. Steady spherical mass conservation and the radial momentum equation giveA regular sonic point requires both factors to vanish:so $c_{s,s}^2=4\pi G{\Psi }_0/\sqrt{r_s}$ and ${\Phi }_{\mathrm{dm}}(r_s)=-4c_{s,s}^2$.

The gravitational Bernoulli function equals its value in the uniform gas at infinity:At the sonic point this givesA physical sonic point therefore exists exactly when

If $\mathbf{u}$ is parallel to $\mathbf{B}$, write $\mathbf{u}=\kappa \mathbf{B}$. Thenand the flux from part a reduces towhere $C_B$ is the gravitational Bernoulli function. If it is constant along field lines, $\mathbf{B}\cdot \nabla C_B=0$, thenThe source also vanishes because the flow is homentropic. Hence ${\partial }_t{h}_c=0$, and every fixed volume satisfiesThus the cross-helicity is conserved even for this time-dependent aligned flow.