CC BY
74
АТОМНОЕ РЕАКТОРОСТРОЕНИЕ, МАШИНЫ, АГРЕГАТЫ И ТЕХНОЛОГИЯ
МАТЕРИАЛОВ АТОМНОЙ ПРОМЫШЛЕННОСТИ
УДК 621.039.517
СТРУКТУРА ТУРБУЛЕНТНОГО ПОТОКА В ПУЧКАХ СТЕРЖНЕЙ ТЕПЛОВЫДЕЛЯЮЩИХ СБОРОК РЕАКТОРНОЙ УСТАНОВКИ БРЕСТ-ОД-300
Д.В. Фомичев1, В.И. Солонин2
1АО “НИКИЭТ”, Москва, Российская Федерация e-mail: [email protected]
2МГТУ им. Н.Э. Баумана, Москва, Российская Федерация e-mail: [email protected]
Приведены результаты численного моделирования течения в 37-стержневых макетах тепловыделяющих сборок реакторной установки БРЕСТ-ОД-300, ди-станционированные двухзаходной трапециевидной навивкой по типу “ребро по ребру”. Показаны распределения избыточного статического давления, составляющих скорости и турбулентных характеристик потока в объеме пучков стержней макетов тепловыделяющих сборок. Проведен анализ структуры турбулентного течения в объеме пучков стержней.
Ключевые слова: тепловыделяющая сборка, структура течения, вычислительная гидродинамика.
STRUCTURE OF TURBULENT FLOW IN THE PIN BUNDLES OF THE BREST-OD-300
D.V. Fomichev1, V.I. Solonin2
1JSC “NIKIET”, Moscow, Russian Federation, e-mail: [email protected]
2Bauman State Technical University, Moscow, Russian Federation e-mail: [email protected]
The paper presents the results of the numerical simulation of 37-pins prototyping fuel assemblies of the reactor plant BREST-OD-300, which are spaced by a bifilar-helix wire wrapper of the “wire to wire” type. It shows the differential static pressure distributions as well as velocity components and the turbulent flow structure in a volume ofpin bundles of the prototyping fuel assemblies.The turbulent flow structure in a volume of pin bundles is analyzed.
Keywords: wire-spaces fuel assembly, turbulent flow structure, computational fluid dynamics.
Структура турбулентного потока в пучках стержней, дистанциони-рованных навитыми на стержни ребрами, определяется формируемой оребрением геометрией трактов течения теплоносителя. Тракты течения повторяют форму с периодом, равным доле продольного шага 4
4 ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение”. 2015. №3
навивки ребер, обратно пропорционально числу ребер на стержне. В области квазистационарного течения (после диссипации турбулентности, сформированной на входе в пучок стержней) турбулентное течение характеризуется периодическим изменением продольной и поперечной компонент скорости, так же как и турбулентных характеристик потока. Последние по сравнению с периодичностью средних компонент скорости в силу конвективного переноса могут иметь периодичность вдоль продольной координаты со сдвигом.
В настоящей работе с использованием современного CFD-кода выполнен анализ структуры течения в пучке стержней с двухзаходной трапециевидной формой оребрения, что может представлять интерес для создания тепловыделяющих сборок (ТВС) перспективных реакторов на быстрых нейтронах с тяжелым жидкометаллическим теплоносителем.
Анализ выполнен для одного значения шага навивки ребер и двух значений высот ребер применительно к исследованию макетов ТВС центральной и периферийной подзон (ЦЗ и ПЗ) активной зоны реакторной установки БРЕСТ-ОД-300 [1,2] при продувке макетов воздухом при числах Рейнольдса 5,02-104 и 5,82-104, соответствующих ранее полученным экспериментальным данным [3].
Объекты исследования. Математическая модель. Объектами исследования являются 37-стержневые макеты ТВС ЦЗ и ПЗ. Стержни расположены в правильной треугольной решетке с шагом 27,57 мм и размещены в шестигранном чехле с внутренним размером “под ключ” 173 мм. Основные геометрические характеристики стержней макетов приведены в таблице.
Расчетные сетки моделей макетов включают в себя проточную часть пучков, а также подводящий тракт конфузорной формы, формирующий поток на входе в макет ТВС. Подводящий тракт состоит из центрального конфузора с диаметром на входе в пучок стержней, равным 56 мм, и кольцевого периферийного конфузора с аналогичным диаметром 176 мм. В центральном и периферийном конфузорах установлены выравнивающие решетки, в том числе и общая для кон-фузоров решетка на входе в макет ТВС. Общая решетка выполнена с отверстиями диаметром 12 мм, центры которых соответствуют расположениям центров гидравлических ячеек пучков стержней ТВС ЦЗ и ПЗ.
Численное моделирование выполнено как для моделей ТВС с подводящим конфузорным трактом, так и без него.
Расчетная сетка моделей — полиэдральная, построенная с использованием пакета ANSYS Meshing [4] и последующей адаптацией в пакете ANSYS Fluent 15.0 [5]. Общее число элементов расчетной сетки равно ~ 108. *
ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение” 2015. № 3 5
Геометрические характеристики макетов ТВС
Параметр Обозначение Значение
Макет ТВС ЦЗ Макет ТВС ПЗ
Наружный диаметр стержня-имитатора твэла, мм d 20,81 22,53
Нижнее основание трапециевидной навивки, мм hi 4,29 3,65
Верхнее основание трапециевидной навивки, мм h-2 2,15
Высота трапециевидной навивки, мм ho 3,38 2,52
Описанный диаметр стержня-имитатора твэла, мм D 27,47
Шаг расположения стержней-имитаторов твэлов, мм S 27,57
Шаг дистанционирующей навивки, мм T 429
Наибольший зазор между двумя смежными навивками дистанци-онирующей проволоки, мм 5i 0,1
Наибольший зазор между чехлом и навивкой дистанционирующей проволоки, мм 52 1,0
Длина стержня-имитатора твэла, мм L 1040
Относительный шаг расположения стержней-имитаторов твэлов S/d 1,33 1,23
Относительный шаг трапециевидной навивки T/d 20,6 19,0
Гидравлический диаметр (проходного сечения пучка стержней без учета чехла), мм dI- 17,37 13,84
Гидравлический диаметр модели (проходного сечения пучка стержней с учетом чехла), мм d ^г.м 14,34 12,28
Для анализа влияния истечения потока воздуха из моделей ТВС на уровень давления в моделях расчетная сетка была продолжена за выходные торцы стержней на расстояние 80 мм (дополнительно анализировалось 180 мм).
Граничные условия на входе для моделей, учитывающих подводящий участок экспериментального стенда, состояли в задании расходов центрального и периферийного потоков воздуха при соответствующих 6
6 ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение”. 2015. №3
температурах, а для моделей без подводящего участка — в задании равномерного распределения расхода воздушного потока по сечению макетов. На выходе из моделей задано равномерное распределение статического давления. Твердые поверхности расчетных моделей предполагаются адиабатическими с выполнением условия прилипания.
Численное моделирование проводилось с помощью программного комплекса ANSYS Fluent 15.0 [5] с использованием стандартной (к-е)-модели турбулентности [6].
Уравнения переноса без учета влияния плавучести и дополнительных источников генерации и порождения диссипации турбулентности записываются в следующем виде:
дрк
~3t
д . , . + дх (рк“’:)
д
дх,
pt
р +---
Ок
дк
дхз
+ Gk — ре
Y к
дре д д
XX + дХ (ре“’:) = дХ,
рЛ де
P + X дХ
3 А
е е2
+ Cu-Gk — С2£р~, кк
где G
Yk =
к ри^и ,
дЩ
дхг
— генерация турбулентной кинетической энергии;
2реМ2 — фактор, учитывающий условия сжимаемости течения
потока воздуха как идеального газа; Mt = W — — турбулентное число
_____ V а2
Маха; а = \/yRT — скорость звука в потоке воздуха.
Турбулентная динамическая вязкость определяется по следующему соотношению:
С к2
Pt = рС^ —.
е
Для данной модели использован стандартный набор констант: Си = 1,44, CW = 1,92, С= 0,09, ок = 1,0, о£ = 1,3.
Свойства воздуха описывались законом идеального газа.
Сходимость полученного численного решения проверялась по равенству расходов на входе и выходе из пучка стержней. Из условия сходимости вычислений для достижения погрешности определения расходов не более 0,1 % требовалось проведение не менее 2000 итераций.
Полученные при численном моделировании результаты приведены для двух характерных областей течения: области регулярного расположения стержней (центральный стержень) и причехловой области (соседний к угловому стержень). Для каждого из стержней на радиусе, соответствующем полувысоте дистанционирующего ребра, выбрано 5 точек диагностики, размещенных относительно друг друга с шагом 30°.
Рассматриваемые характерные области течения и точки диагностики показаны на рис. 1.
ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение” 2015. № 3 7
Рис. 1. Пучок стержней и расположение точек диагностики характеристик потока:
а — пучок стержней (чехол не показан); б — расположение точек диагностики
Точка 1 расположена со стороны ребра, отклоняющего осевой поток в тангенциальном направлении (наветренная сторона), точка 5 — со стороны ребра, освобождающего сечение для потока (подветренная сторона).
Положение точек фиксировано по отношению к дистанциониру-ющим ребрам и изменяется в пространстве при движении потока от входа в пучок стержней к выходу. Отметим, что на расстоянии вдоль оси моделей, равном 1/6 шага навивки ребер, геометрия регулярных ячеек, в которых располагаются контрольные точки, повторяется. В связи с этим можно ожидать периодичности характеристик потока воздушного теплоносителя через 1/6 шага навивки ребер.
Распределение продольной компоненты скорости (близкой по значению к модулю вектора скорости) вдоль потока приведено на рис. 2. При деформировании потока подводящим устройством (рис. 2, а) после входного участка протяженностью ~550 мм (~1,3 шага навивки) значения средних скоростей в точках 1... 5 становятся квазистационарными. Наименьшее значение скорости потока соответствует точке 5, расположенной на подветренной стороне ребра, наибольшее значение скорости — точке 2. Скорости в точках 1 и 3 близки к скорости в точке 2. Скорость в точке 4 занимает промежуточное значение. Таким образом, в области стабилизированного течения осевой поток движется вдоль поверхности стержня между ребрами с продольной скоростью, которая с наветренной стороны на ~25 % больше, чем с подветренной. Значения скорости периодически (через 1/6 шага ди-станционирующей навивки) изменяются по близкой к гармонической
8 ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение”. 2015. №3
Рис. 2. Распределение продольной компоненты скорости в макете ТВС ЦЗ для входа потока в макет из подводящего участка (а) и для равномерного распределения скорости на входе в макет (б)
зависимости с амплитудой 1... 4,5 % среднего значения максимальной для подветренной области потока.
В пределах начального участка значения продольной скорости в выделенных точках зависят от их расположения по отношению к отверстиям входной решетки, перестройки струй из отверстий решетки в потоки в трактах пучка стержней. Периодическая компонента в изменении скорости появляется на расстоянии ^0,7 шага навивки. Значения скорости потока на начальном участке в среднем больше, чем в области квазистабилизированного течения.
Для равномерного распределения скорости на входе в макет (см. рис. 3) протяженность входного участка меньше, чем для варианта с подводящим устройством и составляет ^380 мм (^0,9 шага навивки). Скорость в области квазистационарного течения максимальна в точках 1, 2, 3 и минимальна в точках 4 и 5.
Амплитуда гармонической составляющей пульсации средней скорости составляет 1. . . 5 %. Скорость потока с наветренной стороны на ^25 % больше, чем с подветренной. Изменение скорости в пределах входного участка протекает более плавно, чем в варианте течения за входной решеткой (см. рис. 2, а).
В области стабилизированного течения для двух рассмотренных условий входа потока в макет ТВС наблюдается сдвиг на полпериода
ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение” 2015. №3 9
(1/12 шага дистанционирующей навивки) максимумов скорости потока: для тех продольных координат, в которых скорость в точках 1, 3, 5 имеет максимум, скорость в точках 2 и 4 имеет минимум (и наоборот). Таким образом, изменение продольной скорости потока является периодической функцией не только продольной, но и поперечной координаты.
На рис. 3 и 4 приведены характеристики турбулентности (кинетическая энергия и скорость диссипации) в точках 1... 5 в зависимости от продольной координаты. В пределах длины макета ТВС не наблюдается полной стабилизации характеристик турбулентности (в области ядра турбулентного течения у + ^ 30 для точек 1... 5). Протяженность участка с интенсивным снижением уровня турбулентных пульсаций скорости и уменьшением масштаба турбулентности составляет около
1... 1,5 шага дистанционирующей навивки. 10
Рис. 3. Кинетическая энергия турбулентности в области регулярных ячеек макета ТВС ЦЗ для входа потока в макет из подводящего участка
Рис. 4. Диссипация энергии турбулентности в области регулярных ячеек макета ТВС ЦЗ для входа потока в макет из подводящего участка
10 ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение”. 2015. №3
Рис. 5. Отношение тангенциальной и продольной компонент скорости потока в макете ТВС ЦЗ для входа потока в макет из подводящего участка (а) и для равномерного распределения скорости на входе в макет (б)
На рис. 5 показаны значения отношения тангенциальной и продольной компонент скорости потока в макете ТВС ЦЗ как для входа потока в макет из подводящего участка, так и при равномерном распределении скорости потока на входе в макет.
Векторные поля поперечных компонент скорости в области регулярной решетки стержней и в причехловой области показаны на рис. 6.
Влияние ребер навивки на поток проявляется в генерации вторичного течения, имеющего общее направление вдоль ребер, что обеспечивает переток теплоносителя из одного канала в соседний.
В областях, ограниченных смежными ребрами и поверхностью соседних стержней, имеет место вторичное вихревое течение, которое перемещается вслед за положением дистанционирующих ребер. В областях, где смежные ребра перекрывают переток теплоносителя между соседними ячейками, тангенциальная скорость обращается в ноль. Основной обмен потоков в соседних ячейках осуществляется в областях между поверхностями стержней при максимальном удалении дистанционирующих ребер относительно друг друга. *
ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение” 2015. №3 11
Рис. 6. Векторное поле поперечных скоростей потока в области регулярных (а) и причехловых (б) ячеек (521 мм)
В окрестностях контакта ребер соседних стержней происходит локальное ускорение потока, при этом образуется зона пониженного давления. 12
12 ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение”. 2015. №3
Рис. 7. Распределение статического давления в макете ТВС ЦЗ для входа потока в макет из подводящего участка
В периферийных ячейках пучков стержней поле скорости полностью определяется закруткой потока дистанционирующими ребрами, а также трением теплоносителя о стенку чехла. Кроме того, в периферийной области пучка стержней образуются “открытые” сечения, когда ребра расположены перпендикулярно стенке чехла. На этих участках продольный расход модельного теплоносителя в 1,4 раза превышает расход в ячейках регулярной области. При этом тангенциальная составляющая скорости в указанных зонах стремится к нулю.
Распределение избыточного статического давления в диагностируемых точках 1... 5 для центрального стержня приведено на рис. 7.
Видно, что падение статического давления в пучке стержней содержит периодическую составляющую с шагом 1/6 шага навивки ребер с амплитудой изменения ~100 Па. По представленной картине изменения статического давления начальный участок стабилизации течения может быть определен как ~0,25 шага навивки.
Особенностью распределения статического давления по высоте пучка стержней является наличие на начальном участке зоны пониженного давления в окрестности соприкосновения дистанционирую-щих ребер соседних стержней или ребра и чехла макета ТВС. Зона пониженного давления перестает наблюдаться через 1... 1,5 шага навивки, что может быть связано со стабилизацией течения в объеме пучка стержней. Отмеченная особенность характерна как для моделей с подводящим участком, так и без него.
Заключение. 1. Осредненные по Рейнольдсу значения продольной и поперечных (тангенциальных) компонент скорости достигают квазистационарного значения на расстоянии от входа в макет ТВС порядка 1,3... 0,9 шага навивки ребер в зависимости от степени возмущения потока на входе в макет. Наименьшие протяженности входного участка соответствуют равномерному распределению скорости потока *
ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение” 2015. №3 13
перед макетом. Область существенного изменения характеристик турбулентности составляет 1... 1,5 шага дистанционирующей навивки, однако квазистационарные значения характеристик турбулентности в пределах длины макета ТВС (~2,4 шага дистанционирующей навивки) не достигаются.
2. Осредненные по Рейнольдсу квазистационарные значения продольной скорости в области ядра турбулентного потока у поверхности стержня, расположенного в регулярной решетке, периодически изменяются с шагом, равным 1/6 шага дистанционирующей навивки в направлении вдоль потока. Максимумы продольной скорости в соседних точках вдоль периметра стержня имеют фазовый сдвиг, равный половине углового шага периодичности (180°).
3. Осредненные по Рейнольдсу тангенциальные компоненты скорости достигают стабилизированных значений на расстоянии от входа в макет ТВС порядка 0,9... 1,3 шага навивки, причем наибольшие значения соответствуют области подветренной части потока, где наряду с периодичностью скорости через 1/6 шага навивки ребер четко прослеживается периодичность 1/2 шага навивки ребер. Значения поперечных скоростей составляют 0,04-0,2 от значений продольных компонент, амплитуда их изменения — 0,02-0,05.
4. Векторные распределения тангенциальных компонент скорости показывают, что трапецеидальное оребрение создает вторичное течение вдоль направления завивки ребер, содержащее в своем составе вихревых зоны, что инициирует обмен между ячейками в пучках твэлов с подобной геометрией трактов теплоотвода.
5. Статическое давление в потоке, омывающем оребренные стержни, содержит периодическую составляющую, достигающую ~5 % среднего значения статического давления. Область начала квазилинейного изменения статического давления начинается на расстоянии ~0,25 шага навивки ребер.
6. Анализ полученных данных показывает, что использование стандартной (к-е)-модели турбулентности для расчетов характеристик течения потока в сборках с оребренными стержнями позволяет определить основные особенности структуры течения.
7. Данные анализа указывают на то, что стабилизация течения потока (потерь давления) наступает существенно раньше, чем стабилизация распределений продольной и поперечной компонент скорости. Стабилизация характеристик турбулентности требует протяженности области течения более трех шагов навивки ребер. Эти данные соответствуют результатам анализа развития турбулентного течения в пристенном пограничном слое.
При необходимости уменьшить длину участка стабилизации в сборках стержней с оребрением необходимо уменьшать турбулентные возмущения потока на входе. 14
14 ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение”. 2015. №3
ЛИТЕРАТУРА
1. Белая книга ядерной энергетики / В.В. Орлов, М.М. Селиверстов, В.А. Тищенко, В.В. Ужанова, В.С. Смирнов, И.Х. Ганев, А.В. Лопаткин, С.В. Брюнин,
A. Н. Кархов, С.В. Европин, Г.Е. Шаталов, В.Е. Сытников, П.И. Долгошеев,
B. Б. Козлов, В.П. Фотин / под ред. Е.О. Адамова. М.: Изд-во ГУП НИКИЭТ, 1998. 356 с.
2. Драгунов Ю. Г , Лемехов В. В., Смирнов В. С., Чернецов Н. Г. Технические решения и этапы разработки реакторной установки БРЕСТ-ОД-300 // Атомная энергия. 2012. Т. 113. № 1. С. 58-64.
3. Fomichev D.V., Solonin V.I. A CFD Analysis of Hydraulic Characteristics of the Rod Bundles in the BREST-OD-300 Wire-Spaced Fuel Assemblies // International Conference on Nuclear Energy Science and Technology, Zurich, Switzerland, 2014. Vol. 8. No. 7. Part XII. Р 1514-1519.
4. ANSYS Meshing User’s Guide, Rel. 14.5. ANSYS Inc., 2012.
5. ANSYS FLUENT, Theory Guide, Rel. 14.5. ANSYS Inc., 2012.
6. Launder B.E., Spalding D.B. Lectures in Mathematical Models of Turbulence. London, England: Academic Press, 1972.
REFERENCES
[1] Orlov V.V., Seliverstov M.M., Tishchenko V.A., Zhanova V.U., Smirnov V.S., Ganev I.H., Lopatkin A.V., Bryunin S.V., Karhu A.N., Evropin S.V., Shatalov G.E., Sitnikov V.E., Dolgosheev PI. Kozlov V.B., Photinus V.P., Adamov E.O., ed. The White Book of Nuclear Energy. Moscow, NIKIET Publ, 1998. 356 p.
[2] Dragunov Yu.G., Lemekhov V.V., Smirnov V.S., Chernetsov N.G. The Technical Solutions and Development Phases of the Reactor Plant BREST-OD-300. Atomnaya energiya [Atomic Energy J.], 2012, vol. 113, no. 1, pp. 58-64 (in Russ.).
[3] Fomichev D.V., Solonin VI. A CFD Analysis of Hydraulic Characteristics of the Rod Bundles in the BREST-OD-300 Wire-Spaced Fuel Assemblies. International Conference on Nuclear Energy Science and Technology, Zurich, Switzerland, 2014, vol. 8, no. 7, part XII, pp. 1514-1519.
[4] ANSYS Meshing User’s Guide, Rel. 14.5. ANSYS Inc., 2012.
[5] ANSYS FLUENT, Theory Guide, Rel. 14.5. ANSYS Inc., 2012.
[6] Launder B.E., Spalding D.B. Lectures in Mathematical Models of Turbulence. London, England, Academic Press, 1972.
Статья поступила в редакцию 26.11.2014
Фомичев Дмитрий Вадимович — аспирант кафедры “Ядерные реакторы и установки” МГТУ им. Н.Э. Баумана. Сотрудник отдела “Теплофизика” ОАО “НИКИЭТ”. Автор трех научных работ в области ядерной энергетики.
АО “НИКИЭТ”, Российская Федерация, 107140, Москва, Малая Красносельская ул., д. 2/8.
Fomichev D.V. — postgraduate student of the “Nuclear Reactors and Plants” Department of Bauman Moscow State Technical University. Employee of the “Thermal Physics” Department of JSC “NIKIET”. Author of 3 publications in the field of nuclear technologies.
JSC “NIKIET”, Malaya Krasnosel’skaya ul. 2/8, Moscow, 107140 Russian Federation. *
ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение” 2015. № 3 15
Солонин Владимир Иванович — д-р техн. наук, профессор кафедры “Ядерные реакторы и установки” МГТУ им. Н.Э. Баумана. Автор более 100 научных работ в области ядерных технологий.
МГТУ им. Н.Э. Баумана, Российская Федерация, 105005, Москва, 2-я Бауманская ул., д. 5.
Solonin VI. - Dr. Sci. (Eng.), professor of the “Nuclear Reactors and Plants” Department of Bauman Moscow State Technical University. Author of more than 100 publications in the field of nuclear technologies.
Bauman Moscow State Technical University, 2-ya Baumanskaya ul. 5, Moscow, 105005 Russian Federation.
Просьба ссылаться на эту статью следующим образом:
Фомичев Д.В., Солонин В.И. Структура турбулентного потока в пучках стержней тепловыделяющих сборок реакторной установки БРЕСТ-ОД-300 // Вестник МГТУ им. Н.Э. Баумана. Сер. Машиностроение. 2015. № 3. C. 4-16.
Please cite this article in English as:
Fomichev D.V., Solonin V.I. Structure of Turbulent Flow in the Pin Bundles of the BREST-OD-300. Vestn. Mosk. Gos. Tekh. Univ. im. N.E. Baumana [Herald of the Bauman Moscow State Tech. Univ., Mech. Eng.], 2015, no. 3, pp. 4-16. 16
16 ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение”. 2015. №3
