Енергомашиностроителен факултет Катедра “Топлинна и хладилна техника” маг. инж. Анатоли Венелинов Стойнов МОДЕЛИ НА ТОПЛООБМЕННИТЕ ПРОЦЕСИ В ТЕРМОПОМПЕНИ АГРЕГАТИ АВТОРЕФЕРАТ на дисертация за присъждане на образователна и научна степен “ДОКТОР” Научна специалност 02.06.13 “Промишлена топлотехника” НАУЧЕН РЪКОВОДИТЕЛ: проф. д-р инж. Никола Г. Калоянов РЕЦЕНЗЕНТИ: 1. доц. д-р инж. Цветан Николов Божков 2. проф. д-р инж. Иван Христов Касабов София, 2015г. ТЕХНИЧЕСКИ УНИВЕРСИТЕТ – СОФИЯ
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Енергомашиностроителен факултет
Катедра “Топлинна и хладилна техника”
маг. инж. Анатоли Венелинов Стойнов
МОДЕЛИ НА ТОПЛООБМЕННИТЕ ПРОЦЕСИ В ТЕРМОПОМПЕНИ АГРЕГАТИ
АВТОРЕФЕРАТ
на дисертация за присъждане на образователна и научна степен
“ДОКТОР”
Научна специалност
02.06.13 “Промишлена топлотехника”
НАУЧЕН РЪКОВОДИТЕЛ:
проф. д-р инж. Никола Г. Калоянов
РЕЦЕНЗЕНТИ:
1. доц. д-р инж. Цветан Николов Божков
2. проф. д-р инж. Иван Христов Касабов
София, 2015г.
ТЕХНИЧЕСКИ УНИВЕРСИТЕТ – СОФИЯ
- 2 -
Дисертационният труд е обсъден и насрочен за защита от катедра „Топлинна и хладилна техника” при Технически университет – София на заседание на катедрен съвет, проведено на 09.06.2015г.
Докторантът е зачислен в задочна докторантура към катедра „Топлинна и хладилна техника” при Технически университет – София. Изследванията по дисертационната работа са извършвани в Технически университет – София, зала 2319. Защитата на дисертационния труд ще се състои на 10.11.2015г. от 17.30 ч. в зала 2140 на Технически университет – София. Материалите по защитата се намират в канцеларията на Енергомашиностроителния факултет, стая 2344 на Технически университет – София и са на разположение на интересуващите се. Автор: маг. инж. Анатоли Венелинов Стойнов Заглавие: МОДЕЛИ НА ТОПЛООБМЕННИТЕ ПРОЦЕСИ В ТЕРМОПОМПЕНИ АГРЕГАТИ Тираж: 20 бр.
- 3 -
ВЪВЕДЕНИЕ
Няма съмнение, че термопомпеният принцип на генериране на топлина е една високоефективна алтернатива в системите за отопление, охлаждане, вентилация и загряване на вода за битови нужди. Въпреки своята многогодишна история, реализирането на този принцип в технологичен агрегат все още се свързва с конкретни въпроси, които не са намерили своето решение. Топлообменните процеси в двата основни топлообменни апарати – кондензатора и изпарителя, представляват последователност от еднофазни и двуфазни участъци, които трудно се идентифицират и описват.
ОБЩА ХАРАКТЕРИСТИКА НА ДИСЕРТАЦИЯТА
КАТЕГОРИЗИРАНЕ НА ДИСЕРТАЦИОННИЯ ТРУД
Дисертационният труд се отнася към категорията „Научноизследователски труд”, предварително отпечатан в негови съществени части. ФОРМУЛИРАНЕ НА ПРОБЛЕМА
Изследване на специфични особености при топлообмен в пластинчати топлообменни апарати в термопомпени агрегати “вода-вода”.
АКТУАЛНОСТ И МОТИВИРАНОСТ НА ПРОБЛЕМА
Потребността от адекватни и точни математични модели на топлообмен в пластинчати топлообменни апарати при фазов преход е от съществена важност – както за етапа на проектиране на агрегатите, така и за етапа на анализ на разхода на енергия при работа в условията на конкретна система. Търсенено на детайлно описание на процесите в пластинчати топлообменни апарати на термопомпени агрегати е въпрос, който се изследва от десетилетия. Изследванията най-често показват разнопосочност на получените отклонения при опитите за математическо описание и необходимост от усъвършенстване на съществуващите и търсенето на нови подходи за анализ. Така биха се разширили възможностите както при инженерни приложения, така и в научни изследвания.
ЦЕЛ И ЗАДАЧИ НА ДИСЕРТАЦИОННИЯ ТРУД
Основната цел на дисертационния труд е на основата на експериментално изследване на топлообмена получаване на обобщени математични модели на кондензатора и изпарителя като основни елементи на термопомпен агрегат.
За определеност е избрана термопомпена система “вода - вода” с най-често използвания тип топлообменник – пластинчатият.
Сложните процеси на топлообмен и описанието им налагат оборудване на експериментален стенд с измервателни уреди, планиране и осъществяване на серия от експерименти. Съществен момент в експерименталното изследване е идентификацията на зоните с и без фазов преход и намирането на локалните коефициенти на топлопреминаване.
За реализиране на тази цел се поставят и решават следните основни задачи:
1. Изучаване с техниката на изследване, анализ и изводи от източниците на информация – научни трудове, справочници и др. Това включва базите данни на онлайн библиотеките като ScienceDirect, Elsevier, Pergamon, Taylor & Francis и други – достъпни
- 4 -
библиотеки, както и публикации на фирми по теми в изследваните области.
2. Изграждане на експериментален стенд, оборудване с измервателна апаратура, калибриране на уредите за измерване, разработване на специализиран софтуер за събиране на експериментални данни с контролери,
3. Специализирано експериментално изследване на топлообмена при фазов преход на хладилен агент в пластинчати апарати.
4. Анализ на експерименталните резултати, моделиране на топлообменните процеси и обобщаване на резултатите с нови модели.
АПРОБАЦИЯ НА ДИСЕРТАЦИЯТА
Основните резултати в дисертацията са докладвани и публикувани на: 1. Списание Топлотехника, 2005 г. 2. Научна конференция „ЕМФ’08”, том I, Созопол, 2008г. 3. COFRET 2012, Франкофонски Колоквиум енергия-околна среда-икономия
Созопол, 2012г. 4. Научна конференция „ЕМФ’2015”, том I, Созопол, 2015г.
СЪДЪРЖАНИЕ НА ДИСЕРТАЦИЯТА
2. ГЛАВА 1. ЛИТЕРАТУРЕН ОБЗОР
2.1 Цикъл на Ренкин – теоретичен процес За оценка на ефективността на една термопомпена система в
идеализирания случай определящи са температурите на изпарение и кондензация θE и θC. Отнемането на топлина от нископотенциален топлинен източник и отдаването и към среда с по-висок става по обратния кръгов цикъл на Ренкин (Фиг.1.1).
КОНДЕНЗАТОР,
ИЗПАРИТЕЛ,
M
Течност Пара
Ниско
налягане
Високо
налягане
Cθ
Eθ
1
2
34
5
Фиг. 2.1.
- 5 -
Основният показател за ефективността е безразмерната величина COP (Coefficient Of Performance), отразяваща отношението на получаваната топлинна или студова мощност към необходимата за агрегата. Към необходимата мощност P за агрегата трябва да се включат освен мощността на компресора, така и необходимата мощност за транспортиране на топло/студоносителите в допълнителни елементи (помпи, вентилатори и др.).
( )P
m.ii=
P
iΔ=
P
Q=COP
rE
in,rE
out,rE -
(2.1)
2.2 “Реален” процес
Реалният процес се отличава значително от теоретичния на Ренкин. Първата съществена разлика е в пряка зависимост от топлообменните апарати – увеличаването на тяхната топлообменна повърхност извежда процесите от линията на насищане към еднофазните области, съответно към течна и парна фаза. В кондензатора това води до подохлаждане на течния хладилен агент, а в изпарителя до прегряване на парите му. И в двата случая температурата и налягането не остават постоянни. Същевременно и от компресора има загуби на топлина, както и от останалите елементи на системата към околната среда поради липса и/или невъзможност от поставяне на изолация на външните повърхности. Загуби на налягане има и в отделните линии на системата, които при дълги трасета могат да бъдат съществени. Това прави описанието в log p-i диаграма значително по-сложно, както и детайлното анализиране на процесите с количествени оценки. На Фиг. 2.2 са дадени най-съществените разлики от идеалния процес, като характерните точки са съгласно принципната схема на Фиг. 2.1.
Фиг. 2.2.
2.3 Особености при фазов преход в пластинчат топлообменен апарат (ТОА)
Сложното движение на флуидните потоци в пространствата между пластините, както и фазовия преход на хладилния агент, обуславят трудното описание на топлообмена. Известно е, че ефективни режими се осъществяват при
- 6 -
малки температурни разлики в топлообменните апарати: θext-ΔθE в изпарителя и θint+ΔθC в кондензатора. При определена мощност на компресора това се постига чрез по-голяма топлообменна повърхност или по-голям коефициент на топлопреминаване U съгласно уравнението на топлопреминаване:
(2.2)
където: e топлинната мощност,
А – площта на топлообмен,
е среднологаритмичната температурна разлика.
Горното уравнение е в сила при постоянен коефициент на топлопреминаване U и специфичен топлинен капацитет на флуидите по протежение на топлообменната повърхност, както и липса на топлинни загуби. При условия и потоци различни от дадените в дефиницията, среднологаритмичната температурна разлика се коригира с фактор по-малък от единица.
Коефициентът на топлопреминаване се дава в най-обща форма от (2.3), където hh и hc са коефициентите на топлопредаване от двете страни на стената. Съпротивлението на стената се характеризира с дебелина δ, коефициент на топлопроводност на материала на стената k и съпротивления от повърхностни замърсявания, обединени в параметъра Rfoul.
h
1+R+
k
δ+
h
1
1=U
cfoul
wall
wall
h
(2.3)
Най-често уравненията описващи коефициентите на топлопредаване са в безразмерен критериален вид (2.4), включващи критерия на скоростта Re и критерия на подобието между температурното и скоростното поле Pr . Коефициентите C1, C2 и С3 приемат най-различна форма и стойности, обусловени от типа/геометрията на граничната топлообменна повърхност и параметрите на флуидния поток.
C2 C3Nu C1.Re .Pr ....
(
(2.4)
а връзката с коефициента на топлопредаване се дефинира чрез хидравличния диаметър Dh и коефициента на топлопроводност k:
hD
k.Nuh = (2.5)
2.4 Геометрични особености на пластинчат топлообменен апарат (ТОА)
Изборът на споени пластинчати ТОА за изпарител и кондензатор е продиктувано от високите коефициенти на топлопреминаване, съчетани с ниското тегло, компактни размери, високата надеждност при различните режими и налягания. Също и приемливата цена. Според типа на профила на пластината, ТОА се подразделят основно на лентовопоточни, мрежестопоточни, “рибена кост” и множество под видове. В годините като най-популярна и успешна се е наложила шевронната форма на сечението /или “рибена кост”/ - такива са и двата топлообменни апарата, използвани в настоящата работа. Всяка пластина се отличава с характерни параметри на нагъването (Фиг.2.3).
- 7 -
Фиг. 2.3.
От гледна точка математическото описание на топлообмена в ТОА, най-важният параметър е хидравличният диаметър, фигуриращ в критериите на Нуселт и Рейнолдс. Най-често той се дефинира при две паралелни пластини с дистанция b между тях:
b.2Dh = (2.6)
2.5 Математическо описание на коефициента на топлопредаване при еднофазни потоци в пластинчат ТОА
Известни са множество зависимости с научно и инженерно-приложно значение, описващи еднофазни потоци – основно вода-вода. Основна причина за големия брой уравнения е практически многообразието комбинации от геометричните параметри, от което пък следва, че всяко конкретно изследване извежда решение за специфичен случай. Освен това в повечето изследвания рядко се дават подробни и точни геометрични данни. В настоящата работа е използвана зависимост (2.7) от изледване на пластини на шведски производител, многократно цитирана в литературата:
125,0)30(Pr
4.6
)6(Re
3,0
wall
e.3
1
2C
μ
μ.Pr.Re.1CNu
(2.7)
където константите C1 и C2 са дадени в интервали на изменение на степента на турбулизация на потока чрез критерия на Рейнолдс.
2.6 Зависимости при двуфазни потоци в пластинчат ТОА
- 8 -
Броят на известните зависимости, описващи промяната на агрегатното състояние на фреон в пластинчат ТОА, е ограничен. За разлика от еднофазните потоци, тук топлообменът зависи и от множество други термодинамични параметри като например:
паросъдържание x,
плътност на топлинния поток q,
специфичен масов поток G,
eфекти от смазочното масло за компресора, движещо се с хладилния агент,
структура на повърхността на пластината,
дебелина на граничния слой по дължина на пластината,
ламинарни, преходни и турбулентни локални режими на потока и др.
За двуфазната област при кондензация, използваната формула при анализа е на Yan-Lin:
3/1l
4.0eqph2
C Pr.Re.118,4=Nu
(2.8)
където еквивалентната стойност на Reeq e дефинирана по специфичен начин.
За двуфазната област при изпарение е използвана зависимост (2.9) на Ayub, даваща директно коефициента на топлопредаване. При описанието на експерименталните данни тя даде относително малки отклонения:
(
) *
+
(
)
[Btu/hr-ft2-◦F] (2.9)
Изводи
Сложността на математическото описание на процесите в пластинчатите топлообменни апарати в термопомпен агрегат идва от променливите коефициенти на топлопреминаване по топлообменната повърхност при изпарение и кондензация, описанието на които следва еднофазните и двуфазните области при движението на потоците. Съществуващите зависимости имат ограничена представителност, причинена от базирането им на експерименти при специфични условия, геометрии на ТОА и типове флуиди.
Към днешна дата не съществува общ подход или зависимост, описваща надеждно топлообмена в двуфазна област между хладилен агент и вода /воден разтвор/.
3. ГЛАВА 2. СЪЗДАВАНЕ НА ЕКСПЕРИМЕНТАЛЕН СТЕНД
4.1 Описание на експерименталния стенд
За целите на детайлното изследване на термопомпен агрегат “вода-вода” в лаборатория № 2319 в Технически Университет - София e проектиран и създаден специализиран стенд. Принципната схема, съставните елементи и измервателната апаратура са представени на Фиг.3.1. Основните елементи в схемата са скрол компресор тип ZR48K3E-TFD 551 на фирмата Copeland, пластинчати топлообменни апарати на фирмата Alfa Laval съответно за изпарител
- 9 -
{1} - CB52-20HX и кондензатор {8} -тип CB51-14H. Работните флуиди са вода, показана с непрекъснати цветни линии и фреон 134А.
ком
пресор
конд
енза
тор
изп
ари
тел
kWh
66C
kWh
401
°C
Bar;
°C
kW
h
8
7
5
6
4
3
1
2
9
от
вод
опровод
към
канал
10
11
12
13
14
15
16
17
Фиг. 3.1.
Реализацията на различни работни режими става с промяна на дебитите и входните температури във воден кръг, всички отчитани чрез топломер {15} тип Multical 401 и студомер {14} тип Multical 66C Ultraflow. Контролно-измервателната апаратура във фреоновия кръг се състои от четири глицеринови манометъра на вход/изход на ТОА {2}, два манометъра – вход/изход на компресора и два термометъра на изхода от ТОА част от дигитален анализатор на фреонови системи Testo 560-1 {16}, пет електронни термометъра {13} към общ преобразувател в характерните точки 1 ÷ 5 съгласно Фиг.2.1 и Фиг.2.4, пресостати ниско и високо налягане {9} и наблюдателни стъкла {5}. Данните за потреблението на електричество от компресора се визуализират на дисплея на трифазния
- 10 -
електрически анализатор Поз.{12} на Фиг. 2.4. Данните за измерваните величини постъпват чрез контролери в база данни на персонален компютър. Общият вид на стенда е даден на Фиг. 3.2.
Фиг. 3.2.
4.2 Данни за топлообменните апарати
Двата изследвани ТОА са от пластини-стомана тип AlSl316 с форма на шеврон /рибена кост/, споени с меден припой. Пластините се характеризират с геометрични данни, дадени подробно в Приложение D на дисертацията. 4.3 Точност на уредите
Точността на използваните уреди е отчетена както от техническите им данни, така и от протоколи за калибриране от специализирани лаборатории. 4.4 Измервани величини
Температура на водата вход / изход на кондензатор и
Температура на водата вход / изход на изпарителя и
Дебит на водата през кондензатор / изпарител и
Налягане/температура на кондензация / изпарение
,
, ,
Налягане изход кондензатор и вход изпарител ,
Температура на фреона вход/изход компресор
,
Температура на фреона вход/изход кондензатор ,
Температура на фреона вход/изход изпарител ,
Налягания вход/изход на компресор
,
Работни напрежения/токове на трите фази на компресора
Изводи
- 11 -
1) Окомплектован е експериментален стенд на термопомпен агрегат тип “вода-вода”.
2) Контролно-измервателната апаратура е свързана с контролери и основната част от измерванията се отчита автоматизирано чрез връзка с компютри.
3) Диапазоните на изменение на измерваните величини позволяват търсенето на представителни зависимости.
4. ГЛАВА 3. ЕКСПЕРИМЕНТАЛНО ИЗСЛЕДВАНЕ НА ТОПЛООБМЕНА В ТЕРМОПОМПЕН АГРЕГАТ “ВОДА-ВОДА”
4.1 Последователност на провеждане на експеримента
Провеждането на експеримент съдържа три периода – стартиране, нормална работа и изключване.
Стартирането включва подгряване на маслото в картера на компресора, пуск на водния кръг за стабилизиране на температурата на входящата вода и водата в системата.
Нормалната работа започва с пуск на компресора и всички измервателни уреди, както и проследяване на тяхното нормално функциониране. Изчакват се стационарни режими, с последователни промени на дебитите и температурите. Отчита се коректността на записваните величини и работните интервали.
След приключване на експеримента се изключва първо компресора, проследява се температурата на водата на изход изпарител и се спира водния кръг и записа на данни. 4.2 План на експеримента
Експерименталното изследване е проведено по смесен активно-пасивен план. Активно са променяни параметрите на потоците от страна на водните контури, а пасивно са получавани съответстващите параметри от страната на хладилния агент в съответствие с възможностите на стенда. За всяка настройка на входните параметри (дебити и/или температури на водата в кондензатора и изпарителя) режимът на работа продължава до установяване на стационарен процес. 4.3 Данни от измерванията
Данните от експеримента, отразяващи промените през интервал от 15сек., се групират в обща база от данни. Подробно са изяснени местата и начина на измерване, както и усредняване на получените стойности. Обобщена таблица от експеримента е дадена в приложение Е. Оставени са само стойностите на параметрите от стационарните режими, които са проверени с алгоритъм за достоверност в същото приложение. При математическата обработка автоматично снетите данни са усреднени от няколко последователни измервания за избягване на грешки от наблюдавани колебания, особено при отчитане на работно налягане.
5. ГЛАВА 4. АНАЛИЗ И ОБРАБОТКА НА ДАННИТЕ ОТ ЕКСПЕРИМЕНТАЛНОТО ИЗСЛЕДВАНЕ
- 12 -
5.1. Общ подход на анализ и обработка
Разглеждат се само стационарни режими, които се обработвани независимо един от друг. Построява се процеса в logP/i диаграма на работния хладилен агент – R134A. На Фиг.5.1 е даден произволен процес в диаграмата, като са означени зоните на прегряване на парите (SH), подохлаждане на течността (SQ), отнемане на топлината на прегряване (DSH) и съответните температури на водата, съответстващи на това разделение. За яснота означенията на температурите, съответстващи на състоянието на хладилния агент са удебелени и с долeн индекс “r”. Кореспондиращите им по вода са с долен индекс “w”.
Фиг.5.1
5.1.1. Процеси в топлообменните апарати – обособяване на едно и двуфазни
зони
Изчислението започва с намиране на дебитите
на фреона в изпарителя и кондензатора съгласно алгоритъм, даден в Приложение Е, след което се прави общия топлинен баланс по фреон. Еднофазните зони са описани чрез специфичните топлинни капацитети и температурите, докато двуфазните – чрез енталпиите.
(
) (5.1)
(
) (5.2)
(
) (5.3)
(
) (5.4)
(
) (5.5)
, където: Cp e специфичен топлинен капацитет i e енталпия
in/out – вход изход в съответната зона
- 13 -
1ph/2ph – едно/двуфазна зона на топлообмен.
5.1.2. Изчисление на кореспондиращите температури на водата на тези на фреона по линиите на насищане в диаграмата на процеса
Изчислението се извършва отделно за всяка температура по приложен алгоритъм в Глава 4 на дисертационния труд, съгласно означенията за температурите на водния кръг на Фиг.5.1
5.1.3. Намиране на коефициентите на топлопредаване от страната на водата и фреона за зоните без фазов преход
Тази част от изчисленията с наличните експериментални данни е извършена с критериални уравнения от различни източници, като окончателно бе използвана вече цитираната форма (2.7). В изчислителния алгоритъм даден в Глава 4, са отчетени с последователни приближения температурите от двете страни на топлообменната повърхност, като е използвана и приблизителна оценка на термичните съпротивления от отлагания.
5.1.4. Намиране на коефициентите на топлопредаване от страната на водата и фреона за зоните с фазов преход
Въпреки множеството критериални уравнения в литературата, постигането на задоволителна точност при корелирането на експериментални данни бе проблем и в при тази разработка. Общата схема на изчисление включва съвместното намиране на двата коефициента на топлопредаване с разлика в използваните критериални уравнения от страната на хладилния агент, отговарящи на промяната на агрегатното състояние в двата ТОА. Алгоритъмът следва последователността от т.5.1.3. и е даден подробно в Глава 4 от дисертацията.
5.1.5. Изчисление на произволен работен режим от експерименталните
резултати
За яснота в тази точка от основния труд подробно е дадено изчислението на един работен режим с резултатите от всички междинни изчисления. В края е оценена достоверността на получените резултати чрез сумиране на получените площи на топлообмен за съответните зони спрямо известните общи площи на топлообмен при използваните ТОА, както следва за кондензатор и изпарител:
(5.6)
(5.7)
В изчислителния алгоритъм са залагани последователно множество различни критериални уравнения за коефициента на топлопредаване в двуфазна област на фреона. Отчетени са множество детайли като загуби към околен въздух, замърсявания на топлообменната повърхност, корекция на среднологаритмичната температурна разлика и други. Получени са резултати с големи отклонения, които може да бъдат обяснени със следните причини:
- 14 -
- Разлики в работните диапазони на вариране на входните величини при експериментите въз основа на които са изведени зависимостите и данните от експерименталното изследване,
- Неточности от неясно или различно дефиниране на променливите в математическото описание,
- Липсата на оценка на влиянието на разтвореното във фреона масло при отчитане на топлофизичните характеристики като вискозитет, специфичен топлинен капацитет, коефициент на топлопроводност,
- Неясноти по отношение на формирането на граничния слой на сместа фреон-масло при различните температури и скорости на флуидите, особено предвид сложната форма на обтичаните топлообменни пластини, оттам и флуидния поток,
- Отсъствие в известните зависимости на входни фактори от съществено значение, откъдето търсенето на подходящо математическо описание е обречено на неуспех.
Получените големи отклонения, потвърждават изводите от множество резултати от литературни източници изследващи достоверността на готови критериални уравнения. Това наложи търсенето на нов подход за описание на процесите на топлообмен.
5.2. Подход за поинтервално изчисление в двуфазна област
Приложен е един възможен подход при търсенето на математическо описание - разделянето на двуфазната област на топлообмен на участъци, за всеки от които се намира локален коефициент на топлопреминаване. Прието е линейно разпределение на температурата по дължина на топлообменната повърхност, както е показано на Фиг. 5.2. за изпарителя, съответно за кондензатора.
посока на движение на водата
на
ча
ло н
а з
она
изп
ар
ени
е
кра
й н
а з
она
изп
ар
ени
е
A
посока на движение на фреона
mw
.
mr
.
площ на топлообменния апарат
тем
пе
рату
ра
Фиг. 5.2
Първоначално за всеки от интервалите е извършено изчисление идентично с това за целите двуфазни зони, като бе прието линейно разпределение на
- 15 -
паросъдържанието х. Резултатите от изчислението на отделните работни режими
показаха значителни отклонения, което наложи промяна в поинтервалния подход
Като изходна точка в описанието на интервалите в двуфазна област са приети известните начални условия - температури на водата и фреона по границите на двуфазната област, която се разделя на осем равни по площ интервала. Също частта от общата площ, съответстваща на двуфазен топлообмен, начални и крайни енталпии, дебит и общ топлинен поток. В двата ТОА изчислението започва от паросъдържание х=1.
По алгоритъм с последователни приближения се намират следните параметри: крайна енталпия, топлинен поток по фреон, температура на водата в края на интервала, коефициентите на топлопредаване чрез критериални уравнения, коефициента на топлопреминаване и се проверява топлинния поток за предварително зададен праг на отклонение. Характерно е използването на средноинтервалното паросъдържание в критериалните уравнения за хладилния агент. За начални параметри на всеки следващ интервал се използват граничните стойности от предходния. Стойността на паросъдържанието в последния осми интервал трябва да има стойност 0 при кондензация, а при изпарителя - съответната отчетена стойност при процеса изпарение. Сходна проверка се извършва и при сравнение между получената в последния интервал температура на водата и очакваната от разделянето по зони на едно и двуфазен топлообмен.
Поради малките разлики в поинтервалните температури не се налага корекция на среднологаритмичната температурна разлика.
Резултатите от примерен режим са показани на Фиг.5.3. за кондензация и Фиг.5.4. за изпарение.
Фиг.5.3. Идентифицирана е зависимостта на коефициента на топлопреминаване /показан със синя линия с ромбове/ от паросъдържанието, като тя е в права пропорционалност и при двата двуфазни режима. Разлика се наблюдава в характера на изменение по зони, като това е резултат от конкретните използвани критериални уравнения и формата, под която участва в тях паросъдържанието.
Прилагането на поинтервалния алгоритъм в този му вид върху цялата експериментална база данни още веднъж потвърди лошите оценки за приложимостта на критериалните уравнения в изследвания ТОА и работните диапазони на входните фактори.
- 16 -
Фиг.5.4.
Изводи
1) Съставен е автоматизиран изчислителен алгоритъм в средата Excel, използващ експерименталната база данни и библиотеки с топлофизичните характеристики на използваните флуиди.
2) Доказани са незадоволителните резултати от приложението на известните критериални уравнения за отделните едно и двуфазни зони на топлообмен, когато се третират като една зона. Основен проблем представлява точността при описанието на топлопредаването в двуфазните зони от страна на хладилния агент.
3) Представен е един детайлен подход за поинтервално изчисление на коефициента на топлопреминаване в двуфазните зони.
6. ГЛАВА 5. ПРЕДЛОЖЕНИЕ ЗА ОБОБЩЕНО ОПИСАНИЕ НА ТОПЛООБМЕННИТЕ ПРОЦЕСИ В КОНДЕНЗАТОРА И ИЗПАРИТЕЛЯ
6.1. Коефициенти на топлопреминаване
За определящи температури в кондензатора са приети температурата на
кондензация според налягането на фреона и температурата на
постъпващата вода в него . Избраната математическа форма за
интерпретиране на обобщения топлообменен процес в кондензатора е:
θ
θ
θ
(
)
(6.1)
и съответно за изпарителя:
θ
θ
(
)
(6.2)
Така дефинираните обобщените коефициенти са за “базова” или приведена температурна разлика, затова имат характер на “приведени” инженерно-приложни коефициенти. В таблици 6.1. и 6.2 са илюстрирани част от експерименталните данни , използвани за идентификация на обобщените модели на коефициента на топлопреминаване.
6.2. Математически модел на коефициентите на топлопреминаване
За идентификация на обобщените модели на коефициента на топлопреминаване в кондензатора и изпарителя на термопомпения агрегат е приета следната хипотеза за зависимост:
за кондензация (температурата на фреона в двуфазна област остава постоянна):
(
) (6.3)
за изпарение (поради пад на налягане температурата на фреона в двуфазна област намалява):
(
) (6.4)
От множеството възможни апроксимиращи функции следните две – (6.5) и (6.6) са намерени като най- подходящи за обобщаване на процесите в кондензатора и изпарителя:
(
)
(6.5)
(
)
(6.6)
Оценката за адекватността и точността на моделите е извършена по правилата на регресионния анализ. Резултатите за кондензатора са илюстрирани на Фигури 6.1 и 6.2.
Оценка на резултатите за обобщения коефициент на топлопреминаване в режим изпарение
Y Predicted Y
0%
20%
40%
60%
80%
100%
120%
0
5
10
15
20
25
-8 -6 -4 -2 0 2 4 6 8 More
Че
сто
та
Отклонение
Frequency
Cumulative %
- 21 -
6.3. Оценка на влиянието на получените обобщени модели върху изчислителните мощности на кондензатора и изпарителя
Достоверността е оценена и чрез влиянието на моделираните коефициенти на
топлопреминаване върху топлинната мощност на кондензатора , съответно
хладилната мощност на изпарителя . Това е направено въз основа на приетата в т.6.1. хипотеза и следващите я зависимости. Графичният израз е даден на хистограмите на разпределението на отклоненията – Фигури 6.5 и 6.6.
Фигура 6.5
Фигура 6.6
Изводи
Генерираният метод за определяне на обобщени коефициенти на топлопреминаване, описващи сложния топлообмен при кондензация и изпарение на хладилен агент в пластинчати ТОА е приложим за инженерни пресмятания с достатъчна точност.
0%
20%
40%
60%
80%
100%
120%
0
5
10
15
20
25
30
-8 -6 -4 -2 0 2 4 6 8 More
Че
сто
та
Отклонение
Series1
Series2
0%
20%
40%
60%
80%
100%
120%
0
5
10
15
20
25
-8 -6 -4 -2 0 2 4 6 8
Че
сто
та
Отклонение
Frequency
Cumulative %
- 22 -
7. Математически модел на работата на компресор като част от система
7.1. Постановка при математическо описание на работата на компресор
Стремежът за постигане на оптимални работни режими на термопомпени системи изисква и подробен анализ на енергийните характеристики на компресора. За тази цел е необходимо разполагане с опростен, но с достатъчна точност модел на компресора, който да включва съвместими с използваните в обобщените модели на коефициента на топлопреминаване параметри.
Такъв модел е потърсен въз основа на каталожни данни при различни температури на изпарение и кондензация:
(7.1)
където: C0 ÷ C9 са коефициенти,
Q – хладилна мощност /аналогична е при консумирана мощност P/
– съответно температура на кондензация и изпарение
При такова представяне липсва възможност за отчитане влиянието на подохлаждането на хладилния агент изхода на кондензатора и прегряването на парите на изхода на изпарителя, имащи значително влияние върху общата производителност. В опит да се избегне този недостатък е потърсена корелация, включваща и температурите на подохлаждане и прегряване въз основа на статистическа база данни за използвания в експерименталния стенд скрол компресор Copeland ZR48K3E-TFD 551 в следните интервали на входните фактори и изходния параметър:
Температура на изпарение на хладилния агент -20 до 12.5°С
Температура на кондензация на хладилния агент 30 до 70°С
Температурна разлика при прегряване на парите на хладилния агент в изпарителя
3 до 10 K
Температурна разлика при подохлаждане на кондензата на хладилния агент в кондензатора
0 до 8 K
Хладилна мощност на компресора 2250 до 11750W
7.2. Резултати от моделирането на хладилната мощност
Търсенето на подходящи математически описания се сведе до два подхода. При първия бе потърсена аналогия с уравненията описващи топлообменните апарати, предвид сравнително лесната и бърза употреба. При
втория множество частни зависимости при различни и
бяха обобщени в
една - сложна и обемиста, подходяща за вграждане в изчислителен алгоритъм.
7.2.1. Опростен общ модел
По метода на най-малките квадрати са потърсени коефициентите в опростената, оценена като инженерно приложима форма. Резултатите от регресионния анализ показаха коефициент на множествена корелация R²=0,99. Числените стойности на коефициента и степенните показатели са посочени в полученото уравнение:
(7.2)
- 23 -
Оценката за адекватност и точност e направена чрез процентните отклонения на изчислените стойности от техните оценки съгласно:
(
)
( )
(7.3)
Резултатите от индивидуалните оценки спрямо изходната база данни са обобщени в следната хистограма:
Фигура 7.1
7.2.2. Обобщен модел на основата на множество частни модели
Основната част от общата база данни представлява групирането на данни за търсената хладилна мощност Q при съответни температури на изпарение и
кондензация и
. Всяка от тях е при постоянни стойности за прегряване и
подохлаждане и
. За всеки от коефициентите от до се търси
зависимост от прегряването и подохлаждането. Дори прости линейни уравнения дадоха задоволителна точност. За коефициентите в уравнение (7.1) се получи:
(7.4)
Оценката на достоверността на резултатите спрямо всички изходни данни съгласно уравнение (7.3) , показващо процентното отклонение за хладилната мощност Q, е показана на Фиг. 7.2.
0%
20%
40%
60%
80%
100%
120%
0
200
400
600
800
1000
1200
-8 -6 -4 -2 0 2 4 6 8 More
Че
сто
та
Процентно отклонение
_ Frequency
Cumulative %
- 24 -
Фигура 7.2
Изводи
1) Грешката в рамките на 2% потвърждава възможността за точно описание на характеристиките на компресора чрез сложна зависимост на база множество зависимости при различни стойности на температурите на прегряване в изпарителя и подохлаждане в кондензатора.
2) Лесната и бърза оценка с общия модел (5) даде отклонения в рамките на 5%, достатъчни при случаи, в които не се изисква висока точност.
ЗАКЛЮЧЕНИЕ
В резултат на извършеното експериментално изследване се получиха резултати с
научноприложни приноси, които може да се обобщят както следва:
Формулирана е и е доказана нова хипотеза:
Може да се намери обобщена енергийна характеристика на сложните многосвързани едно и двуфазни топлопреносни процеси в двата основни топлообменни апарата на термопомпен агрегат – кондензатор и изпарител.
Приложен е нов метод на изследване:
За идентифициране на локалните коефициенти на топлопреминаване в двуфазна област е приложен комбиниран подход на експериментално измерване и математическо моделиране на температурното поле с разделяне по интервали.
Получени са потвърдителни факти:
Потвърди се, че общоприетият подход за използване на известните критериални уравнения за среден коефициент на топлопредаване в двуфазна област не отчита спецификата на сложните топлообменни процеси и води до значими отклонения в резултиращата стойност на коефициента на топлопреминаване от реално постигания. По тази причина за преследваните цели е разработен подход за намиране на локални коефициенти на топлопреминаване.
0%
20%
40%
60%
80%
100%
120%
0
500
1000
1500
2000
2500
3000
-8 -6 -4 -2 0 2 4 6 8 More
Че
сто
та
Процентно отклонение
- 25 -
На основата на сравнителен анализ е потвърдена невъзможността за прилагане/използване на известните в литературата критериални зависимости поради факта, че при двуфазни топлообменни процеси с хладилен агент R134A паросъдържанието не е единственият значим фактор, определящ големината на локалните коефициенти на топлопредаване. Съществено влияе и разтвореното във фреона полиолестерно масло, хомогенизирането на потоците с подходящи елементи на входните фланци на ТОА и др.
Резултати, потвърждаващи правомерността, достоверността и точността на хипотезата:
Получени са обобщени стойности за кондензатор и за изпарител както следва:
За кондензатора – обобщен коефициент на топлопреминаване за процесите на отнемане топлината на прегряване, същинска кондензация и подохлаждане на кондензата.
За изпарителя - обобщен коефициент на топлопреминаване за процесите на същинско изпарение и прегряване на парите.
Доказана е строга корелация на обобщените характеристики с режимните параметри както следва:
За кондензатора – дебити на фреона и водата, температурата на същинска кондензация и температурата на водата на входа. Получена е степенна зависимост, която с достатъчна за инженерна приложимост точност ±5% описва експерименталните данни.
За изпарителя - дебити на фреона и водата, температурата на изхода на същинско изпарение и температурата на водата на входа. Получена е степенна зависимост, описваща с достатъчна за инженерна приложимост точност от ±5% експерименталните данни.
Обобщените коефициенти са за “базова” или приведена температурна разлика, затова имат характер на “приведени” инженерно-приложни.
Получен е обобщен модел на енергийните характеристики на компресор в термопомпена система.
Получените резултати може да се ползват при решаването на всяка от двете задачи: проектна и проверочна.
СПИСЪК НА ПУБЛИКАЦИИТЕ ПО ДИСЕРТАЦИОННИЯ ТРУД
1. А. Стойнов. Алгоритъм за моделиране и симулиране на хладилни компресори, Списание Топлотехника, 2005.
2. Калоянов Н., А. Стойнов. Експериментален стенд за изследване на енергийните характеристики на термопомпен агрегат „вода-вода”. Сборник Доклади на XIII-та Научна Конференция с Международно Участие ЕМФ’08, Созопол, Том I, 2008.
3. Kaloyanov N., A. Stoynov. An Approach for developing of generalized models of evaporation and condensation in plate heat exchangers as part of “water-water“ heat pump system. Proceedings of the Conference COFRET 2012, 11-13 june 2012 Sozopol, 2012.
4. Калоянов Н., А. Стойнов Математически модел на работата на хладилен компресор като част от термопомпена система “вода-вода“ – приета за печат
- 26 -
в Сборник Доклади на XX-та Научна конференция с международно участие „ЕМФ 2015”
- 27 -
ABSTRACT OF THE THESIS
Тitle: “Models of the heat exchange in “water to water” heat-pump systems”. The main goal is development of adequate models of the heat exchange in the condenser and the evaporator, based on the experimental study. It’s essential to identify the heat exchange areas with one-phase and two-phase flows and their heat transfer coefficients. Considering this common mathematical models of the heat transfer coefficients have been identified. Chapter 1 presents analysis of the existing techniques of representing the 1 and 2-phase heat exchange in brased plate heat exchangers. Chapter 2 presents the development and establishing of an experimental stand, equipment, and measuring tools necessary for the experimental study. Chapter 3 is dedicated to experimental plan and collected data. Chapter 4 presents developed independent models of the heat transfer in both 1 and 2-phase zones as well as an approach for extended models of the 2-phase zone and their verification. A proposal for obtaining an overall heat transfer coefficients, representing the whole condenser and evaporator is presented in Chapter 5. As well as a simple and detailed models of the compressor have been developed, analyzed and presented.