Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Боланд Дж. Приборы контроля ядерных реакторов (внутризонные)

.pdf
Скачиваний:
22
Добавлен:
24.10.2023
Размер:
10.37 Mб
Скачать

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

Рис. 7.9. Схема ультразвукового расходомера, предназначенного для работы при высоких температурах:

1 — источник звука

с частотой 50

кгц;

2 п е р е д а ю щ а я

линия

(волновод):

3 — устройство д л я

ввода звука;

4

чувствительный

элемент

приемника

 

 

сигнала.

 

 

7.7. РАСХОДОМЕРЫ, ПРИНЦИП ДЕЙСТВИЯ КОТОРЫХ ОСНОВАН НА СТАТИСТИЧЕСКИХ МЕТОДАХ

Статистические методы, первоначально созданные для полу­ чения информации только из сигналов шумов электрических цепей, достаточно обсуждались в литературе и в настоящее время нашли применение во многих областях научных исследо­ ваний. Основные положения теории шумов базируются на сложных математических выкладках, которые можно найти в соответствующей литературе по статистической теории автома­ тического регулирования и теории случайных процессов [22, 23]. Более краткий обзор основных положений теории шумов можно найти в монографии Тая [24].

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

110

зования флуктуации температуры

[25, 26] и электропровод­

ности [27] для измерения расхода

теплоносителя в ядерных

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

Рис. 7.10. Схема расходомера для измерения расхода ме­ тодом анализа временных характеристик переноса тепло­ вых волн:

/ — усилитель; 2 — анализатор; 3 — элементарный объем тепло­ носителя.

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

Если

две

термопары поместить

на расстоянии d друг от

друга по

оси

кассеты (рис. 7.10), то

возмущение температуры

в малом единичном объеме около термопары 1, расположенной выше по направлению движения жидкости, достигнет термо­ пары 2 через некоторое время, если, конечно, траектория эле­ ментарного объема V пройдет через нее. Запаздывание времени в определенных условиях будет прямо пропорционально сред­ нему времени переноса жидкости между термопарами. Тогда объемный расход через кассету можно представить в виде следующего выражения:

где Q.— объемный расход; V — элементарный объем теплоно­ сителя, который движется от термопары 1 к термопаре 2; At — время движения элементарного объема между термопарами; х — функция, которая связывает средний расход в кассете со скоростью движения элементарного объема.

111

Периодические сигналы. Сначала предлагается рассмотреть случай, когда принудительные колебания температуры жидко­ сти около среднего значения вызываются осциллирующим -источником тепла, расположенным выше термопары /. Извест­ но, что периодический сигнал можно разложить в ряд Фурье

у(0 =

^ е х р ( - М 1 ) ,

 

 

(7.11)

в котором амплитуда Фурье при частоте

равна

 

 

 

Р/2

 

 

 

 

 

= у

j

y{t)^(^f-)dt.

 

 

(7.12)

 

-р/2

 

 

 

 

 

Действительные и мнимые части амплитуды Yj можно

записать так:

 

 

 

 

 

 

 

 

Р/2

y(t) cos 2?iL dt;

 

 

Re У - = —

I

 

(7.13)

 

P

J

P

 

 

 

 

-p/2

 

 

 

 

Im Yj

p/2

y(t)sm-=^-dt.

 

-

(7.14)

p

J

 

 

P

 

 

 

 

-p/2

 

 

 

 

Величину и фазу /-го члена

разложения можно

получить по

•формулам

 

 

 

 

 

 

\Yj\ = ( R e F / r - + ( I m F / ; )

 

 

Ф =

arctg Im К,- Re К;-.

J

 

 

Если показания термопар

7"i и 7"2 рассматривать

как сигналы

j / i и у2 с одинаковым началом

отсчета, то

разность

фаз

(в гра­

дусах) для /-го члена разложения сигнала можно в соответ­

ствии с

методом Фурье

вычислить

и

связать со

временем

переноса следующим соотношением:

 

 

 

 

 

t,=

 

p Т 2 ~ ф 1 .

 

 

 

(7.16)

 

1

г

360

 

 

v

'

Если

запаздывание фазы

сигнала

Т2

относительно

сигнала

Гц зависит от скорости теплоносителя, то значение tj не будет зависеть от p/j. Таким образом, чтобы убедиться в том, что расход по формуле (7.10) вычислен правильно, надо рассчи­ тать tj при нескольких частотах возбуждения для определенных значений расхода. Однако у встречающихся на практике систем имеется ряд свойств, из-за которых tj зависит от часто­ ты. К ним относятся следующие:

1) временная характеристика термопар, которая может оказаться гораздо большей, чем наивысшая частота колебаний;

11,2

2)скорость теплопередачи от теплоносителя к стенке кассеты;

3)скорость теплопередачи вдоль кассеты.

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

при известных расходах.

 

 

 

Сторрер [28]

показал,

что фазу и

амплитуду колебаний

температурного

сигнала,

переносимого

потоком

теплоносителя

в необогреваемом канале, можно описать

выражением

 

 

 

 

(7.17)

где У (z) — комплексное

значение температуры

теплоносителя

в точке, расположенной ниже по направлению движения потока па расстоянии z от рассматриваемой точки; У(0) —комплексное

значение температуры в

рассматриваемой точке; v — средняя

скорость теплоносителя; К находится из выражения

 

Х =

т . а + к о т , )

( 7 - 1 8 )

где i = У — 1; т/ = с//Л;

тс = сс//г; с,

и сс — теплоемкости

тепло­

носителя п твэла на единицу длины

твэла; /г — полный

коэффи­

циент теплопередачи от твэла к теплоносителю.

 

Выражение (7.17) переходит в простое уравнение относи­ тельно запаздывания времени при перекосе в том случае, если теплоемкость твэла или коэффициент теплопередачи малы и в пределе стремятся к нулю. Заметим также, что

Я -»• [1 + (VT c)

П Р И ю

о.

Следовательно, для того чтобы

при

измерениях расхода,

уменьшить влияние твэла, необходимо выбрать такую конфигу­

рацию твэла,

которая обеспечила

бы наибольшее

отношение

теплоемкости теплоносителя

к теплоемкости твэла.

 

Сторрер выяснил, что теплопроводность жидкости в формуле

(7.17) не играет роли как

при высоких частотах,

так и при

очень малых скоростях теплоносителя.

 

Бейли (29)

рассчитал

влияние

молекулярной

диффузии

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

S Д ж . В о л а н д

113

нениях температуры теплоносителя от средней на А7" градусов, они также применимы к распространению температурных волн.

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

10°

Ю1

10г

 

Скорость жидкости,

см/сек

Рис. 7.11. Зависимость минимальных расстояний между термопарами от скорости теплоносителя при различных коэффициентах диффузии тепла k в жидкости.

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

Для простоты предположим, что сигналы от термопар 1 и 2 можно изобразить в виде импульсов, как это показано на рис. 7.12. Очевидно, что импульсы на рис. 7.12, б те же самые, что иа рис. 7.12, а, но смещены во времени и каждый второй опущен. С теоретической точки зрения никаких ограничений на размер, форму и расстояние между составляющими темпера­ турных сигналов не накладывается.

Кросскоррелящ-юнную

функцию, связывающую

случайные

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

у\ (t)

и г/2 (0

термопар 1

и 2 можно

записать следующим

образом:

 

 

 

С,,2

(т) =

lim " f

У 1 (t) у2

(t + т) dt.

(7.19)

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

1.14

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

I | 1 г

I '

I

•t=0 /&

t=p-

Qljjрезак бремени, на котором дррврдщся корреляция

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

регистрируются обеими термопарами / и 2). Поскольку нас интересует только измерение времени запаздывания, то относи­ тельная амплитуда температурных сигналов несущественна, поэтому важен факт появления сигналов как больших, так и малых.

Приближенную нормированную кросскорреляционную функ­ цию можно записать в следующем виде:

 

 

 

C l . 2 ( T ) =

Р

 

Ui {t +1) d%

(7.20)

 

 

 

ГI i/i (0

I •

+

Г

 

 

 

 

p

J

 

 

 

 

 

Если импульсы, показанные на рис. 7.12, а,

рассматривать

как сигнал уч, а показанные на рис. 7.12, б как

сигнал г/г, то

очевидно,

что

нормированная

кросскорреляционная функция

С*2

(т) будет

иметь максимум

при т,

равном

t,-. Таким обра­

зом,

кросскорреляционный

анализ

случайных

температурных

сигналов

может дать

теоретически

ту же

информацию, что

метод Фурье для периодических сигналов.

 

 

Экспериментальные результаты. Непрерывный и дискретный способы регистрации данных и кросскорреляционные методы их обработки, если считать указанные данные реакторными шума­ ми, обсуждались в работе Тая [24] с точки зрения относитель-

8*

1 I S

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

Бентлей и Даусон [25] работали с двумя системами такого типа, позволяющими различать пики температурных сигналов величиной 0,0ГС. Рендолл [26, 27] также использовал анало­ говую вычислительную систему для исследования флуктуации электропроводности и температуры потока теплоносителя мето­ дом времени переноса.

Вероятная точность правильно спроектированной системы измерения расхода внутри активной зоны по методу времени переноса пока остается неопределенной, хотя последние экспе­ рименты проводились в весьма неблагоприятных условиях. Из оценок Бентлея и Даусона следует, что при малых ампли­ тудах случайных сигналов погрешность будет не выше 3%. Вполне возможно, что при больших амплитудах периодических сигналов будет достигнута большая точность.

Бентлей и Даусон [25] измерили расход методом времени переноса в петле высокого давления диаметром 12,7 см при давлении 175 атм и температуре 350° С. Термопары диаметром 1,6 мм были изолированы и заключены в оболочку. Они распо­

лагались на входе

и

выходе

модели пучка топливных элемен­

тов на расстоянии

1,7

м друг

от друга. Бентлей и Даусон изме­

рили также расход в натриевой петле диаметром 30,4 см при температуре 400° С. Термопары были подобны использованным в водяной петле, но были запаяны твердым припоем в камеры диаметром 4,75 мм. Расстояние между термопарами составляло 7,55 м. Полученные данные приведены в табл. 7.1.

Приведенные значения

когеренции

рассчитывались ' по

формуле

 

 

 

Скот—' ,_

,t«v=

(7-21)

(Фи Ш-Т1.2Ш)"

 

 

где [ср 1.2 (f) |—абсолютное

значение кросспектральной

плот­

ности мощности процесса в узкой полосе частот около частоты/;

Фи

(/) и фо.2 (f) — автоспектральные

плотности сигналов термо­

пар

У и 2 в узкой

полосе часто около частоты f.

Когерентность

представляет собой

меру корреляций сигналов,

зарегистрированных двумя термопарами. Она была бы равна единице при всех частотах f, если бы сигналы были полиостью коррелированы. К сожалению, для сравнения с результатами других измерений по методу времени переноса получить точные значения измеренного таким образом расхода не удалось, по известно, что они лежат в пределах 10% погрешности. Очевидно,

Т а б л и ц а 7. ] Результаты измерений расходов методом времени переноса

Система

Сигнал, °с

Когеренцня

Время, по которому проводилось усред­ нение (полное время регнстрацнн)

Измеренное время переноса, сек

Измеренный расход, л/мин

Статистическая ошибка (средняя кпадратическая), %

Петля

высокого

0,1

0,05—0,25

15 мин

0,63

692

3

давления

 

 

 

 

 

 

с обычной

 

 

 

 

 

 

водой

 

 

 

 

 

 

Эксперименталь­

0,0

0,05—0,12

6 ч

1,23

23 400

4

ная

петля

 

 

 

 

 

 

с принудитель­

 

 

 

 

 

 

ной

циркуля­

 

 

 

 

 

 

цией

натрия

 

 

 

 

 

 

что условия в этих опытах были весьма далеки от идеальных, потому что расстояние между термопарами было очень боль­ шим, а полученные сигналы»—очень маленькими.

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

пары,

чувствительные

к флуктуациям

температур.

Если

флуктуации

температур,

возникающие в

теплоносителе из-за

теплоотдачи

от

топлива,

слишком

малы

или некоррелированы

п достаточной степени и не обеспечивают

надежные результаты

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

периодически

вносить небольшие изменения мощности

и с малой

частотой

таким

образом,

чтобы

амплитуда

флуктуации

температуры

теплоносителя

все же оказалась

достаточно большой. Термо­

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

скорость потока

теплоносителя оптимальна,

достаточно места

для разнесения

термопар по пространству и

обеспечиваются

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

ходимых для высокоточных

измерений.

 

Результаты

измерений расходов,

полученные

при анализе

флуктуирующих

сигналов

электропроводности

движущейся

воды методом

времени переноса,

показали, что этот метод

117

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

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

 

 

 

 

 

 

 

 

СПИСОК

ЛИТЕРАТУРЫ

 

 

 

 

 

 

 

 

1.

Koch

I. J. е. a. Hazards

Summary Report — Experimental

Breeder

Reactor

 

11 (EBR-II). ANL-5719. Argonne National

Laboratory,

May

1957,

p. 65.

2.

Popper

 

G.

F. Lecture

Notes

 

on

In — Core

Instrumentation

for

Measure­

 

ment

 

 

of

Hydrodynamic

Parameters

in

Water — Cooled

Reactors.

 

ANL-6452, Argonne National

Laboratory, November 1961.

 

 

 

 

 

 

3.

EBWR

 

Project

Group,

EBWR

Test

Reports. ANL-6703,

Argonne

National

 

Laboratory,

January

1964.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4.

Rouhani

Z. Application of the Turbine Type Flowmeters

in

the

Measure­

 

ment

 

of Steam

Quality

and

Void.

In

Proceedings

of

Symposium

on

 

In-Core

Instrumentation,

Oslo.

June

15—19,

1964.

Paper

 

D—6

 

(CONF-640607). Institult for

Atomenergi. Holden,

Norway,

p.

3.

 

 

5.

Schenk

 

K. HBWR Turbine Flow Meters

Developed

Spesially

for

Nuclear

 

Application. In Proceedings of Symposium on

In-Core

Instrumentation,

 

Oslo,

June

15—19,

1964,

Paper

D-7

(CONF-640607), Institult for

Atome­

 

nergi,

Halden, Norway.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6.

Howard

C. L. Development Program on

the

Garigliano

Nuclear

Reactor.

 

USAEC

Report No. GEAP-5190, General Electric Company.

 

 

 

 

 

7.

Oliva

 

D.,

Hench J.

E. Turbine

Flowmeter

Calibration

and

Instrumented

 

Fuel Assembly Pressure Dror>. USAEC

Report

No.

GEAP-4954,

General

 

Electric

Company, October 1965, p. 5.

 

 

 

 

 

 

 

 

 

 

 

8.

Popper

G. F.,

Knox

A.

E.

 

Farel

In — Core

Instrument

Development.

 

USAEC

Report No. ANL-7161,

Argonne National

Laboratory.

 

 

 

 

9.

Seventh

Annual

Report,

O.E.C.D.

Halden

Reactor Project, 1965/66, Report

 

No.

NP-16803,

p. 77,

also

see

K.

Schenk,

Turbine

 

Flowmeters

Tor

 

In-Core

Applications, Report No. HPR-72, p. 15, O.E.C.D. Halden Reactor

 

Project, June 1967.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

10.

Sixth

 

Annual

Repori

O.E.C.D. Halden

Reactor

Project,

1964/65,

Report

 

No.

NP-16230,

Institult

for

Atomenergi,

Halden, Norway,

p.

94.

 

 

 

11.Shercliff J. A. The Theory of Electromagnetic Flow Measurement, Cam­ bridge University Press, N.Y., 1962.

12. Kolin A. An Alternating Field

Induction Flowmeter of High Scnsitivily.

Rev. Sci. Instrum., 1945, vol. 16,

p. 109.

13. Popper G. F., Glass M. C. The Design and Performance of a 1200F Mag-

118

netic

Flowmeter

for

In — Core

Applications

in Sodium

Cooled

Reactors.

I E E E

Trans Nucl. Sci., vol. 14, No. 1, February 1967.

 

 

14. Delwyn D. B. e.

a.

Magnetic

Flowmeter

Colibration,

USAEC

Report

No. IS-1203, Ames Laboratory, August 1966.

 

 

 

15.James W. G. An Induction Flowmeter Design Suitable for Radioactive Liquids. Rev. Scient. lnstrum., 1951, vol. 22, p. 991.

16. Lende

H., Lang W. T. Eddy Current Flowmeter, U . S . Pat., 1948, vol. 2,

435, p.

34.

17.Popper G. F. e. a. Summary Review of Flowmeters Suilable for Measuring Sodium Flow at Temperatures up to 1200 F in the Fast Flux Text Facility (FFTF) . USAEC Report No. ANL-7340, Argonne National Laboratory, December 1967, p. 23.

18.

Evans

P.B.F. e. a. Control

and Instrumentation of

Prototype

Fast

Reac-

 

tor. In

Proceedings

of

Conference

on

Fast

Breeder

Reactors, May

17—19,

 

1966.

British Nuclear Energy Society,

London, p.

754.

 

 

 

 

 

 

J9.

Lion

K. S. Instrumentation

in

Scientific

Research. N. Y.,

McGraw-Hill, 1967.

20. Sawochka S. G. Measurements

of

Forced

Convection

Condensing

Potas-

 

sium

Heat

Transfer. In Proceedings

of

the Conference

on

Application of

 

High

Temperature

Instrumentation

to

Liquid

Metal

Experiments. USAEC

 

Report No. ANL-7100, Argonne National

Laboratory, September

1965, p. 115.

21.

Aeroprojects

Incorporated,

Applications

of

Ultrasonic

Energy:

Ultrasonic

 

Instrumentation

for

Nuclear

Application

Bi — monthly

Progress

Report

 

No.

26.

USAEC

Report

No.

NYO-3622-1, February 1966, p. 5; also see

 

other

progress

reports

issued

under

NYO-3622-series

for

latest develop-

ments.

22.Lee Y. K. Statistical Theory of Communication. N.Y., John Wiley & Sons, Inc., 1960.

23.Bendat J. S. Principles and Applications of Random Noise Theory. John Wiley & Sons, Inc., 1958.

24.Thie J. A. Reactor Noise, N. Y., Rowman and Littlcfield. Inc., 1963.

25.Benlley J. G., Dawson D. G. Fluid Flow Measurements by Transit Time Analysis of Temperature Fluctuations. UKAEA Report No. TRG-1155 (R), 1966, Reactor Engineering Laboratory, Risley.

26.Randall R. L. Flow Measurements Using Noise Analysis Techniques.

 

USAEC

Report No. NAA-SR-Memo 9787,

Atomics

' International,

 

April

1964.

 

 

 

 

 

 

 

27.

Randall

R. L. e. a. Development

of Noise

Analysis Techniques for Measu-

 

ring

Reactor Coolant

Velocities.

USAEC

Report

No. NAA-SR-11193, Ato-

 

mics

International, 1966

 

 

 

 

 

 

2S.

Storrer

F. Temperature Response to Power, Inlet

Coolant

Temperature and

 

Flow

Transients in Solid Fuel Reactors. APDA-132, Atomic

Power

Deve-

 

lopment Associates, Inc., Detroit. Mich., June 1959,

p. 34.

 

 

 

29.

Bailey

R. E. The Effect of the Molecular Diffusion of Heat on Fluid Ve-

 

locity Measurements Accomplished by Thermocouple Cross-Correlation.

 

USAEC

Report No.

NAA-SR-10848, Atomics Internationa],

March

1965.

30.Bentley P. G. e. a. Fluid Flow Measurements, British Patent Application 41034/63; 1963.

Соседние файлы в папке книги из ГПНТБ