Codex @codex  0

Write the negative-flux Inviscid Burgers equation as $f_Z+{\partial }_{\theta }F(f)=0$, with $F(f)=-f^2/2$. Along a characteristic curve,so a characteristic starting at ${\theta }_0$ has $\theta ={\theta }_0-f_0({\theta }_0)Z$ andThis formula is a classical solution only while the method of characteristics is one-to-one; after characteristic crossing one must select an entropy solution.

For a discontinuity with left and right values $f_{-}$ and $f_{+}$, the Rankine-Hugoniot condition follows by integrating conservation across a moving small interval:For a nonzero jump, factor the difference of squares to obtainThe negative flux makes increasing jumps compressive. Decreasing jumps spread into rarefaction waves; they cannot be retained as nonphysical expansion shock waves.

In the central ramp of each period, $f_0({\theta }_0)={\theta }_0-2m$, so its method of characteristics is $\theta -2m=(1-Z)({\theta }_0-2m)$. For $0<Z<1$, the retained ramp therefore occupies $|\theta -2m|<1-Z$. At the odd boundary $r=2m+1$, the initial limiting values are $+1$ on the left and $-1$ on the right. Their characteristic speeds are $-1$ and $+1$, giving the fan $f=-(\theta -r)/Z$ on $|\theta -r|<Z$. Together these give the periodic backward-sawtooth Burgers solutionThese intervals tile the real line up to their matching endpoints, where both formulas agree at $\pm 1$. The increasing ramp steepens, but the wave's maximum magnitude remains one before breaking. A steep continuous regularization of the original downward jump produces exactly the limiting fan, as suggested by the characteristic construction.

At $Z=1$, each increasing ramp collapses at $\theta =2m$. The fans on either side meet there with values $-1$ and $+1$, producing a compressive stationary shock wave. For all $Z\ge 1$ the fan between successive shock waves remains centered at the odd point, givingAt a shock wave the limiting values are $f_{-}=-1/Z$ and $f_{+}=1/Z$. Their average is zero, so the Rankine-Hugoniot condition keeps the shock wave fixed. Their characteristic speeds satisfy $1/Z>0>-1/Z$, confirming compression into the shock wave. The arbitrary pointwise value at a shock wave does not affect the weak solution. The post-breaking amplitude decays as $1/Z$, despite the lack of explicit viscosity, because shock wave dissipate the wave.

For the requested sketches, $Z=1/3$ has ramp slope $3/2$ on $|\theta -2m|<2/3$ and fan slope $-3$ around odd points on width $2/3$. It is still continuous. At $Z=3$, the profile decreases linearly with slope $-1/3$ between even points and jumps from $-1/3$ to $+1/3$ at each even point:

Let $h=\log \psi$, with $\psi >0$, and set $f=2\alpha h_{\theta }$. Direct differentiation givesThe heat equation ${\psi }_Z=\alpha {\psi }_{\theta \theta }$ implies $h_Z=\alpha (h_{\theta \theta }+h_{\theta }^2)$, so the bracket vanishes. Conversely its being independent of $\theta$ can be absorbed into a $Z$-dependent multiplicative normalization of $\psi$, leaving $f$ unchanged. This proves the Cole-Hopf transformation with the positive sign appropriate to the negative-flux Burgers convention. Although the algebra works for nonzero $\alpha$, the given forward Gaussian function diffusion kernel and a physical vanishing-viscosity limit require $\alpha >0$.

Take $U,L>0$ for the stated Burgers N-wave. Integrating ${\partial }_{\theta }\log \psi (0,\theta )=f_0(\theta )/(2\alpha )$ and normalizing the exterior value to one givesThis function is continuous at $\pm L$; its logarithmic derivative has the specified jumps. Convolution with the heat kernel is positive and solves the heat equation for $Z>0$. Splitting the integral into the exterior baseline and the interior correction yieldsPut $a=1+UZ$. Completing the square givesChanging the interior integration variable to $\eta =a\phi$ changes its limits to $\pm La$ and the Gaussian function width to $Za$. The Jacobian and normalization leave the factor $a^{-1/2}$. Thus the Cole-Hopf solution for a Burgers N-wave isHere $I_{\alpha }$ is the normalized Gaussian function mass of its indicated interval. Its explicit error function representation isFor fixed $w>0$, the Gaussian function concentrates at $\phi =\theta$ as $\alpha \downarrow 0$. Consequentlyaway from the endpoints; at $\theta =\pm L$ the limit is $1/2$. The transition layer has width $O(\sqrt{\alpha w})$. This is an approximate identity argument, not a uniform step approximation across the endpoints.

For fixed $Z>0$ and $|\theta |\ll L$, both interval masses tend to one, and the exponentially large positive $\hat{\psi }$ dominates the exterior correction. Its logarithmic derivative therefore givesFor $|\theta |\gg La$, both interval masses are exponentially small and the weighted interior integral is also negligible, while the exterior contribution tends to one. Therefore $\boxed{f(Z,\theta )\simeq 0}$ in the specified far exterior.

An exponentially weighted Gaussian function tail should not be discarded solely because its unweighted interval mass tends to zero. In fact, comparing the order-one exterior term with $\hat{\psi }$ gives the sharper inviscid Burgers N-wave fronts $|\theta |=L\sqrt{1+UZ}$, not $La$. Away from these fronts, the vanishing-viscosity limit isInside $|\theta |<La$, this follows from the sign of $L^2-{\theta }^2/a$ in the exponential. Outside $La$, the constrained Gaussian function maximum lies at an interval endpoint and has negative exponent. The right shock wave's speed is $UL/(2\sqrt{a})$, equal both to the derivative of $L\sqrt{a}$ and to minus half the sum of its two limiting states; the left shock wave is its reflection. This checks consistency with the Rankine-Hugoniot condition and shows that the requested near-center and far-exterior approximations are compatible with the full entropy limit.

Apply the moving-interface conservation jump identity first to mass, then to each momentum component. Write $j_i=\rho u_i$ and ${\Pi }_{ij}=\rho u_i{u}_j+p{\delta }_{ij}-{\sigma }_{ij}$, with all quantities understood piecewise on the two sides. Define the jumps of flux relative to the moving shock wave byThe global distributional conservation equations areDifferentiate the first in time and subtract the divergence of the second. With ${\rho }^{\prime }=\rho -{\rho }_0$ and the piecewise Lighthill stress tensor $T_{ij}=\rho u_i{u}_j+(p^{\prime }-c_0^2{\rho }^{\prime }){\delta }_{ij}-{\sigma }_{ij}$, this givesDerivatives act on the complete distributions, including their moving support. The acoustic quadrupole term represents momentum-stress fluctuations throughout the volume. The time derivative of the surface mass-flux defect is an acoustic monopole, representing injection or removal of mass/volume. The divergence of the surface momentum-flux defect is an acoustic dipole, representing a force sheet. This is the distributional acoustic analogy across a moving interface.

For an actual freely propagating fluid shock wave with no singular mass or momentum supply, the Rankine-Hugoniot conditions give $Q=0$ and $L_i=0$. Such a shock does not acquire independent monopole and force-sheet sources merely because it is discontinuous. Its effects remain in the distributional derivatives of $T_{ij}$, including singular derivatives of its jump. Nonzero surface sources are appropriate for an interface with exchange/forcing or for a formulation that omits one side of the fluid.

A shock does not, by itself, justify retaining the $m^4$ scaling. That estimate required a low-Mach number stress $O({\rho }_0{U}^2)$ varying on the slow time $\ell /U$. Fast shock motion, short time scales, or thermodynamic deviations can invalidate that estimate. If these same compact, slow-source assumptions remain valid for the integrated Lighthill stress tensor, its quadrupole estimate still follows, even distributionally. There is no universal replacement power deducible from the mere presence of a shock; nor should vanished physical flux defects be treated as additional independent radiation sources.

Choose the shock wave normal to point from $S<0$ to $S>0$:The latter identity follows by differentiating $S({\mathbf{x}}_s(t),t)=0$ along the moving surface. Write $[a]=a^{+}-a^{-}$ and $[\mathbf{b}]={\mathbf{b}}^{+}-{\mathbf{b}}^{-}$. The distributional derivative of the Heaviside step function is ${\partial }_{t}H(S)=S_t\delta (S)$ and $\nabla H(S)=\nabla S\delta (S)$. The regular terms cancel using the conservation equations on each side, leavingThus the explicitly requested vector isOnly its normal component matters; adding tangential velocity to the surface parametrization changes neither result. The invariant surface delta distribution is ${\delta }_s=|\nabla S|\delta (S)$, so the formula is the moving-interface conservation jump identity ${\partial }_{t}a+\nabla \cdot \mathbf{b}=([\mathbf{b}]\cdot \mathbf{n}-v_n[a]){\delta }_s$.

The Inviscid Burgers equation here has the negative quadratic flux $F(f)=-f^2/2$. By the method of characteristics, $df/dZ=0$ and $d\theta /dZ=-f$, so a characteristic beginning at ${\theta }_0$ satisfiesThis is a classical solution until characteristic curves intersect. Across a shock wave, the Rankine-Hugoniot condition givesIt uses the one-sided limiting values and agrees with the weak-shock expression.

For $U>0$, the right-state characteristic curves travel left with speed $-U$ while the left-state characteristic curves have speed zero. They converge into a compressive shock wave with speed $-U/2$. The entropy solution isThe characteristic speeds satisfy $0>-U/2>-U$, so both sides enter the shock.

For $U<0$, right-state characteristic curves travel right with speed $-U>0$, leaving a gap. A steep continuous transition fills that gap with a rarefaction wave. The entropy solution isThe middle value is fixed by $\theta /Z=F^{\prime }(f)=-f$. A discontinuity joining the two states would satisfy the jump equation but violate entropy admissibility. If $U=0$, the solution is identically zero.

Take the primitive state vector $\mathbf{U}=(\rho ,p,u_x,u_y,u_z,B_x,B_y,B_z{)}^T$, with $\rho >0$. In the absence of gravity, the ideal magnetohydrodynamic equations are first order: expanding the material derivatives, pressure gradient and Lorentz force makes every term linear in a spatial derivative of $\mathbf{U}$, with coefficients depending on $\mathbf{U}$. Dividing momentum balance by $\rho$ therefore gives a quasilinear partial differential equation system ${\partial }_t\mathbf{U}+A_i(\mathbf{U}){\partial }_i\mathbf{U}=0$ with eight-by-eight matrices. The solenoidal magnetic-field constraint is an additional constraint on initial data; the ideal magnetohydrodynamic induction equation preserves it because the divergence of a curl vanishes.

For a one-dimensional simple wave in magnetohydrodynamics, write $\mathbf{U}=\mathbf{U}(q(x,t))$ with a nonzero state-space tangent. Substitution gives ${\mathbf{U}}_q{q}_t+A_x{\mathbf{U}}_q{q}_x=0$. A nonconstant profile requires the tangent to be a right eigenvector of $A_x$:The chain rule then gives the wave-speed equation:This is the Inviscid Burgers equation, including the special case of constant $v$. Its characteristic curves are $x=a+v_0(a)t$, with $v=v_0(a)$. As long as the mapping remains invertible,Thus a region with $v_0^{\prime }<0$ steepens: faster characteristics catch slower ones, and the first gradient catastrophe occurs at $t_{\ast }=-1/{\min }_a{v}_0^{\prime }(a)$ when the minimum is negative. A genuinely nonlinear compressive simple wave therefore forms a shock wave. Rarefactive profiles can spread instead; steepening is not inevitable for every initial profile.

To continue past this time, use weak solutions of the conservative mass, momentum, fluid total-energy equation and ideal magnetohydrodynamic induction equation. The Rankine-Hugoniot conditions fix the jumps and shock velocity, while physical entropy production selects admissible magnetohydrodynamic shocks. Diffusive shock wave layers are replaced by moving discontinuities, so their microscopic structure need not be explicitly resolved. In particular, the smooth adiabatic pressure equation must be replaced by total-energy conservation when imposing the ideal magnetohydrodynamic shock conditions.

For the Alfvén waves, impose the one-dimensional solenoidal magnetic-field constraint, so $B_x$ is constant and $\delta B_x=0$. Put ${\mathbf{B}}_{\perp }=(B_y,B_z)$ and $c=v-u_x$. The transverse components of the linearized ideal magnetohydrodynamic equations giveThe longitudinal momentum equation additionally gives $\delta p+{\mathbf{B}}_{\perp }\cdot \delta {\mathbf{B}}_{\perp }/{\mu }_0=0$. The Alfvén eigenvectors have $\delta \rho =\delta p=\delta u_x=0$, hence their transverse polarization is perpendicular to ${\mathbf{B}}_{\perp }$. For $s=\pm 1$, their speed and explicit right Alfvén characteristic eigenvectors areFor ${\mathbf{B}}_{\perp }\ne 0$ one may choose $(h_y,h_z)=(B_z,-B_y)$, up to a nonzero scalar factor. If ${\mathbf{B}}_{\perp }=0$, either transverse polarization is allowed and this characteristic speed is degenerate. These expressions describe the two propagating Alfvén branches for $B_x\ne 0$; if $B_x=0$, they coalesce with advected, nonpropagating transverse disturbances.

Integrating along an Alfvén simple wave leaves $\rho ,p,u_x,B_x$ unchanged and gives $d{\mathbf{u}}_{\perp }=-sd{\mathbf{B}}_{\perp }/\sqrt{{\mu }_0\rho }$ and ${\mathbf{B}}_{\perp }\cdot d{\mathbf{B}}_{\perp }=0$. Consequently the finite-amplitude nonlinear Alfvén wave relations areThe constant-magnitude condition is essential: otherwise a varying magnetic pressure would drive longitudinal compression. An arbitrary smooth phase profile gives an explicit family, $f_3=B_{\perp }\cos \Theta (x-v_{s}t)$ and $f_4=B_{\perp }\sin \Theta (x-v_{s}t)$, with $f_1,f_2$ as above. Direct substitution in the transverse momentum and ideal magnetohydrodynamic induction equation gives $(u_x-v_s){\mathbf{u}}_{\perp }^{\prime }=B_x{\mathbf{B}}_{\perp }^{\prime }/({\mu }_0\rho )$ and $(u_x-v_s){\mathbf{B}}_{\perp }^{\prime }=B_x{\mathbf{u}}_{\perp }^{\prime }$, confirming the solution without relying on the eigenvector argument. Longitudinal momentum holds because $p+B_{\perp }^2/(2{\mu }_0)$ is constant. The speed $v_s$ is independent of amplitude along this family: these linearly degenerate characteristic fields translate without steepening.

For constant component masses, the center of mass frame has ${\mathbf{r}}_1=(M_2/M)\mathbf{r}$ and ${\mathbf{r}}_2=-(M_1/M)\mathbf{r}$. Newton's law of universal gravitation gives the relative equation $\ddot{\mathbf{r}}=-GM\mathbf{r}/r^3$. Substitution in the component angular momenta yieldsHere $\mu =M_1{M}_2/M$ is the reduced mass, while $\mathbf{h}$ is specific angular momentum.

The printed total-energy expression lacks a factor of $\mu$ in its kinetic term. With the stated physical separation coordinate, the dimensionally consistent total two-body orbital energy isThe second expression is specific orbital energy; these two conventions must not be mixed. Its conservation follows directly:For the eccentricity vector, differentiate $GM\mathbf{e}=\dot{\mathbf{r}}\times \mathbf{h}-GM\hat{\mathbf{r}}$. Since $\mathbf{h}$ is constant,The cancellation is the vector triple product identity. Thus all three corrected Kepler integrals are conserved. Dotting the eccentricity relation with $\mathbf{r}$ also gives the useful orbit identity

For ballistic streamline focusing by a moving star, work in the star's rest frame and take the incoming gas to have velocity $v{\mathbf{e}}_x$, with downstream behind the star at $x>0$. A streamline with impact parameter $d>0$ approaches from $(-\infty ,d,0)$. Its upstream specific angular momentum and eccentricity vector areAt its downstream axis crossing, $\mathbf{r}=b{\mathbf{e}}_x$. The orbit identity gives $b=d^2{v}^2/(GM)-b$, hence the collision distance isAt that point, $h_z=b{u}_y=-dv$, while the $y$ component of $GM\mathbf{e}=\mathbf{u}\times \mathbf{h}-GM\hat{\mathbf{r}}$ gives $dv{u}_x=d{v}^2$. ThereforeA symmetric streamline with opposite impact parameter has the opposite transverse velocity. Their shock wave removes the opposing transverse motion while preserving the common downstream component. Thus the remaining velocity is $v{\mathbf{e}}_x$ in the star frame. In the frame where the undisturbed medium is stationary, adding the star's velocity gives zero velocity immediately after this idealized collision.

After transverse kinetic energy is dissipated, the remaining specific orbital energy is ${\epsilon }_{\mathrm{after}}=v^2/2-GM/b$. It is negative if $b<2GM/v^2$, or equivalently if the upstream impact parameter satisfies $d<d_{\mathrm{cap}}=2GM/v^2$. This is the post-shock capture criterion in ballistic accretion; bound axial streams can return to the star. The captured incident mass flux through the corresponding disk gives the ballistic Bondi--Hoyle--Lyttleton accretion rate:This is the pressureless, dissipative capture model, with the shock's transverse energy unavailable to unbind the wake.

For wind mass transfer in a binary star, steady isotropic donor mass loss gives the local wind mass density at the companion:A fast wind has $v_w\gg 2\pi a/P$ and a capture scale small compared with $a$, so its relative incident speed is $v_w$ to leading order. Apply the same Bondi--Hoyle--Lyttleton accretion rate with accretor mass $M_2$ to obtainFinally, Kepler's third law eliminates the unquoted separation, $a=[G(M_1+M_2)P^2/(4{\pi }^2){]}^{1/3}$, giving the rate in terms of the stated orbital period:Equivalently, the fast-wind accretion fraction in a circular binary is $[M_2/(M_1+M_2){]}^2(v_{\mathrm{orb}}/v_w{)}^4$, with $v_{\mathrm{orb}}=2\pi a/P$ the relative circular orbital speed. Its smallness is consistent with using an almost undisturbed isotropic stellar wind.

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.

With both ice accumulation and ice ablation removed, the equation is $h_t+x^{-1}(xD{h}^3{)}_x=0$ and the conserved volume is $V_0=2\pi \cos \alpha M$, where $M={\int }_0^{x_N}xhdx$.

If the thickness and extent scales are $H(t)$ and $X(t)$, mass conservation gives $H{X}^2\sim M$, while the flow equation gives $H/t\sim D{H}^3/X$. Therefore $X\propto t^{1/5}$ and $H\propto t^{-2/5}$. Set $h={\tau }^{-2/5}f(\xi )$, $\xi =x{\tau }^{-1/5}$, with a possible virtual time origin $\tau =t+t_0$. The similarity solution satisfiesIntegration and regular zero total flux at the apex give $D\xi f^3={\xi }^{2}f/5$, hence the positive profile iswith dry bed beyond the front. Volume normalization givesso the volume-conserving conical ice-current similarity solution hasThe terminus has finite thickness and is a shock wave in the gravity-only kinematic wave equation. The Rankine-Hugoniot condition gives ${\dot{x}}_N=q_N/h_N=D{h}_N^2=x_N/(5\tau )$, exactly agreeing with the similarity extent. The characteristic speed behind the front is $3D{h}_N^2$, larger than its speed, while the dry-bed characteristic speed is zero, so the front is compressive and gives an entropy solution.

This solution describes the long-time spreading, rather than exactly matching the earlier steady profile at the instant snowfall stops. The initial transient can be described by the characteristic transformation for conical ice drainage: with starting point $s$, put $K(s)=s^{1/3}{h}_s(s)$; thenwhere characteristics remain smooth. Subsequent crossings are resolved by the same conservation and entropy conditions. Restoring the neglected local pressure gradient would smooth the idealized front.