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

электронное научно-техническое издание

НАУКА и ОБРАЗОВАНИЕ

________Эл № ФС 77 - 30569. Государственная регистрация №0420900025. ТЭБЫ 1994-0408_

Теплофизические аспекты повышения мощности активных зон реакторов с водой под давлением

# 05, май 2009 автор: Сатин А. А.

Первые

публикации

Программы и программные системы

Учебные

программы

Студенческая

Общие

проблемы

инженерного

образования

Инженер в современной России

Экобионика

Зарубежное

образование

История

технического

прогресса

Будущий

инженер

Вне рубрик

Расширеный поиск Подписаться на новости

ПОИСК

Ред. совет Специальности Рецензентам Авторам Архив

[email protected] МГТУ им. Н.Э. Баумана

Введение

Повышение тепловой мощности активных зон является важным для повышения конкурентоспособности действующих и проектируемых АЭС. Над данной задачей работают многие коллективы организаций атомной отрасли. Предложены и находятся в большой степени проработанности четыре направления повышения тепловой мощности.

■ Увеличение загрузки топлива в активной зоне. Например, высота топливного столба в ТВС-2М по сравнению с ТВС-2 увеличена на 150 мм за счет изменения конструкции хвостовика и головки [1].

■ Снятие избыточного консерватизма в теплогидравлических и нейтронно-физических расчетах позволило обосновать возможность повышения тепловой мощности до 104% за счет уменьшения коэффициента запаса по линейной мощности [2].

■ Введение в кассеты интенсификаторов тепломассообмена с тем, чтобы повысить значения критического теплового потока, увеличить теплосъем с ТВС [1].

■ Изменение диаметров твэлов, что требует большого объема НИОКР [2].

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

Кризис кипения — явление характеризуемое внезапным падением коэффициента теплоотдачи вследствие изменения механизма кипения, что приводит к росту температуры поверхности нагрева. Тепловой поток непосредственно перед наступлением кризиса называется критическим тепловым потоком [3].

На кризис кипения при вынужденном движении потока жидкости влияние оказывают [4] пузырьковый слой, движущийся параллельно поверхности нагрева, который так близко примы

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

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

1. Теоретические подходы к изучению явления кризиса кипения

Одной из первых моделей, позволяющих описать механизм возникновения кризиса кипения, стала модель баланса сил, действующих на пузырь предложенная П. Ченгом [4].

СОБЫТИЯ

Международная

научно-практическая

конференция

"Использование технологий дистанционного обучения в высшей школе "

III Международная

научно-практическая

конференция

"Информационная среда вуза XXI века"

НОВОСТНАЯ ЛЕНТА

24.08.2009

Комсомольская правда >

24.08.2009 > Евгений СТРИГУНОВ > Толковым студентам -платникам государство будет помогать

13.07.2009

РИА "Новости" > 13.07.2009 > Около 28 тысяч студентов вузов переведены с платных на бюджетные места в 2009 году

8.07.2009

Письмо Рособразования от

29.05.2009 N20-353 "О Конкурсе 2009 года на соискание медалей Российской академии наук с премиями для молодых ученых РАН... "

1.07.2009

Интерфакс > 01.07.2009 > Глава Рособрнадзора: Процедура аккредитации вузов претерпит изменения

24.06.2009

Интерфакс > 23.06.2009 > Путин: при подготовке нового бюджета расходы на науку и образование должны быть оптимизированы

весна

ВХОД

регистрация забыли пароль?

Пресс-релизы

Библиотека

Конференции

Выставки

Доска объявлений

Архив

Ассоциация технических Университетов Информация о проекте Авторы

Координационный совет

Рис. 1 . Баланс сил, действующих на прикрепленный к стенке пузырь

На пузырь, прикрепленный к вертикальной стенке (рис.1) действуют следующие силы:

4 лгі

объем и радиус пузыря в момент отрыва, Ал.— плотность насыщенного пара, £ — ускорение свободного падения

сила тяжести

о V Г

э , где мл. — масса насыщенного пара внутри пузыря, к? / кр —

• сила Архимеда жидкости;

4 лг1

= Рж%

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

Р =С(е)а1т = С(в)2лг,„сг

г V } у? V / ^ ' где ^ — I

коэффициент, зависящий от краевого угла

сила поверхностного натяжения 1,7 / ‘ V где — длина контакта пузыря и стенки в

момент отрыва, коэффициент поверхностного натяжения, смачивания;

сила инерции, с учетом поправки на неинерциальность системы пузырь-поток жидкости,

Рцх = Рг + Ря =тжаг+тж й;

Алгі

Рг =тжа, = рж — = — ржг и

и , где ?— тангенциальная составляющая силы инерции

Плит 4 Л" 2 2

нормальная составляющая силы инерции

=тжаъ =

4л-

2 2

Д1, — тангенциальная и нормальная составляющие ускорения, и “ш — составляющие относительной скорости

между жидкостью и пузырем, нормальная и параллельная к стенке соответственно.

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

Ра + ^-Рт = О.

= О-

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

(1)

Ченг [4] представил критический тепловой поток как сумму количества теплоты, переносимой за счет конвекции жидкости, и скрытой теплоты парообразования, переносимой пузырями

= + РаРшжУжжГж!*, (2)

скс3_ ' '

^ %идО!^ = % ~ Т

: — температура насыщения, 7 — температура в месте кризиса, — теплота парообразования. Скорости ^ус и

В формуле (2) — критический тепловой поток,

постоянные, нормальные к стенке скорости в

момент отделения пузыря жидкости и насыщенного пара соответственно, млЭдгр - _ температура недогрева,

Т.

V.

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

* ♦ +

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

' \| п

пузырем. Последняя определяется через число Вебера и радиус пузыря в момент отрыва и имеет вид \ .

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

Важным этапом для описания механизма возникновения кризиса стало появление модели критического перегрева в пристенном слое жидкости Л. Тонга [4]. Тонг рассмотрел баланс энергии для слоя перегретой жидкости, отделенной от основного потока пузырьковым слоем (рис. 2) и представил его в виде

(3)

В уравнении (3) Р — средняя плотность слоя жидкости, средняя скорость слоя жидкости, — периметр

нагревателя, — толщина слоя жидкости, — средняя энтальпия слоя жидкости, ^ — энталь

жидкости, С£— коэффициент теплоотдачи от слоя жидкости пузырьковому слою, — средняя удельная теплоемкость

перегретой жидкости, # — тепловой поток от поверхности, г — расстояние.

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

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

Решение уравнения (3) для состояния, близкого к кризису кипения при заданных параметрах (т. е. давлении, массовом расходе, паросодержании и эквивалентном диаметре) и дополнительных упрощениях для однородного (4) и неоднородного (5) теплового потока имеют вид

а

(4)

(5)

Из предположения, что начало кризиса кипения определяется предельным значением энтальпии перегретого слоя, Тонг получил соотношение для эквивалентного однородного теплового потока при кризисе кипения

с

(6)

с = -

а

Здесь коэффициент

Р^ш.ж'^агж'- р зависит от режима течения и паросодержания, что в принципе можно свести к

зависимости от скорости а.ж, которую можно определить из зависимостей Ченга о связи отрывного радиуса пузыря и скорости течения жидкости, 2=1— положение точки начала кризиса кипения, изме

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

С

$ га .т?дн

-а

е ) а

(7)

Поправочный коэффициент в общем случае также является функцией паросодержания. Соотношения для его определения приведены в книге Тонга [4], а также в работах наших исследователей В. С. Осмачкина [5], В.Н. Смолина [6,7].

2. Эмпирические соотношения для расчета критического теплового потока в реальных конструкциях. Табличный метод ФЭИ.

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

■ сложность геометрии реальных конструкций ТВС, представляющих собой пучки стержней в треугольной (ВВЭР) либо

или

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

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

разных ячейках;

■ неравномерное энерговыделение по сечению и длине активной зоны, разный расход, перетечки теплоносителя из ячейки в ячейку;

■ влияние дистанционирующих и перемешивающих решеток, различающихся по конструкции и имеющих разные гидравлические характеристики.

Рис. 3 . Пучок с треугольной упаковкой твэлов.

1- ячейка образована твэлами, 2 - ячейка образована твэлами и направляющим каналом (НК), 3 - периферийная ячейка

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

паросодержание. Например, зависимость, предложенная в 1951 г. Аргонской Национальной лабораторией (А1\11_) для

АТ

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

массовой скоростью С* [4] , имеет вид

^-=с(— І Д7°'п 10е I 1

(8)

Значения параметров С и УУ1, при которых данная зависимость верна, приведены в таблице 1. Таблица 1. Значения параметров с и т в зависимости АМЬ

* МПа 3,5 7 14

0,817 0,626 0,445

т 0,160 0,275 0,500

Подобные зависимости были получены и для расчета активных зон реактора ВВЭР. Например в [9] предложена зависимость

_ (21,20-0,44^)0^-4,4 1(Г3СЦ Чщ (0,898 - 0, ОЗІЗр)-1 - 84,5 +1 /

(9)

где Р — давление, X— массовое паросодержание, & — массовая скорость, Ь — длина пучка, ^Т ^ ^ — тепловой диаметр сборки, Р — площадь сечения сборки; И — обогреваемый периметр сечения сборки. Как и все эмпирические зависимости, зависимость (9) применима только в определенном диапазоне значений входящих в нее параметров, представленных в таблице 2.

Таблица 2. Диапазоны изменения параметров в соотношении (9).

Р , МПа Ст , кг/м2 час X ^ , мм Ь , м

10,3 — 18, 6 530 — 3230 - 0,07— 0,4 9 1— 3,5

В расчетах ОАО ОКБ «Гидропресс» применяется соотношение, предложенное Ю.А. Безруковым, и приведенное в [9]

^ =0,795(1- хПО)*1 (1-0,0185р)

! = 0,311(1 - х) - 0,127 ^ п= 0,105;; - 0,5 _

Таблица 3. Диапазоны изменения параметров в формуле (10)

(10)

Р , МПа Сг , кг/м2 час X сі, мм Ь , м

7,45—16,7 700 —3800 - 0,07— 0,4 9 ,5 3, —3 ,7

В результате экспериментов по определению кризиса кипения в различных конструкциях были накоплены данные для широкого диапазона изменения параметров, а с вводом соответствующих поправок удается получить величины критического потока для любых вариантов этих параметров. Данная работа велась в 70-х гг. прошлого века в США. В настоящее время такие же работы проводятся в России. Недавно опубликован табличный метод, представляющий собой массивы (таблицы) экспериментальных данных для сборок из базы данных ФЭИ [12]. Эти массивы характеризуются следующими параметрами:

■ число экспериментальных точек 963;

■ диаметр твэлов от 10,2 до 15 мм;

• относительный шаг з!от 1,15 до 1,88;

■ число твэлов в сборках 4, 9, 16 и 20;

■ геометрия сечения обечайки — прямоугольник;

• относительная длина сборки (отношение длины пучка твэлов Ь к тепловому диаметру ^т) 40—960;

■ все твэлы обогреваемые;

■ обогрев равномерный по длине и сечению;

■ диапазон давления Р =4 — 15,6 МПа;

О =115— 4810 кг/(м2с);

диапазон массовой скорости диапазон паросодержания от X = - 0,17 до 0,87.

Опубликованные в [12] таблицы содержат значение критического теплового потока , для которого

■ тепловой диаметр регулярной ячейки “ I- ■

=11лш

относительный шаг 1 относительное расстояние от места кризиса до входа в сборку / > 3300Л4М ;

коэффициент гидравлического сопротивления дистанционирующих решеток ~ ^ ; расстояние от места кризиса до ближайшей дистанционирующей решетки 2 = 0;

отношение максимального теплового потока к среднему по сечению ^Ям ^ Яс ~ ^ .

Для расчета критического теплового потока по табличному методу ФЭИ при равномерном по длине энерговыделения используется соотношение

В соотношении (11)

. ^ = ^и) = №1/11)-1/3_

(11)

поправочная функция на тепловой диаметр ячейки;

поправочная функция на относительный шаг;

где

Къ = Р(Г) = 1 + 0,6 ехр(-0,01// ¿у _

до входа в сборку;

, Л^ЭЯ^ГОЛООО)0"

решетки (переменная ' / 4 У,

— поправочная функция на относитель

— поправочная функция, учитывающая влияние дистанциони

, Л=0,Э^“-5С0/1000)°':!

(переменная r f ' /

■ ^-5 ^к) ^,б ■' ^а) — поправочная функция на теплогидравлическую неравноцен

сборки, равная единице в случае одинаковых ячеек;

К6 = 1 + 5(АГ- 1)(0,5-0,04р)

сборки;

К*

Поправочная функция <5 сформирована и оптимизирована авто

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

3. Анализ путей повышения величины критического теплового потока с использованием табличного метода ФЭИ.

Повышение мощности активной зоны требует повышения величины критического теплового потока. Конструкции ТВС

активных зон отечественных и зарубежных водо-водяных реакторов имеют определенные размеры твэлов, шаг их

ТУ' Vе тг

расположения, тепловые диаметры ячеек, однозначно определяя величины 1, л2, л3, из зависимости (11).

Единственной возможностью повышения величины критического теплового потока остается варьирование параметрами К

функции 1 4 — шагом расположения решеток и коэффициентом гидравлического сопротивления.

Например, для квадратной ячейки ТВС, типичной для реактора PWR (рис. 4), средний равномерный тепловой поток с q„ = 522 кВт/ м2 п „ „

поверхности твэлов . Допускается наличие максимально напряженной ячеики с равномерным по

высоте тепловым потоком (ТП)

q = 1,35# = 704,7 кВт/м2

При повышении тепловой мощности активной зоны на

30%, тепловой поток в рассматриваемой ячейке станет равен .2

■7 = 1,3-1,35,? =916 кВт / м

. Массовая скорость

^ — 2000 КГ(м • С) не меняется, Тогда при наличии только обычных дистанционирующих решеток, расположенных с

К* = 0,65

шагом 500 мм и имеющих

. наблюдается кризис кипения на расстоянии 3480 мм от низа активной зоны.

Рис. 4 . Квадратная ячейка PWR Возможны два варианта повышения величины КТП по всей высоте активной зоны:

1. Варьирование шагом решеток.

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

Рис. 5 . График зависимости КТП от шага расположения решеток

Уменьшение шага с 500 мм до 250 мм (рис.5) при ~ позволяет поднять КТП на 9% по сравнению с

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

2. Варьирование коэффициентом гидравлического сопротивления.

В этом случае количество решеток не меняется, не увеличивается доля металла в активной зоне, но меняется величина гидравлических потерь.

£

Рис. 6 . График зависимости КТП от коэффициента гидравлического сопротивления /

К.

Увеличение / с 0,65 до 1,2 при шаге 250 мм (рис.6) позволяет поднять КТП еще на 3,4% по отношению к варианту с

шагом расположения 250 мм и ~ ^ . Хотя количество решеток в этом случае не увеличивается, но возрастает их

гидравлическое сопротивление, а, следовательно, и гидравлическое сопротивление активной зоны.

Заключение.

Одним из ограничений увеличения мощности как работающих, так и проектируемых ТВС служит кризис теплообмена. Кризис теплообмена в пучке ТВС имеет локальный характер и в настоящее время рассчитывается по моделям, обобщающим экспериментальные данные, полученные на электрообогреваемых имитаторах ТВС, учитывающих особенности конкретных пучков твэлов. Но как правильно заметили авторы [12], в ряде случаев полученные результаты стоит рассматривать как оценочные и нуждающиеся в подтверждении экспериментом. Тем не менее, было установлено, что дистанционирующие, перемешивающие решетки позволяют существенно увеличить значение критического теплового потока. Выбором положения решетки, их характеристик и количества можно повысить тепловую мощность ТВС принятых для ВВЭР-1000 размеров на 15-20%.

Список используемой литературы.

1. Драгунов Ю.Г., Рыжов С.Б., Васильченко И.Н., Кобелев С.Н., Вьялицын В.В. Проект активной зоны для РУ АЭС-2006. Материалы конференции: Обеспечение безопасности АЭС с ВВЭР. ФГУП ОКБ ГИДРОПРЕСС, 2007.

2. Асмолов В.Г. Атомная энергетика - второе рождение. Научный совет РНЦ «Курчатовский институт», М, 2007.

3. Петухов Б.С., Генин Л. Г., Ковалев С. А., Соловьев С. Л. Теплообмен в ядерных энергетических установках. М., Изд-во МЭИ, 2003. 548с.

4. Тонг Л. Теплоотдача при кипении и двухфазное течение. М., «Мир», 1960. 280с.

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

Препринт ИАЭ — 2014. М., 1970.

6. Смолин В. Н. Модель механизма кризиса теплоотдачи при движении пароводяной смеси и методика расчета кризисных условий в трубчатых твэлах.— В кн.: Исследование критических тепловых потоков в пучках стержней. Семинар СЭВ ТФ—74. М., 1974, с. 209—224.

7. Смолин В. Н., Поляков В. К. Критический тепловой поток при продоль Теплоэнергетика, 1967, № 4, с., 54-58.

8. Дубровский И. С., Югай Т., Пащенко М. П., Чалых А. Ф. Кризис теплообмена при вынужденном течении пароводяной смеси в сборке стержней в стационарных и нестационарных режимах. В кн.: Исследование критических теп пучках стержней. Семинар СЭВ ТФ—74. М., 1974, с. 67—78.

9. Безруков Ю. А., Астахов В. И., Салий Л. А. и др. Исследование критических тепловых потоков в пучках стержней приме

стержней. М., 1974, с. 57—66.

10. Емельянов И.Я., Солонин В. И., Михан В.И., Рекшня Н.Ф., Конструирование ядерных реакторов. М. Энергомиздат.1982.

11. Кириллов П.Л. Справочник по теплогидравлическим расчетам (ядерные реакторы, теплообменники, парогенераторы).

М. Энергоатомиздат. 1984.

12. Бобков В. П., Ефанов А. Д. и др. Таблица для критических тепловых потоков в сборках твэлов в квадратной упаковке. М.

Атомная энергия: том 98, вып.1: 2005.

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

Публикации со словами: активная зона ядерного реактора, реактор с водой под давлением, кризис кипения. критический тепловой поток

Тематические рубрики:

Наука в образовании: Электронное научное издание

Публикации с ключевыми словами: активная зона ядерного реактора. реактор с водой под давлением. кризис кипения. критический тепловой поток

Публикации со словами: активная зона ядерного реактора, реактор с водой под давлением. кризис кипения. критический тепловой поток

Тематические рубрики:

• Наука в образовании: Электронное научное издание

Ассоциация технических Университетов Вузы

Информационное агентство

Координационный совет Новости УМО Вузов

телефон (8499) 263-68-67 ; 1*55

□ 5ТИСК ОЯОиР

© 2003-2009 «Наука и образование: электронное научно-техническое издание»

т