Цель лабораторной работы - изучение процесса охлаждения изолированной жилы. Задачей лабораторной работы является исследование влияния параметров технологического режима охлаждения изолированной жилы на процесс охлаждения с применением метода математического моделирования и численных методов.
Цель данной лабораторной работы заключается в определении температурного поля в сечениях проводника и изоляции с целью выбора рационального режима охлаждения т. е. правильного выбора скорости изолирования V0, длины охлаждающей ванны и температуры воды в секциях охлаждающей ванны.
Для определения температурного поля в сечении проводника и изоляции существуют несколько моделей, рассмотрим некоторые из них.
Модель №1
Геометрия изолированной жилы представлена на рис.1.
Рис.1. Геометрия изолированной жилы
Допущения данной модели:
– процесс стационарный;
– теплофизические параметры постоянны;
– диффузией тепла в направлении оси z можно пренебречь;
– температура жилы не изменяется.
Уравнение энергии будет иметь вид:
(1)
Граничные условия:
где Тпр – температура жилы провода; Тср – температура среды (воды) в охлаждающей ванне.
Начальные условия:
где Твых – температура изоляции на выходе из кабельной головки.
Для решения данной задачи воспользуемся методом конечных разностей. Выберем равномерную сетку по z и r:
Аппроксимируем производные и получим:
(2)
С учётом выражений (2) выражение (1), граничные и начальные условия будут иметь вид:
(3)
Преобразуем выражения (3):
, (4)
. (5)
Алгоритм решения данной задачи:
– с учётом начальных условий рассчитываем температуру Тi,1;
– по выражению (4) рассчитываем поле температур для шага по длине j+1;
– по выражению (5) рассчитываем температуры в точка на границе изоляции;
– переходим к следующему шагу по длине z и повторяем расчёт сначала.
Модель №2
Данная модель охлаждения изолированной жилы отличается от предыдущей тем, что
температура провода изменяется по длине охлаждающей ванны, а все остальные допущения остаются в силе. Схема разбиения области конечным числом узлов представлена на рис.2.
Уравнения энергии для изоляции и для жилы соответственно будут иметь вид:
(6)
Рис.2. Схема разбиения
Граничные условия:
Начальные условия:
Также как и в предыдущем случае для решения системы дифференциальных уравнений (4) используем метод конечных разностей. Разностные уравнения (использована явная разностная схема)имеют вид:
(7)
(8)
(9)
Алгоритм решения данной задачи:
– с учётом начальных условий рассчитываем температурное поле по выражениям (7) и (8) во внутренних точках областей I и II;
– по выражениям (9), рассчитывается температура в граничных точках;
– переходим к следующему шагу по длине z и повторяем расчёт сначала.
В случае решения данной задачи по неявной разностной схеме для каждой из областей на каждом шаге по длине решается система алгебраических уравнений, пересчитываются значения температур на границах и переходят к следующему шагу по длине.
Модель №3
Данная модель охлаждения изолированной жилы отличается от предыдущей тем, что отбрасывается допущение о постоянстве теплофизических характеристик материала, а все остальные допущения остаются в силе.
Рис.3. Зависимость теплофизических характеристик от температуры
Теплофизические характеристики для текущей точки можно определить (ниже приведён пример для теплоемкости с):
При решении задачи методом конечных разностей можно воспользоваться постановкой задачи для модели 2.
Алгоритм решения данной задачи:
– с учётом начальных условий рассчитываем температуру Тi,1;
– по выражениям (5) и (6) рассчитываем температурное поле во внутренних точках областей I и II;
– по выражениям (7), (8) и (9) рассчитывается температура в точках на границе изоляции и между изоляцией и жилой а также в точке, лежащей в центре жилы;
– пересчитываем значения теплофизических характеристик для каждой точки сечения;
– пересчитываем поле температур в соответствии с выражениями (5), (6), (7), (8) и (9) для этого же сечения и пересчитываем теплофизических характеристики, повторяем расчёт до тех пор пока разность между предыдущими значениями и значениями полученными в текущем расчёте не будет меньше некоторой величины, взятой за погрешность;
– переходим к следующему шагу по длине z и повторяем расчёт сначала.
Задание
В соответствии со своим вариантом задания выполнить следующие расчеты.
1) Для заданного номинального технологического режима и заданных теплофизических характеристик материалов:
– разработать алгоритм и расчетную программу;
– рассчитать температурное поле в сечении изолированной жилы на выходе из каждой секции охлаждения;
2) Исследовать влияние на распределение температуры по сечению изолировыанной жилы следующих факторов:
– теплофизических характеристик материалов;
– линейной скорости изолирования жилы;
– температур по секциям охлаждения;
– толщины изоляции жилы;
3) Провести анализ полученных закономерностей.
4) Построить графики полученных зависимостей.
5) Оформить отчет.
Варианты заданий
Для всех вариантов заданий:
мм, мм,
Lвоз = 0.6м, м, м, м
Дж/(м2 ∙с∙ ○ С)Дж/(м2 ∙с∙ ○ С)Дж/(м2 ∙с∙ ○ С)
Дж/(м2 ∙с∙ ○ С)
Технологические и теплофизические характеристки
№
V0,
м/мин
lи/ lпДж/(м∙с∙ ○ С)
rи / rпкг/м3
си / сп
Дж/кг×°С
Tвых1°С
Tвых2°С
0.24 / 440
800 / 8840
330 / 2030
0.18 / 370
770 / 9210
305 / 1940
0.32 / 410
730 / 10190
310 / 2010
0.27 / 500
690 / 8350
320 / 2000
0.25 / 280
810 / 8070
330 / 1990
0.20 / 310
750 / 12430
315 / 2020
0.34 / 390
650 / 10830
325 / 2040
0.30 / 300
710 / 8550
335 / 2100
0.38 / 470
790 / 7330
300 / 2150
0.24 / 470
800 / 14140
340 / 1900
0.18 / 300
630 / 8840
350 / 1950
0.32 / 390
700 / 9210
360 / 1800
0.27 / 310
600 / 10190
300 / 2030
0.25 / 280
880 / 8350
350 / 1940
0.20 / 500
610 / 8070
345 / 2010
0.34 / 410
820 / 12430
340 / 2000
0.30 / 370
620 / 10830
305 / 1990
0.38 / 440
810 / 8550
360 / 2020
0.34 / 470
630 / 7330
355 / 2040
0.20 / 300
700 / 14140
315 / 2100
0.25 / 390
880 / 12350
365 / 2150
0.27 / 310
610 / 10270
350 / 1900
0.32 / 280
820 / 8430
340 / 1950
0.28 / 500
620 / 7830
295 / 1800
0.30 / 410
810 / 12550
290 / 2100
6.3 Контрольные вопросы
Какие допущения принимаются при расчёте температурного поля в сечении изолированной жилы для каждой модели?
В чем отличия моделей расчёта температурного поля?
Сформулировать основные положения метода конечных разностей.
Какие виды сеток существуют?
В чём суть аппроксимации производных?
ПРИЛОЖЕНИЕ
Теплофизические характеристики полимерных материалов приведены на рис.1.
а б
Рис. 1. Зависимость теплоёмкости и плотности полимеров от температуры и давления: а – полиэтилен высокой плотности; б – полиэтилен низкой плотности алатон ТМ 1540.
Реологические характеристики полимерных материалов приведены на рис.2.