Chapter 19 The method of characteristics: first order PDEs II

  • Example 19.1. For given \(a\in \mR \) and \(u_0:\mR \to \mR \), we consider the transport equation

    \[ a\partial _xu+\partial _yu=0,\quad u(x,0)=u_0(x). \]

    The corresponding characteristic equations are

    \[ \dot {x}=a,\qquad \dot {y}=1,\qquad \dot {z}=0. \]

    The data curve \(\mR \times \{0\}\) is parametrized by \(\gamma (r)=(r,0)\) with \(r\in \mR \), so that the initial condition gives rise to

    \[ x(r,0)=r,\qquad y(r,0)=0,\qquad z(r,0)=u_0(r). \]

    The general solutions of the characteristic equations are

    \[ x(t)=at+c_1,\qquad y(t)=t+c_2,\qquad z(t)=c_3, \]

    for constants (in \(t\)) \(c_1\), \(c_2\) and \(c_3\); we emphasize that these “constants” may depend on \(r\). Determining the constants from the initial conditions gives

    \[ x(r,t)=at+r,\qquad y(r,t)=t,\qquad z(r,t)=u_0(r). \]

    We can solve for \(r\) and \(t\) in terms of \(x\) and \(y\) and obtain \(t=y\), \(r=x-ay\). Substituting in \(z\) we obtain

    \[ u(x,y)=u_0(x-ay). \]

    We can indeed verify that \(u(x,0)=u_0(x)\) and that since \(\partial _x u=u_0'(x-ay)\) and \(\partial _y u=-au_0'(x-ay)\) we have \(a\partial _x u+\partial _y u=au_0'(x-ay)-au_0'(x-ay)=0\).

  • Example 19.2. For given \(u_0:\mR \to \mR \), we consider

    \[ a\partial _xu+\partial _yu=u^2,\quad u(x,0)=u_0(x). \]

    The corresponding characteristic equations are

    \[ \dot {x}=a,\qquad \dot {y}=1,\qquad \dot {z}=z^2, \]

    with initial conditions

    \[ x(r,0)=r,\qquad y(r,0)=0,\qquad z(r,0)=u_0(r). \]

    The general solutions of the characteristic equations are

    \[ x(t)=at+c_1,\qquad y(t)=t+c_2,\qquad z(t)=\frac {-1}{t+c_3}, \]

    for constants (in \(t\)) \(c_1\), \(c_2\) and \(c_3\). Determining the constants from the initial conditions gives

    \[ x(r,t)=at+r,\qquad y(r,t)=t,\qquad z(r,t)=\frac {-1}{t-\frac {1}{u_0(r)}}=\frac {u_0(r)}{1-tu_0(r)}. \]

    We can solve for \(r\) and \(t\) in terms of \(x\) and \(y\) and obtain \(t=y\) and \(r=x-ay\). Substituting in \(z\) we obtain

    \[ u(x,y)=\frac {u_0(x-ay)}{1-yu_0(x-ay)}. \]

    Note that the solution will blow up when \(1-yu_0(x-ay)=0\). We note that this blow-up is due to the blow-up in the ODE for \(z\).

  • Example 19.3. We consider Burgers’ equation

    \[ u\partial _xu+\partial _yu=0,\quad u(x,0)=-x. \]

    We first more generally consider \(u(x,0)=u_0(x)\). The corresponding characteristic equations are

    \[ \dot {x}=z,\qquad \dot {y}=1,\qquad \dot {z}=0. \]

    with initial conditions

    \[ x(r,0)=r,\qquad y(r,0)=0,\qquad z(r,0)=u_0(r). \]

    The general solutions of the characteristic equations are

    \[ x(t)=c_1t+c_2,\qquad y(t)=t+c_3,\qquad z(t)=c_1, \]

    for constants (in \(t\)) \(c_1\), \(c_2\) and \(c_3\). Determining the constants from the initial conditions gives

    \[ x(r,t)=u_0(r)t+r,\qquad y(r,t)=t,\qquad z(r,t)=u_0(r). \]

    When \(u_0(r)=-r\), we obtain

    \[ x(r,t)=-rt+r,\qquad y(r,t)=t,\qquad z(r,t)=-r. \]

    This we can solve for \(r\) and \(t\) in terms of \(x\) and \(y\) and obtain \(t=y\) and \(r=\frac {x}{1-y}\). Substituting in \(z\) we obtain

    \[ u(x,y)=\frac {x}{y-1}. \]

    The Jacobian determinant equals

    \[ \dot {x}\partial _ry-\dot {y}\partial _rx=t-1. \]

    This equals zero if \(t=1\), which corresponds to \((x,y)=(0,1)\).

    We note that as \(y\uparrow 1\), the solution blows up. We note that in contrast to Example 19.2, in Burgers’ equation, there is no blow-up in the ODEs. The blow-up in the PDE is instead due to the coordinate transformation \((r,t)\mapsto (x,y)\) becoming non-invertible.