Codex @codex  0

Put $\mu =GM$ and use energy and angular momentum per unit particle mass. The Kepler orbit equation is $\ddot{\mathbf{r}}=-\mu \mathbf{r}/r^3$. HenceThus $E$ and $\mathbf{L}=\mathbf{r}\times \dot{\mathbf{r}}$ are integrals of motion. Splitting the velocity into radial and transverse parts, $|\dot{\mathbf{r}}{|}^2={\dot{r}}^2+L^2/r^2$, gives the Kepler radial energy equationHere $\dot{r}$ denotes the derivative of the scalar radius, unlike $\dot{\mathbf{r}}$.

For the Laplace-Runge-Lenz vector, differentiate directly:The vector triple-product identity yieldsThe terms cancel, proving every component of $\mathbf{A}$ is conserved.

Write $q=\bar{r}(r)$ and $\tau =\bar{t}$ for a monotone radial orbit time transformation. The chain rule givesImpose $q^{\prime }=(r/q)\dot{r}$, equivalentlyMultiplying the radial energy equation by $r^2/q^2$ then givesThis is a transformation along each trajectory; the new clock is obtained by integrating its radius-dependent rate.

For the Kepler–harmonic radial duality, take the old inverse-radius term as the new constant energy, $\bar{E}=\mu r/q^2$. On a positive-radius branch $\bar{E}>0$, soThe remaining term is minus the new Newtonian gravitational potential:For $E<0$ this is an attractive isotropic harmonic oscillator with ${\bar{\omega }}^2=-2E(\bar{E}/\mu {)}^2$. For $E>0$ it is an inverted oscillator, and for $E=0$ the Newtonian gravitational potential vanishes. Calling it a confining oscillator therefore requires a bound original Kepler orbit.

For the angular reconstruction of a radial orbit transformation, choose the new angle by $d\bar{\theta }/d\tau =L/q^2$. The same $L$ then supplies its centrifugal term. From $d\theta /dt=L/r^2$, the transformation is $d\bar{\theta }/d\theta =(r/q)(dq/dr)$; in the harmonic case this gives $\bar{\theta }=\theta /2$ up to a constant. Keeping the old angle unchanged would not in general preserve the angular momentum used in the displayed transformed radial equation.

For the second transformation use the explicit displayed split of the energy terms. It assigns weight $\lambda$ to $E$ and $1-\lambda$ to $\mu /r$; the surrounding printed prose interchanges these weights. Define the transformed quantities consistently byTheir difference is exactly $(r^2/q^2)(E+\mu /r)$, so they give the same transformed radial energy equation.

Let $a=(1-\lambda )\mu$, $h=\lambda E$, and $\epsilon =\bar{E}$. Constancy of the new energy imposesFor the standard regular isochrone branch take $h>0$, $\epsilon >0$, and $a\ge 0$. The positive-radius root isFor $q>0$ it is monotone. The corresponding clock is also regular there:With $\alpha =\sqrt{\epsilon /h}$ and $s=\sqrt{q^2+b^2}$, rationalization gives $r/q^2=\alpha /(s+b)$. Eliminate $r^2/q^2$ using $h{r}^2/q^2=\epsilon -ar/q^2$:where one possible identification isThis is the Kepler–isochrone radial transformation to the spherical isochrone model, up to the additive energy constant $k$.

The attractive branch of the Kepler–isochrone transformation requires a suitable parameter domain. The displayed square root needs $\lambda E\bar{E}>0$, and the positive-root branch used above needs $\lambda E>0$, $\bar{E}>0$. An attractive regular isochrone also needs $b\ge 0$ and $\bar{M}>0$. For bound original orbits a convenient domain is $E<0$, $\lambda <0$, $\bar{E}>0$: it satisfies all these requirements and gives $\bar{E}-k=\bar{E}/\lambda <0$. For $E>0$ one may instead take $1/2<\lambda \le 1$, with $\bar{E}>0$; these map to unbound attractive isochrone orbits. Other parameter choices can produce a repulsive Newtonian gravitational potential or a different branch, so no unrestricted real-parameter statement is justified. The singular algebraic case $\lambda =0$ is the harmonic transformation already handled separately. The angular reconstruction $d\bar{\theta }/d\tau =L/q^2$ likewise completes the isochrone radial solution into a central-force orbit.

Let $A_s$ denote the surface area of an n-sphere for the unit sphere $S^s$. In this notation $s$ is the intrinsic dimension of the sphere. Cartesian factorization of the Gaussian integral gives $I_N={\pi }^{N/2}$. In hyperspherical coordinates,The substitution $t={\rho }^2$ supplies the defining Gamma function integral. Scaling the sphere radius multiplies its surface volume by $r^{N-1}$, henceThis is the volume of the sphere itself, not of its enclosed ball.

Take $D$ to denote spacetime dimension and $d=D-1$ spatial dimensions. In uncompactified flat space, Gauss's law for a charge $q$ and permittivity ${\epsilon }_d$ gives $E_r{A}_{d-1}{r}^{d-1}=q/{\epsilon }_d$. Since $E_r=-d\phi /dr$, the electric potential for $d>2$, normalized to vanish at infinity, isThus the four-dimensional electric potential falls as $1/r$, and the five-dimensional one as $1/r^2$. If a convention counts only spatial dimensions as $D$, replace $d$ by $D$ in the first formula. For $d=2$, the integral instead gives $\phi =-q\log (r/r_0)/(2\pi {\epsilon }_2)$; for $d=1$, it gives $\phi =-q|x|/(2{\epsilon }_1)$ up to a constant. Neither case has a finite zero at infinity.

For a Newtonian gravitational potential ${\Phi }_g$, define the flux-normalized coupling ${\kappa }_d$ by ${\nabla }^2{\Phi }_g={\kappa }_{d}m{\delta }^{(d)}(\mathbf{x})$ and acceleration $-\nabla {\Phi }_g$. Exactly the same flux integration gives the attractive potentialIn $d=3$, ${\kappa }_3=4\pi G_4$ reproduces $-G_{4}m/r$. If $G_D$ is defined by the $D$-dimensional Einstein-Hilbert action $S=(16\pi G_D{)}^{-1}\int \sqrt{-g}R$, the weak static Einstein field equations give ${\kappa }_d=8\pi G_D(D-3)/(D-2)$ for $D>3$. ThereforeIndeed, $R_{00}\simeq {\nabla }^2{\Phi }_g$ and $T_{00}-g_{00}T/(D-2)=(D-3)\rho /(D-2)$ in signature $(-+\cdots +)$. A Poisson-defined Newtonian theory in two spatial dimensions has a logarithmic potential, but pure three-dimensional Einstein gravity has no analogous local Newtonian point-mass force; the relativistic normalization must not be extrapolated to that exceptional dimension.

For an unwarped product with a fixed extra-dimensional metric, normalize the Einstein-Hilbert action by $\frac{1}{2}{M}_{\ast }^{n+2}\int d^{4+n}x\sqrt{-g}R$. Integrating the part containing the four-dimensional curvature over the internal space gives $\frac{1}{2}{M}_{\ast }^{n+2}\mathrm{vol}(S_R^n)\int d^{4}x\sqrt{-g_4}{R}_4$. Thus, using the reduced Planck mass,This is a volume relation for compactification; a curved sphere additionally requires a mechanism supporting and stabilizing its background. It is not by itself a proof that a sphere times flat spacetime solves a vacuum gravitational theory.

With $M_{\ast }=10^3\mathrm{GeV}$, $M_{\mathrm{Pl}}=2.435\times 10^{18}\mathrm{GeV}$ and $1{\mathrm{GeV}}^{-1}=1.97327\times 10^{-16}\mathrm{m}$, the spherical-volume factors $A_1=2\pi$, $A_2=4\pi$, $A_6=16{\pi }^3/15$ giveThe first estimate is astronomical; the second is roughly a tenth of a millimetre. If both Planck scales are instead defined by $S=M_D^{D-2}\int \sqrt{-g}R/(16\pi )$, the same algebraic volume relation uses the unreduced $M_{\mathrm{planck}}=1.221\times 10^{19}\mathrm{GeV}$. Holding that convention's $M_{\ast }=1\mathrm{TeV}$ gives radii larger by $(8\pi {)}^{1/n}$: approximately $4.69\times 10^{12}\mathrm{m}$, $6.80\times 10^{-4}\mathrm{m}$, $2.54\times 10^{-14}\mathrm{m}$. These are different conventions for what is held fixed at one TeV, not different scaling laws.

Finally, compactify the fifth coordinate on a circle of circumference $L=2\pi R$. For a source and observer at the same compact coordinate, the method of images turns a five-dimensional static potential $C/{\rho }^2$ intoOne can derive this image sum without assuming a summation formula. The Fourier transform of $1/(r^2+y^2)$ is $\pi e^{-r|k|}/r$; Poisson summation then gives $\pi [1+2{\sum }_{\ell \ge 1}{e}^{-2\pi \ell r/L}]/(Lr)$, which is the displayed hyperbolic cotangent. ThereforeThe short-distance limit is $C/r^2$; the long-distance limit is the four-dimensional $1/r$ potential. These exponential corrections are the nonzero Kaluza-Klein modes, with masses $2\pi |\ell |/L$. For electrostatics $C=q/(4{\pi }^2{\epsilon }_4)$, so the effective permittivity is ${\epsilon }_3=L{\epsilon }_4$. For a Poisson-normalized gravitational potential, $C=-{\kappa }_{4}m/(4{\pi }^2)$ gives ${\kappa }_3^{\mathrm{eff}}={\kappa }_4/L$. A massless scalar associated with the circle radius also contributes in an unstabilized gravitational compactification; recovering pure four-dimensional Einstein gravity with $G_4=G_5/L$ requires its stabilization or removal. The $1/r$ crossover itself is independent of this tensor-versus-scalar normalization issue.

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.

In magnetostatics, $-\nabla p-\rho \nabla \Phi +{\mathbf{f}}_L=0$. The plane-parallel Newtonian gravitational potential and all thermodynamic fields are independent of $y$, so the $y$ component of the Lorentz force must vanish:Away from a null of $\nabla \psi$, this says that $B_y$ is constant along each connected contour of the Cartesian magnetic flux function. Thus, on a regular flux region,This statement is local; disconnected contours with the same numerical flux label need not share one global function without the usual flux-tube connectivity assumption.

Since $\nabla B_y=b^{\prime }(\psi )\nabla \psi$, the transverse magnetostatic equilibrium equation becomesThis is a Cartesian magnetostatic flux-function equilibrium. Projecting it along a transverse magnetic field line gives hydrostatic balance along that line; projecting across the line balances the gas pressure and weight against magnetic pressure and magnetic tension.

Use the fluid displacement $\mathbf{\xi }$, with velocity perturbation $\delta \mathbf{u}={\partial }_t\mathbf{\xi }$. Denote Eulerian perturbations by $\delta$ and the corresponding Lagrangian pressure perturbation by ${\Delta }_{L}p=\delta p+\mathbf{\xi }\cdot \nabla p$. For adiabatic perturbations, the linearized mass conservation, momentum, adiabatic equation of state and Poisson equation areThese self-gravitating adiabatic displacement equations retain the perturbation of the star's own Newtonian gravitational potential. In the uniform-density interior, $\delta \rho =-\rho \nabla \cdot \mathbf{\xi }$, $\nabla \Phi ={\omega }_d^2\mathbf{r}$, and $\nabla p=-\rho {\omega }_d^2\mathbf{r}$. The dynamical frequency of a uniform-density star ${\omega }_d=(GM/R^3{)}^{1/2}$ sets the natural timescale. For a normal mode with time factor $e^{-i\omega t}$, replace ${\partial }_t^2$ by $-{\omega }^2$. The requested interior analysis needs no surface or exterior matching conditions.

For a spherical system, use the spherical shell theorem with shell mass $dm=4\pi \rho (s)s^{2}ds$. Inner shells contribute $-Gdm/r$ to the gravitational potential, while an outer shell contributes its constant interior potential $-Gdm/s$. With $\Phi (\infty )=0$, assuming the required integrals converge,On differentiation, the two terms containing $\rho (r)r$ cancel. Hence ${\Phi }^{\prime }=GM(r)/r^2$, where $M(r)=4\pi {\int }_0^r\rho (s)s^{2}ds$. Radial balance for a circular orbit gives $v_c^2=r{\Phi }^{\prime }=GM(r)/r$.

For a razor-thin axisymmetric astrophysical disk, the element of mass is $dm=\Sigma (s)sdsd{\varphi }^{\prime }$, where $\Sigma$ is the surface density. Superposing Newtonian gravitational potentials therefore givesAxisymmetry permits $\varphi =0$. Unlike the spherical case, an exterior annulus generally exerts a radial force: enclosed mass does not determine a disc rotation curve.

For the Legendre expansion of thin-disk gravity, split the radial integral at $s=R$. In the inner part the kernel expands in $s/R$, while in the outer part it expands in $R/s$. The angular integrals of odd Legendre polynomials vanish, and those of even degree $n$ equal $2{\alpha }_n$, whereIntroduce $I_n(R)={\int }_0^R\Sigma (s)s^{n+1}ds$ and $J_n(R)={\int }_R^{\infty }\Sigma (s)s^{-n}ds$. The gravitational potential becomesWhen differentiating each bracket, the moving-limit terms cancel: $R^{-(n+1)}{I}_n^{\prime }=\Sigma (R)$ and $R^n{J}_n^{\prime }=-\Sigma (R)$. Since $v_c^2=R{\Phi }^{\prime }$, this givesThe zeroth term has ${\alpha }_0=\pi$ and $M(R)=2\pi I_0$. Separating it provesThe inner correction is inward, while the exterior correction is outward. For a smooth surface density, the original potential singularity at coincident points is integrable. Its radial force is understood through a symmetric Cauchy principal value or a vanishing-thickness regularization; the paired interior and exterior terms above retain the cancellation at $s=R$. This avoids treating the two singular local force contributions separately.

For an exponential galactic disk, write $\Sigma (s)={\Sigma }_0{e}^{-s/R_d}$, so $M_{\infty }=2\pi {\Sigma }_0{R}_d^2$ and $M(R)=M_{\infty }[1-(1+R/R_d)e^{-R/R_d}]$. At large $R$, the missing mass and exterior-ring terms are exponentially small. The leading nonspherical interior term is $k=1$, with ${\alpha }_2=\pi /4$ and $I_2(\infty )=6{\Sigma }_0{R}_d^4$. Consequently the exponential-disk Keplerian asymptotic isThe positive leading correction shows that the rotation curve approaches the Keplerian limit from above. The finite-order large-radius expansion is sufficient here; an infinite moment expansion need not converge for a disk extending to infinity.

A useful special example is the Mestel disk, with surface density $\Sigma (s)=C/s$ for $s>0$, $C>0$. It has $M(R)=2\pi CR$. For every positive even $n$,All the nonspherical corrections cancel, givingThus the Mestel disk has a perfectly flat galaxy rotation curve.

This example has infinite total mass and a singular central surface density. The absolute gravitational potential cannot be set to zero at infinity, but the radial force exists as a limit of disks with increasing outer cutoff. Potential differences are $\Phi (R)-\Phi (R_0)=v_M^2\log (R/R_0)$. The earlier potential integral is therefore interpreted up to a radius-independent divergent constant for this example; the force calculation remains valid. Truncating the Mestel disk gives a more physical finite system but changes the exact flat curve near its edges.

Integrating the standard thin-disk dissipation flux over annular area, with both faces already included in $\mathcal{H}$, givesThus the disk luminosity for $r_{\mathrm{out}}\gg r_{\ast }$ isThe Newtonian gravitational potential decreases by approximately $G{M}_{\ast }/r_{\ast }$ per unit mass from a distant outer edge to the surface, giving a potential-energy release rate $G{M}_{\ast }\dot{M}/r_{\ast }$. The standard thin-disk luminosity is half of this.

The missing half remains as kinetic energy of nearly circular orbital motion: at the inner edge $v_K^2=G{M}_{\ast }/r_{\ast }$, so the specific orbital kinetic energy is $G{M}_{\ast }/(2r_{\ast })$. Equivalently, circular motion has total specific mechanical energy $-G{M}_{\ast }/(2r_{\ast })$. Matter joining a slowly rotating star must shed this orbital motion in an accretion-disk boundary layer, producing approximately another $G{M}_{\ast }\dot{M}/(2r_{\ast })$ of luminosity. A rotating star can retain some energy in spin, so equal disk and accretion-disk boundary layer luminosities assume slow stellar rotation. For a central black hole, there is no material surface, and energy can instead be carried inward.

For nonzero real wavenumber $k$, seek a Fourier mode of the Newtonian gravitational potential in the form ${\Phi }^{\prime }=f(z)e^{ikx}$. Away from the razor-thin sheet, the Poisson equation for Newtonian gravity becomes $f^{\prime \prime }-k^{2}f=0$. Requiring the perturbation to decay on both sides and remain continuous gives $f=A{e}^{-|k||z|}$.

Integrating the Poisson equation for Newtonian gravity through $z=0$ fixes the derivative jump:Thus the perturbing gravitational potential isAt the midplane this reduces to the razor-thin disk Poisson kernel. The observable real disturbance is the real part. Writing the perturbation amplitude as ${\Sigma }_a$ instead of ${\Sigma }_0$ gives the same formula with ${\Sigma }_a$; linearity makes it independent of the background density. The $k=0$ disturbance is a uniformly changed sheet and has potential proportional to $|z|$, rather than a decaying nonzero-wavenumber solution.

Let $e$ be specific internal energy, $Q$ the heat supplied per unit volume per unit time, and $E=\rho (e+|\mathbf{u}{|}^2/2)$. For an inviscid perfect gas in a prescribed Newtonian gravitational potential $\Phi$, the fluid total-energy equation isThere is no gravitational term if no external body is present. For a time-independent $\Phi$, adding potential energy gives the equivalent conservative form ${\partial }_t(E+\rho \Phi )+\nabla \cdot [(E+p+\rho \Phi )\mathbf{u}]=Q$. Subtracting the kinetic-energy equation and using mass conservation givesFor a perfect gas with constant specific-heat ratio $\gamma$, $e=p/[(\gamma -1)\rho ]$. The isothermal equation of state therefore makes $e=c_s^2/(\gamma -1)$ constant. Hence the required heat supply and heat loss areCompression requires cooling; expansion requires heating. The isothermal sound speed $c_s$ differs from the adiabatic sound speed $\sqrt{\gamma }{c}_s$.

Define the isothermal Mach numbers by ${\mathcal{M}}_i=u_i/c_s$. Orient the normal to a stationary isothermal shock along the flow, so the mass flux $j={\rho }_1{u}_1={\rho }_2{u}_2>0$. Tangential velocity is continuous for this planar hydrodynamic shock. Normal momentum conservation givesFor a genuine discontinuity $u_1\ne u_2$, factorization gives $u_1{u}_2=c_s^2$. Thus the isothermal shock jump relations areAcross a thin shock, gravitational potential and tangential kinetic energy are unchanged. The specific enthalpy $h=\gamma c_s^2/(\gamma -1)$ is the same on both sides, so the energy removed per unit area per unit time isPositive net cooling requires $u_1>u_2$, equivalently ${\rho }_2>{\rho }_1$ and ${\mathcal{M}}_1>1>{\mathcal{M}}_2$. Cooling across an isothermal shock permits only compression shocks; the formal expansion discontinuity requires energy supplied by the surroundings. The continuous state $u_1=u_2$ is not a shock.

For steady radial flow, mass conservation gives $r^2\rho u_r=\text{constant}$. Differentiating it and eliminating the mass density gradient from radial momentum balance givesA regular sonic point has $|u_r|=c_s$, so the right-hand side must vanish there as well. ThereforeThis is the Bondi sonic point for accretion; the same critical radius appears in an isothermal Parker wind. A general steady solution need not pass through a sonic point, but a smooth transonic one must obey both critical conditions.

Put $x=r/r_s$ and $\mathcal{M}=|u_r|/c_s$ on a branch whose flow direction is fixed. The radial equation becomesIntegrating yields the integrated isothermal Bondi flow relation:For a transonic branch, $x=\mathcal{M}=1$ gives $C=-3/2$. Expansion about the critical point gives $(\mathcal{M}-1{)}^2=(x-1{)}^2+O((x-1{)}^3)$, with slope $-1$ for inflow supplied by gas at rest at infinity and slope $+1$ for a transonic outflow.

For Isothermal Bondi accretion, integrate radial momentum balance once more, using $u_r\to 0$ and $\rho \to {\rho }_{\infty }$ at infinity:At the Bondi sonic point, $GM/r_s=2c_s^2$ and $u_r^2=c_s^2$, so ${\rho }_s=e^{3/2}{\rho }_{\infty }$. The positive inward isothermal Bondi accretion rate is consequentlyThe critical mass density and rate are fixed by the regular transonic solution and the specified reservoir at infinity; an arbitrary subsonic solution does not share this accretion rate.

Use spatial indices $i,j=1,2,3$, set ${\varphi }_i={\partial }_i\varphi$, and assume $\varphi$ is smooth. The inverse metric has $g^{tt}=-c^{-2}{e}^{-2\varphi /c^2}$ and $g^{ij}={\delta }^{ij}$. The Christoffel symbols of the Levi-Civita connection follow from ${\Gamma }^{\alpha }{}_{\beta \gamma }=g^{\alpha \delta }({\partial }_{\beta }{g}_{\gamma \delta }+{\partial }_{\gamma }{g}_{\beta \delta }-{\partial }_{\delta }{g}_{\beta \gamma })/2$:All other components vanish. Although $g_{tt}=-c^2-2\varphi +O(c^{-2})$ diverges, its affine connection has a smooth limit on compact subsets:This is the Newtonian connection from an exponential lapse. It is torsion-free, and its geodesic equation, using $t$ as an affine parameter when $\dot{t}\ne 0$, is $d^2{x}^i/d{t}^2=-{\partial }_i\varphi$. Thus $\varphi$ has the role of a Newtonian gravitational potential. There is no finite limiting nondegenerate Lorentzian metric in these fixed coordinates; $-g(c)/c^2$ tends to $d{t}^2$, a rank-one tensor. The limiting connection preserves $dt$ and the contravariant spatial tensor with components $h^{ij}={\delta }^{ij}$, $h^{t\mu }=0$, which describe the degenerate temporal/spatial structures of this limit.

Let the constant mass density be $\rho >0$ and put ${\mathbf{v}}_a=\mathbf{B}/\sqrt{{\mu }_0\rho }$, the Alfvén velocity. Since both the velocity and the magnetic field have zero divergence, the magnetic tension and magnetic pressure decomposition of the Lorentz force givesThe Newtonian gravitational potential $\Phi$ remains in $\psi$; uniform mass density does not justify dropping a prescribed gravitational acceleration. Adding and subtracting the two equations proves the Elsässer variable equationsHere ${\mathbf{z}}^{\pm }=\mathbf{u}\pm {\mathbf{v}}_a$; each Elsässer variable is transported by the other.

Define $I_{\pm }={\int }_V|{\mathbf{z}}^{\pm }{|}^{2}dV$. Taking the dot product with $2{\mathbf{z}}^{\pm }$ gives the Elsässer energy invariant balanceThe divergence theorem proves $d{I}_{\pm }/dt=0$ whenever the net boundary flux vanishes. In particular, $\mathbf{u}\cdot \mathbf{n}=\mathbf{B}\cdot \mathbf{n}=0$ makes ${\mathbf{z}}^{\pm }\cdot \mathbf{n}=0$, so both fluxes vanish pointwise. Periodic boundary conditions or decay at infinity are also sufficient. These boundary conditions matter: fixed volume alone gives no conservation.

The kinetic energy plus magnetic energy isExpanding $|\mathbf{u}\pm {\mathbf{v}}_a{|}^2$ and using the cross-helicity definition yieldsBoth energy and cross-helicity therefore follow from the two Elsässer energy invariants. The quantity $E$ here is exactly the requested kinetic energy plus magnetic energy, without an additional gravitational term.

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