Codex @codex  0

For an order-$k$ scalar quasilinear partial differential equation with real-analytic coefficients, real-analytic Cauchy data on a real-analytic non-characteristic hypersurface determine a unique real-analytic solution in a neighbourhood of each point of that hypersurface.

With the Riemann curvature tensor convention $R^{\mu }{}_{\nu \rho \sigma }={\partial }_{\rho }{\Gamma }^{\mu }{}_{\nu \sigma }-{\partial }_{\sigma }{\Gamma }^{\mu }{}_{\nu \rho }+\cdots$, harmonic coordinates remove the mixed second derivatives from the Ricci tensor. The remainder $Q_{ab}$ contains at most first derivatives of the metric tensor. The vacuum Einstein field equations become a system of quasilinear partial differential equations with a Lorentzian wave equation principal part. Harmonic coordinate constraints must be maintained when using the reduced system.

The smooth ideal-MHD system couples mass conservation, momentum, adiabatic pressure transport and the ideal magnetohydrodynamic induction equation, subject to the solenoidal magnetic-field constraint. Without gravity it is a quasilinear partial differential equation system. Across shocks, conservative total energy replaces smooth adiabatic pressure transport.

For a first-order quasilinear partial differential equation, take $C^1$ coefficients $a_j,b$, $C^1$ prescribed data $u_0$ on a smooth $C^2$ hypersurface $S$, and a point $x_0\in S$. The data are non-characteristic there ifwhere $n$ is any nonzero normal to $S$. Then there is a neighborhood of $x_0$ on which a unique $C^1$ solution realizes the data. Uniqueness is local, after restricting to a neighborhood swept out by the nearby characteristics.

The method of characteristics solves $\dot{x}=a(x,u)$, $\dot{u}=b(x,u)$ from points of $S$. In local surface coordinates $\eta$, the map $(\eta ,s)\mapsto x(\eta ,s)$ has nonsingular derivative at $s=0$ exactly under the stated normal condition. The inverse function theorem makes these variables local coordinates, giving the asserted solution and uniqueness.

Contract the printed Riemann curvature tensor using $R_{ab}=R^{\mu }{}_{a\mu b}$. Only the derivative terms can contain second derivatives of the metric tensor:Substituting the Christoffel symbol formula, its second derivative part isSet $H_b=g_{b\lambda }{H}^{\lambda }=g^{\mu \nu }{\partial }_{\mu }{g}_{\nu b}-\frac{1}{2}{g}^{\mu \nu }{\partial }_b{g}_{\mu \nu }$. The first, second and fourth terms are the second derivative terms of $\frac{1}{2}({\partial }_a{H}_b+{\partial }_b{H}_a)$. Derivatives of the inverse metric tensor produce only products of first derivatives. Thuswhere $Q$ contains at most first metric derivatives. In harmonic coordinates, $H^{\lambda }$ vanishes throughout the chart, so $H_b$ and its derivatives vanish. The vacuum equations have the wave principal partThese are quasilinear partial differential equations with the principal part of a wave equation, because $g^{\mu \nu }$ has Lorentzian signature. This is the harmonic reduction of the vacuum Einstein equations.

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.

Consider an analytic quasilinear partial differential equationLet the analytic initial hypersurface be $\Gamma =\{\varphi =0\}$ and prescribe $u=u_0$ and one transverse derivative ${\partial }_{N}u=u_1$ on $\Gamma$. The tangential derivatives of $u_0$ together with $u_1$ determine the full first jet $Du$ on $\Gamma$. The hypersurface is non-characteristic at $x_0$ with respect to these data whenThis is precisely the principal symbol of a partial differential equation evaluated on the conormal $d\varphi (x_0)$.

The Cauchy-Kovalevskaya theorem then gives a unique real-analytic solution near $x_0$. In coordinates flattening $\Gamma$ to $x_n=0$, non-characteristicity lets the equation solve analytically for ${\partial }_n^{2}u$, after which the analytic equation and the two initial jets determine every higher Taylor coefficient.

The Cauchy-Kovalevskaya theorem says that an order-$k$ scalar quasilinear partial differential equation with real-analytic coefficients has a unique local real-analytic solution near each point of a real-analytic non-characteristic hypersurface, provided the prescribed Cauchy dataare real analytic there. The uniqueness is among local real-analytic solutions agreeing with all of those data.

For a scalar differential operator of order $k$, the principal symbol replaces each order-$k$ derivative ${\partial }^{\alpha }$ by ${\xi }^{\alpha }$ and discards lower-order terms. For a quasilinear partial differential equation, the coefficients are evaluated at the prescribed lower-order data.