Transient Two-Dimensional Heat Conduction Using Chebyshev Collocation

Initializing live version
Download to Desktop

Requires a Wolfram Notebook System

Interact on desktop, mobile and cloud with the free Wolfram Player or other Wolfram Language products.

Consider the two-dimensional heat equation given by


, with and ,

which represents heat conduction in a two-dimensional domain. The boundary conditions are such that the temperature, , is equal to 0 on all the edges of the domain:

and for

and for

Without loss of generality, one can take the thermal diffusivity, , equal to . The initial condition is given by

for .

The dimensionless temperature, , can be found using either NDSolve or the Chebyshev collocation technique. As shown in the table of data given at , the two methods give the same results. Also, three contours of the dimensionless temperature, , (i.e., 0.25, 0.5, and 0.75) at are shown using the red curves and dashed black curves for the solution obtained with NDSolve and Chebyshev collocation, respectively. Again, perfect agreement is observed. Finally, you can set the value of the dimensionless time, and the contour plot tab lets you display the contour plot of the solution obtained using the Chebyshev collocation method.


Contributed by: Housam Binous, Brian G. Higgins, and Ahmed Bellagi (April 2013)
Open content licensed under CC BY-NC-SA



In the discrete Chebyshev–Gauss–Lobatto case, the interior points are given by . These points are the extremums of the Chebyshev polynomial of the first kind .

The Chebyshev derivative matrix at the quadrature points, , is an matrix given by

, , , for , and for and ,

where for and .

The discrete Laplacian is given by , where is the identity matrix, is the Kronecker product operator, , and is without the first row and first column.

An affine transformation, , allows shifting from the interval to .


[1] L. N. Trefethen, Spectral Methods in MATLAB, Philadelphia: SIAM, 2000.

Feedback (field required)
Email (field required) Name
Occupation Organization
Note: Your message & contact information may be shared with the author of any specific Demonstration for which you give feedback.