Размещено на http://www.allbest.ru/

Предварительные результаты расчета международной тестовой задачи v1000ct-1 c использованием интегральных теплогидравлических программных комплексов

Г.В.Алехин, С.А.Курбаев, М.А.Быков

ОКБ «ГИДРОПРЕСС», Подольск, Россия

Аннотация

Результаты испытаний режима подключения одного главного циркуляционного насоса (ГЦН) к трем работавшим, полученные во время пусконалодочных работ на 6 блоке АЭС «Козлодуй» в Болгарии и приведенные в спецификации OECD/NEA VVER-1000 (V1000 CT-1), были использованы для проведения расчетов по интегральным теплогидравлическим кодам, которые применяются в ОКБ «Гидропресс» и включают модели пространственной кинетики и межпетлевого перемешивания в реакторе.

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

Результаты проведенных расчетов сопоставляются с результатами проведенных испытаний и результатами расчетов по другим теплогидравлическим кодам (DIN3D/ATHLET, RELAP5-3D и др.).

Введение

Для расчета переходных процессов в наилучшем приближении в ОКБ «Гидропресс» в настоящее время применяются программные комплексы (ПК) ТРАП-КС [1]-[2], ПК ДКМ [3] и ПК КОРСАР/ГП [4], включающие модели пространственной кинетики и неполного перемешивания в реакторе.

ПК ТРАП-КС включает ранее разработанный ПК ТРАП-97 (включает программы ДИНАМИКА-97 и ТЕЧЬ-М-97) и программные модули (ПМ) КАМАЗ и КАМЕРА-В2.

КАМАЗ - ПМ для расчета теплогидравлических и нейтронно-физических процессов в активной зоне реактора. Модель теплогидравлических процессов основана на разработанной ранее в комплексе ТРАП-97 модели КАНАЛ-97 и позволяет проводить расчеты активных зон с числом каналов до 500 и количеством расчетных узлов по высоте до 30. Трехмерный расчет нейтронного поля в активной зоне осуществляется программным модулем (ПМ) КАРТА, использующим ПК САПФИР [5] для подготовки констант и расчета стационарного состояния реактора в процессе выгорания (разработаны в НИТИ им. А.П.Александрова). ПМ КАРТА предназначен для расчета нейтронно-физических процессов в активной зоне реакторов типа ВВЭР в трехмерном двухгрупповом диффузионном приближении. Для решения уравнений кинетики используется конечноразностный метод с 1, 6 или 24 расчетными узлами на ТВС в плане.

Для уточненного описания потоков вблизи границ областей с различными нейтронно-физическими свойствами в ПМ КАРТА реализована модификация алгоритма Аскью-Такеды или метод коррекции коэффициентов диффузии (МККД) [6].

КАМЕРА-В2 - расчет процессов перемешивания теплоносителя в напорной и сборной камерах реактора. Методика расчета процессов в камерах реактора основана на уравнениях сохранения массы и энергии теплоносителя. Турбулентный массобмен между расчетными ячейками моделируется введением в уравнение энергии и концентрации борного раствора коэффициента турбулентного массообмена.

ПК ДКМ включает код ДИНАМИКА-97, входящий в состав ПК ТРАП-97, ПМ КАМЕРА-В2 и ПМ МАЗ-3. ПМ МАЗ-3 (разработан ОКБ «Гидропресс») предназначен для исследования теплогидравлической обстановки в активной зоне и реакторе типа ВВЭР в нестационарных режимах с учетом трехмерной пространственной кинетики в покассетном приближении. Для решения уравнений кинетики применяется библиотека коэффициентов аппроксимации кода БИПР-7. Для решения уравнений кинетики используется метод конечных разностей с одной точкой на кассету и итерационный метод последовательной верхней релаксации. Для повышения точности расчета кода применяется МККД [6].

ПК КОРСАР/ГП является развитием ранее разработанной НИТИ им. А.П.Александрова первой базовой версии кода КОРСАР, аттестованной в Ростехнадзоре. Дополнительно в ПК КОРСАР/ГП включен ПМ КАРТА для расчета пространственной кинетики, а также ПМ КАМЕРА-В2 (расчет процессов межпетлевого перемешивания в реакторе) и ПМ ТВЭЛ-2 и ТВЭЛ-3 (расчет термомеханических свойств твэлов). ПМ КАМЕРА-В2, ТВЭЛ-2 и ТВЭЛ-3 разработаны в ОКБ «Гидропресс».

Ранее были проведены тестовые расчеты по этим программным комплексам переходного режима с подключением ГЦН петли к трем работавшим для 6 блока АЭС «Козлодуй» на момент кампании 172 эффективных суток [7].

В стандартной проблеме V1000 CT-1 [8] рассматривается аналогичный режим с подключением ГЦН для 6 блока АЭС «Козлодуй» на мощности 27,46% (824 МВт), но на момент кампании 30,7 эффективных суток. В исходном положении регулирующая 10-я группа органов регулирования находится в положении 36% от низа активной зоны. В ходе переходного процесса ее положение не изменяется.

Определены три задачи для теста V1000CT-1:

· Задача 1 - расчет переходного процесса в реакторной установке с использованием точечной модели нейтронной кинетики.

· Задача 2-совместный трехмерный нейтронно-физический и теплогидравлический расчет переходного процесса в реакторе и активной зоне с использованием заданных граничных условий.

· Задача 3 - расчет переходного процесса в реакторной установке в наилучшем приближении.

Для задачи 3 предлагается рассмотреть основной и дополнительно экстремальный сценарий. В экстремальном сценарии дополнительно рассматривается выброс органа регулирования из регулирующей группы в секторе, находящимся вблизи петли с подключаемым ГЦН. Результаты решения стандартной проблемы V1000CT-1 приведены, в частности, в [9]-[10].

В настоящей работе приведены результаты расчета основного сценария по ПК ТРАП-КС и некоторые результаты расчета, полученные по ПК ДКМ. В тесте необходимо провести сравнение результатов расчета с полученными результатами испытаний.

Дополнительно с целью проверки расчета процессов межпетлевого перемешивания в реакторе проведен расчет по ПК ТРАП-КС режима закрытия отсечного клапана (БЗОК) на одном парогенераторе (стандартная проблема V1000-CT2 (задача 1) [11]. Это испытание проведено на 6 блоке АЭС «Козлодуй» для первой топливной загрузки в начале кампании на мощности 281 МВт (концентрация борной кислоты - 7,2 г/кг; группы 9-10 полностью погружены в активную зону, 8-я группа находится в положении 84% от низа активной зоны, остальные группы полностью извлечены из активной зоны, положение групп регулирования не меняется от времени). После закрытия БЗОК парогенератора температура теплоносителя в изолированной петле значительно возрастает, а в других петлях изменяется незначительно, что приводит к несимметричному распределению температуры теплоносителя на входе в активную зону.

Принятые допущения

В таблицах 1-2 приведены принятые в тексте и на графиках условные обозначения и сокращения, на рис. 1-4 используемые расчетные схемы.

Для расчета режимов закрытия БЗОК ПГ и подключения ГЦН к трем работавшим по ПК ТРАП-КС используется следующее разбиение реакторной установки. Система теплоносителя первого контура представлена четырьмя расчетными циркуляционными петлями, на одной из которых расположен КД.

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

Каждая петля разбивалась на следующее количество расчетных объемов: горячий трубопровод - два расчетных объема; горячий коллектор ПГ - один расчетный объем; трубчатка ПГ - пять расчетных объемов; холодный коллектор ПГ - один расчетный объем; холодный трубопровод ГЦТ - пять расчетных объемов; соединительный трубопровод между КД и ГЦТ - три расчетных объема.

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

Активная зона реактора представлена многоканальной моделью, включающей 163 канала, характеризующих поведение всех ТВС активной зоны и 18 каналов протечек. Нодализационная схема активной зоны приведена на рис. 3. В расчете использовалась расчетная схема с шестью точками на кассету и МККД. Периферийные кассеты окружены рядом шестигранных каналов (всего 54 канала), в которых проводится расчет плотности нейтронного потока в шести точках на один канал в плане.

Подготовка библиотеки нейтронно-физических констант для ТВС и условных каналов отражателя для ПК ТРАП-КС проведена по комплексу программ САПФИР. В качестве узловых значений параметров обратных связей НФХ приняты следующие значения: плотность теплоносителя - 0,7816; 0,7144; 0,6093 и 0,4501 г/см³; температура топлива - 543, 873, 1773 и 2600 К; концентрация бора в теплоносителе - 0; 0,5; 1,20; 2,00 г В/кг. При нахождении текущих значений теплогидравлических параметров между используемыми расчетными узлами значения НФХ интерполируются, а при выходе за пределы экстраполируются при помощи интерполяционного полинома, приведенного в [2]. Учет эффектов гетерогенности на границе с отражателем учитывается неявно при подготовке констант отражателя.

В расчетах моделируются показания комплектов АКНП, в каждом из которых моделируются показания 3 каналов, расположенных равномерно по азимуту в соответствии с проектной документацией. Каждый из этих каналов включает два детектора, расположенных симметрично относительно середины активной зоны. К первому комплекту АКНП относятся каналы 1,10 и 20, ко второму- 5, 15 и 25 (рис.4).

При расчете НКР с использованием ПМ КАМЕРА-В2 (рис. 1) рассматривается опускной и подъемный участок НКР. При этом принято следующее разбиение на расчетные элементы: количество участков разбиения по сечению для опускного участка - 40; количество участков разбиения по высоте опускного участка - 20; количество шестигранных каналов по сечению для подъемного участка напорной камеры - 163; количество участков по высоте в подъемном участке напорной камеры - 10.

Расположение входных патрубков реактора задано на первом участке по ходу движения теплоносителя.

Принято допущение о выравнивании поля аксиальных скоростей на последнем участке опускной части НКР.

На участках опускной части и последнем участке подъемной части НКР на основе уравнений сохранения массы и минимизации потерь давления по замкнутому контуру определяются расходы теплоносителя в поперечном направлении. На остальных участках подъемной части принято допущение о равенстве аксиальных скоростей и отсутствии поперечного течения.

На первом участке СКР выделена область, в которой выделяются шестигранные каналы (их количество соответствует количеству кассет в активной зоне). Принято допущение о выравнивании полей аксиальных скоростей на выходе из этих каналов. В этой области происходит расчет поперечных расходов на основе уравнений сохранения массы и минимизации потерь давления по замкнутому контуру.

При расчете СКР принято следующее разбиение на расчетные элементы:

? количество участков разбиения по сечению - 40;

? количество участков разбиения по высоте - 10.

На участках с первый по десятый проводится расчет поперечных расходов между выделенными секторами из уравнений сохранения массы и минимизации потерь давления по замкнутому контуру. Область в центральной части СКР (рис. 2), в которой принимается идеальное перемешивание теплоносителя, в расчете не учитывается.

ПМ КАМЕРА-В2 использует моделирование турбулентного массообмена в реакторе, основанное на использовании коэффициента турбулентного массообмена, который задается в зависимости от числа Рейнольдса (Re) пропорционально Re^0,875 c параметрическим коэффициентом K_Т, который определяется по результатам сопоставления расчетных и экспериментальных значений параметров на входе в активную зону при испытаниях по межпетлевому перемешиванию на действующих энергоблоках и на экспериментальных стендах. Кроме того, в модели имеется возможность задать угол закрутки потока в опускном участке НКР (б_аз), определяющего дополнительное смещение потоков теплоносителя с различной температурой или концентрацией бора относительно ожидаемого. При проведении расчетов можно найти оптимальные значения Кт и б_аз на основании определения минимальных среднеквадратичных отклонений между расчетом и экспериментом по температуре теплоносителя на входе в активную зону.

Для режима закрытия БЗОК ПГ проведен расчет двух вариантов. В первом случае расчет проведен по ПМ КАМЕРА-В2 с заданными граничными условиями в петлях РУ (температура и расход теплоносителя в соответствии с таблицей 3). В этом случае определяется поле температур только на входе в активную зону. Путем поиска оптимальных значений Кт и б_аз, определяются настроечные параметры ПМ КАМЕРА-В2. Во втором варианте расчет проведен по ПК ТРАП-КС с использованием полноконтурной схемы и полученных ранее по ПМ КАМЕРА-В2 значений Кт и б_аз. Результаты расчета основных параметров РУ сравниваются с результатами испытаний.

При расчете режима с закрытием БЗОК ПГ по ПК ТРАП-КС приняты также следующие предположения:

· исходная мощность реактора, давление первого и второго контура, положение ОР СУЗ приняты в соответствии с данными испытаний;

· ЭГСР турбины работает в режиме поддержания давления в ГПК 5,13 МПа;

· гидравлические характеристики петель в пределах погрешности их определения были скорректированы таким образом, чтобы расчетное значение расхода теплоносителя через активную зону в стационарном состоянии примерно соответствовало данным испытания;

· частота вращения ГЦН 1-4 была скорректирована таким образом, чтобы расчетные значения расходов теплоносителя в петлях примерно соответствовали данным испытания;

· регулятор уровня в КД поддерживает исходное значение уровня в КД;

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

При решении стандартной проблемы V1000-CT1 приняты следующие основные допущения:

· исходная мощность реактора, давление первого и второго контура, положение ОР СУЗ приняты в соответствии с данными испытаний;

· ЭГСР турбины работает в режиме поддержания давления в ГПК 6,05 МПа;

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

· уставки на включение (отключение) ТЭН КД приняты в соответствии с [8];

· расход питательной воды в ПГ1-3 в течение первых 130 с процесса приняты по данным испытаний, а в дальнейшем учитывается работа регуляторов уровня в ПГ режиме поддержания заданного уровня;

· в исходном состоянии значения физических уровней в ПГ1-4 были выставлены в соответствии с данными испытаний путем задания характеристик регуляторов уровня в ПГ

· изменение частоты вращения подключаемого ГЦН принято в соответствии с [8];

· значение настроечного коэффициента турбулентного массообмена (K_Т) для ПМ КАМЕРА-В2 принято по результатам его определения в режиме закрытия БЗОК ПГ, закрутка потока в опускной части НКР в расчете не учитывается;

· регулятор уровня в КД поддерживает исходное значение уровня в КД.

В проведенных расчетах не учитываются отклонения НФХ от ожидаемых. Расчет пространственного поля остаточных тепловыделений проведен в соответствии со стандартом MS ISO 10645-92 без учета погрешности расчета, исходное значение остаточного тепловыделения соответствует стационарному значению для исходной мощности реактора.

В [8] упоминается о том, что перед режимом подключения ГЦН к трем работавшим энергоблок работал на мощности 75% от номинальной, затем мощность была снижена до 21% от номинальной и поддерживалась несколько часов в диапазоне 21-30% от номинальной.

Поскольку точная предыстория изменения мощности и концентрации борной кислоты в [8] не приведена для стандартной проблемы V1000 -CT1 проведен расчет трех вариантов.

Во всех вариантах принято допущение, что отравление реактора соответствует уровню мощности 75% от номинальной. В варианте 1 определялась критическая концентрация борной кислоты, в варианте 2 она принята в соответствии с [8] равной 5,95 г/кг. Изменение отравления в ходе переходного процесса в вариантах 1 и 2 не учитывается.

В варианте 3 проведена оценка влияния нестационарного отравления на ход переходного процесса. В исходном состоянии на уровне мощности 75% от номинальной принято ожидаемое положение регулирующей группы ОР СУЗ ( 80% от низа активной зоны). Принимается стационарное отравление реактора, соответствующее исходному уровню мощности. В начальный период происходит снижение уровня мощности до уровня менее 30% от номинальной путем ввода регулирующей группы до положения 36% и изменения концентрации борной кислоты. Рассматривается период перед подключением ГЦН к трем работавшим в течение 3023 с. По истечении этого времени при мощности реактора 824 МВт было проведено подключение ГЦН. В течение всего рассматриваемого процесса учитывается нестационарное отравление реактора.

Расчет по ПК ДКМ проведен при тех же условиях, что и расчет варианта 1, проведенного по ПК ТРАП-КС.

При проведении расчетов первые 100-200 с просчитывались без возмущений, чтобы убедиться в достижении стабилизации параметров.

Результаты выполненных расчетов и их обсуждение

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

Результаты расчета режима с закрытием БЗОК ПГ с использованием ПК ТРАП-КС и ПМ КАМЕРА-В2 приведены на рис. 5-18.

Сопоставление основных расчетных (по ПК ТРАП-КС) и экспериментальных параметров в исходном и конечном состояниях приведено в таблице 4. Сравнение показывает, что основные сравниваемые параметры в исходном и конечном состояниях находятся в диапазоне погрешности измерений. В отличие от эксперимента в расчете происходит незначительное снижение мощности примерно до уровня 8,9% от номинальной. Исходное значение критической концентрации борной кислоты отличается от полученной в испытании незначительно.

Распределение температуры теплоносителя на входе в ТВС, а также разницы температуры холодной нитки петли 1 и температурами на входе в ТВС (dTk), полученное по ПМ КАМЕРА-В2 приведены на рис. 12-14, 18. По результатам проведенных вариантных расчетов получено значение угла закрутки потоков теплоносителя в НКР, при котором имеет место минимальное отклонение результатов расчета от эксперимента, равное 12є. Это значение принято и при расчете по ПК ТРАП-КС.

Качественно картина распределения полей температур на входе в ТВС, полученная в расчете по ПК ТРАП-КС и ПМ КАМЕРА-В2, совпадает с данными испытания (рис. 13-14). В активной зоне образуется сектор с повышенной температурой, расположенный вблизи петли с закрытым БЗОК ПГ. Имеются переходные области, в которых температура теплоносителя на входе в ТВС соответствует промежуточным значениям между температурой теплоносителя в холодной нитке отсеченной петли и соседних петель. В этих областях наблюдаются максимальные отклонения расчетных и экспериментальных данных. В большинстве ТВС отклонения расчетных и экспериментальных значений входной температуры теплоносителя незначительны.

Результаты проведенного сопоставления показывают, что ПК ТРАП-КС ПМ КАМЕРА-В2 может моделировать рассмотренный процесс межпетлевого перемешивания с приемлемой точностью.

Результаты расчета режима подключения ГЦН к трем работавшим и их сопоставление с данными испытания и результатами расчета по другим кодам приведены на рис. 19-41.

Сопоставление основных параметров в исходном состоянии, приведенное в таблице 5, показывает, что основные теплогидравлические параметры находятся в диапазоне погрешности измерения. При этом ПК ТРАП-КС для варианта 1 дает завышенное значение критической концентрации борной кислоты на 0,5 г/кг по сравнению с данными [8]. Для варианта 3 это отклонение составляет 0,3 г/кг. Это отклонение находится в пределах погрешности определения критической концентрации борной кислоты по ПК ТРАП-КС. Изменение мощности реактора перед подключением ГЦН для варианта 3 приведено на рис. 30.

Критическая концентрация, определенная по ПК ДКМ, составила в исходном состоянии 6,25 г В/кг.

После подключения ГЦН со стороны петли №3 происходит снижение температуры теплоносителя примерно на 3,2єС (рис.26). Из-за наличия отрицательных обратных связей по плотности теплоносителя и температуре топлива происходит несимметричное увеличение мощности реактора от 24 до 42 МВт (рис.19-20) . Наибольшее увеличение мощности имеет место для варианта 2. В варианте 3 после 30 с процесса начинается снижение мощности реактора, связанное с влиянием нестационарного отравления. При расчете по ПК ДКМ приращение мощности составило 30,7 МВт, т.е. немного меньше полученного в варианте 1 по ПК ТРАП-КС (35,5 МВт).

Сопоставление расчетного изменения мощности реактора и ее приращения по ПК ТРАП-КС с результатами расчета по другим кодам (DIN3D/ATHLET, RELAP5-3D, TRAC-PF1/NEM, RELAP5/PARKS) приведено на рис. 19-20. Результаты сопоставления показывают, что приращение мощности для варианта 2 в ПК ТРАП-КС больше, чем по другим кодам, для варианта 1 меньше, чем по коду RELAP5-3D, но больше, чем по другим кодам.

Сопоставление расчетного исходного усредненного по сечению активной зоны аксиального профиля энерговыделения с профилями, полученными по кодам BIPR и TRAC/NEM (рис. 31) показывает на их удовлетворительное согласие для вариантов 1 и 2. Для варианта 3 в нижней части активной зоны достигаются более высокие значения энерговыделения, чем в вариантах 1 и 2.

Картограммы распределения энерговыделения между ТВС (Kq) в различные моменты времени приведены на рис. 36-40. На рис.38-40 приведено сопоставление результатов расчета между ПК ТРАП-КС и DIN3D/ATHLET в момент времени 15 с. Максимальные отклонения составляют для вариантов 1-3 соответственно 3.58, 3.94 и 5.05%.

Распределение полей температур на входе в ТВС в зависимости от времени после подключения ГЦН приведено на рис. 35.

В стационарном состоянии температура теплоносителя в холодной нитке петли №2 меньше, чем в петлях 1 и 4, что связано с влиянием холодного теплоносителя, попадающего через СКР в петлю №2 из горячей нитки петли №3. Расчетные и экспериментальные значения температуры теплоносителя в петлях РУ качественно согласуются между собой (рис. 26-27). В результате на входе в ТВС образуется сектор с пониженной температурой теплоносителя, примыкающий к петле №2 (рис 35). После подключения ГЦН на 11 с процесса в активную зону начинает поступать более холодный теплоноситель из петли №3 и в активной зоне начинает образовываться «холодный» сектор со стороны подключаемой петли, который постепенно расширяется и становится максимальным на 15 с. При этом сектор с пониженной температурой, примыкающий к петле №2, постепенно уменьшается. К 20 с процесса на входе в ТВС практически происходит выравнивание поля температур.

Из полученных результатов расчета следует удовлетворительное согласие с экспериментальными данными основных параметров (давление первого и второго контура, температура теплоносителя в петлях РУ). Отклонение расчетной температуры теплоносителя в петле №3 связано с отсутствием учета инерционности термопар, которое значительно влияет на их показания.

Заключение

Программный комплекс ТРАП-КС был использованы для проведения расчетов стандартных проблем V1000-CT1 (задача 3) и V1000-CT2 (задача 1) и сопоставления с данными испытаний и результатами расчета по другим кодам. Дополнительно для решения стандартной проблемы V1000-CT1 (задача 3) был применен программный комплекс ДКМ.

Результаты проведенных расчетов показывают, что рассматриваемые комплексы могут моделировать рассматриваемые режимы с приемлемой по отношению к экспериментальным данным точностью.

Список литературы

1. М.А.Быков, С.И.Зайцев, Ю.В.Беляев, Г.В.Алехин (ФГУП ОКБ «ГИДРОПРЕСС», Подольск, Россия), А.П.Егоров, В.И.Гусев (НИТИ им. Александрова, Сосновый Бор, Россия). Совершенствование программно-расчётного комплекса ТРАП-97. Учёт пространственных эффектов в реакторе. ТЕПЛОЭНЕРГЕТИКА №1, 2006 г (стр 48-52).

2. М.А.Быков, С.И.Зайцев, Г.В.Алехин, Курбаев С.А. (ФГУП ОКБ «ГИДРОПРЕСС», Подольск, Россия), А.П.Егоров, В.И.Гусев (НИТИ им. Александрова, Сосновый Бор, Россия). Развитие программного кмплекса ТРАП-97. Анализ чувствительности. ТЕПЛОЭНЕРГЕТИКА №1, 2006 г (стр 53-61).

3. G.V.Alekhin, Yu.V.Belyaev, S.I.Zaitsev, M.A.Bykov, Yu.N.Nadinsky, O.V.Kudryavtsev. Experimental&Design organization “Gidropress”, Podolsk, Russian Federation. Development of local parameters in WWER core considering 3-D kinetics. Proceedings of the thirtennth Symposium of AER. Dresden. Germany. 22-26 September 2003.

4. Драгунов Ю.Г., Быков М.А., Василенко В.А., Мигров Ю.А.. ОКБ «Гидропресс», НИТИ имени А.П.Александрова. Опыт применения и развитие расчетного кода КОРСАР для обоснования безопасности АЭС с ВВЭР. Теплоэнергетика №1, 2006, стр. 43-47.

5. G.V.Artyomov, A.V.Elshin, A.S.Ivanov et al. “Devepopment of neutron-physics models of varies types of reactors on the basis if unified algoritms of applied code package SAPFIR”. Proceedings of the 10-th International Seminar on Reactor Physics, Moscow, 2-6 September, 1997, p.34.

6. T.Takeda, Y.Komano. Extension of Askew's Coarse Mesh Method to Few-Group Problems for Calculating Two-Dimensional Power Distribution in Fast Breeder Reactors. Jounal of Nuclear Science and Technology, 15(7), pp. 523-532,1978.

7. Г.В.Алехин, С.А.Курбаев, М.А.Быков, С.И.Зайцев, Ю.В.Беляев. ФГУП ОКБ «Гидропресс», Подольск, Россия. Ю.А.Мигров, В.Г.Коротаев, О.В.Кувшинова. ФГУП НИТИ им. Александрова, Сосновый Бор, Россия. Кросс-верификация расчетных комплексов ТРАП-КС, ДКМ и Корсар/ГП по результатам динамических испытаний на действующих энергоблоках с ВВЭР-1000. The 5-th International conference. FSUE OKB “Gidropress”, Podolsk, Russia, 29 May - 1 June, 2007.

8. B. Ivanov, K. Ivanov, P. Groudev, M. Pavlova, V Hadjiev. VVER-1000 coolant transient benchmark. Phase 1 (V1000CT-1). Volume 1: Final Specifications (Revision 4). Main coolant pump (MCP) switching on. NEA/OECD, NEA/NSC/DOC (2002)6.

9. S. Kliem, T.Hohne, U.Rohde. FZ Dresden-Rossendorf; Institute of Safety Research, Germany. Y. Kozmenkov. SSC RF - IPPE, Obninsk, Russian Federation. DYN3D/ATHLET and ANSYS CFX calculations of the OECD VVER-1000 coolant transient benchmark. “Safety assurance of NPP with WWER”. The 5-th International conference. FSUE OKB “Gidropress”, Podolsk, Russia, 29 May - 1 June, 2007.

10. B. Ivanov, K. Ivanov, E. Royer, S.Aniel, U.Bieder, N.Kolev, P.Groudev. Overview and status of the OECD/DOE/CEA VVER-1000 coolant transient (V1000CT) benchmark. 14-th Symposium AER, 2004.

11. N.Kolev, S.Aniel, S.Royer, U.Bieder, D.Popov, Ts.Topalov. VVER-1000 Coolant Transient Benchmark (V1000CT), Volume II: Specification of the VVER-1000 Vessel Mixing Problems, OECD/NEA, March 2004.

Приложение

Таблица 1 -Принятые условные обозначения

Таблица 2 - Принятые сокращения

Таблица 3 - Граничные условия для расчета по ПМ КАМЕРА-В2 режима закрытия БЗОК ПГ

Таблица 4. Основные параметры в исходном и конечном состоянии. Закрытие БЗОК ПГ1. Эксперимент/ТРАП-КС

Таблица 5. Стационарное состояние. Сопоставление основных расчетных и экспериментальных параметров. Подключение ГЦН к трем работавшим

Подключение ГЦН

Подключение ГЦН

Подключение ГЦН

Подключение ГЦН

Размещено на Allbest.ru