CC BY
193
Наука и Образование
МГТУ им. Н.Э. Баумана
Сетевое научное издание
Наука и Образование. МГТУ им. Н.Э. Баумана. Электрон. журн. 2014. № 10. С. 398-412.
ISSN 1994-0448
DOI: 10.7463/1014.0727220
Представлена в редакцию: Исправлена:
© МГТУ им. Н.Э. Баумана
УДК 621.039.577-182.3
Моделирование теплообмена в змеевиковом теплообменнике применительно к реакторной установке "УНИТЕРМ"
Солонин В. И.1, Сатин А. А.1*
02.09.2014 21.09.2014
asatin9Qgmail.ru :МГТУ им. Н.Э. Баумана, Москва, Россия
В данной статье рассматривается моделирование процесса теплообмена в промежуточном теплообменнике змеевикового типа, входящего в состав контура теплоотвода реакторной установки "УНИТЕРМ", выполненное с применением гидрогазодинамического пакета ANSYS CFX. Проводится сравнение и анализ результатов, полученных с использованием моделей турбулентности k-e RNG, Shear Stress Transport (SST) на "грубой сетке". Показаны особенности течения потоков в межтрубном пространстве, их влияние на теплообмен и температур теплоносителей первого и промежуточных контуров. Представлены картины распределения теплогидравличеких параметров: температуры, перепада давления, скорости.
Ключевые слова: Реакторная установка "УНИТЕРМ", контур теплоотвода, активная зона, змеевиковый теплообменник, вычислительная гидродинамика и теплообмен
1 Введение
Проект РУ «УНИТЕРМ», разработанный «НИКИЭТ» [1-3], является эффективной основой для создания автономных ядерных энергоисточников. Высокая надежность, ядерная и радиационная безопасность, самозащищенность реакторной установки определяется эффектами реактивности водоохлаждаемой активной зоны, естественной циркуляцией теплоносителя, принятыми конструкторско-технологическими решениями. Мощность реактора регулируется за счет температурных эффектов реактивности, что предъявляет высокие требования к точности определения теплогидравлических параметров контуров теплоотвода в расчетных обоснованиях стационарных и нестационарных режимов РУ. Наряду с использованием методик расчета, базирующихся на критериальных соотношениях для гидродинамики и теплоотдачи [4], целесообразно разрабатывать и расчетные методы определения теплогидравлических параметров контура теплоотвода, использующие современные CFD коды. Полнота описания процессов такими кодами дает возможность выполнить оптимизацию конструкторских и технологических решений, обоснованно выстроить объекты физического моделирования и программы экспериментов, необходимые для валидации расчетных кодов.
В настоящей статье приводятся результаты использования гидрогазодинамического пакета ANSYS CFX для моделирования одного из элементов теплоотвода РУ "УНИТЕРМ" - промежуточного теплообменника, передающего теплоту активной зоны к теплоносителю промежуточного контура, который в свою очередь отдает теплоту поверхностям парогенератора [2,3]. Использование промежуточного контура практически исключает попадание радиоактивности в контур потребителей пара, но приводит к дополнительному термическому сопротивлению при теплопередаче от твэлов к контуру потребителя РУ.
Необходимые для моделирования размеры поверхностей теплообмена, высоты их расположения в корпусе реактора были получены при предварительных расчетах с использованием критериальных соотношений для реактора с тепловой мощностью 20 МВт. Ограниченные вычислительные ресурсы привели к необходимости на этапе отработки технологии моделирования ограничить число змеевиков промежуточного теплообменника двумя, сохранив их геометрические размеры близкими к необходимым для реактора тепловой мощностью 20 МВт. Тепловая мощность двух змеевиков промежуточного теплообменника равна 300 кВт.
Рисунок 1 - Модель промежуточного теплообменника 1 - внутренняя обечайка; 2 - наружная обечайка; 3 - змеевиковые трубы внутреннего ряда; 4 - змеевиковые
трубы наружного ряда
На рисунке 1 представлена трехмерная модель промежуточного теплообменника, состоящая из внутренней и внешней обечаек, змеевиковых труб. Промежуточный теплообменник состоит из двух рядов змеевиков по 5 параллельных труб диаметром 13х1,5 мм в каждом ряде. Продольный шаг навивки змеевика - 14,5 мм; поперечный -15,5 мм. В радиальном направлении зазоры между змеевиками составляют 2,5 мм, между змеевиками и кожухами - 0,2 мм. Высота змеевиков составляет 1609,5 мм. Вид сверху на модель показан на рисунке 2.
Рисунок 2. Модель промежуточного теплообменника (Вид сверху)
Расчетная сетка для модели строилась с использованием программ Icem CFD и ANSYS Meshing. Расчетная область разбивалась на две подобласти: первая - межтрубное пространство (МТП) - вода между трубами и обечайками; вторая - внутритрубное пространство (ВТП) - вода, текущая в змеевиках. Стенка моделировалась с помощью функции Interface Model - Thin material - Steell при задании контакта сред. Элементы сетки - тетраэдры. Общее число элементов в сетке модели промежуточного теплообменника - 13 970 098.
На рисунке 3 представлена расчетная сетка для модели промежуточного теплообменника.
Необходимые для расчетов расходы, температуры теплоносителей первого и промежуточного контуров приведены в таблице 1. Давление теплоносителя принято равным 16,5 МПа Движение воды в МТП - нисходящее, в ВТП - восходящее.
Рисунок 3 - Расчетная сетка для модели промежуточного теплообменника
2 Выбор моделей турбулентности
В расчете использованы модели турбулентности: k-S RNG и Shear Stress Transport (SST) [4].
Выбор модели турбулентности k-S связан с ее успешным использованием для решения различных инженерных задач, меньшей требовательностью к качеству сетки. Турбулентная вязкость определяется выражением:
к2
14 = СцР —
где C^ - константа, С^ = 0 , 0 9; k - кинетическая энергия;
£ - скорость диссипации кинетической энергии турбулентности.
Для кинетической энергии турбулентности и диссипации используются уравнения:
+ h (р UJ£) =-t + - С*Р£+ СМ
г»/-<£|
где - константы,
С£1=1,44; С£2=1,92; ак=1; <г£=1,3
Рк - генерация турбулентности из-за действия сил вязкости,
Р— ь, Реь - влияние сил плавучести,
D _ Vt dp ь = -^ ^
Р£Ь = С3т ах (О ,/\ь)5 1 п (<р),
С3 = 1 ; бгр=1,
- угол между вектором скорости и вектором ускорения свободного падения. Для модели к-е КЫО уточнены значения коэффициентов:
4.38,
С£1=1,42-£ь /„ = ; п = С„ = 0,08 5; С£2=1,68;
IV - 1+/J3
// = 0, 0 1 2 ; о.=0,7179; а£=0,7179
Модель турбулентности SST - модель переноса касательных напряжений, применяется, когда требуется хорошее решение в пристеночном слое, что важно для расчета теплообмена. Модель SST использует модель к-Ю - для пристеночной области и модель к-е для внешнего течения.
Уравнения для кинетической энергии к и частоты турбулентности Ю имеют вид:
^+(р //' ^+1»-
асрш) д
(д + + ( 1 - Fi ) 2р—1—+ «2 ■- - //2 Р *>2 + Р. ь,
F1 = tan/ ( ar,gi 4) ,
. , , 2Vfc 500v4 4pk ч
ar^i = mm (rnax fe^ ,
CD—w=rn dx ( 0 - 1 о ) ,
где y - расстояние от стенки, м; «2 , //2, °"k2, 0"ш2 - константы,
«2 = 0,44; /2 = 0,082 8; о— 2 = 1; а-2 = 1/0,856. При моделировании течения вблизи стенки использованы пристеночные функции: стандартная (standart), масштабируемая (scalable) и автоматическая (automatic).
Толщина первой ячейки, для получения безразмерной координаты от стенки
10
оценивалась по [5,6]. В [6] используется формула для приближенного расчета у , которая дает значение толщины первой ячейки в первом приближении. После расчета и получения конкретного значения , необходимо изменить сетку в соответствии с рекомендациями к моделям турбулентности [5].
В пристеночной области тепловой поток вычислялся по формуле:
_ рсри
— т+ (т и> ~ т /) . При использовании пристеночной масштабируемой функции:
т + — о, 1 ы (у*) + р.
где .
При использовании автоматической пристеночной функции: Т+ — Рг у* е(-г) + [0, 1 1п( у*) + ¡3] е- х/г,
0,01 {Рг у*)*
где
1+5 Рг3 у*
3 Результаты моделирования
В качестве граничных условий задавались: давление, температура на входе и расход сред (таблица 1). Модель теплопереноса - Thermal Energy. Сходимость решения уравнений порядка 10-4. Свойства воды брались из библиотеки IAPWS-IF97 [7].
Таблица 1 - Задание граничных условий для модели промежуточного теплообменника
Межтрубное пространство (МТП)
Давление на входе, МПа 16,5
Температура на входе, °С 335,0
Шероховатость, мкм 10,0
Расход, кг/с 0,751
Внутритрубное пространство (ВТП)
Давление на входе, МПа 16,5
Температура на входе, °С 194,9
Шероховатость, мкм 10,0
Расход, кг/с 0,4515
Сводка балансных результатов расчетов промежуточного теплообменника, выполненных с использованием моделей к-е КЫО, приведена в таблице 2.
Промежуточный теплоноситель на выходе из труб (ВТП) нагревается с входной температуры 194,4 °С до средних значений температур 325,7 °С и 320,2 °С. Средние скорости его движения в змеевиках составляют 0,91 и 0,88 м/с. Гидравлические потери в змеевиках различается существенно и составляют 9,8 кПа для к-е NN0 и 5,6 кПа для 88Т модели. Теплоноситель первого контура (МТП) на выходе из межтрубного пространства теплообменника охлаждается с 335 °С до 266,8 °С и 272,3 °С, движется со средней
скоростью 0,33 и 0,45 м/с. Гидравлические потери при обтекании змеевиков составляют 1,26 кПа и 1,77 кПа для моделей к-е NN0, ££Т.
Различие в количестве переданной в теплообменнике теплоты к теплоносителю промежуточного контура от теплоносителя первого контура составляет ~7%. Тепловой баланс теплообменника лучше выполняется для расчетов с использованием 88Т модели. Полученный результат обусловлен с одной стороны достаточно "грубо"расчетной сеткой, с другой - подтверждает потенциальные преимущества модели ББТ для расчета теплообмена за счет лучшего описания пристенной области.
Таблица 2 - Результаты теплового расчета промежуточного теплообменника
Параметр Модель к-е ШО Модель 88Т
На выходе из ВТП
Температура, °С 325,68 320,19
Скорость, м/с 0,908 0,884
Давление, МПа 16,4902 16,4944
Перепад давления, Па 9781 5600
У+ 207,86 149,32
На выходе из МТП
Температура, °С 266,80 272,31
Скорость, м/с 0,331 0,453
Давление, МПа 16,4991 16,4982
Перепад давления, Па 1264 1774
У+ 172,19 174,46
Анализ распределений скоростей, температур в теплоносителях первого и промежуточного контуров модели позволил установить особенности гидродинамики и теплообмена в моделируемом змеевиковом теплообменнике. При этом учитывалось, что при выбранных геометрических размерах обечаек и змеевиковых труб, гидравлические диаметры трактов теплоносителя, движущегося между трубами, равны 9,02 мм для трактов между змеевиками и 3,53 мм для трактов у кожухов. Поэтому при неучете перемешивания между трактами скорость потока в трактах у кожухов [8] составляет ~0,58 от скорости потока в тракте между трубами, а расходы в трактах у кожухов ~0,18 от расхода в тракте между трубами. Перемешивание потоков, развивающихся в
параллельных трактах, уменьшает разницу в расходах. Тем не менее, меньший расход в прикожуховой области теплообменника приводит к большему охлаждению теплоносителя в этой области, что и наблюдается в расчетах.
Указанная особенность геометрии области течения теплоносителя первого контура в межтрубном пространстве промежуточного теплообменника может привести к разверке температуры трубного пучка по периметру трубы, неоднородности температуры теплоносителя промежуточного контура на выходе из отдельных труб, что и определило направление анализа результатов моделирования.
Во-первых установлено, что движение теплоносителя первого контура в межтрубном пространстве (рисунок 4) происходит преимущественно вдоль змеевиков труб. Не содержащие векторов скорости области эллипсной формы на рисунке 4 соответствуют поперечным сечениям змеевиков труб (5 труб в ряду).
Рисунок 4 - Векторы скорости в поперечном сечении межтрубного пространства промежуточного теплообменника: а) расчеты по модели к-£ КЫС'., б) расчеты по модели 88Т
В тракте между змеевиками, где скорости максимальны, угол ориентации части потока, удаленной от стенки, к горизонту составляет ~27,1° угловых (угол подъема спирали труб внутреннего и периферийного рядов составляет соответственно ~11,5° угловых и ~9,1° угловых). В трактах, прилегающих к кожухам, угол ориентации удаленной от стенки части потока к горизонту, составляет соответственно ~21,6° угловых и ~23,7° угловых. Приведенные значения соответствуют результатам моделирования с использованием модели ББТ. Сводка данных о геометрии течения в пространстве между трубами приведена в таблице 3 и на рисунке 5.
змеевик
внешнего ряда
змеевик
внутреннего ряда
линии тока
Рисунок 5. Линии тока в продольном сечении межтрубного пространства Таблица 3 - Характеристика геометрии трактов и течения потока в межтрубном пространстве (МТП)
Угол ориентации к горизонту Модель к-е 88Т
Змеевиков внутреннего ряда 11,5
Змеевиков наружного ряда 9,1
Потока у внутренней стенки 21,6
Потока в тракте между змеевиками 27,1
Потока у наружной стенки 23,7
Таким образом, потоки у обечайки имеют меньшие углы ориентации к горизонту, что увеличивает из охлаждение, равно как и меньше гидравлические диаметры трактов.
Распределение скорости теплоносителя промежуточного контура, движущегося внутри труб, приведено на рисунке 6. Влияние массовых сил, связанных с кривизной тракта, малосущественно: распределение скорости симметричны относительно оси трубы.
Рисунок 6 - Векторы скорости в поперечном сечении внутритрубного пространства промежуточного
теплообменника
Расчеты показали, что максимальные величины температуры теплоносителя первого контура на выходе из поверхности теплообменника имеют место в тракте между змеевиками (рисунок 7), а наименьшее значение этих температур - в окрестности кожухов. Этот результат согласуется с гидродинамикой течения в межтрубном пространстве.
Рисунок 7 - Распределение температуры на выходе из межтрубного пространства промежуточного теплообменника: а) расчеты по модели к-£ КЫС'., б) расчеты по модели 88Т
Температуры нагреваемого теплоносителя промежуточного контура на выходе из труб (рисунок 8) практически одинаковы для всех пяти параллельных труб каждого из
змеевиков, но различны для змеевиков внутреннего первого ряда и внешнего ряда. Для обеих моделей расчета более высокие температуры наблюдаются для змеевиков периферийного ряда. Это может быть связано с большей площадью его трубных поверхностей (на 27%) . Однако разница в температурах на выходе из змеевиков невелика и составляет ~2,4 °С для модели к-е КМО, ~4,1 °С для модели ££Т, что составляет менее 3% от величины подогрева теплоносителя промежуточного контура, равного ~128 °С.
Рисунок 8 - Распределение температуры на выходе из внутритрубного пространства промежуточного теплообменника: а) расчеты по модели к-е КЫС'., б) расчеты по модели 88Т
Распределение температуры на наружной поверхности змеевиков промежуточного теплообменника приведено на рисунке 9. Максимальное различие температур поверхности составляет ~5 °С, меньшие значения температур поверхности получаются при использовании 88Т модели расчета. Значения температур поверхности промежуточного теплообменника лежат в пределах 195...328 °С.
Рисунок 9 - Распределение температуры на поверхности змеевика: а) расчеты по модели к-е RNG; б)
расчеты по модели SST
4 Анализ результатов моделирования
Выполненное моделирование с использованием программного кода ANSYS и моделей к-е RNG, SST характеристик турбулентности позволило установить важные особенности течения и теплообмена в змеевиковых поверхностях модельного промежуточного теплообменника применительно к РУ "УНИТЕРМ". Показано, что течение в межтрубном пространстве змеевиковой поверхности формируется как продольное вдоль осей труб спиральное течение, а не как поперечное обтекание коридорного пучка труб. Оценки показывают, что гидравлические потери при продольном течении потока вдоль змеевиков оказываются меньшими, чем при традиционно рассматриваемом поперечном обтекании пучков прямых, ориентированных нормально потоку труб. Теплоотдача от отдельных областей течений в трактах между змеевиками и кожухами к теплоносителю промежуточного контура оказывается различной, что объясняется гидродинамикой течения и выбранными размерами зазоров между трубами рядов змеевиков, змеевиками и кожухами. Несмотря на отмеченное, разверка температур промежуточного теплоносителя на выходе из отдельных змеевиков, различие в температурах теплоносителя первого контура лежит в пределах 5...9 °С. Эта разверка может быть уменьшена при оптимизации геометрии зазоров змеевиковой поверхности и кожуха.
Количественные различия в результатах расчетов с использованием моделей к-е RNG, SST связаны не столько с особенностями замыкающих соотношений, сколько с использованием достаточно грубой расчетной сетки, в которой первые пристенные конечные элементы соответствовали значениям y+—150...200. При этом каждый из выполненных расчетов потребовал около 40 часов с использованием компьютера с процессором intel core i7, оперативной памятью - 16 Гб.
Список литературы
1. Драгунов Ю.Г., Шишкин В.А., Гречко Г.И., Гольцов Е.Н. Малая ядерная энергетика: задачи и ответы // Атомная энергия. 2011. Т. 111, вып. 5. С. 294-297.
2. Ачкасов А.Н., Гольцов Е.В., Гречко Г.И., Кузнецов Ю.Н. Реакторные установки для атомных станций малой мощности // Атомная энергия. 2012. Т. 113, вып. 1. С. 43-48.
3. Солонин В.И., Шишкин В.А. Атомные станции малой мощности - перспективный источник локального электро-теплоснабжения регионов Севера и Востока России // Международная конференция «Региональные проблемы в стратегии устойчивого развития России»: сб. статей. М.: «Аванта+», 1999. С. 300-304.
4. ANSYS CFX-Solver Theory Guide. SAS IP, Inc., 2011. 402 p.
5. Волкинд Д.К. Оценка качества результатов численного моделирования // Пользовательская конференция ANSYS 2013 «Направления развития инженерного анализа» (Москва, 4-5 июня 2013 г.): материалы. ANSYS, Inc.; «Делкам Урал» - «ПЛМ Урал», 2013. С. 56.
6. Compute Grid Spacing for a Given Y+. Pointwise, Inc.: website. Available at: http://www.pointwise.com/yplus , accessed 01.09.2014.
7. IAPWS-IF97. Release on the IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam. The International Association for Properties of Water and Steam, Erlangen, September 1997. 48 p.
8. Емельянов И.Я., Михан В.И., Солонин В.И. Конструирование ядерных реакторов: учеб. пособие для вузов / под общ. ред. акад. Н.А. Доллежаля. М.: Энергоиздат, 1982. 400 с.
Science and Education of the Bauman MSTU, 2014, no. 10, pp. 398-412.
DOI: 10.7463/1014.0727220
Received: Revised:
02.09.2014 21.09.2014
Science ^Education
of the Bauman MSTU
ISSN 1994-0448 © Bauman Moscow State Technical Unversity
Modeling of Heat Transfer in the Helical-Coil Heat Exchanger for the Reactor Facility "UNITERM"
J,* a satin9 0 'gmail.ru
:Bauman Moscow State Technical University, Moscow, Russia
V.I. Solonin1 , A.A. Satin1
Keywords: reactor facility UNITERM, heat sink circuit, core, helical-coil heat exchange,
computational hydrodynamic and heat exchange
Circuit heat sink plays an important role in the reactor system. Therefore it imposes high requirements for quality of determining thermal-hydraulic parameters. This article is aimed at modeling of heat exchange process of the helical-coil heat exchanger, which is part of the heat sink circuit of the reactor facility "UNITERM."
The simulation was performed using hydro-gas-dynamic software package ANSYS CFX. Computational fluid dynamics of this package allows us to perform calculations in a three-dimensional setting, giving an idea of the fluid flow nature. The purpose of the simulation was to determine the parameters of the helical-coil heat exchanger (temperature, velocity at the outlet of the pipe and inter-tubular space, pressure drop), and the nature of the fluid flow of primary and intermediate coolants. Geometric parameters of the model were determined using the preliminary calculations performed by the criterion equations. In calculations Turbulence models k-s RNG, Shear Stress Transport (SST) are used. The article describes selected turbulence models, and considers relationship with wall function.
The calculation results allow us to give the values obtained for thermal-hydraulic parameters, to compare selected turbulence models, as well as to show distribution patterns of the coolant temperature, pressure, and velocity at the outlet of the intermediate cooler.
Calculations have shown that:
- maximum values of primary coolant temperature at the outlet of the heat exchanger surface are encountered in the space between the helical-coil tubes;
- higher temperatures of intermediate coolant at the outlet of the coils (in space of helical-coil tubes) are observed for the peripheral row;
- primary coolant movement in the inter-tubular space of helical-coil surface is formed as a spiral flow, rather than as a in-line tube bank cross flow.
References
1. Dragunov Iu.G., Shishkin V.A., Grechko G.I., Gol'tsov E.N. Low-capacity nuclear power: questions and answers. Atomnaia energiia, 2011, vol. 111, no. 5, pp. 294-297. (English translation: Atomic Energy, 2012, vol. 111, iss. 5, pp. 366-369. DOI: 10.1007/s10512-012-9505-x ).
2. Achkasov A.N., Gol'tsov E.V., Grechko G.I., Kuznetsov Iu.N. Reactor facilities for small nuclear power plants. Atomnaia energiia, 2012, vol. 113, no. 1, pp. 43-48. (English translation: Atomic Energy, 2012, vol. 113, iss. 1, pp. 51-56. DOI: 10.1007/s10512-012-9594-6 ).
3. Solonin V.I., Shishkin V.A. Small nuclear power plants are promising source of local electricity and heat supply of the North and East of Russia. Mezhdunarodnaia konferentsiia "Regional'nye problemy v strategii ustoichivogo razvitiia Rossii": sb. statei [International conference "Regional problems in the sustainable development strategy of Russia": collect. of papers]. Moscow, "Avanta+" Publ., 1999, pp. 300-304. (in Russian).
4. ANSYS CFX-Solver Theory Guide. SAS IP, Inc., 2011. 402 p.
5. Volkind D.K. Quality assessment of the results of numerical modeling. Pol'zovatel'skaia konferentsiia ANSYS 2013 "Napravleniia razvitiia inzhenernogo analiza": materialy [ANSYS User Conference 2013 "Current trends in engineering analysis": proc.]. Moscow, 4-June 2013. Publ. of ANSYS, Inc.; "Delkam Ural" - "PLM Ural", 2013, p. 56. (in Russian).
6. Compute Grid Spacing for a Given Y+. Pointwise, Inc.: website. Available at: http://www.pointwise.com/yplus , accessed 01.09.2014.
7. IAPWS-IF97. Release on the IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam. The International Association for Properties of Water and Steam, Erlangen, September 1997. 48 p.
8. Dollezhal' N.A., ed., Emel'ianov I.Ia., Mikhan V.I., Solonin V.I. Konstruirovanie iadernykh reaktorov [Design of nuclear reactors]. Moscow, Energoizdat Publ., 1982. 400 p. (in Russian).
