Теплоэнергетика, 2021, № 3, стр. 36-44

Исследование гидродинамики закрученных потоков в области направляющего канала за перемешивающими дистанционирующими решетками ТВС реактора типа PWR

С. М. Дмитриев a, А. В. Герасимов a*, А. А. Добров a, Д. В. Доронков a, А. П. Живодеров a, Е. В. Рубцова a, А. В. Рязанов a, Д. Н. Солнцев a, А. Н. Пронин a, А. Е. Хробостов a

a Нижегородский государственный технический университет им. Р.Е. Алексеева
603155 г. Нижний Новгород, ул. Минина, д. 24, Россия

* E-mail: poorvenom@yandex.ru

Поступила в редакцию 08.08.2019
После доработки 18.06.2020
Принята к публикации 24.06.2020

Полный текст (PDF)

Аннотация

Приведены результаты исследований течения теплоносителя в тепловыделяющих сборках за перемешивающими дистанционирующими решетками реактора типа PWR. Цель работы – оценка эффективности применения перемешивающих решеток различных конструкций. Для достижения поставленной цели были проведены экспериментальные исследования на аэродинамическом стенде с масштабными моделями фрагментов топливных кассет с перемешивающими дистанционирующими решетками разных конструкций. В качестве показательной области исследований были выбраны две смежные ячейки области направляющего канала, конструктивной особенностью которых является различная пространственная ориентация турбулизаторов, установленных на решетках-интенсификаторах. Общая картина течения представлена векторными полями тангенциальных скоростей, а также распределением относительных скоростей в зазорах между твэлами и направляющим каналом. Для оценки эффективности применения конструкций решеток относительно перемешивания потока теплоносителя был проведен анализ коэффициентов внутриячейкового вихреобразования и межъячейкового перемешивания. В результате этого анализа показано, что в альтернативной конструкции изменение направления пары дефлекторов ячейки позволило увеличить коэффициент межъячейкового обмена в 1.13 раза относительно базовой конструкции, а коэффициент внутриячейкового вихреобразования – в 2.2 раза. Следовательно, использование решетки альтернативной конструкции в целях лучшего перемешивания потока теплоносителя является более предпочтительным. Накопленная база данных по течению теплоносителя в ТВС-Квадрат служит для инженерного обоснования конструкций активных зон реакторов PWR. Результаты экспериментальных исследований используются для верификации CFD-кодов как зарубежной, так и отечественной разработки, а также программ детального поячеечного расчета активных зон для уменьшения консерватизма при обосновании их теплотехнической надежности .

Ключевые слова: перемешивающие решетки, дефлекторы, направляющий канал, аэродинамический стенд, внутриячейковое вихреобразование, межъячейковое перемешивание

С учетом результатов исследовательских работ, проведенных в целях обоснования работоспособности тепловыделяющих сборок (ТВС), а также опыта конструирования и эксплуатации топливных сборок [1] были предложены перемешивающие дистанционирующие решетки (ПДР) двух типов. Особенность конструкций данных решеток заключается в установке дополнительных турбулизаторов-дефлекторов, имеющих различную форму и пространственную ориентацию, которые создают интенсивную вихревую структуру внутри ячейки и поперечное перетекание теплоносителя между соседними ячейками, расположенными в области направляющего канала (НК). Учет этих процессов очень важен при конструировании дистанционирующих и перемешивающих элементов, а также при оценке эффективности использования данных конструкций для перемешивания потока теплоносителя [2]. Перемешивание потока способствует интенсификации процесса конвективного теплообмена на поверхности тепловыделяющих элементов и увеличивает поперечное межъячейковое перетекание, приводящее к выравниванию неоднородностей поля температур в поперечном сечении тепловой сборки, что способствует увеличению запаса до кризиса теплоотдачи.

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

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

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

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

ЭКСПЕРИМЕНТАЛЬНЫЙ СТЕНД И МЕТОДИКА ИССЛЕДОВАНИЙ

Комплекс экспериментальных исследований, направленных на изучение течения потока в ТВС, выполнялся на аэродинамическом стенде, включающем в себя экспериментальную модель – канал квадратного сечения, в котором размещена изучаемая перемешивающая дистанционирующая решетка [1320]. Использование воздуха в качестве рабочей среды обусловлено необходимостью выхода режима движения потока в экспериментальной модели в область автомодельного течения для соблюдения условий теории гидродинамического подобия [3]. При проведении экспериментальных исследований учитывалось геометрическое и кинематическое подобие процессов.

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

По длине экспериментальной модели за ПДР выбраны характерные сечения, в которых проводились измерения (рис. 1) для изучения течения потока. Для измерения полного вектора скорости в показательной области направляющего канала использовался специальный пневмометрический зонд [21].

Рис. 1.

Общий вид пучка твэлов модели с установленной ПДР и схема расположения сечений измерения гидродинамических характеристик потока. 1 – твэл; 2 – направляющий канал; 3 – поперечные плоскости; 4 – дефлекторы; 5 – ПДР

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

Рис. 2.

Поперечное сечение экспериментальной модели области ячеек НК базовой конструкции (а), области ячеек НК альтернативной конструкции (б) и общий вид модели ПДР (в)

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

Рис. 3.

Профили исследуемых дефлекторов в ячейках области НК. а – дефлекторы, применяемые в ячейках Type 1 базовой и альтернативной конструкций ПДР; б – дефлекторы, применяемые в ячейках Type 2 базовой конструкций ПДР; в – вид сбоку

в альтернативной конструкции ячейка Type 1 идентична аналогичной ячейке базовой конструкции, в ячейке Type 2 изменено направление отгиба дефлекторов на обратное и увеличена площадь перекрытия затесненного дефлектора профилем НК до аналогичного значения для дефлектора в ячейке Type 1 конструкций обоих типов (см. рис. 2, б).

Угол отгиба β = 25° (см. рис. 3, в) и угол наклона линии отгиба α = 15° (см. рис. 3, а) остаются одинаковыми для всех представленных конструкций.

Анализ интенсивности течения теплоносителя производится с использованием коэффициентов межъячейкового перетекания ${{F}_{m}}$ (интегральная интенсивность скоростей в зазорах) и внутриячейкового вихреобразования ${{F}_{{{v}or\_mom}}}$ (интегральный параметр закрутки потока для переменного поля скорости) [2]:

(1)
${{F}_{m}} = \frac{1}{{{{l}_{{rod}}}}}\int\limits_0^l {\frac{{{{w}_{p}}}}{{\bar {w}}}} \,\,dl;$
(2)
${{F}_{{{v}or\_mom}}} = \frac{{\int\limits_{ - r}^r {{{r}^{2}}{{V}_{{tan}}}{{U}_{z}}dr} }}{{R\int\limits_{ - r}^r {rU_{z}^{2}dr} }},$
где $l$ – расстояние между твэлами в зазоре; ${{l}_{{rod}}}$ – длина зазора между твэлами или направляющим каналом и твэлом, м/с; wp – составляющая вектора скорости в точке измерения составляющих вектора скорости потока, м/с; $\bar {w}$ – среднерасходная скорость, м/с; Vtan – тангенциальная скорость, м/с; Uz – аксиальная составляющая абсолютной скорости, м/с; r – расстояние от центра ячейки до измеряемой точки в зазоре, м; R – минимальное расстояние от центра ячейки до твэла или направляющего канала, м.

РЕЗУЛЬТАТЫ ЭКСПЕРИМЕНТАЛЬНЫХ ИССЛЕДОВАНИЙ

Исследования коэффициента гидравлического сопротивления (КГС) решетки-интенсификатора проводились в диапазоне чисел Re от 20 000 до 85 000. Анализ полученных результатов показал, что изменение ориентации турбулизаторов повлияло на КГС незначительно, значение КГС решетки изменилось менее чем на 1% и составило 1.25 при числе Re = 85 000.

Локальные гидродинамические характеристики потока измерялись при числах Re = 90 000. Данные исследования были направлены на определение проекций вектора абсолютной скорости в области НК, на основе которых строились векторные поля скоростей ${{W}_{x}}$ и ${{W}_{y}}$ в поперечных сечениях исследуемой модели на безразмерных расстояниях ${L \mathord{\left/ {\vphantom {L {{{d}_{{\text{г}}}}}}} \right. \kern-0em} {{{d}_{{\text{г}}}}}}$ ($L$ – расстояние от края перемешивающей дистанционирующей решетки до точки измерения гидродинамических характеристик потока, т.е. до кончика чувствительной части пневмометрического зонда; ${{d}_{{\text{г}}}}$ – гидравлический диаметр всей модели, т.е. квадратного канала с пучком твэлов), а также выявлялась графическая зависимость аксиальных и тангенциальных скоростей от этого безразмерного расстояния.

Векторные поля в зазоре решетки базовой конструкции для ячеек обоих типов показаны на рис. 4, а распределения относительных скоростей в этих ячейках – на рис. 5. Аналогичные характеристики для решетки альтернативной конструкции представлены на рис. 6 и 7.

Рис. 4.

Векторное поле скоростей за решеткой базовой конструкции. Расстояние ${L \mathord{\left/ {\vphantom {L {{{d}_{{\text{г}}}}}}} \right. \kern-0em} {{{d}_{{\text{г}}}}}}{\text{:}}$ а – 2.5; б – 22.5

Рис. 5.

Распределение относительных скоростей (${{{{W}_{i}}} \mathord{\left/ {\vphantom {{{{W}_{i}}} {\bar {W}}}} \right. \kern-0em} {\bar {W}}}$) в зазорах между твэлами за решеткой базовой конструкции ячейки Type 1 (а) и Type 2 (б) ($\bar {W}$ = 33.29 м/с). Зазор: 1 – левый; 2 – правый; 3 – верхний; 4 – нижний

Рис. 6.

Векторное поле скоростей альтернативной конструкции. Расстояние ${L \mathord{\left/ {\vphantom {L {{{d}_{{\text{г}}}}}}} \right. \kern-0em} {{{d}_{{\text{г}}}}}}{\text{:}}$ а –3.0; б – 13.6

Рис. 7.

Распределение относительной скорости (${{{{W}_{i}}} \mathord{\left/ {\vphantom {{{{W}_{i}}} {\bar {W}}}} \right. \kern-0em} {\bar {W}}}$) в зазорах ячейки альтернативной конструкции Type 1 (а) и Type 2 (б) ($\bar {W}$ = 33.29 м/с). Обозначения см. рис. 5

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

поток в ячейке Type 1 имеет интенсивную вихревую структуру, которая не рассеивается по всей исследуемой длине; вихрь изменяет свое положение в поперечной плоскости в сторону зазора, смежного с ячейкой Type 2, на расстоянии ${L \mathord{\left/ {\vphantom {L {{{d}_{{\text{г}}}}}}} \right. \kern-0em} {{{d}_{{\text{г}}}}}}$ = = 22.5 (см. рис. 4), в свою очередь в этой ячейке вихревая структура отсутствует;

в нижнем и левом зазорах ячейки Type 1 поперечное течение потока имеет интенсивный характер по всей длине исследуемой области и поперечные скорости в них составляют 0.2 и 0.3 относительно среднерасходной скорости соответственно (см. рис. 5). В ячейке Type 1 теплоноситель меняет свое направление на обратное через верхний зазор на расстоянии ${L \mathord{\left/ {\vphantom {L {{{d}_{{\text{г}}}}}}} \right. \kern-0em} {{{d}_{{\text{г}}}}}}$ = 15, что обусловлено смещением вихря в сторону правого зазора (см. рис. 4). Межъячейковое перетекание потока через зазор между ячейками Type 1 и Type 2 меняет свое направление на обратное на расстоянии ${L \mathord{\left/ {\vphantom {L {{{d}_{{\text{г}}}}}}} \right. \kern-0em} {{{d}_{{\text{г}}}}}}$ = 7.5;

в правом и нижнем зазорах ячейки Type 2 поперечное течение теплоносителя не изменяет своей направленности по всей исследуемой длине, тангенциальные скорости составляют 0.1 и 0.2 соответственно (см. рис. 5).

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

в обеих ячейках области НК образуется интенсивная вихревая структура (см. рис. 6). В ячейке Type 1 безразмерная длина, на которой поток сохраняет вихревую структуру течения, равна ${L \mathord{\left/ {\vphantom {L {{{d}_{{\text{г}}}}}}} \right. \kern-0em} {{{d}_{{\text{г}}}}}}$ = = 13.6 (см. рис. 6), а в ячейке Type 2 длина жизни вихря составляет ${L \mathord{\left/ {\vphantom {L {{{d}_{{\text{г}}}}}}} \right. \kern-0em} {{{d}_{{\text{г}}}}}}$ = 7.5, что свидетельствует о более симметричном образовании вихревого течения потока, чем в ячейке базовой конструкции (в аналогичной ячейке Type 2 вихрь отсутствовал);

в ячейке Type 1 в нижнем зазоре меняется направление течения потока на обратное на расстоянии ${L \mathord{\left/ {\vphantom {L {{{d}_{{\text{г}}}}}}} \right. \kern-0em} {{{d}_{{\text{г}}}}}}$ = 13.6 (рис. 6, б и 7, а), а в правом зазоре, смежном с ячейкой Type 2, направление течения меняется на расстоянии ${L \mathord{\left/ {\vphantom {L {{{d}_{{\text{г}}}}}}} \right. \kern-0em} {{{d}_{{\text{г}}}}}}$ = 7.5. В левом зазоре ячейки Type 1 межъячейковое перетекание происходит очень интенсивно, на расстоянии ${L \mathord{\left/ {\vphantom {L {{{d}_{{\text{г}}}}}}} \right. \kern-0em} {{{d}_{{\text{г}}}}}}$ = 25 за решеткой тангенциальная скорость составляет 0.2;

в ячейке Type 2 в верхнем зазоре межъячейковое перетекание потока теряет свою интенсивность на расстоянии ${L \mathord{\left/ {\vphantom {L {{{d}_{{\text{г}}}}}}} \right. \kern-0em} {{{d}_{{\text{г}}}}}}$ = 16.2. В верхнем зазоре ячейки Type 1 и в правом зазоре ячейки Type 2 межъячейковое перетекание довольно интенсивное; на расстоянии ${L \mathord{\left/ {\vphantom {L {{{d}_{{\text{г}}}}}}} \right. \kern-0em} {{{d}_{{\text{г}}}}}}$ = 25 за решеткой тангенциальная скорость составляет 0.2.

ОЦЕНКА ЭФФЕКТИВНОСТИ ИСПОЛЬЗОВАНИЯ КОНСТРУКЦИЙ ПДР

Оценка эффективности использования конструкций перемешивающих дистанционирующих решеток производится на основе анализа коэффициентов внутриячейкового вихреобразования ${{F}_{{{v}or\_mom}}}$ и межъячейкового перетекания ${{F}_{m}},$ которые усреднялись по зазорам и центральным диагоналям каждой характерной ячейки, далее строились графики зависимости данных коэффициентов относительно безразмерной координаты L/dг и проводилось их сравнение (рис. 8).

Рис. 8.

Распределение коэффициентов ${{F}_{{{v}or\_mom}}}$ (а) и ${{F}_{m}}$ (б). 1 – конструкция А, Type 1; 2 – конструкция А, Type 2; 3 – конструкция Б, Type 1; 4 – конструкция Б, Type 2

На основе результатов расчета коэффициентов ${{F}_{m}}$ и ${{F}_{{{v}or\_mom}}}$ можно заключить, что для решетки альтернативной конструкции в ячейках Type 1 и Type 2 коэффициент ${{F}_{m}}$ увеличился в 1.25 и 1.07 раза соответственно по сравнению с базовой конструкцией, а коэффициент внутриячейкового вихреобразования ${{F}_{{{v}or\_mom}}}$ увеличился в 1.4 и 3.0 раза для ячеек обоих типов по сравнению с базовой решеткой соответственно (см. таблицу).

Усредненные значения ${{F}_{m}}$ и ${{F}_{{{v}or\_mom}}}$ в решетках базовой и альтернативной конструкций

Показатель Ячейка Type 1 Ячейка Type 2
базовая конструкция (Б) альтернативная конструкция (А) базовая конструкция (Б) альтернативная конструкция (А)
Коэффициент межъячейкового перетекания ${{F}_{m}}$ 0.16 0.20 0.15 0.16
Отношение  $\frac{{{{F}_{m}}{\text{(А)}}}}{{{{F}_{m}}{\text{(Б)}}}}$ 1.25 1.07
Коэффициент внутриячейкового вихреобразования ${{F}_{{vor\_mom}}}$ 0.09 0.11 0.02 0.06
Отношение  $\frac{{{{F}_{{vor\_mom}}}(А)}}{{{{F}_{{vor\_mom}}}(Б)}}$ 1.4 3

ВЫВОДЫ

1. Анализ результатов экспериментальных исследований течения теплоносителя за решетками‑интенсификаторами показал, что изменение формы и взаимного расположения дефлекторов привело к симметричному и более интенсивному образованию вихревой структуры в ячейках обоих типов в области НК по сравнению с базовой конструкцией, а также позволило увеличить интенсивность межъячейкового перетекания, что приводит к увеличению среднего значения относительных тангенциальных скоростей в зазорах до 0.2.

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

Список литературы

  1. Основное оборудование АЭС с корпусными реакторами на тепловых нейтронах / С.М. Дмитриев, О.А. Бых, Ю.К. Панов, Н.М. Сорокин, Д.Л. Зверев, В.А. Фарафонов. М.: Машиностроение, 2013.

  2. Митрофанова О.В. Гидродинамика и теплообмен закрученных потоков в каналах ядерно-энергетических установок. М.: Физматлит, 2010.

  3. Гухман А.А. Введение в теорию подобия. М.: Высшая школа, 1973.

  4. Trupp A.C., Azad R.S. The structure of turbulent flow in triangular array rod bundles // Nucl. Eng. Des. 1975. V. 32. Is. 1. P. 47–84. https://doi.org/10.1016/0029-5493(75)90090-4

  5. Vonka V. Measurement of secondary flow vortices in a rod bundle // Nucl. Eng. Des. 1988. V. 106. Is. 2. P. 191–207. https://doi.org/10.1016/0029-5493(88)90277-4

  6. Eifler W., Nijsing R. Experimental investigation of velocity distribution and flow resistance in a triangular array of parallel rods // Nucl. Eng. Des. 1967. V. 5. Is. 1. P. 22–42. https://doi.org/10.1016/0029-5493(67)90075-1

  7. Rowe D.S., Johnson B.M., Knudsen J.G. Implications concerning rod bundle crossflow mixing based on measurements of turbulent flow structure // Int. J. Heat Mass Transfer. 1974. V. 17. Is. 3. P. 407–419. https://doi.org/10.1016/0017-9310(74)90012-X

  8. Caraghiaur D. Experimental study and modelling of spacer grid influence on flow in nuclear fuel assemblies: Technical report. KTH Stockholm, 2009.

  9. Carajilescov P. Experimental and analytical study of axial turbulent flows in an interior subchannel of a bare rod bundle: PhD thesis. Massachusetts Institute of Technology, 1976.

  10. Large eddy simulations of turbulent flow for grid-to-rod fretting in nuclear reactors / J. Bakosi, M.A. Christon, R.B. Lowrie, L.A. Pritchett-Sheats, R.R. Nourgaliev // Nucl. Eng. Des. 2013. V. 262. P. 544–561. https://doi.org/10.1016/j.nucengdes.2013.06.007

  11. Peña-Monferrer C., Muñoz-Cobo J.L., Chiva S. CFD turbulence study of PWR spacer-grids in a rod bundle // Sci. Technol. Nucl. Install. 2014. V. 2014. P. 15. https://www.hindawi.com/journals/stni/2014/635651/

  12. Sergeev D.A., Kandaurov A.A., Troitskaya Yu.I. The particular use of PIV methods for the modeling of heat and hydrophysical processes in the nuclear power plants // J. Phys.: Conf. Ser. 2017. V. 891. P. 012088. https://doi.org/10.1088/1742-6596/891/1/012088

  13. Гидродинамические особенности течения теплоносителя за перемешивающей дистанционирующей решеткой ТВС-КВАДРАТ реактора PWR / О.Б. Самойлов, А.С. Носков, Д.Л. Шипов, С.М. Дмитриев, А.А. Добров, Д.В. Доронков, М.А. Легчанов, А.Н. Пронин, Д.Н. Солнцев, В.Д. Сорокин, А.Е. Хробостов // Теплоэнергетика. 2019. № 4. С. 32–38. https://doi.org/10.1134/S0040363619040076

  14. Расчетно-экспериментальные исследования локальной гидродинамики и массообмена потока теплоносителя в ТВС-Квадрат реакторов PWR с перемешивающими решетками / С.М. Дмитриев, О.Б. Самойлов, А.Е. Хробостов, А.В. Варенцов, А.А. Добров, Д.В. Доронков, В.Д. Сорокин // Теплоэнергетика. 2014. № 8. С. 20–27. https://doi.org/10.1134/S0040363614080050

  15. Экспериментальные исследования локального массообмена и эффективности перемешивания теплоносителя дистанционирующими решетками в ТВС реактора КЛТ-40С / А.В. Варенцов, Д.В. Доронков, Е.С. Купричева, Д.Н. Солнцев, В.Д. Сорокин // Труды НГТУ им. Р.Е. Алексеева. 2012. № 1. С. 107–113.

  16. Исследования локальной гидродинамики и межъячеечного массообмена потока теплоносителя в районе направляющих каналов тепловыделяющих сборок реакторов PWR / С.М. Дмитриев, М.А. Легчанов, А.Е. Хробостов, А.В. Варенцов, Д.В. Доронков, А.А. Добров // Промышленная энергетика. 2013. № 12. С. 45–50.

  17. Особенности течения теплоносителя в ТВС-КВАДРАТ реактора PWR при постановке перемешивающих дистанционирующих решеток с различными типами дефлекторов / А.В. Варенцов, Д.В. Доронков, Е.М. Илютина, И.В. Каратушина, В.Д. Сорокин, А.Е. Хробостов // Труды НГТУ им. Р.Е. Алексеева. 2015. № 3. С. 134–143.

  18. Исследование массообмена теплоносителя за перемешивающими решетками ТВС реакторов ВБЭР-300 с целью обоснования их эффективности / С.М. Дмитриев, А.В. Варенцов, А.А. Добров, Д.В. Доронков, М.А. Легчанов, А.Е. Хробостов // Труды НГТУ им. Р.Е. Алексеева. 2013. № 5. С. 197–205.

  19. Исследования гидродинамических характеристик и особенностей течения теплоносителя за дистанционирующей решеткой тепловыделяющей сборки реактора плавучего энергоблока / С.М. Дмитриев, Д.В. Доронков, М.А. Легчанов, А.Н. Пронин, Д.Н. Солнцев, В.Д. Сорокин, А.Е. Хробостов // Теплофизика и аэромеханика. 2016. Т. 23. № 3. С. 385–394.

  20. Изучение гидродинамических процессов течения теплоносителя в ТВС-КВАДРАТ реактора PWR с различными перемешивающими дистанционирующими решетками / С.М. Дмитриев, А.А. Добров, Д.В. Доронков, А.Н. Пронин, Д.Н. Солнцев, В.Д. Сорокин, А.Е. Хробостов // Теплофизика и аэромеханика. 2018. Т. 25. № 5. С. 725–734.

  21. Применение многоканального пневмометрического зонда для исследования профиля скорости теплоносителя в моделях топливных кассет ядерных реакторов / С.М. Дмитриев, А.А. Добров, М.А. Легчанов, А.Е. Хробостов // Приборы и методы измерений. 2015. Т. 6. № 2. С. 188–195.

Дополнительные материалы отсутствуют.