Processing math: 100%

2.1 HW 1

PDF (letter size)

PDF (legal size)

2.1.1 Problem 1.8a

Find all quadratic polynomial solutions of the 3D Laplace equation $\frac{\partial ^{2}u}{\partial x^{2}}+\frac{\partial ^{2}u}{\partial y^{2}}+\frac{\partial ^{2}u}{\partial z^{2}}=0$

Solution

A quadratic polynomial in variables $x,y,z$ is $$\begin{equation} u=a_{1}+a_{2}x+a_{3}y+a_{4}z+a_{5}x^{2}+a_{6}y^{2}+a_{7}z^{2}+a_{8}xy+a_{9}xz+a_{10}yz \tag{1} \end{equation}$$ Hence $u_{x}=a_{2}+2a_{5}x+a_{8}y+a_{9}z$ which implies that $u_{xx}=2a_{5}$. Similarly $u_{y}=a_{3}+2a_{6}y+a_{8}x+a_{10}z$, therefore $u_{yy}=2a_{6}$. And finally $u_{z}=a_{4}+2a_{7}z+a_{9}x+a_{10}y$ and $u_{zz}=2a_{7}$. Substituting these results in the Laplace equation gives above result in $$\begin{align*} 2a_{5}+2a_{6}+2a_{7} & =0\\ a_{5}+a_{6}+a_{7} & =0 \end{align*}$$

Therefore $a_{5}=-\left ( a_{6}+a_{7}\right )$. Using this relation back in (1) gives $$\begin{align*} u & =a_{1}+a_{2}x+a_{3}y+a_{4}z-\left ( a_{6}+a_{7}\right ) x^{2}+a_{6}y^{2}+a_{7}z^{2}+a_{8}xy+a_{9}xz+a_{10}yz\\ & =a_{1}+a_{2}x+a_{3}y+a_{4}z+a_{6}\left ( -x^{2}+y^{2}\right ) +a_{7}\left ( -x^{2}+z^{2}\right ) +a_{8}xy+a_{9}xz+a_{10}yz \end{align*}$$

Which can be written as $$u\left ( x,y,z\right ) =A_{1}+A_{2}x+A_{3}y+A_{4}z+A_{5}\left ( y^{2}-x^{2}\right ) +A_{6}\left ( z^{2}-x^{2}\right ) +A_{7}xy+A_{8}xz+A_{9}yz$$

2.1.2 Problem 1.7

Find all real solutions to 2D Laplace equation $u_{xx}+u_{yy}=0$ of the form $u=\log \left ( p\left ( x,y\right ) \right )$ where $p\left ( x,y\right )$ is a quadratic polynomial.

Solution

A quadratic polynomial $p\left ( x,y\right )$ in variables $x,y$ is $$p\left ( x,y\right ) =a_{1}+a_{2}x+a_{3}y+a_{4}x^{2}+a_{5}y^{2}+a_{6}xy$$ Therefore $$u\left ( x,y\right ) =\log \left ( a_{1}+a_{2}x+a_{3}y+a_{4}x^{2}+a_{5}y^{2}+a_{6}xy\right )$$ Hence $$u_{x}=\frac{a_{2}+2a_{4}x+a_{6}y}{p\left ( x,y\right ) }$$ and $$\begin{equation} u_{xx}=\frac{2a_{4}}{p\left ( x,y\right ) }-\frac{\left ( a_{2}+2a_{4}x+a_{6}y\right ) ^{2}}{p\left ( x,y\right ) ^{2}} \tag{1} \end{equation}$$ Similarly $$u_{y}=\frac{a_{3}+2a_{5}y+a_{6}x}{p\left ( x,y\right ) }$$ And $$\begin{equation} u_{yy}=\frac{2a_{5}}{p\left ( x,y\right ) }-\frac{\left ( a_{3}+2a_{5}y+a_{6}x\right ) ^{2}}{p\left ( x,y\right ) ^{2}} \tag{2} \end{equation}$$ Substituting (1,2) into $u_{xx}+u_{yy}=0$ gives $$\begin{align*} \left ( \frac{2a_{4}}{p\left ( x,y\right ) }-\frac{\left ( a_{2}+2a_{4}x+a_{6}y\right ) ^{2}}{p\left ( x,y\right ) ^{2}}\right ) +\left ( \frac{2a_{5}}{p\left ( x,y\right ) }-\frac{\left ( a_{3}+2a_{5}y+a_{6}x\right ) ^{2}}{p\left ( x,y\right ) ^{2}}\right ) & =0\\ 2a_{4}-\frac{\left ( a_{2}+2a_{4}x+a_{6}y\right ) ^{2}}{p\left ( x,y\right ) }+2a_{5}-\frac{\left ( a_{3}+2a_{5}y+a_{6}x\right ) ^{2}}{p\left ( x,y\right ) } & =0\\ 2a_{4}+2a_{5}-\frac{\left ( a_{2}+2a_{4}x+a_{6}y\right ) ^{2}+\left ( a_{3}+2a_{5}y+a_{6}x\right ) ^{2}}{p\left ( x,y\right ) } & =0 \end{align*}$$

Or $$\left ( 2a_{4}+2a_{5}\right ) p\left ( x,y\right ) =\left ( a_{2}+2a_{4}x+a_{6}y\right ) ^{2}+\left ( a_{3}+2a_{5}y+a_{6}x\right ) ^{2}$$ But $p\left ( x,y\right ) =a_{1}+a_{2}x+a_{3}y+a_{4}x^{2}+a_{5}y^{2}+a_{6}xy$. Hence the above becomes $$\left ( 2a_{4}+2a_{5}\right ) \left ( a_{1}+a_{2}x+a_{3}y+a_{4}x^{2}+a_{5}y^{2}+a_{6}xy\right ) =\left ( a_{2}+2a_{4}x+a_{6}y\right ) ^{2}+\left ( a_{3}+2a_{5}y+a_{6}x\right ) ^{2}$$ Expanding and comparing coefficients gives

$$\begin{multline*} 2x^{2}a_{4}^{2}+2x^{2}a_{4}a_{5}+2a_{6}a_{4}xy+2a_{6}a_{5}xy+2a_{2}a_{4}x+2a_{2}xa_{5}+2y^{2}a_{4}a_{5}+2y^{2}a_{5}^{2}+2a_{3}a_{4}y+2a_{3}a_{5}y+2a_{1}a_{4}+2a_{1}a_{5}=\\ 4x^{2}a_{4}^{2}+x^{2}a_{6}^{2}+4a_{4}a_{6}xy+4a_{5}a_{6}xy+4xa_{2}a_{4}+2a_{3}a_{6}x+4y^{2}a_{5}^{2}+y^{2}a_{6}^{2}+2a_{2}a_{6}y+4a_{3}a_{5}y+a_{2}^{2}+a_{3}^{2} \end{multline*}$$ Simplifying $$\begin{multline*} 2a_{4}a_{5}x^{2}+2a_{2}a_{5}x+2a_{4}a_{5}y^{2}+2a_{3}a_{4}y+2a_{1}a_{4}+2a_{1}a_{5}=\\ 2x^{2}a_{4}^{2}+a_{6}^{2}x^{2}+2a_{4}a_{6}xy+2a_{5}a_{6}xy+2a_{2}a_{4}x+2a_{3}a_{6}x+2a_{5}^{2}y^{2}+a_{6}^{2}y^{2}+2a_{2}a_{6}y+2a_{3}a_{5}y+a_{2}^{2}+a_{3}^{2} \end{multline*}$$ Comparing coefficients of terms that contain no $x,y$ and coefficients of $x,y,xy,x^{2},y^{2}$ gives the following equations in order $$\begin{align*} 2a_{1}a_{4}+2a_{1}a_{5} & =a_{2}^{2}+a_{3}^{2}\\ 2a_{2}a_{5} & =2a_{2}a_{4}+2a_{3}a_{6}\\ 2a_{3}a_{4} & =2a_{2}a_{6}+2a_{3}a_{5}\\ 0 & =4a_{4}a_{6}\\ 2a_{4}a_{5} & =2a_{4}^{2}+a_{6}^{2}\\ 2a_{4}a_{5} & =2a_{5}^{2}+a_{6}^{2} \end{align*}$$

Equation $0=4a_{4}a_{6}$ above implies that $a_{4}=0$ or $a_{6}=0$ or both are zero. But if both are zero, there is no solution. On the other hand, if $a_{4}=0,$ then this also leads to no solution as all equations reduce to $0=0$. Therefore only choice left is $a_{6}=0$. Now the above equations become $$\begin{align*} 2a_{1}a_{4}+2a_{1}a_{5} & =a_{2}^{2}+a_{3}^{2}\\ 2a_{2}a_{5} & =2a_{2}a_{4}\\ 2a_{3}a_{4} & =2a_{3}a_{5}\\ 0 & =0\\ 2a_{4}a_{5} & =2a_{4}^{2}\\ 2a_{4}a_{5} & =2a_{5}^{2} \end{align*}$$

Or $$\begin{align*} 2a_{1}a_{4}+2a_{1}a_{5} & =a_{2}^{2}+a_{3}^{2}\\ a_{5} & =a_{4}\\ a_{4} & =a_{5}\\ 0 & =0\\ a_{5} & =a_{4}\\ a_{4} & =a_{5} \end{align*}$$

Hence $$\begin{align} a_{4} & =a_{5}\tag{3}\\ a_{6} & =0\tag{4}\\ 2a_{1}a_{4}+2a_{1}a_{5} & =a_{2}^{2}+a_{3}^{2}\nonumber \end{align}$$

Since $a_{4}=a_{5}$ then $$\begin{align} 2a_{1}a_{5}+2a_{1}a_{5} & =a_{2}^{2}+a_{3}^{2}\nonumber \\ a_{5} & =\frac{a_{2}^{2}+a_{3}^{2}}{2a_{1}} \tag{5} \end{align}$$

Using (3,4,5) in $p\left ( x,y\right ) =a_{1}+a_{2}x+a_{3}y+a_{4}x^{2}+a_{5}y^{2}+a_{6}xy$ gives $$\begin{align*} p\left ( x,y\right ) & =a_{1}+a_{2}x+a_{3}y+a_{5}x^{2}+a_{5}y^{2}\\ & =a_{1}+a_{2}x+a_{3}y+a_{5}\left ( x^{2}+y^{2}\right ) \\ & =a_{1}+a_{2}x+a_{3}y+\frac{a_{2}^{2}+a_{3}^{2}}{2a_{1}}\left ( x^{2}+y^{2}\right ) \end{align*}$$

Only three arbitrary constants are needed. Let $a_{1}=a,a_{2}=b,a_{3}=c$ the above becomes $$p\left ( x,y\right ) =a+bx+cy+\frac{b^{2}+c^{2}}{2a}\left ( x^{2}+y^{2}\right )$$ And the solution becomes $$u\left ( x,y\right ) =\log \left ( a+bx+cy+\frac{b^{2}+c^{2}}{2a}\left ( x^{2}+y^{2}\right ) \right )$$

2.1.3 Problem 1.13

Find all solutions $u=f\left ( r\right )$ of the 3D Laplace equation $u_{xx}+u_{yy}+u_{zz}=0$ that depends only on radial coordinates $r=\sqrt{x^{2}+y^{2}+z^{2}}$

Solution

The Laplacian in 3D in spherical coordinates is $$\nabla ^{2}u\left ( r,\theta ,\phi \right ) =u_{rr}+\frac{2}{r}u_{r}+\frac{1}{r^{2}}\left ( \frac{\cos \theta }{\sin \theta }u_{\theta }+u_{\theta \theta }\right ) +\frac{1}{r^{2}\sin ^{2}\theta }u_{\phi \phi }$$ The above shows that the terms that depend only on $r$ makes the laplacian $$\nabla ^{2}u\left ( r\right ) =u_{rr}+\frac{2}{r}u_{r}$$ Hence the PDE $\nabla ^{2}u\left ( r\right ) =0$ becomes an ODE now since there is only one dependent variable giving $$u^{\prime \prime }\left ( r\right ) +\frac{2}{r}u^{\prime }\left ( r\right ) =0$$ Let $v=u^{\prime }\left ( r\right )$ and the above becomes $$v^{\prime }\left ( r\right ) +\frac{2}{r}v\left ( r\right ) =0$$ This is linear first order ODE. The integrating factor is $I=e^{\int \frac{2}{r}dr}=e^{2\ln r}=r^{2}$. Therefore the above becomes $\frac{d}{dr}\left ( vr^{2}\right ) =0$ or $vr^{2}=C_{1}$ or $v\left ( r\right ) =\frac{C_{1}}{r^{2}}$. Therefore $$\begin{align*} u^{\prime } & =\frac{C_{1}}{r^{2}}\\ du & =\frac{C_{1}}{r^{2}}dr \end{align*}$$

Integrating gives the solution $$u=-\frac{C_{1}}{r}+C_{2}$$ The above is the required solution. Hence $$\fbox{$f\left ( r\right ) =-\frac{C_1}{r}+C_2$}$$ Where $C_{1},C_{2}$ are arbitrary constants.

2.1.4 Problem 1.20

The displacement $u\left ( t,x\right )$ of a forced violin string is modeled by the PDE $u_{tt}=4u_{xx}+F\left ( t,x\right )$. When the string is subjected to the external force $F\left ( t,x\right ) =\cos x$, the solution is $u\left ( t,x\right ) =\cos \left ( x-2t\right ) +\frac{1}{4}\cos x$, while when $F\left ( t,x\right ) =\sin x$, the solution is $u\left ( t,x\right ) =\sin \left ( x-2t\right ) +\frac{1}{4}\sin x$. Find a solution when the forcing function is (a) $\cos x-5\sin x$, (b) $\sin \left ( x-3\right )$

Solution

2.1.4.1 Part (a)

Since the PDE is linear, superposition can be used. When the input is $F\left ( t,x\right ) =\cos x-5\sin x$ then the solution is $$\begin{align*} u\left ( t,x\right ) & =\left ( \cos \left ( x-2t\right ) +\frac{1}{4}\cos x\right ) -5\left ( \sin \left ( x-2t\right ) +\frac{1}{4}\sin x\right ) \\ & =\cos \left ( x-2t\right ) +\frac{1}{4}\cos x-5\sin \left ( x-2t\right ) -\frac{5}{4}\sin x \end{align*}$$

2.1.4.2 Part (b)

Since the PDE is linear, superposition can be used. When the input is $F\left ( t,x\right ) =\sin \left ( x-3\right )$ then the solution same as when the input is $\sin x$ but shifted by $3$. Hence $$u\left ( t,x\right ) =\sin \left ( \left ( x-3\right ) -2t\right ) +\frac{1}{4}\sin \left ( x-3\right )$$

2.1.5 Problem 1.27b

Solve the following inhomogeneous linear ODE $5u^{\prime \prime }-4u^{\prime }+4u=e^{x}\cos x$

Solution

First the homogeneous solution $u_{h}$ is found, then a particular solution $u_{p}$ is found. The general solution will be the sum of both $u=u_{h}+u_{p}$. Since this is a constant coefficient ODE, the characteristic equation is $5\lambda ^{2}-4\lambda +4=0$. The roots are $\lambda _{1}=\frac{2}{5}+\frac{4}{5}i,\lambda _{1}=\frac{2}{5}-\frac{4}{5}i$, which implies the solution is $$u_{h}\left ( x\right ) =e^{\frac{2}{5}x}\left ( c_{1}\cos \left ( \frac{4}{5}x\right ) +c_{2}\sin \left ( \frac{4}{5}x\right ) \right )$$ Using the method of undetermined coefficients, and since the forcing function is $e^{x}\cos x$, then let $$\begin{equation} u_{p}=Ae^{x}\left ( B\cos x+C\sin x\right ) \tag{1} \end{equation}$$ Hence $$\begin{align} u_{p}^{\prime } & =Ae^{x}\left ( B\cos x+C\sin x\right ) +Ae^{x}\left ( -B\sin x+C\cos x\right ) \tag{2}\\ u_{p}^{\prime \prime } & =Ae^{x}\left ( B\cos x+C\sin x\right ) +Ae^{x}\left ( -B\sin x+C\cos x\right ) +Ae^{x}\left ( -B\sin x+C\cos x\right ) +Ae^{x}\left ( -B\cos x-C\sin x\right ) \nonumber \\ & =Ae^{x}\left ( B\cos x+C\sin x-B\sin x+C\cos x-B\sin x+C\cos x-B\cos x-C\sin x\right ) \nonumber \\ & =Ae^{x}\left ( -B\sin x+C\cos x-B\sin x+C\cos x\right ) \nonumber \\ & =Ae^{x}\left ( -2B\sin x+2C\cos x\right ) \tag{3} \end{align}$$

Substituting (1,2,3) back into the original ODE gives $$\begin{align*} 5Ae^{x}\left ( -2B\sin x+2C\cos x\right ) -4\left ( Ae^{x}\left ( B\cos x+C\sin x\right ) +Ae^{x}\left ( -B\sin x+C\cos x\right ) \right ) +4Ae^{x}\left ( B\cos x+C\sin x\right ) & =e^{x}\cos x\\ Ae^{x}\left ( -10B\sin x+10C\cos x\right ) -Ae^{x}\left ( 4B\cos x+4C\sin x\right ) -Ae^{x}\left ( -4B\sin x+4C\cos x\right ) +Ae^{x}\left ( 4B\cos x+4C\sin x\right ) & =e^{x}\cos x\\ Ae^{x}\left ( -10B\sin x+10C\cos x-4B\cos x-4C\sin x+4B\sin x-4C\cos x+4B\cos x+4C\sin x\right ) & =e^{x}\cos x \end{align*}$$

Hence $$Ae^{x}\left ( 6C\cos x-6B\sin x\right ) =e^{x}\cos x$$ Comparing coefficients shows that $$\begin{align*} A & =1\\ B & =0\\ C & =\frac{1}{6} \end{align*}$$

Hence from (1) $$u_{p}=e^{x}\frac{\sin x}{6}$$ Therefore the general solution is $$\begin{align*} u\left ( x\right ) & =u_{h}\left ( x\right ) +u_{p}\left ( x\right ) \\ & =e^{\frac{2}{5}x}\left ( c_{1}\cos \left ( \frac{4}{5}x\right ) +c_{2}\sin \left ( \frac{4}{5}x\right ) \right ) +e^{x}\frac{\sin x}{6} \end{align*}$$

2.1.6 Problem 2.1.6

Solve the PDE $\frac{\partial ^{2}u}{\partial x\partial y}=0$ for $u\left ( x,y\right )$

Solution

Integrating once w.r.t $x$ gives $$\frac{\partial u}{\partial y}=F\left ( y\right )$$ Where $F\left ( y\right )$ acts as the constant of integration, but since this is a PDE, it becomes an arbitrary function of $y$ only. Integrating the above again w.r.t. $y$ gives $$u=\int F\left ( y\right ) dy+G\left ( x\right )$$ Where $G\left ( x\right )$ is an arbitrary function of $x$ only. If we let $\int F\left ( y\right ) dy=H\left ( y\right )$ where $H\left ( y\right )$ is the antiderivative for the indefinite integral which depends on $y$ only. Then the above can be written as $$\fbox{$u\left ( x,y\right ) =H\left ( y\right ) +G\left ( x\right ) $}$$ To verify, from the above $\frac{\partial u}{\partial y}=H^{\prime }\left ( y\right )$ and hence $$\begin{align*} \frac{\partial ^{2}u}{\partial x\partial y} & =\frac{d}{dx}\left ( H^{\prime }\left ( y\right ) \right ) \\ & =0 \end{align*}$$

2.1.7 Problem 2.2.2

Solve the following initial value problems and graph the solutions at $t=1,2,3$

a

$u_{t}-3u_{x}=0,u\left ( 0,x\right ) =e^{-x^{2}}$

b

$u_{t}+2u_{x}=0,u\left ( -1,x\right ) =\frac{x}{1+x^{2}}$

c

$u_{t}+u_{x}+\frac{1}{2}u=0,u\left ( 0,x\right ) =\arctan \left ( x\right )$

d

$u_{t}-4u_{x}+u=0,u\left ( 0,x\right ) =\frac{1}{1+x^{2}}$

Solution

2.1.7.1 Part a

Let $\xi$ be the characteristic variable defined such that $\xi =x-ct$. Where characteristic lines are given by $x=x_{0}+ct$. But $c=-3$ in this problem. Hence characteristic lines are $$x=x_{0}-3t$$ Where $x_{0}$ means the same as $x\left ( 0\right )$, i.e. $x\left ( t\right )$ at time $t=0$. Since $c=-3$ then $$\xi =x+3t$$ Let $$u\left ( t,x\right ) \equiv v\left ( t,\xi \right )$$ $u_{t}-3u_{x}=0$ is now transformed to $v\left ( t,\xi \right )$ as follows $$\begin{align} \frac{\partial u}{\partial t} & =\frac{\partial v}{\partial t}\frac{\partial t}{\partial t}+\frac{\partial v}{\partial \xi }\frac{\partial \xi }{\partial t}\nonumber \\ & =\frac{\partial v}{\partial t}+3\frac{\partial v}{\partial \xi } \tag{1} \end{align}$$

And $$\begin{align} \frac{\partial u}{\partial x} & =\frac{\partial v}{\partial t}\frac{\partial t}{\partial x}+\frac{\partial v}{\partial \xi }\frac{\partial \xi }{\partial x}\nonumber \\ & =0+\frac{\partial v}{\partial \xi }\nonumber \\ & =\frac{\partial v}{\partial \xi } \tag{2} \end{align}$$

Substituting (1,2) in $u_{t}-3u_{x}=0$ gives the transformed PDE as $$\begin{align*} \frac{\partial v}{\partial t}+3\frac{\partial v}{\partial \xi }-3\frac{\partial v}{\partial \xi } & =0\\ \frac{\partial v}{\partial t} & =0 \end{align*}$$

Integrating w.r.t $\xi$ gives the solution in $v\left ( t,\xi \right )$ space as $$v\left ( t,\xi \right ) =F\left ( \xi \right )$$ Where $F\left ( \xi \right )$ is an arbitrary continuous function of $\xi$. Transforming back to $u\left ( t,x\right )$ gives $$\begin{equation} u\left ( t,x\right ) =F\left ( x+3t\right ) \tag{3} \end{equation}$$ At $t=0$ the above becomes $$e^{-x_{0}^{2}}=F\left ( x_{0}\right )$$ This means that (3) becomes (since $x=x_{0}+ct$ or $x=x_{0}-3t$ or $x_{0}=x+3t$) $$u\left ( t,x\right ) =e^{-\left ( x+3t\right ) ^{2}}$$

2.1.7.2 Part b

$$\begin{align*} u_{t}+2u_{x} & =0\\ u\left ( -1,x\right ) & =\frac{x}{1+x^{2}} \end{align*}$$

Let $\xi$ be the characteristic variable defined such that $\xi =x-ct$. Where characteristic lines are given by $x=x_{0}+ct$. But $c=2$ in this problem. Hence characteristic lines are $$x=x_{0}+2t$$ And $$\xi =x-2t$$ Let $u\left ( t,x\right ) \equiv v\left ( t,\xi \right )$. Then $u_{t}+2u_{x}=0$ is transformed to $v\left ( t,\xi \right )$ as was done in part (a) (will not be repeated) which results in $$\frac{\partial v}{\partial t}=0$$ Integrating w.r.t $\xi$ gives the solution $$v\left ( t,\xi \right ) =F\left ( \xi \right )$$ Where $F\left ( \xi \right )$ is an arbitrary continuous function of $\xi$. Transforming back to $u\left ( t,x\right )$ results in $$\begin{equation} u\left ( t,x\right ) =F\left ( x-2t\right ) \tag{3} \end{equation}$$ At $t=-1$ the above becomes $$\frac{x_{0}}{1+x_{0}^{2}}=F\left ( x_{0}+2\right )$$ Let $x_{0}+2=z$. Then $x_{0}=z-2$. And the above becomes $$\frac{z-2}{1+\left ( z-2\right ) ^{2}}=F\left ( z\right )$$ This means that (3) becomes $$\begin{align*} u\left ( t,x\right ) & =\frac{\left ( x-2t\right ) -2}{1+\left ( \left ( x-2t\right ) -2\right ) ^{2}}\\ & =\frac{x-2t-2}{1+\left ( x-2t-2\right ) ^{2}} \end{align*}$$

2.1.7.3 Part c

$$\begin{align} u_{t}+u_{x}+\frac{1}{2}u & =0\tag{1}\\ u\left ( 0,x\right ) & =\arctan \left ( x\right ) \nonumber \end{align}$$

Let $\xi$ be the characteristic variable defined such that $\xi =x-ct$. Where characteristic lines are given by $x=x_{0}+ct$. But $c=1$ in this problem. Hence characteristic lines are given by solution to $$\begin{align*} \frac{dx}{dt} & =1\\ x\left ( t\right ) & =x_{0}+t \end{align*}$$

And $$\begin{align*} \xi & =x-ct\\ & =x-t \end{align*}$$

Then $u_{t}+u_{x}$ are transformed to $v\left ( t,\xi \right )$ as was done in part (a) (will not be repeated) which results in $$u_{t}+u_{x}=\frac{\partial v}{\partial t}$$ Substituting the above into (1) gives (where now $v$ is used in place of $u$). $$\frac{\partial v}{\partial t}+\frac{1}{2}v=0$$ This is now first order ODE since it only depends on $t$. Therefore $v^{\prime }+\frac{1}{2}v=0$. This is linear in $v$. Hence the solution is $\frac{d}{dt}\left ( ve^{\int \frac{1}{2}dt}\right ) =0$ or $ve^{\frac{1}{2}t}=F\left ( \xi \right )$ where $F$ is arbitrary function of $\xi$. Hence $$v\left ( t,\xi \right ) =e^{\frac{-1}{2}t}F\left ( \xi \right )$$ Converting back to $u\left ( t,x\right )$ gives $$\begin{equation} u\left ( t,x\right ) =e^{\frac{-t}{2}}F\left ( x-t\right ) \tag{2} \end{equation}$$ At $t=0$ the above becomes $$\arctan \left ( x_{0}\right ) =F\left ( x_{0}\right )$$ From the above then (2) can be written as $$u\left ( t,x\right ) =e^{\frac{-t}{2}}\arctan \left ( x-t\right )$$

2.1.7.4 Part d

$$\begin{align*} u_{t}-4u_{x}+u & =0\\ u\left ( 0,x\right ) & =\frac{1}{1+x^{2}} \end{align*}$$

Let $\xi$ be the characteristic variable defined such that $\xi =x-ct$. Where characteristic lines are given by $x=x_{0}+ct$. But $c=-4$ in this problem. Hence characteristic lines are $$x=x_{0}-4t$$ And $$\xi =x+4t$$ Then $u_{t}-4u_{x}$ are transformed to $v\left ( t,\xi \right )$ as was done in part (a) (will not be repeated) which results in $$u_{t}-4u_{x}=\frac{\partial v}{\partial t}$$ Substituting the above into (1) gives (where now $v$ is used in place of $u$). $$\frac{\partial v}{\partial t}+v=0$$ This is now first order ODE since it only depends on $t$. Therefore $v^{\prime }+v=0$. This is linear in $v$. Hence the solution is $\frac{d}{dt}\left ( ve^{\int dt}\right ) =0$ or $ve^{t}=F\left ( \xi \right )$ where $F$ is arbitrary function of $\xi$. Hence $$v\left ( t,\xi \right ) =e^{-t}F\left ( \xi \right )$$ Converting to $u\left ( t,x\right )$ gives $$\begin{equation} u\left ( t,x\right ) =e^{-t}F\left ( x+4t\right ) \tag{2} \end{equation}$$ At $u\left ( 0,x\right ) =\frac{1}{1+x^{2}}$ the above becomes $$\frac{1}{1+x_{0}^{2}}=F\left ( x_{0}\right )$$ From the above then (2) can be written as $$u\left ( t,x\right ) =\frac{e^{-t}}{1+\left ( x+4t\right ) ^{2}}$$

2.1.8 Problem 2.2.3

Graph some of the characteristic lines for the following equation and write down the formula for the general solution

(b) $u_{t}+5u_{x}=0\,$, (d) $u_{t}-4u_{x}+u=0$

Solution

2.1.8.1 Part b

$$u_{t}+5u_{x}=0$$

Let $\xi$ be the characteristic variable defined such that $\xi =x-ct$. Where characteristic lines are given by $x=x_{0}+ct$. But $c=5$ in this problem. Hence characteristic lines are $$\begin{equation} \fbox{$x\left ( t\right ) =x_0+5t$} \tag{1} \end{equation}$$ And $$\xi =x-5t$$ Then $u_{t}-5u_{x}=0$ is transformed to $v\left ( t,\xi \right )$ as was done in earlier (will not be repeated) which results in $$u_{t}-5u_{x}=\frac{\partial v}{\partial t}$$ Therefore $\frac{\partial v}{\partial t}=0$ which has the general solution $v\left ( t,\xi \right ) =F\left ( \xi \right )$ where $F$ is arbitrary function of $\xi$. Transforming back to $u\left ( t,x\right )$ gives $$\fbox{$u\left ( t,x\right ) =F\left ( x-5t\right ) $}$$ On the characteristic lines given by (1) the solution $u\left ( t,x\right )$ is constant. The slope of the characteristic lines is $5$ and intercept is $x_{0}$. The following is a plot of few lines using different values of $x_{0}$.

Figure 2.1:Showing some characteristic lines for part b

2.1.8.2 Part d

$$u_{t}-4u_{x}+u=0$$ Let $\xi$ be the characteristic variable defined such that $\xi =x-ct$. Where characteristic lines are given by $x=x_{0}+ct$. But $c=-4$ in this problem. Hence characteristic lines are $$\begin{equation} \fbox{$x\left ( t\right ) =x_0-4t$} \tag{1} \end{equation}$$ And $$\xi =x+4t$$ Then $u_{t}-4u_{x}$ is transformed to $v\left ( t,\xi \right )$ as was done in earlier (will not be repeated) which results in $$u_{t}-4u_{x}=\frac{\partial v}{\partial t}$$ Therefore the original PDE becomes $\frac{\partial v}{\partial t}+v=0$, where $u$ is replaced by $v$. This is linear first order ODE which has the solution $v\left ( t,\xi \right ) =e^{-t}F\left ( \xi \right )$ where $F$ is arbitrary function of $\xi$. Transforming back to $u\left ( t,x\right )$ gives the general solution as $$u\left ( t,x\right ) =e^{-t}F\left ( x+4t\right )$$ The following is a plot of few characteristic lines $x=x_{0}-4t$ using different values of $x_{0}$.

Figure 2.2:Showing some characteristic lines for part d

2.1.9 Problem 2.2.5

Solve $u_{t}+2u_{x}=\sin x,u\left ( 0,x\right ) =\sin x$

Solution

Let $\xi$ be the characteristic variable defined such that $\xi =x-ct$. Where characteristic lines are given by $x=x_{0}+ct$. But $c=2$ in this problem. Hence characteristic lines are $$\begin{equation} \fbox{$x=x_0+2t$} \tag{1} \end{equation}$$ And $$\xi =x-2t$$ Then $u_{t}+2u_{x}$ is transformed to $v\left ( t,\xi \right )$ as was done in earlier (will not be repeated) which results in $$u_{t}+2u_{x}=\frac{\partial v}{\partial t}$$ Substituting this into the original PDE gives $$\frac{\partial v\left ( t,\xi \right ) }{\partial t}=\sin \left ( \xi +2t\right )$$ Integrating w.r.t $t$ gives $$\begin{align*} v\left ( t,\xi \right ) & =\int \sin \left ( \xi +2t\right ) dt+F\left ( \xi \right ) \\ & =-\frac{\cos \left ( \xi +2t\right ) }{2}+F\left ( \xi \right ) \end{align*}$$

Transforming back to $u\left ( t,x\right )$ gives $$\begin{align} u\left ( t,x\right ) & =-\frac{\cos \left ( x-2t+2t\right ) }{2}+F\left ( x-2t\right ) \nonumber \\ & =\frac{-1}{2}\cos \left ( x\right ) +F\left ( x-2t\right ) \tag{1} \end{align}$$

When $t=0$, $u\left ( 0,x\right ) =\sin x$, therefore the above becomes $$\begin{align*} \sin x_{0} & =F\left ( x_{0}\right ) -\frac{1}{2}\cos x_{0}\\ F\left ( x_{0}\right ) & =\sin x_{0}+\frac{1}{2}\cos x_{0} \end{align*}$$

Therefore the solution (1) becomes $$\begin{align*} u\left ( t,x\right ) & =\left ( \sin \left ( x-2t\right ) +\frac{1}{2}\cos \left ( x-2t\right ) \right ) -\frac{1}{2}\cos x\\ & =\sin \left ( x-2t\right ) +\frac{1}{2}\cos \left ( x-2t\right ) -\frac{1}{2}\cos x \end{align*}$$

2.1.10 Problem 2.2.9

(a) Prove that if the initial data is bounded, $\left \vert f\left ( x\right ) \right \vert \leq M$ for all $x\in \mathbb{R}$, then the solution to the damped transport equation (2.14) $u_{t}+cu_{x}+au=0$ with $a>0$ satisfies $u\left ( t,x\right ) \rightarrow 0$ as $t\rightarrow \infty$. (b) Find a solution to (2.14) that is defined for all $\left ( t,x\right )$ but does not satisfy $u\left ( t,x\right ) \rightarrow 0$ as $t\rightarrow \infty$.

Solution

2.1.10.1 Part(a)

$u_{t}+cu_{x}+au=0$ is solved to show what is required. Let $\xi$ be the characteristic variable defined such that $\xi =x-ct$. Where characteristic lines are given by $x=x_{0}+ct$. Hence characteristic lines are $$\begin{equation} x=x_{0}+ct \tag{1} \end{equation}$$ And $$\xi =x-ct$$ Then $u_{t}+cu_{x}$ is transformed to $v\left ( t,\xi \right )$ as was done in earlier (will not be repeated) which results in $$u_{t}+cu_{x}=\frac{\partial v}{\partial t}$$ Substituting this into the original PDE gives $$\frac{\partial v}{\partial t}+av=0$$ Where $u$ is replaced by $v$. This can be viewed as first order linear ODE since it depends on $t$ only. Its solution is $v\left ( t,\xi \right ) =e^{-at}F\left ( \xi \right ) \,$ where $F$ is arbitrary function of $\xi$. Transforming back to $u\left ( t,x\right )$ gives $$\begin{equation} u\left ( t,x\right ) =e^{-at}F\left ( x-ct\right ) \tag{1} \end{equation}$$ At $t=0$ initial data is $f\left ( x\right )$. Hence the above becomes at $t=0$ $$f\left ( x\right ) =F\left ( x\right )$$ Hence (1) now becomes $$\begin{equation} u\left ( t,x\right ) =e^{-at}f\left ( x-ct\right ) \tag{2} \end{equation}$$ But since $\left \vert f\left ( x\right ) \right \vert$ is bounded, and since $a>0$ then $e^{-at}\rightarrow 0$ as $t\rightarrow \infty$. Which implies the solution itself $u\left ( t,x\right )$ goes to zero as well. This is the reason why initial data needed to be bounded for this to happen.

2.1.10.2 Part(b)

Keeping $a>0$. If initial data have the form $f\left ( x\right ) e^{-bx}$ where $\left \vert b\right \vert >a$, then at $t=0$ the solution found in (1) becomes $$f\left ( x_{0}\right ) e^{-bx_{0}}=F\left ( x_{0}\right )$$ Then the solution (2) now becomes, after replacing $x_{0}$ by $x-ct$ $$\begin{align*} u\left ( t,x\right ) & =e^{-at}e^{-b\left ( x-ct\right ) }f\left ( x-ct\right ) \\ & =e^{-at+bct}e^{-bx}f\left ( x-ct\right ) \\ & =e^{\left ( bc-a\right ) t}e^{-bx}f\left ( x-ct\right ) \end{align*}$$

The problem is asking to show that this does not go to zero for all $x\in \mathbb{R}$ as $t\rightarrow \infty$. Since $\left \vert b\right \vert >a$ then $bc-a$ is positive quantity ($c$ is assumed positive)¹ .

Therefore $e^{\left ( bc-a\right ) t}$ will blow up as $t\rightarrow \infty$. And therefore the whole solution will not go to zero. For any $x$, no matter how large $x$ is, a large enough $t$ can be found to make the product $e^{\left ( bc-a\right ) t}e^{-bx}$ blow up.

2.1.11 Key solution for HW 1

PDF