One-Dimensional Heat Conduction with Temperature-Dependent Conductivity

Consider the following nonlinear boundary value problem , with , , and . This equation with the boundary conditions (BCs) describes the steady-state behavior of the temperature of a slab with a temperature-dependent heat conductivity given by .
Introducing nondimensional temperature and position using and , the governing equation and BCs become , with , , and . Without loss of generality, set . Then this problem admits an analytical solution given by .
We use Galerkin's method to find an approximate solution in the form . The unknown coefficients of the trial solution are determined using the residual and setting the integrals for . You can vary the order of the trial solution, . The Demonstration plots the analytical solution (in gray) as well as the approximate solution (in dashed cyan). The polynomial expression of the trial solution is displayed on the same plot. In addition, the difference between the exact and approximate solutions, , is plotted. This error is seen to decrease as the order of the trial solution, , is increased.


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


[1] W. F. Ramirez, Computational Methods for Process Simulation, 2nd ed., Oxford: Butterworth–Heinemann, 1997.
    • 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.