Codex @codex  0

The growth bound is $|F(t,x)|\le K(1+|x|)$. A characteristic curve through $(s,a)$ solvesThe $C^1$ field is locally Lipschitz in space, giving local existence and uniqueness. On either direction of a bounded time interval the Gronwall inequality gives $1+|X(t;s,a)|\le (1+|a|)e^{K|t-s|}$. Thus no trajectory can escape to infinity at finite time. On the resulting compact space-time region, local ordinary differential equation existence continues it, proving global definition for all $t\in \mathbb{R}$.

Uniqueness gives the flow composition law and inverse $X(s;t,\cdot )$. Differentiation with respect to $a$ solves the variational equation, making each map a $C^1$ diffeomorphism. Its positive Jacobian determinant isThese facts constitute the global characteristic flow under linear growth. A uniform global bound on $D_{x}F$ is not required; its bounds on each compact characteristic tube suffice.

Write $X(t,a)=X(t;0,a)$ and $A(t,x)=X(0;t,x)$. If $u_0\in C^1({\mathbb{R}}^d)$, the unique global classical solution isIt is $C^1$, and the chain rule makes its value constant along every characteristic curve. The same calculation proves classical uniqueness. No boundedness of $u_0$ is needed for this classical statement.

For $u_0\in L^{\infty }$, a bounded weak solution on the positive half-line is defined byfor every compactly supported $C^1$ test function on $[0,\infty )\times {\mathbb{R}}^d$. The analogous backward-time identity defines the solution on all real time. The $u\mathrm{div}F$ term is essential: this is transport, not the conservative equation $u_t+\mathrm{div}(Fu)=0$.

The same inverse-flow formula defines the unique weak solution, with $\Vert u(t){\Vert }_{\infty }=\Vert u_0{\Vert }_{\infty }$, equality almost everywhere, and a representative continuous into $L^1$ on each compact spatial set. For a direct uniqueness argument, change variables $x=X(t,a)$ in the weak identity and put $w(t,a)=u(t,X(t,a))$. Using $J_t=(\mathrm{div}F)(t,X)J$ givesinitially for $C^1$ compactly supported $\zeta$. Since $J>0$ and $J,J_t$ are continuous on compact sets, approximation in the spatial parameter permits $\zeta =\eta /J$: both $\zeta$ and ${\zeta }_t$ are uniformly approximable there, which are the only derivatives left in this identity. Hence $\int w{\eta }_t+\int u_0\eta (0)=0$ and $w=u_0$ almost everywhere. This also verifies existence by reversing the argument, and avoids assuming second derivatives of the $C^1$ field. Local $L^1$ continuity follows by approximating $u_0$ by continuous functions on compact flow tubes.

This is the weak transport solution under characteristic coordinates. The weak-strong uniqueness principle says that a weak solution in the specified class agrees with a classical solution having the same data for as long as that classical solution exists. Here it holds globally in the bounded class whenever the classical initial data are also bounded.

Along a curve $s\mapsto (x(s),y(s))$, the chain rule givesTogether with the second-order partial differential equation, these form a linear system for $(u_{xx},u_{xy},u_{yy})$ whose coefficient matrix isA characteristic curve is precisely one along which this system is singular, equivalently where the principal symbol vanishes on a normal covector. Taking the determinant yields

Expanding the equation givesFor the hyperbolic partial differential equation discriminant, $A=1$, $B=\rho$, $C={\rho }^2-1$, soFor a curve $\rho =\rho (t)$, the characteristic hypersurface equation is $({\rho }^{\prime }{)}^2-2\rho {\rho }^{\prime }+{\rho }^2-1=0$, or ${\rho }^{\prime }=\rho \pm 1$. Therefore the two families of characteristic curves areThe values $C=0$ include the characteristic lines $\rho =-1$ and $\rho =1$. Every curve in the first family approaches $\rho =-1$ as $t\to -\infty$; every curve in the second approaches $\rho =1$. The two curves through $(\rho ,t)=(0,0)$ are $\rho =e^t-1$ and $\rho =1-e^t$, shown as the thick curves below. Equivalently, the characteristic coordinates are $e^{-t}(\rho +1)$ and $e^{-t}(\rho -1)$.

Use the shallow-water approximation, a hydrostatic vertically mixed layer, a deep motionless ambient, and no ambient fluid entrainment. Particle loss changes reduced gravity but changes total layer volume only at the discarded dilute-particle order. Likewise thermal expansion changes the equation of state while the leading Boussinesq approximation retains volume conservation. DefineDepth integration of mass conservation, horizontal momentum balance, particle transport, and the heating law givesHere the integrated excess hydrostatic pressure is $b{h}^2/2$. There is no particle source from the bed. Bed drag and the particle contribution to inertia are neglected at this order. The equivalent variable-buoyancy shallow water equations areThe term $h{b}_x/2$ is necessary when heating or particle loss makes the reduced gravity vary horizontally.

For $v=(h,u,\varphi ,\theta {)}^T$, the quasilinear system $v_t+A(v)v_x=S(v)$ hasIts characteristic polynomial factors asThus the characteristic curves have slopesOn the two repeated contact characteristic curves, the particle and heating laws are $D\varphi =-W_s\varphi /h$ and $D\theta =\beta \varphi$. Writing $c_g=\sqrt{bh}$ and $D_{\pm }={\partial }_t+(u\pm c_g){\partial }_x$, the gravity-wave compatibility relations areFor $h>0$ and $b>0$, gravity-wave right eigenvectors can be taken as $(1,\pm c_g/h,0,0{)}^T$. The contact eigenspace has dimension two: $\delta u=0$ and $b\delta h+(h/2)\delta b=0$, with independent $\delta \varphi ,\delta \theta$. Therefore $A$ has a full set of real eigenvectors.

The system is hyperbolic throughout strict static stability, but not strictly hyperbolic because the contact speed is repeated. At $b=0$ the eigenvalues coalesce and this diagonalizability is lost; when $b<0$ the gravity-wave speeds become complex, reflecting the loss of a statically stable lower layer.

Expanding the equation givesso its principal symbol isIf a characteristic curve is locally a graph $y=y(x)$, its conormal is proportional to $(-y^{\prime },1)$. The characteristic equation is thereforeOn each region separated by $y=\pm 1$, separation of variables givesThus all the characteristic curves areThe inner curves approach the horizontal characteristics $y=\pm 1$ as $x\to \pm \infty$; the outer hyperbolic-cotangent branches have vertical asymptotes and also approach $y=\pm 1$. This describes the requested sketch.

The initial line $x=0$ has conormal $(1,0)$, and $p(y,1,0)=1$, so it is a non-characteristic hypersurface at every point. Since the coefficients and prescribed data are real analytic, the Cauchy-Kovalevskaya theorem gives a unique real analytic solution in a neighborhood of each point of $\{x=0\}$, hence in a neighborhood of that line.

The shallow-water approximation requires a small depth-to-horizontal-scale ratio, slow streamwise variation of the channel and interface, small vertical acceleration, and rapid enough cross-sectional mixing for $u$ and $\varphi$ to be representative averages. Near the channel apex, or during the initial collapse at a sharp front, these conditions need not hold. The high Reynolds number permits neglect of leading viscous resistance; weak stirring mixes the particles without appreciable ambient fluid entrainment.

Put $\lambda =\beta /H$. At elevation $z$, the channel width is $2\lambda xz$, so the occupied area and interface width areVolume conservation and suspended-particle conservation giveThe excess hydrostatic pressure is ${\rho }_{0}d\varphi (h-z)$. Its streamwise gradient averaged over the triangular section gives acceleration $-d\varphi h_x-(dh/3){\varphi }_x$: the weighted mean of $h-z$ is $h/3$. Thus the triangular-channel shallow water equations areThe quasilinear coefficient matrix in $(h,u,\varphi )$ isIts three distinct real eigenvalues, for $h,d\varphi >0$, establish strict hyperbolicity:On the material characteristic curve, $d\varphi /dt=-2W_s\varphi /h$. On each of the other two characteristic curves, multiplying by the left eigenvector $(\pm 2c/h,1,\pm dh/(3c))$ gives the compatibility ordinary differential equationThe derivatives are along the indicated characteristic, with $g^{\prime }=d\varphi$. At vanishing concentration, the gravity-wave speeds coalesce and strict hyperbolicity is lost.

The one-dimensional Isentropic Euler equations areFor the polytropic equation of state,The continuity equation therefore implieswhile the momentum equation becomesAdding and subtracting these equations givesThus the Riemann invariantsare constant along the characteristic curves

The one-dimensional Isentropic Euler equations areDefine the Riemann invariants for one-dimensional isentropic flowCombining the two equations directly givesso $R_{\pm }$ is constant on the corresponding characteristic curve $dx/dt=u\pm c$.

Initially both invariant disturbances are supported in $0\le x\le L$. Once the two characteristic families have separated, the left-moving packet is a simple wave with $R_{+}$ at its equilibrium value, while the right-moving packet is a simple wave with $R_{-}$ at equilibrium. Between and outside them, both invariants retain their equilibrium values. The schematic spacetime diagram below uses the equilibrium slopes $\pm c_0$ and shows separation after $t=L/(2c_0)$.

Let $\tau$ be the time at which a right-running characteristic curve leaves the piston. The boundary velocity and the characteristic it emits areWiththese relations becomeBecause $d\varphi /d\theta =1-ϵ\cos \theta >0$ and $\varphi (\theta +2\pi )=\varphi (\theta )+2\pi$, $U$ is single-valued and $2\pi$-periodic behind the wavefront. Ahead of the first characteristic, where $\varphi <0$, the fluid remains undisturbed and $U=0$. The implicit relation is Kepler equation; unless $ϵ=0$, its inverse is not a pure sinusoid.

Series reversion givesExpanding the cosine then givesand hence

The principal symbol ofis $p(x,t;\xi ,\tau )=(1-x^2){\xi }^2$. For a regular curve $s\mapsto (x(s),t(s))$, the characteristic curve equation isAway from $x=\pm 1$ this gives $t=\text{constant}$. The two degenerate lines $x=1$ and $x=-1$ are also characteristic. These are the characteristic curves, apart from reparametrization and pieces joined at the degenerate lines.

Define the concentration-dependent reduced gravityUnder the Boussinesq approximation, the three shallow water equations can be written in advective form asThe coefficient matrix of this quasilinear system has characteristic speedsThus this is a hyperbolic system when $g^{\prime }h>0$. Along the intermediate characteristic curve $dx/dt=u$, the concentration obeysIt decreases because ambient entrainment dilutes the chemical; detrainment does not change the concentration of a well-mixed parcel.

For $u_0>u_A$, define the Alfvénic Mach-cone slopeThe equation in part (b) is the hyperbolic partial differential equationwhose characteristic curves are $x\pm mz=\text{constant}$. Let $W(s)$ be one for $|s|<L/2$ and zero otherwise. Selecting the characteristics that trail downstream, toward negative $x$, and imposing part (c) givesThe stationary induction and mass conservation equations then givewith $\delta u_z=\delta u_y=\delta B_y=0$. In the strongly super-Alfvénic limit, $\delta u_x\simeq u_{A}b/B_0$ and $\delta B_z\simeq b{u}_A/u_0$ inside the disturbed region.

The perturbations have support only whereBecause the wing is infinite in $y$, this is the union of two inclined slabs bounded by the characteristic planes $x+m|z|=\pm L/2$. These two trailing slabs are the Alfvén wings generated by the conductor.

Now the mass per unit radius satisfies the nonlinear scalar conservation lawIts characteristic curves have speedBecause $v_m<0$ and $|q^{\prime }|$ increases with $M_r$, the high central part of a localized pulse travels inward faster than its low-density edges. On the inner, rising side, faster high-density characteristics catch slower low-density ones, causing characteristic crossing and an inward-facing shock wave. On the outer, falling side, the characteristics separate and form a rarefaction wave. The initial bump therefore develops a sharp inner front and a broad outer tail while moving toward the star.

For $s\ge 0$ and $\xi \in \mathbb{R}$, define the characteristic curve $X(t;s,\xi )$ byThe bounded derivative $F_x$ makes $F(t,\cdot )$ globally Lipschitz, uniformly in $t$. On each finite time interval, $|F(t,x)|\le |F(t,0)|+L|x|$, so Gronwall inequality prevents finite-time escape. The Picard-Lindelöf theorem therefore gives a unique trajectory for every finite $t\ge 0$. Differentiation in $\xi$ givesso the characteristic flow map is a $C^1$ increasing diffeomorphism.

Along a characteristic, the chain rule changes the equation intoTracing $(t,x)$ backward by the flow therefore givesThe regularity of the flow makes this a $C^1$ classical solution. Conversely, every classical solution obeys the same ordinary differential equation along every characteristic, so the formula also proves uniqueness.

The characteristic curves are $x=\xi +at$. Along one of them, $U(t)=u(t,\xi +at)$ obeys the separable ordinary differential equationThereforeFor $0\le t<1$, the denominator is positive everywhere. At $t=1$, it first vanishes where $\cos \xi =1$, namely at $\xi =2\pi k$. The solution consequently has finite-time blowup at timealong the points $x=a+2\pi k$. No finite classical solution can continue through that time.

The characteristic flow map solves the ordinary differential equationand henceAlong this characteristic curve, the chain rule givesThe value is therefore constant, and tracing $(t,x)$ back to time zero gives the classical solutionDirect differentiation verifies both the linear transport equation and its initial value.