Файл: Козырев, А. П. Теория тепловых и гидродинамических процессов в атомных энергетических установках учеб. пособие.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 15.10.2024

Просмотров: 190

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

тем, что теоретический анализ проводился без учета ак­ сиальной теплопроводности вдоль потока жидкого металла и наличия посторонних примесей в опытах. При Р е = 30 300 опытные данные обобщаются простой по струк­

туре зависимостью

[ 15 ]

Ни = 4 , 3 6 + 0,016Ре .

( б . 7 б )

А/и

Р е

Рис. 6.24. Теплоотдача при течении жидких металло» в трубах:

I - А/и = 7+0,025Я ?°»8j 2 - А/и = 5 + + 0,025 Ре"*8

При уменьшении числа Ре зависимость (6.76) дает пе­ реход допредельному случаю теплообмена при ламинарном

режиме ( N u = 4 ,3 6 ).

В турбулентном режиме ( /?е>1СП) при числах Ре = = 300fI0 ООО теплоотдачу можно рассчитывать по эмпири­

ческой зависимоети[15J

Ни = 1,5 + 0,005Ре .

(6.77)

у тяжелых и щелочных металлов физические параметры слабо зависят от температуры. Кроме того, благодаря высокой интенсивности теплообмена температурный напор стенка-жидкость обычно очень мал. Поэтому расчетные зависимости не включают поправку, учитывающую измене­ ние физических параметров по сечению канала.

253

Влияние примесей на

теплоотдачу

Зависимости (б.7б) и (б .77)

получены при условиях,

обеспечивающих малое содержание неметаллических приме­ сей, однозначно характеризующих чистоту жидкого металла. Однако опыт эксплуатации контуров [ 7 ] , особенно при высоких рабочих температурах, показывает, что в жидких металлах содержатся металлические и неметалли­ ческие примеси, существенно влияющие на теплообмен. Металлические примеси ( Fe , Сх , Ni и т .д .) по­ являются в результате коррозии контура при частичном растворении стенок и вымывания отдельных фракций из конструкционных материалов.

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

Присутствие в потоке расплавленных металлов взвешен­ ных частиц окислов металлов приводит к появлению так называемого термического контактного сопротивления и

254


снижению теплоотдачи. Природа термического контактного сопротивления до конца не выяснена. Наиболее вероятной причиной его появления является рост концентрации не­ растворимых частиц в пристенном вязком слое и образо­ вание дополнительной прослойки из примесей и окислов на границе раздела жидкий металл-стенка.

Термическое контактное сопротивление является ре­ зультатом сложных физико-химических и теплогидродина­ мических процессов на поверхности теплообмена. Нали­ чие контактного сопротивления проявляется в снижении коэффициента теплоотдачи по сравнению с теоретическим значением, а также в нестабильности теплообмена во вре­ мени и усиленной пульсации температуры стенки. При сильно загрязненном теплоносителе коэффициенты тепло­ отдачи могут снизиться в 3 - 5 раз. При таком резком ухудшении теплоотдачи невозможна нормальная эксплуа­ тация установки и требуется специальная очистка тепло­ носителя.

Нижний уровень теплоотдачи на загрязненном приме­

сями

металле при

Ре = IOO ~ 20 000 можно опреде­

лить

по формуле

[15 ]

 

 

Ни =

4,36 + 0,00SiPe .

(6.78)

Учет влияния окислов в потоке на теплоотдачу к на­ трию при турбулентном движении в вертикальном обогре­ ваемом канале приближенно можно определить по зави­

симости

[61 ]

 

 

 

N u

(6.79)

 

Ни,

 

 

t + fNu

 

где Ни

-

число Нуссельта для чистого

металла#

Nuj

-

число Нуссельта при течении загрязненного

 

 

 

255


металла}

-поправка на наличие термического контакт­ ного сопротивления,

(п

0 , 5 ’ М Ъ Ъ +16-ЮгС01?- 30

____

/

 

0.7

; (6*80у

 

 

R e 0’7

 

С0& = 0,02 т

0 , I# -

весовое содержание

кислорода в

жидком металле.

 

 

 

С ростом числа Re

и уменьшением содержания кисло­

рода в теплоносителе термическое контактное сопротивле­ ние снижается, теплоотдача растет. При изменении весо­ вого содержания кислорода в натрии от 0 , 1 до 0 , 0 2 # при

Ре

=

200

теплоотдача увеличивается примерно на 30#.

П ри/?е>2 •

Ю5 поправка становится малой. Влияние числа

Re

 

обусловлено тем,

что

с ростом скорости потока

уменьшается

толщина ламинарного

подслоя и большее

количество окислов уносится из пристенного слоя.

 

Влияние

неравномерного

тепловыделения

по

 

 

 

 

 

высоте

канала

 

 

Существенная разница между теплоотдачей при

=

const и tCT = const

, как это

видно

из формул

(б .74)

и (6 .75),

свидетельствует о сильном влиянии гра­

ничных условий на теплоотдачу.

Это влияние

наблюдается

только

при низких

числах

Я г

• При

Я г >

0,7 и

R e

>

10 ^

влияние указанных граничных условий состав­

ляет

всего

несколько процентов,

и ими можно пренебре­

гать .

Причина

заключается в сильной зависимости профи­

ля температур

от

числа

р г

.

При высоких

числах Я г

термическое

 

сопротивление сосредоточено вблизи стен­

ки, а профили температур не

чувствительны к аксиальным

изменениям

температуры и плотности теплового потока на

256


стейке. При низких числах р х высокая теплопровод­ ность жидкости делает форму профиля температуры, рас­ пределение термического сопротивления, а следователь­ но, и теплоотдачу чувствительными к граничным условиям.

Указанные выводы имеют важное практическое значение, в частности при расчете системы охлаждения и темпера­ турного состояния активных зон реакторов. Если в каче­ стве теплоносителя используется газ или вода под дав­ лением ( Р г ^ 0 ,7 ), то аксиальное изменение плотности теплового потока по высоте активной зоны слабо влияет на теплоотдачу. Б этом случае для расчета местного тем­ пературного напора можно пользоваться традиционным ме­ тодом, связанным с применением коэффициента теплоотда­ чи Ы. по местной плотности теплового потока в дан­ ном сечении. Однако при использовании в качестве теп­ лоносителя жидкого металла такой метод может дать су­ щественную ошибку в определении температуры оболочки тепловыделяющего элемента. Причина ошибки состоит в том, что при переменном тепловыделении ( ^ ф const ) распределение температур в данном поперечном сечении канала будет отличаться от условий при ^ = const

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

я реакторных системах с жидкометаллическим тепло­ носителем влияние локальных изменений теплового потока

на

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

лежащих ниже по пото­

ку,

проявляется особенно сильно,

из-за малого вклада

турбулентного перемешивания в общий процесс переноса

тепла. ,Глина затухания

локальных возмущений темпера-

17, зак. 7д

257

туры соизмерима с длиной тепловыделяющего канала. Для учета этого физического факта при определении темпера­ тур в зоне можно использовать метод суперпозиции [48.J,

основанный на принципе независимого действия источника тепла на температурном поле. Сущность принципа в следу­ ющем. Так как дифференциальное уравнение энергии линей­ ное и однородное, то сумма его ранений также является решением, произвольное тепловыделение по высоте канала можно аппроксимировать ступенчатым и определить общее решение путем суммирования (суперпозиции) решений для каждой ступени. Применение принципа суперпозиции приво­ дит к вычислению интеграла Дюамеля

V

’ > ■

%

d *. (6.81)

 

 

 

 

где

t* (%)

- температура, определенная при постоянном

 

 

равном ^

. Способы вычисления ин­

 

 

тегралов типа Х б.81) приводятся

в работе

На рис.

В.Е. Минашина и др .[4 8 ].

опре­

6.25 приведены температуры стенки,

деленные различными методами при синусоидальном тепло­ выделении по высоте ка­ нала. При неравномерном

 

тепловыделении, близком

 

по своему характеру к

 

тепловыделению в активной

 

зоне реактора,температура

 

стенки больше и макси­

 

мум смещен. В подобных

 

задачах проще использо­

 

вать коэффициент

тепло­

 

отдачи

ос t но

только

 

до стадии эскизного про­

_

ектирования. При принятии

технических решений ин-

Рис.6.25. К расчету теплообмениматьЛ°какНвпипр?°ПопяК

на при синусоидальном тепло-

 

HepaB-

выделении по высоте канала:

нияеия°£1ипопа^0выделе-

I - решение методом суперпо-

??ояние

зоныРи вРкяки5 °~

зиции; 2 - приближенное ре-

стояние

зоны и в каких

шение

условиях

какая допускает­

ся ошибка при использовании традиционных методов расче­ та теплоотдачи.

258


Влияние входных у с л о в и й

Приведенные выше зависимости (б .74) - (6.77) полу­ чены в условиях стабилизированного теплообмена жидких металлов в длинных каналах с ^ / d > 30. Однако даже при условии £ = const процесс стабилизации профиля температур в теплообменных аппаратах происходит на оп­ ределенном начальном участке тепловой стабилизации, где коэффициент теплоотдачи меняется по длине. В этом случае расчет локальных температур поверхности нагре­ ва по формулам стабилизированного теплообмена может

привести к ошибочным результатам. Теплоотдачу на вход­

ном участке и в коротких трубах

с

< 30

приближен­

но можно определять по формулам

стабилизированной теп­

лоотдачи, вводя поправочный множитель

в ^

, учитыва­

ющий влияние входных условий:

 

 

 

6

 

 

(6.82)

Следует отметить, что в общем случае термический начальный участок может иметь значительную длину и оп­ ределяться не только проходным сечением канала, но и другими факторами. Например, в реакторе имеют место по­ перечные потоки тепла, вызванные различной мощностью соседних тепловыделяющих сборок, неоднородностью раз­ множающей среды, присутствием регулирующих органов и т .д . В отдельных случаях большая часть зоны может на­

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

259