CC BY
62
УДК 621.039.542:[532.546.6+532.5.031]
НЕРАВНОМЕРНОСТИ РАСПРЕДЕЛЕНИЯ ПАРОВОДЯНОГО ТЕПЛОНОСИТЕЛЯ В ЗАСЫПКЕ ТЕПЛОВЫДЕЛЯЮЩИХ ЧАСТИЦ Канд. техн. наук, доц. СОРОКИН В. В.
Белорусский национальный технический университет E-mail: [email protected]
Существенное повышение безопасности атомной энергетики может быть достигнуто за счет внедрения микротвэльного топлива в реакторных установках с водяным теплоносителем. Микротвэл представляет собой шарик диаметром порядка миллиметра, состоящий из топливного ядра, покрытого защитной оболочкой. Оболочки разделяют топливо и теплоноситель и обеспечивают удержание продуктов деления ядер. Благоприятные теплофизиче-ские свойства засыпки микротвэлов в сочетании с водяным теплоносителем позволяют избежать аварийных последствий перегрева топлива. Известны предложения по использованию микротвэлов в вВэР, прямоточном реакторе с перегревом пара, реакторе со спектральным регулированием за счет изменения паросодержания.
В статье расчетным методом анализируется возможность возникновения режимов с неравномерностью распределения теплоносителя в засыпке. Установлена неустойчивость, связанная с неоднозначностью зависимости перепада давлений от расходного паросодержания. Неустойчивость уменьшается с увеличением массового расхода, размера частиц, со снижением абсолютного давления, а влияние пористости на эффект мало. Наиболее опасный диапазон - при паросодержании порядка единицы. Возникают поры, заполненные только паром, и поры с водопаровой фазой при паросодержании, меньшем среднего по засыпке. В прямоточном реакторе с перегревом пара в таком режиме на выходе теплоноситель будет содержать долю влаги в капельной форме. Найдено, что в части режимов с недогревом на входе движение теплоносителя неустойчиво против сильного возмущения течения. В реакторе радиального типа при прочих равных условиях неравномерность распределения двухфазного теплоносителя по засыпке выше по сравнению с однофазным теплоносителем.
Ключевые слова: пароводяной теплоноситель, засыпка, тепловыделяющие частицы, ядерный реактор.
Ил. 3. Библиогр.: 11 назв.
NONUNIFORMITIES OF TWO-PHASE COOLANT DISTRIBUTION IN A HEAT GENERATING PARTICLES BED
SOROKIN V. V.
Belorussian National Technical University
Sufficient atomic power generation safety increase may be done with microfuel adapting to reactor plants with water coolant. Microfuel particle is a millimeter size grain containing fission material core in a protecting coverage. The coverage protects fuel contact with coolant and provides isolation of fission products inside. Well thermophysical properties of microfuel bed in a direct contact with water coolant excludes fuel overheating when accidents. Microfuel use was suggested for a VVER, а direct flow reactor for superheat steam generation, a reactor with neutron spectra adjustment by the steam partial content varying in the coolant.
Nonuniformities of two-phase coolant distribution in a heat generating particles bed are predicted by calculations in this text. The one is due to multiple-valuedness of pressure drop across the bed on the steam quality dependency. The nonuniformity decreases with flow rate and particle size growths absolute
pressure diminishing while porosity effect is weak. The worse case is for pressure quality of order of one. Some pure steam filled pores appears parallel to steam water mixture filled pores, latter steam quality is less than the mean of the bed. Considering this regime for the direct flow reactor for superheat steam generation we predict some water drops at the exit flow. The two-phase coolant filtration with subcooled water feed is unstable to strong disturbance effects are found. Uniformity of two-phase coolant distribution is worse than for one-phase in the same radial type reactor.
Keywords: steam-water coolant, fillings, heat-emitting particles, nuclear reactor.
Fig. 3. Ref.: 11 titles.
Задача расчета возможной неравномерности распределения пароводяного теплоносителя в засыпке тепловыделяющих частиц возникает в связи с созданием насыпных активных зон ядерных реакторов с микротвэльным топливом. Такая техника разрабатывается с целью кардинального улучшения безопасности атомной энергетики.
Совокупность требований к ядерному топливу, теплоносителю и другим материалам активной зоны настолько ограничена, что на сегодняшний день системе «вода - цирконий - оксид урана (плутония)» в активной зоне мощной коммерческой атомной электростанции трудно предложить альтернативу. К сожалению, эта система уязвима при авариях с нарушением отвода теплоты от топлива. Уже при 900 °С прочность циркониевого сплава снижается, заметно увеличивается выход водорода вследствие реакции циркония с водой, прогрессирует течь радиоактивных веществ. Рабочая температура в центре некоторых таблеток топлива превышает 2000 °С, так что угроза реализации опасного перегрева оболочки не умозрительна. Причин, которые могут привести к нарушению отвода теплоты от топлива, много: течи теплоносителя, отключение электроэнергии, кризис теплоотдачи. Ситуация усугубляется невозможностью мгновенного прекращения генерации теплоты в активной зоне, в течение времени порядка минуты мощность можно снизить на порядок, остаточная тепловая энергия на уровне 40 МВт выделяется месяцами. Сценарии инцидентов с возможным нарушением отвода теплоты от топлива учитываются уже на стадии проекта современного энергоблока, для каждого инцидента предусматриваются меры и средства, гарантирующие непопадание радиоактивности в окружающую среду. Тем не менее возможность нарушения отвода теплоты от топлива объективна и вызывает обоснованные опасения специалистов и общественности в достаточности уровня безопасности атомной электростанции, что является постоянным аргументом против выбора вариантов с ядерной компонентой энергетики.
Существенно повысить надежность при нарушении отвода теплоты от топлива можно при использовании топлива в форме микротвэлов, непосредственно охлаждаемых водяным теплоносителем. Микротвэл представляет собой шарик, состоящий из топливного ядра, покрытого защитной оболочкой. Оболочки разделяют топливо и теплоноситель и обеспечивают удержание продуктов деления ядер. Микротвэлы, свободно насыпанные в корзины, образуют засыпки, которые являются структурными элементами ТВС. Теплоноситель входит сквозь стенку корзинки, движется сквозь поры между микротвэлами, нагревается и выходит через противоположную стенку корзинки. Температура топлива в микротвэле низкая, в центре зерна - порядка 350 °С. Это обеспечивает высокий уровень самозащищен-
ности, даже небольшое повышение температуры топлива выключает цепную реакцию без использования стержней управления и защиты. Тепловые переходные процессы в засыпках микротвэлов протекают быстро. Оболочки мик-ротвэлов не имеют швов, прочные и стойкие они надежно удерживают продукты деления, не содержат металлического циркония, могут работать при больших степенях выгорания топлива до 10-15 % тяжелых ядер, выдерживают нагрев до 1600 °С. Засыпка характеризуется большой площадью поверхности теплообмена и низкой мощностью теплового потока. Благоприятные теплофизические свойства системы микротвэлов и водяного теплоносителя позволяют избежать аварийных последствий перегрева топлива.
Известны предложения по использованию микротвэлов в ВВЭР [1], прямоточном реакторе с перегревом пара [2], реакторе со спектральным регулированием за счет изменения паросодержания [3] и в других [4-6]. Двухфазные течения в энергетическом оборудовании могут быть неустойчивыми [7], в частности известны примеры неустойчивых течений через пористые среды и зернистые засыпки [8]. Эффект может проявляться как автоколебание и/или расслоение теплоносителя по паросодержанию. С появлением уравнений движения, адаптированных для приложения к течениям двухфазной жидкости в засыпке микротвэлов при параметрах, характерных для атомной энергетики [9], возникла возможность расчетного определения неравномерностей распределения и неустойчивых режимов движения теплоносителя. Такая информация позволяет избежать нежелательного режима или предусмотреть меры для ослабления его проявления. Рассмотрим некоторые типы неустойчивости.
Для устойчивости состояния равномерного распределения пара по параллельным поровым каналам увеличение перепада давления Ар в зависимости от массового расходного паросодержания х должно происходить быстрее, чем по линейному закону (вторая производная перепада давления по расходному паросодержанию больше нуля). Тогда любая локальная разница расходов в соседних порах будет приводить к возрастанию местного сопротивления и появлению сил, стремящихся ее устранить. Режимы на рис. 1 получены расчетом по зависимостям, предложенным в [9], Ик - толщина засыпки. Варианты 1, 2, 4 устойчивы, режим 3 при х ~ 1 не устойчив. Неустойчивость уменьшается с увеличением массового расхода и размера частиц, со снижением абсолютного давления. Влияние пористости на эффект мало.
6
4
2
,0
Рис. 1. Зависимость перепада давления от паросодержания: 1 - (рм)/ = 55 кг/(м2 с), абсолютное давление р = 0,9 МПа; 2 - 55, 2,0; 3 - 55, 3,0; 4 - 22, 0,9; диаметр микротвэла << = 2 мм; пористость 8 = 0,37, «вода + пар», расчет по [9]
Режимы с более заметным проявлением эффекта подобраны на рис. 2. Течения, соответствующие линии 3 на рис. 1 и линиям 1, 3 на рис. 2, отличаются массовыми расходами. Снижение расхода приводит к усилению неустойчивости и возникновению многозначности в зависимости Дp(x). По мере снижения расхода на графике возникает точка перегиба при х = 1, затем она сдвигается в область более низких паросодержаний, параллельно формируется максимум перепада давления при x < 1 (х = 0,95-0,99), который также продвигается в область меньших х. Кривая 1 в совокупности с прямой 2, проведенной через точки Ap(0) и Ap(1) на рис. 2, ограничивают область неустойчивости (выделена заливкой).
0 0,1
0,3
0,5
0,7
0,9 х
Рис. 2. Зависимость перепада давления от паросодержания: 1, 2 - (рн>/ = 22 кг/(м2с); 3 - 11; p = 3,0 МПа; d = 2 мм; 8 = 0,37, «вода + пар», расчет по [9]
Наличие максимума давления при х ~ 1 должно приводить к появлению пор, заполненных только паром, и пор с водопаровой фазой при паросо-держании, меньшем среднего по засыпке. Следует ожидать, что в прямоточных реакторах с перегревом пара на выходе теплоноситель будет содержать долю влаги в капельной форме, если не предусмотреть меры по подавлению неустойчивости.
Оценим влияние фактора двухфазности теплоносителя на равномерность распределения потока в засыпке плоского реактора с радиальной раздачей. Схема устройства приведена на рис. 3. Теплоноситель при температуре насыщения входит через участок 1 в раздающий коллектор, поступает в насадки 7, 8, выходит в собирающий коллектор и покидает реактор через участок 6. В засыпке теплоноситель нагревается и обогащается паровой фазой. В собирающем коллекторе жидкость двухфазная, в раздающем - однофазная.
Р1
Н
Р4
¡111 ||Ш 1»
а ш ••••
Р2
р 5
ь
Р1
¡¡¡II ¡¡III ¡¡III ¡111
ИИ! .....) ":: !;!:
¡¡||:|::;|||
¡¡¡:| ¡¡Ж ¡¡Ж ¡Ж
Р6
Рис. 3. Схема реактора
2
3
И
1
6
4
5
Примем размер раздающего коллектора к, собирающего - И, длину реактора - ь. Отличие расходов по участкам 7, 8 определяет неравномерность распределения теплоносителя по засыпке. Пусть скорости по участкам отличаются незначительно: и = ^(1 + 88), 88 << 1.
В засыпке паросодержание растет линейно по длине, если энерговыделение равномерно по объему. По рис. 1 и 2 при х < (0,15-0,30) приближенно можно принять Лр(х) ~ (1 + кх), где к - коэффициент. Тогда перепад давления на слое засыпки представится выражением рки2(1 + кхь/2), где р - плотность жидкости; к - проницаемость засыпки; и - скорость фильтрации; хь - паросодержание на выходе из засыпки.
Динамические условия по участкам реактора (для простоты коллекторы считаем короткими и их гидравлическое сопротивление не учитываем)
определятся по формулам:
Р1 + ри:2 = р2 + ри22 = рз; рб + р иб2 = р5 + р и42 = р4; (1)
(Р1 + р)/2 - (р4 + Р5)/2 = ркит2( 1 + кхь/2); (2)
(рз + р)/2 - (рб + р5)/2 = рки82( 1 + кхь/2), (3)
где р - плотность двухфазной жидкости. Условия сохранения массы:
Ш12 = Ьи7; Ииб = р/р Шт. (4)
Из уравнений (1)-(4) следует оценка величины неравномерности 88:
88дф = [(р/ р )(Ь/И)2+(Ь/к)2]/[4к( 1 + кхь/2)]; (5)
88оф = [(Ь/И)2+(Ь/к)2]/4к; (6)
88дф/88оф = [(р/ р )(Ь/И)2+(Ь/к)2]/[((Ь/И)2+(Ь/к)2)(1 + кхь/2)], (7)
где 88дф, 88оф - различие скоростей по частям засыпки в случае двух- и однофазного теплоносителей соответственно.
Согласно (6), неравномерность распределения теплоносителя по засыпке в устройстве по рис. 3 неизбежна. Причина - во влиянии скоростного напора в коллекторах. Неравномерность возрастает с длиной коллекторов (реактора) и убывает с увеличением сопротивления засыпки. Подбором параметров устройства неравномерность может быть ограничена приемлемым уровнем. Несмотря на простоту уравнений (1)-(4), выражение (6) адекватно представляет наблюдаемые на практике эффекты [7, 10, 11].
Неравномерность при наличии фазового перехода в засыпке, согласно (5), дополнительно ослабляется за счет роста сопротивления насадки и увеличивается за счет уменьшения плотности потока в собирающем коллекторе. Итоговое изменение равномерности по сравнению с однофазным случаем дает формула (7).
Выражение (7) можно упростить, приняв гомогенную модель для теплоносителя в собирающем коллекторе р = рр"/((1 - хь)р"+ хьр), равенство
размеров коллекторов и пренебрегая величинами 1, (1 - хь)р" по сравне-
нию с хьр. Гомогенизации течения в коллекторе должен способствовать выброс из засыпки распределенных по длине парогазовых струек. С учетом упрощений из (7) получим
88дф/88оф = р/р"к. (8)
Поскольку при небольших хь можно записать к = ©р/р"; © < 1 (рис. 1, 2), то
88дф/88оф = 1/© > 1, (9)
т. е. в том же реакторе, в рамках принятых упрощений, равномерность распределения двухфазного теплоносителя по засыпке ухудшится по сравнению с однофазным теплоносителем.
Рассмотрим пример движения теплоносителя через засыпку с возможностью развития большой неравномерности расхода пара на выходе. Внутри засыпки однородно по объему зерен выделяется тепловая энергия. От зерен она передается жидкости. На вход подается теплоноситель (жидкость) с недогревом, внутри засыпки происходит рост его температуры до насыщения, далее идет движение с увеличением паросодержания по длине до выхода.
Изучим сначала случай равномерного движения. Воспользуемся упрощенным уравнением движения двухфазной среды в засыпке. Для случая равномерного по сечению паросодержания перепад давления Лр на засыпке выразится формулой
Лр = и-2 (1 + k Г xdz I; (10)
(1+k И
х = (х0 + Ь)2/м> - Ь; 20 = мЬ/(Ь + х0),
м - расходная скорость теплоносителя; 2 - координата по толщине засыпки; 20 - координата достижения теплоносителем температуры насыщения; х0 - паросодержание на выходе из засыпки; Ь - паросодержание на входе в засыпку (недогрев).
Тогда, принимая сечение засыпки за 1, получим, что расход пара через нее составляет х0, расходная скорость теплоносителя м (далее м положим равной 1), значение перепада давления
Лр = 1 + Ш2(х0 + Ь). (11)
Рассмотрим случай неравномерного движения. Разделим сечение засыпки на две половины. Предположим, что через первую течет теплоноситель со скоростью, превосходящей м настолько, что насыщение достигается только на выходе из засыпки. Скорость составит величину (х0 + Ь)/Ь, а перепад давления
Лр = [(х0 + Ь)/Ь]2. (12)
Если через вторую половину течет теплоноситель со скоростью, меньшей м и равной (Ь - х0)/Ь, то суммарный расход теплоносителя и пара по двум половинам окажется эквивалентным случаю равномерной фильтрации. Перепад давления по второй половине составит
Ар = [(Ь - хс)/Ь]2(1 + 2кЪх0/(Ъ - Х0)). (13)
Рассматривая величину к как параметр, определим его, приравнивая перепады давлений в (12) и (13) по половинам засыпки:
[(х0 + Ь)/Ь]2 = [(Ь - Х0)/Ь]2(1 + 2кЬх0/(Ь - Х0)),
откуда
к = 2/(Ь - х0). (14)
Подставив (14) в (11), найдем величины Ь и Х0, при которых перепад давления по формуле (12) больше перепада, вычисленного по (11)
1 + Х02/(Ь - Х0)(Х0 + Ь) > [(х0 + Ь)/Ь]2. (15)
Неравенство справедливо при р = хэ/Ь е [0,84, 1]. Для этих режимов смена равномерной фильтрации на сильно неравномерную будет сопровождаться снижением сопротивления засыпки. Таким образом, система в части режимов фильтрации оказывается неустойчивой против сильного возмущения.
В Ы В О Д
В прямоточных реакторах с паросодержанием около единицы или с перегретым паром на выходе теплоноситель может содержать долю влаги в капельной форме, если не предусмотреть меры по нейтрализации неустойчивости равномерного паросодержания. В реакторе радиального типа выбранной конструкции равномерность распределения двухфазного теплоносителя по засыпке ухудшится по сравнению с однофазным теплоносителем. Засыпка с двухфазным теплоносителем в части режимов фильтрации с недогревом на входе неустойчива против сильного возмущения течения.
Л И Т Е Р А Т У Р А
1. П е р с п е к т и в ы применения микротвэлов в ВВЭР/ Н. Н. Пономарев-Степной [и др.] // Атомная энергия. - 1999. - Т. 85, вып. 6. - С. 443-449.
2. Ф и л и п п о в, Г. А. Перспективы создания прямоточных микротвэльных ядерных реакторов с перегревом пара / Г. А. Филиппов, Р. Г. Богоявленский, А. А. Авдеев // Тяжелое машиностроение. - 2002. - № 1. - С. 7-11.
3. Я д е р н ы й реактор с водой под давлением с активной зоной на основе микротвэлов и способ осуществления его работы: пат. 2475869 РФ, МПК7 G21C1/00 / Е. И. Гришанин, П. Н. Алексеев; заявитель Федеральное государственное бюджетное учреждение «Национальный исследовательский центр "Курчатовский институт"». - № 2012105061/07; заявл. 15.02.2012; опубл. 20.02.2013 // Бюл. № 5.
4. И с а е в, А. Н. Реактор с водным теплоносителем малой мощности с топливом на базе шариковых микротвэлов / А. Н. Исаев // Атомная техника за рубежом. - 2007. - № 8. -С. 14-20.
5. И с а е в, А. Н. Применение шариковых микротвэлов реакторов с газовым теплоносителем в малых водо-водяных реакторах / А. Н. Исаев // Атомная техника за рубежом. -2007. - № 10. - С. 12-18.
6. О в о з м о ж н о с т и создания безопасных водо-водяных реакторов с тепловыделяющими сборками на основе шаровых твэлов / Е. Д. Домашев [и др.] // Промышленная теплотехника. - 2000. - Т. 22, № 4. - С. 51-55.
7. Р у к о в о д я щ и й документ по стандартизации Министерства тяжелого, энергетического и транспортного машиностроения СССР. Тепловой и гидравлический расчет тепло-обменного оборудования АЭС. Методические указания: РД 24.035.05-89. - Введ. 01.07.90. -Л.: НПО ЦКТИ, 1991. - 211 с.
8. K a v i a n y, M. Principles of Heat Transfer in Porous Media / M. Kaviany. - NY: Springer-Verlag, 1991. - 626 p.
9. С о р о к и н, В. В. Гидравлика и теплообмен шаровых засыпок в условиях активной зоны водо-водяных ядерных реакторов с микротвэлами / В. В. Сорокин. - Минск: Беларус. навука, 2010. - 191 с.
10. Б ы с т р о в, П. И. Гидродинамика коллекторных теплообменных аппаратов / П. И. Быстров, В. С. Михайлов. - М.: Энергоиздат, 1982. - 223 с.
11. Л е б е д е в, А. В. Влияние геометрии зернистого слоя на его гидродинамические и тепловые характеристики: дис. ... канд. техн. наук: 01.04.14 / А. В. Лебедев. - Новосибирск: ИТФ, 1988. - 156 с.
R E F E R E N C E S
1. P o n o m a r e v-S t e p n o y, N. N., Kukharkin, N. E., Khrulev, A. A., Degal'tsev, Iu. G., Glushkov, E. S., Filippov, G. A., Grishanin, E. I., & Fal'kovskii, L. N. (1999) Future Aspects of Coated Particles Fuel Application in LWR. Atomnaia Energiia [Nuclear Energy], 85 (6), 443-449 (in Russian).
2. F i l i p p o v, G. A., Bogoyavlensky, A. A., & Avdeev, A. A. (2002) Perspectives of Parallel-Flow Coated Particles Fuel of Nuclear Reactors Creation with Overheat of Steam. Tiazheloe Mashinostroenie [Heavy Engineering], 1, 7-11 (in Russian).
3. G r i s h a n i n, E. I., & Alekseev, P. N. (2013) Nuclear Reactor with Water Under Pressure with Active Zone on the Base of Coated Particles Fuel and the Manner of its Operation. Patent of the Russian Federation No 2475869 (in Russian).
4. I s s a e v, A. N. (2007) Reactor with Water Coolant of Small Capacity with Fuel on the Base of Spherical Coated fuel Particles. Atomnaia Tekhnika za Rubezhom [Nuclear Equipment Abroad], 8, 14-20 (in Russian).
5. I s s a e v, A. I. (2007) Application of Spherical Coated Fuel Particles of Reactors with Gas Heat Coolant in Small Two-Phase Reactors Atomnaia Tekhnika za Rubezhom [Nuclear Equipment Abroad], 10, 12-18 (in Russian).
6. D o m a s h e v, E. D., Domashev, V. E., Kolykhan, L. I., Fal'kovskii, L. N., & Shraiber, A. A. (2000) About Possibilities of Creation Safety Two-Phase Reactors with Heat-Emitting Adjustment on the Base of Spherical Fuel Elements. Promyshlennaia Teplotekhnika [Industrial Heat Engineering] 22 (4), 51-55 (in Russian).
7. RD 24.035.05-89. Guiding Document on Standardization of Ministry of Heavy, Energetic and Transport Engineering SSSR. Heat and Hydraulic Calculation of Heat Exchange Equipment Nuclear Power Station. Methodology Instructions. Leningrad: Scientific-Industrial Association of Boiler-and-Turbine Institute, 1991. 211 p. (in Russian).
8. K a v i a n y, M. (1991) Principles of Heat Transfer in Porous Media. NY, Springer-Verlag. 626 p.
9. S o r o k I n, V. V. (2010) Hydraulics and Heat Exchange of Spherical Fillings Under Conditions of Active Zone of Water-Cooled Nuclear Reactors with Coated Fuel Particles. Minsk: Belorussian Science. 191 p. (in Russian).
10. B y s t r o v, P. I., & Mikhailov, V. S. (1982) Hydrodynamics of Collector Heat Exchange Devices. Мoscow, Energoizdat. 223 p. (in Russian).
11. L e b e d e v, A. V. (1988) Vliianie Geometrii Zernistogo Sloia na Ego Gidrodina-micheskie i Teplovye Kharakteristiki. Dis. kand. tekhn. nauk [Influence of Geometry of Granular Bed on its Hydrodynamic and Heat Characteristics. Dissertation of Candidate of Technical Sciences]. Novosibirsk, ITF. 156 p. (in Russian).
Представлена кафедрой ТЭС Поступила 14.04.2014
