Codex @codex  0

A characteristic caustic is a locus where the coordinate map from initial-data and evolution parameters to physical coordinates has zero Jacobian determinant. The inverse function theorem then fails, and different characteristic branches can reach the same point. A classical solution constructed from initial data must be restricted to its appropriate invertible branch before such a fold. In geometrical optics the corresponding ray-envelope phenomenon is a caustic; in nonlinear first-order equations a fold can signal the onset of multiple gradients or shock formation by characteristic intersection.

The Burgers' equation with constant linear damping has the displayed form. For zero viscosity and $\gamma >0$, characteristics are $u=u_0(a)e^{-\gamma t}$ and $x=a+u_0(a)(1-e^{-\gamma t})/\gamma$. A negative initial slope can make this coordinate map cease to be invertible, producing shock formation by characteristic intersection. The transformation $u=e^{-\gamma t}v$, $\tau =(1-e^{-\gamma t})/\gamma$ converts the inviscid equation to the ordinary inviscid Burgers' equation. Across a thin viscous layer the damping contributes only at lower order to the leading Rankine-Hugoniot condition.

With $B_x=0$, a fast simple wave in magnetohydrodynamics has constant $u_y$ and $B_y/\rho =K$. Its longitudinal velocity satisfies $d{u}_x/d\rho =\pm c_f/\rho$, and ${\rho }_t+(u_x\pm c_f){\rho }_x=0$. Along this state curve the speed derivative is $\pm [c_f/\rho +K^2/(2{\mu }_0{c}_f)]$, which is nonzero. Compressive profiles therefore exhibit shock formation by characteristic intersection, unlike a linearly degenerate characteristic field.

For zero viscosity, the method of characteristics gives $u=u_0(a)e^{-t}$ and $x=a+u_0(a)(1-e^{-t})$. On the sloping part of the initial profile, $u_0(a)=-2a$, so $x=a(2e^{-t}-1)$. While this map is increasing, inversion givesThe sloping interval shrinks from width one to zero while its height decreases from two to one. Its spatial derivative is $-2e^{-t}/(2e^{-t}-1)$ and diverges as $2e^{-t}-1\downarrow 0$. Thus the first shock formation by characteristic intersection occurs atFor $0\le t<t_{\ast }$ the profile consists of a left plateau, a straight descending ramp ending at zero, and a right zero plateau. At $t_{\ast }$ the ramp becomes a vertical jump. This is a damped Burgers equation; viscosity smooths that jump into a thin transition.

Write $u=u_x$, $w=u_y$ and $B=B_y$. The isothermal equation of state replaces the adiabatic thermal equation by $p=c_s^2\rho$. Since the fields depend only on $x,t$, the solenoidal constraint gives ${\partial }_x{B}_x=0$ and the $x$ component of the ideal magnetohydrodynamic induction equation gives ${\partial }_t{B}_x=0$. Hence $B_x$ is constant in both space and time.

The continuity, two momentum and transverse induction equations areHere $B$ is the transverse magnetic component and the constant $B_x$ is the longitudinal component. The longitudinal magnetic force is the gradient of transverse magnetic pressure; the transverse force is magnetic tension. In primitive variables the isothermal coplanar magnetohydrodynamic characteristic matrix isThe imposed $z$ components stay zero because their corresponding equations have zero right-hand sides under this reduction.

Let $c=v-u$ be a characteristic speed relative to the fluid and define the Alfvén velocity ${\mathbf{v}}_a=\mathbf{B}/\sqrt{{\mu }_0\rho }$. To derive the characteristic polynomial, an eigenvector with components $(\delta \rho ,\delta u,\delta w,\delta B)$ obeysEliminating $\delta u$ and $\delta w$ on a nondegenerate branch givesTheir consistency yields $(c^2-c_s^2)(c^2-v_{ax}^2)=v_{ay}^2{c}^2$. Equivalently, direct expansion of $\det (\mathbf{A}-vI)$ gives the same polynomial without dividing by $c$ or a possibly vanishing field component:ThusThe discriminant is $(c_s^2-v_a^2{)}^2+4c_s^2{v}_{ay}^2\ge 0$, and both squared speeds are nonnegative. These are the fast magnetosonic wave and slow magnetosonic wave families, transported with the mean flow. The equations are hyperbolic; on a generic state the four distinct real eigenvalues give four characteristic directions. Special field alignments can make speeds coincide or the slow speed vanish. There is no separate entropy characteristic because the thermal model is isothermal, and the independent out-of-plane Alfvén wave polarization is excluded by $u_z=B_z=0$.

To construct a simple wave in magnetohydrodynamics, take $\mathbf{U}=\mathbf{U}(a)$ along a state-space curve whose tangent is a right eigenvector of one selected branch:Substitution into the partial differential equations then gives the scalar transport equationSuch curves exist locally by integrating the eigenvector field away from degeneracies. For example, parameterizing by $a=\rho$ where its eigenvector density component is nonzero gives the explicit state-curve equationswith signed $c$ equal to one of the four relative wave speeds. These formulas are used only where their denominators are nonzero.

A particularly transparent nonlinear example is an isothermal perpendicular magnetosonic simple wave. Set the allowed constant $B_x=0$, choose $B=K\rho$ and keep $w$ constant. For the fast branch,These relations satisfy the eigenvector equations and therefore all four original equations wheneverThe amplitude-dependent characteristic speed has derivativeThis already proves the existence of genuinely nonlinear waves in the stated system.

More generally, the method of characteristics applied to initial data $a(x,0)=a_0(x)$ givesIf the initial speed decreases towards increasing $x$ somewhere, the denominator reaches zero in finite time. When the negative minimum is attained, the first steepening time isFaster trailing states catch slower leading states, giving shock formation by characteristic intersection. For the right-moving perpendicular fast example, a negative density gradient provides such a compression. Beyond $t_{\ast }$ a smooth simple wave is no longer sufficient and an appropriate shock solution is required. Rarefactive profiles instead spread. A linearly degenerate characteristic field has $v^{\prime }=0$ along its state curve and avoids this mechanism; the conclusion is that compressive nonlinear magnetosonic simple waves steepen, not that every possible MHD waveform must do so.

Put $A=(\gamma +1)/2$. A shock formation by characteristic intersection first occurs when the map from emission time to position loses monotonicity:ThusFor $X=c_0{t}^3/T^2$, write $s=\tau /T$ to obtainMinimizing over $s>0$ gives $s^2=[3(3A-1){]}^{-1}$ and hence