Партнерка на США и Канаду по недвижимости, выплаты в крипто

30% recurring commission
Выплаты в USDT
Вывод каждую неделю
Комиссия до 5 лет за каждого referral

;

; ;

;

где: Pi - давление на i-ом участке (Р1= Рбс ), кгс/м2 ;

DPМ - местная потеря давления ü

DPТ - потеря давления на трение ý на i-ом участке, кгс/м2;

DPН - потеря за счет перепада высот þ

z - коэффициент местного сопротивления;

yМ; yТ - поправочные коэффициенты для расчета местных потерь и потерь на трение;

G - расход через канал, кг/час;

Si - проходное сечение участка, м2;

g² ; g' - удельный вес пара и воды, кг/м3 ü на линии насыщения;

u'', u' - удельный объем пара и воды, м3/кг ý при давлении Рi ;

i'' , i' - энтальпия пара и воды, ккал/кг þ

g = 9.81 м/с2 – ускорение свободного падения;

R - коэффициент двухфазного потока;

l - коэффициент трения;

li , di , hi - длина участка, гидравлический диаметр и перепад высот, м;

g - удельный вес двухфазной смеси, кг/м3;

j, b, X – истинное объемное, расходное объемное (без учета скольжения фаз) и весовое паросодержание двухфазной смеси;

- средняя шероховатость канала на участке;

Fr, Rе - критерии Фруда и Рейнольдса;

w0 - скорость циркуляции, м/с;

n' [м2/с]; m¢ [кг×с /м2] - кинетическая и динамическая вязкость воды на линии насыщения, (находится интерполяцией таблицы m'T (p) );

КПР - коэффициент проскальзывания пара;

РКР - критическое давление = 225.65·104 кг/м2;

i - энтальпия теплоносителя на участке, ккал/кг;

НЕ нашли? Не то? Что вы ищете?

WК – мощность канала, кВт; 860 ккал/(кВт·ч) – переводной коэффициент.

Т.4.3 Для обогреваемого участка (активная зона) энтальпия рассчитывается пропорционально интегралу от энерговыделения по высоте:

;

где - высотное распределение энерговыделения по данной половине реактора (интерполяция по ВРД-В).

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

;

где:

qi - тепловой поток на участке, ккал/ч. м2;

Fтепл - поверхность теплообмена ТВС, м2;

Кяч – коэффициент неравномерности тепловыделения по сечению ТВС;

hгр – относительная доля тепла, уходящего в графит;

- теплота парообразования, ккал/кг;

Фср – среднее по высоте значение Ф(z).

Т.4.3.1 В случае, если на данном участке i > iнк, рассчитывается энтальпия жидкой фазы

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

после чего применяются данные выше зависимости для двухфазного потока, за исключением:

;

где: c - поправка на неподобие пучка и трубы;

e - доля сечения канала, занятая твэлами и несущей трубой ;

Sтрубы, Sпу - площади поверхности трубы (канала) и пучка твэлов, м2 ;

d, D - диаметры твэла и трубы, м.

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

;

где:

- условный удельный объем.

Т.4.3.2 В случае, если на данном участке i < iнк , расчет участка ведется по формулам однофазного потока, где по сравнению с Т.4.2 :

; R = 1 ; ;

j = 0 ; b = 0 ; Fr = 0 ;

m - динамическая вязкость недогретой воды (находится интерполяцией таблицы mT(i)), кг·с/м2 ;

u - удельный объем недогретой воды (находится интерполяцией таблицы uТ(i)), м3 /кг.

Т.5 Теплогидравлический расчет по п. Т.4 производится для условных каналов левой и правой половины реактора (с усредненными геометрическими характеристиками), в исходном и конечном (возмущенном) состояниях. В результате расчетов вычисляются распределения объемного паросодержания по высоте:

- в левой половине, в исходном состоянии;

- в правой половине, в исходном состоянии;

- в левой половине, в конечном состоянии;

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

Далее вычисляются высотные функции изменения паросодержания в опыте по половинам реактора:

Т.6 Ввиду неравномерности изменения паросодержания по высоте реактора, вклад каждого высотного участка в изменение реактивности, согласно [17], составит

где Y (z) - текущая форм-функция;

- исходная ценность нейтронов.

Интегрируя по высоте и производя общепринятую замену:

получаем:

Коэффициент пропорциональности aj, описывающий свойства всего реактора, должен стоять за знаком интеграла, тогда значение парового коэффициента реактивности равно:

где знаменатель рассчитывается по формуле :

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

Т.7 Изменение реактивности в опыте рассчитывается по формуле

где rв - «вес» взвешенного по реактиметру участка стержней АР;

z1, z2 - исходное и конечное средние положения стержней АР в опыте с возмущением Gпв , м;

z3, z4 - то же, в опыте по «взвешиванию», м;

- градуировочная кривая, рассчитанная по сигналам ВРД-В;

DW - изменение мощности в зоне нечувствительности АР, МВт;

az - допплеровская составляющая быстрого мощностного коэффициента реактивности, b/МВт (оценка по данным измерений на Nном):

значение может быть уточнено по результатам расчётов.

Учет составляющей - az DW обязателен, так как при малых aj отработка АР (z1 ® z 2 ) может быть слабой или отсутствовать вовсе.

Расчет aj в соответствии с формулами Т.2 - Т.7 реализован в программе ТРАКТ [16,22].

Т.8 Определение «веса» взвешенного участка АР и быстрого мощностного коэффициента реактивности aw производится путем обработки сигнала изменения нейтронной мощности в процессе её самостабилизации после ввода отрицательной реактивности стержнями АР на участке, близком к отрабатывавшим возмущения по питательной воде. (Методикой предусмотрены и отдельные опыты по определению aw, см. приложение Р).

В реакторах РБМК структура мощностного эффекта реактивности достаточно сложна. При aj < 1b мощностной эффект обусловлен в большей степени изменением температуры топлива (допплер-эффект). При "взвешивании" стержней АР переходной процесс в этом случае носит устойчивый характер, и значение aw определяется из соотношения:

где DW1- установившееся значение отклонения нейтронной мощности от исходной (через 40-60 с после возмущения).

Ввиду искажения действительного «веса» участка АР действием отрицательного мощностного эффекта реактивности за время ввода реактивности стержнями , в измеренное реактиметром значение «веса» стержней должна вводиться поправка:

Здесь: DWТ - изменение нейтронной мощности за время ;

- относительное изменение тепловой мощности по отношению к нейтронной за время , зависящее от соотношения времени и постоянной времени твэлов t при линейном законе изменения нейтронной мощности. Экспериментальное значение t = 8 с (для мощностей, близких к номинальной).

Расчет rв , aw в соответствии с формулами Т.8 реализован в программе РЕЛЬС [21, 22].

Т.9 Итоговая величина aj по серии опытов рассчитывается как средняя взвешенная, с весами, пропорциональными суммарным возмущениям GПВ в каждом опыте, по фактически зарегистрированным значениям:

;

где j – номер опыта.

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

Т.10 Оценка погрешности определения aj по данным измерений

Т.10.1. Для детального учета всех составляющих погрешности определения aj следует, помимо знания источников ошибки, иметь совокупность весовых коэффициентов, с которыми эти источники влияют на результат.

Тогда оценка погрешности aj по данным измерений выразится в виде:

где:

d (aj ) - погрешность определения aj ;

d (DGпв ), drв - погрешности измерения исходных данных;

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

Для определения вклада в погрешность aj различных составляющих в исходные данные одного из опытов по определению aj (3 блок ЧАЭС, 25.12.1987 г., 81 ДП, aj ~ 0.9 b, DGпв ~ 400 т/ч) вводились вариации, соответствующие абсолютным погрешностям измерения технологических параметров в соответствии с оценками [17]. Результаты расчетов представлены ниже в виде таблицы Т.1, где буквенные обозначения – те же, что и в приведенных выше алгоритмах.

1- Составляющие погрешности определения aj по данным измерений

Группа параметров	Параметр	Абсолютная погрешность	di (aj) , b
Исходного состояния		212 м3/ч	0,002
	33 т/ч	0,001
	1,55 кгс/см2	-0,006
	°C	-0,004
	кгс/см2	0,001
	°С	0,005
	40 МВт (т)	0.007
при скачке Gпв		33 т/ч	-0,030
	0,41 кгс/см2	0,002
	3 МВт (т)	0,011
z1 –z2	0,05 м	0,083
При “взвешивании” АР	rв	с	0,021
z3 –z4	м	-0,060
aw	0,4·10-4 b/МВт (т)	0,051

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

Наибольшее влияние на погрешность aj оказывают погрешности измеренных в опыте возмущений ( DGпв) и их отработки (z1 – z2 , DW), а также параметры "взвешивания" АР.

При повторении циклов возмущения Gпв n раз составляющая погрешности за счет параметров второй группы уменьшается пропорционально 1/ в то время как составляющая за счет параметров третьей группы остается неизменной.

При проведении трех циклов возмущения Gпв (+200/-400 т) оценка приборной погрешности определения средней величины aj составляет ± 0.1 b, что составляет 11% для aj =0.9 b.

Приведенные расчетно-аналитические оценки следует считать минимальной границей погрешности определения aj. Практически дополнительные источники ошибок (субъективные ошибки считывания, неполная стационарность состояний, неучтенный гистерезис в показаниях приборов и пр.) приводят к увеличению фактического разброса значений aj, получаемых в опытах. Поэтому для средневзвешенного значения aj, рассчитанного по формуле из Т.9, следует рассчитывать опытно-статистическую оценку погрешности:

, , ,

где:

D () – погрешность средневзвешенного значения aj при доверительной вероятности e, b ;

- среднеквадратичная ошибка , b;

t (e ,n) – коэффициент Стьюдента, при доверительной вероятности e = 0.95:

t = 3.2 при n = 4, t = 2.6 при n = 6;

n - число опытов;

j - порядковый номер опыта;

- значение aj в j-ом опыте;

В протокол результатов измерении aj включается средневзвешенное значение aj с опытно-статической оценкой погрешности, если она равна или превосходит минимальную расчетно-аналитическую оценку, в противном случае приводится минимальная расчетно-аналитическая оценка ± 0.1 bэф.

Т.11 Программы РЕЛЬС, ТРАКТ, реализующие алгоритмы Т.1 – Т.10, включены в состав прикладного программного обеспечения оперативной рабочей станции физических измерений (ОРС-Ф), входящей в состав ВИ СКУЗ. Эти программы могут вызываться из меню пользователя как самостоятельно, так и под управлением головной программы alfafimetr. Полное описание работы с этими программами содержится в [22].

Перед запуском программы в директории /home/alfafi должны находиться файлы регистрации параметров с расширениями *.rgg, *.rg0, созданные программой orsf_reg в ходе опытов по возмущениям Gпв и "взвешиванию" участка стержней АР (см. приложения В, П). По умолчанию, последние файлы регистрации (скажем, .rgg, .rg0) дублированы с именами alfafi.rgg, alfafi.rg0.

Приложение У

(обязательное)

Алгоритмы обработки данных для определения величины парового и мощностного коэффициентов реактивности

(метод наименьших квадратов)

У.1 Расчетные зависимости

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

Баланс реактивности в этих процессах описывается единого вида уравнениями, аналогичными формулам пп. Т.5 - Т.6:

az · DWl + aj · Dj l + rсуз, l = rк, l - rн, l ,

где l = 1....L - номер уравнения баланса реактивности,

L - количество уравнений баланса реактивности реактора (по числу опытов, включая опыт (опыты) по "взвешиванию участка стержней АР;

(dT / dW) - производная температуры топлива по мощности; aТ - коэффициент реактивности по температуре топлива, (bэф /°С);

aj - паровой коэффициент реактивности, (bэф / %) ;

rсуз, l - реактивность, введенная стержнями АР в l - ом опыте, (bэф);

Dj l - изменение среднего по реактору объёмного паросодержания, взвешенное с квадратом высотного поля для l - ого опыта, % (рассчитывается согласно алгоритмам пп. Т.3 - Т.6 , с использованием программы ТРАКТ;

DWl - изменение мощности реактора для l - ого опыта, МВт(т);

rк, l , rн, l - конечная и начальная реактивности реактора для l - ого опыта, рассчитанные по программе РЕЛЬС (bэф) .

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

У.2 Аналитические МНК - оценки дисперсии параметров

В матричном виде уравнения баланса реактивности для событий с изменением и перемещением стержней АР (ЛАР) записываются как

= , (У.1)

МНК - оценки вектора параметров определяется из решения соответствующей системы нормальных уравнений

= (У.2)

В явном виде (на примере одной из форм балансных уравнений реактивности для событий с изменением и перемещением стержней АР (ЛАР)) (У.1) записывается как