Codex @codex  0

Assume the impulse point is in the interior, $0<a<l$. Expand the displacement in the fixed-end Fourier sine seriesOrthogonality gives the initial velocity coefficient by applying the Dirac delta to the sine function:Substitution in the wave equation yields $T_n^{\prime \prime }+2k{T}_n^{\prime }+{\alpha }_n^2{T}_n=0$. For $k<{\alpha }_n$, its roots are $-k\pm i{\omega }_n$, where ${\omega }_n=\sqrt{{\alpha }_n^2-k^2}$. The two initial conditions determineHence the normalized velocity-impulse Green function for a damped string isFor a mode with critical damping, use the limit $e^{-kt}t$ in place of the time factor. For $k>{\alpha }_n$, put ${\eta }_n=\sqrt{k^2-{\alpha }_n^2}$ and use $e^{-kt}\sinh ({\eta }_{n}t)/{\eta }_n$. These choices solve the same modal ordinary differential equation and initial conditions, so the formula covers arbitrary $k>0$. The impulse initial velocity is a distribution, so initial differentiation of the series is interpreted distributionally.

The printed formula requires correction: the factor depending on $n$ must lie inside the sum, and sine-mode normalization requires $2/l$, rather than $2$. The printed slash after the final sine is also stray. Direct differentiation at $t=0$ verifies the corrected normalization: the initial sine coefficients are exactly those of $\delta (x-a)$.

Suppose $J$ is a Jacobi field vanishing at both endpoints of a nonconstant geodesic segment $[0,t_0]$. Its tangential component $h(t)=⟨J(t),T(t)⟩$ satisfiesbecause $D_{t}T=0$ and the Riemann curvature tensor has the usual skew-adjoint symmetry. The two zero endpoint values imply $h\equiv 0$. Both $J$ and $D_{t}J$ are therefore perpendicular to $T$.

The linear map from these Jacobi fields to $T(0{)}^{\perp }$ given by $J\mapsto D_{t}J(0)$ is injective: $J(0)=0$ together with $D_{t}J(0)=0$ forces $J\equiv 0$ by uniqueness for the linear ordinary differential equation defining a Jacobi field. Its target has dimension $n-1$, soThis proves the endpoint-vanishing normal Jacobi fields have dimension at most n-1 bound. For a constant geodesic the equation is $J^{\prime \prime }=0$, and two zero values instead force $J=0$, so there is no conjugate point.

For a connected smooth Riemannian manifold $(M,g)$ without boundary, let $I(M)$ be the group of all global diffeomorphisms $f$ with $f^{\ast }g=g$, endowed with its compact-open topology. Equivalently these maps preserve the intrinsic distance; the Myers-Steenrod theorem supplies smoothness of distance-preserving bijections and makes $I(M)$ a finite-dimensional Lie group acting smoothly. The group is metric-dependent: a topological manifold can carry a very symmetric metric or one with no continuous symmetries. Connectedness matters here. With infinitely many identical disconnected components, permutations of those components can produce a group which is not a finite-dimensional Lie group.

An isometry is determined by its value and derivative at one point. To prove this, suppose $f(p)=p$ and $d{f}_p=I$. Preservation of the Levi-Civita connection giveswhere the exponential map is defined locally. Thus $f$ is the identity in a normal neighborhood of $p$. The set where its value and derivative agree with the identity is both closed and, by the same argument, open. Connectedness makes it all of $M$. No completeness is needed for this determination argument.

The stabilizer $I(M{)}_p$ consequently acts faithfully and orthogonally on $T_{p}M$, so embeds in $O(n)$. More generally an isometry takes a chosen orthonormal frame at $p$ to one at its image point, providing an injective map into the orthonormal frame bundle. It is a smooth embedding, and the full isometry action is proper. The orbit-stabilizer dimension formula therefore givesIsotropy is compact. If $M$ is compact then $I(M)$ itself is compact, by the same frame description or Arzela-Ascoli applied to isometries and their inverses. Properness also means that the quotient by isometries has well-controlled local orbit structure, even though orbit dimensions may vary.

The infinitesimal equation is the Killing equation. Differentiating $f_t^{\ast }g=g$ at zero givesConversely a vector field satisfying this equation has local flows preserving the metric. A global one-parameter subgroup requires the field to be complete. On a geodesically complete manifold every Killing field is complete: its length is constant along each of its own trajectories, so a finite-time trajectory has bounded speed and stays in a closed bounded ball. The Hopf-Rinow theorem makes that ball compact, and the ordinary differential equation can then be continued. The caveat is real: on the Euclidean interval $(0,1)$, the field ${\partial }_x$ is Killing but its translations are not globally defined for all time. This illustrates Killing fields on incomplete manifolds need not generate global isometries.

The Killing transport identities determine a Killing field from $X(p)$ and the skew-adjoint map $(\nabla X{)}_p$. Along a geodesic the field is a Jacobi field, so uniqueness for its second-order differential equation propagates these initial data. They have $n+n(n-1)/2$ components, reproducing the dimension bound infinitesimally. The important distinction on an incomplete manifold is between all local Killing solutions and the Lie algebra of complete global isometry flows.

In dimension three the bound is six, and rotational isotropy gives stronger information. The Lie subalgebras of $\mathfrak{so}(3)$ have dimensions zero, one or three: identifying the bracket with the cross product, a two-plane is not closed under brackets. Suppose $\dim I(M)=5$. Since an orbit has dimension at most three, its stabilizer would have dimension at least two, and hence exactly three. Its connected isotropy acts as all of $SO(3)$ on the tangent space. But the orbit tangent space is invariant under that action and would have dimension $5-3=2$, impossible because the standard rotational representation is irreducible. Therefore dimension five cannot occur.

If the dimension is six, each orbit has dimension three and each stabilizer dimension three. All orbits are open, so connectedness makes the action transitive. At each point the isotropy is transitive on tangent two-planes, forcing sectional curvature to be independent of the plane; transitivity makes its value independent of the point. Thus a maximally symmetric connected three-manifold has constant sectional curvature. Homogeneous Riemannian manifolds are complete, so its simply connected cover is the round sphere, Euclidean space or hyperbolic space after normalization. Conversely these three simply connected metrics attain dimension six:Here $O^{+}(3,1)$ preserves the chosen sheet of the hyperboloid; it includes orientation-reversing hyperbolic isometries. Constant curvature alone does not imply that a quotient retains six-dimensional global symmetry.

Dimension four likewise forces transitivity: stabilizer dimension three would give a one-dimensional invariant orbit tangent space, again impossible for full rotational isotropy. The stabilizer thus has dimension one and the orbits are three-dimensional. Typical four-dimensional isometry groups occur for $S^2\times \mathbb{R}$ and ${\mathbb{H}}^2\times \mathbb{R}$, where the factor groups have dimensions three and one. The usual Nil metric has three translation symmetries and one rotational isotropy symmetry. In centred Heisenberg coordinates it iswhich exhibits the horizontal rotations. Its distinct vertical and horizontal Ricci eigenvalues bound continuous isotropy by that rotation group, so its full isometry dimension is four. Standard metrics on $\widetilde{S{L}_2(\mathbb{R})}$ and nonround Berger spheres also have dimension four. The standard Sol geometry has a three-dimensional translation group and only finite isotropy, giving dimension three. These examples show how a preferred vertical direction reduces rotational symmetry.

The isometry dimensions in dimension three are precisely $0,1,2,3,4,6$, all realizable. A compact hyperbolic three-manifold gives zero; the product of a circle and a closed hyperbolic surface gives one. For dimension two take $T^2\times (0,1)$ withIts scalar curvature is the constant $-14$, but its Ricci eigenvalues are the distinct constants $-3,-6,-5$ in the two horizontal directions and the vertical direction. Hence connected isotropy is trivial and possible isometries preserve these line distributions. A local isometry preserving them has $t^{\prime }=t+c$ and would require constant dilations $e^{-c}$ and $e^{-2c}$ of the two circle coordinates. A global circle diffeomorphism with such a constant dilation has degree $\pm 1$, forcing $c=0$. The identity component is therefore exactly the two circle translations. A flat three-torus gives dimension three; $S^2\times \mathbb{R}$ gives four; a simply connected constant-curvature model gives six. This construction distinguishes local homogeneity from the smaller global isometry group.

For quotients, let a complete simply connected Riemannian manifold $X$ cover $M=X/\Gamma$. Every isometry of $M$ lifts to an isometry of $X$ which normalizes the deck transformation group. Conversely, every normalizing isometry descends. Two lifts differ by a deck transformation, giving the isometry group of a Riemannian quotient formulaThus one should not replace the isometry group of a quotient by that of its model geometry. For a flat torus, for example,The second factor is finite, whereas translations give the three-dimensional identity component.

There is a useful curvature test for compact manifolds. Contracting and differentiating the Killing equation gives ${\nabla }^{\ast }\nabla X=\mathrm{Ric}(X)$, with the nonnegative rough-Laplacian convention. Integration yields the integrated Bochner identity for Killing fieldsIf Ricci curvature is negative definite, no nonzero Killing field exists. The compact isometry group then has dimension zero and is finite. This explains why closed hyperbolic manifolds have finite global symmetry, despite the six-dimensional isometry group of their universal cover. Compactness is essential to this conclusion.

For a compact flat Riemannian manifold, the same identity makes every Killing field parallel. Parallel fields are exactly the vectors fixed by linear holonomy. ThereforeApplied to the six flat types of Question 1, this gives dimensions three for the torus, one for each of the four nontrivial cyclic types, and zero for the Hantzsche-Wendt manifold. These are identity-component dimensions; disconnected finite symmetries can still be present. Local curvature symmetry, global deck-group constraints, and completeness thus play different roles in determining the isometry group.

Write the perturbation Stokes streamfunction as $\psi =\varphi (r)e^{ikz-i\omega t}$, with $\omega =kc$ and nonzero $k$. Its radial and axial velocity amplitudes are $u=-ik\varphi /r$ and $w={\varphi }^{\prime }/r$. These automatically obey incompressibility. Linearizing the radial and axial Euler equations, with pressure amplitude $P$, givesThusDifferentiating the first expression, the terms involving $U^{\prime }{\varphi }^{\prime }$ cancel. Equating the two expressions for $P^{\prime }$ proves the axisymmetric inviscid pipe stability equation:At the axis, finite smooth axisymmetric velocity requires $\varphi =O(r^2)$, after choosing the irrelevant constant in the streamfunction. At the wall, impermeability gives $\varphi (a)=0$. There is no tangential no-slip condition in this inviscid problem. A decoupled azimuthal disturbance satisfies $ik(U-c)u_{\theta }=0$ and does not introduce a growing discrete mode when $U\ne c$.

For real $k$, suppose $c=c_r+i{c}_i$ with $c_i\ne 0$, so division by $U-c$ is legitimate. Multiply the ordinary differential equation by ${\varphi }^{\ast }/[r(U-c)]$, integrate, and apply integration by parts. The endpoint terms vanish by the preceding regularity and wall conditions, leavingTaking its imaginary part yieldsA growing temporal mode therefore requires a sign change of $(U^{\prime }/r{)}^{\prime }$, rather than a sign change of $U^{\prime \prime }$ alone. This is the cylindrical counterpart of Rayleigh's inflection-point theorem. If this derivative vanishes identically, the positive first integral rules out a nontrivial growing real-wavenumber mode as well. Isolated zeros without a sign change do not make the second integral vanish for a nonzero regular mode.

For the parabolic Poiseuille flow, $U^{\prime }/r=-2U_0/a^2$ is constant. Under the specified exclusion of critical levels, division givesSetting $\varphi =rf$, or equivalently choosing $n=1$ in $r^{-n}\varphi$, turns this intoIt is the Modified Bessel differential equation of order one. Hence $\varphi =r[A{I}_1(kr)+B{K}_1(kr)]$. Since $r{K}_1(kr)$ tends to a nonzero constant at the axis, it gives a singular radial velocity and must be excluded. The admissible mode shape is thereforeFor $k=0$, direct solution gives $\varphi =C{r}^2$ after axis regularity, and the wall condition again forces the trivial solution. For nonzero $k$, the nonzero roots of $I_1$ are imaginary. Writing $j_{1,l}>0$ for the zeros of the ordinary Bessel function $J_1$, the evanescent potential modes of inviscid Poiseuille flow areTheir perturbation vorticity vanishes: ${\partial }_{z}u-{\partial }_{r}w=k^2\varphi /r-({\varphi }^{\prime }/r{)}^{\prime }=0$. These modes grow exponentially in one axial direction and decay in the other. They can describe end-forced or spatially localized potential fields, but neither sign gives a bounded normal mode on a whole infinite pipe. Moreover, the reduced radial equation contains no $c$: it supplies no temporal eigenvalue condition. An arbitrary imaginary part assigned to $\omega =kc$ here is consequently not evidence for a physical temporal instability.

The real-wavenumber calculation does exclude growing discrete axisymmetric modes of this inviscid profile. It does not establish every form of stability: singular neutral disturbances with $c$ in the range of $U$, transient dynamics, non-axisymmetric disturbances, and finite-amplitude transition require separate treatment. The exclusion $U\ne c$ specifically removes the inviscid continuous spectrum.

With viscosity, linearize the Navier-Stokes equations instead and eliminate pressure to obtain a fourth-order cylindrical counterpart of the Orr-Sommerfeld equation. Impose both $\varphi (a)=0$ and ${\varphi }^{\prime }(a)=0$, with appropriate regularity at the axis. The resulting eigenvalue problem involves the Reynolds number and determines complex frequencies, rather than leaving $c$ arbitrary. Viscosity supplies diffusion, no-slip boundary layers and critical-layer regularization, so the inviscid sign argument alone no longer decides its spectrum. That observation does not imply that viscous pipe flow must have a growing linear eigenmode; finite-amplitude pipe transition is a different question.

The positive drift places the boundary layer at the left endpoint. Set $y_{\mathrm{out}}=y_0+ϵy_1+\cdots$ and impose the right boundary condition order by order. The leading ordinary differential equation is $[(1+x)y_0{]}^{\prime }=0$, so $y_0=2/(1+x)$. At the next order,which integrates toThe outer limit at the left endpoint is two, so it cannot satisfy the left boundary value without a rapidly varying correction.

For the inner expansion use $X=x/ϵ$ and $y=Y_0(X)+ϵY_1(X)+\cdots$. Its equations and boundary values areMatching the leading constant gives $Y_0=2-e^{-X}$. The first-order forcing is then $-2+(1-X)e^{-X}$. A particular solution is $-2X+\frac{1}{2}{X}^2{e}^{-X}$, and matching determines the constant to be $3/2$. Enforcing $Y_1(0)=0$ givesIn the overlap region the common expansion is $2+ϵ(3/2-2X)$. Adding the inner and outer approximations and subtracting that common part gives the uniform asymptotic approximationIt has error $O({ϵ}^2)$ uniformly for $0\le x\le 1$, apart from an exponentially small endpoint adjustment of the same asymptotic irrelevance. At $x=0$ it equals one exactly, and at $x=1$ it differs from one only exponentially. The polynomial multiplying the exponential remains bounded on the inner scale.

There is an independent check from a first integral. The equation is $[ϵy^{\prime }+(1+x)y{]}^{\prime }=0$, so its exact solution has the formwhere $C$ is selected by the right boundary value. Endpoint expansion of this integral gives $C=2-ϵ/2+O({ϵ}^2)$, the same outer coefficients and the same inner expansion. This check also establishes the uniform order of the composite approximation without treating an inner residual as an outer estimate.

Along the exact ordinary differential equation solution, the chain rule gives $g(y)=f^{\prime }(y)f(y)=y^{\prime \prime }$. Insert the exact solution into the symmetric two-step two-derivative formula and expand about the middle time $t=t_{n+1}$. Its unnormalized local truncation error isThe coefficients of $h{y}^{\prime }$, $h^3{y}^{\prime \prime \prime }$ and $h^5{y}^{(5)}$ are respectivelyThe centered residual is an odd function of $h$, so every even power cancels. The next coefficient isThus the method has order six. Its zero-stability roots at $h=0$ are the simple roots $1,-1$. With sufficiently accurate starts and smooth derivative evaluation maps, convergence of a zero-stable multiderivative method gives sixth-order error on fixed time intervals. The first Dahlquist barrier for ordinary linear multistep methods does not apply, because this is a multiderivative multistep method.

A constant-step linear multistep method has the formwith ${\alpha }_k\ne 0$. The characteristic polynomials of a linear multistep method encode both accuracy and error propagation. The method has order of a numerical method $p$ when insertion of a sufficiently smooth exact solution leaves an unnormalized local truncation error $O(h^{p+1})$. The exponential-symbol order criterion for a multistep method expresses this as $\rho (e^z)-z\sigma (e^z)=O(z^{p+1})$. In particular, first-order consistency of a numerical method requires $\rho (1)=0$ and ${\rho }^{\prime }(1)=\sigma (1)$.

Accuracy alone is insufficient. For $f=0$, perturbations obey the recurrence $\rho (E)e_n=0$. The root condition for a multistep method requires every root of a polynomial $\rho$ to lie in the closed unit disk, with each unit-circle root simple. This is zero-stability. A root outside the unit disk causes geometric amplification as $n\sim T/h$ increases; a repeated unit-circle root permits polynomial growth in $n$. Either can destroy convergence even for starting perturbations tending to zero.

The Dahlquist equivalence theorem states that, for a fixed-coefficient method and a well-posed ordinary differential equation with a suitably Lipschitz continuous right-hand side, consistency of a numerical method and zero-stability are equivalent to convergence for all convergent starting data. If the starting errors are $O(h^p)$ and the exact solution is sufficiently smooth, the global error is $O(h^p)$ on a fixed time interval. A stability estimate bounds the propagated sum of $O(h^{p+1})$ local defects over $O(1/h)$ steps; the discrete Gronwall inequality controls the dependence of $f$ on the accumulated error. Implicit methods also require the locally consistent, uniquely solvable update branch for sufficiently small $h$.

For example, the explicit Euler method has order one and $\rho (\zeta )=\zeta -1$, satisfying the root condition for a multistep method. The two-step Adams-Bashforth method has order two and $\rho (\zeta )=\zeta (\zeta -1)$, again satisfying it. The trapezoidal rule has order two with one step. The two-step Simpson multistep method, $Y_{n+2}-Y_n=h(f_{n+2}+4f_{n+1}+f_n)/3$, has order four and the simple roots $1,-1$: it is zero-stable, although its parasitic amplification root can spoil long-time decay. These examples distinguish finite-time numerical convergence from absolute stability.

The first Dahlquist barrier gives the maximum possible order of a convergent ordinary real $k$-step method:These bounds assume only first-derivative evaluations, fixed coefficients and zero-stability; the sixth-order multiderivative multistep method in solution is outside that class.

To obtain arbitrarily high convergent order within the ordinary class, increase the number of steps. The Adams-Bashforth method integrates a polynomial interpolation of $k$ past derivative values and has order $k$. The Adams–Moulton method includes the new endpoint and has order $k+1$ for $k$ steps. In both families $\rho (\zeta )={\zeta }^{k-1}(\zeta -1)$, so the root condition for a multistep method holds for every fixed $k$. Starting values can be generated by a Runge-Kutta method of sufficient accuracy. The backward differentiation formula instead differentiates an interpolating polynomial; its order-$k$ members are zero-stable only for $1\le k\le 6$.

Finally, high order and favorable stiff stability are separate design goals. By the Second Dahlquist barrier, an ordinary A-stable linear multistep method has order at most two. One may accept a smaller linear stability domain, use suitable backward differentiation formula members for stiff problems, or switch to high-order implicit Runge-Kutta methods or multiderivative multistep methods. A convergent high-order formula still needs time steps appropriate to its stability properties and the problem's eigenvalues.

A linear, or “three-thirds”, boundary heat transfer closure needs its own dimensional conductance $L>0$, so $q=L\Delta$. One cannot simply keep the dimensions of the previous nonlinear coefficient $K$. With equal conductances at both boundaries, the mean-temperature equation becomesIntegrating this linear ordinary differential equation givesThe relaxation time is $\rho c_{p}h/(2L)$. If the conductances differed, conservation of energy would instead select their weighted boundary temperature; the arithmetic mean follows from the assumed upper/lower symmetry.

Put $G=g({\rho }_p-{\rho }_0)/{\rho }_0$ and $b=G\varphi >0$, the reduced gravity. Use the Boussinesq approximation, hydrostatic pressure, uniform vertical velocity and concentration profiles, negligible fluid entrainment and drag, and a flat bed. Integrating volume conservation, momentum conservation and the particle budget over the layer givesThe hydrostatic pressure resultant is ${\rho }_{0}b{h}^2/2$; the last right-hand side is the deposition flux. With $D={\partial }_t+u{\partial }_x$, these areThe quasilinear coefficient matrix in variables $(h,u,b)$ isIts eigenvalues are $u,u\pm c$, where $c=\sqrt{bh}$. Therefore the characteristic curves and their compatibility ordinary differential equations areFor the last equation the left eigenvectors are $(\pm c/h,1,\pm h/(2c))$, and multiplying the quasilinear equations by them gives the displayed source terms. This is the sedimenting rectangular-channel characteristic compatibility relation; it supplies one differential relation on each wave family, rather than three independently closed scalar equations there. The contact eigenvector has $du=0$ and $d(b{h}^2)=0$. For reference, an equivalent form isOnly for constant $b$ without settling do $u\pm 2c$ become genuine Riemann invariants. With settling, the spatial buoyancy gradient cannot be silently omitted.

A freely swimming body has neither an imposed end couple nor an imposed end transverse force. Its end bending moment is therefore zero, and the transverse force $S=G_x$ is also zero at both ends. These are the four free-end bending boundary conditions; they do not require the end displacement or slope to vanish.

Put $R=m_b{h}_{tt}+\mathcal{D}(m\mathcal{D}h)$, using the physically consistent sign. Integrating $G_{xx}=R$ with the two leading-end conditions givesThe two trailing-end conditions consequently requireThese free-end compatibility for a swimming beam conditions are total transverse force and torque balances. They explain why an arbitrary plausible prescribed $h$ does not suffice: two integration constants can satisfy only two end conditions unless the forcing already satisfies both compatibility integrals.

The missing freedom is the swimming recoil correction. An observed or prescribed deformation $h_d$ relative to the body must be supplemented by a lateral translation $Y(t)$ and a small yaw $\Theta (t)$:To identify $Y$ with the actual centre of mass displacement, first remove the body-mass-weighted mean from $h_d$. One may also remove its body-mass-weighted first moment about $x_c$ to define the yaw reference uniquely. This changes the decomposition of the observed motion, not the physical displacement.

For $\xi =x-x_c$, the recoil contribution to the forcing isInsert $R=R_d+R_{\mathrm{rec}}$ into the two compatibility integrals, using $1$ and $\xi$ as weights. This gives two coupled ordinary differential equations for the two recoil motions. Their principal acceleration matrix isFor a nondegenerate positive mass distribution it is a positive-definite matrix: its quadratic form at $(a,b)$ is ${\int }_0^L(m_b+m)(a+b\xi {)}^{2}dx>0$ unless $a=b=0$. Thus initial translation and yaw determine the recoil through these equations, and the corrected $h$ supplies the two compatibility conditions. The four free-end conditions determine the bending moment together with two rigid recoil motions; they do not overdetermine a physically complete swimming motion.

A spatially uniform pressure is a useful leading approximation when the sound speed makes the acoustic crossing time much shorter than the radiative cooling and thermal conduction times. The motions must remain at low Mach number, and gravity and other stresses must be negligible. Sound then redistributes pressure before appreciable thermal evolution occurs. Small pressure corrections can still drive the slow motion; imposing exactly zero pressure gradient in the full momentum equation at every order would be a stronger assumption.

The equation of state gives the specific internal energy of an ideal gas as $e=p/[(\gamma -1)\rho ]$. With $D/Dt={\partial }_t+v{\partial }_x$, the continuity equation is $D\rho /Dt=-\rho v_x$. Differentiating $e$ therefore yieldsInsert this into the internal energy equation, whose compressional term is $-p{v}_x$, and use spatially uniform pressure, so $Dp/Dt=p_t$. The result isFor temporally constant pressure put $A=\mu p/\mathcal{R}$, so the mass density is $\rho =A/T$. The continuity equation now gives $v_x=T^{-1}DT/Dt$. Multiplication of the last equation by $(\gamma -1)T/(\gamma p)$, with $\lambda ={\lambda }_0{T}^{\alpha }$, givesThis is the isobaric conductive cooling in mass coordinates transformation before changing coordinates.

There is a boundary assumption implicit in the proposed Lagrangian coordinate. For a mass integral from a fixed origin, integration of the continuity equation givesThus $m$ is material when the base mass flux vanishes, as appropriate for the stationary cold reservoir below. If the origin permits mass flux, replace $m$ by $m-{\int }^t\rho (0,s)v(0,s)ds$; this corrected label has zero material derivative. In a material mass coordinate, $D/Dt={\partial }_t{|}_m$ and ${\partial }_x=(A/T){\partial }_m$. In particularThe conduction coefficient therefore becomes $C{\lambda }_0$, with the same $C$ as the cooling coefficient. With $\tau =Ct$, the equation is

In the hot gas there is no volumetric radiative cooling. The initial temperature is uniform and the cold boundary introduces no finite mass length. The remaining nonlinear diffusion equation is invariant under $m\mapsto bm$, $\tau \mapsto b^2\tau$ with the same temperature scale. Its natural similarity solution is thereforeThe far-field condition also recovers the initial temperature at every fixed $m>0$ as $\tau \downarrow 0$. Substitution gives the ordinary differential equationFor $\alpha =1$, this reduces to $f^{\prime \prime }+\xi f^{\prime }/2=0$. Integrating first for $f^{\prime }$ gives $f^{\prime }=K\exp (-{\xi }^2/4)$. The two boundary conditions fix $K=1/\sqrt{\pi }$, hence the error-function isobaric cooling front isHere the error function has normalization $2/\sqrt{\pi }$, as printed in the original PDF; the converted TeX's $2/\pi$ is incorrect. The following original sketch shows the temperature profiles. The cooling layer broadens as $\sqrt{\tau }$, and at fixed positive mass label the temperature falls as time increases, while the remote hot gas remains at $T_0$.

For $\alpha =1$, the magnitude of the heat flux into the cold reservoir isAlthough the thermal conductivity tends to zero there, the physical temperature gradient diverges: their product has this finite limit. The stationary reservoir has no incoming advective energy flux and radiates the conducted energy immediately. Since $T_m(0,\tau )=T_0/\sqrt{\pi {\lambda }_0\tau }$, the radiated power per unit area isThus the requested exponent is $\boxed{k=1/2}$. The initial divergence is the usual idealization of an instantaneous temperature discontinuity in a diffusion equation.

Set $a=2\pi G{\rho }_0$. The axisymmetric Poisson equation for Newtonian gravity integrates to $(R{\Phi }_R{)}^{\prime }=4\pi G{\rho }_{0}R$. Regularity on the axis fixes the integration constant, so ${\Phi }_R=aR$. Radial force balance in the rotating fluid givesWith the stated central pressure, $p_c={\rho }_0{a}^2/(4k^2)$, this is a perfect square:The free surface lies at the first zero of pressure. For the source's positive rotation parameter, $R_0=\sqrt{2\pi G{\rho }_0}/k$; allowing either rotation direction requires $|k|$ in the denominator. The inward effective gravity is ${\Phi }_R-R{\Omega }^2=R(a-k^2{R}^2)$, so $\boxed{g_{\mathrm{eff}}(R_0)=0}$. In particular $p_R(R_0)=0$. This is the critically rotating self-gravitating cylinder.

For a normal mode proportional to $e^{i\omega t+im\varphi }$, the background material derivative becomes $i\sigma$, where $\sigma =\omega +m\Omega (R)$. Linearizing the centrifugal term gives $-2\Omega u_{\varphi }$ in the radial equation. In the azimuthal equation the background velocity gradient and cylindrical geometric term together give $(2\Omega +R{\Omega }^{\prime })u_R=3\Omega u_R$. Since the interior mass density is constant, write $W=p^{\prime }/{\rho }_0+{\Phi }^{\prime }$. The two momentum equations and incompressible flow condition are consequentlyFor $m\ne 0$, the last two equations express the other amplitudes in terms of $u=u_R$:Differentiate the second expression using $R{\Omega }^{\prime }=\Omega$ and $R{\sigma }^{\prime }=m\Omega$, then substitute into the radial momentum equation. Before any division by $\sigma$, the result isAway from corotation this gives the required ordinary differential equationAt corotation the original coupled equations, or the undivided equation, must be used. A regular axisymmetric radial velocity would vanish by the incompressible flow condition; the requested dipolar calculation has $m=1$.

For $m=1$, use the amplitude normalization $u=A(\omega +kR)$, which remains useful when $\omega$ tends to zero. Substitution into the undivided equation gives zero identically. For nonzero $\omega$ this is the source's $1+kR/\omega$ after choosing $A=1/\omega$. The other amplitudes areThese formulas also satisfy the three original equations directly. By the regularity assumption in the question, this is the interior velocity family to which the free surface conditions must be applied.

The gravitational potential of a displaced cylindrical interface supplies an essential part of those conditions. For surface fluid displacement $\eta e^{im\varphi }$, the bulk mass density perturbation vanishes but the displaced mass density step contributes ${\rho }_0\eta \delta (R-R_0)$. The interior and exterior gravitational potential amplitudes solve Laplace's equation, are proportional to $R^m$ and $R^{-m}$ for $m\ge 1$, and are continuous at the surface. Integrating the Poisson equation for Newtonian gravity across the mass density step givesSince the two radial derivatives at the surface are $-m{\Phi }_s^{\prime }/R_0$ and $m{\Phi }_s^{\prime }/R_0$, respectively,It would be incorrect to set the entire potential perturbation to zero merely because the bulk mass density is unchanged.

Let ${\Omega }_s=k{R}_0$, so $a={\Omega }_s^2$. Zero fluid Lagrangian perturbation of surface pressure gives $p_s^{\prime }+\eta p_R(R_0)=0$, hence $p_s^{\prime }=0$. The surface momentum potential and kinematic condition for $m=1$ are thereforeIf $\omega +{\Omega }_s\ne 0$, the kinematic condition fixes $\eta =-iA$. The same value follows if the radial fluid Lagrangian displacement $u_R/[i(\omega +kR)]=-iA$ is required to extend continuously to the free surface, including a possible surface-corotation limit. Comparing the gravitational condition with the interior value of $W_s$ then givesfor a nontrivial perturbation. These dipolar perturbations of a critically rotating cylinder describe a displacement of the whole cylinder's axis. Translating an isolated self-gravitating system gives another equilibrium, so there is no external restoring force and the translation is neutral. A uniform drift of the centre also follows from invariance under a Galilean transformation, explaining the double zero of the frequency equation. The singular $1/\omega$ in the originally chosen velocity normalization does not obstruct the finite-amplitude limiting solution.

There is an admissibility qualification to the frequency conclusion as printed. Regular interior fluid Eulerian perturbation of velocity alone does not justify cancelling the surface factor $\omega +{\Omega }_s$. Eliminating $\eta$ without that cancellation instead givesThe additional surface-corotation case $\omega =-{\Omega }_s$ formally satisfies every displayed velocity equation and the usual linearized surface conditions withHere $u_R(R_0)=p^{\prime }(R_0)=0$, and the velocity is finite on the axis. However, its interior radial fluid Lagrangian displacement tends to $-iA$, whereas the undisplaced surface has $\eta =0$. Thus a regular single-frequency fluid Lagrangian displacement continuous up to the surface excludes this degenerate branch. The intended ${\omega }^2=0$ conclusion holds with that regular-displacement condition, or with surface corotation excluded; the printed velocity regularity assumption by itself leaves this formal exception.