Оценка надежности экранных труб топки котла к блоку 660 МВт ТЭС «Бар» в стационарных режимах

УДК 621.184.54

ОЦЕНКА НАДЕЖНОСТИ ЭКРАННЫХ ТРУБ ТОПКИ КОТЛА К БЛОКУ 660 МВт ТЭС «БАР» В СТАЦИОНАРНЫХ РЕЖИМАХ

© 2013 г. А.А. Белов , А.Н. Безгрешное , А.Н. Озеров , А.Б. Баранников , \В.Г. Лебсак

*Южно-Российский государственный *South-Russian State

технический университет

ОАО «Энергомаш»

Technical University

**OJSC «Energomash»

Проведен расчетный анализ надежности гидродинамики котла к энергоблоку 660 МВт ТЭС «Бар». Показано, что выбранная гидравлическая схема экранных поверхностей нагрева удовлетворяет критериям надежности как в прямоточных режимах при сверх- и докритическом давлениях, так и в сепараторном режиме.

Ключевые слова: гидравлическая схема; гидродинамика котла; анализ надежности; экранные поверхности нагрева; сепараторный режим; прямоточный котёл; сверхкритическое давление.

Conducted a computer analysis of the reliability of hydrodynamics of the boiler of the power unit of 660 mW power plant «Bar». It is shown that the hydraulic circuit-screen surfaces of heating meets the criteria of reliability as an in-through modes for super - and subcritical pressure, and in сепараторном mode.

Keywords: hydraulic scheme; hydrodynamics boiler; reliability analysis; on-screen surface heating; separator mode; direct-flow boiler; supercritical pressure.

При проектировании прямоточных котлов центральное место занимает проблема надежной работы экранных поверхностей топки. Особенно остро эта проблема стоит для прямоточных котлов, предназначенных для работы в широком диапазоне изменения параметров пара, как это предусмотрено в техническом задании пылеугольного парового котла для блока 660 МВт ТЭС «Бар», который должен работать как при сверхкритическом давлении (СКД) в диапазоне нагрузок 90 - 100 % _Оном, так и при докритическом давлении (ДКД) в скользящих (прямоточных) режимах и в сепараторном режиме при пусках котла. Оценка надежности экранов при этом должна производиться как в соответствии с требованиями руководящих материалов по проектированию паровых котлов СКД [1], так и по критериям надежности их работы при докритических давлениях [2, 3].

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

Тепловой расчет котла проводится с целью определения тепловосприятий поверхностей нагрева [Qj]. Для получения более надежных результатов этот расчет был выполнен по двум различным программам «Boiler-designer» и «Тракт», показавшим хорошую сходимость.

[Qj]

С целью более точного выполнения гидравлических расчетов и расчетов показателей надежности обогреваемая часть панелей (блоки) топочных экранов по высоте была разбита на 11 участков (рис. 2) в соответствии с геометрией топки и разделением экранов на нижнюю (НРЧ) и верхнюю (ВРЧ) радиационные части. Участок 3 (между осями ярусов горелок) соответствует зоне активного горения топлива. Для адекватного моделирования температурного режима в районе конца НРЧ и начала ВРЧ (т.е. в зоне сочетания максимальных тепловых потоков и достаточно высоких температур рабочей среды), а также в районе аэродинамического выступа выделены участки 5, 6 и 8 высотой 1 м.

Ввиду двойной симметрии, тепловосприятия рассчитывались для одной четвертой части топочных экранов. Найти на каждом участке панели (а таких участков - 99) тепловой поток с помощью программы «Тракт» или «Boiler-designer» практически нереально. Для этой цели была использована программа распределения нагрузок «PZNAG», разработанная на кафедре «Парогене-раторостроение» ЮРГТУ (НПИ) и успешно, в течение многих лет, применяемая на ОАО ТКЗ «Красный котельщик». С помощью этой программы были найдены тепловосприятия каждого участка всех панелей с учетом нормативных значений коэффициентов неравномерности по высоте ув [3] и ширине [2] топочных экранов, общее тепловосприятие которых принималось из теплового расчета по программе «Тракт».

Тепловой расчет котла

[%1

Gl

[Aj.

Гидравлический расчет

[Gjl

Расчет показателей надежности

Рис. 1. Этапы оценки надежности

На рис. 2 отмечены значения коэффициента ув в соответствии с зонами, рекомендуемыми в [3, табл. 8.3], а для участков экранов эти коэффициенты определялись линейной интерполяцией между величинами на границах зон.

Верхний

62.75

х 3,7

>в = 0.9

- 1,0.

42.0

©-

-21,45

>>„ = 1.3

6

^—

1,0_ 1,0_

6,4

ярус горелок 4 у = 4

о 3"

[-4 1—1

¡1 0-20,1

Нижний

ярус горелок

IX

4 = 1-2 Q___3,3

0-

4,8

\

Ув = 0.6

59.05

ВРЧ

36.6

28.2

НРЧ

4.8

0.0

Гидравлический расчет выполнялся с помощью программы «Гидравлика», разработанной на кафедре «Парогенераторостроение» ЮРГТУ (НПИ) и зарегистрированной в отраслевом фонде алгоритмов и программ [4]. Расчетная схема в этой программе представляется в виде ориентированного графа произвольной степени сложности, для которого решаются основные уравнения сохранения.

Математическая модель гидравлической схемы котла представляет собой систему уравнений, состоящей из уравнений сохранения массы для смешивающего узла (1), сохранения импульса для замкнутого контура (2), энергии для узлов (3) и дуг графа (4):

[a0iJ][G}] +[ G,r] = [0];

[bj][Apj] = [0];

[-aj ][ hjKGj+] - j h/Gj] - [WG4 = [hG];1

[hjK] - [hjH] = [Qj / Gj],

гзф 2ф 1ф

Рис. 2. Схема разбиения экранов топки на участки: О _ ©-номера участков

Оценка надежности, а соответственно, тепловые и гидравлические расчеты проводились для пяти нагрузок: 8, 15, 32, 50 и 100 % от номинальной производительности Dном. При этом котел работает при различных параметрах (ДКД и СКД) в трех качественно различных ситуациях (режимах):

1) сепараторный режим при ДКД (нагрузка 8, 15 и 32 %);

2) прямоточный режим при ДКД (нагрузка 50 %);

3) прямоточный режим при СКД (нагрузка 100 %).

Основной целью гидравлических расчетов, со-

гласно рис. 1, для всех трех режимов является определение расходов [Gj•], энтальпий и давлений [р], необходимых для расчета показателей надежности поверхностей нагрева котла. Кроме этого, для прямоточных режимов важным показателем является общий перепад давления в пароводяном тракте, который регламентируется по условиям контракта. А для сепараторного режима дополнительно проверялись соответствия характеристик насоса (напор и кавитацион-ный запас) данным условиям эксплуатации.

где [а0/] - сокращенная матрица инцидентности; [GгГ] - вектор-столбцы поперечных величин графа (массовый расход среды в дугах) и граничных расходов; [Ьк/] - матрица контуров; [к] [квх/г] - вектор-столбцы среднерасходных энтальпий в узлах и граничных энтальпий на входе в систему; [к/], [к/] -вектор-столбцы энтальпий на входе и выходе из дуг; [а/+], [а/], [G/], [Gj~] - положительная и отрицательная части матрицы [а/ = [а/+] + [а/] и вектор-столбца расходов [Gj] = [Gj+] + [Gj~]; [G/S] = =[-а// ] [G/] -[а/] [Gj~]-[G/Г~] - вектор-столбец суммарных расходов, входящих в узел; [Огт~] - отрицательная часть вектор-столбца граничных расходов [^-г] = [G/г+] + [^-г_]; Qj - тепловосприятие /-го компонента, кДж/с.

Если Gj > 0 или Gj = 0 и Qj = 0, то вектор энтальпий на входах в дуги определяется по формуле

[/ = [а,/+к/+]т[к/],

где [к/] - вектор-столбец кодов, учитывающих знак массового расхода в дуге.

Вектор [к/] состоит из нулей и единиц (1 - если расход в дуге больше нуля). Энтальпия на выходе из дуги [к/] определяется из уравнения (4).

Если Gj < 0, то вектор энтальпий на входах в дуги определяется по формуле

[/ = [~а~к-]Ш

где [к/] - вектор-столбец кодов, учитывающих знак массового расхода в дуге.

Вектор [к[] состоит из нулей и единиц (1 - если расход в дуге меньше нуля). Энтальпия на выходе из дуг [к/] определяется из уравнения (4).

Если энтальпии на выходе из смешивающего узла не одинаковы, например в полнопроходном сепарато-

1 Запись вида [а/Ь/] означает вектор, /-ми компонентами которого служат элементы а/Ь/.

ре, то к вышеприведенным уравнениям добавляются соотношения, отражающие закон сохранения энергии для потоков, исходящих из такого узла [5]. Возможность моделирования неравномерности энтальпий на выходе из узла является особенностью данной модели и соответствующей программы [4]. Кроме этого, как видно из приведенной модели, в программе предусмотрено наличие отрицательных расходов. Эта особенность оказалась очень полезной при моделировании узла рециркуляции (сепараторный режим), где при определенных условиях в некоторых компонентах возникал отрицательный расход.

При выполнении гидравлических расчетов использовались расчетные схемы с различной степенью детализации. Такой подход позволил уменьшить трудозатраты, повысить наглядность схем, причем, как показали специальные тесты, точность гидравлического расчета не пострадала. В общей расчетной гидравлической схеме пароводяного тракта до полнопроходного сепаратора, включая растопочный узел (рис. 3), четыре панели 1-го хода НРЧ объединены в один элемент, пять панелей 2-го хода НРЧ объедине-

ны также в один элемент, а девять панелей 3-го хода ВРЧ представлены в виде двух элементов

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

С применением межпанельной расчетной схемы (рис. 4) определялись действительные расходы в отдельных панелях и их разверенных трубах.

Для общей расчетной схемы использовались граничные условия, представленные в табл. 1.

Рис. 3. Общая расчетная гидравлическая схема 47

Таблица 1

Граничные условия для общей расчетной гидравлической схемы

Величина Вариант нагрузки, %

8 15 32 50 100

Давление на входе P, бар 35,4 75,6 132,3 145,3 295,5

Энтальпия на входе 1, кДж/кг 438,5 441,5 445,7 1057,0 1278,0

Расход среды на входе в котел G, кг/с

Гидравлическое сопротивление байпаса 1-го хода НРЧ было установлено таким образом, чтобы расход через него при нагрузке 100 % составил 20 % от общего расхода. При расчетах на пониженные нагрузки геометрические характеристики байпаса не изменялись. В сепараторных режимах использовался насос фирмы Хайвард Тайлер с известными характеристиками. Рас-

ход в экранных поверхностях в этих режимах принимался одинаковым при всех нагрузках и равным 700 т/ч (194,4 кг/с). Это постоянство расхода обеспечивалось изменением местного сопротивления регулирующего клапана на выходе из циркуляционного насоса. Результаты расчета массовых скоростей р^ в основных элементах общей расчетной схемы приведены в табл. 2.

Как видно из табл. 2, потоки рабочей среды по параллельным ветвям распределяются достаточно равномерно.

Общий перепад давления от входа до собирающего коллектора за растопочным узлом в зависимости от нагрузки составил: для нагрузки 8 % - 9,41 бар (9,60 ат), 15 % - 6,56 бар (6,69 ат), 32 % - 4,65 бар (4,74 ат), 50 % - 15,25 бар (15,55 ат), 100 % - 29, 84 бар (30,44 ат).

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

Снижение общего перепада давления с ростом нагрузки в сепараторных режимах обусловлено следующими причинами:

1) при увеличении нагрузки уменьшается расход рециркулирующей среды, имеющей температуру насыщения при давлении в сепараторе, и увеличивается прямоточный расход с температурой среды после экономайзера. В результате на входе в НРЧ энтальпия среды уменьшается, что приводит к росту экономай-

зерных и уменьшению паросодержащих участков труб панелей;

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

В итоге, несмотря на увеличение массового паро-содержания с ростом нагрузки, общий перепад давления от входа до собирающего коллектора за растопочным узлом снижается. Соответственно этому, гидравлическое сопротивление регулирующего клапана, который стоит на напорном трубопроводе циркуляционного насоса, увеличивается от 9,65 бар при нагрузке 8 % до 15,31 бар при нагрузке 32 %. При этом, в соответствии с характеристиками насоса фирмы Хайвард Тайлер, минимальный запас по напору равен 9,65 бар.

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

Таблица 2

Основные результаты расчетов расходов по общей расчетной гидравлической схеме

Наименование элемента № дуги (НПИ) Массовые скорости рм>, кг/(м2-с), при нагрузках, %

8 15 32 50 100

Экономайзер (Эк) П3 304,5 304,5 304,5 445,2 890,4

Нижняя радиационная часть (НРЧ 1) -ход 1 П14 682,0 662,5 626,8 879,9 2019,0

Нижняя радиационная часть (НРЧ 2) -ход 2 П17 687,7 687,7 687,7 1006 2011,0

Боковые экраны верхней радиационной части (ВРЧ Б) - ход 3 П19 453,6 452,6 457,1 685,9 1370,0

Фронтовые экраны верхней радиационной части (ВРЧ Ф) - ход 3 П44 479,5 480,6 473,7 665,7 1335,0

Фестон 1 (Ф1) П21 446,3 446,2 446,2 652,4 1305,0

Потолочный экран (Пот) П22 511,6 511,5 511,5 747,9 1496,0

Фестон 2 (Ф2) П25 427,1 429,8 433,3 619,3 1239,0

Горизонтальный газоход 1 (ГорГх 1) П46 341,0 339,7 340,3 503,4 1007,0

Горизонтальный газоход 2 (ГорГх2) П109 361,7 359,1 354,8 532,2 1064,0

Фестон 3 (Ф3) П36 281,2 288,2 299,5 396,8 791,2

Внутренний боковой экран конвективной шахты (Вн БЭ КШ) П37 321,0 312,5 283,9 480,1 961,1

Фронтовой и задний экраны конвективной шахты (ФиЗЭ КШ) П33 338,1 341,3 361,4 493,5 987,5

Наружный боковой экран конвективной шахты (Нар БЭ КШ) П30 350,3 354,4 370,5 507,1 1014,0

Граничные условия для межпанельной расчетной гидравлической схемы

Величина Вариант нагрузки, %

8 15 32 50 100

Давление на входе Р, бар 39,63 79,70 136,80 148,40 296,50

Энтальпия на входе i, кДж/кг 859,30 944,20 671,50 1316,00 1493,0

Расход среды на входе в НРЧ G, кг/с

Таблица 3

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

Таблица 4

Минимальные значения коэффициентов межпанельных гидравлических разверок ргт|п в экранах топки

Экран Панель ргт1П при нагрузках, %

8 15 32 50 100

НРЧ ход 1 1б 0,901 0,874 0,828 0,866 0,945

НРЧ ход 2 2б 0,888 0,878 0,848 0,939 0,957

ВРЧ 5б 0,900 0,905 0,933 0,932 0,943

Как видно из табл. 4, распределение расходов по панелям экранных поверхностей нагрева достаточно равномерное.

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

В результате расчета надежности при СКД (нагрузка 100 %) получено, что температура наружной стенки лобовой части разверенных труб лежит в пределах от 402 °С (НРЧ 1) до 549 °С (ВРЧ Б). Эти тем-

пературы значительно ниже предельной по условию окалинообразования [2].

Результаты оценки температурного режима при ДКД показаны в табл. 5.

В табл. 5 приняты следующие обозначения: /нар, 4н - температуры наружной и внутренней поверхности трубы в ее лобовой части; р^, р^пул - массовые скорости в трубе и минимальная по условию возникновения пульсаций; х, хпр - действительное и предельное массовые паросодержания.

Как видно из табл. 5, в прямоточном режиме (нагрузка 50 %) и в сепараторном режиме при нагрузке 32 % в некоторых экранах х > хпр. Эта ситуация типична для прямоточных режимов, а при нагрузке 32 % режим близок к прямоточному. Однако в этом случае увеличение температуры металла из-за ухудшения теплоотдачи незначительно и лежит в пределах 5 ... 15 °С.

Расчет коэффициентов запаса по застою и опрокидыванию потока показал, что все экраны по этим показателям удовлетворяют требованиям норм [2]. Минимальный коэффициент запаса по застою равен 1,3 и наблюдается в панели 1б ВРЧ Б, а наименьший запас по опрокидыванию равен 2,1 в панели 1ф ВРЧ Ф.

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

Таблица 5

Основные результаты расчета показателей надежности при ДКД для разверенных труб

Нагрузка, % Экран (панель) t °С "-нар? ^ ^вн, С pwTp, кг/(м2-с) Р^пуЛ, Кг/(м4) x хпр

8 НРЧ 1 (3ф) 254 248 686 355 0,009 0,847

НРЧ 2 (6б) 256 249 793 426 0,082 0,795

ВРЧ Ф (3ф) 254 247 447 * 0,069 1,0

ВРЧ Б (5б) 253 247 416 * 0,069 1,0

15 НРЧ 1 (3ф) 309 299 728 320 -0,01 0,824

НРЧ 2 (6б) 309 297 806 550 0,155 0,789

ВРЧ Ф (3ф) 308 296 457 * 0,125 1,0

ВРЧ Б (5б) 308 296 421 * 0,126 1,0

32 НРЧ 1 (3ф) 338 319 741 ** -0,170 0,60

НРЧ 2 (6б) 366 342 847 469 0,224 0,550

ВРЧ Ф (3ф) 404 390 477 * 0,827 0,724

ВРЧ Б (5б) 405 391 462 * 0,841 0,736

50 НРЧ 1 (3ф) 378 351 962 628 0,384 0,467

НРЧ 2 (6б) 427 401 950 451 0,865 0,479

ВРЧ Ф (3ф) 496 480 609 385 1,493 -

ВРЧ Б (5б) 493 476 626 372 1,474 -

Примечание:*, ** - в панели отсутствует экономайзерный или испарительный участки соответственно.

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

В работе исследовались гидравлические характеристики ветвей, образующих экранную систему первого и второго хода НРЧ и экранов ВРЧ в режимах 100, 50, 32, 15, 8 % от номинального расхода пара Dном и гидравлические характеристики средних и раз-веренных труб панелей.

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

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

Выводы

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

2. Гидравлические характеристики экранов топки однозначны и обеспечивают устойчивую гидродинамику в экранах котла.

3. Общий перепад давления от входа до собирающего коллектора за растопочным узлом для 100 % нагрузки составил 30,4 ат.

Поступила в редакцию

4. Циркуляционный насос фирмы Хайвард Тайлер имеет большой запас по напору.

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

Литература

1. Проектирование котлов сверхкритического давления на скользящем давлении (дополнение к Нормативному методу гидравлического расчета котельных агрегатов). Руководящие указания. Вып. 56. Л., 1988. 17 с.

2. Гидравлический расчёт котельных агрегатов : нормативный метод / О.М. Балдина [и др.]. М., 1978. 256 с.

3. Тепловой расчет котлов (Нормативный метод): 3-е изд., перераб. и доп. СПб., 1998. 256 с.

4. Баранников А.Б., Белов А.А., Федоров В.С. Программа поверочного гидравлического расчета «Гидравлика» : Свидетельство об отраслевой регистрации разработки № 5251 / Федер. агентство по образованию; Гос. коорд. центр информ. технол.; Отрасл. фонд алгоритмов и программ. Зарег. 04.10.2005.

5. Белов А.А., Баранников А.Б. Математическая модель гидравлической схемы котельного агрегата с учетом неравенства входной энтальпии в потоках, исходящих из одного узла // Кибернетика электрических систем : материалы XXIII сессии семинара «Энергоснабжение промышленных предприятий», Новочеркасск, 25 - 28 сент. 2001 г. / Юж.-Рос. гос. техн. ун-т.: ред. журн. «Изв. вузов. Электромеханика», Новочеркасск 2002. С. 86 - 88.

6 ноября 2012 г.

Белов Александр Алексеевич - канд. техн. наук, заведующий кафедрой «Парогенераторостроение», Южно-Российский государственный технический университет (Новочеркасский политехнический институт). Тел. (863-52) 55-644. E-mail: warme@npi-tu.ru

Безгрешнов Александр Николаевич - канд. техн. наук, профессор, кафедра «Парогенераторостроение», Южно-Российский государственный технический университет (Новочеркасский политехнический институт). Тел. (863-52) 55-644. E-mail: warme@npi-tu.ru

Озеров Александр Николаевич - канд. техн. наук, профессор, кафедра «Парогенераторостроение», Южно-Российский государственный технический университет (Новочеркасский политехнический институт). Тел. (863-52) 55-644. E-mail: warme@npi-tu.ru

Баранников Алексей Борисович - канд. техн. наук, доцент «Парогенераторостроение», Южно-Российский государственный технический университет (Новочеркасский политехнический институт). Тел. (863-52) 55-644. E-mail:warme@npi-tu.ru_

Лебсак Владимир Геннадиевич! - главный конструктор проекта ОАО «Энергомаш».

BelovAleksandrAleksandrovich - Candidate of Technical Sciences, head of department «Steam-Generating Buildings», South-Russia State Technical University (Novocherkassk Polytechnic Institute). Ph. (863-52) 55-644. E-mail: warme@npi-tu.ru

Bezgreshnov Alexander Nikolaevich - Candidate of Technical Sciences, professor, department «Steam-Generating Buildings», South-Russia State Technical University (Novocherkassk Polytechnic Institute). Ph. (863-52) 55-644. E-mail: warme@npi-tu.ru

Ozerov Alexander Nikolaevich - Candidate of Technical Sciences, professor, department «Steam-Generating Buildings», South-Russia State Technical University (Novocherkassk Polytechnic Institute). Ph. (863-52) 55-644. E-mail: warme@npi-tu.ru

Barannikov Alexey Borisovich - Candidate of Technical Sciences, assistant professor, department «Steam-Generating Buildings», South-Russia State Technical University (Novocherkassk Polytechnic Institute). Ph. (863-52) 55-644. E-mail: warme@npi-tu.ru

Lebsak Vladimir Leonidovich| - chief designer of the project of OJSC «Energomash».