УДК (519.9 + 518.5); 532.54
ИССЛЕДОВАНИЕ ГИДРАВЛИЧЕСКИХ РЕЖИМОВ РАБОТЫ ЦИРКУЛЯЦИОННЫХ СИСТЕМ ТЭЦ НА КОМПЬЮТЕРНЫХ МОДЕЛЯХ
С.В. КОЛЕСНИКОВ*, И.В. КУДИНОВ**, А.В. ЕРЕМИН,* А.С. КОЛЕСНИКОВА*, А.Н. БРАНФИЛЕВА*
* Самарский государственный технический университет Самарский государственный архитектурно-строительный университет
Даны теоретические основы построения компьютерных моделей сложных многокольцевых разветвленных трубопроводных систем с использованием теории графов и двух законов Кирхгофа. Разработана компьютерная модель циркуляционной системы ТЭЦ Волжского автомобильного завода и проведены многочисленные исследования её текущего состояния, а также разработаны рекомендации по изменению режимов работы и по реконструкции. Полученные результаты могут быть использованы и для любых других циркуляционных систем подобного типа.
Ключевые слова: компьютерная модель, теория графов, законы Кирхгофа, цирксистемы ТЭЦ, эпюры давлений.
Введение
При гидравлическом расчете сложных многокольцевых трубопроводных систем эффективным оказывается применение компьютерных моделей.
В основу построения компьютерной модели положены два закона Кирхгофа, применяемые при расчетах электрических сетей. Использование этих законов, применительно к расчетам гидравлических сетей, обусловлено полной аналогией процессов протекания тока в проводниках и жидкости в трубопроводных гидравлических системах [1].
Найдем распределение расходов в сети, состоящей из одного кольца (рис. 1), имеющего три ответвления. Расходы по участкам кольца а, Ь, с, ё обозначим через ,
Яъ, Qc , Яё , а по ответвлениям - через Ql, Я2 , Qз . Требуется найти распределение
расходов по участкам а, Ь, с, ё при заданном расходе Я на входе в кольцо. Расходы
воды по ответвлениям от кольца , Я2 , Яз известны, а их сумма равна расходу Я на
входе в кольцевую сеть: Я = + Я2 + Яз .
При выполнении расчетов примем следующие допущения: 1) приток воды в узел будем считать положительным, а отток из узла - отрицательным; 2) потерю напора для расхода воды, протекающей в контуре по часовой стрелке, считаем положительной, а против часовой стрелки - отрицательной.
© С.В. Колесников, И.В. Куликов, А.В.Еремин, А. С. Колесникова, А.Н. Бранфилева Проблемы энергетики, 2013, № 7-8
Рис. 1. Схема кольцевой сети
Первый закон Кирхгофа устанавливает равенство притока и оттока среды в каждом узле сети, то есть требуется выполнение уравнения баланса расходов:
^ = 0, (1) г =1
где п - число трубопроводов, соединяющихся в узле; бг (г = 1,п) - расходы среды по всем трубопроводам, соединяющимся в данном узле.
Согласно второму закону Кирхгофа сумма напоров для любого замкнутого контура равна нулю:
= £^б6 = о, (2)
г=1 г =1
где Б г (г = 1, п) - гидравлическое сопротивление г-го участка; бг 0 = 1,п) - расходы
среды на г-м участке.
Используя уравнения (1), (2), на основе итеративного метода расчета можно найти распределение расходов по всем участкам сети при известном расходе б, заданном на входе в кольцо. На первом шаге итерации задается произвольное распределение расходов среды на каждом участке кольца, то есть задаются значения ба , бь , бс, ба . Тогда для узлов 0, 1, 2 из первого закона Кирхгофа имеем:
ба = б-ба ; ба = б\ + бь ; бь = 62 + бс . В записи уравнения Кирхгофа для узла под номером 3 нет необходимости, т. к. расход 6з = ба + бс может быть найден, если будут известны значения расходов по всем другим участкам сети.
По второму закону Кирхгофа на основе принятых расходов по всем участкам кольца найдем величину невязки напоров
8Н = £ Бб = Баб1 + Бьб1 + БсбС - Баб2 . (3)
г =1
Если величина невязки 8Н положительна, то, следовательно, перегружены участки, в которых направление движения жидкости осуществляется по часовой стрелке, и недогружены участки, где направление движения противоположное. Для приближения невязки напоров 8Н к нулю введем в расчет поправочный (увязочный) расход 86 . Этот расход должен вычитаться из величины расхода на перегруженных и добавляться - на недогруженных участках. Величину увязочного расхода 86 можно найти из соотношения (3), положив 8Н = 0. Полагая, что невязка напоров, полученная из (3), положительна, величину 86 найдем из соотношения
Ба (ба - 8б)2 + БЬ (бь - 8б)2 + Бс (бс - 8б)2 - Ба (ба + 8б)2 = 0. (4)
© Проблемы энергетики, 2013, № 7-8
Если пренебречь членами, содержащими (8<2)2, как величинами относительно малыми, то соотношение (4) относительно увязочного расхода SQ будет представлять алгебраическое линейное уравнение. Его решение
8Q = ?>н,
( п л SiQi
(5)
'I
\ i=1 у
После нахождения SQ уточняются расходы на участках кольца, и указанная последовательность расчетов вновь повторяется. Такой итеративный расчет выполняется до тех пор, пока получаемые из двух последних итераций расходы не будут отличаться на заданную (достаточно малую) величину.
При разработке компьютерной модели определяются гидравлические характеристики трубопроводов теплосети. Потери напора в трубопроводе складываются из потерь на трение (линейные) и потерь в местных сопротивлениях:
, 2 2 l w ^ w
ЛЬ = + > (6) d 2 g 2 g
где ЛЬ - потери напора, м; X - коэффициент трения; / - длина трубопровода, м; й -
внутренний диаметр, м; w - средняя скорость, м/с; сумма коэффициентов
местных сопротивлений на участке.
Вводя понятие эквивалентной длины, потери напора на местные сопротивления можно свести к линейным потерям [2 - 9]. Эквивалентная длина местных сопротивлений находится из соотношения
!э = й X Ук. (7)
С учетом (7) формула (6) примет вид
w2 ( / 1э Л kw2 , ч
ЛЬ =-| к— + к— 1 =-((+ 1э ). (8)
2g ^ й й у 2dg
Если расход жидкости через участок Q м3/с, то
w = 4Q/пd2. (9)
Подставляя (9) в (8), находим
8к(/ + 1э ) 2 ЛЬ = 2 э Q2. п 2 gd 5
Таким образом, гидравлическая характеристика участка-трубы принимает вид
ЛЬ = sQ2,
8Х(/ + /э) 2 5 где 5 =- - гидравлическое сопротивление участка, с /м .
п 2 gd 5
Гидравлическая характеристика участка-задвижки имеет вид
ЛЬ = sQ 2,
где 5 - коэффициент сопротивления задвижки, зависящий в основном от степени ее открытия. Для полностью открытой задвижки коэффициент местного сопротивления принимается равным 0,07.
Участки-насосы в расчетной схеме гидравлической сети представлены аналитическими характеристиками, связывающими между собой напор, развиваемый насосом, и подачу. Характеристика насоса в координатах Q-H с достаточной для практики точностью может быть выражена уравнением вида
Я = Я0 - 02н Б,
(10)
где Н о - напор, развиваемый насосом при закрытой на выходе задвижке (0н = 0), м;
Рассмотрим применение компьютерных моделей для анализа работы циркуляционных систем тепловых электрических станций, которые предназначены для охлаждения циркуляционной воды в градирнях с целью использования ее для конденсации пара в конденсаторах паровых турбин.
Для нормального функционирования цирксистемы важной проблемой является поддержание заданных уровней воды в стояках и чашах градирен и аванкамерах (см. рис. 2). От уровня жидкости в стояках зависит работа разбрызгивающих устройств (сопел), а в аванкамере - устойчивая работа насосов. От этих двух факторов в существенной степени зависят производительность (расход воды) цирксистемы и охлаждающие способности градирен.
При работе циркуляционных систем в ряде случаев наблюдается рассогласование уровней жидкости в системе аванкамера-стояк градирни-чаша градирни. Это приводит к малому напору жидкости на соплах, переливу воды через борта чаш градирен или к недостатку воды в них и проч.
Теоретическая зависимость уровня жидкости в стояке градирни (рис. 2 - 4) от различных параметров цирксистемы имеет вид [9]
где Н о - напор, развиваемый насосом при закрытой на выходе задвижке; (0 = 0); ро - высота уровня жидкости в стояке градирни; ргр - высота уровня жидкости в чаше градирни; рвх - высота уровня жидкости в аванкамере; рвых - высота расположения разбрызгивающих устройств (сопел градирни); к = кн + г + ^ ; кн - коэффициент гидравлического сопротивления насоса; к - сопротивление сети на участке от части градирен до аванкамеры; к1 - сопротивление сети от насоса до сопел градирни; к2 -сопротивление сопел градирни.
Рис. 2. Схема системы градирня-аванкамера-насос: 1 - насос; 2 - аванкамера; 3 - градирня; 4 - стояк градирни; 5 - чаша градирни; 6 - разбрызгивающие устройства; 7 - трубопроводы, соединяющие чашу градирни с бетонными каналами; 8 - бетонные каналы
0Н - подача насоса, м3/с; Б - гидравлическое сопротивление насоса.
(10)
Рис. 3. Упрощенная схема цирксистемы (левая подсистема ТЭЦВАЗ): I - напорный трубопровод; II - сбросной трубопровод; А - аванкамера; 1, 2, 3, 4, 5, 6 - насосы ЦНС-1; К-1, К-2, К-3 - конденсаторы; ГР-1, ГР-2, ГР-3, ГР-4 - градирни
Рис. 4. Упрощенная схема цирксистемы с аванкамерой: 1 - емкость (чаша) градирни;
2 - аванкамера; 3 - насос; 4,5 - задвижки; б - конденсатор; 7 - стояк градирни;
8 - трубопроводы с соплами
Анализируя формулу (10), приходим к заключению, что величина ро будет возрастать при увеличении &2 (сопротивление разбрызгивающих устройств) и при уменьшении к (при кн - const , r - const). В самом деле, насос создает разность потенциалов между стояком и аванкамерой, равную величине H ст . Чем больше сопротивление в сети, тем меньше Hст . В предельном случае, когда насос не р^отает Hст = 0 , наблюдается равенство уровней в чаше градирни, аванкамере и в стояке.
Из сказанного выше следует, что управлять величиной уровня в стояке можно лишь изменением сопротивления к = кн + r + ki (регулировать величину к2 не представляется возможным), в основном - это регулированием задвижками на выходе из насосов или на входе и выходе в конденсаторы, градирни, а также включением или отключением оборудования.
Ввиду того, что уровни жидкости в стояке и аванкамере взаимосвязаны, формулу (10) можно применить и для анализа факторов, влияющих на высоту уровня в аванкамере, характеризуемого величиной рвх (рис. 4). Зависимость уровней здесь обратная. Если, например, уровень жидкости в стояке уменьшается, что происходит при увеличении к1 и уменьшении к2 , то уровень в аванкамере будет возрастать. При увеличении уровня жидкости в стояке (увеличение H ст) уровень жидкости в аванкамере будет падать. Это обосновывается тем, что количество жидкости, циркулирующей в системе, неизменно, и насос может создать разность потенциалов только за счет уменьшения одного уровня и повышения другого.
Для нормального функционирования цирксистемы весьма важным является понимание взаимосвязи между уровнями в чашах градирен и аванкамере и, по возможности, управления ими. Для вывода необходимых соотношений, с целью выявления основных закономерностей и проведения качественного анализа, рассмотрим упрощенную схему цирксистемы, включающей чашу градирни, аванкамеру и насос (рис. 5). Будем считать, что насос подает жидкость в чашу градирни, которая соединена с аванкамерой каналом, имеющим гидравлическое сопротивление г . Площадь сечения чаши градирни равна 5^, аванкамеры - £2.
Рис. 5. Чаша градирни, аванкамера,насос. 1 - чаша градирни; 2 - аванкамера; 3 - насос; г - сопротивление канала между аванкамерой и чашей градирни
Составим уравнение равновесия системы аванкамера - насос:
2
Н2 + Но - кнО, = Ратм . Уравнение взаимосвязи между уровнями в чаше градирни и аванкамере имеет
(11)
вид
Н1 + гО2 + Н2 = 0. (12)
Для исключения из уравнения (12) одного из уровней (Н1 или Н2) привлечем уравнение материального баланса:
НI + Н2 52 = М ,
где М - масса жидкости в чаше градирни и аванкамере. Выразим Н1 из (13):
М - Н 2 5 2
51
М 5 2 -- Н 2 — 51 51
Подставляя (14) в (12), находим
5, - 5 -
+ О2 -
51 - 52
Подставим (15) в (11):
О2
2
2 ( г51 - кн (51 - 52 )
= Р
51 - 52
Выражая из (16) О и подставляя в (15), получаем М г51
Н2 =-
(Ратм - Н0) ■
М
51 - 5 2
(16)
(17)
5! -52 г51 -кн(5! -52) ч , _
Из анализа формулы (17) следует, что с увеличением сопротивления кн уровень жидкости в аванкамере увеличивается и, соответственно, исходя из формулы (14) - в чаше градирни уменьшается. Количественные величины изменения уровней определяются параметрами, входящими в формулы (14) и (17).
В случае, когда с аванкамерой соединены чаши нескольких градирен, получение аналитических зависимостей весьма затруднительно. В связи с этим расчеты по определению изменения уровней в чашах градирен и аванкамере были выполнены на компьютерной модели. Результаты таких исследований применительно к цирксистеме ТЭЦВАЗ приведены ниже. © Проблемы энергетики, 2013, № 7-8
0
51 - 52
В цирксистеме ТЭЦВАЗ имеется 7 градирен, из которых 6 - брызгального типа и одна (ГР-7) - пленочного (рис. 6). На рис. 6 указаны высоты стояков градирен, отметки врезки сливных трубопроводов к стоякам и отметки высот расположения сопел для различных ярусов.
Зависимость расхода воды на соплах от высоты уровня в стояке, для какой-либо из градирен 1,2,3, показана на рис. 7 (2 = /(И)). Анализируя эту зависимость, можно заключить, что до момента, пока вода не поднимется до отметки 6,4 м (отметка врезки сливных трубопроводов к стоякам), расход на соплах градирен равен нулю. После достижения этой отметки расход за достаточно короткое время становится равным порядка 13000 т/ч. Расчеты показали, что скорость движения жидкости в стояке при 2 = 13000 т/ч составляет 0,35 м/с. После того, как жидкость достигает сливных трубопроводов (16 шт. 0 420 мм), скорость поднятия жидкости в стояке снижается до 0,125 м/с. Таким образом, примерно через 4-5 секунд сливные трубопроводы будут полностью закрыты водой, и расход на соплах возрастет за это время от нуля до 13000 т/ч. Падая с высоты 6,4 м до высоты расположения сопел 4,2 м, большие массы жидкости создают условия для возникновения гидравлического удара в распределительной системе сопел и на соплах. Достаточно сказать, что при общем расходе воды в стояке 4,44 т/с на каждую из 16 сливных труб будет приходиться по 250 кг/с жидкости, падающей с высоты более 2 м.
М
1
се! ✓ ,— и
Рис. 6. Геометрические отметки водораспределительной системы градирен ТЭЦВАЗ.
Градирни № 1, 2, 3.. .6 - брызгального типа; № 7 - пленочного типа
^_
Рис. 7. Зависимость расЕодайды через сопла от высоГй уровня в стояках градиреЕй,2{3 Сплошная линия - реальная зависимость 2 = /(Н); штриховая линия - рекомендуемая зависимость после реконструкции стояка 2 = /1(Н)
Гр №4
© Проблемы энергетики, 2013, № 7-8
118
Анализируя зависимость Q = f (к), можно также заметить, что нормальная работа градирни, когда вся жидкость из стояка сливается в чашу градирни через сопла, находится в диапазоне 6,4 - 8,5 м, где 8,5 м - верхняя отметка стояка. При достижении этой отметки начинается перелив жидкости через верх стояка, и любое дальнейшее прибавление расхода не будет иметь эффекта в смысле охладительной способности градирни. Таким образом, рабочий диапазон стояка градирни составляют лишь 8,5 -6,4 = 2,1 м. При этом давление на соплах верхнего яруса, даже если жидкость заполняет стояк доверху, составляет 8,5 - 4,2 = 4,3 м. Если стояк заполнен лишь до половины рабочей зоны (7,4 м), то давление на соплах становится еще меньше (около 3 м). Давления в 3 - 4 м и меньше явно недостаточно для нормальной работы сопел по эффективности разбрызгивания жидкости (требуется давление около 6 м). К тому же расстояние от нижнего яруса сопел до поверхности воды в её чаше составляет практически для всех градирен (исключая пленочную градирню №7) лишь 3 - 4 м. Учитывая небольшую высоту фонтанчика (при напоре на соплах 4 м - около 1,6 м (см. рис. 8) и радиуса факела разбрызгивания, такой высоты падения жидкости (около 5 м) явно недостаточно для достижения высоких паспортных показателей охладительной способности градирен. Высота врезки сливных трубопроводов в стояки и высота расположения сопел градирен в данных градирнях недостаточны. Такая конструкция не позволяет достигать паспортных давлений на соплах из-за малой высоты стояка, а также небольшой высоты факела разбрызгивания. Высота падения капель тоже оказывается заниженной. В этой связи можно сформулировать следующие предложения по реконструкции градирен (см. рис. 9): 1 - увеличить высоту стояка; 2 -выполнить врезку сливных труб на уровне сопел; 3 - принять меры к поднятию уровня жидкости в стояке увеличенной высоты.
h, м R, м
Рис. 8. Зависимости: 1 - радиуса Я; 2 - высоты к факела разбрызгивания - от напора Н
Эти мероприятия приведут к увеличению давления на соплах и, как следствие, к увеличению высоты факела разбрызгивания, что в итоге увеличит охладительную способность градирен. Кроме того, отпадут проблемы, связанные с гидравлическим ударом, приводящим к срыву сопел.
Разработанная компьютерная модель была использована для расчетов следующих вариантов работы цирксистемы ТЭЦВАЗ:
1. Все насосы и турбогенераторы в работе, оборудование с паспортными характеристиками.
2. Расчет идентифицированной цирксистемы с реальными (действительными) характеристиками.
3. Расчет цирксистемы как единой целой - без разделения на две подсистемы. Характеристики системы реальные.
Применительно к первому варианту исследований в работе находились 10 насосов, 7 градирен и 11 конденсаторов турбин. Система разделена между ТГ-7
Д м
и ТГ-8. Характеристики оборудования паспортные. Результаты исследований показаны на рис. 9. И, в частности, здесь представлены данные о высотах уровней жидкости в стояках и чашах градирен, а также расходы воды через каждую градирню и суммарный расход в цирксистеме.
Предлагаемые варианты реконструкции
и
Рис. 9. Варианты реконструкции градирен (все градирни брызгальные)
Анализируя полученные результаты, можно отметить, что высота стояков явно недостаточна для работы при таком составе оборудования, и если её увеличить, то давление на соплах составит около 7 - 8 м, что превышает рекомендуемую высоту 6 м. Этот режим работы будем считать как бы эталоном, к которому следует, по возможности, приближаться путем каких-либо реконструкций и усовершенствований. Суммарный расход воды через все насосы цирксистемы для данного варианта исследований составляет 145719 т/ч.
По второму варианту исследований (идентифицированная цирксистема) в работе находились 4 насоса ЦНС-1; 3 насоса ЦНС-2; 7 градирен, 10 конденсаторов (отключен конденсатор ТГ-6), цирксистема разделена на две подсистемы между турбогенераторами ТГ-7 и ТГ-8.
Анализ полученных результатов позволяет заключить, что суммарный расход через насосы составляет 53115 т/ч - левая подсистема, 29148 т/ч - правая. Расход через градирни - 82105 т/ч (52202 т/ч - левая и 29903 т/ч фравая подсистемы). Следует отметить низкий напор на соплах градирен (3-4 м), при этом верхний ярус сопел ГР-4 не работает. Отмечается существенная неравномерность уррвня воды в чашах градирен левой подсистемы. Так уровень воды в чаше ГР-1 ндх^дится на грани перелива 15-10 см от верхней кромки борта, тогда как в градирне 4 уровень воды в чаше на 60-70 см ниже. В градирне 1 имеет место большое разлйяйе (до 50 см) в уровнях воды? ^ чаши и приемного канала, а также чаши и аванкамеры (до 1,4 м). В правой подсистеме уровень в чашах градирен достаточно высок: 50-60 см от верхней кромки. Различие уровней в чашах и аванкамере около 0,5-0,6 м.
Следует особо отметить существенно отличающиеся подпоры насосов в левой (Нн = 5,5 м) и правой (Нн = 3,8 м) подсистемах. Значительно меньший подпор в правой части создает предпосылки для неустойчивой работы насосов. Анализ состояния конденсаторов позволяет заключить, что расходы через конденсаторы ТГ-10 и ТГ-11 существенно разнятся с паспортными их значениями. Так в конденсаторе ТГ-10 расход, по сравнению с паспортным, меньше на 11400-5340=6060 т/ч, а в ТГ-11 на 13904-7973=4931 т/ч, что свидетельствует либо о большой их загрязненности, либо о прикрытости задвижек, через которые подается вода на эти конденсаторы. То же самое можно отметить и для конденсаторов ТГ-1 (7468-5123=2345 т/ч) и ТГ-5. © Проблемы энергетики, 2013, № 7-8
Гр 4
Анализ результатов расчётов позволяет заключить, что расход, по сравнению с предыдущим вариантом, практически не изменился. Это вполне объяснимо тем, что система сеть — насосы является самонастраивающейся, т.е. изменения происходят лишь внутри системы при неизменных общих показателях. Что касается изменений внутри системы, то они весьма значительны. Так уровни воды в чашах градирен левой подсистемы значительно выше кромок бортов: в градирне 1- на 70, ГР-2 - 40, ГР-3 -30, ГР-4 - 20 см. То есть во всех этих градирнях, если не увеличивать высоту бортов чаш, будет наблюдаться перелив в больших количествах.
В градирнях правой подсистемы уровень в чашах градирен и аванкамере существенно уменьшается (до 0,5—0,7 м в чашах градирен). В результате подпор на насосах ЦНС-2 уменьшается до 1,6 м, что недопустимо по условиям их устойчивой работы. Подпор на соплах практически не изменяется по сравнению с предыдущим вариантом и составляет 3 — 4 м.
В результате такого перераспределения уровней в чашах градирен расход на насосах правой подсистемы уменьшается с 33105 до 32349 т/ч, а в левой увеличивается с 52202 до 52970 т/ч по сравнению с предыдущим вариантом. Из приведенного анализа следует, что объединение левой и правой подсистем цирксистемы ТЭЦВАЗ недопустимо.
Литература
1. Соколов Е.Я. Теплофикация и тепловые сети. М. - Л.: Госэнергоиздат, 1963. 360 с.
2. Панамарев Ю.С., Кудинов В.А. и др. Расчет гидравлических и температурных режимов работы теплосети Самарской ТЭЦ с помощью компьютерной модели // Теплоэнергетика. 2005. № 5. С. 35 - 39.
3. Кудинов И.В. Использование компьютерных моделей для проектирования тепловых сетей // Вестник Самарского государственного технического университета. Сер. Технические науки. 2010. №3 (27).
4. Кудинов В.А., Коваленко А.Г., Колесников С.В., Панамарев Ю.С. Разработка компьютерной модели и исследование работы циркуляционной системы Новокуйбышевской ТЭЦ-2 // Изв. АН. Энергетика. 2001. №6. С. 118-124.
5. Колесников С.В., Дикоп В.В., Томкин С.В., Кудинов В.А. Исследование гидравлических режимов работы цирксистемы Тольяттинской ТЭЦ на компьютерной модели // Изв. вузов СНГ. Энергетика. 2002. № 6. С. 90-95.
6. Зройчиков Н.А., Кудинов В.А., Коваленко А.Г., Колесников С.В., и др. Разработка компьютерной модели и расчет оптимальных режимов работы циркуляционной системы ТЭЦ—23 ОАО «Мосэнерго» // Теплоэнергетика. 2007. №12. С. 7 - 15.
7. Кудинов И.В. Построение компьютерных моделей систем теплоснабжения больших городов // Вестник Самарского государственного технического университета. Сер. Технические науки 2011. №1 (29). С. 212 - 219.
8. Дикоп В.В., Бухтияров А.В., Коваленко А.Г., и др. Исследование гидравлических режимов работы цирксистемы ТЭЦ Волжского автомобильного завода на компьютерной модели // Известия вузов СНГ. Энергетика. 2005. №1. С. 69—76.
9. Дикоп В.В., Бухтияров А.В., Коваленко А.Г., и др. Разработка компьютерной модели и исследование гидравлических режимов работы цирксистемы ТЭЦВАЗ // Известия вузов СНГ. Энергетика. 2005. №3. С. 60 — 65.
Поступила в редакцию 28 мая 2013г
Колесников Сергей Владимирович - канд. техн. наук, докторант кафедры «Теоретические основы теплотехники и гидромеханика» Самарского государственного технологического университета (СамГТУ).
Кудинов Игорь Васильевич - канд. техн. наук, доцент кафедры «Гидравлика и теплотехника» Самарского государственного архитектурно-строительного универститета. Тел.: 8 (846) 332-42-35.
Еремин Антон Владимирович - аспирант кафедры «Теоретические основы теплотехники и гидромеханика» Самарского государственного технологического университета (СамГТУ). Тел. 8 (917) 111 07 97. E-mail: a.v.eremin@list.ru.
Колесникова Анна Сергеевна - аспирант кафедры «Теоретические основы теплотехники и гидромеханика» Самарского государственного технологического университета (СамГТУ).
Бранфилева Анастасия Николаевна - аспирант кафедры «Теоретические основы теплотехники и гидромеханика» Самарского государственного технологического университета (СамГТУ).