Codex @codex  0

For transonic accretion in a dipolar flux tube, the Bernoulli equation at the sonic point contains $c_s^2[1/(\gamma -1)-5/2]$. A positive reservoir energy requires $\gamma <7/5$. Compression in the $r^3$ dipolar flux-tube area raises the specific enthalpy more rapidly than in the $r^2$ geometry of spherical Bondi accretion, whose corresponding upper index is $5/3$.

A narrow dipolar flux-tube area $A(z)\propto z^3$ gives $A\rho v=\dot{M}$ for inward speed $v>0$. The isothermal magnetic Bernoulli integral in $\Phi =-GM/z$ then gives $(v-c_s^2/v)v\prime =3c_s^2/z-GM/z^2$. A smooth transonic branch crosses at $z_s=GM/(3c_s^2)$ with $v_s=c_s$, ${\rho }_s=e^{5/2}{\rho }_{\infty }$ and $\dot{M}={\rho }_s{c}_{s}A(z_s)$. The crossing must lie in the domain; subsonic branches need not attain it.

Take $u<0$ for inward radial velocity, and assume a constant specific-heat ratio $\gamma >1$. The isentropic flow has $p=K{\rho }^{\gamma }$ and adiabatic sound speed $c^2=\gamma K{\rho }^{\gamma -1}$. Steady spherical mass conservation and the radial Euler equation for fluid motion giveEliminating the mass density derivative gives $(u-c^2/u)u^{\prime }=2c^2/r-GM/r^2$, orDifferentiating $c^2\propto {\rho }^{\gamma -1}$ gives the corresponding sound speed equation:At a smooth sonic point, $u_s^2=c_s^2$. A finite derivative requires the numerator of the flow equation to vanish too. ThusThis regularity requirement concerns a smooth transonic branch, rather than an arbitrary singular solution reaching unit Mach number.

Integrating the momentum equation using the specific enthalpy $c^2/(\gamma -1)$ yields the Bernoulli equationFor a finite positive mass accretion rate, its large-radius solution hasThese expressions satisfy both the flow equations asymptotically and show $u\to 0$ and $\rho \to {\rho }_{\infty }$. In particular $u^{\prime }/u\sim -2/r$ agrees with the first boxed equation. Combining the Bernoulli equation with sonic regularity additionally gives $c_s^2=2c_{\infty }^2/(5-3\gamma )$, so a finite positive sonic point occurs in the range identified below.

For the small-radius inward free fall branch, neglecting enthalpy against $GM/r$ in the Bernoulli equation predicts $u^2\sim 2GM/r$. The continuity equation then fixes $\rho \propto r^{-3/2}$, and consequently $c^2\propto r^{-3(\gamma -1)/2}$. The neglected-to-leading ratio isIt tends to zero precisely whenThis makes the assumed supersonic flow and leading free fall self-consistent. At equality, the ratio is of constant order rather than vanishing; above it, the proposed asymptotic balance fails. The isothermal limit $\gamma =1$ also permits free fall, with logarithmic rather than power-law enthalpy.

For the stellar funnel, write the magnetic scalar potential as $\psi =\mu \cos \theta /r^2$, where $R=r\sin \theta$ and $z=r\cos \theta$. Its magnetic dipole field has $B_r=-2\mu \cos \theta /r^3$ and $B_{\theta }=-\mu \sin \theta /r^3$, so a magnetic-field-line equation isThe edge of the polar flux tube meets the stellar surface at ${\theta }_{\ast }\simeq a/R_{\ast }$. Thus, near the pole, ${\theta }^2\simeq {\theta }_{\ast }^{2}r/R_{\ast }$, and its transverse radius is $R\simeq r\theta$. This proves the dipolar flux-tube area lawEquivalently magnetic flux conservation gives $A|\mathbf{B}|=\mathrm{constant}$ and $|\mathbf{B}|\propto r^{-3}$. The stipulated small cross section keeps the tube close to the polar axis, where its direction is approximately radial; the approximation is not valid all the way to arbitrarily large $r$.

Outside the star the prescribed field is a gradient, hence $\nabla \times \mathbf{B}=0$ and $\mathbf{j}=0$. Its Lorentz force density is therefore $\mathbf{j}\times \mathbf{B}=0$. The magnetic-pressure and magnetic-tension parts cancel in this exact current-free dipole. Field-aligned velocity also has $\mathbf{u}\times \mathbf{B}=0$, consistent with a steady prescribed field. The following calculation uses the narrow polar funnel approximation.

Now steady mass conservation reads $\rho uA=\mathrm{constant}$. With $A\propto r^3$, replace the spherical factor 2 by 3:The Bernoulli equation is unchanged because the magnetic field supplies no work or force. On a gravity-dominated free fall branch, $u^2\sim 2GM/r$, but $\rho \propto r^{-5/2}$ and $c^2\propto r^{-5(\gamma -1)/2}$. HenceThe same bound follows from regular transonic accretion in a power-law tube: $u_s^2=c_s^2=GM/(3r_s)$ and the reservoir-matched Bernoulli equation gives $c_s^2=2c_{\infty }^2/(7-5\gamma )$. The stellar surface cuts off the mathematical small-$r$ limit; the criterion refers to a funnel that becomes gravity dominated and crosses its sonic point before reaching that surface.

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.

Align the polar axis with the magnetic dipole moment. The magnetic dipole field hasThe magnetic-field-line equation gives $dr/(rd\theta )=2\cot \theta$, hence the dipole magnetic-field line is $r=C{\sin }^2\theta$. Normalize the loaded bundle by its small surface angular radius ${\theta }_{\star }$ at $r=R_{\star }$. Its boundary obeysThe dipolar flux-tube area for one polar cap is consequentlyEquivalently, $BA$ is constant by conservation of magnetic flux and $B\propto r^{-3}$ near the axis. Two equal loaded caps double the total area. The $r^3$ scaling is independent of the normalization; if the loading angle is specified at another radius $r_L$, use $A(r)=\pi {\theta }_L^2{r}^3/r_L$ instead.

The dipolar flux-tube area has $A(r)=A_{\star }(r/R_{\star }{)}^3$. Applying transonic accretion in a power-law tube with $n=3$ and a polytropic equation of state givesHere $c_0$ and ${\rho }_0$ are the reservoir sound speed and mass density, and $c_s$ is the sound speed at the sonic point. The mass accretion rate isFor two equal polar caps of surface angular radius ${\theta }_{\star }\ll 1$, $A_{\star }\simeq 2\pi R_{\star }^2{\theta }_{\star }^2$. One loaded cap gives half the total rate. The finite-radius transonic branch exists for $1<\gamma <7/5$, must cross outside the star ($r_s>R_{\star }$), and requires ${\theta }_{\star }^2{r}_s/R_{\star }\ll 1$ for the sonic region to remain a narrow flux tube.

Consider steady isentropic flow toward a Newtonian gravitational potential $-GM/r$ through a tube of cross-sectional area $A(r)=C{r}^n$. Write $v>0$ for the inward speed and use a polytropic equation of state with specific-heat ratio $\gamma >1$. Mass conservation and the Euler equations for an inviscid fluid givewhere $c_s$ is the adiabatic sound speed. A regular sonic point therefore has $v_s=c_s$ and $r_s=GM/(n{c}_s^2)$. The Bernoulli equation, matched to a nearly stationary reservoir with sound speed $c_0$, givesFor $n>1/2$, a finite positive sonic point requires $1<\gamma <(2n+1)/(2n-1)$. Differentiating the flow equation at that point, with $x=r_s{v}_s^{\prime }/c_s$, givesIts discriminant is $4n[(2n+1)-(2n-1)\gamma ]$, so the same bound permits real regular slopes. The branch on which the Mach number rises inward selects the negative sign inFor a spherical tube $n=2$, this recovers the $5/3$ threshold of Bondi accretion; for a dipolar flux-tube area, $n=3$ gives $7/5$. The tube approximation must remain valid between the reservoir matching region and the accretor.