ФИЗИКО–ХИМИЧЕСКИЕ ПРОЦЕССЫ ПРИ ПОЛИМЕРИЗАЦИИ

На правах рукописи
Иванова Евгения Владимировна
ФИЗИКО–ХИМИЧЕСКИЕ ПРОЦЕССЫ
ПРИ ПОЛИМЕРИЗАЦИИ ИЗОЛЯЦИОННЫХ СЛОЕВ
КАБЕЛЬНЫХ ИЗДЕЛИЙ В УСЛОВИЯХ
НЕСТАЦИОНАРНОГО СОПРЯЖЕННОГО
ТЕПЛООБМЕНА
01.04.17 – Химическая физика, горение и взрыв, физика экстремальных
состояний вещества
01.04.14 – Теплофизика и теоретическая теплотехника
АВТОРЕФЕРАТ
диссертации на соискание ученой степени кандидата
физико-математических наук
Томск – 2013
Работа выполнена на кафедре автоматизации теплоэнергетических процессов
ФГБОУ ВПО «Национальный исследовательский Томский политехнический
университет».
Научные руководители:
доктор физико-математических наук,
профессор
Кузнецов Гений Владимирович
доктор физико-математических наук
Стрижак Павел Александрович
Официальные оппоненты:
Ханефт Александр Вилливич
доктор физико-математических наук,
профессор,
ФГБОУ ВПО «Кемеровский
государственный университет»
Мамонтов Геннадий Яковлевич
доктор физико-математических наук,
профессор,
ФГБОУ ВПО «Томский государственный
архитектурно-строительный университет»
Ведущая организация:
ФГБОУ ВПО «Пермский национальный
исследовательский политехнический
университет»
Защита состоится 27 сентября 2013 г. в 10.00 часов на заседании
диссертационного совета Д 212.269.13 при Национальном исследовательском
Томском политехническом университете в аудитории 217 учебного корпуса
№ 8 по адресу: 634050, г. Томск, ул. Усова, 7.
С диссертацией можно ознакомиться в Научно-технической библиотеке
Национального исследовательского Томского политехнического университета.
Автореферат разослан: «15» августа 2013 г.
Ученый секретарь диссертационного совета
Д 212.269.13
кандидат технических наук
А.С. Матвеев
Общая характеристика работы
Актуальность темы. Процессы полимеризации изоляционных слоев типичных кабельных изделий предполагают движение последних через специализированные нагревательные камеры. Основной характеристикой процесса является полнота завершения полимеризации в изоляционном слое. Наиболее распространенная технология полимеризации изоляции предполагает воздействие
на образец беспаровой среды (например, воздуха) с температурой 180–200 0С.
Нагрев до таких достаточно высоких температур обеспечивает снижение доли
влаги и газовых включений в изоляционном покрытии (при этом реализуется
защита от коррозии, эрозии и других негативных факторов). Однако, вместе с
тем возникает опасность перегрева изделия и деформации его поверхности в
результате плавления и термического разложения. Достаточные и предельно
допустимые температуры, характерные времена прогрева, а также способы полимеризации могут существенным образом меняться в зависимости от компонентного состава изоляционного материала и его структуры.
Полимеризация может реализовываться несколькими способами, отличающимися по агрегатному состоянию полимеризуемой системы. Однако, вне зависимости от способа полимеризации и агрегатного состояния изоляционного
материала рассматриваемые физико–химические процессы являются энергозатратными. Для выбора оптимальных по масштабам привлекаемых материальных ресурсов и качеству изделий технологических параметров целесообразно
прогностическое моделирование теплового состояния изделия в процессе его
полимеризации.
Вопросам исследования механических, теплофизических, термохимических и других характеристик типичных изоляционных материалов, макроскопических закономерностей процессов изготовления и эксплуатации последних
посвящено достаточно много работ. Можно выделить труды Лукомской А.И.,
Баденкова П.Ф., Кеперша Л.М., Панфиловой О.А., Аваева А.А., Осипова Ю.Р.,
Осипова С.Ю., Короткова В.Н., Цаплина А.И., Милехина Ю.М., Янкова В.И.,
Труфановой Н.М., Ковригина Л.А., Силуянова А.Ю., Сырчикова И.Л., Татарникова А.А., Андыка В.С., Кима В.С., Леонова А.П., Бутакова А.А., Седова И.В.,
Халтуринского Н.А., Левина П.П., Худяева С.И., Берлина А.А., Штейнберга
А.С., Шкадинского К.Г., Зеленского А.Д., Холодного С.Д., Назарова Э.С., Михайлова Ю.М., Куличихина С.Г., Захватова А.С. и других исследователей.
Однако, несмотря на достаточно многочисленные, на первый взгляд, научные работы и многолетний производственный опыт, технологические процессы
изготовления изоляции не проработаны в полной мере теоретически. Во многом это можно объяснить тем, что ранее не исследовались процессы полимеризации с учетом реальной геометрии изделий, а также сопряженного конвективного, кондуктивного и радиационного теплообмена в нагревательных камерах.
Отсутствие теории и, как следствие, прогностических моделей, сказывается как
на эффективности производства (брак, неэффективное расходование материальных и временных ресурсов), так и безопасности технологических процессов
(перегрев, плавление, тление и возгорание изоляционных материалов).
3
Разработка математических моделей и методов решения соответствующих
задач полимеризации в условиях сложного теплообмена может позволить существенно повысить качество выпускаемых кабельных изделий и снизить энергоемкость технологических процессов их изготовления.
Цель диссертационной работы – теоретический анализ макроскопических закономерностей полимеризации изоляции кабельных изделий при помощи математического моделирования комплекса нестационарных взаимосвязанных физико–химических процессов, протекающих при их изготовлении.
Для достижения поставленной цели решались следующие задачи:
1.
Разработка физических и математических моделей процессов полимеризации в условиях нестационарного сопряженного конвективного, кондуктивного и радиационного теплопереноса в системе «кабельное изделие – воздух –
нагревательная камера».
2.
Выбор методов решения нестационарных дифференциальных уравнений
теплопереноса в частных производных. Разработка алгоритмов решения задач.
Тестирование полученных моделей на адекватность. Выполнение эксперимента
и сопоставление полученных теоретических и экспериментальных результатов.
3.
Численное исследование макроскопических закономерностей полимеризации изоляционных слоев кабельных изделий с учетом влияния таких факторов и процессов, как конвективный, кондуктивный и радиационный теплообмен, параметры нагревательной камеры, скорости движения и времена прогрева
изделий, характерные толщины изоляции и другие.
4.
Анализ влияния термохимических и теплофизических характеристик материала оболочки кабельного изделия на времена его полной полимеризации.
5.
Исследование влияния внутренней структуры кабельного изделия на времена полимеризации его внешней оболочки.
6.
Разработка рекомендаций по повышению качества изоляции кабельных
изделий и энергоэффективности технологических процессов их изготовления.
Научная новизна работы. Впервые поставлена и решена группа задач математического моделирования комплекса физико–химических процессов, протекающих в типичных кабельных изделиях при полимеризации в рамках моделей, учитывающих двумерный нестационарный сопряженный конвективный,
кондуктивный и радиационный теплоперенос. Задачи не имеют аналогов по постановке, алгоритму решения и полученным результатам.
Практическая значимость работы. Разработанные физические и математические модели, алгоритмы численного решения задач теплопереноса могут
быть использованы для анализа качества изоляции кабельных изделий и повышения энергоэффективности процессов полимеризации при их изготовлении.
Сформулированные в тексте рукописи рекомендации позволят существенно
повысить ресурсоэффективность (время, энергия, сырье и т.д.) рассматриваемых в диссертации процессов.
Степень достоверности результатов численных исследований. Оценка
достоверности полученных в ходе исследований результатов проводилась проверкой консервативности используемых разностных схем. Также выполнено
тестирование выбранных численных методов и разработанного алгоритма ре4
шения основной задачи теплопереноса на ряде менее сложных нестационарных
нелинейных задач теплопроводности. Проведены экспериментальные исследования, результаты которых являются основанием для выводов об удовлетворительной достоверности результатов проведенных численных исследований.
Научные положения, выносимые на защиту:
1.
Физическая и математическая модели процесса полимеризации изоляционной оболочки кабельного изделия, отличающиеся от известных учетом двумерного нестационарного сопряженного конвективного, кондуктивного и радиационного теплопереноса.
2.
Результаты численных исследований влияния основных факторов (размеры и температура нагревательной камеры, скорость движения кабеля, внутренняя структура изделия, радиационный, конвективный и кондуктивный механизмы теплопереноса, теплофизические и термохимические параметры изоляции) на интегральные характеристики физико–химических процессов, протекающих в кабельном изделии при полимеризации.
3.
Рекомендации по повышению качества изоляционных оболочек типичных кабельных изделий и ресурсоэффективности процессов полимеризации.
Личный вклад автора состоит в постановке задач полимеризации изоляционных слоев кабельных изделий, выборе методов и разработке алгоритмов
решения поставленных при выполнении работы задач, установлении основных
закономерностей протекания физико-химических процессов при полимеризации оболочек кабелей, проведении экспериментальных исследований в реальных производственных условиях, обработке и анализе полученных результатов,
разработке рекомендаций для повышения энергоэффективности производств
кабелей, формулировке основных выводов диссертационной работы.
Апробация работы. Основные положения и результаты диссертации докладывались и обсуждались на ХIV Международной научно–практической
конференции «Современная техника и технологии» (Томск, 2008), ХV Международной научно–практической конференции «Современная техника и технологии» (Томск, 2009), VII Международной конференции «Перспективы развития фундаментальных наук» (Томск, 2010), VIII Всероссийской научно–
практической конференции «Молодежь и современные информационные технологии» (Томск, 2010), Всероссийской научно–практической конференции
«Теплофизические основы энергетических технологий» (Томск, 2010), XVI
Международной научно–практической конференции «Современные техника и
технологии» (Томск, 2010), II Международной конференции «Техническая химия. От теории к практике» (Пермь, 2010), V Российской национальной конференции по теплообмену (Москва, 2010), XVII Международной научно–
практической конференции «Современные техника и технологии» (Томск,
2011), VII всероссийском семинаре вузов по теплофизике и энергетике (Кемерово, 2011), II Всероссийской научно–практической конференции «Теплофизические основы энергетических технологий» (Томск, 2011), VII Всероссийской
научно–технической конференции «Энергетика: эффективность, надежность,
безопасность» (Томск, 2011), III Всероссийской научно–практической конференции «Теплофизические основы энергетических технологий» (Томск, 2012),
5
Всероссийской конференции «Химическая физика и актуальные проблемы
энергетики» (Томск, 2012), XIX Международной научно–практической конференции «Современные техника и технологии» (Томск, 2013).
Публикации. Основные результаты диссертации представлены в 23 печатных работах (5 статей в журналах из списка, рекомендованного ВАК). Список основных публикаций приведён в конце автореферата.
Структура и объём работы. Диссертация состоит из введения, 3 глав, заключения, списка литературы, включающего 170 наименований, содержит 33
рисунка, 7 таблиц, 129 страниц.
Краткое содержание работы
Во введении обоснована актуальность темы диссертационной работы,
сформулирована цель и задачи исследований, отражена практическая значимость и научная новизна полученных результатов, представлены защищаемые
автором положения.
Первая глава отражает современное состояние теоретических и экспериментальных исследований физико–химических процессов, протекающих при
изготовлении различных изоляционных материалов. Установлено отсутствие
результатов исследований процессов полимеризации с учетом реальных производственных факторов и процессов, в частности, нестационарного конвективного, кондуктивного и радиационного теплообмена.
Во второй главе представлена постановка задачи нестационарного сопряженного конвективного, кондуктивного и радиационного теплопереноса, позволяющая выполнять численные исследования физико–химических процессов,
протекающих в процессе полимеризации изоляции кабельных изделий. Приведено описание методов решения всех уравнений, используемых при моделировании исследуемых процессов. Рассмотрен алгоритм проверки консервативности разностной схемы, примененный для оценки достоверности полученных результатов при решении системы нестационарных дифференциальных уравнений в частных производных. Описан алгоритм решения задачи. Приведены результаты численного решения тестовых задач для верификации используемого
алгоритма и выбранных численных методов решения основной задачи. Представлены результаты решения нестационарных и нелинейных задач теплопроводности с учётом зависимости теплофизических свойств от температуры и излучения на границе, а также задачи тепловой конвекции в замкнутой области.
Принята следующая наиболее общая физическая модель исследуемого
процесса (рис. 1). Кабель состоит из жилы (1) и оболочки (2). Изделие с
начальной температурой Θ0 и постоянной скоростью Wc движется через нагревательную камеру. Оболочка нагревается при высокой температуре воздуха в
камере (3). Температура свода камеры (начальная температура воздуха в камере) Θv принималась значительно больше Θ0. Учитывалась негерметичность
нагревательной камеры на входе (R2<R<R3, Z=0) и выходе (R2<R<R3, Z=Z1) –
приток окружающего воздуха с температурой Θh<<Θv и скоростью Wc. В результате прогрева изоляционный слой кабеля полимеризуется. Полнота завер6
шения процесса определяется степенью полимеризации изоляционного материала φ. Время завершения полимеризации (φ≈1 по всей толщине (R1<R<R2) оболочки) τp является основной интегральной характеристикой процесса.
При
численном
моделировании
рассматривались кондуктивный и радиационный механизмы теплопереноса, а также смешанная конвекция воздуха в камере.
Решена
осесимметричная задача теплопереноса и химического реагирования в
безразмерной
постановке в цилиндрической системе координат (рис. 1).
Численный анализ
комплекса
физикоРис. 1. Схема области решения задачи при 0≤τ≤τp:
химических процессов
1 – жила, 2 – оболочка (изоляция), 3 – воздух в камере;
выполнен при следуU, W – компоненты скорости конвективных потоков
ющих допущениях:
1. Контакт между жилой и слоем изоляционного материала идеален. На производствах этому фактору уделяется особое внимание. Выполняются специальные испытания изделий (допустимые размеры зазоров, как правило, не
превышают 0,01 мм).
2. Кабель имеет правильную цилиндрическую форму. Из всего многообразия
возможных конфигураций кабельных изделий наиболее типичным можно
считать вытянутый вдоль оси симметрии цилиндр.
3. Теплофизические и термохимические характеристики взаимодействующих
веществ не зависят от температуры. Для рассматриваемых условий теплофизические характеристики (λ, C, ρ) элементов системы (рис. 1) изменяются
в пределах 5÷7 %. При высокой инерционности полимеризации типичных
кабельных изделий такие изменения λ, C и ρ будут несущественно влиять на
интегральные характеристики основных физико–химических процессов.
Система нелинейных нестационарных дифференциальных уравнений в
частных производных, соответствующая сформулированной постановке задачи,
при 0≤τ≤τp имеет следующий вид:
2
2
1 1  1 1 1  1



;
(1)
0<R<R1, 0<Z<Z1,
Fo1 
R 2 R R Z 2
7
1 2  22 1 2  22 q2 w2 z12




,
 2 T
Fо2 
R 2 R R
Z 2
 E 
d
d
w2  2 2 ,  2  2k20 exp   a 2  ;
dt
dt
 RT2 
R1<R<R2, 0<Z<Z1,
  2
3
3
3
1
1 3  23 
3
U
W



R2<R<R4, 0<Z<Z1,

,

R
Z
Ra 3 Pr3  R 2 R R Z 2 
(2)
(3)
(4)
 2 1   2
(5)


  R ,
R 2 R R Z 2
Pr3   2 1   2   3




U
W
U 


 
.
(6)


R
Z
R
Ra 3  R2 R R Z 2 R2  R
Здесь τ – безразмерное время; τp – безразмерное время полимеризации; R, Z
– безразмерные координаты цилиндрической системы; Fo – число Фурье; Θ –
безразмерная температура; q2 – теплота химической реакции полимеризации,
Дж/кг; w2 – массовая скорость полимеризации, кг/(м3с); z1 – продольный размер
камеры, м; λ2 – теплопроводность изоляции, Вт/(мК); ΔT – разность температур
(ΔT=Tv–T0), К; Tv – температура свода камеры (начальная температура воздуха),
К; T0 – начальная температура изделия, К; ρ2 – плотность, кг/м3; φ2 – степень завершенности процесса полимеризации; t – время, с; k20 – предэкспонент химической реакции, с–1; Еa2 – энергия активации химической реакции, Дж/моль; T2
– температура в изоляционном слое, К; U и W – безразмерные составляющие
скорости конвекции; Ra – число Рэлея; Pr – число Прандтля; Ψ – безразмерный
аналог функции тока; Ω – безразмерный аналог вектора вихря.
Начальные (τ=0) условия: 0<R<R1, 0<Z<Z1, Θ=Θ0; R1<R<R2, 0<Z<Z1, Θ=Θ0,
φ=φ0; R2<R<R4, 0<Z<Z1, Θ=Θv, Ψ=0, Ω=0.
Граничные условия при 0≤τ≤τp: на оси симметрии (R=0, 0<Z<Z1) и границах Z=0, Z=Z1, 0<R<R2, R3<R<R4 для уравнений энергии выставлялись условия
равенства нулю градиента температуры; на границах Z=0, Z=Z1, R2<R<R3 принималось условие Θ=Θh; на границах R=R1, R=R2, 0<Z<Z1 выставлялись граничные условия IV рода для уравнений теплопроводности; на границе R=R4,
0<Z<Z1 принималось условие Θ=Θv; для уравнений движения и неразрывности
на стенках камеры (R=R4, 0<Z<Z1; Z=0, Z=Z1, R3<R<R4) и поверхности кабеля
(R=R2, 0<Z<Z1) выставлялись условия прилипания, на входе и выходе из камеры
(Z=0, Z=Z1, R2<R<R3) использовались краевые условия, соответствующие смешанной конвекции в нагревательной камере.
Для решения системы дифференциальных уравнений (1)–(6) с соответствующими краевыми условиями применён метод конечных разностей. Разностные аналоги дифференциальных уравнений решены локально-одномерным
методом и методом переменных направлений. Система одномерных разностных уравнений решена методом прогонки с использованием неявной четырёхточечной разностной схемы. Для решения нелинейных одномерных уравнений
применён метод итераций.
8
Оценка достоверности полученных результатов проводилась проверкой
консервативности разностной схемы. Для верификации используемых методов
и разработанного алгоритма решены тестовые задачи теплопроводности с учётом химической реакции в материале, гидродинамики и конвективного теплообмена. Также проведен эксперимент, показавший удовлетворительное соответствие полученных значений температуры в кабельном изделии результатам
численных исследований.
Третья глава содержит результаты (представлены в размерных и безразмерных переменных) численного анализа основных макроскопических закономерностей протекания физико–химических процессов при полимеризации типичных кабельных изделий.
Показано (рис. 2),
что
использование
упрощенных постановок,
предполагающих
применение декартовой
системы координат, при
численных исследованиях приводит к отклонениям температур в
изделии до 10 К относиРис. 2. Распределения температуры в момент времени тельно моделей в циt=100 с: 1, 2 – в бесконечной пластине (с учетом и без линдрической системе
излучения), 3, 4, 5 – в цилиндре (с учетом и без излуче- координат.
ния, с учетом эффекта полимеризации)
Установлено, что
радиационный теплообмен достаточно существенно (рис. 2) влияет
на распределение температуры в изделии (отклонения составили более 10 К). Также выявлено, что непосредственно тепловой эффект процесса полимеризации незначительно
влияет на температурное
Рис. 3. Изменение степени полимеризации в окрестно- поле изоляции кабеля.
сти границы «изоляция – жила» во времени при разных
Показано (рис. 3),
температурах в камере: 1 – с учетом теплообмена излу- что процесс полимеричением, 2 – без учета излучения
зации при учете теплообмена излучением заканчивается на 300 с раньше, чем в условиях только конвективного теплообмена. В результате численных исследований получены типичные распределения степени полимеризации в процессе прогрева изделия
(рис. 4).
9
Установлено,
что
при
малых
временах
нагрева зависимость φ(r)
имеет экспоненциальный
вид
(кривая
1)
с
максимумом на границе с
горячим воздухом (r=r2).
Это обусловлено тем, что
оболочка не успевает
прогреться до температур,
при
которых
полимеризация
завершается по всему
слою (r1<r<r2). С ростом
Рис. 4. Распределения степени полимеризации по времени зависимость φ(r)
толщине оболочки кабеля в разные моменты времени: становится
практически
1 – t=1 с, 2 – t=50 c, 3 – t=100 c, 4 – t=200 c, линейной (кривые 2–5).
5 – t=300 c, 6 – t=tp
Оболочка прогревается по
толщине и разница между φ(r1) и φ(r2) уменьшаются с увеличением времени.
Процесс завершается практически полной полимеризацией (φ≈0,99, r1<r<r2) при
времени около 500 с (кривая 6).
Определены зависимости времени полимеризации от температуры в
нагревательной камере Тv (рис. 5) и ее характерных размеров rv (рис. 6).
Показано (рис. 5), что с
понижением температуры
Тv время полимеризации tp
увеличивается нелинейно.
Это обусловлено экспоненциальной
зависимостью
скорости химического реагирования при полимеризации оболочки кабеля от
температуры. При понижении Тv уменьшается энергия, подводимая к поверхРис. 5. Зависимости времени завершения полимери- ности оболочки изделия.
зации tp от температуры воздуха в камере Тv: Это приводит к снижению
1 – rv=0,02 м, 2 – rv=0,01 м, 3 – rv=0,005 м
скорости прогрева приповерхностного слоя оболочки и, как следствие, уменьшению скорости химической реакции. При повышении Тv увеличивается теплоприход в изоляционном
слое кабельного изделия. Возрастает температура в этом слое кабеля. В таких
условиях процесс вулканизации проходит быстрее (рис. 5).
На рис. 5 показано, что, например, при Тv=350 К время полной
полимеризации составляет tp=911 с. При дальнейшем снижении температуры
воздуха время tp экспоненциально увеличивается и составляет более 1200
10
секунд при Тv=300 К. При повышении Тv до 500 К значения tp существенно
(более чем в три раза) уменьшаются (рис. 5), хотя на реализацию процесса
расходуется на 40 % больше энергии в единицу времени. В таких условиях
важную роль играет баланс между временем вулканизации одного кабельного
изделия, затратами на разогрев и поддержание определенной температуры
поверхности камеры, а также производительностью технологического
процесса.
При
увеличении
размеров нагревательной
камеры
(повышении
поперечного размера rv в
частности) интенсивность
исследуемых процессов
теплопереноса снижается
(рис.
6)
и
время
tp
полимеризации
возрастает. Это можно
объяснить тем, что при
увеличении rv возрастает
Рис. 6. Зависимости времени завершения полимериза- объем воздуха в камере,
ции tp от характерного размера камеры rv: на разогрев
которого
1 – Tv=450 K, 2 – Tv=465 K, 3 – Tv=480 K
расходуется энергия. При
увеличении этого объема больше энергии затрачивается на прогрев воздуха и
меньше передается непосредственно оболочке кабеля. Особенно этот фактор
играет заметную роль при увеличении rv в несколько раз (рис. 6) относительно
значения 0,01 м и умеренной продолжительности процесса.
Выполнено исследование влияния вынужденной конвекции в
нагревательной камере вследствие притока относительно холодного воздуха на
ее входе и выходе (рис. 7).
Установлено, например, что в условиях отсутствия поступления в
камеру холодного воздуха (Тh=300 K) процесс полимеризации завершается
за время tp=284,6 с. При
этом степень полимеризации по толщине оболочки изделия составляет
φ≈0,99 (рис. 7). В условиРис. 7. Распределения степени полимеризации по ях притока в камеру
толщине оболочки: 1 – при t=tp без учета вдува окружающего воздуха с
холодного воздуха, 2 – с учетом (при t=498,2 с)
температурой Th (существенно ниже температуры Тv) даже через отверстие относительно небольшого
размера процесс вулканизации существенно замедляется и завершается только
11
при времени более 500 с. Показано, что с ростом скорости движения кабеля и,
как следствие, притока воздуха, влияние этого фактора усиливается. Можно
сделать вывод о том, что энергоэффективность процесса вулканизации существенно снижается (инерционность прогрева значительно возрастает). Полученный результат обусловлен значительным падением температуры непосредственно вблизи оболочки по сравнению с температурой верхней части камеры,
в которую холодный воздух не поступает (рис. 7).
Также выполнено сравнение интегральных характеристик полимеризации одножильных и многожильных кабелей.
Установлено,
что,
например, для многожильного
кабеля время полной полимеризации его резиновой оболочки при температуре в камере
Тv=500 К составляет tp≈775 с.
Для одножильного кабеля при
а
идентичных условиях (температура в камере, конструктивные размеры проходного сечения и т.д.) и параметрах (толщина и материал внешнего
изоляционного слоя) процесс
вулканизации завершается при
tp≈497 с. При уменьшении температуры Тv до 450 К отклонения времен tp для одножильного и многожильного кабелей
б
достигают 45 %. При увеличении Тv до 550 К эти отклонения
несколько уменьшаются (до 30
%). Можно сделать вывод о
существенном влиянии внутренней структуры кабеля на
времена полимеризации его
внешней оболочки.
В тоже время поля температур, функции тока и вихрей
скорости (рис. 8) в нагревав
Рис. 8. Температурное поле (а), изолинии функ- тельных камерах с одно– и
ции тока (б) и вихря скорости (в) при полной по- многожильными кабелями отлимеризации оболочки многожильного кабель- личаются незначительно.
ного изделия (tp=774,6 с)
Проведен анализ влияния
на условия полимеризации теплофизических и термохимических характеристик
12
материала оболочки кабельного изделия. Сформулирована группа соответствующих аппроксимационных выражений для интегральной характеристики
процесса – времени полной полимеризации tp, например:
tp(ρ)=0, 071ρ+413,26 при 1000<ρ<1400 кг/м3;
tp(C)=0, 063C+411,81 при 1100<C<1600 Дж/(кгК);
tp(λ)= –116573λ3+24865λ2–3001,5λ+659,95 при 0,12<λ<0,19 Вт/(мК).
Установлено, что при выборе режимов полимеризации изоляционных материалов на производствах особое внимание целесообразно уделять точности
определения теплофизических и особенно термохимических характеристик
изоляционных материалов. Отклонения этих характеристик в пределах 15 %
приводит к изменению интегральных параметров процесса полимеризации более чем на 30 %.
Выполнено сравнение результатов экспериментов (рис. 9) и численных исследований, проведенных для гибкого трубчатого электронагревательного элемента с использованием двух моделей теплопереноса: упрощенная постановка
(№ 1) – в декартовых координатах; наиболее полная постановка (№ 2) – в цилиндрических координатах с учетом нестационарного сопряженного кондуктивного, радиационного и конвективного теплообмена.
Установлено удовлетворительное соответствие
результатов численных исследований и экспериментов для наиболее полной
Рис. 9. Условная схема экспериментальной установ- постановки задачи полимеки: 1 – экструдер, 2 – вулканизационная камера, 3 – ризации (табл. 1). Как
охлаждающая ванна, 4 – намоточное устройство, 5 – следствие, физическую постановку, математическую
термоэлектрический преобразователь
модель, методы и алгоритм
решения задач теплопереноса, рассмотренных в диссертации, можно считать
приемлемыми и обеспечивающими достоверность полученных результатов. В
тоже время установлено, что использование упрощенных моделей (в частности,
постановки № 1 в декартовых координатах без учета радиационного и конвективного теплообмена) приводит к достаточно существенным отклонениям от
эксперимента (до 18 %).
Табл. 1. Сравнение результатов численных и экспериментальных исследований температуры изоляции кабельного изделия в момент времени t=144 с
Параметр
Эксперимент Постановка № 1 Постановка № 2
Температура поверхности оболочки, К
363
306
355
Абсолютная погрешность, К
–
57
8
Относительная погрешность, %
–
18
2
По результатам выполненных исследований сформулированы следующие
рекомендации для производителей кабельных изделий:
13
При выборе технологических параметров процесса полимеризации необ1.
ходимо учитывать реальную форму изделия, конвективный и радиационный
механизмы теплообмена в камере, а также смешанную конвекцию при возможных притоках окружающего холодного воздуха в окрестности участков входа и
выхода изделий из камеры.
Для повышения энергоэффективности процессов полимеризации целесо2.
образно минимизировать размеры входных и выходных отверстий в вулканизационных камерах с целью снижения влияния относительно холодного окружающего воздуха на условия прогрева изоляционных материалов.
Размеры нагревательных камер вулканизационных печей следует мини3.
мизировать и выполнять сопоставимыми с характерными поперечными размерами протягиваемых кабелей.
Режимы полимеризации изоляционных слоев кабельных изделий следует
4.
выбирать с учетом внутренней структуры кабелей. В частности, при наличии
несколько жил изделия целесообразно пропускать через камеры с меньшей
скоростью относительно одножильных.
При выборе режимов полимеризации изоляционных слоев кабельных из5.
делий целесообразно учитывать возможные предельные отклонения теплофизических и химических характеристик изоляционных материалов от номинальных.
Также в тексте рукописи сформулированы рекомендации по дальнейшему
развитию разработанного в диссертации подхода.
В заключении подведены основные итоги выполненных исследований и
обобщены теоретические и экспериментальные следствия.
Основные результаты и выводы
1.
Впервые разработана группа физических и математических моделей для
анализа макроскопических закономерностей полимеризации изоляционных
слоев типичных кабельных изделий, отличающаяся от известных учетом комплекса основных физико–химических процессов и факторов.
Выполнена верификация моделей полимеризации изоляции кабельных
2.
изделий на базе разработанного алгоритма оценки консервативности используемых разностных схем, решения группы тестовых задач и проведения экспериментальных исследований в производственных условиях. Установлено, что отклонения результатов теоретических и экспериментальных исследований не
превышают 2 %.
Определены максимальные отклонения интегральных характеристик
3.
процесса, полученные при использовании моделей одномерного и двумерного
приближения (до 30 %).
Численно исследованы масштабы влияния на интегральные характери4.
стики полимеризации (в частности, времена полной полимеризации) конвективного и радиационного теплопереноса, скорости движения изделия через
нагревательную камеру, размеров камеры, теплофизических и термохимиче-
14
ских характеристик изоляционных материалов, наличия нескольких жил в кабеле и других факторов.
Проведен анализ предельных температур нагревательной поверхности,
5.
при которых полимеризация оболочек кабелей завершается за минимальный
интервал времени с обеспечением требуемого качества изделий и высокой
энергоэффективностью.
Установлено, что наличие притока относительно холодного воздуха на
6.
входе и выходе из камеры при протяжке кабеля существенно снижает температуру его изоляционной оболочки. Это приводит к значительному (более 40 %)
увеличению времени полной полимеризации.
Показано, что для многожильных кабельных изделий целесообразно уве7.
личивать характерные времена нагрева в специализированных камерах относительно одножильных при прочих равных условиях.
Сформулированы апроксимационные выражения tp(ρ), tp(C), tp(λ), tp(Тv),
8.
tp(rv) для прогностической оценки влияния большой группы факторов на условия полимеризации.
Представленные в работе математические модели, а также сформулиро9.
ванные рекомендации можно использовать при выборе температурных режимов полимеризации изоляционных слоев типичных кабельных изделий, а также
прогностической оценке интегральных характеристик процессов полимеризации при известных параметрах работы нагревательных камер.
Основные публикации по теме диссертации
1.
Иванова, Е.В. Об одном подходе к выбору технологических параметров
процесса вулканизации при изготовлении кабелей / Иванова Е.В., Кузнецов
Г.В. // В кн.: Труды ХV Международной научно–практической конференции
«Современная техника и технологии». – Томск, 2009.– С. 274–275.
Иванова, Е.В. Математическое моделирование температурных полей в
2.
процессе вулканизации типичных кабельных изделий / Иванова Е.В., Кузнецов
Г.В. // «Известия Томского политехнического университета», 2010. – Т. 316, №
4. – С. 38–41.
Иванова, Е.В. Математическое моделирование температурного поля ци3.
линдрического изделия в процессе изготовления / Иванова Е.В., Кузнецов Г.В.
// В кн.: Труды VII Международной конференции «Перспективы развития фундаментальных наук». – Томск, 2010. – С. 472–474.
Иванова, Е.В. Математическое моделирование процесса вулканизации
4.
кабельных изделий с учетом влияния условия естественной конвекции / Иванова Е.В., Стрижак П.А. // В кн.: Труды Всероссийской научно-практической
конференции «Теплофизические основы энергетических технологий». – Томск,
2010. – С. 105–108.
Иванова, Е.В. Влияние теплообмена излучением на процесс вулканиза5.
ции эластомерных материалов / Иванова Е.В., Кузнецов Г.В. // В кн.: Труды II
Международной конференции «Техническая химия. От теории к практике». –
Пермь, 2010. – С. 259–263.
15
Иванова, Е.В. Численное моделирование теплопереноса при вулканиза6.
ции типичных кабельных изделий / Иванова Е.В., Кузнецов Г.В. // В кн.: Труды
пятой российской национальной конференции по теплообмену. – Москва, 2010.
– С. 89–91.
Иванова, Е.В. Математическое моделирование физико-химических про7.
цессов при вулканизации кабельных изделий / Иванова Е.В., Красильникова
В.О. // В кн.: Труды Всероссийской научно-практической конференции «Теплофизические основы энергетических технологий». – Томск, 2011. – с. 22-26.
Иванова, Е.В. Математическое моделирование теплофизических и тер8.
мохимических процессов при вулканизации кабельных изделий / Иванова Е.В.
// В кн.: Труды VII всероссийского семинара вузов по теплофизике и энергетике. – Кемерово, 2011. – С. 31.
Иванова, Е.В. Численное моделирование процессов тепломассопереноса
9.
при вулканизации кабельных изделий / Иванова Е.В., Аманжолова Н.А. // В
кн.: Труды XVII Всероссийской научно-технической конференции «Энергетика: эффективность, надежность, безопасность». – Томск, 2011. – С. 85–89.
10. Иванова, Е.В. Численное моделирование температурных полей цилиндрических изделий при вулканизации / Иванова Е.В., Кузнецов Г.В. // «Химическое и нефтегазовое машиностроение», 2011. – № 7. – С. 10–11.
11. Иванова, Е.В. Численный анализ потенциала энергосбережения технологий производства кабельных изделий / Иванова Е.В., Стрижак П.А. // «Известия
вузов. Проблемы энергетики», 2012. – № 7–8. – С. 72–79.
12. Иванова, Е.В. Численное моделирование комплекса теплофизических и
термохимических процессов при вулканизации кабельных изделий / Иванова
Е.В., Стрижак П.А. // «Тепловые процессы в технике», 2012. – Т. 4, № 4. – С.
187–192.
13. Иванова, Е.В. Влияние параметров вулканизационной камеры на процесс полимеризации резиновой оболочки кабеля / Иванова Е.В. // В кн.: Труды
Всероссийской конференции «Химическая физика и актуальные проблемы
энергетики». – Томск, 2012. – С. 120–121.
14. Иванова, Е.В. Численное моделирование физико–химических процессов
при вулканизации многожильного кабельного изделия / Иванова Е.В., Стрижак
П.А. // Бутлеровские сообщения, 2013. – Т. 34, №. 6. – С. 22–27.
15. Иванова, Е.В. Программа расчета температурных полей при изготовлении кабелей с учетом влияния потоков воздуха на входе и выходе из вулканизационной печи / Иванова Е.В., Стрижак П.А. // Свидетельство о государственной регистрации программы для ЭВМ № 2012617107. – Бюлл. прогр. № 2, 2012.
16. Иванова, Е.В. Программа расчета температурных полей при изготовлении кабелей с учетом их реальной конфигурации и теплообмена излучением /
Иванова Е.В., Стрижак П.А. // Свидетельство о государственной регистрации
программы для ЭВМ № 2012617285. – Бюлл. прогр. № 2, 2012.
17. Иванова, Е.В. Программа расчета температурных полей при изготовлении кабелей с учетом реальных условий процесса вулканизации / Иванова Е.В.,
Стрижак П.А. // Свидетельство о государственной регистрации программы для
ЭВМ № 2012617239. – Бюлл. прогр. № 2, 2012.
16