Рассмотрим неограниченный цилиндр радиуса R, температура поверхности которого остается неизменной на протяжении всего процесса теплообмена. Радиальное распределение температур в начальный момент задано в виде некоторой функции Т(r). Необходимо найти распределение температур. Такие задачи встречаются при расчете процессов охлаждения полимерного волокна, затвердевания литников литьевых форм и т. п.

Дифференциальное уравнение теплопроводности для цилиндра имеет вид:

(1.1)

Краевые условия: (1.2)

(1.3)

(1.4)

Решение, полученное методом разделения переменных, имеет сложный вид потому задачей данной работы является найти численное его решение.