Расчёт распределения температуры в многослойной цилиндрической стенке.

Рассчитать стационарное температурное поле осесимметричной многослойной цилиндрической стенки, в одном из слоёв которой от r1 до r2 равномерно распределены внутренние источники теплоты удельной мощностью q=10^7 Вт/м^3.
Определить тепловые потоки при r=r2 и r=r3, а также построить график изменения температуры по толщине стенки.
Геометрические размеры, мм: r1=20, r2=28, r3=38, r4=40.
Теплопроводности материалов, Вт/(м*К): k1=15, k2=3, k3=80.
Коэффициент теплоотдачи с внутренней стороны стенки h1=100 Вт/(м2*К).
Температура газа внутри трубы tг1=80 град С.
Линейный тепловой поток при r=r4 qст4=5 кВт/м.

Линейный тепловой поток через r2, r3 и r4 будет один и тот же:
Q2=5000 Вт/м. Направлен он от оси, то есть его проекция на радиус-вектор положительна.
Проекция линейного теплового потока через r1 равна:
Q1=Q2-q*Pi*(r2^2-r1^2)=-7064 Вт/м - отрицательная величина.
А температура в начале первой зоны:
t1=t0-Q1/(2*Pi*r1*h)=642.113 град С
Далее решаем для первой зоны уравнение
qист=div q, где q - плотность теплового потока.

Вычисление дивергенции в полярных координатах.



Программа вычислений на Delphi, которая строит графики