Transient Two-Dimensional Heat Conduction Using Chebyshev Collocation

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.


  • [Snapshot]
  • [Snapshot]
  • [Snapshot]
  • [Snapshot]
  • [Snapshot]


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.
    • Share:

Embed Interactive Demonstration New!

Just copy and paste this snippet of JavaScript code into your website or blog to put the live Demonstration on your site. More details »

Files require Wolfram CDF Player or Mathematica.

Mathematica »
The #1 tool for creating Demonstrations
and anything technical.
Wolfram|Alpha »
Explore anything with the first
computational knowledge engine.
MathWorld »
The web's most extensive
mathematics resource.
Course Assistant Apps »
An app for every course—
right in the palm of your hand.
Wolfram Blog »
Read our views on math,
science, and technology.
Computable Document Format »
The format that makes Demonstrations
(and any information) easy to share and
interact with.
STEM Initiative »
Programs & resources for
educators, schools & students.
Computerbasedmath.org »
Join the initiative for modernizing
math education.
Step-by-Step Solutions »
Walk through homework problems one step at a time, with hints to help along the way.
Wolfram Problem Generator »
Unlimited random practice problems and answers with built-in step-by-step solutions. Practice online or make a printable study sheet.
Wolfram Language »
Knowledge-based programming for everyone.
Powered by Wolfram Mathematica © 2018 Wolfram Demonstrations Project & Contributors  |  Terms of Use  |  Privacy Policy  |  RSS Give us your feedback
Note: To run this Demonstration you need Mathematica 7+ or the free Mathematica Player 7EX
Download or upgrade to Mathematica Player 7EX
I already have Mathematica Player or Mathematica 7+