CC BY
11
HEAT PUMP GAS COOLER CONTROL USING CRITERION OF MINIMUM OF EXERGY LOSSES
Sit M.L., Sit B.M., Burciu V.I., Ioi§er A.M.
Institute of Power Engineering of the Academy of Sciences of Moldova
Abstract. This paper focuses on the development of the criterion of optimality of transients of the control system, based on the minimum of exergy losses in the gas cooler of carbon dioxide heat pump. It is noted that the exergy quality criterion has a clear physical meaning, as compared with the integral quadratic criterion in which the choice of the coefficients in the integrand is not justified. Mathematic model of heat exchanger is obtained using the method of solving differential equations, without going to the irrational transfer functions. The model is reduced to transfer functions of the first and second order with the delay. The continuous temperature control system of heat pump gas cooler is considered. It is shown, that one of the versions of the control system for the minimization of the proposed criterion can be a combined control system using both the principle of the negative feedback and the principle of the invariance related to a number of disturbances affecting the processes of heat transfer in the heat exchanger.
Keywords: heat pump, carbon dioxide, control system, gas cooler, exergy, heat exchanger.
DIRIJAREA RACITORULUI DE GAZE AL POMPEI DE CALDURA DUPA CRITERIUL DE MINIMUM AL PIERDERILOR DE EXERGIE §it M.L., §it B.M., Burciu V.I., Ioi^er A.M.
Institutul de Energetica al Academiei de §tiinte a Moldovei Rezumat. Se propune de a utiliza ca criteriu de optimizare a proceselor tranzitorii in sistemul de dirijare valoarea minimala a pierderilor exergiei in procesul tranzitoriu. Aceasta se prezinta de catre functia de sub integrala a relatiei integrale patrate a criteriului de calitate a procesului tranzitoriu din sistemul de racire cu gaze a pompei de caldura cu dioxid de carbon. S-a demonstrat, ca pentru minimizarea in baza criteriului propus, ca varianta de realizare a sistemului de dirijare poate servi sistemul combinat de dirijare, care include concomitent principiul reactiei inverse negative §i principiul de invarianta a reactiei pentru unele perturbatii, care au un impact asupra proceselor de schimb de caldura in schimbatorul de caldura. Se mentioneaza, ca criteriul exergetic al calitatii poseda un sens evident fizic in comparatie cu criteriul integral patratic, deoarece in ultimul criteriu selectarea valorilor coeficientilor din expresia de sub semnul integralei nu este argumentata. Prin utilizarea metodei de obtinere a solutiilor sistemului de ecuatii cu derivate partiale §i fara a utiliza procedura de trecere la functii de transfer irationale, s-a obtinut modelul matematic al sistemului de racire cu gaze a pompei termice cu dioxid de carbon pentru ciclul supercritic de functionare. Modelul este prezentat de functii de transfer de primul §i al doilea ordin cu intarziere.
Cuvinte-cheie: pompa de caldura, dioxid de carbon, sistemul de dirijare, exergia, schimbator de caldura.
УПРАВЛЕНИЕ ГАЗООХЛАДИТЕЛЕМ ТЕПЛОВОГО НАСОСА ПО КРИТЕРИЮ МИНИМУМА ПОТЕРЬ ЭКСЕРГИИ Шит М.Л., Шит Б.М., Бурчиу В.И., Иойшер А.М.
Институт Энергетики АН Молдовы Аннотация. Предложено использовать минимум потерь эксергии в качестве подынтегрального выражения интегрального квадратического критерия качества переходного процесса системы управления газоохладителем теплового насоса на диоксиде углерода. Показано, что для минимизации предложенного критерия одним из вариантов системы управления может быть комбинированная система управления, использующая как принцип отрицательной обратной связи, так и принцип инвариантности по отношению к некоторым возмущениям, действующим на процессы теплопередачи в теплообменнике. Отмечено, что эксергетический критерий качества имеет ясный физический смысл, по сравнению с интегральным квадратическим критерием, где часто выбор коэффициентов в подынтегральном выражении не обоснован. С использованием метода решения системы дифференциальных уравнений в частных производных, без перехода к иррациональным передаточным функциям, получена модель газоохладителя теплового
насоса на диоксиде углерода, работающего в сверхкритическом цикле. Модель приведена к передаточным функциям первого и второго порядка с запаздыванием.
Ключевые слова: тепловой насос, диоксид углерода, система управления, эксергия, теплообменник.
Условные обозначения: Т - температура газа после газоохладителя, Т2 - температура воды после газоохладителя, Тс - температура стенки, Л - интервал времени дискретизации, Лх — шаг дискретизации по длине газоохладителя (ГО), с — теплоемкость газа, шх — массовый расход газа, V] — скорость газа, аг — коэффициент теплоотдачи от газа к стенке, Тм — задание по температуре нагреваемого теплового агента.
1. Введение
Эксергетический метод является теоретической базой анализа тепловых процессов и установок. Применению этого метода посвящена обширная литература, например [1,2,3,4], обзор литературы [5]. В работе [4] показано, какие условия должны быть выполнены для оптимальной организации теплообмена для условий линейного закона те-плопереноса. Предложена структурная схема системы автоматического управления теплообменником, при которой производство энтропии в такой системе будет минимальным. Система стабилизирует температуру горячего потока изменением температуры на входе первичного потока и поддержанием отношения водяных эквивалентов первичного и вторичного потоков [4]. Примерами работ, где исследуется деструкция эк-сергии в статических режимах работы теплообменников, являются работы [6-8]. В работах [9,10] рассмотрены вопросы применения эксергии в исследовании динамики систем управления дистилляционными колонами. Вопросы, связанные с получением эк-сергетического критерия качества для оценки качества переходного процесса системы управления теплообменником и обоснования структуры системы управления им, в частности, газоохладителем теплового насоса, до настоящего времени не рассматривались. В работе [11], были предложены структурные схемы тепловых насосов на диоксиде углерода. Был рассчитан газоохладитель, выполненный на базе кожухозмеевикого противоточного теплообменника и была решена задача определения переходных характеристик газоохладителя на диоксиде углерода на основании численного решения системы уравнений динамики с переменными коэффициентами. В данной работе показана связь между моделью газоохладителя, как объекта управления в системе управления температурой нагреваемой среды и влиянием эксергетического критерия качества на структуру системы управления.
2. Математическая модель газоохладителя
Модели теплообменников «газ-жидкость», как объектов с распределенными параметрами исследовалась в работах В.Н. Девятова [12], А.А.Шевякова [13], однако, в их работах рассматривались теплообменники с параметрами теплоносителя, которые практически постоянны. Нами рассматривается работа теплообменника в сверхкритическом цикле, когда его параметры изменяются в зависимости от температур, давлений и расходов теплоносителей. Для определения статического режима работы газоохладителя была использована методика, разработанная И.В. Деревичем и Е.Г. Смирновой
[14].
Для определения инерционности рассматриваемого газоохладителя (противо-точного теплообменника «газ-жидкость» без учета тепловой инерционности стенки) использованы уравнения [15,16].
lt + V1 ^ = ^12(V1' V2' C1'A, С2' р2) ~ T1)
T dt
dT1 dx
= K13(vb ^ cb P ^ P2)-(T1 - T2 )
(1) (2)
Газоохладитель ТНУ при работе в сверхкритическом цикле представляет собой динамическое звено с переменными параметрами. При необходимости учета теплоемкости стенки газоохладитель, может быть описан системой из трех дифференциальных уравнений в частных производных:
dT
dT
G1 ~T+cm ~T =a1f1 iTc- T);
dt
dx
dT
CcGc~& = af2 ' (T2 - TC )- ^ ' (TC - T1 ) '
C1G2 ^ - C2m2 dT2 = a2f2 ■ (T2 - TC )
(3)
при граничных условиях Т (0, ^) = То, Т ^) = То и нулевых начальных условиях. Выходными величинами объекта управления являются температура воды после газоохладителя - Т2 и температура газа после ГО - Т • Управляющими воздействиями являются температура рабочего тела Т (0, t), и скорость его движения у1 . Зависимости коэффициента теплоотдачи от длины теплообменника приведены на рис.1.
^ 2000 -И
я 1800
и 1600 &
-ГАЗ-СТЕНКА ■'СТЕНКА-ВОДА
10 15 20
Длина теплообменника, м
Рис. 1. Зависимость коэффициентов теплоотдачи на единицу длины трубки от длины
теплообменника
c
2200
S 1400
200
н 1000
800
600
0
5
25
30
Отличием от известных работ является то, что в рассматриваемом случае в уравнениях (1) и (2) коэффициент а зависит от значений величин c, m, V • Решение этой системы уравнений может быть найдено методами приближенного интегрирования дифференциальных уравнений в частных производных. Были сделаны следующие предпосылки при решении системы уравнений (1):
1) потерь тепла в окружающее пространство нет; 2) скорость теплопередачи по всей длине аппарата пропорциональна разности первых степеней температур сред; 3) периметр поперечного сечения поверхности раздела сред постоянен по всей длине аппарата; 4) смешение теплоносителя в направлении движения сред незначительно и не учитывается; 5) тепловая емкость стенки мала по сравнению с тепловой емкостью движущихся масс обменивающихся теплотой сред и не учитывается; 6) свойства диоксида углерода изменяются по длине аппарата в зависимости от температуры и давления газа.
Именно п.6 является «камнем преткновения» при определении динамических характеристик теплообменника. В этом случае «не работают» методы передаточных функций, которыми обычно аппроксимируют систему (1) и (2). Решение задачи осложняется также и тем, что скорость газа намного превосходит скорость жидкости, и система уравнений (1) и (2) становится «жесткой». При этом способ решения приходится искать среди множества методов решения гиперболических систем дифференциальных уравнений в частных производных. Эта система нами решалась с использованием одного из явных методов Эйлера. Здесь частная производная по пространству аппроксимирована с помощью центральной разности первого порядка, а частная производная по времени аппроксимирована разности вперед.
Эта схема устойчива и имеет погрешность аппроксимации O(At, (Дх2)) .
= Tj - dt
1,i , 2 • dx
T/+1 = Ti, + dt
2,i , 2 • dx
Vi(j-1,0[j -Tj] + dt• Kl2(j-1,0(Tj -Tj)
Tg = Tj +—• v2(j -1,0 [j - Tj-1 ]+dt • Ku(j -1, о (T1j - Tj)
(4)
Структурная схема газоохладителя ТНУ, как объекта управления представлена на рис.2.
Рис. 2. Структурная схема газоохладителя теплового насоса как объекта управления
Вид передаточных функций и числовые значения коэффициентов для кожухоз-меевикового газоохладителя, состоящего из семи параллельно включенных по газу и
3 3
воде блоков, при общем расходе газа 5,5 м/с и расходе воды 5, 84 м /час и давлении газа 9 МПа, температуре газа от 30 оС до 80оС в каждом из которых длина трубки (материал - сталь) составляет 30м., внутренний диаметр трубки 0,01м., толщина стенки трубки 0, 001 м, диаметр навивки трубки - 0,25м. представлены в таблице 1. Необходимо сразу оговориться, что значения коэффициентов существенно зависят от режима работы ТНУ и значительно изменяются в зависимости от параметров цикла ТНУ.
Таблица 1. Вид упрощенных передаточных функций газоохладителя
Передаточная функция Значения коэффициентов.
к р"рг / wii(P) = V /тир+1 Ти = 3,1;г = 0,8с; / = 0,01
W12( Р) = /4 Р + 1)(Тз1 Р + 1) Тп = 1,1; Т12 = 0,22; кх = 0,368.
W21( р) = ¡/(Т1р + 1)(Т 2 р +1) к = 0,95; Т1 = 3,81; Т2 = 0,79.
к\е рт / W22( р) = 1 /(тп р + 1)(Т21 р +1) / = 0,33; г = 6,4Т = 6,5;Т2 = 2,3.
W31( р) = р + 1)(Тз2 р + 1) / = 2,5;Т31 = 16,8;Т32 = 0,74.
W32( р ) =% р + 1 кх = 5,04; Г32 = 2,2.
W41( р) = /{«р + 1)(Т42 р +1) к41 = 0,82; Т41 = 6,26;Т42 = 6,4.
W42( р) = /{Т4з р + 1)(Т44 р + 1) к42 = 11,17;Т43 = 5,35;Т44 = 0,76.
3. Описание закона управления
Критерием качества работы теплонасосной установки является обеспечение максимального значения коэффициента тепловой эффективности (СОР), а критерием качества контура управления температурой нагреваемого теплового агента является минимум деструкции эксергии в теплообменнике в процессе регулирования .
Для системы управления газоохладителем в динамике это означает требование апериодического переходного процесса сигнала температуры воды при минимальном времени переходного процесса. Управлениями ГО являются: следующие параметры: температура и давление газа на входе в теплообменник, расход газа через теплообменник. Возмущения - температура и расход воды на входе в теплообменник. Из-за переменности параметров передаточных функций ГО в зависимости от режима работы система выполнена комбинированной, т.е. использует принцип обратной связи и принцип инвариантности.
Результаты численного решения уравнений динамики газоохладителя, записанных в виде системы дифференциальных уравнений с распределенными параметрами, полностью учитывают динамику объекта. При этом подходе, который мы назовем агрегированным, предварительный выбор коэффициентов регулятора является более обоснованным, чем при использовании других методов, так как отсутствуют промежуточные преобразования передаточных функций, их упрощение и неизбежные при этом погрешности. В полученной данным методом передаточной функции коэффициенты
зависят от свойств газа, которые являются переменными (рассматриваемый перегреватель газа работает в широком диапазоне давлений).
Рассмотрим динамику потерь эксергии в теплообменнике. Уравнение эксергетического баланса теплообменника в статике имеет вид [17]:
(в11 " в21 )• О1 = ( в22 " в12 )• О2 0 (5) или ^ 0 = ( в11 - в21 )• О1 — (в22 - в12 )• О2 (6)
В уравнениях (5), (6) индексы при расходах сред 1,2 относятся соответственно к хладагенту и нагреваемому теплоносителю, а цифровые индексы - к выходам и входам хладагента и теплоносителя (хладагент охлаждается от состояния ви до состояния в21, а тепловой агент нагревается от состояния в12 до состояния в22). Зависимость эксергии от температуры для рассматриваемого примера имеет вид:
в1(2)1 = ^1(2)1 — — Та (^1(2)1 — ^вхг) , (7)
где. Нех,,— энтальпия и энтропия окружающей среды, 1 = 1,2. Так как, энтропия ^ и энтальпия \ диоксида углерода на выходе из газоохладителя, при малых отклонениях от установившегося режима работы теплообменника, линейно зависят от температуры и давления хладагента, то можно написать, что потери эксергии:
0(Г) = к1 • О10 •&Т1 Г) + к2 • О20 • ЛТ2 Г) + (вцо — в210 ) • АО1 (т) + (в120 — в2м ) • ЛО2(т) + к3 • О • Щ(Т), (8)
где, приращения температур хладагента АТ (т) и нагреваемой среды ЛТ2 (т), давления хладагента Ар (т) рассматривают в зависимости от времени, к1,к2 — коэффициенты, характеризующие линеаризованную зависимость энтропии от температуры. Индексы «0» при переменных, входящих в (8), обозначают установившийся режим работы. Эксерге-тическими потерями из-за потерь давления в теплообменнике пренебрегаем. Интегральный квадратический критерий качества для оптимизации САУ газоохладителем по критерию минимума деструкции эксергии после ряда упрощений (пренебрежения членами выражения (Л^ О(т))2), кроме квадратических) будет иметь вид:
!
1 = 1
(к, • Ою • АТ1(т))2 +(к2 • Ом • ЛТ2(т))2 +(кз • Ою • Ар(т))2 +... ... + ((«110 — в210 )-А°1(г)) +((в 120 в220 ) • ЛО2 (т))
• йт-
->Ш1И
(9)
0
Из вида уравнения (8) следует, что для минимизации потерь эксергии при регулировании процесса теплообмена в теплообменнике, система управления должна минимизировать колебания температуры и давления хладагента в газоохладителе, а также и расхода нагреваемой среды (при регулировании температуры нагреваемой среды изменением расхода хладагента). Из вида уравнения (8) следует, что система управления теплообменником, например, газоохладителем теплового насоса, спроектированная по критерию минимума потерь эксергии, может включать в себя как контуры компенсации возмущений, так и главный контур управления, построенный на основании использо-
вания принципа обратной связи. Значение величины давления на входе в газоохладитель теплового насоса стабилизируется посредством изменения перепада давлений на регулирующих вентилях теплового насоса (так учитывается параметр Д^ (г) в формуле (6)), температура хладагента ДГ (г) регулируется изменением режима работы теплообменника, установленного перед компрессором. Изменение расхода хладагента G (г) регулируется в зависимости от нагрузки теплового насоса ДG2 (г). В приведенных на рис.3 структурных схемах контуров системы регулирования температуры на выходе газоохладителя управлениями являются T (г), P (г) - давление и температура газа перед газоохладителем, т.е. - управления, которые обычно присутствуют в классических критериях оптимальности систем управления. Преимуществом эксергетического критерия качества перед «классическим» квадратическим критерием вида J = J u 2(z)dz является
ясный смысл коэффициентов перед величинами управлений, входящих в подынтегральное выражение, что позволяет обоснованно выработать требования к качеству переходных процессов в системе управления.
G 2
IN
Рис. 3. Структурная схема контуров регулирования температуры теплоносителя
Для системы управления (рис.3) регулятор в контуре с обратной связью может быть выбран из класса ПИД-регуляторов с переменными коэффициентами. На рис.3 T2m,T2OUT,T2out pr - температура воды на входе в ГО, на выходе из ГО, заданная температура воды, РО1- регулирующие органы (система регулирования давления компрессора, система регулирования температуры газа перед компрессором) передаточная функция перегревателя газа, С1 - контроллер, ОУ- объект управления, WK1 - компенсатор, РО2- регулирующий орган расхода хладагента (регулируемый привод компрессора, регулирующий вентиль). Можно показать, что, что при постоянных настройках ПИД-регулятора внутреннего контура и при переменных характеристиках объекта управления в переходном процессе появляется перерегулирование, что является нежелательным. Для устранения этого необходимо использовать ПИД-регулятор с параметрами настройки, изменяемыми в зависимости от параметров объекта управления. Также можно показать, что компенсатор является реализуемым. Компенсация запаздывания обеспечивается путем установки измерителя температуры воды перед ГО на необходимом расстоянии с учетом максимального расхода воды и введении в схему регулируемого запаздывания. Структура блока компенсации возмущений по расходу воды аналогична структуре схемы, рис. 3, и отличается от нее только тем, что управлением является расход хладагента - диоксида углерода. В целом, структура системы аналогична, описанной в [4] с коррекцией на условия работы газоохладителя теплового насоса.
4. Заключение
Предложена методика проектирования системы управления газоохладителем теплового насоса, оптимальной по критерию минимума потерь эксергии. Отмечено, что одним из вариантов такой системы может быть комбинированная система с тремя управлениями, использующая как принцип обратной связи, так и принцип инвариантности по отношению к возмущениям. Отмечено, что эксергетический критерий качества имеет ясный физический смысл, по сравнению с интегральным энергетическим квадрати-ческим критерием, где выбор коэффициентов в подынтегральном выражении часто трудно обосновать.
Литература
[1] Rozonoer L.I., Tsirlin А.М. "Optimalinoe upravlenie termodinamicheskimi protse-sami. I. Avtomatika I telemehanika, Nr.1, 1983, s. 70-79. (in Russian)
[2] Rozonoer L.I., Tsirlin А.М. "Optimalinoe upravlenie termodinamicheskimi protse-sami. II. Avtomatika I telemehanika, Nr.2, 1983, s. 88-101. (in Russian)
[3] Rozonoer L.I., Tsirlin А.М. "Optimalinoe upravlenie termodinamicheskimi protse-sami. III. Avtomatika I telemehanika, Nr.3, 1983, s. 50-64. (in Russian)
[4] Tsirlin А.М. " Optimalinoe upravlenie protsesami neobratimogo teplo - i massope-renosa". Tehnichescaia kibernetika, №2, 1991, s.171-179. (in Russian)
[5] Orlov V.N., Rudenco А.У. " Optimalinoe upravlenie v zadachah o predelinyh voz-mojnosteah neobratimyh termodinamicheskih processov". Avtomatika I telemehanika Nr.5, 1985, s. 7-41. (in Russian)
[6] Kazakov V.G., Lukanin G.P., Smirnova О.Б. Eksergeticheskie metody оцotsenki teplotehnologhicheskih ustanovok: Uchebnoe posobie. Sankt-Peterburg, 2013.- 62s. (in Russian)
[8] [9]
McGovern J., Tirca-Dragomirescu G., Feidt M., Dobrovicescu A. A critical Exploration of the Usefulness of Rational Efficiency as a Performance Parameter for Heat Exchangers. Proceedings of ECOS 2012 - the 25th international conference on efficiency, cost, optimization, simulation and environmental impact of energy systems June 26-29, 2012, vol.2, pp.14-23, Perugia, Italy.
Madara§an T., Ungure§an P. Contributions regarding the exergetic analzsis of heat exchangers. Termotehnica, Nr.1-2, 2006, pp.25-27.
Munir M.T., Yu W. and Young B.R. Control loop configuration and eco-efficiency.
http://www.google.md/url?sa=t&rct=i&q=&esrc=s&source=web&cd=10&cad=ria &uact=8&ved=0CF4QFjAJ&url=http%3A%2F%2Fwww.nt.ntnu.no%2Fusers%2F skoge%2Fprost%2Fproceedings%2Fcpc8-focapo-
2012%2Fdata%2Fpapers%2F025.pdf&ei=GHBcVM7VEcriaNmbgKgM&usg=AF QiCNFdBGWia0eZK1JAQ IGEF20fQIyZQ&bvm=bv.79184187,d.d2s (дата обращения 22.10.2014)
Munir M.T., Yu W., Young B.R. Eco - efficiency of control configurations using exergy. UKACC International Conference on Control 2012 Cardiff, UK, 3-5 September 2012, pp. 160-165. DOI 10.1109/C0NTR0L2012.6334623. Sit B.M., Juravliov A.A., Sit M.L., Povyshenie energheticheskoi effectivnosti tep-lonasosnoi ustanovki na dioxide ugleroda v sisteme kombinirovannogo teplosnab-jenia za schet optimizatsii eio sistmy upravlenia. Problemy reghionalinoi energetiki, N1, 2008, http://iournal.ie.asm.md/assets/files/m71 2 61.doc. (in Russian) Deviatov B.N. «Teoria perehodnyh protsesov v tehnologhicheskih apparatah s tochki zrenia zadach upravlenia». Novosibirsk, Nauka, 1964. (in Russian) Sheveakov А.А., Iakovleva R.V. Upravlenie teplobymi obiectami s raspredelenny-mi parametrami.- М: Energoatomizdat, 1986.-208 s. (in Russian) Derevich I.V., Smirnova E.G. Metod rascheta teploobmena pri protivotochnom dvi-jenii teplonositelei s peremennymi teplofizicheskimi svoistvami. Teoreticheskie os-novy himicheskoi tehnologhii, 2002, т.36, Nr.4, s.376-380. (in Russian) Abu-Hamdeh N.H. Control of a liquid-liquid heat exchanger. Heat and Mass Transfer, 38 (2001), pp.687-693.
Malinowski L. Equations for transient behavior of parallel flow multichannel heat exchangers. Heat and Mass Transfer, 39 (2003), pp.321-325. Eksergeticheskie raschety tehnicheskih sistem: Spravochnoe posobie. Brodeanskii V.M., Verhivker G.P., Karchev Ia.Ia. i dr. Kiev: Naukova Dumka, 1991.- 360s. (in Russian)
Moroziuc T.V. Novyi etap v razvitii eksergeticheskogo analiza. Holodilinaia tehnika I tehnologhia, Nr.1, 2014, s.13-17. (in Russian)
Сведения об авторах.
Шит Михаил Львович -к.т.н., в.н.с., лаборатории, «Энергетической эффективности и возобновляемых источников энергии» института энергетики АНМ. Область научных интересов: тепловые насосы, автоматическое управление технологическими процессами в
Иойшер Анатолий
Матусович - к.ф-м.н., зав. лабораторией микропроводов и нитевидных наноструктур Института "ЕИК1" (Кишинев). Область
научных интересов:
нанотехноло-гия, сенсоры, производ-ство тепла и холода с помощью вихревых
£} V
mm
энергетике, промышленности, сельском хозяйстве. E-mail: [email protected]
труб и термоэлектрических теп-ловых насосов. E-mail: amyo [email protected].
Бурчу Виталий, д.т.н., Академик РАЕН. В настоящее время является заведующим лабораторией «Энергетической эффектив-ности и возобновляемых источников энергии», Инс-титута
энергетики Акаде-мии Наук Молдовы. E-mail: [email protected]
Шит Борис
Михайлович - инженер-программист Института энергетики АНМ. Область научных интересов:
тепловые насосы, автоматическое управление технологическими процессами в энергетике, промышленности, сельс-ком хозяйстве. E-mail: [email protected].
