Научная статья на тему 'Уравнения режимных задач шестифазных турбогенераторов предельной мощности'

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

CC BY
151
25
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ТУРБОГЕНЕРАТОРЫ ПРЕДЕЛЬНОЙ МОЩНОСТИ / ШЕСТИФАЗНЫЕ ТУРБОГЕНЕРАТОРЫ / РЕЖИМНЫЕ ЗАДАЧИ / УРАВНЕНИЯ СОСТОЯНИЯ / ЛИНЕЙНЫЕ ПРЕОБРАЗОВАНИЯ. / TURBINE GENERATORS OF MAXIMUM CAPACITY / SIX-PHASE TURBINE GENERATORS / OPERATION CONDITIONS PROBLEMS / OPERATION CONDIONS EQUATIONS / LINEAR TRANSFORMATION.

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

В целях увеличения единичной мощности сверхмощных турбогенераторов могут применяться многофазные обмотки статора, в том числе шестифазные, которые состоят из двух сдвинутых в пространстве на 30 электрических градусов трехфазных обмоток. Мощные шестифазные турбогенераторы находят широкое применение на строящихся в настоящее время энергоблоках атомных электростанций. Такие турбогенераторы за счёт наличия индуктивной связи между трехфазными системами обмотки статора обладают рядом специфических особенностей по сравнению с трехфазными. Ряд вопросов, касающихся шестифазных генераторов предельной мощности, нашёл освещение в литературе. Вместе с тем режимные задачи и связанные с ними аспекты изучены недостаточно. В статье, используя обычные допущения, применяемые при анализе трехфазных машин, получены дифференциальные уравнения шестифазного турбогенератора в фазных координатах. Определены соотношения между индуктивностями самои взаимоиндукции контуров статора и ротора. Индуктивности рассеяния контуров статора представлены в виде суммы индуктивностей рассеяния по путям самоиндукции и взаимоиндукции между трехфазными системами обмотки статора. Произведено преобразование уравнений к вращающимся синхронно с ротором координатам, что обеспечивает большую наглядность и меньшую ресурсоемкость вычислений. Выявлены отличия от известных выражений для трехфазных машин. Определено уравнение движения ротора шестифазного турбогенератора. Произведен выбор относительных единиц, позволяющий упростить полученные выражения и проводить сопоставление параметров с параметрами обычных трехфазных машин. В относительных единицах получены уравнения состояния шестифазного турбогенератора аналогичные уравнениям Парка–Горева трехфазных машин. Это позволяет по аналогии с известными методами трехфазных турбогенераторов проводить исследования режимов работы и переходных процессов шестифазного турбогенератора.

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

EQUATIONS FOR OPERATION CONDITIONS PROBLEMS OF SIX PHASE TURBINE GENERATORS OF MAXIMUM CAPACITY

Polyphase stator winding can be used for increasing the unit capacity of high-power turbine generators. This includes six-phase winding consisting of two separate three-phase windings displaced by 30 electrical degrees. The modern high-power six-phase turbine generators are widely used in the construction of nuclear power plant units. Six-phase turbine generators have some specifics compared with three-phase generators due to mutual inductive coupling between three-phase stator winding systems. Some issues regarding six-phase generators have been described in literature; however, the operation conditions and related aspects have not been studied enough. This article has obtained differential equations defining the operation conditions of a six-phase turbine generator in phase coordinates using the common assumptions from the theory of three-phase machines. The relationship between self and mutual induction of rotor and stator circuits is defined. Leakage reactance of stator circuits is presented as a sum of leakage reactance for the self and mutual induction paths between three-phase stator winding systems. The equations were converted to rotating coordinates synchronous with rotor to provide higher clarity and reduce the amount of resources required for the calculations. We have explored the differences from the well-known equations for three-phase machines. The equation for the rotor movement of the six-phase generator has been determined. We have chosen the relative units allowing to simplify the obtained expressions and compare the parameters with those of ordinary three-phase generators. By this equations of six-phase generator operation are obtained. These relative units enables us to find the equations similar to the Park-Gorev-form equations for three-phase turbine generators. Solving these equations by traditional methods known for three-phase machines allow studying the operation modes and transient processes of six-phase turbine generators.

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

ЭНЕРГЕТИКА

DOI 10.5862/JEST.243.1 УДК 536.7

В.А. Барилович, Ю.А. Смирнов

ПАРОСТРУЙНЫЙ ИНЖЕКТОР: РАСЧЕТ И ПРИМЕНЕНИЕ В ОТОПИТЕЛЬНЫХ СИСТЕМАХ

V.A. Barilovich, Yu.A. Smirnov

STEAM-EJECTING INJECTOR: СALCULATЮN AND APPLICATION IN HEATING SYSTEMS

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

ИНЖЕКТОР; СОПЛО; КАМЕРА СМЕШЕНИЯ; ДИФФУЗОР; ПОДОГРЕВАТЕЛЬ СЕТЕВОЙ ВОДЫ; НАСОС; ЭКОНОМИЯ ТОПЛИВА.

On the basis of the system of one-dimensional stationary differential equations of conservation laws а calculation of the injector working at water vapor is executed. The main parts of the injector are considered: an accelerating Laval nozzle, a mixing camera and a diffuser with a cylindrical entrance. Calculation allows determining the characteristic sizes of a design of the injector and change of the key physical parameters of water and steam by its length. The calculation results are compared to author's own experimental data. The scheme with an injector for heating of network water is suggested and compared to existing ones. The injector usage instead of ordinary water heater with a pump leads to significant saving of fuel and electric energy.

INJECTOR; NOZZLE; MIXTURE CAMERA; DIFFUSER; HEATER OF NETWORK WATER; PUMP; ECONOMY OF FUEL.

Введение

Пароструйные инжекторы находили и находят применение в самых различных областях техники. Изобретенный в 1858 году Г. Джиффар-дом инжектор [1] использовался для подачи воды в котел паровоза, но был вытеснен появивши-

мися позднее более эффективными устройствами — насосами. В период бурного развития космонавтики возникла необходимость в компактных энергетических установках — жид-кометаллических магнитогидродинамических (МГД) преобразователях для космических объектов. С этого времени вновь начались исследо-

вания как МГД установок в целом, так и инжекторов, являющихся одним из основных элементов этих устройств. В этой области работали коллективы Института высоких температур АН СССР, Энергетического научно-исследовательского института имени Г.М. Кржижановского, Московского энергетического института, Ленинградского политехнического института, Казанского авиационного института. Основные результаты расчетно-опытных исследований изложены в [2—6], [8—10]. Значительный вклад в исследование струйных аппаратов и систем на их основе внесли и зарубежные ученые [1, 7].

В настоящее время инжектор вновь привлек внимание исследователей. Этот интерес обусловлен тем, что профессор В.В. Фисенко нашел новую область для его применения: он предложил с целью экономии топлива заменить существующие системы отопления зданий, включающие теплообменные аппараты и насосы, инжекторами [11]. Актуальность развития энергосберегающих технологий в современных условиях не вызывает сомнений.

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

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

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

Разгонное сопло

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

(di dp

dw v dp 2cfw — + —— + —f— = 0; dz w dz D

dv v dw v dF = o_

dz w dz F dz

dT , „ л T dp

--(v -v )--= 0;

dz r dz

dv dw

(9)

dp + ( di dz

dv

P

— + w — = 0. dz dz

Разрешая эту систему уравнений относительно dp/dz, получим

dp 1

dz

Ъ, (

D

д_ dp

d__ dv

- v -

Jy

д i

w2 l3vyP

x

л

+ w P J

F i dv

dF_ P dz

Для двухфазной области

di 1 = T—= r • dv P dT v'- v''

f -v • л d i

dp Тогда

=c

Jv

P dT

v— + v = dp

cvT i " r\ . —— (v - v ) + v . r

dp dz

w

>(M2 -1)

2c;

(

D

1 + W- (v''-v') vr

_L L

w

>(m 2 -1)

2c,

D

1 + -

w

_L L

Окончательно получаем

dp dz

w

>(M2 -1)

2c,

(

dw dz dv dz

D

v w

1 + -

w

V

лdp dz

2cfw.

1 L

D

\

dw + v dF; dz F dz'

dT = v'-v') TdP, dz r dz

где M = w/a;

а =

rv

1

F =

4

v = (1 -x)v' + xv'; cv = (1 -x)cv + xc'v; c'v и c"v -удельные изохорные теплоемкости соответ-

ственно насыщенной жидкости и пара со стороны двухфазной области.

Для сопла, состоящего, например, из двух усеченных конусов с цилиндрической проставкой:

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

D = Dвx-Аг)г/1ск, если 0<г</ск;

D = ^ если /ск <г<(1ск +1сц); (7)

D = DГ + (Дср - Dг) (г - (/ск + /Сц))/!Сд, если (/ск + 1сц + 1сд ) — г — (/ск + 1сц ) •

Таким образом, на основе прямой задачи получена замкнутая система уравнений, позволяющая выполнить расчет сопла, работающего на сухом насыщенном паре. Численный эксперимент показывает: чтобы процесс пошел по нижней ветке кривой (dp/dz<0), достаточно задать расход выше критического не более чем на 1 %.

Камера смешения

В работе [2] выполнено опытное исследование пароводяного инжектора, геометрические размеры которого приводятся. Достигнутая степень повышения давления составила всего около 1,4 , что обусловлено значительной потерей кинетической энергии двухфазного потока из-за несовершенной проточной части, в частности из-за большой длины (100—200 мм) камеры смешения (КС). С увеличением коэффициента ин-жекции и = 0^/0п у инжектора с длинной камерой создаются условия для отрыва потока (в этом случае протяженность конденсирующейся струи меньше длины КС) и в еще большей степени возрастают гидравлические потери. При малых значениях и невозможна полная конденсация пара, что в конечном итоге не позволило авторам [2] получить высокое давление за инжектором.

Наши опыты с инжекторами, работающими на перегретой воде, близкой к состоянию насыщения (недогрев составлял 3—5°) при рвх = 0,6—0,8 МПа, показали [8], что из трех конических камер длиной 0,15, 0,1 и 0,05 м наибольшую степень повышения давления я* = = 1,8—2 имел инжектор с самой короткой камерой, при этом максимум я* с уменьшением I /d смещается в сторону больших значений

кс' мин г ^

и. Максимальное значение я* было получено при

I /d = 12,5 и диаметре среза сопла Лаваля,

кс мин

равном 15 мм.

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

'т/тг =

6 (1 - х

кр

'V гк г

пА-кРг к ^ ™

и холодная жидкость.

Предполагалось, что в начальном сечении КС холодные и горячие капли равномерно распределены в сухом насыщеннном паре, но в целом поток термически и механически неравновесен. При его движении в камере и наличии обменных процессов он будет переходить в равновесное состояние: значения температуры и скорости фаз будут близки между собой. Переход в это состояние или близкое к нему и определит длину камеры смешения. Поскольку ввод холодной жидкости в КС осуществляется, как правило, при незначительной скорости, то считается, что в результате силового взаимодействия между несущим потоком и холодными каплями происходит дробление последних, а диаметр мелких холодных капель Dх к определяется по формуле [12]

А,

ч0,25

А,

1,9

Же,

0,25

+ 0,315

0х к

_Р0л Р0х:

\1,5

сЖД251п

А0х к Ахк

Ро'п (^п -№хк)0 ^Охк где Ше0 =--критерий Вебера,

■'Охк

определяемый условиями на входе в КС;

К - ^хк )0 А0хкР0П

"х0

= /^); Reo =

М-оп

При

этом делалось допущение, что разрушение капель происходит на участке, протяженность которого много меньше длины КС.

Ниже приводится система уравнений, позволяющая выполнить расчет КС на основе прямой задачи:

dp dz

w2 ° гк Г ^п

dp/dрП) С п 1 к

х [ Г1 ™пр п V ^г к Ргк

V " Л2 Рп

С хк к А3

С п V ™ х к ) ^Рх к )

( П "

1— ^пРп о'* л

1 ^х кРХк у

2 А -1

+п!к,

1 +

М — М W Р п гк пгп

Мп М Р

п г к г к

1 + »пРп (»п — »х к ) — »пРп

»х кРх к — »пРп

» х к Рх

2

Агк«икД/(гк—п) ■

Аос «ухк Д/(п—х

^кс + У пDкс 1 .

»п

dw.r

1 dp 6Д

-+

*

х к

dz dz пD3 р ^

К— ^хк) йг;„—хк.

dz '

гк

1 йр + 6^

*

гк

Лх ^ГКРГК & Л^кРгк^Гк

г

f ± + т пБ | 6Дх к^х к п . + V ПБкс + „3 ^ пБх к Рх к ^

х к х к х к

Л

ГК ГК +(^п — )- ~гк-п

пБ3 Р ^ У dz

гкггк'гк

йх к = йг

= [(С — 4к)4/(п— хк) — (С (ТХк)— 4к)4/(хК—п) + %—хк) (Тп — ТХк)]х

6

йг

й/^ к ^^хк

^хк

йг

™хкРхк-°хк

4 с

йг

рхк

V //

где D* = сх^Ок2Рп (»п — »к )|»п — »к| /8 -сопротивления единичной капли;

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

Сх = / (Reк); Reк =

^ |Кп — Кк |Рк

^п

А/п

1 ( 2х А р8 рхк(Тхк)

х=

35

0,56 '

- = п!)т к п¥ т к А/Ш-ПГКС;

ск

а = (2 + 0,6ReК'5 Ргпп/БК.

Система из шести дифференциальных уравнений позволяет произвести расчет КС на основе прямой задачи, то есть определить изменение давления, температуру пара и холодных капель, скорости пара, холодных и горячих капель по длине камеры.

Температура горячих капель определяется из уравнения энергии для адиабатического потока

(

Ог

'гк +-

w

(

+

+-

wп

(

+Ох

1хк +-

w

= I.

(14)

сила

см0'

(15)

У*

где I см0 — полная энтальпия смеси на входе в камеру смешения.

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

ГОп =0хвёрхв [(^ (Рп ) — ^о ) — Мн ] (16)

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

Диффузор

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

Приведем дифференциальные уравнения, на основе которых выполнялся расчет диффузора:

dI dFд

—=р~г-V ; dz dz

dDтт

(17)

= жСРжРж К — Т I (тж — Т ) Хк/((гРп )2 ^ );

(18)

::

+

1

|

А

(

w

= 0,

(19)

где Гд = Дг); I = Gw + pFg — полный импульс потока; тж — касательное напряжение на стенке канала; ск — коэффициент, учитывающий реальность процесса конденсации.

Определив величину полного импульса по координате г, находим текущее значение статического давления р = ( I -

Результаты расчета инжектора

На рис. 1 показано изменения основных физических величин вдоль продольной оси инжектора, полученных в результате численного решения систем уравнений.

Из этого рисунка видно, что давление горячей воды на выходе из инжектора в пять раз превышает давление сухого насыщенного пара на входе в сопло. Между тем наши опыты показали (рис. 2), что максимальное превышение давления — около 3,6. Такое расхождение можно объяснить упрощенной моделью процесса, в которой не учитывались дробление и коагуляция капель, полидисперсность капельного потока, наличие жидкой пленки и срывные явления. Однако расчеты, выполненные на основе пред-

2

4

6

8

Рис. 2. Степень повышения давления в инжекторе в зависимости от коэффициента инжекции и

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

и

V Сопло Лаваля А Камера смешения Диффузор

-- —^^

¡,, ми.. м/с Г-Ю2 К

1

У 'л

\ Зкк

-----р "'хк / "'"ем.

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

Применение инжектора в отопительных системах

Инжектор сочетает в себе компактный высокоэффективный теплообменный аппарат, где в процессе конденсации пара на холодных каплях пассивной жидкости коэффициент теплоотдачи достигает сотен тысяч Вт/(м2К) (при конденсации пара на холодных стенках а = = 7000—9000 Вт/(м2К)), и насос, позволяющий в несколько раз повысить давление по сравнению с давлением активного агента (пар) на входе в инжектор. Поэтому естественно, что им можно заменить подогреватель сетевой воды, используемый для отопления зданий, и насос, обеспечивающий прокачку этой воды (что, как уже отмечалось, было предложено профессором В.В. Фисенко).

Далее в данной статье рассматривается использование инжектора в котельной СПбГПУ с целью экономии топлива и электрической энергии.

Последовательно определяем: необходимую производительность парового котла (парогенератора) (исходя из уравнения теплового баланса при заданном расходе сетевой воды Gс в)

(п = (св 0'в5 1в8 )/(1П1 1в2 );

тепловой поток, получаемый потребителем,

бпот = Gс в(гВ5 — гв8);

расход топлива в парогенераторе

В = ^п (*п1 1в2)КН1 Ппг);

среднюю температуру воды в водоподогре-вателе

Т ср = (Т5 Т )/ 1п (Т5/ Т8);

среднюю температуру пара в водоподогрева-теле

Тп ср = 0п1 - 'в2 У( 5п1 _ 5в2 );

температурный напор в водоподогревателе АТ = Т - Т .

в п п ср в ср

Так как инжектор представляет собой открытую адиабатическую систему, то можно написать

(пИп1 81в8 = °СВ1В5 , где Опи, О8 - расход пара и

воды через инжектор.

Вводя коэффициент инжекции и = G8/ Gп и =

G в = G8 + Gn и, определим расход пара через инжектор Gn и = Gc в/(1 + и) и расход топлива в схеме с инжектором Вти = BTJfiпи jGп .

Сравнение тепловой схемы с инжектором и обычной — с водоподогревателем сетевой воды

В существующей тепловой схеме перегретый пар после парогенератора дросселируется (/п1 = = const), а затем поступает в водоподогреватель (температура пара снижается с Тп1 = 280 °С до Тп1д = 270 °С, а давление — с 29 до 19 бар.). Средняя температура пара в водоподогревателе понизится и станет равной

Тп ср д = (гп1 - 'в2д V(^п1д _ ^в2д )'

что приведет к снижению температурного напора AT = Т - Т .

вп д п ср д в ср

Если считать, что при этом коэффициент теплопередачи не изменится, то это приведет к увеличению расхода пара и топлива при неизменном Q

^пот

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

G = GAT /AT ; В = G В/G ,

п д п вп вп д' т д п д т' и'

п д т п'

(1п1 1в5)

0в5 1в8)

и учитывая, что расход сетевой воды

где АТВп д = Тп1д ТВср.

Расход топлива для схемы с инжектором определим по формулам Вти = Ви0п, Втид = Втдрпи /Gпд. Экономия пара и топлива составит

А(п = (п -(пи; А(пд = (пд -(пи;

АВТ = Вт -Вти; Щд = ВТд -Вти.

Расчет тепловых схем был выполнен при следующих данных:

G = 70 кг/с; Т , = 553 К; Т , = 543 К;

св ' ' п1 ' п1д

/п1 = 2940 кДж/(кг- К); = 994,4 кДж/(кг К); 4д = 896,6 кДж/(кг- К)"^ = 398 кДж/(кг К);

iв8 = 209 кДж/(кг- К); Лн = 0,7; Я = 48819 кДж/кг; Арсв = 3 • 105 Па; Т8 = 323 К; Т5 = 368 К; = 0,833; = 16,34 кг/кмоль; р = 1,0133 • 105 Па- Т = 273 К;

^ ос ' 'ос '

5п1 = 6,46 кДж/(кг • К); 5п1д = 6,66 кДж/(кг • К); 5вю = 2,63 кДж/(кг • К); ^ = 2,42 кДж/(кг • К);

В результате расчета имеем: О = 6,798 кг/с; Т = 345 К; N = 30759,67 Вт;

п ' ' ' вср ' н

б = 13230 Вт; В = 0,325 кг/с; В = 0,232 кг/с;

пот т ти

О = 7,14 кг/с; АО = 0,342 кг/с; и = 13,45;

пд пд

Оп и=4,84 кг/с; АОп = 1,955 кг/с; АВт = 0,0935 кг/с; ™АВ = 336,76 кг/ч; АВ = 8082,2 кг/сут.

т ч ' ' ' т с ' / •>

Потребитель

Рис. 3. Тепловая схема водоподогревателя сетевой воды с насосом

Новая схема

Потребитель орячей воды

Потребитель теплоты

Линия запуска инжектора

Бак для слива воды

Линия подпитки холодной водой

Насос

Рис. 4. Тепловая схема подогрева сетевой воды с инжектором

При нагреве сетевой воды от 50 до 95 °С с расходом 70 кг/с и установке инжектора экономится пара ДGп = 6,798 - 4,843 = 1,95 кг/с и топлива ДВт = 0,3253 - 0,2317 = 0,0936 кг/с (ДВт -100/Вт = 28,7 %). Таким образом, в сутки (даже при отсутствии дросселирования) экономия топлива по массе составит 8082 кг, а по объему газа, приведенного к нормальным условиям, ДВ = ДВ и Т /Р /р = 0,128 м3/с или

' ту ту г к ос' ос' г г ' '

11079 м3/сутки.

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

соса в старой схеме Nн = GсвАРсв / Рсв / пн = 30759 Вт. Экономия в сутки составит 2,66-109 Дж.

Рассмотрим процессы в водоподогревателе и инжекторе с позиций энтропийного и эксер-гетического методов. Для водоподогревателя применим энтропийный метод, как более удобный для анализа, а для инжектора - эксергети-ческий. На основании теоремы Гюи-Стодолы определим потерю потока эксергии для водоподогревателя: ( (т \ \

Ат а =То с = То с

= 3531,6 кВт.

°с вСр 1п

Т.

—G п ( — )

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

D = Ех. + Ех2 - Ех5 = 2816 кВт,

и 1 2 5 '

где Ел, = G ((/ , — I ) — Т (у , — я )); Ех8 = G8

^ 1 п и ^ п1 о с о с ^ п1 ос''' 8 8

(('в8 - 'о с) - То Х -Яос)); Щ = ^ в (('в5 - 'о с) -

- Т (у 5 -я )).

о с в5 ос

Таким образом, применение инжектора уменьшает потерю потока эксергии в процессе нагрева воды на ДD = Dт а — Dи = 715,2 кВт, или на 20 %.

Выводы

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

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

Термодинамический анализ показывает, что с точки зрения экономии топлива и электрической

энергии применение тепловых схем с инжекторами в отопительных системах значительно эффективнее по сравнению с существующими.

Приложение

Принятые обозначения: ав — термодинамическая скорость звука, м/с; ср — удельная изобарная теплоемкость, Дж/(кгК); с{ - коэффициент трения; сх - коэффициент сопротивления; F, f — площадь поперечного сечения, м2; G - расход, кг/с; I — удельная энтальпия Дж/кг; j — плотность потока массы, кг/(м2с); Д - результирующая плотность потока пара от испаряющейся капли, кг/(м2с); I - длина, м; М - число Маха; N — мощность, Вт; п — поток капель, 1/с; пу - объемная концентрация капель, 1/м3; р - статическое давление, Па; Q — тепловой поток, Вт; R - радиус, м; R п - газовая постоянная пара, Дж/(кгК); г — удельная теплота парообразования, Дж/кг; Т — температура, К; t—температура, °С; V — удельный объем, м3/кг; w — скорость, м/с; х — степень сухости пара; г — продольная координата, м; а — коэффициент теплоотдачи, Вт/(м2К); I — коэффициент теплопроводности, Вт/(мК); р — коэффициент динамической вязкости, Па-с; р — плотность, кг/м3; t — время, с; — касательное напряжение трения на стенке, Па.

Подстрочные индексы: г - горло; г к - горячие капли; д - диффузор; ж - жидкость; к - капля; кр - критический; КС - камера смешения; п -пар; с - сопло; с д - дивергентная часть сопла; с к - конвергентная часть; см - смесь; ср - срез; с ц - цилиндрическая часть сопла; х к - холодные капли.

СПИСОК ЛИТЕРАТУРЫ

1. Эллиот Д. Исследование применения насосов с несущим газом в ракетных двигателях // Вопросы ракетной техники. 1964. №4. С. 88—112.

2. Данилин В.С., Дейч М.Е., Севостьянов А.П. [и др.] Экспериментальное исследование пароводяного инжектора // Симпозиум по производству электроэнергии с помощью МГД- генераторов, Варшава, 24-30 июля 1968. Доклад SM-107/143,

3. Богомолов Б.Г., Духовлинов С.Д., Черных Е.В. [и др.] Некоторые результаты изучения однокомпо-нетной схемы жидкометаллического МГД-преобразования энергии // Симпозиум по производству электроэнергии с помощью

МГД-генераторов, Варшава. 24—30 июля 1968. Доклад SM-107/135.

4. Шпильрайн Э.Э., Якимович К.А. Термодинамика энергетической МГД- установки с парожид-костным инжектором // ТВТ. 1965. №5.

5. Калафати Д.Д., Козлов В.Б. Термодинамика циклов со смешением жидкого металла перед МГД-преобразователем // ТВТ. 1963. №3.

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

7. Radebold [etc]. Energy conversation with liquid-metal working-fluids in the MHD-Staustrahlrohr // Electricity from MHD. IAEA. Vienna. 1968. Vol. 3. P. 1439-1461.

8. Баранов Г.А., Барилович В.А., Батуев Ш-Б.Б., Зысин В.А. Исследование процессов в разгонном устройстве, работающем на вскипающей жидкости // ТВТ. 1972. Т. 10. №3. С. 629-634.

9. Зысин В.А. [и др.] Вскипающие адиабатные потоки. М.: Атомиздат, 1976. 152 с.

10. Барилович В.А. [и др.] Расчет жидкостного

инжектора // Труды ЛПИ. 1980. №370. С. 87-92.

11. Новожилов И.А., Фисенко В.В. Новая энергоресурсосберегающая технология // Энергетик. 1996. №3. С. 4-5.

12. Кириллов И.И., Яблоник Р.М. Основы теории влажнопаровых турбин. Л.: Машиностроение, 1968. 264 с.

REFERENCES

1. Elliot D. Issledovaniye primeneniya nasosov s nesushchim gazom v raketnykh dvigatelyakh [Research of use of pumps with the bearing gas in rocket engines.]. Voprosy raketnoy tekhniki. 1964. №4. S. 88-112.(rus)

2. Danilin V.S., Deych M.Ye., Sevostyanov A.P. [i dr.] Eksperimentalnoye issledovaniye parovodyanogo inzhektora [Experimental investigation of vapor-water injector]. Simpozium po proizvodstvu elektroenergii s pomoshchyu MGD- generatorov. Varshava, 24-30 iyulya 1968. SM-107/143 (rus)

3. Bogomolov B.G., Dukhovlinov S.D., Chernykh Ye.V. [i dr.] Nekotoryye rezultaty izucheniya odnokomponetnoy skhemy zhidkometallicheskogo MGD-preobrazovaniya energii [Some results of investigation of one-component scheme of liquid-metal MHD energy conversation]. Simpozium po proizvodstvu elektroenergii s pomoshchyu MGD- generatorov. Varshava, 24-30 iyulya 1968. SM-107/135. (rus)

4. Shpilrayn E.E., Yakimovich K.A. Termodinamika energeticheskoy MGD- ustanovki s parozhidkostnym inzhektorom [Thermodynamics of MHD power installation with vapor-liquid injector]. TVT. 1965. №5. (rus)

5. Kalafati D.D., Kozlov V.B. Termodinamika tsiklov so smesheniyem zhidkogo metalla pered MGD-preobrazovatelem. [Thermodynamics of cycles with

liqud-metal mixing before MHD-convertor]. TVT. 1963. №3. (rus)

6. Sokolov Ye.Ya., Zinger N.M. Struynyye apparaty [Jet apparatus]. M.: Energoatomizdat, 1989. 351 s. (rus)

7. Radebold [etc]. Energy conversation with liquid-metal working-fluids in the MHD-Staustrahlrohr. Electricity from MHD,IAEA, Vienna. 1968. Vol. 3. P. 1439-1461.

8. Baranov G.A., Barilovich V.A., Batuyev Sh-B.B., Zysin V.A. Issledovaniye protsessov v razgonnom ustroystve, rabotayushchem na vskipayushchey zhidkosti [Research of processes in the accelerating device working at the boiling liquid.] TVT. 1972. T.10. №3. S. 629-634. (rus)

9. Zysin V.A. [i dr.] Vskipayushchiye adiabatnyye potoki [Boiling adiabatic streams]. M.: Atomizdat 1976. 152 s. (rus)

10. Barilovich V.A. [i dr.] Raschet zhidkostnogo inzhektora [Calculation of liquid injector]. Trudy LPI. 1980. №370. S. 87-92. (rus)

11. Novozhilov I.A., Fisenko V.V. Novaya energo resurs osberegayushchaya tekhnologiya [New power resource-saving technology]. Energetik. 1996. № 3. S. 4-5. (rus)

12. Kirillov I.I., Yablonik R.M. Osnovy teorii vlazhnoparovykh turbin [Basics of wet vapor turbines]. L.: Mashinostroyeniye, 1968. 264 s. (rus)

СВЕДЕНИЯ ОБ ABTOPAX/AUTHORS

БАРИЛОВИЧ Владимир Антонович — доктор технических наук профессор Санкт-Петербургского

политехнического университета Петра Великого.

195251, Россия, г. Санкт-Петербург, Политехническая ул., 29.

E-mail: vba_exerg@mail.ru

BARILOVICH Vladimir A. — Peter the Great St. Petersburg Polytechnic University. 29 Politechnicheskaya St., St. Petersburg, 195251, Russia. E-mail: vba_exerg@mail.ru

СМИРНОВ Юрий Александрович — кандидат технических наук доцент Санкт-Петербургского

политехнического университета Петра Великого.

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

195251, Россия, г. Санкт-Петербург, Политехническая ул., 29.

E-mail: y_a_s@mail.ru

SMIRNOV Yurii A. — Peter the Great St. Petersburg Polytechnic University. 29 Politechnicheskaya St., St. Petersburg, 195251, Russia. E-mail: y_a_s@mail.ru

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