% This is a sample file for the Undergraduate Faculty Program at
% PCMI, containing one lecture 
%
%   What is a Partial Differential Equation?
%
%  from the Park City Lectures of Andrew J. Bernoff.
%
% The sample file illustrates the use of epsf.tex to include postscript 
% graphics, as well as various "AmS-LaTeX" constructions from the amsmath
% package (which is automatically loaded by the pcms-l document class).
% To run this file you need the files 
%
%        pcms-l.cls and pcmslmod.tex
%


\documentclass[lecture,12pt,]{pcms-l}
\input pcmslmod.tex  % v.1.2
%\input epsf.tex

%\usepackage{amssymb} % this command would have loaded all the extra symbols,



% authors should not define these, they will be defined by the volume editors
%\def\currentvolume{3}
%\def\currentyear{1993}


% EQUATION NUMBERING AND THEOREM SETUP

\numberwithin{section}{chapter}
\numberwithin{equation}{chapter}

\theoremstyle{plain}
\newtheorem{theorem}[equation]{Theorem}
\newtheorem{lemma}[equation]{Lemma}

\theoremstyle{definition}
\newtheorem{definition}[equation]{Definition}

\newtheorem{exercise}{Exercise}
\newtheorem{problem}{Problem}
\newtheorem*{remark}{Remark}

% Set enumerate to use letters, not numbers for problem parts.

\renewcommand{\theenumi}{\alph{enumi}}
\renewcommand{\labelenumi}{(\theenumi)}

% AUTHOR-DEFINED MACROS:

\newcommand{\cA}{\mathcal{A}}
\newcommand{\cB}{\mathcal{B}}
\newcommand{\C}{\mathbb{C}}
\newcommand{\cC}{\mathcal{C}}
\newcommand{\cD}{\mathcal{D}}
\newcommand{\cE}{\mathcal{E}}
\newcommand{\bU}{{\bar U}}

\begin{document}

\mainmatter
\setcounter{page}{1}

%\LogoOn

\lectureseries[An Introduction to PDE]{An Introduction to Partial Differential Equations}


\auth[Andrew J. Bernoff]{Andrew J. Bernoff} 

%\address{Department of Mathematics, Harvey Mudd College,Claremont, CA 91711}
\email{ajb@hmc.edu}

%The following items will become first page footnotes; they are optional.
%\date{\today}
%\setaddress


% the following hack starts the lecture numbering at 2
\setcounter{lecture}{1}
\setcounter{chapter}{1}

\lecture{Cooling of a Hot Bar:\break The Diffusion Equation}


\section{Outline of Lecture}
\begin{itemize}
\item An Introduction to Heat Flow
\item Derivation of the Diffusion Equation
\item Examples of Solution to the Diffusion Equation
\item The Maximum Principle
\item Energy Dissipation and Uniqueness
\end{itemize}


\section{An Introduction to Heat Flow}

A classical example of the application of ordinary differential equations
is Newton's Law of Cooling which, basically, answers the question {\it ``How does a cup of coffee cool?''}  Newton hypothesized that the rate at which the temperature, $U(t)$,
changes is proportional to the difference with the ambient temperature, which we call $\bU$,
\begin{equation}
\frac{dU}{dt} = - \kappa (U- \bU) .
\end{equation}
Here $\kappa$ is a positive rate constant (with units of inverse time) that measures how fast heat is lost from the coffee cup to the ambient environment. If we specify 
the initial temperature,
\begin{equation}
U(0) =U_0,
\end{equation}
we can solve for the evolution of the temperature,
\begin{equation}
U(t) = \bU +(U_0 -\bU) e^{-\kappa t}.
\end{equation}


\vskip 3 in

\begin{center}
Figure 2.1: (a) A coffee cup  (b) Its temperature as a function of time.\\
 (draw your own figure).
\end{center}
If we graph the temperature as a function of time, we see that it decays exponentially to the ambient temperature, $\bU$, at a rate governed by $\kappa$. 

When we derived Newton's Law of cooling we made several assumptions -- most importantly that the temperature in the coffee cup did not vary with location. If we account for the variation of temperature with location, we can derive a PDE called the {\bf heat equation} or,  more generally, the  {\bf diffusion  equation}. If the temperature, $U(x,t)$ is a function of a single spatial variable, $x$, we will show that it satisfies the diffusion equation,
\begin{equation*}
U_t = DU_{xx},
\end{equation*}
where $D$ is a constant known as the thermal diffusivity. In higher dimensions, the equation can be written 
\begin{equation*}
U_t = D \nabla^2 U ,
\end{equation*}
where $\nabla^2$ is the {\bf Laplacian}.

\section{Derivation of the Diffusion Equation}

The diffusion equation will be our second example of a conservation law; we can derive the equation by accounting for the flow of thermal energy. Suppose we consider a metal bar, with a uniform cross-sectional area, $A$, whose temperature, $U(x,t)$, is a function  of time, $t$, and the position, $x$, along the bar (that is we assume the temperature is uniform in every cross-section).



\vskip 3 in

\begin{center}
Figure 2.2: Conservation of heat in a metal bar\\ of cross-sectional area $A$.\\
 (draw your own figure).
\end{center}

Let the thermal energy in the region $a<x<b$ is given by
\begin{equation}
\cE =   \rho_0 c_v A  \int_a^b U(x,t) ~ dx
\end{equation}
The important term in the integral is the temperature, $U(x,t)$, measured in degrees. The remaining constants, $A$, the cross-sectional area (with units of [(length)$^2$]), $\rho_0$, the density [mass/(length)$^3$] and $c_v$, the heat capacity [energy/(degree $\cdot$ mass)] are physical properties of the material -- think of them as being obligatory for making the units work out.

We wish to equate the change in thermal energy to the heat flux out of the bar through the planes at $x=a$ and $x=b$. To do this we use {\bf Fourier's heat law} 
which states that the flux density  of thermal energy, $q(x,t)$ is proportional to the temperature gradient,
\begin{equation}
q(x,t) =  - k U_x  ,
\end{equation}
where the negative sign reflects the fact that heat flows from hot to cold,
just as in Newton's law of cooling, with a constant of proportionality, $k$, called the thermal conductivity [(energy$\cdot$length)/(degrees$\cdot$time)].
\eject
\phantom{redrum}


\vskip 3 in

\begin{center}
Figure 2.3: Heat flux is  from hot to cold!!\\
 (draw your own figure).
\end{center}

Now, the total flux of thermal energy into the 
{\bf into} the region $a<x<b$ is given by
\begin{equation}
Q = A[q(a,t)-q(b,t)],
\end{equation}
where we multiply by the area $A$ to get the total flux through the cross-section.

By {\bf conservation of energy}, the rate of change of the energy between $a$ and $b$ is given by the flux into the region,
\begin{equation}
\frac{d\cE}{dt} = Q .
\label{CofE}
\end{equation}
Once again we can rewrite the flux by a clever application of the fundamental theorem of calculus,
\begin{eqnarray}
Q &=& A[ q(a,t)-q(b,t)] = - \left. Aq(x,t) \right  |_{x=a}^{x=b} \\
    &=&- A\int_a^b q_x ~dx .
\end{eqnarray} 
We now rewrite the conservation of energy equation as  
\begin{equation}
\frac{d\cE}{dt}  = \frac{d}{dt} \left [ \rho c_v A \int_a^b U ~ dx \right ]    
                        =  \int_a^b  \rho c_v A U_t  ~ dx    
                        =  Q  =- A\int_a^b ~ q_x dx ,
\label{CofM2}
\end{equation}
or, rearranging
\begin{equation}
  \int_a^b \rho c_v A ~U_t   +Aq_x ~ dx =0 .   
\end{equation}
Since this is true for {\bf every} interval $a<x<b$, the integrand must vanish identically. So
\begin{equation}
 \rho c_v A  U_t   +Aq_x  =0 .   
\end{equation}
Substituting for the flux function $q(x,t) = -kU_x$ yields
\begin{equation}
 \rho c_v A  U_t   -kA(U_x)_x  =0 .   
\end{equation}
Rearranging the equation yields the diffusion equation,
 \begin{equation}
\boxed{ U_t   = DU_{xx}    ,}
\end{equation}
where the {\bf diffusivity}, $D = k/(\rho c_v)$, is a constant which
is determined by the geometry and physical properties of the metal
bar.

To complete the description of the problem, we need to supplement the
diffusion equation with boundary conditions and initial conditions.
Suppose we consider a bar of finite length $L$, occupying the region
$0<x<L$. At the boundaries of the metal bar we can specify 
a fixed temperature,
\begin{equation}
U(0,t) = U_0 \qquad U(L,t)=U_1 ,
\end{equation}
which are usually referred to as {\bf Dirichlet} boundary conditions.
Alternatively, we could specify a heat flux,
\begin{equation}
q_0= q(0,t) = -kU_x(0,t) \qquad q_1= q(L,t) = -kU_x(L,t) .
\end{equation}
Specifying the gradient across the boundary is referred to as {\bf Neumann} boundary conditions.

Finally,  we also need to specify the initial temperature distribution,
\begin{equation}
U(x,0) = f(x) \qquad 0 < x <L .
\end{equation}
We will demonstrate below that the solution to this problem (if it
exists) is unique; later in this course we will solve this problem
using the method of separation variables.

For completeness, we also comment here that the problem can be posed on the infinite line,
$-\infty < x < \infty$
sometime called the {\bf Cauchy} problem -- in this case one usually replaces the boundary condition with the specification that the temperature remains bounded as we approach infinity, 
\begin{equation}
\lim_{x \to \pm \infty}|U(x,t)| < C,
\end{equation}
for some constant $C$.  This condition may seem superfluous at first
glance, but actually is necessary to stop heat from leaking in from
infinity (speaking very, very informally an infinite source of heat
infinitely far away can have a finite effect in  a short amount of
time). If you are interested in details, look for the examples of
Tychonov in a PDE's text\footnote{See, for example, T. W. K\"orner, ``Fourier Analysis," Cambridge University Press, p. 338.}.


\section{Examples of Solution to the Diffusion Equation}

We can summarize the last section by restating a well-posed 
problem for the diffusion equation on the interval $0 < x <L$ with Dirichlet boundary conditions,
\begin{center}
{\sc The Dirichlet Problem for the Diffusion Equation\\
(Non-homogeneous Boundary Conditions)}
\end{center}
\begin{eqnarray*}
U_t=DU_{xx} \quad & 0<x < L, t>0 \quad  & {\bf \rm DE} \\
 U(0,t) =U_0 \quad U(L,t) =U_1 \quad & t>0 & {\bf \rm BC} \\
U(x,0) = f(x) & 0 < x < L . &  {\bf \rm IC} \\
\end{eqnarray*}
Solving the general problem will have to wait, but we can find some
specific solutions to the problem using the ideas of {\bf Separation
 of Variables}.  For the moment, we will restrict ourselves to 
 homogeneous boundary conditions,
 \begin{center}
{\sc The Dirichlet Problem for the Diffusion Equation\\
(Homogeneous Boundary Conditions)}
\end{center}
\begin{eqnarray*}
U_t=DU_{xx} \quad & 0<x < L, t>0 \quad  & {\bf \rm DE} \\
 U(0,t) =0 \quad U(L,t) =0 \quad & t>0 & {\bf \rm BC} \\
U(x,0) = f(x) & 0 < x < L , &  {\bf \rm IC} \\
\end{eqnarray*}
If you want, you can skip the derivation for the moment and
jump ahead to Exercise 1, if you don't mind the solution appearing {\it deus ex machina}
( a fancy term for ``out of thin air").

\subsection{A Solution to the Homogeneous Dirichlet Problem}  
Let us look for solutions to the homogeneous Dirichlet problem of the form
\begin{equation}
U(x,t)=X(x)T(t)
\end{equation}
we find from the differential equation (DE) that
\begin{equation}
XT_t = DX_{xx}T
\end{equation}
and dividing by $XT$ we find
\begin{equation}
\frac {T_t}{DT} = \frac {X_{xx}}{X} = - \lambda .
\end{equation}
where $\lambda$ is to be determined.  Now because $T_t/DT$ is {\bf
only} a function of $t$ and $X_{xx}/X$ is {\bf only} a function of
$x$ we know that $\lambda$ must be independent of $x$ and $t$
respectively, and therefore must be a constant -- consequently it is
known as the {\bf separation constant}. 
We can now solve the resulting ODE for $T(t)$
\begin{equation}
 T_t = -\lambda DT  \quad \Rightarrow \quad T(t) = e^{-\lambda t},
\end{equation}
or some constant multiple of it.

We now look for a solution for the $X(x)$ equation that also satisfies the homogeneous boundary conditions.
From the boundary conditions (BC), we know that 
\begin{eqnarray}
U(0,t) &=& X(0)T(t)=0 \quad \Rightarrow \quad  X(0)=0 \\
U(0,t) &=& X(L)T(t)=0 \quad \Rightarrow \quad X(L)=0
\end{eqnarray}
So finally we conclude that we are looking for solutions to the {\bf Boundary Value Problem} for 
$X(x)$,
\begin{equation}
\boxed{
X_{xx} + \lambda X =0, \qquad   
X(0) =0 \quad X(L) =0.} 
\end{equation}
Solving the DE, we find that 
\begin{equation}
X(x) = B\cos(\sqrt{\lambda}x) + C \sin(\sqrt{\lambda} x)
\end{equation}
and applying the boundary conditions we see that $X(0)=0$ implies that $B=0$, and that 
\begin{equation}
C\sin(\sqrt{\lambda}L)=0.
\end{equation}
Consequently, a non-trivial solution (that is a solution for which $X(x) \ne 0$) for $X(x)$
can be found if and only if
\begin{equation}
\boxed {\lambda = \lambda_n \equiv  \left(\frac {n \pi }{L} \right )^2  \quad {\rm for} \quad n =1,2,3 \ldots} 
\end{equation}
for which we find
\begin{equation}
\boxed{
X(x)=X_n(x) \equiv  \sin\left ( \frac{ n \pi x}{L} \right )     \quad {\rm for} \quad n =1,2,3 \ldots ,}
\end{equation}
or some constant multiple of it. These special values of $\lambda$ are called {\bf eigenvalues} 
and the associated functions, $X_n(x)$, are known as {\bf eigenfunctions}.

Multiplying the solution for $X_n(x)$ and $T(t)$ together finally yields a solution for $U_n(x,t)$,
\begin{equation}
\boxed{U(x,t)=U_n(x,t) \equiv  \sin\left ( \frac{ n \pi x}{L} \right ) 
e^{- \left ( \frac{n \pi }{L} \right )^2 Dt } ~ {\rm for} ~ n =1,2,3 \ldots .}
\end{equation}
The method of separation of variables is very powerful -- it will be one of our primary tools for finding solutions to PDE's in the coming lectures.
\vfill

\begin{exercise}
Verify that 
\begin{equation*}
U_n(x,t) \equiv  \sin\left ( \frac{ n \pi x}{L} \right ) 
e^{- \left ( \frac{n \pi }{L} \right )^2 t } \quad {\rm for} \quad n =1,2,3 \ldots ,
\end{equation*}
satisfies the diffusion equation $U_t = DU_{xx}$ and the homogeneous boundary conditions $U(0,t)=U(L,t)=0$. Explain why any linear combination of $U_n$,
\begin{equation*}
U(x,t) = \sum_{n=1}^{\infty}
a_n U_n(x,t)
=\sum_{n=1}^{\infty} a_n \sin\left ( \frac{ n \pi x}{L} \right ) 
e^{- \left ( \frac{n \pi }{L} \right )^2D t }
\end{equation*}
also satisfies the diffusion equation and the homogeneous boundary condition. Does it worry you that this is an infinite sum? What initial condition, $U(x,0)$, does this correspond to? 
\end{exercise}
\vfill
\eject

\subsection{A Solution to the Cauchy Problem}
We can also consider a solution to the Cauchy problem for the diffusion equation, which you hopefully remember is the problem posed on the entire real line,
 \begin{center}
{\sc The Cauchy Problem for the Diffusion Equation}
\end{center}
\begin{eqnarray*}
U_t=DU_{xx} \quad & - \infty <x < \infty , t>0 \quad  & {\bf \rm DE} \\
 \lim_{x \to \pm \infty}|U(x,t)| < C  & t>0, & {\bf \rm BC} \\
U(x,0) = f(x) & -\infty < x < \infty . &  {\bf \rm IC} \\
\end{eqnarray*}
While there are many clever derivation for the solution to this problem, for the moment I will simply give you the most important solution, usually called the 
{\bf fundamental solution} or the {\bf diffusion kernel},
\begin{equation}
U(x,t) =G(x,t+\tau) \equiv \frac{1}{\sqrt{4 \pi D (t+\tau)}} e^{-\frac{x^2}{4D(t+\tau)}}.
\end{equation}
where $\tau$ is a constant (which we will assume is positive).  This solution can be used to construct a general solution of the diffusion equation for an arbitrary initial condition, $f(x)$.
\vfill
\begin{exercise}
Verify that 
\begin{equation*}
G(x,t+\tau) \equiv \frac{1}{\sqrt{4 \pi D (t+\tau)}} e^{-\frac{x^2}{4D(t+\tau)}}.
\end{equation*}
satisfies the diffusion equation and the boundary conditions for the Cauchy problem
when $\tau >0$ . Show that this solution corresponds to a Gaussian with time varying width and height. How does the Gaussian's  width, height and area vary in time?
\end{exercise}
\vfill
\eject





\section{The Maximum Principle}

Looking at solutions to the heat equation, we note that they tend to average out maximums and minimums. We can develop some intuition for this by considering 
what the equation says. Basically, $U_t = DU_{xx}$ means: {\it The temperature is decreasing when the profile is convex down and the temperature is increasing when the profile is convex up}.
 
 

\vskip 3 in

\begin{center}
Figure 2.4: The heat equation interperted graphically\\
 (draw your own figure).
\end{center}
From which we conclude that interior maximums in temperature are decreasing and interior minimums of temperature are increasing. This reasoning is not quite airtight
(how to make it tighter is a good question to ponder). We can give a rigorous statement (without proof) of the maximum principle:
\begin{theorem}[Maximum Principle for the Diffusion Equation]
If \\ $u(x,t)$ satisfies the Dirichlet problem for the diffusion equation in the 
semi-infinite strip $0<x<L$, $0<t$, then it assumes its maximum value (as a function of $x$ and $t$) either initially (when $t=0$) or on the lateral boundaries (where
$x=0$ or $x=l$).
\end{theorem} 

The same  is also true of the minimum of $u(x,t)$. A proof can be found in most advanced PDE texts.

\begin{exercise}
Interpret the solutions we have found for the diffusion equation in terms of the maximum principle. Show examples where the maximum value of $u(x,t)$ occur 
in the initial condition and on the lateral boundaries.
\end{exercise}

\section{Energy Dissipation and Uniqueness}

By looking at what is normally known as energy for the diffusion equation, we can show that the solution for the Dirichlet problem is unique. Note this energy is a mathematical construct, not to be confused with the thermal energy discussed in the derivation of the diffusion equation.

First, suppose that $U(x,t)$ is a solution to the homogeneous Dirichlet problem,
\begin{eqnarray*}
U_t=DU_{xx} \quad & 0<x < L, t>0 \quad  & {\bf \rm DE} \\
 U(0,t) =0 \quad U(L,t) =0 \quad & t>0 & {\bf \rm BC} \\
U(x,0) = f(x) & 0 < x < L , &  {\bf \rm IC} \\
\end{eqnarray*}

Let's define, the energy,
\begin{equation}
W = \frac{1}{2} \int_0^L ~ U^2 ~ dx ,
\end{equation}
which is a function of $t$ dependent on the particular solution $U(x,t)$ (technically it is a function of $t$ and a {\bf functional} of $U(x,t)$). Note that $W \ge 0$ with 
$W=0$ only for the trivial solution $U(x,t)=0$.

If we differentiate the energy with respect to time, we find
\begin{eqnarray*}
\frac{dW}{dt}  &=& \int_0^L ~ UU_t ~ dx ,\\
                      &=& D \int_0^L ~ UU_{xx} ~ dx ,\\
                      &=&  - \int_0^L ~ (U_x)^2 ~ dx  + UU_{xx} \left. \right |_{x=0}^{x=L}, 
\end{eqnarray*}
where we have substituted the DE and used integration by parts. Now, applying the BC's, we fine that the boundary terms from the integration by parts vanish, so,
\begin{eqnarray*}
\frac{dW}{dt}  =  - \int_0^L ~ (U_x)^2 ~ dx  \le 0
\end{eqnarray*}
Now, we can conclude that $W$ is decreasing (that is energy is dissipated!!){\bf unless} $U_x=0$, that is to say 
that $U$ is constant. As the only constant solution satisfying the boundary conditions
is $U=0$, we might be tempted to conclude that the solution always decays to this trivial state. This turns out to be true, although one must invest some analysis to show it rigorously.
\vfill
\eject
\phantom{shrdlu}
\vskip 3 in

\begin{center}
Figure 2.4: (a) A solution to the homogeneous Dirichlet problem for the heat  equation.
(b) The corresponding energy, W, which is decreasing to zero.\\
 (draw your own figure).
\end{center}


A second conclusion one can reach is that if $f(x)=0$, that $U(x,t)=0$
for all $t>0$.  This follows quickly because $W=0$ at $t=0$, it is
non-increasing and non-negative.  While this seems like a trivial
result, it has a very powerful consequence. Suppose we had two
solutions to the non-homogeneous Dirichlet problem, call them $V_1$
and $V_2$. You should be able to convince yourself that there
difference $U=V_1-V_2$ satisfies the homogeneous Dirichlet problem
with $f(x)=0$. Consequently, we know that $U(x,t)=0$ for all $t>0$,
which implies $V_1=V_2$. From this we conclude that {\it The solution
  to the non-homogeneous Dirichlet problem is unique}, a powerful
result indeed.

\vfill

\begin{exercise}
Convince yourself the energy argument for uniqueness of solutions in
the previous paragraph is correct. Show that a similar argument can
be made for the Neumann problem.
\end{exercise}



\vfill
\eject

\section{Challenge Problems for Lecture 2}



\begin{problem}
Consider the diffusion equation with homogeneous Neumann boundary conditions.
\begin{eqnarray*}
U_t=DU_{xx} \quad & 0<x < L, t>0, \quad  & {\bf \rm DE} \\
 U_x(0,t) =0 \quad U_x(L,t) =0 \quad & t>0 ,& {\bf \rm BC} \\
U(x,0) = f(x) & 0 < x < L .&  {\bf \rm IC} \\
\end{eqnarray*}
\end{problem}
\begin{enumerate}
\item Explain physically why this corresponds to the diffusion of heat in a metal bar with insulated ends.
Make sure you understand what each of the equations corresponds to.

\item Show that  
\begin{eqnarray*}
&{\rm  (i)}& ~ U_0(x,t)=1 \\
&{\rm  (ii)}& ~ U_n(x,t)=\cos  \left ( \frac{n \pi x}{L} \right )
e^{-\left ( \frac{n \pi}{L} \right )^2 Dt} \qquad n = 1,2,3, \ldots
\end{eqnarray*}
satisfy both the diffusion equation (DE) and the homogeneous Neumann boundary conditions (BC).

\item Write down a general solution as a linear combination of the solutions you found in part (b). 
What does this say about $f(x)$ if we assume that this solution also satisfies the initial condition (IC)?
\end{enumerate}


\vfill

\begin{problem}
In this problem, we will argue that for the homogeneous Neumann problem discussed in Problem 1,
that the solution approaches a constant temperature, given by the average of the initial temperature.

\begin{enumerate}
\item
Suppose we define the total heat energy in the bar as
\begin{equation*}
Q(t) = \int_0^L ~U(x,t)~ dx .
\end{equation*}
Show that $Q$ is {\bf conserved}, that is that it is independent of time (Hint: compute $\frac{dQ}{dt}$).

\item
Use the initial condition to compute $Q$ in terms of $f(x)$.


\item 
Modify the energy argument in the previous section show that the energy is decreasing unless $U(x,t)$ is constant. Use this to argue that $U(x,t)$ approaches a constant solution as $t \to \infty$.

\item 
Finally, use parts (a) and (b) of the problem to show that there is only one possible
constant solution for $U$ that is consistent with the conservation of $Q$.  Show that  solution corresponds to the bar approaching the average temperature of the initial condition.

\end{enumerate}

\end{problem}
\vfill
\eject

\begin{problem}
Previously we showed that the diffusion kernel,
\begin{equation*}
U(x,t)=G(x,t+\tau) \equiv \frac{1}{\sqrt{4 \pi D (t+\tau)}} e^{-\frac{x^2}{4D(t+\tau)}},
\end{equation*}
satisfies the diffusion equation with an initial condition
\begin{equation*}
U(x,0)= G(x,\tau) .
\end{equation*}

\begin{enumerate}
\item Show that the total heat ,
\begin{equation*}
Q(t)=\int_{-\infty}^{\infty} ~ U(x,t) ~ dx ,
\end{equation*}
is conserved, and in fact $Q(t)=1$, for the heat kernel.


\item
Show that as $\tau \to 0$, that $G(x,\tau) \to 0 $ for $x \ne 0$ and that 
$G(0,\tau) \to \infty$.

\item Explain why the solution $G(x,t)$ (i.e. with $\tau=0$) corresponds to introducing a unit amount of heat concentrated at the origin when $t=0$. This is 
called a {\it $\delta$-function} initial condition.

\end{enumerate}
\end{problem}

\vfill

\begin{problem}
A curious property of the diffusion equation is that both derivatives and integrals of 
solution satisfy the diffusion equation also. This can be used to generate new 
solutions from existing ones.

\begin{enumerate}

\item
Show that if $U(x,t)$ satisfies the diffusion equation then $\psi(x,t)=U_x$ satisfies the diffusion equation also (Hint: differentiate both sides of the diffusion equation with respect to $x$).

\item 
Generate a new solution for the heat equation by differentiating the {heat kernel}
with respect to $x$. Graph this solution (MAPLE may be useful here) -- what does it look like?

\item
Show that if $\psi(x,t)$ satisfies the heat equation, then so does
$U(x,t) = \int_{x_0}^x \psi(\xi,t) ~ d \xi$.

\item
Show that 
\begin{equation*}
U(x,t) = \frac{1}{2} + \frac{1}{\sqrt{\pi}} \int_0^{x/\sqrt{4Dt}} e^{-z^2} ~ dz
\end{equation*}
satisfies the heat equation, by showing its derivative is the heat kernel.
Graph the solution at various times -- what physical problem does this solution 
correspond to ?
\end{enumerate}

\end{problem}

\vfill
\eject





\end{document}