WWW.KONF.X-PDF.RU
БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА - Авторефераты, диссертации, конференции
 

Pages:   || 2 | 3 | 4 |

«СОЗДАНИЕ НЕЙТРОННО-ФИЗИЧЕСКОГО КОДА НА ОСНОВЕ DSn-СХЕМ И НЕСТРУКТУРИРОВАННОЙ СЕТКИ ИЗ ПРЯМЫХ ПРИЗМ ДЛЯ УЧЕТА ПРОСТРАНСТВЕННЫХ НЕОДНОРОДНОСТЕЙ В НЕТВЭЛЬНОЙ ЧАСТИ АКТИВНОЙ ЗОНЫ РЕАКТОРНЫХ ...»

-- [ Страница 1 ] --

АКЦИОНЕРНОЕ ОБЩЕСТВО ОПЫТНО-КОНСТРУКТОРСКОЕ БЮРО «ГИДРОПРЕСС»

На правах рукописи

Николаев Александр Александрович

СОЗДАНИЕ НЕЙТРОННО-ФИЗИЧЕСКОГО КОДА НА ОСНОВЕ DSn-СХЕМ

И НЕСТРУКТУРИРОВАННОЙ СЕТКИ ИЗ ПРЯМЫХ ПРИЗМ ДЛЯ УЧЕТА

ПРОСТРАНСТВЕННЫХ НЕОДНОРОДНОСТЕЙ В НЕТВЭЛЬНОЙ ЧАСТИ



АКТИВНОЙ ЗОНЫ РЕАКТОРНЫХ УСТАНОВОК СО СВИНЦОВОВИСМУТОВЫМ ТЕПЛОНОСИТЕЛЕМ

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

Научный руководитель кандидат технических наук, А.В. Дедуль Москва - 2015

ОГЛАВЛЕНИЕ

Введение Глава 1 Постановка задачи и обзор литературы 15

1.1 Разработка подхода к сеточной аппроксимации активных зон 15 реакторных установок со свинцово-висмутовым теплоносителем и выбор совместимого с геометрией способа аппроксимации оператора переноса

1.2 Решение проблемы пре- и постпроцессинга в расчетном 27 обосновании проектных характеристик активных зон, выполняемом с использованием трехмерных неструктурированных сеток из прямых призм

1.3 Заключение по главе 1 Глава 2 "PMSNSYS" – нейтронно-физический код на основе DSn-схем 36 и неструктурированной комбинированной сетки из прямых призм для расчетов характеристик реакторов на быстрых нейтронах со свинцововисмутовым теплоносителем

2.1 Общие сведения о программе PMSNSYS 36

2.2 Математическая модель расчета уравнения переноса нейтронов и 37 гамма-квантов

2.3 Точность выбранных алмазных разностных схем 46

2.4 Дополнительные возможности сеточной аппроксимации 48

2.5 Заключение по главе 2 Глава 3 "REBEL" - программа пре- и постпроцессинга расчетов 50 нейтронно-физических характеристик реакторов на быстрых нейтронах со свинцово-висмутовым теплоносителем

3.1 Общее описание программы REBEL 50

3.2 Элементы твердотельного проектирования в REBEL 5

3.3 Создание расчетных сеток в REBEL 59

3.4 Постпроцессинг нейтронно-физических расчетов 66

3.5 Заключение по главе 3 Глава 4 Результаты применения программ PMSNSYS и REBEL для 69 расчета проблемно-ориентированных задач

4.1 Расчет двумерной модели тепловыделяющей сборки реакторной 69 установки со свинцово-висмутовым теплоносителем с поглощающим стержнем в центре

4.2 Расчет трехмерная модели тепловыделяющей сборки быстрого 74 реактора со свинцово-висмутовым теплоносителем с поглощающим стержнем в центре

4.3 Расчет двумерной модели полномасштабной активной зоны 95 быстрого реактора со свинцово-висмутовым теплоносителем

4.4 Заключение по главе 4 Заключение Приложение А. Дополнительные сведения о программе PMSNSYS 110 Список литературы

ОСНОВНЫЕ ОБОЗНАЧЕНИЯ И СОКРАЩЕНИЯ

–  –  –

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

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

Таким образом, продемонстрирована «блеск и нищета» гексагональных сеток применительно к расчету активных зон РУ с СВТ. Возникает вопрос, как можно продолжать использовать преимущества гексагональных сеток для моделирования твэлов, но при этом минимизировать модельные погрешности при описании нетвэльной части активной зоны (и бокового отражателя)?

Решение этого вопроса позволило бы снизить уровень консерватизма при проектировании, повысить качество расчетного обоснования, и, соответственно, конкурентоспособность проектов РУ с СВТ, что является важной и актуальной научно-технической задачей. При этом исчезла бы необходимость: 1) в разработке методик оценки соответствующих погрешностей расчета, 2) в непосредственном выполнении таких оценок,





3) в разработке и внедрении (громоздких) методик компенсации погрешностей,

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

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

Достижение цели исследования требует решения следующих задач:

1) разработки подхода к сеточной аппроксимации активных зон РУ с СВТ, которая обеспечивала бы возможность описания твэлов прямыми призмами с правильным шестиугольным основанием, и, одновременно, обеспечивала бы, при необходимости, возможность осуществления детальной сеточной аппроксимации элементов конструкций в нетвэльной части активной зоны. При этом согласованный с такой комбинированной сеточной аппроксимацией метод решения уравнения переноса не должен иметь существенных методических ограничений в расчете гетерогенных сред1;

2) разработки нейтронно-физического кода, обеспечивающего возможность выполнения расчетов с использованием указанной сеточной аппроксимации и метода решения уравнения переноса нейтронов и гаммаквантов;

3) обеспечения замыкания (комплексности) процесса расчетного обоснования с использованием указанного НФ-кода за счет разработки пре- и постобработчика для подготовки исходных данных к нейтронно-физическому расчету активной зоны РУ с СВТ и обработки полученных результатов;

4) апробация разработанного подхода к сеточной аппроксимации применительно к расчету проблемно-ориентированных задач (типовых конструкций активных зон РУ с СВТ).

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

С точностью до погрешности сечений взаимодействия излучения с веществом

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

Для решения попутно возникающей проблемы обеспечения системного сопровождения вычислений на неструктурированных сетках из прямых призм необходимо решение следующих задач: 1) автоматизации процесса создания сеточных моделей для НФ-расчетов, 2) генерации файлов входных данных для НФ-кода, 3) обработки результатов НФ-расчета. Для решения этих задач в диссертационной работе выполнен анализ современной ситуации в указанной предметной области, и с учетом проектных потребностей (опыта промышленной эксплуатации физических кодов в конструкторской организации ОКБ «ГИДРОПРЕСС») создан соответствующий пре- и постпроцессор.

Научная новизна результатов исследования 1 Автором впервые предложено при аппроксимации активных зон РУ с СВТ применять комбинированные сетки (в общем случае неструктурированные), состоящие из прямых призм с правильным шестиугольным, произвольным четырехугольным и произвольным треугольным основанием, в сочетании с аппроксимацией оператора переноса на основе алмазных (DD-схем) и алмазоподобных разностных схем (DDL-схем GQ-метода) DSn-приближения МДО.

2 Автором создана программа PMSNSYS, в которой впервые реализовано решение стационарного уравнения переноса нейтронов и гамма-квантов с использованием указанных численных схем на указанных комбинированных сетках.

3 Автором создан оригинальный интерактивный трехмерный графический пре- и постпроцессор REBEL, обеспечивающий автоматизацию процедур процесса создания трехмерных сеточных моделей для нейтронно-физических расчетов РУ с СВТ, генерацию файлов исходных данных для выполнения таких расчетов по программе PMSNSYS и постобработку результатов ее работы.

4 Автором лично и в соавторстве впервые получены результаты применения созданных программных средств к расчету проблемноориентированных задач.

На защиту выносится 1 Подход к аппроксимации активных зон РУ с СВТ, в соответствии с которым предлагается применять в совокупном сочетании: комбинированные сетки, состоящие из прямых призм с правильным шестиугольным, произвольным четырехугольным и произвольным треугольным основанием, с аппроксимацией оператора переноса на основе сочетания алмазных (DD-схем) и алмазоподобных (DDL-схем GQ-метода) разностных схем DSn-приближения МДО.

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

3 Интерактивный трехмерный графический пре- и постпроцессор REBEL, обеспечивающий создание сеточных моделей активных зон РУ с СВТ, подготовку файлов исходных данных и постобработку результатов расчетов.

4 Результаты применения программ PMSNSYS и REBEL к расчету проблемно-ориентированных задач.

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

Практическая направленность 1 Созданы программы PMSNSYS и REBEL, при создании которых учтен опыт промышленной эксплуатации нейтронно-физических кодов в конструкторской организации ОКБ «ГИДРОПРЕСС».

2 В 2011 г. на программу PMSNSYS в соответствии с ЕСПД выпущен комплект программной документации [8-15], включая отчет о верификации применительно к РУ СВБР-100 [8]. Программа PMSNSYS введена в опытную эксплуатацию в ОКБ «ГИДРОПРЕСС» для использования в расчетном обосновании РУ с ТЖМТ по приказу № 101 от 05.07.2011 г.;

3 В период с 2011 по 2013 гг. (REBEL c 2009 г.) программы REBEL и PMSNSYS были применены для обоснования проектных характеристик и безопасности активной зоны РУ СВБР-100, что отражено в материалах технического проекта [16-30].

4 По заказу ИБРАЭ РАН с учетом полученного в диссертационном исследовании опыта сформулированы технические требования [31, 32] на разработку нейтронно-физического решателя, пре- и постпроцессора для частного проекта «Коды нового поколения».

Апробация работы Результаты работы докладывались и получили положительную оценку на российских и международных научных семинарах и конференциях (Семинар Нейтроника (ГНЦ РФ-ФЭИ, 2011-2014 гг.), МНТК «50 лет БФС» (ГНЦ РФФЭИ, 2012 г.), МНТК «Инновационные проекты и технологии ядерной энергетики» (НИКИЭТ, 2012 г.), конференция молодых специалистов ОКБ «ГИДРОПРЕСС» (2010-2013 гг.). Отдельные части представленной работы неоднократно отмечались на различных научных конкурсах на лучшую научную работу, в том числе в 2011 г. отмечены дипломом победителя конкурса научных работ молодых ученых семинара Нейтроника-2011. За достигнутые успехи в развитии науки постановлением Главы города Подольска №94-П от 31.01.2014 автору выдан диплом и присвоено звание «Лучший работник науки 2013 г.».

Автор в соответствии с распоряжением Губернатора Московской области № 429-РГ от 22.12.2014 г. является лауреатом «Ежегодной премии Губернатора Московской области в сфере науки и инноваций для молодых ученых и специалистов» за 2014 г. за представленные на конкурсе результаты диссертационного исследования.

Публикации Основное содержание диссертации отражено в двух публикациях в рецензируемом журнале «Тяжелое машиностроение», в 13 научных работах в виде докладов на научно-технических конференциях.

Личный вклад автора Автором лично сформулировано выносимое на защиту предложение применять комбинированные расчетные сетки в сочетании с аппроксимацией оператора переноса на основе алмазных (DD-схем) и алмазоподобных (DDL-схем GQ-метода) разностных схем DSn-приближения МДО для расчетов активных зон РУ с СВТ. Автором лично разработаны программы REBEL и PMSNSYS.

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

Структура диссертации Диссертация состоит из четырех глав, введения, заключения, одного приложения и списка литературы. Общий объем диссертации составляет 130 страниц, в том числе 43 рисунка и 8 таблиц. Список литературы содержит 126 наименований.

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

В главе 2 рассмотрены основные особенности созданной в процессе диссертационного исследования МДО-программы PMSNSYS. Представлены конечно-разностные уравнения DD- и DDL-схем, используемые в программе для аппроксимации оператора переноса, в зависимости от типа освещенности прямых призм с выпуклым многоугольным основанием. Рассмотрена природа методических погрешностей DDL-схем на сетках, содержащих существенно несимметричные ячейки, способы ее оценки и снижения ее величины.

Представлены полезные особенности в строении сеток, поддерживаемые PMSNSYS.

В главе 3 выполнен обзор функциональных возможностей разработанной автором пре- и постпроцессинговой программы REBEL. Рассмотрены возможности REBEL в части создания трехмерных твердотельных моделей и дискретных пространственных сеток. Рассмотрен демонстрационный пример применения программы REBEL к расчету фрагмента активной зоны РУ с СВТ по коду PMSNSYS, включая визуализацию результатов расчета.

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

Благодарности Автор выражает благодарность научному руководителю работы А.В. Дедулю за доброту и терпение в научном руководстве автора, А.В. Воронкову (ИПМ им. М.В. Келдыша РАН) за постановку исходной задачи и драгоценное время, уделенное автору во время сотрудничества, преподавательскому коллективу физического факультета СГУ им. Н.Г. Чернышевского (в особенности А.И. Жбанову, а также Б.Е. Железовскому, А.Е. Васильеву, В.Н. Рачкову, В.П. Вешневу, Н.Г. Недогреевой, Н.В. Романовой), Г.Н. Мантурову (ГНЦ РФ-ФЭИ) за помощь и ценные пояснения по ряду проблемных вопросов, коллегам по работе (в особенности В.В. Кальченко, П.Б. Афанасьеву и В.В. Усенкову). Автор благодарит уважаемых Г.Л. Пономаренко (ОКБ «ГИДРОПРЕСС»), Л.П. Басса и О.В. Николаеву (ИПМ им. М.В. Келдыша РАН) за прочтение диссертационной работы и ценные советы по улучшению ее качества.

Особо автор благодарит своих коллег А.Н. Скобелева и М.В. Потапова за неоценимый вклад в освоение, тестирование, активное использование и рекомендации по совершенствованию разработанных автором программ, способствовавший их успешному ускоренному внедрению в практику проектных расчетов в ОКБ «ГИДРОПРЕСС».

Автор благодарит руководство ОКБ «ГИДРОПРЕСС», и, в особенности, В.С. Степанова, А.С. Зубченко и А.А. Сошникова, за постоянное внимание к работе, помощь и создание условий для ее выполнения.

В заключение автор считает наиболее уместным процитировать своего коллегу Гленна Сьдена, разработчика нейтронно-физического кода PENTRAN [33,34]: «Finally, and not the least, I thank God for wisdom to develop this code».

ГЛАВА 1 ПОСТАНОВКА ЗАДАЧИ И ОБЗОР ЛИТЕРАТУРЫ

1.1 Разработка подхода к сеточной аппроксимации активных зон реакторных установок со свинцово-висмутовым теплоносителем и выбор совместимого с геометрией способа аппроксимации оператора переноса В расчетном обосновании БР с ЖМТ активно используются программы, решающие уравнение переноса на регулярных гексагональных сетках [35-39].

Не является исключением и расчетное обоснование активных зон РУ с СВТ [6], где традиционно1 используется потвэльное приближение (ввиду достаточно небольших размеров таких активных зон). Гексагональная сеточная аппроксимация устанавливает удобную и практичную взаимосвязь между результатами расчета и каждым отдельно взятым твэлом. При этом методической погрешностью, связанной с неучетом мелкомасштабной пространственной гетерогенности, благодаря быстрому спектру нейтронов в активной зоне, вполне можно пренебречь [7].

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

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

В случае расчета интегральных величин (эффективность ПС, усредненная по объему плотность групповых потоков и т.п.) погрешностями подобных упрощений пространственной аппроксимации во многих случаях можно пренебречь. В особенности, применение таких упрощений может быть целесообразным на начальной стадии разработки проекта при выборе основных Определяющую роль в развитии практических подходов и программ для расчетого обоснования нейтроннофизических характеристик активных зон РУ с СВТ сыграли сотрудники ФГУП «ГНЦ РФ-ФЭИ»

(Н.Н. Новикова, Е.А. Земсков, и др.) и ИПМ им. М.В. Келдыша РАН (коллектив авторов ППП РЕАКТОР во главе с А.В. Воронковым) проектных решений, требующих многовариантного поиска. Если целью расчета является не выбор, а уточнение параметров установки, то вопрос согласованности сеточной аппроксимации с расчетной областью обостряется.

Особую остроту проблема погрешности гексагональной сеточной аппроксимации приобретает при выполнении поверочных расчетов, при оценке возможности снижения инженерных запасов, или уточнения величины запасов для закладывания в проект, и т.п.. Например, когда для этого необходимо рассчитать детальный ход НФ-функционала (например, тепловыделения) по радиусу одиночного ПС (рис. В.2, б), или по радиусу элементов ПЭЛ кластерной конструкции (рис. 1.1). Ясно, что применение гексагональной аппроксимации с размером ячеек, равному шагу твэлов в активной зоне (1,36 см, см. Введение), для решения этой задачи малоэффективно.

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

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

Загрузка...

Использование диффузионной теории для решения проблемы учета эффектов пространственной гетерогенности зачастую может носить формальный характер из-за возникающих случаев грубого нарушения условий ее применимости (наличие сильнопоглощающей среды или источников нейтронов, выраженная пространственная гетерогенность, необходимость расчета транспорта гамма-квантов для расчета их вклада в пространственное распределение тепловыделения). Практика расчетов показывает, что из наиболее распространенных на сегодняшний день методов решения уравнения переноса (ММК, МДО, МСГ, МПГ) лучше всего приспособлены для решения подобных задач (расчеты активных зон без пространственной гомогенизации) именно ММК и МДО [40-42].

Платой за предоставляемые в этом случае методические преимущества является вычислительное удорожание расчетов, из-за чего использование только вероятностных МК-кодов (т.е. полное замещение детерминистических кодов) в массовых проектных расчетах на сегодняшний день представляется затруднительным даже с учетом возможностей современных ЭВМ [43-47].

Хотя надо признать, что в этом направлении все же достигнуты определенные успехи [48-50]. В то же время, производительность ЭВМ создает предпосылки для широкого внедрения сеточных МДО-кодов в расчетное обоснование НФХ активных зон РУ с СВТ, во многом, конечно, благодаря малым размерам таких активных зон. Это важное обстоятельство – как с точки зрения небольшого количества ячеек в расчетной модели, так и с точки зрения скорости сходимости решения – жесткосвязные задачи сходятся быстрее. Все это во многом обусловливает остановку выбора на этом методе при разработке перспективных кодов для активных зон РУ с СВТ.

Этот вывод подкрепляется практическим опытом автора, как одного из основных исполнителей в процессе расчетного обоснования активных зон РУ с СВТ в ОКБ «ГИДРОПРЕСС». Фактически начало активного внедрения кинетического приближения в расчеты именно активных зон РУ с СВТ в ОКБ «ГИДРОПРЕСС» датируется концом 2009 г.. Тогда для решения проектных задач применялся кинетический DSn-модуль KIN3D6 (автор Е.П. Сычугова) комплекса программ РЕАКТОР-ГП [51]. Проведение отдельных проектных расчетов с использованием KIN3D6 в потвэльной постановке в те годы было возможно даже на одноядерных ПЭВМ благодаря реализованным в программе эффективным алгоритмам ускорения сходимости внутренних и внешних итераций [52]. Распараллеливание расчетных модулей комплекса, предназначенных как для расчета радиационной защиты [53], так и для расчета активной зоны, позволило в продолжение 2010 г. при использовании многоядерной ПЭВМ (16 ядер) и OpenMP-версии KIN3D6 выполнить цикл расчетов [54,55] в обоснование кампании активной зоны РУ СВБР-100, в обоснование условий работы твэлов, характеристик ПС и условий их работы.

К сожалению, комплекс программ РЕАКТОР-ГП не включал в свой состав МДО-программы, решающие уравнения переноса на сетках, в требуемой мере согласованных с расчетной областью. Этот факт отчасти послужил толчком к выполнению настоящего диссертационного исследования.

В настоящее время разработаны множество трехмерных МДО-программ (например MCCG3D [56], САТУРН-3 [57], CRONOS2-SN [41], APPOLLO2 [58], РАДУГА-Т [59], ATTILA [60], SNAC [61], MOCK-3D [62], NEWTRNX [63]), возможности которых обеспечивают решение проблемы учета пространственной неоднородности в активных зонах «в лоб» – путем явного описания расчетной области (рис. 1.2).

Рис. 1.2 Детальное сеточное описание фрагмента нетвэльной части активной зоны Последние пять из приведенных программ выполняют решение в полном смысле в произвольной трехмерной геометрии, а в первых четырех рассматривается не полностью нерегулярная трехмерная геометрия – здесь регулярная геометрическая структура сетки повторяется вдоль оси Z («SmaRT-Z» сетка, [61]).

Такой упрощенный подход представляется наиболее целесообразным для решения проектных задач, т.к. явный учет этой регулярной структуры позволяет существенно снизить трудоемкость расчетов [64]. Это важное преимущество для проектного кода, рассчитанного на интенсивную эксплуатацию. Полностью трехмерные программы с численными схемами повышенного порядка точности более востребованы для решения задач радиационной защиты.

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

Лишь ограниченное количество МДО-кодов предоставляют возможность расчета на подобных сетках. Например, это отечественная программа MCCG3D.

Программа MCCG3D реализует метод длинных характеристик, и хорошо известна в научном сообществе. Реализованные в программе возможности позволяют получать с ее использованием референсные решения (с точностью до используемых групповых констант) для очень широкого класса реакторных и защитных задач [41,42,64] на сетках, составленных из прямых призм с произвольным многоугольным основанием. Более того, геометрический блок программы не ограничивается полигональной (в плоскости) аппроксимацией – он поддерживает криволинейную форму поверхности моделируемых объектов.

Специально для расчетов активных зон в MCCG3D, кроме SC-схемы, также предусмотрена возможность применения алмазной разностной схемы вдоль характеристических направлений из-за определенных преимуществ в точности.

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

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

Вычислительной привлекательностью [65] для расчетов активных зон обладает широкоиспользуемая во многих отраслевых отечественных и зарубежных программах линейная алмазная разностная схема (DD – diamond difference). В соответствии с этой схемой поток в данной точке определяется в виде среднего арифметического между значениями потока в соседних счетных точках [66]. Эффективный размер расчетных ячеек DD-схемы составляет вплоть до двух длин свободного пробега [65] (для БР с ЖМТ примерно 3-8 см, в зависимости от свойств конкретной среды). DD-схема имеет второй порядок точности расчета интегральных величин, и обеспечивает достаточно высокую аппроксимацию решения при умеренных расчетных затратах [67].

Наиболее известный пример подобной МДО-программы, решающей уравнение переноса на неструктурированных сетках (но к сожалению, двумерных), состоящих из произвольных выпуклых многоугольников, с привлечением разностных схем – программа ЭКРАН DSn-метода (ВНИИЭФ) [68].

Развитием этой программы является трехмерный кинетический модуль комплекса САТУРН-3 (предшественник кинетического модуля комплекса SERENA [69]). Его особенностью является применение для конечно-разностной аппроксимации уравнения переноса разностных линейных схем, отличительной чертой которых является то, что наряду с традиционно используемыми в DD-схемах средними значениями искомой функции на гранях многогранника и в самой ячейке привлекаются величины на ее боковых ребрах. Такой расширенный шаблон позволяет обеспечить точность, сравнимую с точностью традиционных разностных схем на сетках из ячеек правильной формы. К сожалению, в отличие от своего двумерного предшественника, поддерживаемая в САТУРН-3 геометрия может быть составлена из прямых призм, имеющих только произвольное четырехугольное основание [56].

Анализ литературы показывает, что в 1982 г. Т.Р. Хиллом и Р.Р. Патерностером [70] были предложены двумерные разностные схемы для сеток из произвольных выпуклых четырехугольников и треугольников (DDL-схемы GQ-метода, в терминологии [71]), родственные реализованным в программе ЭКРАН. Схемы GQ-метода существенно проще последних, как в вычислительном плане, так и в плане реализации. Принципиальное отличие состоит в том, что в GQ-методе для приближения решения не привлекаются значения в вершинах четырехугольника и не решаются соответствующие дополнительные уравнения. Но вместе с тем (и поэтому) схемы GQ-метода несколько проигрывают схемам ЭКРАН (или САТУРН-3, что, в общем, то же) в методической точности, величина которой определяется степенью отличия ячеек от правильной формы [72]. Известно, что понижение порядка точности метода в зависимости от степени искажения ячеек вообще присуща численным схемам, и этот вопрос решается измельчением ячеек в проблемной области. Не являются исключением и схемы GQ-метода.

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

Схемы Т.Р. Хилла и Р.Р. Патерностера реализованы ими в программе LaMEDOC. Р.Е. Алькуфом для этих схем разработан и опробован DSA-метод ускорения сходимости [73], после чего схемы GQ-метода были реализованы в программе входящей в состав МДО-комплекса TWODANT-GQ [74], DANTSYS [72]. В [71] подчеркнуто, что GQ-метод простой и быстрый («The GQ method is simple, fast, and is compatible with popular diffusion synthetic acceleration techniques»), в [70] – что вычислительно-эффективный и быстрый «computational-efficient and fast».

Следует также отметить, что двумерные разностные DDL-схемы для произвольных выпуклых четырехугольников, аналогичные [70], предложены в 2000 г. в работе [78]. Основное различие между [70] и [78] состоит в записи DDL-соотношений для случая, когда освещены три стороны четырехугольника.

Как показывают результаты численных расчетов, выполненных автором, DDL-соотношения, предложенные Т.Р. Хиллом и Р.Р. Патерностером для такого варианта освещенности ячейки, несколько точнее предложенных в [78].

Этот результат объясняется тем, что в [70] DDL-соотношение для выходящего излучения выводится из соотношения общего баланса, а в [78] для его определения используется просто DD-соотношение.

Далее возникает вопрос, как двумерные схемы GQ-метода можно использовать в целях работы применительно к трехмерным прямым призмам с произвольным четырехугольным и треугольным основанием? Ответ на этот вопрос следующий. В соответствии с [67] алмазные разностные схемы допускают простое обобщение на требуемый трехмерный случай прямых призм с произвольным многоугольным основанием. Для этого достаточно связать входящее и выходящее через торцевые грани прямой призмы излучение с помощью алмазного разностного соотношения, как это, например, сделано в [75] для прямых призм с треугольным основанием. Таким образом, рассмотренные алмазные разностные схемы (DD- и DDL-схемы) вместе с соотношением баланса [76] (нулевого порядка) в ячейке I для направления полета m с освещенными (i) и неосвещенными (j) гранями (сторонами) Ai Aj VI VI QI i j tI I i j образуют замкнутую система уравнений, которая легко разрешается относительно выходящего излучения на неосвещенных гранях призмы.

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

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

Указанный подход, в соответствии с которым предлагается для сеточной аппроксимации активных зон РУ с СВТ использовать сетки, состоящие одновременно из прямых призм с правильным шестиугольным и произвольным четырехугольным (и треугольным) основанием, а для аппроксимации оператора переноса использовать алмазные разностные схемы (DD-схемы, DDL-схемы GQ-метода), является новым. Он впервые предложен автором в 2011 г. в [77].

Рис. 1.3 Пример применения комбинированной сетки для моделирования ТВС (гомогенно шестиугольниками представлены твэлы и гетерогенно представлены поглощающие элементы в центре ТВС) Оригинальность предлагаемого в диссертационной работе решения, в отличие от рассмотренных ЭКРАН, САТУРН-3 и других MCCG3D, аналогичных, заключается:

во-первых, в отказе от поддержки НФ-кодом произвольного полигонального основания прямых призм (в отличие от MCCG3D) и реализация комбинированных сеток из ячеек только требуемых типов (на базе только правильных шестиугольников и произвольных четырехугольников1), и именно для трехмерного, а не двумерного случая комбинированных сеток (как это реализовано в ЭКРАН);

- во-вторых, в отличие от ЭКРАН, САТУРН-3, в применении к именно таким комбинированным сеткам вычислительно более дешевых разновидностей численных схем (двумерных DDL-схем GQ-метода для произвольных четырехугольников и для правильного DD-схем шестиугольника). При этом для получения решения в трехмерном случае указанные схемы применяются для получения решения на боковых гранях, а для связи входящего и выходящего излучения на торцевых гранях призмы используется только DD-схема.

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

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

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

–  –  –

в а – гексагональная аппроксимация, б – аппроксимация прямыми призмами с произвольным четырехугольным основанием, в – аппроксимация с использованием предлагаемых в работе комбинированных сеток Рис. 1.4 Преимущество применения комбинированных сеток для сеточной аппроксимации типовых конструкций активной зоны РУ с СВТ Решение проблемы пре- и постпроцессинга в расчетном 1.

обосновании проектных характеристик активных зон, выполняемом с использованием трехмерных неструктурированных сеток из прямых призм Практическй комфорт в применении проектного кода, решающего задачи на столь нетривиальных (в сравнении с классическими регулярными структурированными сетками) сеточных моделях, требует решения задачи согласованной по геометрии подготовки файлов исходных данных и постобработки результатов НФ-расчетов. Эта задача не может быть решена вручную с разумными трудозатратами [45]. Более того, столь важные для ядерной и радиационной безопасности результаты нейтронно-физических расчетов активной зоны должны быть наглядны и легко проверяемы (в том числе в части подготовки расчетных моделей [79], тем более, когда речь идет о нетривиальных сетках, разработка и внедрение которых выполняется впервые).

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

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

Фактически для создания сеточной модели необходимо (рис. 1.5):

1) создание набора регулярно расположенных гексагональных ячеек для аппроксимации твэлов; 2) создание типовых наборов ячеек (из произвольных четырехугольников и треугольников) для моделирования каналов СУЗ;

3) создание ячеек для аппроксимации бокового отражателя.

Рис. 1.5 Рассмотренный вариант сеточной аппроксимации в области границ активной зоны и бокового отражателя (с детализацией канала СУЗ) Ясно, что создание подобных сеток не требует применения сложных сеточных генераторов. Этот вопрос на хорошем практическом уровне вполне может быть решен в рамках элементарного препроцессора, единственное требование к которому – визуализация создаваемой сетки. Создание сетки для НФ-расчетов БР подразумевает также наполнение ее информацией о пространственном распределении материалов и их ядерно-физическом составе.

Этот процесс очень трудоемкий в случае, если его выполнять вручную.

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

Еще большие трудности возникают в случае, когда нужно получить гомогенизированные ядерно-физические составы конструкций, не согласованных с сеткой. Здесь также необходима автоматизация. Решение этой проблемы можно осуществлять разными способами, например как это предложено в [45] с использованием лучевого трассирования. Однако указанный метод требует пропускания достаточного количества лучей для обеспечения выхода решения на насыщение. Возможен более примитивный вариант решения проблемы (и, в ряде случаев более времязатратный, что, впрочем, зависит от типов трехмерных тел – громоздкие многогранные тела обрабатываются дольше, чем тела, поверхности которых заданы аналитически, например цилиндры или трубы), при котором объемные доли конструкций, входящих в ячейки сетки, вычисляются «в лоб» - аналитически (рис. 1.6).

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

Трехмерные фрагменты пересечения трубы с ячейками раскрашены в разные цвета Следующий этап – генерация файлов исходных данных для программы нейтронно-физического расчета из подготовленной сеточной модели. По сути это вопрос чисто технический. Для его разрешения достаточно знать особенности синтаксиса файла начальных данных конкретной МДО-программы. После этого легко решается задача записи формализованного (во внутреннем представлении препроцессора) представления сетки в текстовый (или бинарный) файл входных данных к МДО-программе. Такая подпрограмма конвертации, например, для модулей DIFF3D и KIN3D комплекса РЕАКТОР-ГП, может содержать менее 200 строк исходного кода [80] на языке фортран.

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

Естественным образом возникает вопрос, как обстоят дела в данной предметной области в России и мире? Какие в этих вопросах реализуются подходы, традиции и тренды?

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

- комплекс программ РЕАКТОР-ГП [81,82]: подготовка моделей расчета активной зоны – программа G-Core, подготовка моделей расчета защиты методами комбинаторной геометрии – модуль GNCOMP (G-Comp), постобработка – программа ARCSERVICE («Редактор Архива»);

- программная оболочка ShIPR [83] (программы DP3 и DP3T и редактор расчетных сеток VIZAGIST);

- программа MCCG3D [84] и «MCCG3D Geomery Editor»;

- программный комплекс ModExSys для математического обеспечения программ TRIGEX, MMKKENO [85,86].

Интересный современный программный продукт, разработанный в АО «НИКИЭТ» (г. Москва), представлен в [87]. Он представляет собой графический псевдотрехмерный препроцессор для разработки моделей экспериментальных сборок на стенде БФС и обеспечивает подготовку исходных данных для ряда нейтронно-физических программ (преимущественно для программ ММК). Данный препроцессор успешно используется для верификации НФ-кодов проекта БРЕСТ-ОД-300. Функции постпроцессинга в описанной версии программы не реализованы. Однако ее дальнейшее развитие предполагает реализацию таких функций, и, кроме того, предполагает расширение списка кодов, для которых будет осуществляться препроцессинг.

Другая разработка АО «НИКИЭТ» – интегральная система кодов «Платформа-БРЕСТ» Кроме прочего, «Платформа-БРЕСТ» с [88].

использованием графического интерфейса пользователя обеспечивает разработку расчетных моделей (в псевдотрехмерной постановке и только в гексагональной геометрии) проектируемых в АО «НИКИЭТ» БР, и генерацию файлов исходных данных для ряда используемых в организации программ нейтронно-физического расчета.

За рубежом ситуация в настоящее время в целом выглядит несколько иначе. Анализ литературы [41, 89-97] показывает, что создание модели и постобработка результатов здесь преимущественно ориентирована на продвинутые коммерческие пре-, постпроцессоры (BOT3P 5.1, SolidWorks, ANSYS ICEM-CFD, Catia, и других), а также специальные программывизуализаторы (Tecplot, VisIt, ParaView, Surfer). При этом связь (передача данных) между ними осуществляется через специальные интерфейсные файлы (стандартов IGES, STEP, BREP, ACIS и т.п.), чтение/запись которых осуществляется, в том числе, посредством специально разрабатываемых утилит. Отметим также, что комплекс BOT3P активно используется и в России (для МДО-кодов КАТРИН и TORT). Важное отличие отечественных пре- и постпроцессинговых программ от используемых за рубежом графических средств пре- и постобработки то, что последние, за редким исключением, полностью трехмерные (а не псевдотрехмерные) – таков тренд.

Таким образом, можно констатировать, что в соответствии с тенденциями современности НФ-код должен быть обеспечен пре- и постпроцессором с графическим интерфейсом пользователя. Без этого эффективная разработка, тестирование и внедрение в проектную деятельность МДО-программ на неструктурированных сетках на сегодняшний день практически немыслима.

Тем более что развитие ЭВМ и графических библиотек предоставляют для этого все необходимые возможности.

В соответствии с традицией, сложившейся в отечественной отрасли, такой пре- и постпроцессор, как правило, разрабатывается автором НФ-кода (или при его плотном участии, либо коллективом активных пользователей НФ-кодов).

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

Приспособление уже имеющихся (отраслевых или зарубежных) пре- и постпроцессоров (в том числе открытых, типа зарубежных SALOME, легитимного сопряжено с трудностями их получения и NURESIM) использования, необходимостью ориентации разрабатываемой Что также касается и открытых программных продуктов, т.к. при этом возникают понятные трудности получения надежных гарантий по обеспечению технической поддержки и прав на использование на долгосрочную перспективу МДО-программы на их конкретные возможности, зачастую фиксированные, что является отрицательной стороной привлечения готовых программных продуктов. В любом случае, применение «чужих» программных средств требует выполнения действий по взаимоувязке с НФ-программой.

Указанная общая ситуация определяет общий фон требований к решению проблемы пре- и постпроцессинга при расчете активных зон РУ с СВТ, которые в рамках настоящего диссертационного исследования приняты следующими:

создание «собственного» пре- и постпроцессора в единой пользовательской среде (по аналогии с BOT3P 5.1, SolidWorks, ANSYS ICEM-CFD) – без разделения на отдельные программы для препроцессинга и для постпроцессинга (как это, например, реализовано в комплексе РЕАКТОР-ГП);

- реализация трехмерного интерактивного графического интерфейса пользователя, что достаточно ново для «собственных» пре- и постпроцессоров детерминистических кодов отечественной отрасли (исключение составляют, например, сервисные программы НФ-кодов ВНИИЭФ в комплексе SERENA);

- обеспечение возможности формирования универсальных расчетных моделей – отсутствие привязки к типу реактора, экспериментальных критических сборок, тестовых задач. При этом необходимо обеспечение возможности заполнения сеток из заранее подготовленных твердотельных моделей (гомогенизация материалов в ячейках сетки пропорционально объемным долям пересекающихся с ней конструкций – достаточно уникальная возможность в сравнении с другими пре- и постпроцессорами). Все это обеспечивает большую гибкость и определенные преимущества перед многими отечественными пре- и постпроцессорами (РЕАКТОР-ГП, ModExSys, Платформа-БРЕСТ, ПC [87] и др.), ориентированными либо на отдельно взятые проекты БР, либо на решение преимущественно отдельно взятого круга задач (проектное обоснование, либо верификация кодов).

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

Получение опыта и освоение технологии таких разработок открывает возможности по созданию и эксплуатации «экспериментальной площадки» для отработки перспективных научно-технических решений в области расчетного обоснования проектов (разработка оригинальных способов визуализации и представления данных, поддержка разработки или внедрения численных схем для неструктурированных многогранных сеток [98], разработка взаимосвязей «чертеж» - «расчет», разработка и реализация интерфейсной взаимосвязи между дисциплинами в области мультифизичных расчетов, и т.д.).

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



Pages:   || 2 | 3 | 4 |
Похожие работы:

«СМИРНОВ АНДРЕЙ НИКОЛАЕВИЧ ВЕРИФИКАЦИЯ ЦИФРОВЫХ ДИНАМИЧЕСКИХ МОДЕЛЕЙ КРУПНЫХ ЭНЕРГООБЪЕДИНЕНИЙ ПО ДАННЫМ СМПР Специальность 05.14.02 – Электрические станции и электроэнергетические системы Диссертация на соискание ученой степени кандидата технических наук Научный руководитель – кандидат технических наук, доцент А.Х. Есипович Санкт-Петербург – 2013 ОГЛАВЛЕНИЕ ВВЕДЕНИЕ ГЛАВА 1...»

«Егоров Денис Эдуардович СОВЕРШЕНСТВОВАНИЕ МЕТОДОВ РАСЧЕТА МНОГОФУНКЦИОНАЛЬНЫХ ФИЛЬТРОКОМПЕНСИРУЮЩИХ УСТРОЙСТВ ДЛЯ СЕТЕЙ 10 0,4 КВ Специальность: 05.14.02 «Электрические станции и электроэнергетические системы» Диссертация на соискание ученой степени кандидата технических наук Научный руководитель – Доктор технических наук, профессор В. П. Довгун Красноярск 2015 ОГЛАВЛЕНИЕ ОГЛАВЛЕНИЕ ВВЕДЕНИЕ ГЛАВА 1 Проблемы обеспечения качества...»

«УДК 621.039.5 Жемков Игорь Юрьевич НАУЧНО-МЕТОДИЧЕСКОЕ СОПРОВОЖДЕНИЕ ЭКСПЛУАТАЦИИ ИССЛЕДОВАТЕЛЬСКОГО РЕАКТОРА НА БЫСТРЫХ НЕЙТРОНАХ Специальность: 05.14.03. Ядерные энергетические установки, включая проектирование, эксплуатацию и вывод из эксплуатации Диссертация на соискание ученой степени доктора технических наук Научный консультант: доктор технических наук,...»

«САУШИН Илья Ирекович ТУРБУЛЕНТНОСТЬ В ПОГРАНИЧНОМ СЛОЕ ПУЛЬСИРУЮЩЕГО ПОТОКА Специальность 01.02.05 Механика жидкости, газа и плазмы Диссертация на соискание...»

«НИПРУК ОКСАНА ВАЛЕНТИНОВНА ХИМИЯ СЛОЖНЫХ ОКСИДНО-СОЛЕВЫХ СОЕДИНЕНИЙ УРАНА (VI). СИНТЕЗ, СТРОЕНИЕ, СОСТОЯНИЕ В ВОДНЫХ РАСТВОРАХ 02.00.01 – Неорганическая химия химические науки ДИССЕРТАЦИЯ на соискание ученой степени доктора химических наук Научный консультант: доктор химических наук, профессор Н.Г. Черноруков...»

«ПЕТРОВ ИЛИЯН ИВАНОВ Эволюция структур мировых и европейских энергетических рынков и перспективы развития газотранспортных сетей в Юго-Восточной Европе с участием Болгарии и России Специальность 08.00.14 „Мировая экономика Диссертация на...»

«Паянен Рейно Игоревич ИССЛЕДОВАНИЕ И РАЗРАБОТКА РЕГЕНЕРАТИВНЫХ ЦИКЛОВ НА ЭЛЕГАЗЕ Специальность 01.04.14 – Теплофизика и теоретическая теплотехника ДИССЕРТАЦИЯ на соискание учёной степени кандидата технических наук Научный руководитель: Доктор технических наук Мазурин И.М. МОСКВА – 2015 ОГЛАВЛЕНИЕ: ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ Глава 1. НАКОПЛЕНИЕ ПРИМЕСЕЙ В ЭЛЕГАЗЕ В ВЫСОКОВОЛЬТНОМ...»

«Гавриленко Сергей Сергеевич ИССЛЕДОВАНИЕ И РАЗРАБОТКА ПОДХОДОВ К ПРОЕКТИРОВАНИЮ ВОДОПОДГОТОВИТЕЛЬНЫХ УСТАНОВОК ДЛЯ ПАРОГАЗОВЫХ ТЭС Специальность 05.14.14 – Тепловые электрические станции, их энергетические системы и агрегаты Диссертация на соискание ученой степени кандидата технических наук Научный руководитель: доктор...»

«ПЕТРОВ ИЛИЯН ИВАНОВ Эволюция структур мировых и европейских энергетических рынков и перспективы развития газотранспортных сетей в Юго-Восточной Европе с участием Болгарии и России Специальность 08.00.14 „Мировая экономика Диссертация на...»

«ЕВДОКИМОВА НАТАЛЬЯ ГЕОРГИЕВНА РАЗРАБОТКА НАУЧНО-ТЕХНОЛОГИЧЕСКИХ ОСНОВ ПРОИЗВОДСТВА СОВРЕМЕННЫХ БИТУМНЫХ МАТЕРИАЛОВ...»

«ЧУВАРАЯН Александра Асватуровна ЭНЕРГЕТИЧЕСКИЙ ФАКТОР В ПОЛИТИКЕ РОССИИ НА БЛИЖНЕМ И СРЕДНЕМ ВОСТОКЕ Специальность 23.00.04 политические проблемы международных отношений, глобального и регионального развития Диссертация на соискание ученой степени кандидата политических наук Научный руководитель Почётный работник науки и техники РФ, Доктор военных наук, профессор Анненков В.И. Научный...»

«УРАЗОВ ДМИТРИЙ ЮРЬЕВИЧ НАУЧНОЕ ОБЕСПЕЧЕНИЕ ПРОЦЕССА ТЕРМОВЛАЖНОСТНОЙ ОБРАБОТКИ КОЛБАСНЫХ ИЗДЕЛИЙ СПЕЦИАЛЬНОСТЬ 05.18.12 ПРОЦЕССЫ И АППАРАТЫ ПИЩЕВЫХ ПРОИЗВОДСТВ ДИССЕРТАЦИЯ на соискание ученой степени кандидата технических наук Научный руководитель: доктор технических наук,...»

«Суворова Ирина Александровна ВЫБОР СЕЧЕНИЙ ПРОВОДНИКОВ И РАЦИОНАЛЬНЫХ НАПРЯЖЕНИЙ РАСПРЕДЕЛИТЕЛЬНЫХ ЭЛЕКТРИЧЕСКИХ СЕТЕЙ В СОВРЕМЕННЫХ УСЛОВИЯХ Специальность 05.14.02 Электростанции и электроэнергетические системы ДИССЕРТАЦИЯ на соискание ученой степени кандидата технических наук Научный руководитель доктор технических наук, профессор В.В.Черепанов Киров, 2015 2 Содержание СОДЕРЖАНИЕ Введение Глава 1 Анализ состояния распределительных...»

«Варков Артем Александрович РАЗРАБОТКА И ИССЛЕДОВАНИЕ СИСТЕМЫ УПРАВЛЕНИЯ МАНИПУЛЯЦИОННЫМ ПРОМЫШЛЕННЫМ РОБОТОМ НА БАЗЕ КОНТРОЛЛЕРА ДВИЖЕНИЯ Специальность 05.09.03 – «Электротехнические комплексы и системы» ДИССЕРТАЦИЯ на соискание ученой степени кандидата технических наук Научный руководитель: доктор...»

«ТРУФАНОВ Виктор Васильевич МЕТОДИЧЕСКОЕ ОБЕСПЕЧЕНИЕ ПЕРСПЕКТИВНОГО РАЗВИТИЯ ЭЛЕКТРОЭНЕРГЕТИЧЕСКИХ СИСТЕМ РОССИИ В СОВРЕМЕННЫХ УСЛОВИЯХ Специальность 05.14.02 Электрические станции и электроэнергетические системы Диссертация на соискание ученой степени доктора технических наук Научный консультант: Воропай Николай Иванович,...»

«Хуршудян Смбат Размикович Оптимизация режимов ПГУ при участии ее в регулировании мощности и частоты в энергосистеме (на примере ПГУ-450) Специальность 05.13.06 Автоматизация и управление технологическими процессами и...»

«СИМОНЯН Левон Ашотович Экономические проблемы развития электроэнергетического хозяйства стран Совета сотрудничества арабских государств Персидского залива (1970-2013 гг.) (08.00.14 – мировая экономика) Диссертация на соискание ученой степени кандидата экономических наук Научный руководитель – кандидат экономических наук Смирнова Галина Ивановна Москва – ОГЛАВЛЕНИЕ Введение Глава I....»

«Дубоносов Антон Юрьевич ГИДРОДИНАМИКА ВХОДНЫХ ЦИЛИНДРИЧЕСКИХ КОЛЛЕКТОРОВ ТЕПЛООБМЕННЫХ АППАРАТОВ ТЕПЛОЭНЕРГЕТИЧЕСКИХ УСТАНОВОК Специальность: 05.14.14 «Тепловые электрические станции, их энергетические системы и агрегаты » Диссертация на соискание ученой степени кандидата технических наук Научный руководитель: д-р технических наук, профессор А.М. Гапоненко г. Краснодар 2015 ОГЛАВЛЕНИЕ ВВЕДЕНИЕ ГЛАВА 1 ОБЗОР...»

«Марьяндышев Павел Андреевич СОВЕРШЕНСТВОВАНИЕ ТЕХНОЛОГИИ ЭНЕРГЕТИЧЕСКОГО ИСПОЛЬЗОВАНИЯ ДРЕВЕСНОГО БИОТОПЛИВА Специальность 05.14.04 «Промышленная теплоэнергетика» Диссертация на соискание ученой степени кандидата технических наук Научный руководитель: д.т.н, профессор...»

«Николаев Александр Александрович УЧЕТ ЭФФЕКТОВ ПРОСТРАНСТВЕННОЙ ГЕТЕРОГЕННОСТИ В МНОГОМЕРНЫХ РАСЧЕТАХ РЕАКТОРОВ С ЖИДКОМЕТАЛЛИЧЕСКИМ ТЕПЛОНОСИТЕЛЕМ СВИНЕЦ-ВИСМУТ Специальность 05.14.03 Ядерные энергетические установки, включая проектирование, эксплуатацию и вывод из эксплуатации Диссертация на соискание учной степени кандидата технических наук Научный руководитель кандидат технических наук, А.В. Дедуль...»









 
2016 www.konf.x-pdf.ru - «Бесплатная электронная библиотека - Авторефераты, диссертации, конференции»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.