Diffusion Examples
William O. Bray
University of Maine
The purpose of this notebook is to illustrate graphically the solution to the heat equation on the real line. The intial value problem has the form
IC: u(x,0)=f(x)
The Gaussian
In[29]:=
The above defines the heat kernel. In what follows are several
examples of solutions to the Cauchy problem for the heat equation.
These illustrate ideas presented in class.
In[30]:=
Out[30]=
Example 1
In[31]:=
Out[33]=
Notice that the above formula is the same as obtained in the text obtained by hand calculation and using the definition of the error function.
In[34]:=
Out[34]=
A three dimensional rendition of the solution is given as follows.
In[35]:=
Out[35]=
Example 2
In[37]:=
Out[39]=
In[40]:=
Out[40]=
The initial data above is an odd function. The following plot illustrates that our solution also gives the solution to an IBVP on the half-line with homogeneous Dirichlet boundary condition at x=0. This foreshadows the method of reflection discussed in Chapter 4 of the text.
In[41]:=
Out[41]=
Example 3
In[42]:=
Out[44]=
In[45]:=
Out[45]=
In[46]:=
In[47]:=
Out[47]=
In[48]:=
Out[48]=
In[49]:=
Out[49]=