ТЕПЛОМАССОПЕРЕНОС В ЯДЕРНО-ЭНЕРГЕТИЧЕСКИХ УСТАНОВКАХ

Август 31, 2022

Главная
Физика
ТЕПЛОМАССОПЕРЕНОС В ЯДЕРНО-ЭНЕРГЕТИЧЕСКИХ УСТАНОВКАХ

Содержание

2. ТЕПЛОМАССОПЕРЕНОС В ЯДЕРНО-ЭНЕРГЕТИЧЕСКИХ УСТАНОВКАХ Тема №12 Особенности процессов гидродинамики и теплообмена в активных зонах реакторов
3. Методы теплофизического расчёта стержневых ТВС 03
4. 04
5. 05
6. Термогидродинамический расчёт ТВС проводится с целью определения параметров потока теплоносителя, охлаждающего сборку. На рисунках 1 и
7. 07
8. 08
9. Рисунок 3 – Поперечные сечения сборок стержней с квадратной (А) и гексагональной (Б) решётками 09
10. 10
11. 11
12. 12
13. 13
14. Исходная (а) и трансформированная (рассчитываемая) (б) области в методе согласованного с формой границы преобразования координат 14
15. Пример дискретизации пучка стержней в методе, основанном на приближении пористого тела 15
16. При выборе способа дискретизации межстержневого пространства пучка на субканалы следует руководствоваться следующими соображениями: 1) взаимодействие через
17. Таким образом, субканал – часть занятого теплоносителем пространства сборки, ограниченная (рисунки 5а, 5б и 6а и
18. Рисунок 4 – К определению термина "меж-твэльная ячейка» 18
19. Рисунок 5а – Фрагмент из участков трёх твэлов, расположенных в гексагональной упаковке, и заключённого между ними
20. Рисунок 5б – Фрагмент из участков трёх твэлов, расположенных в гексагональной упаковке, и заключённого между ними
21. Рисунок 6а – Центральная, боковая и угловая межтвэльные ячейки 21
22. 22
23. Рисунок 6б – Центральная, боковая и угловая межтвэльные ячейки 23
24. Такой подход вполне корректен, если ограничиваться ● расчётом параметров потока, усредненных по поперечному сечению субканала, и
25. Второе условие, с учётом сделанного выше предположения о неизменности плотности теплоносителя в пределах элементарной ячейки, справедливо
26. Рисунок 7 – Поперечные сечения контрольных объёмов, используемых при решении субканальных уравнений баланса массы, энергии и
27. Между субканалами происходит обмен массой, импульсом и энергией – так называемое "межканальное перемешивание". Перемешивание обусловлено различными
28. 28
29. По полноте и точности описания процессов гидродинамики и тепломассообмена субканальные коды можно классифицировать на основании анализа
30. Рисунок 9 – Этапы разработки компьютерного кода, реализующего выбранную математическую модель 30
31. Математические модели теплогидравлических процессов в пучках стержней будем классифицировать согласно схеме, представленной на рисунке 10. Пучки
32. Рисунок 10 – Схема классификации субканальных моделей 32
33. В топливных сборках реакторов на быстрых нейтронах (БН, LMFBR) твэлы дистанционируются, как правило, спиральными проволочными навивками
34. Количество "комплектов" решаемых уравнений сохранения зависит от того, рассматривается ли теплоноситель как одна фаза (жидкость) или
35. Коэффициент межканального перемешивания субстанцией Ψ между субканалом i и смежным с ним субканалом j определим, аналогично
36. Основные тенденции в развитии субканальных моделей Cуществует большое количество субканальных моделей. Это разнообразие объясняется, с одной
37. Очевидно, что даже при анализе одного типа аппаратов допустимо иметь набор математических моделей различной степени сложности
38. 38
39. ОСОБЕННОСТИ ПРОЦЕССОВ ТЕПЛООБМЕНА В РАЗЛИЧНЫХ РЕЖИМАХ РАБОТЫ РЕАКТОРА ЭНЕРГОВЫДЕЛЕНИЕ В РАЗЛИЧНЫХ РЕЖИМАХ 39
40. В переходных процессах наибольшие деформации поля энерговыделения возникают при умеренном снижении мощности (40÷50%) и последующем восстановлении
41. Рисунок 2 – Изменение распределения энерговыделения в реакторе ВВЭР-1000 в режимах разной мощности: 1 – 100%;
42. Рисунок 3 – Изменение потока нейронов (Ф) и плотности теплового потока (q) в ВВЭР-440 при удалении
43. 43
44. ● Высокая температура, ● низкая теплопроводность, ● достаточно высокая теплоёмкость диоксида урана приводят к тому, что
45. Если реактор до выключения работал в течение времени τс на мощности N0 , то мощность остаточного
46. Приведенные в таблице значения мощности малы по сравнению с номинальными значениями. Но они значительны в абсолютных
47. ИЗМЕНЕНИЕ ТЕМПЕРАТУРЫ В ПЕРЕХОДНЫХ ПРОЦЕССАХ Под термином "переходные процессы" подразумеваются процессы перехода элементов реактора с одного
48. Плотность теплового потока на поверхности: Средняя по объёму температура цилиндра: (1) (2) 48
49. Подставив (1) в (2), получим: Среднее термическое сопротивление цилиндрического твэла без оболочки: Уравнение теплового баланса для
50. 50
51. С другой стороны, плотность теплового потока на поверхности цилиндрического твэла без оболочки можно выразить через термическое
52. – постоянная времени, определяющая динамические характеристики цилиндрического твэла. Можем записать τ0=τвн+τн , где стоящие в правой
53. Критериальное число Био – критерий краевого подобия, характеризующий связь между полем температур в твёрдом теле и
54. 54
55. В чём физический смысл постоянной времени? Если начальное состояние твэла стационарное (то есть F(0)=1; dt/dτ=0), то
56. Аналогично, для мгновенного изменения энерговыделения на ΔF имеем Таким образом, τ0 есть постоянная экспоненты, описывающей изменение
57. Для линейного по времени изменения температуры теплоносителя Δtж=mτ изменение температуры твэла Максимальное отставание температуры твэла от
58. Рисунок 6 – Изменение средней температуры цилиндрического твэла без оболочки при линейном изменении температуры теплоносителя 58
59. Постоянную времени для тел любой формы можно записать в общем виде (10) где V и F
60. Рисунок 7 – Распределение температуры в твэле с газовым зазором между топливом и оболочкой 60
61. Рассмотрим теперь задачу охлаждения разогретой без учёта ее внутреннего тепловыделения. Пусть АЗ охлаждается жидкостью с расходом
62. В качестве примера рассмотрим характерные времена запаздывания для ВВЭР-440. Основные параметры: объём топлива VT=4.5 м3; радиус
63. Время запаздывания ‒ τвн ‒ вследствие внутреннего термического сопротивления твэла можно оценить так: Время запаздывания, обусловленное
64. Время прохода теплоносителя через активную зону: τтр = H/w = 2.5/3.4 ≈ 0.735 с. Теплоёмкость топлива:
65. ТЕПЛОВЫЕ УДАРЫ И ПЕРИОДИЧЕСКИЕ КОЛЕБАНИЯ ТЕМПЕРАТУРЫ В переходных режимах температура в разных частях конструкционных элементов изменяется
66. 66
67. При быстром снижении температуры теплоносителя внутренние слои трубы охлаждаются быстрее, чем наружные. Распределение температуры приобретает вид,
68. (15) Максимальные напряжения возникают на поверхностях конструкционных элементов. 68
69. Под термином "тепловой удар" понимается однократное изменение температуры поверхности. Однако в потоке теплоносителя существуют колебания температуры,
70. Для определения температурных напряжений необходимо найти распределение температуры в стенке толщиной h, на одной из поверхностей
71. ОСОБЕННОСТИ ТЕМПЕРАТУРНЫХ РЕЖИМОВ ПРИ РАССЛОЕНИИ ТЕПЛОНОСИТЕЛЯ В переходных, а особенно в аварийных режимах практически всегда наблюдается
72. Опасность стратифицированных течений – большой градиент температуры в узкой области, что может вызвать чрезмерные термические напряжения
73. Рисунок 9 – Расслоение теплоносителя в верхней камере реактора: 1 – активная зона; 2 – уровень
74. Тепло, передаваемое через поверхность раздела между горячей и холодной областями, быстро размывается турбулентными вихрями, и градиент
75. Числами подобия, определяющими характер течений являются также числа Фруда и Рейнольдса: Другим примером расслоенных течений является
76. 76
77. Вопросы, выносимые на зачёт 1. Характер распределения температуры в цилиндрическом твэле с зазором между топливом и
79. Скачать презентацию

Слайд 2

ТЕПЛОМАССОПЕРЕНОС В ЯДЕРНО-ЭНЕРГЕТИЧЕСКИХ УСТАНОВКАХ
Тема №12
Особенности процессов
гидродинамики и теплообмена
в активных зонах реакторов

Слайд 3

Методы
теплофизического расчёта
стержневых ТВС
03

Слайд 4

04

Слайд 5

05

Слайд 6

Термогидродинамический расчёт ТВС проводится с целью определения параметров потока теплоносителя, охлаждающего

сборку.
На рисунках 1 и 2 представлены тепловыделяющие сборки ректоров ВВЭР-440 и ВВЭР-1000, соответственно.
Широкий спектр задач теплогидравлических расчётов стержневых тепловыделяющих сборок решается с помощью методов, основанных на субканальном приближении.
Идея его сводится к следующему. Пространство сборки (рисунок 3), занятое теплоносителем, представляет собой многосвязную область, которую можно рассматривать как совокупность взаимосвязанных субканалов сложной формы.
Субканалом называем часть занятого теплоносителем пространства сборки.

Слайд 7

07

Слайд 8

08

Слайд 9

Рисунок 3 – Поперечные сечения сборок стержней с квадратной (А) и

гексагональной (Б) решётками

Слайд 10

10

Слайд 11

11

Слайд 12

12

Слайд 13

13

Слайд 14

Исходная (а) и трансформированная (рассчитываемая) (б) области в методе согласованного с

формой границы преобразования координат

Слайд 15

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

Слайд 16

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

следующими соображениями:
1) взаимодействие через границы между субканалами должно быть минимальным;
2) должна быть возможной разумная аппроксимация функций, описывающих характеристики теплоносителя внутри субканала.
Заметим, что пучки твэлов являются гидравлически "тесными": относительный шаг t решётки стержней, определяемый как отношение расстояния s между центрами соседних стержней (шага решётки) к диаметру d стержня, много меньше 2,0, т.е.:
(1.1)
Как правило, t≤1.4.
Для ТВС реакторов на быстрых нейтрона, t≤1.2. Шаг решётки в сборках активной зоны меньше, чем в сборках зоны воспроизводства.
Применительно к таким пучкам стержней вышеуказанным условиям соответствует разбиение на субканалы, совпадающие в поперечном сечении с межтвэльными ячейками (рисунок 4).

Слайд 17

Таким образом, субканал – часть занятого теплоносителем пространства сборки, ограниченная (рисунки

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

Слайд 18

Рисунок 4 –
К определению термина "меж-твэльная ячейка»
18

Слайд 19

Рисунок 5а – Фрагмент из участков трёх твэлов, расположенных в гексагональной

упаковке, и заключённого между ними субканала

Слайд 20

Рисунок 5б – Фрагмент из участков трёх твэлов, расположенных в гексагональной

упаковке, и заключённого между ними субканала

Слайд 21

Рисунок 6а – Центральная, боковая и угловая межтвэльные ячейки
21

Слайд 22

22

Слайд 23

Рисунок 6б – Центральная, боковая и угловая межтвэльные ячейки
23

Слайд 24

Такой подход вполне корректен, если ограничиваться
● расчётом параметров потока, усредненных по

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

Слайд 25

Второе условие, с учётом сделанного выше предположения о неизменности плотности теплоносителя

в пределах элементарной ячейки, справедливо лишь, если аксиальная составляющая скорости потока существенно превосходит поперечные.
Это позволяет, в свою очередь, провести разделение переменных в уравнении движения и рассматривать отдельно уравнения баланса продольной и поперечных составляющих импульса.
Как следствие, существенно упрощается решение задачи. Необходимо отметить, что для аксиальной и поперечной составляющих импульса приходится использовать различные контрольные объёмы (рисунок 7).

Слайд 26

Рисунок 7 – Поперечные сечения контрольных объёмов, используемых при решении субканальных

уравнений баланса массы, энергии и аксиальной составляющей импульса (a) и поперечных составляющих импульса (b)

Слайд 27

Между субканалами происходит обмен массой, импульсом и энергией – так называемое

"межканальное перемешивание".
Перемешивание обусловлено различными факторами конструктивного и режимного характера (рисунок 8).

Слайд 28

28

Слайд 29

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

можно классифицировать на основании анализа реализованных в них математических моделей.
Причину этого легко понять, рассмотрев рисунок 9, на котором представлена схема разработки программы для ЭВМ, реализующей некоторую математическую модель.
Как видно, при разработке кода могут возникнуть обстоятельства, ограничивающие, в сравнении с используемой моделью, область применения. С другой стороны, использование некоторых дополнительных математических процедур или моделей, применение более совершенных численных схем и алгоритмов, использование более прогрессивных языков программирования, учёт возможностей новых трансляторов и операционных систем способны значительно расширить диапазон решаемых с помощью данной программы задач, при этом реализуемая в ней модель не претерпит никаких изменений.

Слайд 30

Рисунок 9 – Этапы разработки компьютерного кода, реализующего выбранную математическую модель
30

Слайд 31

Математические модели теплогидравлических процессов в пучках стержней будем классифицировать согласно схеме,

представленной на рисунке 10.
Пучки твэлов можно разделить по способу дистанционирования.
В ТВС водоводяных реакторов (ВВЭР, PWR, BWR) в качестве дистанционирующих элементов обычно используют решётки, расположенные перпендикулярно пучку.
В этом случае твэльную сборку рассматривают как решётку стержней с дополнительными местными сопротивлениями, а для многофазного теплоносителя ещё и дополнительными поверхностями осаждения и конденсации.

Слайд 32

Рисунок 10 – Схема классификации субканальных моделей
32

Слайд 33

В топливных сборках реакторов на быстрых нейтронах (БН, LMFBR) твэлы дистанционируются,

как правило, спиральными проволочными навивками (реже – спиральными рёбрами). Корректное описание течения теплоносителя в областях с такой геометрией – сложная задача. Наличие же спиральной навивки – главный фактор, влияющий на распределение параметров теплоносителя.
Математическая модель разрабатывается для исследования конкретного типа сборки. Это объясняется хотя бы тем, что кроме решаемых уравнений модель включает в себя, в частности, соотношения, описывающие теплофизические свойства теплоносителя.
Очевидно, что изначально следует разделить модели на два класса: пучки со стержнями без навивки и пучки с твэлами с навивкой.

Слайд 34

Количество "комплектов" решаемых уравнений сохранения зависит от того, рассматривается ли теплоноситель

как одна фаза (жидкость) или же как многофазная (многожидкостная) среда. Поэтому в рамках каждого из выбранных классов выделим семейства "одножидкостных" и "многожидкостных" моделей.
В названных семействах следует определить по два рода моделей.
"Одножидкостные" модели описывают процессы в однофазном теплоносителе или/и в двухфазном, рассматриваемом как гомогенная среда.
Модели "многожидкостные" используют одно из приближений:
1) "дрейф фаз" ("скорость дрейфа" );
2) "раздельное течение фаз". Каждую же из групп "одножидкостного" семейства классифицируем по типу решаемых систем уравнений на стационарные и нестационарные.
Стационарные модели разделим на три группы:
1) использующие приближение «изобарного сечения» и коэффициенты межканального перемешивания,
2) использующие приближение «изобарного сечения»,
3) свободные от перечисленных ограничений и включающие в себя уравнение баланса поперечных составляющих импульса.

Слайд 35

Коэффициент межканального перемешивания субстанцией Ψ между субканалом i и смежным с

ним субканалом j определим, аналогично [17, 18], как отношение поперечного потока субстанции Ψ через имеющий единичную длину участок межстержневого зазора к потоку субстанции Ψ в осевом направлении через поперечное сечение субканала i , а именно
где ‒ коэффициент межканального перемешивания субстанцией Ψ между субканалом i и смежным с ним субканалом j;
‒ поперечный поток субстанции Ψ через имеющий единичную длину
участок межстержневого зазора; ‒ поток субстанции Ψ в осевом направлении через поперечное сечение субканала i.
Вопросу корректного определения коэффициентов межканального перемешивания на основе обработки данных натурных экспериментов уделено большое внимание в работах Физико-энергетического института (ГНЦ "ФЭИ", г. Обнинск, Россия) [18‒29].

Слайд 36

Основные тенденции в развитии субканальных моделей
Cуществует большое количество субканальных моделей.
Это разнообразие

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

Слайд 37

Очевидно, что даже при анализе одного типа аппаратов допустимо иметь набор

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

Слайд 38

38

Слайд 39

ОСОБЕННОСТИ ПРОЦЕССОВ ТЕПЛООБМЕНА
В РАЗЛИЧНЫХ РЕЖИМАХ РАБОТЫ РЕАКТОРА
ЭНЕРГОВЫДЕЛЕНИЕ В РАЗЛИЧНЫХ РЕЖИМАХ
39

Слайд 40

В переходных процессах наибольшие деформации поля энерговыделения возникают при умеренном снижении

мощности (40÷50%) и последующем восстановлении до номинальной.
Пример такого распределения для ВВЭР-1000 показан на рисунке 2. Это довольно медленные процессы, в которых тепловые характеристики следуют за физическими параметрами реактора.

Слайд 41

Рисунок 2 – Изменение распределения энерговыделения в реакторе ВВЭР-1000 в режимах разной

мощности:
1 – 100%;
2 – 0 ч, 50%;
3 – 4.62 ч, 50%; работа 5.3 ч на 50% мощности;
4 – 5.3 ч, 100%; 5 – 10.6 ч, 100%; 6 – 15.8ч, 100%.
ζ = Z/Н, F(ζ) = q(ζ)/q.

Слайд 42

Рисунок 3 – Изменение потока нейронов (Ф) и плотности теплового потока

(q) в ВВЭР-440
при удалении одной ТВС за время 0.1 с

При быстрых процессах принципиальным моментом является тот факт, что тепловые характеристики активной зоны изменяются медленнее, чем физические.
Так, при выбросе из активной зоны одной ТВС мощность ВВЭР-440 возрастает в 15 раз за время 0.1 с, а тепловой поток увеличивается всего в 1.34 раза (рисунок 3).

Слайд 43

43

Слайд 44

● Высокая температура,
● низкая теплопроводность,
● достаточно высокая теплоёмкость
диоксида урана приводят

к тому, что в топливе аккумулируется большое количество тепла.
При аварийных остановах реактора, связанных с ухудшением теплоотвода, это тепло может повысить температуру оболочки твэла до недопустимого значения.
При использовании металлического урана или керметного топлива, обладающих более высокой теплопроводностью, чем диоксид урана, количество тепла, запасённое в твэле, значительно меньше, и температура оболочки твэла при отключении реактора оказывается в пределах допустимой.
Остаточное тепловыделение в реакторе после его останова обусловлено двумя процессами.
● В начальные моменты времени (до 100 с) после останова продолжаются процессы деления на запаздывающих нейтронах.
● В дальнейшем основной вклад вносит тепловыделение от распадающихся радиоактивных осколков деления. Периоды полураспада осколков составляют от нескольких секунд до тысяч лет, поэтому в остановленном реакторе и в отработавшем топливе после прекращения деления ядер будет выделяться энергия распада осколков. Это вызывает необходимость постоянного охлаждения реактора после останова, а также охлаждения отработавшего топлива в хранилищах.

Слайд 45

Если реактор до выключения работал в течение времени τс на мощности

N0 , то мощность остаточного тепловыделения N уменьшается по закону
В таблице показано, как с течением времени после останова спадает мощность остаточного тепловыделения.

Слайд 46

Приведенные в таблице значения мощности малы по сравнению с номинальными значениями.

Но они значительны в абсолютных величинах.
Так, в реакторе с электрической мощностью 1000 МВт, который долго работал на полной мощности, мощность остаточного тепловыделения через час после остановки составит 48 МВт, через сутки около 15 МВт, через год около 0,8 МВт.
Отсюда следует необходимость постоянного охлаждения остановленного реактора.

Слайд 47

ИЗМЕНЕНИЕ ТЕМПЕРАТУРЫ В ПЕРЕХОДНЫХ ПРОЦЕССАХ
Под термином "переходные процессы" подразумеваются процессы перехода

элементов реактора с одного температурного уровня на другой.
Расчёт (оценочный) большинства таких процессов, например, с уменьшением расхода, может быть выполнен на основе квазистационарных подходов.
Для таких расчётов необходимо знать тепловой поток и расход в каждый момент времени. Коэффициенты теплообмена рассчитываются по зависимостям для соответствующего стационарного режима, что не приводит к большим ошибкам, так как изменения коэффициентов теплообмена обычно относительно невелики.
Далее приведены выражения для оценки основных параметров таких процессов.
Большинство твэлов ядерных реакторов имеют вид цилиндрических стержней. Распределение температуры в цилиндре с плотностью объёмных источников тепла qV [Вт/м3] имеет вид

Слайд 48

Плотность теплового потока на поверхности:
Средняя по объёму температура цилиндра:
(1)
(2)
48

Слайд 49

Подставив (1) в (2), получим:
Среднее термическое сопротивление цилиндрического твэла без оболочки:
Уравнение

теплового баланса для элемента с внутренним тепловыделением (для нестационарного процесса), усредненное по объему, запишем в виде:

(5)

(3)

(4)

Слайд 50

50

Слайд 51

С другой стороны, плотность теплового потока на поверхности цилиндрического твэла без

оболочки можно выразить через термическое сопротивление:

Сопоставив эти два уравнения, получим

(6)

Подставив (6) в (5), получим

(7)

Слайд 52

– постоянная времени, определяющая
динамические характеристики
цилиндрического твэла.
Можем записать τ0=τвн+τн ,

где стоящие в правой части слагаемые это, соответственно, постоянные времени для процессов внутри твэла и процессов теплообмена с теплоносителем («наружных процессов»):

(9)

(8)

Слайд 53

Критериальное число Био – критерий краевого подобия, характеризующий связь между полем

температур в твёрдом теле и условиями теплоотдачи на его поверхности, являясь мерой внутреннего и внешнего термического сопротивления (термического сопротивления стенки и термического сопротивления на границе стенки с потоком):
где L – характерный размер.

Слайд 54

54

Слайд 55

В чём физический смысл постоянной времени?
Если начальное состояние твэла стационарное (то

есть F(0)=1; dt/dτ=0), то из (7) следует
Пусть в момент τ=0 температура скачком увеличивается на Δtж.
Температура твэла при этом будет изменяться по закону

Слайд 56

Аналогично, для мгновенного изменения энерговыделения на ΔF имеем
Таким образом, τ0 есть

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

Слайд 57

Для линейного по времени изменения температуры теплоносителя Δtж=mτ изменение температуры твэла
Максимальное

отставание температуры твэла от температуры теплоносителя при τ→∞ составляет m·τ0, то есть это отставание равно изменению температуры теплоносителя за время τ0 (см. рисунок 6).

Слайд 58

Рисунок 6 – Изменение средней температуры цилиндрического твэла без оболочки при

линейном изменении температуры теплоносителя

Слайд 59

Постоянную времени для тел любой формы можно записать в общем виде

(10)
где V и F — объём и площадь поверхности теплообмена тела, RТ – суммарное термическое сопротивление переносу тепла от тела к теплоносителю, включающее среднее значение внутреннего термического сопротивления тела.
Постоянная времени реального цилиндрического твэла, имеющего оболочку и газовый зазор между топливом и оболочкой, учитывает их термические сопротивления. Приближенное выражение имеет вид:
(11)

Слайд 60

Рисунок 7 – Распределение
температуры в твэле с газовым зазором между топливом

и оболочкой

Слайд 61

Рассмотрим теперь задачу охлаждения разогретой без учёта ее внутреннего тепловыделения.
Пусть АЗ

охлаждается жидкостью с расходом M и постоянной температурой t. Начальная температура топлива t0 , теплоёмкость топлива CT=VTρTcT; теплоёмкость теплоносителя Cж=Vρc=SHρc.
Здесь VT ‒ объём топлива; S, H ‒ площадь проходного сечения и высота ТВС, соответственно.
Из уравнения теплового баланса можно найти изменение температуры за интервал времени dτ:
‒(CT + Cж) dt = Mctdτ , (12)
Заметим, что Mc=wSρc=(wC)/H, поэтому выражение (12) может быть переписано в виде
(12а)
где H/w=τтр ‒ время транспорта теплоносителя через активную зону.
Выполнив интегрирование, получим
(13)

Слайд 62

В качестве примера рассмотрим характерные времена запаздывания для ВВЭР-440.
Основные параметры:
объём топлива

VT=4.5 м3;
радиус топливной таблетки R=3.775 мм;
высота топливного столба Н=2.5 м;
объём теплоносителя V=7.56 м3;
скорость теплоносителя w=3.4 м/с;
коэффициент теплообмена αср=2260 Вт/(м2·К);
свойства топлива:
ρТ=1·104 кг/м3; сТ=330 Дж/(кг·К); λТ=2.3 Вт/(м·К);
свойства теплоносителя: ρ=750 кг/м3; с=5460 Дж/(кг·К).

Слайд 63

Время запаздывания ‒ τвн ‒ вследствие внутреннего термического сопротивления твэла можно

оценить так:
Время запаздывания, обусловленное наружным термическим сопротивлением (граница раздела «топливо‒теплоноситель»), можно оценить следующим образом:
Общее время запаздывания

Слайд 64

Время прохода теплоносителя через активную зону:
τтр = H/w = 2.5/3.4 ≈

0.735 с.
Теплоёмкость топлива:
CT = cT ρT VT = 330·104·4.5 = 1.485·107 Дж/К.
Теплоёмкость теплоносителя:
C = cρV = 5460·750·7.56 ≈ 3.096·107 Дж/К.
Отношение теплоёмкостей равно
(C+CT )/C = 1+CT /C = 1+1.485·107/3.096·107 ≈ 1.48.
С учётом отношения теплоёмкостей время запаздывания получим такую оценку времени запаздывания:
τз = τтр·(C+CT)/C = 0.735·1.48 ≈ 1.09 с.

Слайд 65

ТЕПЛОВЫЕ УДАРЫ И
ПЕРИОДИЧЕСКИЕ КОЛЕБАНИЯ ТЕМПЕРАТУРЫ
В переходных режимах температура в разных частях

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

Слайд 66

66

Слайд 67

При быстром снижении температуры теплоносителя внутренние слои трубы охлаждаются быстрее, чем

наружные. Распределение температуры приобретает вид, показанный на рисунке 8.
Если средняя температура стенки , то слои, имеющие температуру выше средней, оказываются сжатыми, а ниже средней растянутыми.
Значение напряжений оценивается по формуле

(14)

где β ‒ коэффициент линейного расширения;
Е ‒ модуль упругости;
μ ‒ коэффициент Пуассона;
‒ средняя температура стенки:

Слайд 68

(15)
Максимальные напряжения возникают на поверхностях конструкционных элементов.
68

Слайд 69

Под термином "тепловой удар" понимается однократное изменение температуры поверхности.
Однако в потоке

теплоносителя существуют колебания температуры, которые необходимо учитывать, особенно при больших тепловых потоках.
Колебания температуры (пульсации) оказывают влияние на коррозионные процессы, поскольку приводят к разрушению защитных пленок на поверхности и могут вызвать усталостные напряжения.
Причинами пульсаций температуры могут быть: ●колебания мощности; ●турбулентные пульсации, которые определяются внутренней структурой потока; ●неустановившиеся конвективные течения, особенно в патрубках, тройниках и т.д. при наличии значительных градиентов температуры; ●смена режимов кипения (в месте кризиса, при переходе к пленочному кипению).

Слайд 70

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

h, на одной из поверхностей которой задают пульсации температуры, на другой – теплообмен с греющим теплоносителем при постоянном коэффициенте теплообмена.
Термоупругие нестационарные напряжения на поверхности можно определить по известным соотношениям из теории упругости.

Слайд 71

ОСОБЕННОСТИ ТЕМПЕРАТУРНЫХ РЕЖИМОВ ПРИ РАССЛОЕНИИ ТЕПЛОНОСИТЕЛЯ
В переходных, а особенно в аварийных

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

Слайд 72

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

может вызвать чрезмерные термические напряжения в конструкциях.
Расслоение теплоносителя может происходить в верхней камере у реакторов различного типа (ВВЭР, БН) в режимах с уменьшением энерговыделения.
Нагретые слои оказываются в верхней части камеры, более холодные – в нижней.
Устойчивая стратификация препятствует развитию естественной конвекции.
Между горячей и холодной областями имеется тонкий горизонтальный слой, в котором реализуется весь перепад температуры между верхней и нижней практически изотермическими областями (рисунок 9).
Это обстоятельство необходимо учитывать при установке датчиков измерения температуры, иначе возможны грубые ошибки в оценке её уровней в переходных режимах.

Слайд 73

Рисунок 9 – Расслоение теплоносителя в верхней камере реактора:
1 – активная

зона;
2 – уровень теплоносителя;
3 – граница горячей и холодной зон.

Слайд 74

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

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

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

(16)

Слайд 75

Числами подобия, определяющими характер течений являются также числа Фруда и Рейнольдса:
Другим

примером расслоенных течений является режим подачи холодной воды в трубопровод одной из петли первого контура ВВЭР в аварийном режиме при утечке теплоносителя.
Структура течения теплоносителя в этом случае состоит из двух встречных потоков (рисунок 10).
Горячая жидкость распространяется на длину L в одну сторону, холодная течёт в другую.

Слайд 76

76

Слайд 77

Вопросы, выносимые на зачёт
1. Характер распределения температуры в цилиндрическом твэле
с

зазором между топливом и оболочкой (изобразить графически).

ТЕПЛОМАССОПЕРЕНОС В ЯДЕРНО-ЭНЕРГЕТИЧЕСКИХ УСТАНОВКАХ

Содержание

ТЕПЛОМАССОПЕРЕНОС В ЯДЕРНО-ЭНЕРГЕТИЧЕСКИХ УСТАНОВКАХТема №12Особенности процессовгидродинамики и теплообменав активных зонах реакторов

Методытеплофизического расчётастержневых ТВС03

04

05

Термогидродинамический расчёт ТВС проводится с целью определения параметров потока теплоносителя, охлаждающего

07

08

Рисунок 3 – Поперечные сечения сборок стержней с квадратной (А) и

10

11

12

13

Исходная (а) и трансформированная (рассчитываемая) (б) области в методе согласованного с

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

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

Таким образом, субканал – часть занятого теплоносителем пространства сборки, ограниченная (рисунки

Рисунок 4 –К определению термина "меж-твэльная ячейка»18

Рисунок 5а – Фрагмент из участков трёх твэлов, расположенных в гексагональной

Рисунок 5б – Фрагмент из участков трёх твэлов, расположенных в гексагональной

Рисунок 6а – Центральная, боковая и угловая межтвэльные ячейки21

22

Рисунок 6б – Центральная, боковая и угловая межтвэльные ячейки23

Такой подход вполне корректен, если ограничиваться ● расчётом параметров потока, усредненных по