Научная статья на тему 'Математическое моделирование процессов теплопереноса при охлаждении экструдированной кабельной жилы с учетом фазовых превращений полимерной изоляции'

Математическое моделирование процессов теплопереноса при охлаждении экструдированной кабельной жилы с учетом фазовых превращений полимерной изоляции Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

CC BY
28
21
i Надоели баннеры? Вы всегда можете отключить рекламу.

Аннотация научной статьи по электротехнике, электронной технике, информационным технологиям, автор научной работы — Митрошин В. Н.

Предлагается уточненная модель для расчета температурных полей экструдированных кабельных жил с изоляцией из полиэтилена высокого давления низкой плотности в процессе охлаждения проводов при их изготовлении.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «Математическое моделирование процессов теплопереноса при охлаждении экструдированной кабельной жилы с учетом фазовых превращений полимерной изоляции»

По данным «Регион-газа» цена природного газа среднего давления (без НДС) с учётом роста тарифа в настоящее время составляет 726 руб. за 1000 м^. С учётом НДС - 871.2 руб. за 1000 м3. Удельная теплота сгорания газа - 8500 ккал/м^. Следовательно, энергетический потенциал, выделяющийся при сжигании 1000 м^ природного газа, равен 8500000 ккал или 8.5 Гкал. Это количество энергии используется на нагрев воды с кпд 0.9, т.е. 8,5 Гкал х 0.9 = 7.65 Гкал. Значит, на получение 7.65 Гкал израсходовано 871,2 рубля. Отсюда стоимость газа в себестоимости выработки тепла составит 871,2 : 7.65 = 113.9 руб ./Г кал. Результирующая стоимость тепла будет определяться заработной платой персонала (она минимальна ввиду работы модуля в автоматическом режиме), накладными расходами, величиной заложенной прибыли, налогами. Можно ожидать её на уровне 140-150 руб./Гкал,

Микропроцессорная система управления модулем обеспечивает качественно* количественное регулирование процесса теплоснабжения, т.е. в зависимости от температуры наружного воздуха меняется температура горячей воды и её расход. Этим достигается большая комфортность температурного режима помещений на всех этажах отапливаемого здания. Температура в служебных помещениях в ночные часы, в выходные дни поддерживается на существенно более низком уровне, обеспечивая снижение эксплуатационных расходов. Нормальный температурный режим восстанавливается в течение часа перед началом рабочего дня. Локальное построение системы позволяет снять привычное ограничение на температуру обратной воды в 70°С Используя соответствующий алгоритм управления, можно стабилизировать ее на уровне 55-60°С, допуская снижение необходимого расхода горячей воды. В итоге возникает возможность при проектировании уменьшать диаметры труб в теплотрассах, экономя металл.

Что касается работы тепловых сетей, то ГМ НПЦ ИНТЕХ разработал технологию внешнего и внутреннего эмалирования труб с защитой сварных швов теплотрасс от коррозии в полевых условиях, что увеличивает стоимость труб вдвое. Этим гарантируется нормальная их работа в течение 40-50 лет вместо 4-6-ти. С целью обеспечения гидро- и теплоустойчивости тепловых сетей разработана методика их наладки путём моделирования процессов теплоснабжения на ЭВМ. определения необходимых величин гидравлических сопротивлений для всех узловых точек сети с одновременной их установкой на теплотрассах.

Дополнительную информацию можно получить: тел.: (8-8482) 35-01-09 (E-mail: pliss(«'iavtograd,rii) или тел, (8-8462) 37-09-87. E-mail: zolotovvpiVimail.ru.

ЬИЬЛИОГРАФИЧКС'КИЙ список

I 11Л1.:юфикании 11 тепловые сети: Учеб, для вузов / Е.Я. Соколов: М.: Энергоатомизлат. 1982. 360 с.

2. Соколов Е.Я.. Зингер II. М. Струйные аппараты. М.: Эн ер го атом из дат, 1989, 352 с.

.V Гообатис К Ф. Нол та ре и кий А.Н. Справочник по котельным установкам малой производительности - Пол рсл К Ф Ролдатиса. М.: ')пер[ оатомиздат. 1989. 488 с.

Статья поступила « редакцию 2 ’ августа 2004 г.

УДК 621.315 В.Н. Митрошип

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ПРОЦЕССОВ ТЕПЛОПЕРЕНОСА ПРИ ОХЛАЖДЕНИИ ЭКСТРУДИРОВАННОЙ КАБЕЛЬНОЙ ЖИЛЫ С УЧЕТОМ ФАЗОВЫХ ПРЕВРАЩЕНИЙ ПОЛИМЕРНОЙ ИЗОЛЯЦИИ

Предлагаете» уточненная модель для расчета температурных тлей экструдированных кабельных жги с изоляцией из полиэтилена высокого давления низкой плотности (ПЭНП) в процессе охлаждения проводов при их изготовлении. Разработанная модель учитывает преобразование фаз. имеющее место при охлаждении пластмассы, и связанные с uu.\t значительные тепловыделения при фазовом переходе, обусловленные кристаллизацией полимера.

Одной из важнейших технологических операций производства кабелей с пластмассовой изоляцией является операция изолирования. При изготовлении кабелей с полиэтиленовой изоляцией ее наложение осуществляется на экструзионных линиях, содержащих помимо червяч-

ных экструдеров охлаждающие ванны, обеспечивающие постепенное охлаждение полиэтиленовой изоляции. По существующей в настоящее время технологии изолирования кабельных жил с пластмассовой изоляцией применяют так называемый режим воздушно-водяного охлаждения.

При пропускании кабельной жилы через ванну охлаждения в изоляции возникает значительный градиент температурного поля по радиусу изоляции, обусловленный, в частности, тонкостенностью изоляции (составляющей единицы миллиметров), а также существенной разницей температур расплавленной изоляции и охлаждающей среды - воды в ванне.

Наличие значительного градиента температуры может привести к неравномерной деформации (усадке) полимера в различных ее слоях по радиусу, что вызовет появление механических напряжений в изоляции; может явиться причиной образования внутренних пустот в ней и даже вызвать эффект так называемой “закалки” изоляции - быстрого охлаждения и кристаллизации ее.

Появление вы'неперечисленных эффектов неизбежно приведет к ухудшению электрических характеристик изоляции, а также к ее преждевременному старению и даже механическому разрушению. Для уменьшения градиента температуры используют многосекционные ванны охлаждения с различной температурой охлаждающей воды в них.

Регулирование процесса охлаждения изоляции можно осуществлять как изменением длины участков воздушного и водяного охлаждения, так и температурой охлаждающей среды (например, воды в ваннах охлаждения).

Для математического описания и температурного поля в изолированной кабельной жиле принимаем следующие допущения.

1. Изолированная жила считается симметричной относительно своей оси,

2. Жила движется с постоянной скоростью.

3. Материал жилы и изоляции является изотропным.

4. Не учитываются изменения размеров провода, вызванные усадкой изоляции,

5. Процесс считаем стационарным.

Зависимость температуры любой точки изоляции кабеля от времени его охлаждения, радиальной и осевой координат при отсутствии внутренних источников тепла есть решение уравнения нестационарной теплопроводности [1]:

где Т2 ~ температура изоляции; г - текущий радиус; г - продольная координата,

/ - текущее время;

а2 - коэффициент температуропроводности материала изоляции;

^2 - коэффициент теплопроводности материала изоляции;

V - скорость изолирования (вытяжки); г0 - радиус внутреннего проводника;

Яиз -■ радиус жилы по изоляции;

£ - общая длина участка охлаждения.

Обычно в качестве изоляции радиочастотных кабелей используется полиэтилен высокого давления низкой плотности (ПЭНП). Для ПЭНП при комнатной температуре коэффициент теплопроводности А = 0,335 Дж/м-с-К, а теплоемкость С = 2,3 кДж/(кг-К).

Для внутреннего металлического проводника зависимость температуры любой точки Г1(г.г,/)есть решение уравнения нестационарной теплопроводности [2]:

дг2 г дг А. 2 ог дг дг2

с!

(1)

^гМг.2,1)

Го* г 0<:<Ь;

дт^(л-,2г3/) 1 эт;(г,г,/),

--- -—а 1 --------------------г-н г г ^ л

---------—а] ---------г------1-------------------г----------

дг дг г дг 8г~

(2)

О < г <г„ ; 0 < ( <= ™;

и ^ у

здесь Т] - температура проводника;

а; - коэффициент температуропроводности материала проводника.

С учетом процесса выделения тепла при кристаллизации полимера при охлаждении кабельной изоляции в уравнение (1) необходимо добавить член, учитывающий наличие внутренних источников тепла Е(г.-,/) [3]:

-/ \ ¥Ят дг\ ,

а(г,г,/)=-(3)

где ()т - экспериментально определяемое максимальное тепловыделение при кристаллиза-

ции полимера;

г[ - относительная величина, принимающая значения от 0 до 1 и характеризующая степень завершенности процесса кристаллизации;

С’2 - удельная теплоемкость полиэтилена.

Р

Ч = ^: (4)

ИР

р -- степень кристалличности полимера;

[3^ - равновесная степень кристалличности полимера (при комнатной температуре для ПЭ равна 0.7 [4]).

Для полиэтилена ()щ ~ [377 кДж/кг] [3].

Математическая модель неизотермической кристаллизации, описывающая изменение степени кристаллизации полимера в процессе охлаждения изоляции, описывается кинетическим \ равнением [3]:

ар(г,г.г)

дг

\ ко ■ехР

Я-Т(г,и) т" -Т(г,2,г)

В уравнении (5) К(), и, Сд , Т - кинетические константы, полученные экспериментально:

Ко = 0.04 [с_1], и = - 1000 [Дж/моль]; С# =1,8; Т = 43 [К];

I? =8,31441 [Дж/моль-К] - универсальная газовая постоянная.

Зависимость коэффициента теплопроводности ПЭНП /(т) от температуры приведена в [5]. Данная кривая может быть аппроксимирована линейной зависимостью:

Х2 = 0.488-0,0017-7’. (6)

Так как толщина кабельной изоляции, как правило, невелика и составляет единицы миллиметров. разность температур наружного и поверхностного слоев изоляции не превышает при охлаждении 40-60°С. то в первом приближении можно считать, что в указанном интервале температур коэффициент теплопроводности ПЭНП постоянен. Тогда уравнение (1) с учетом сказанного и уравнения (3) принимает вид

с72(>\=.!)

■ = а-

с/

у.дя, сіі(г,,л).

уЗТ2(г,г,{) |

г0 < Г <Яуз ; 0 < 2 < І.

На поверхности контакта изоляции и проводника выполняется граничное условие четвертого рода, что соответствует теплообмену тел, находящихся в тепловом контакте (температуры сопрягающихся поверхностей одинаковы) [2]:

= ^Ъ^Л. 7;(го,/) = Г,(го,0- (§)

Ог дг

На внешней поверхности изоляции, охлаждаемой теплоносителем (водой - в ванне охлаждения, либо воздухом ~ на участке воздушного охлаждения), выполняется граничное условие второго рода [2];

— (9)

где ап - коэффициент теплоотдачи на поверхности изоляции, [Вт/м2 К],

Тв - температура теплоносителя (воды или воздуха - для участков водяного или воздушного охлаждения соответственно).

Для определения коэффициента теплоотдачи на поверхности изоляции ап используются аналитические зависимости для случая продольного обтекания цилиндра с учетом нелинейности свойств воды [6].

Начальные условия для металлического проводника в общем случае:

7'1(/-,0,/) = У,(г,г). (10)

Обычно для лучшей адгезии (прилипаемости) расплавленной изоляции к внутреннему проводнику его предварительно подогревают до температуры порядка 90’С. Справедливо допущение:

/1М = Т„ , (И)

где Тп - температура предварительного подогрева проводника.

Условия в центре жилы могут быть записаны следующим образом:

ЭГ1(0,г,г) = (); {12)

дг

Начальное условие для изоляции:

Г2(лО,0 = /2(г,0- 03)

При хорошей гомогенизации расплава изоляции в зоне дозирования экструдера справедливо допущение:

/2 М=7>. (14)

где Тр - температура расплава изоляции.

Переходя к безразмерной записи уравнений (2), (7) и (5), с учетом (4) получим

аЭ1Су,/,т)_ а2Э|(д:,/,т) [ 1 эе)(х,лт)| 2.5Ч0ил). (15)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

5т дх2 х Зх 1 д12 ' '

/е[0,1]; лге [О,.*:,];

ае,(*./,х) э202(х,;,х), 1 502(х:/,т) 2 а2е2(.г,/,х) ^ ае2(х,/,х)

дт = ал-2 5л- т ' 812 * ' Ы

51

/ е [од]; д:е[л1;1];

(16)

д! V

•[1 + С0-рр-л(я,/,т)];

и ¥

ЯТ\{х,1,т) Г*{1-е2(х,/, т)}

(П)

здесь

ш 2 (Я,л2 У-1

01=4; е3=Ч; / = 7; * = Т = 4т’ Т2= т-р л =

т т I йцз2 V I )

г - А_ ■ о* л1 “ Т ! У»! -

(18)

К/ г-С,

где 0[ и02- безразмерная температура проводника и изоляции соответственно; Т" - температура приведения (температура плавления полимера); I - безразмерная осевая (продольная) координата; х - безразмерная радиальная координата; х1 - граница сопряжения двух неоднородных сред (изоляции и металлического проводника); т - безразмерное время; Ре - число Пекле;

Qm - безразмерное максимальное тепловыделение при кристаллизации полимера.

Граничные условия:

дв2(хь1,х)_ ЗО^./.т).

п. ~ a j А, “ э Vi5v

ОХ ОХ А-2

e2(x1,/,x)=0i(jci)/,t); (20)

- а (',т) ■ ■ [е2(1.Лг)-е,(/,т)], (21)

где 0е (/,т)- относительная температура охлаждающей среды (воды, например).

(22)

В безразмерной форме начальные условия записываются следующим образом:

01(jc,/,O)=91o(x,/); 92(jc,/,0)= 02О(л,/);

Т Т (23)

е1(хДт)-01(х,т) = ^- = 0*о; 02(*,О,т) = 02{л,т)=^- = 021 :

Э01(О,/,т)=О;

дх

Таким образом, уравнения (15)-(17) с граничными условиями (19)-(22) и начальными условиями (23) являются математической моделью процессов теплопереноса и фазовых превращений при охлаждении кабельной жилы с экструдированной полимерной изоляцией.

БИБЛИОГРАФИЧЕСКИ СПИСОК

1. MtunpouiUH В.Н. Математическая модель для расчета температурных полей экструдированных изоляционных покрытий кабельных жил Н Вести. Самар, гос. техн, ун-та. Сер. Физико-математические науки. 2004. Вып. 27. С. 8890.

2. Карташов ЭМ. Аналитические методы в теории теплопроводности твердых тел. М.: Высш. шк., 2001. 550 с.

3. Зиннатуллин P.P., Труфанова И.М., Шилинг А,А. Исследование процессов теплопереноса и фазовых превращений при охлаждении провода с полимерной изоляцией Н V Минский между нар. форум по тепло- и массообмену; Тез. докл. и сообщ. Т. 2. 24-28 мая 2004 г, Минск, 2004. С. 130-131.

4. Полиэтилен. Справочное руководство / Под ред. М.И. Гарбара. Л.: Госхимиздат, 1955.

5. Оганесян К.Г. Расчет времени вулканизации при изолировании жил кабелей вулканизирующимся полиэтиленом // Электротехническая промышленность. Серия: Кабельная техника, 1975. № 12. С. 7-10.

6. Теплотехнический справочник. / Под ред. В.Н.Юренева Т. 2. М.: Энергия, 1976, 896 с.

Статья поступила в редакцию 18 ноября 2004 г.

УДК 621.311.22.002.56; 62! .317 В.И. Немченко, С.Ю. Абрамов

СИСТЕМНЫЙ ПОДХОД К ОГАНИЗАЦИИ УЧЕТА ТЕПЛОВОЙ ЭНЕРГИИ НА МУНИЦИПАЛЬНОМ УРОВНЕ

Представлены результаты системного анализа основных факторов, определяющих выбор приборов систем коммерческого учета тепловой энергии и теплоносителя. Изложены основные результаты внедрения приборного учета на муниципальном уровне.

Рассмотрены результаты системного подхода к решению многокритериальной технической задачи анализа работы системы централизованного теплоснабжения на примере г, Ново-куйбышевска Самарской области [1]. Организационно-технические мероприятия проводились в несколько этапов.

i Надоели баннеры? Вы всегда можете отключить рекламу.