WWW.DISSERS.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА

   Добро пожаловать!


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

ТИХОНЕНКО АЛЕКСЕЙ ВИТАЛЬЕВИЧ

КОМПЛЕКСЫ ПРОБЛЕМНО-ОРИЕНТИРОВАННЫХ ПРОГРАММ В СИСТЕМАХ СИМВОЛЬНОЙ МАТЕМАТИКИ

Специальность: 05.13.18 – Математическое моделирование, численные методы и комплексы программ

Автореферат диссертации на соискание ученой степени доктора физико-математических наук

Обнинск – 2009

Работа выполнена в Обнинском государственном техническом университете атомной энергетики (г. Обнинск).

Официальные оппоненты доктор физико-математических наук, профессор Зингерман Константин Моисеевич доктор физико-математических наук, профессор Кудряшов Николай Алексеевич доктор физико-математических наук Якобовский Михаил Владимирович Ведущая организация Объединенный институт ядерных исследований, Лаборатория информационных технологий, г. Дубна, Московская обл.

Защита состоится 27 ноября 2009 г. в 1400 час. на заседании диссертационного совета Д212.263.04 в Тверском государственном университете по адресу:

170000, г. Тверь, ул. Желябова, 33, ауд. 52.

С диссертацией можно ознакомиться в библиотеке ТвГУ по адресу: 170000, г. Тверь, ул. Володарского, 44а.

Автореферат разослан « » _______________ 2009 г.

Ученый секретарь диссертационного совета Д212.263.доктор технических наук, профессор В.Н. Михно

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

АКТУАЛЬНОСТЬ ТЕМЫ Современную фундаментальную и прикладную науку уже невозможно представить без математического моделирования, а само математическое моделирование – без повсеместного использования вычислительной техники, компьютеров и математического обеспечения разного типа. Являясь научным методом познания и проектирования, в основе которого лежит работа не с самим объектом или процессом, а с его моделью, математическое моделирование сочетает в себе многие достоинства теории и эксперимента: вычислительные эксперименты с моделями дают возможность достаточно полного изучения объектов, что недоступно при чисто теоретических подходах. При этом возможности моделирования и сложность решаемых задач зависят как от применяемых вычислительных методов и мощности вычислительной техники, так и уровня математического обеспечения.

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

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

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

ССМ – интерактивные многофункциональные компьютерные системы высокого ин теллектуального уровня, сочетающие достаточную простоту использования с мощью самых современных математических инструментов: они построены на технологиях символьных и численных вычислений с произвольной точностью и продвинутых математических алгоритмах. К наиболее мощным ССМ на сегодняшний день относятся MAPLE, MATHEMATICA, MATLAB и MATHCAD. Интерактивный тип диалога пользователя с ССМ на функциональном языке позволяет программировать ход решения задач аналогично аналитическому, когда каждый этап исследований можно верифицировать и, постепенно выполняя ряд операций, получить решение задачи. Такое программирование, осуществляемое на входном языке ССМ, требует не только задание операторов и функций, но и программное описание параметров, процедур, правил преобразования и т.п. в другом – функциональном – виде и соответствует основным принципам традиционных форм программирования.

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

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

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

Интерактивные многофункциональные компьютерные ССМ стали универсальными математико-информационными средами; причем ССМ возможно использовать на основе как встроенных инструментов и функций, так и операторов и расширений, задаваемых самим пользователем, наращивая, тем самым мощь математических операций. Поэтому разработка теоретических основ создания комплексов программ на основе аналитико-численных алгоритмов в ССМ (MAPLE, MATHEMATICA, MATHCAD) и их практическое воплощение для решения научных и прикладных задач ядерной энергетики и физики частиц представляет собой актуальную комплексную задачу.

АКТУАЛЬНОСТЬ ПРИКЛАДНЫХ ЗАДАЧ Одним из важных направлений использования ядерной энергии является использование реакторов малой и сверхмалой мощности (до 100 МВт и 300 кВт соответственно). Такие реакторы могут найти широкий спектр применения для автономного теплоснабжения и электроснабжения, в удаленных и труднодоступных районах и т.п. Создание комплекса программ в среде ССМ и проведение вычислительного эксперимента с целью обоснования возможности разработки таких реакторов с саморегулированием является актуальной научной проблемой.

Актуальной комплексной проблемой является также создание библиотек нейтронных активационных ядерных данных по сечениям нейтронных реакций с ядрами мишеней от Z = 184 в энергетическом диапазоне от 150 МэВ до 1 ГэВ, поскольку такие данные необходимы для изучения свойств конструкционных материалов ядерной техники. Среди большого числа задач, решаемых при создании библиотек, важным этапом является сравнительный анализ моделей ядерных данных по их соответствию с экспериментом. Поскольку исходные экспериментальные данные сильно различаются как по количеству для разных реакций, так и по степени их однородности внутри каждого диапазона (причем не все исходные расчетные данные покрывают соответствующие области экспериментальных данных), появляется необходимость в создании комплексов программ в ССМ для комплексного анализа данных и ранжирования моделей на основе как численных расчетов с применением внутренних ресурсов ССМ, так и дополнительных процедур по определению и визуализации разных типов промежуточных и конечных результатов.

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

Многочисленные данные о сечениях деления свидетельствуют наличии двух и более локальных максимумах потенциального барьера деления для ядер элементов c массовым числом A = 229247. Решение задачи о вычислении коэффициента прохождения для таких потенциалов представляет важную научную и прикладную проблему в области ядерной энергетики и физики частиц. Такая задача является сложной в математическом плане и для решения требует разработки соответствующих алгоритмов и программ. Кроме того, необходимы значительные вычислительные ресурсы для визуализации и получения численных значений характерных физических параметров ядер, что естественным образом, возможно реализовать с помощью комплексов программ в ССМ.

Задачи теории поля в пространстве де Ситтера длительное время привлекают внимание исследователей, а в последние годы интерес к таким задачам был актуализирован появлением теории инфляционной Вселенной. Оказывается, что задача исследования полей в пространстве де Ситтера является более ясной при использовании статической системы координат и позволяет рассчитывать квантовые процессы аналитически. Поэтому создание комплекса программ в ССМ для использования формализма Ньюмена-Пенроуза как инструмента вычислений в искривленном пространстве-времени и его применение для решения уравнений безмассовых полей со спином в пространстве де Ситтера является актуальной задачей физики частиц.

ЦЕЛЬ РАБОТЫ Разработка комплексов проблемно-ориентированных программ на основе аналитико-численных алгоритмов в системах символьной математики и их практическое применение для решения прикладных задач ядерной энергетики и физики частиц.

ТЕОРЕТИЧЕСКИЕ ПРОБЛЕМЫ создания комплексов проблемно-ориентированных программ – многопрофильная программная задача, заключающаяся в:

– алгоритмизации математически поставленных задач средствами встроенных математических ресурсов ССМ в сочетании со специфическим для них программным заданием физических и математических данных и дополнительных условий;

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

– сочетании аналитических вычислений и технологии вычислительного эксперимента;

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

ПРИКЛАДНЫЕ НАУЧНЫЕ ПРОБЛЕМЫ – использование комплексов программ для решения прикладных задач ядерной энергетики и физики частиц, связанных с:

– разработкой и реализацией моделей физических и энергетических объектов на основе аналитико-численных алгоритмов в ССМ;

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

ПРАКТИЧЕСКАЯ ЦЕННОСТЬ РЕЗУЛЬТАТОВ РАБОТЫ 1. Разработаны комплексы программ для практического применения компьютерных технологий при решении прикладных задач ядерной энергетики и физики частиц:

– показана принципиальная возможность разработки реактора малой мощности с саморегулированием;

– обосновано использование моделей CASCADE, ISABEL, CEM2K, INCL4, Dresner, Bertini для расчета активационных нейтронных данных для нуклидов с Z = 684 в области энергий 1501000 МэВ;

– обоснованы рекомендации по изотопному составу тепловыделяющего элемента и условиям хранения многослойных устройств с ТВЭЛ;

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

2. Обосновано использование ССМ для обеспечения учебного процесса информационно-инновационного типа и внедрены в ИАТЭ учебно-исследовательские методики учебного процесса информационного типа, основанного на использовании современного программного обеспечения и обеспечивающего преемственность и фундаментальность образования. Разработанные программные комплексы используются на лекциях, практических и лабораторных занятиях (в виде компьютерного практикума моделирующего типа с использование ССМ) по курсам элементарной, общей, теоретической и математической физики; они являются универсальными учебноисследовательскими ресурсами, применяемыми студентами и аспирантами при выполнении учебно-исследовательских курсовых и дипломных работ и диссертаций.

НАУЧНАЯ НОВИЗНА РАБОТЫ 1. Показана возможность использования ССМ для решения теоретических и прикладных задач различного математического уровня и создания на их основе комплексов программ научного и учебно-исследовательского типов.

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

3. На основе символьных и аналитико-численных алгоритмов разработаны комплексы прикладных программ для:

– реализации, вычислительного эксперимента и проверки адекватности математической модели реактора теплоснабжения малой мощности с саморегулированием;

– обработки и статистического анализа данных для нейтронной библиотеки активационных файлов «IEAF-2005» с целью выбора наилучшей модели взаимодействия протонов с ядрами мишеней;

– реализации и расчета серии моделей, учитывающих различные физикотехнические характеристики, многослойного устройства с ТВЭЛ: линейная и нелинейная модели теплопроводности; учет контактного термического сопротивления, излучения и различных граничных условий;

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

– использования формализма Ньюмена-Пенроуза в искривленном пространстве-времени с целью исследование безмассовых полей.

4. Получен ряд новых аналитических и численных результатов при решении прикладных задач ядерной энергетики и физики частиц:

– выявлены оптимальные режимы работы реактора теплоснабжения малой мощности для нормального и «тлеющего» режимов;

– выполнено ранжирование моделей взаимодействия протонов с ядрами мишеней для нейтронной библиотеки активационных файлов «IEAF-2005»;

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

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

ДОСТОВЕРНОСТЬ РЕЗУЛЬТАТОВ Достоверность результатов работы обеспечивается:

– использованием строгих математических выкладок и фундаментальных физических законов при решении прикладных задач;

– многократным тестированием разработанных алгоритмов и элементов программных комплексов на известных классических задачах теоретической и математической физики;

– сопоставлением с известными промежуточными результатами и соответствием частным результатам, полученным в работах других авторов;

– применением встроенных в ССМ средств проверки промежуточных и конечных результатов;

– полным соответствием результатов решения ряда задач, полученных использованием разных ССМ;

– возможностью проверки ряда аналитических результатов, допускающих ясную физическую и математическую интерпретацию;

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

АПРОБАЦИЯ РАБОТЫ Основные результаты диссертации докладывались на 14 российских и между народных научных конференций:

– 6-ой Всесоюзной конференции «Современные теоретические и экспериментальные проблемы теории относительности и гравитации». МГПИ, Москва, 1985;

– сессии Ядерного отделения АН СССР, ФИАН, Москва, 1985;

– Всесоюзном симпозиуме «Гравитация и объединение фундаментальных полей», Киев, 1985;

– межреспубликанской школе «Квантовые эффекты в сильных полях», Кишинев, 1985;

– 7-ой Всесоюзной конференции «Современные теоретические и экспериментальные проблемы теории относительности и гравитации». ЕГУ, Ереван, 1988;

– международной школе «Manufacturing Opportunities through Science and Technology», USA, ORNL, Oak Ridge, Tennessee, 1994;

– международной конференции «International Science Partners Program», USA LLNL, Livermore, California, 1994;

– международной конференции «International Science Partners Program», Los Alamos, New Mexico, 1995;

– 3-ей международной конференции «Математические идеи П.Л. Чебышёва и их приложение к современным проблемам естествознания». Обнинск, 2006 г.;

– международной конференции «Физика в системе инженерного образования стран ЕврАзЭС». Москва, 2006.

– X международной конференции «Безопасность АЭС и подготовка кадров». Обнинск, 2007 г.;

– 4-ой международной конференции «Математические идеи П.Л. Чебышёва и их приложение к современным проблемам естествознания». Обнинск, 2008 г.;

– международной конференции «On Nuclear Data for Science and Technology ND 2007», Nice, France, 2007 г.;

– международном семинаре « First Workshop on Accelerator Radiation Induced Activation », Paul Scherrer Institut, Switzerland, 2008.

В 2002 и 2004 г.г. работы автора по использованию прикладных математических пакетов были признаны лучшими, а автор признан победителем конкурса, проводимого Компанией SoftLine. Шесть учебных пособий автора, разработанных по материалам диссертации, прошли экспертизу учебно-методического управления по физике и им присвоен гриф УМО.

ЛИЧНЫЙ ВКЛАД СОИСКАТЕЛЯ Результаты работы, выносимые на защиту, получены лично соискателем. Из работ, выполненных в соавторстве, в диссертацию включены результаты, которые получены лично соискателем или при его непосредственном участии.

СПИСОК НАУЧНЫХ ПУБЛИКАЦИЙ Основные результаты диссертации опубликованы в 26 статьях (включая 7 работ в журналах из списка ВАК) и докладах и представлены в 2 научных монографиях, 14 учебных пособиях, а также в 13 электронных работах, опубликованных на официальных сайтах производителей систем символьной математики.

Автором опубликовано также 56 работ в электронном виде в ИНТЕРНЕТ (в форме рабочих листов прикладных пакетов и электронных книг).

СТРУКТУРА И ОБЪЕМ ДИССЕРТАЦИИ Диссертационная работа изложена на 286 страницах машинописного текста и состоит из введения, шести глав, заключения и списка литературы из 218 наименований работ.

СОДЕРЖАНИЕ РАБОТЫ

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

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

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

Проводится сравнение реализаций программных решений задач в стандартных подходах и подходах, основанных на использовании ССМ (рис. 1.1).

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

Рис. 1.2. Концепция создания комплексов в ССМ Таблица 1.Основные типы операций при реализации функциональных алгоритмов в ССМ Тип операции Содержание Выход Уравнение в формате 1. Ввод уравнений Программная запись уравнений ССМ.

Программный ввод параметров и данных.

2. Ввод данных Программная запись начальных условий.

Параметры и уравнения Программная запись дополнительных условий. в формате ССМ.

Преобразования переменных.

3. Использование Вывод преобразованинструментов Преобразования координат.

ных уравнений и вырапреобразования Преобразования систем уравнений. жений.

Аналитическое решение алгебраических или 4. Использование дифференциальных уравнений или систем инструментов реЧисленное решение алгебраических или диффе- Вывод решений уравшения уравнений ренциальных уравнений или систем уравнений нений.

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

6. Использование Ввод визуализируемых функций.

инструментов виВывод результатов в виде графиков.

зуализации Вывод графиков Рис. 1.3. Классическая триада математического моделирования и ее реализация средствами ССМ К типам информации относятся: основная (информация для инициализации ресурсов пакетов и о решаемой задаче); дополнительная (информация для управления процессом вычислений и вспомогательная информация о системе); рабочая (информация об использовании конкретной ССМ и данные о выводе результатов).

Показывается адекватность использования ССМ классической триаде математического моделирования (рис. 1.3). Обосновывается возможность использования ССМ для обеспечения учебного процесса информационно-инновационного типа.

В разделе 1.3 формулируются программно-информационные основы создания комплекса программ для решения уравнений в частных производных и разработаны два программных комплекса, которые состоят из трех программ (рис. 1.4).

1. Программный комплекс «Уравнения параболического и гиперболического типов»;

2. Программный комплекс «Уравнения эллиптического типа», Рис. 1.4. Структура программных комплексов Программа 1 предназначена для построения формальных решений уравнений в частных производных различными методами. Выходные данные этой программы представляются в виде формальных рядов или интегралов, представляющих формальное решение, и используются затем в программе 2. Программа 2 предназначена для решения краевых задач с использованием полученных формальных решений.

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

В частности с помощью программного блока «Уравнения эллиптического типа» решается следующие краевые задачи:

Краевые условия Краевая задача Первое граничное условие Второе граничное условие На границе (радиуса R1) круга:

I. Первая краевая зада(А) В задаче внутри круга: ограu r, f ча для круга ниченность в начале координат:

rRu r, На границе круга:

rII. Вторая краевая заru r, g дача для круга (Б) В задаче вне круга: ограниченrRность на бесконечности:

На границе круга:

III. Третья краевая заu r, r u r, ru r, rR h дача для круга На границе (радиуса R1) кольца: На границе (радиуса R2) кольца:

IV. Первая краевая заu r, f1 u r, f2 дача для кольца rR1 rR В разделе 1.4 формулируются программно-информационные основы создания комплекса программ для использования векторного анализа, тензорной алгебры и тензорного анализа и разработаны три программных комплекса.

1. Программный комплекс «Векторный анализ», который состоит из двух программ (рис. 1.5).

Программа 1 предназначена для введения необходимых координатных систем, задания скалярных и векторных операторов, а также программного вычисления градиента, дивергенции, ротора и оператора Лапласа для этих полей с использованием ресурсов пакетов MAPLE или MATHEMATICA. Программа 2 предназначена для задания пользователем необходимых координатных систем (например, невстроенных в системы) и вычисления градиента, дивергенции, ротора и оператора Лапласа для этих полей. Программа представляет собой рабочие листы пакетов MAPLE, MATHEMATICA или MATHCAD.

Рис. 1.5. Структура программного комплекса «Векторный анализ» 2. Программный комплекс «Тензорная алгебра», который состоит из двух программ (рис. 1.6).

Рис. 1.6. Структура программного комплекса «Тензорная алгебра» Программа 1 предназначена для программного задания векторов и тензоров в четырехмерном псевдоевклидовом пространстве-времени. Выходные данные этой программы представляются в виде векторных и тензорных объектов пакета и могут использоваться в программе 2. Программа 2 предназначена для выполнения операций с векторами и тензорами в четырехмерном псевдоевклидовом пространствевремени. Выходные данные этой программы представляются в виде преобразованных векторных и тензорных объектов и их функций.

3. Программный комплекс «Тензорный анализ», который состоит из двух программ (рис. 1.7).

Программа 1 предназначена для задания тензорных объектов в 4-мерном искривленном пространстве-времени. Выходные данные этой программы представляются в виде тензорных объектов пакета и могут использоваться в программе 2. Программа 2 предназначена для вычисления тензорных объектов в конкретных физически важных пространствах и выполнения над ними операций. Выходные данные программы представляются в виде тензорных объектов и их функций.

Рис. 1.7. Структура программного комплекса «Тензорный анализ» В разделе 1.5 формулируются программно-информационные основы создания комплексов для решения задач электродинамики и разработаны четыре программных комплекса.

1. Программный комплекс «Общие компьютерные методы в электродинамике», который состоит из двух программ (рис. 1.8).

Рис. 1.8. Структура программного комплекса «Общие компьютерные методы в электродинамике» Программа 1 «Поля движущегося электрического заряда» предназначена для вычисления 4-потенциала и поля движущегося электрического заряда. Выходные данные этой программы представляются в виде формул и компонент векторов и могут использоваться в программе 2. Программа 2 «Визуализация электрического и магнитного полей» предназначена для визуализации электрического и магнитного полей с помощью различных типов графиков. Выходные данные этой программы представляются в виде двумерных и трехмерных графиков (в том числе, поверхностных и векторных). Заметим, что выходные данные программы 2 имеют универсальный характер и могут быть использованы во всех случаях, где необходимо выполнить визуализацию данных в виде различных графиков.

2. Программный комплекс «Тензоры в релятивистской электродинамике», который состоит из двух программ (рис. 1.9).

Рис. 1.9. Структура программного блока «Тензоры в релятивистской электродинамике» Программа 1 «Программное задание векторов и тензоров электродинамики» предназначена для задания векторов и тензоров релятивистской электродинамики в четырехмерном псевдоевклидовом пространстве-времени. Выходные данные этой программы представляются в виде векторных и тензорных объектов пакета и могут использоваться в программе 2. Программа 2 «Преобразования векторов и тензоров электродинамики» предназначена для выполнения преобразований векторов и тен зоров релятивистской электродинамики в четырехмерном псевдоевклидовом пространстве-времени при переходе от одной инерциальной системы отсчета к другой.

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

3. Программный комплекс «Уравнения в частных производных в электродинамике», который состоит из четырех программ (рис. 1.10).

Программа 1 «Стационарные поля: двумерные эллиптические уравнения в электродинамике» предназначена для решения двумерных (или с цилиндрической симметрией) стационарных краевых задач электродинамики на основе построенных формальных решений эллиптических уравнений. Выходные данные этой программы представляются в виде формул и функциональных соотношений, представляющих физические величины, описывающие поля и заряды. Программа 2 «Стационарные поля: трехмерные эллиптические уравнения в электродинамике» предназначена для решения трехмерных краевых стационарных задач электродинамики на основе программно построенных формальных решений эллиптических уравнений. Выходные данные этой программы представляются в виде формул и функциональных соотношений, представляющих физические величины, описывающие поля и заряды.

Программа 3 «Квазистационарные поля: параболические уравнения в электродинамике» предназначена для решения краевых квазистационарных задач электродинамики на основе программно построенных формальных решений параболических уравнений. Выходные данные этой программы представляются в виде формул и функциональных соотношений, представляющих физические величины, описывающие поля и токи. Программа 4 «Численное интегрирование уравнений движения заряженных частиц» предназначена для решения численного решения уравнений движения заряженных частиц в электромагнитных полях. Выходные данные этой программы представляют собой наборы численных данных и/или графики, представляющие зависимости координат и скоростей от времени, а также траекторию движения.

Рис. 1.10. Структура программного комплекса «Уравнения в частных производных в электродинамике» 4. Программный комплекс «Интегрирование уравнений движения заряженных частиц», который состоит из двух программ (рис. 1.11).

Программы 1 и 2 («Интегрирование уравнения Гамильтона-Якоби») предназначены для интегрирования уравнений Гамильтона-Якоби, описывающих движения классических и релятивистских заряженных частиц соответственно (S – действие):

22 S S S r r rU r,, S, 2 m 2 m r2 2 m r2 sin2 t 22 S 1 S 1 S m2 c4 c2 U r,, S, t r r2 r r2 sin2 r Выходные данные этой программы представляются в виде формул и функциональных соотношений, представляющих законы движения и траектории частиц.

Программы 3 и 4 («Интегрирование уравнений Лагранжа») предназначены для интегрирования уравнений Лагранжа, описывающих движения классических и релятивистских заряженных частиц в постоянных однородных электромагнитных полях соответственно:

dv e0 em e0 E v H, e0 E v H.

dp dt c dt c Рис. 1.11. Структура программного комплекса «Интегрирование уравнений движения заряженных частиц» Выходные данные этой программы представляются в виде формул и функциональных соотношений, описывающих законы движения и траектории частиц.

Рис. 1.12. Структура программного комплекса «Теория поля в искривленном пространстве-времени» В разделе 1.6 формулируются программно-информационные основы создания комплексов для решения прикладных задач теории поля в искривленном простран стве-времени и разработан:

Программный комплекс «Теория поля в искривленном пространстве-времени» (рис. 1.12).

Программа 1 «Инварианты тензора кривизны» предназначена для вывода инвариантов тензора Римана 1 I1 Riklm Riklm Riklm Riklm, I2 Riklm Rlmpr Rprik i Riklm Rlmpr Rprik 48 и их вычисления для конкретных пространств. Выходные данные этой программы представляются в виде функциональных соотношений для инвариантов.

Программа 2 «Решение уравнений гравитационного поля» предназначена для решения уравнений гравитационного поля для некоторых физически важных пространств. Выходные данные этой программы представляются в виде функциональных соотношений для компонентов метрического тензора:

ds2 e t,r c2 dt2 et,r dr2 r2 d r2 sin d или 2.

ds2 c2 dt2 a t sin d r2 sin d d 2 r Программа 3 «Вектор поля гравитационно-инерциальных сил и визуализация искривленного пространства» предназначена для визуализации физически важных искривленных пространств на основе вычисления вектора гравитационноинерциальных сил (G – ненулевые символы Кристоффеля) Gi,F Fi Fj hij, Fi gи параметров эргосферы черной дыры. Программа представляет собой рабочие листы пакетов MAPLE и MATHCAD. Выходные данные этой программы представляются в виде разных типов графиков.

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

1. Программный комплекс «Одномерное квантовое рассеяние в MAPLE и MATHEMATICA» (рис. 1.13).

Программа 1 «Решение уравнения Шредингера» предназначена для решения уравнения Шредингера с заданными потенциальными функциями, его анализа и выделения асимптотик волновых функций. Выходные данные этой программы представляются в виде функциональных соотношений для волновых функций и используются в программе 2. Программа 2 «Потоки вероятности, коэффициента прохождения и отражения» предназначена для вычисления физических параметров рассеяния.

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

Рис. 1.13. Структура программного блока «Одномерное квантовое рассеяние в MAPLE и MATHEMATICA» 2. Программный комплекс «Одномерное квантовое рассеяние на прямоугольных барьерах в MATHCAD», который состоят из трех программ (рис. 1.14).

Рис. 1.14. Структура программного блока «Одномерное квантовое рассеяние на прямоугольных барьерах в MATHCAD» Программа 1 «Сшивание волновых функций» предназначена для «сшивания» волновых функций, являющихся решением уравнения Шредингера с постоянным значением потенциальной функции. Выходные данные этой программы представляются в виде функциональных соотношений для коэффициентов сшивки и используются в программе 2. Программа 2 «Потоки вероятности, коэффициента прохождения и отражения» предназначена для вычисления физических параметров рассеяния.

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

3. Программный комплекс «Квантовые связанные системы» (рис. 1.15);

Программа 1 «Решение уравнения Шредингера и анализ квантовых состояний» предназначена для решения уравнения Шредингера с заданной потенциальной функцией. Программа представляет собой рабочие листы в пакетах MAPLE и MATHEMATICA. Выходные данные программы представляются в виде функциональных соотношений для волновых функций и используются в программе 2. Программа 2 «Анализ волновых функций, квантовые состояния» предназначена для анализа волновых функций и определения квантовых состояний. Выходные данные этой программы представляются в виде формул и функциональных соотношений, представляющих нормированные волновые функции и энергетические уровни и используются в программе 3.

Рис. 1.15. Структура программного комплекса «Квантовые связанные системы» Программа 3 «Численные расчеты и визуализация полученных результатов» предназначена для получения численных значений физических квантовых состояний и их визуализации. Выходные данные этой программы представляются в виде наборов численных данных и графиков различных типов.

ВО ВТОРОЙ ГЛАВЕ разрабатывается комплекс программ для расчета и выполнения на его основе вычислительного эксперимента в модели реактора теплоснабжения малой мощности.

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

В разделе 2.2 приводится программная реализация комплекса программ в системе MATHCAD (рис. 2.1); описываются стационарные состояния работы реактора;

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

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

Точность вычисления переменных: TOL=0.001 (Convergence Tolerance). Численное решение систем: (блок «Given–Find», критерий точности для дополнительных условий (Constraint Tolerance) CTOL=0.001.

Система уравнения нестационарной теплогидравлики – жесткая система дифференциальных уравнений; ее решение осуществляется в MATHCAD с помощью встроенных функций (блок «Given–Find») по алгоритму RADAUS с аналогичными параметрами точности. Точность вычисления производных (метод Риддера) TOL=0.001; точность решения до 10-6.

Комплекс программ для расчета модели реактора состоит из трех программ и двух программных блоков (рис. 2.2).

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

Рис. 2.2. Структура комплекса программ для расчета модели реактора Программа 3 предназначена для численного решения системы уравнений нестационарной теплогидравлики реактора. Математически эта система представляет собой систему нелинейных алгебраических и дифференциальных уравнений. Выходные данные этой программы, полученные при заданных начальных условиях работы реактора (характеризующих некоторое стационарное состояние), представляются операторами MATHCAD, возвращающими решения системы уравнений; они могут быть представлены в виде наборов численных данных для параметров реактора (для определенных моментов времени) и в виде графиков зависимостей параметров от времени. Программный блок 1 предназначен для моделирования переходных процессов, возникающих в системе в результате воздействия внешних факторов:

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

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

1. Система уравнений стационарной теплогидравлики реактора.

2. Система уравнений нестационарной теплогидравлики реактора.

Математиче- 3. Функциональные соотношения, которые моделируют переходные процессы, ская модель воздействующие на систему.

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

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

2. Изменение температуры воды на входе во второй контур (моделирование сеИсследование зонного изменения температуры) и возвращение в исходное состояние.

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

4. Изменение начальной реактивности реактора (моделирование факторов, меняющих начальную реактивность) и возвращение в исходное состояние.

Метод функционального программирования в среде прикладного пакета Вычислитель- MATHCAD с использованием численно-аналитических алгоритмов решения ный алгоритм системы нелинейных алгебраических уравнений и системы нелинейных дифференциально-алгебраических уравнений.

Компьютерная Комплекс программ (рабочие листы) в среде MATHCAD на входном языке.

программа Компьютер Персональный компьютер.

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

В разделе 2.3 реализуется вычислительный эксперимент (рис. 2.3, табл. 2.1) и проверяется адекватность математической модели реактора теплоснабжения.

Проведенный вычислительный эксперимент (по изучению саморегулируемых режимов работы реактора и переходных процессов, возникающих в системе в результате воздействия внешних факторов), интерактивный анализ режимов работы и визуализация базовых параметров реактора позволяют выявить оптимальные режимы работы реактора (в том числе и для некоторых критических ситуаций). Тем самым обосновывается принципиальная возможность разработки реактора малой мощности с саморегулированием В ТРЕТЬЕЙ ГЛАВЕ разрабатывается комплекс программ для сравнительного анализа моделей ядерных данных (Bertini/ABLA; Bertini/Dresner; CEM2K;

NCL4/ABLA; INCL4/Dresner; ISABEL/ABLA; ISABEL/Dresner; CASCADE), представленных в нейтронной библиотеке активационных файлов «IEAF-2005» (Intermediate Energy Activation File) для 682 нуклидов с данными по сечениям нейтронных реакций с ядрами мишеней с Z = 184 в энергетическом диапазоне от 150 МэВ1 ГэВ.

В разделе 3.1 разрабатываются программные процедуры подготовки данных для сравнительного анализа, основанные на экспертных оценках, количественно характеризуемых факторами согласия D1 0.4, D2 0.8, где icalc iexp Dil , iexp expil и calcil – экспериментальные и расчетные сечения для i-ого экспериментального значения и реакции l; и с помощью полиномиальной регрессии. В результате экспериментальные (примерно для тысячи реакций) и теоретические данные представлялись в виде:

exp exp E1l 1l 1l X1l Y1exp Y1exp l l ..................

, EXPl Eil iexp iexp Xil Yilexp Yilexp l l ..................

exp exp X YIexp YIexp l I l l EI l I l I l где l = 1, …, L – номер реакции, I = 1, …, I – количество экспериментальных точек соответствующих ей, Eil – энергия (в МэВ) налетающего протона, expil и expil - сечение (в мб) и погрешность сечения выхода остаточного ядра;

Yj 1l Y0 1l Y0 1l X1l X ..., Y0 ... Y0 ..., Yj ..., , X Y0 I l Y0 I l Yj I l I l где X – вектор энергий, Y0 и Y0 – векторы экспериментальных данных с сечениями протонных реакций и их погрешностей соответственно, Yj (j нумерует анализируе мые модели: Bertini/ABLA (j = 1); Bertini/Dresner (j = 2); CEM2K (j = 3);

INCL4/ABLA (j = 4); INCL4/Dresner (j = 5); ISABEL/ABLA (j = 6); ISABEL/Dresner (j = 7), CASCADE (j = 8)).

Всего было выбрано 3999 экспериментальных значений сечений реакций (p, xnypza) с ядрами Z=684 для более 1000 реакций в энергетическом диапазоне Е=1501000 МэВ. Наряду с независимыми выходами ядер в рассмотрение были включены некоторые реакции типа (p, xnypza). Данные ранжировались в порядке возрастания массового числа ядра-мишени с А=6210 и разбиты на 9 поддиапазонов.

В каждом поддиапазоне экспериментальные данным были сопоставлены данные, полученные с использованием программы CASCADE/INPE и многоцелевой программы MCNPX.

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

Анализ основывается на использовании векторов imax cvar Y0L,YjL MMQL Y0L YjL , CORL , j j ii i1 stdev Y0L stdev YjL MMQ _ CORL MMQL CORL, j jj где MMQ(L)j – вектор обратных значений сумм квадратов разностей экспериментальных и теоретических данных для каждой реакции; imax – число значений энергий для реакции L, j – номер модели (j = 1, 2,…8); COR(L)j – вектор корреляций экспериментальных и теоретических данных для каждой реакции; cvar(Y0,Yj) и stdev(Yj) – операторы для вычисления ковариации и стандартного отклонения соответственно.

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

Комплекс программ для анализа моделей ядерных данных состоит из пяти программ (рис. 3.1).

Программы 1 и 2 предназначены предварительной обработки экспериментальных и соответствующих им расчетных данных на основе факторов согласия D = 0.4 и D = 0.8. Программы написаны на языке Fortran; они представляют данные в виде единого EXCEL-файла как для всей реакций, так и в виде наборов файлов для всех реакций. Выходные данные этих программ используются в программе 4 для проведения сравнительного анализа данных для каждой реакции. Программа 3 (набор однотипных программ для каждой реакции) предназначена для предварительной обработки экспериментальных и соответствующих им расчетных данных на основе метода полиномиальной регрессии в прикладном математическом пакете MATHCAD.

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

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

Рис.3.1. Структура комплекса программ для анализа моделей ядерных данных Программы 3 – 5 представляют собой рабочие листы (и наборы однотипных рабочих листов для каждой реакции) пакета MATHCAD.

а б Рис.3.2. Визуализация ранжирования моделей с помощью трехмерных гистограмм для всего диапазона (а) и для поддиапазонов (б) массового числа ядер-мишеней В результате проведено ранжирование (для всего диапазона (рис. 3.2а) и для поддиапазонов (рис. 3.2б) массового числа ядер-мишеней) моделей взаимодействия протонов с ядрами мишеней для нейтронной библиотеки активационных файлов «IEAF-2005» для данных 2006 и 2007 г.г. и обосновано использование моделей CASCADE, ISABEL, CEM2K, INCL4, Dresner, Bertini для расчета активационных нейтронных данных для нуклидов с Z = 684 в области энергий 150–1000 МэВ.

В ЧЕТВЕРТОЙ ГЛАВЕ разрабатывается комплекс программ для реализации моделей теплопроводности и анализа температурных полей в многослойном устройстве со сферическим тепловыделяющим элементом.

В разделе 4.1 описываются устройство, его технологические реализации (отличающиеся по геометрии и мощности ТВЭЛ) и решаемые краевые задачи (табл. 4. и рис. 4.1 и 4.2); приводится структура комплекса программ в системах MAPLE и MATHEMATICA для расчета температурных полей.

Таблица 4.Модели теплопроводности в многослойном устройстве со сферическим ТВЭЛ Модель Учитываемые факторы I, Базовая модель (БМ) Стандартная модель стационарной теплопроводности.

II, Модель с учетом контактного Учет контактного термического сопротивления между термического сопротивления (МКС) оболочками.

III, Ограниченная нелинейная мо- Учет зависимости коэффициента теплопроводности от дель (ОНМ) температуры в IV оболочке.

IV, Полная нелинейная модель Учет зависимости коэффициентов теплопроводности от (ПНМ) температуры во всех оболочках.

V, Полная нелинейная модель с Учет зависимости коэффициентов теплопроводности от учетом контактного термического температуры во всех оболочках и контактного термичесопротивления (ПНК) ского сопротивления.

VI, Модель, с учетом излучения на Учет условия излучения на внешней границе устройства.

внешней границе (МИГ) VII, Модель с учетом излучения на Учет условия излучения в прослойке между сферическим ТВЭЛ (МИТ) ТВЭЛ и отражателем.

VIII, Модель с учетом внешнего Учет внешнего воздействия (алюминиевые стержни для воздействия (МВВ) усиления отвода тепла).

Рис. 4.2. Геометрические характеристики устройства для трех классов технологий (радиусы оболочек указаны в см) Рис. 4.1. Объект исследования и модели Расчет температурных полей – важная задача для определения условий хранения ядерных устройств., которые должны удовлетворять требованиям «нераспространения» и технической реализуемости (не слишком высокая температура). Подход к решению – денатурирование ТВЭЛа: изменение его изотопного состава (в т.ч. как фактор защиты от несанкционированного распространения ядерных материалов).

Комплекс состоит из совокупности программ, схематически представленных на рис. 4.3. Более подробно на рис. 4.4. приведена часть комплекса для расчета базовой (БМ) и ограниченной нелинейной (ОНМ) моделей.

Рис. 4.3. Структура комплекса программ Рис. 4.4. Часть структуры комплекса для расчета базовой и ограниченной нелинейной моделей Программы 1 и 2 предназначены для решения первой и третьей краевых задач для температурного поля в рамках линейной модели. Выходные данные этих программ представляются в виде функциональных зависимостей температуры от радиальной координаты и соответствующих двумерных и трехмерных графиков и используются затем в программе 5.

Программы 3 и 4 предназначены для получения решений первой и третьей краевых задач для температурного поля в рамках нелинейной модели. Выходные данные этих программ представляются в виде функциональных зависимостей температуры от радиальной координаты и соответствующих двумерных и трехмерных графиков и используются затем в программе 5.

Программа 5 предназначена для сравнения результатов, полученных в рамках линейной и нелинейной моделей: результаты представляются в виде графиков и численных значений температур (в том числе и для характерных точек).

В работе рассматриваются два типа краевых условий:

1) на границе устройства поддерживается постоянная температура Text (первая краевая задача):

u5 r Text, rrгде Text – постоянная температура, поддерживаемая на границе устройства.

2) на границе устройства происходит теплообмен (с постоянным тепловым потоком) с внешней средой, температура которой равна Text (третья краевая задача):

ru5 r h Text u5 r, rrrrгде h – коэффициент теплообмена.

В разделе 4.2 рассчитываются тепловые поля в рамках линейных моделей.

1) Реализуется базовая линейная (с постоянными коэффициентами теплопроводности) модель стационарной теплопроводности. Уравнения теплопроводности для сферически-симметричных I ( –объемная тепловая плотностью ТВЭЛа) и II–V областей записываются (i = 2, 3, 4, 5) соответственно:

2u1 r u1 r 2ui r ui r 2 2 k1 0, 0.

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

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

В разделе 4.3 рассчитываются тепловые поля в рамках нелинейных моделей.

1) Реализуется ограниченная нелинейная (предполагающая зависимость коэффициента теплопроводности от температуры в IV слое) модель стационарной теплопроводности. При этом предполагается линейная зависимость коэффициента теплопроводности k4 материала в IV слое от температуры u4:

K4 4 u4 r k4, где 4 и k4 – постоянные параметры. При этом температурное поле u4 будет удовлетворять нелинейному уравнению второго порядка:

2u4 r u4 r u4 r 2 a 0.

4 u4 r k4 r2 r r r Полученные аналитические решения представляют собой кусочно-гладкие функции, описывающие зависимость температуры от радиальной координаты r (в рамках обеих моделей) для первой и третьей краевых задач соответственно. Решения визуализируется с помощью двумерных и трехмерных графиков (рис. 4.5).

Рис. 4.5. Двух- и трехмерная визуализации температурных полей в ограниченной нелинейной модели 2) Реализуется полная нелинейная модель (предполагающая зависимость коэффициента теплопроводности от температуры во всех слоях) модель стационарной теплопроводности. При этом предполагается линейная зависимость коэффициентов теплопроводности ki материалов от температуры ui:

Ki i ui r ki, где i и ki – постоянные параметры. При этом температурное поле u1 будет удовлетворять нелинейному уравнению второго порядка:

2u1 r u1 r u1 r 2 1 u1 r k1 0, 1 r2 r r r а остальные температурные поля будут удовлетворять нелинейному уравнению, приведенному для ограниченной нелинейной модели (как в IV слое). Для определения коэффициентов линейной зависимости от температуры были использованы экспериментальные справочные данные и программа в системе MATHCAD, реализующая метод линейной регрессии.

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

В разделе 4.4 рассчитываются тепловые поля в рамках моделей, учитывающих тепловое излучение и внешнее воздействие 1) Реализуется полная нелинейная модель с учетом теплового условия излучения на внешней границе устройства:

u5 r 4 h Text u5 r H text u5 r, rrrrr rrгде text –температура окружающего помещения, где H – коэффициент, характеризующий излучательные свойства поверхности устройства.

2) Реализуется модель, построенная на основе базовой и учитывающая условие теплового излучения в прослойке между сферическим ТВЭЛ и отражателем.

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

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

k4,i k4 i k3, i – номер подоболочки, i – коэффициент, зависящий от соотношения объема стержней и подоболочки:

n Vi i , V4,i n Vi где V4,i – объем подоболочки; Vi – объем стержня, проходящего через эту подоболочку, n – число стержней.

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

Полученные результаты в виде аналитических зависимостей температурных полей визуализируются, численно анализируются; и на основании полученных данных делаются выводы по изотопному составу ТВЭЛ и физически приемлемым условиям хранения устройств для трех классов технологий и удовлетворяющих требованиям нераспространения делящихся материалов.

В ПЯТОЙ ГЛАВЕ разрабатывается комплекс программ для реализации моделей потенциальных барьеров деления тяжелых ядер; аналитического решения уравнения Шредингера; вычисления коэффициентов прохождения (s-состояния), энергии резонансных уровней и квазистационарных состояний для различных реализаций потенциалов.

Рис. 5.1. Структура комплекса программ в рамках модели I В разделе 5.1 разрабатывается модель I барьера многопрофильными прямоугольными потенциалами:

0 x 0, x a3,V01 0 x a VI x , V a1 x a2, V03 a2 x a где V0i – постоянные, определяющие «высоты» прямоугольных барьеров, ai – размерные параметры, определяющие их границы.

Для построения моделей потенциальных барьеров и решения задачи рассеяния в рамках модели I используется комплекс программ в среде прикладного математического пакета MAPLE. Комплекс программ состоит из трех программ (рис. 5.1).

Программа 1 предназначена для «сшивания» волновых функций, которые являются решением уравнения Шредингера с постоянной или равной нулю потенци альной энергией. Выходные данные этой программы представляют собой набор постоянных коэффициентов, получаемых при решении системы уравнений и используются затем в программе 2.

Программа 2 предназначена для вычисления потоков вероятности i j x x x x x x x 2 для падающей и прошедшей волн и определения коэффициента прохождения T jtr x / jin x с использованием алгебраических и дифференциальных соотношений. Выходные данные этой программы представляют собой функциональные зависимости и используются затем в программе 3. Программа 3 предназначена для визуализации и получения численных значений физических параметров.

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

0, x 0, x a , V x i V 1 1 i2 x xi, ai1 x ai; i 1,2,3; a0 0i где V01 и V03 – локальные максимумы (в точках x1 и x2 соответственно), а V02 – ло2 кальный минимум (в точке x2) потенциальной энергии; величины 1, 2 и 3 – описывают кривизны парабол; – эффективная масса частиц, которая определяется массовым числом A: 0.0540 A5/3 2 / Мэв.

Для построения моделей потенциальных барьеров и решения задачи рассеяния в рамках модели II используется комплекс программ в среде прикладного математического пакета MAPLE. Комплекс программ состоит из пяти программ (рис. 5.4).

В разделе 5.3 разрабатывается модель III барьера пятью квадратичными функциями (рис. 5.3):

0, x 0, x a , V x i V 1 1 i2 x xi, ai1 x ai; i 1,2,3,4,5; a0 0;V01 V05 0i где V02 и V04 – локальные максимумы (в точках x2 и x4 соответственно), а V03 – локальный минимум (в точке x3) потенциальной энергии.

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

Рис. 5.4. Структура комплекса программ в рамках моделей II и III Рис. 5.5. Визуализации семейства барьеров и коэффициента прохождения в модели III Программа 1 предназначена для моделирования потенциала распада с помощью трех (для модели II) или пяти (для модели III) квадратичных функций. Выходные данные этой программы представляют собой набор параметров парабол и координат точек их гладкого сшивания и используются затем в программах 2 и 5. Программа 2 предназначена для программного решения уравнения Шредингера для двух типов областей (с нулем и квадратичным потенциалами). Выходные данные этой программы представляют собой функциональные выражения для волновых функций и используются затем в программе 3.

Рис. 5.6. Зависимость ln(T(E)) Рис. 5.7. Зависимость ln(T(E)) (подпороговые резонансы) (квазистационарные уровни) Программа 3 предназначена для вычисления потоков вероятности и определения коэффициента прохождения с использованием алгебраических и дифференциальных соотношений. Выходные данные этой программы представляют функциональные зависимости и используются затем в программе 4. Программа 4 предназначена для вычисления потоков вероятности и определения коэффициента прохожде ния с использованием алгебраических и дифференциальных соотношений. Выходные данные этой программы представляют функциональные зависимости и используются затем в программе 5. Программа 5 предназначена для визуализации и получения численных значений физических параметров (рис. 5.5).

В разделе 5.4 проводится исследование результатов, полученных в рамках трех моделей, получены численные значения физических параметров на примере ядра с массовым числом A = 240: определяются уровни резонансов и соответствующие значения коэффициента прохождения для несимметричных барьеров; определяются уровни квазистационарных состояний для симметричных барьеров; исследуется влияние изменения параметров потенциала на результат (рис. 5.6 и 5.7).

В ШЕСТОЙ ГЛАВЕ разрабатывается комплекс программ (рис. 6.1) по программной реализации формализма Ньюмена-Пенроуза и применению его для исследования безмассовых полей со спином (s = 0, 1/2, 1, 3/2, 2), распространяющихся в искривленном пространстве-времени.

Комплекс состоит из двух программ (рис. 2.1).

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

В разделе 6.1 программно реализуется формализм Ньюмена-Пенроуза в искривленном пространстве-времени; строятся комплексные световые тетрады и вычисляются спиновые коэффициенты для ряда аксиально-симметричных 4-мерных искривленных пространств.

В разделе 6.2 производится разделение переменных и программными методами получены аналитические решения уравнений безмассовых полей со спином в статическом пространстве де Ситтера:

1) скалярное поле (s = 0):

Dl Dn Dn ' ' ' ' Dl Dm Dm Dm ' ' ' ' Dm R 0, где – волновая функция скалярного поля, R – скалярная кривизна;

2) электромагнитное поле (s = 1):

Dl 2 1 Dm 2 ' ' 2 '0 0, Dm 2 1 Dn 2 ' ' 0 '0 0, где 0, 1, '0 – тетрадные проекции тензора электромагнитного поля;

3) гравитационные возмущения метрики (s = 2):

Dl 2 4 1 Dm 4 ' ' 0 0, Dm 2 4 1 Dn 4 ' ' 0 0, где 0, 1 – возмущения вейлевских скаляров;

4) безмассовое спинорное поле (s = 1/2):

Dl E1 Dm ' ' E2 0, Dm E2 Dn ' ' E1 0, где E1, E2 – компоненты 4-спиноров поля;

5) поле гравитино (s = 3/2):

Dl 2 Dn 3 ' 2 ' Dm 2 ' ' Dm 3 ' 2 ' a2 0, Dn 2 ' ' ' Dl 3 2 Dm 2 ' 2 Dm 3 2 ' a1 0, где a1, a2 – компоненты гравитино.

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

Таблица 6.Элемент программного кода в MAPLE решения уравнений безмассовых полей со спином в статическом пространстве де Ситтера 3.6. Решение радиальных уравнения с помощью оператора «dsolve» Ввод > beta1=1/2*I*omega*r0; beta2=-1/2*I*omega*r0+s;

delta1=1/2*s-1/2*(s+I*omega*r0); delta2=1/2*s+1/2*(s+I*omega*r0);

EQN:=(x^2-1)*diff(f(x),`$`(x,2)) +(-2*x*s+2*delta*x+2*x+2*beta*x+*delta-2*beta)* diff(f(x),x)+(-l-2*s*beta-s-l^2+ 2*beta*delta+s^2+ delta+beta^2-2*s*delta+beta+delta^2)*f(x); dsolve(EQN);

1 x Ввод f(x)_C1 hypergeom [sl, s1l ], [s12 ], 2 ( s2 ) 1 x 1 x _C2 hypergeom [l, 1l ], [1s2 ], 2 2 2 3.7. Вывод набора решений f(x) 3.8. Вывод набора решений R(x) > R1(x):=subs([beta=1/2*I*omega*r0,delta=-1/2*I*omega*r0], Ввод x^alpha*(1+x)^beta*(1-x)^ delta*f1(x));

R2(x):=subs([beta=-1/2*I*omega*r0+s,delta=s+1/2*I*omega*r0], x^alpha*(1+x)^beta*(1-x)^delta*f1(x));

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

В ЗАКЛЮЧЕНИИ формулируются основные результаты работы, выносимые на защиту.

ОСНОВНЫЕ РЕЗУЛЬТАТЫ РАБОТЫ, ВЫНОСИМЫЕ НА ЗАЩИТУ 1. Теоретические основы создания комплексов проблемно-ориентированных программ на основе аналитико-численных алгоритмов в ССМ для решения прикладных научных задач.

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

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

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

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

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

3. Комплекс программ для обработки и статистического анализа данных нейтронной библиотеки активационных файлов «IEAF-2005».

КП разработан с целью выбора наилучшей модели взаимодействия протонов с ядрами мишеней и предназначен для обработки и статистического анализа активационных нейтронных данных для нуклидов с Z = 684 в области энергий 1501000 МэВ.

Проведено ранжирование моделей (для всего диапазона и для поддиапазонов массового числа ядер-мишеней) для данных 2006 и 2007 г.г.

4. Комплекс программ для реализации и расчета серии моделей, учитывающих различные физико-технические характеристики, многослойного устройства с ТВЭЛ.

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

различные граничные условия и специальные внешние воздействия.

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

5. Комплекс программ для реализации моделей многогорбых потенциалов деления ядер тяжелых элементов.

КП предназначен для реализации моделей потенциальных барьеров деления тяжелых ядер (А = 229 247); аналитического решения уравнения Шредингера и вычисления коэффициентов прохождения.

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

6. Комплекс программ по применению формализма Ньюмена-Пенроуза для исследования безмассовых полей в искривленном пространстве-времени.

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

ОСНОВНЫЕ РЕЗУЛЬТАТЫ ДИССЕРТАЦИИ ОПУБЛИКОВАНЫ:

в периодических научных изданиях, выпускаемых в Российской Федерации:

1. Гальцов Д.В. Морозов М.Ю., Тихоненко А.В. Безмассовые поля в статическом пространстве де Ситтера: точные решения и выбор вакуумных состояний // Теор. и мат. физика. 1988, т. 77, № 3, с. 190–203.

2. Казанский Ю.А., Левченко В.А., Матусевич Е.С., Юрьев Ю.С., Балакин И.П., Белугин В.А., Дорохович С.Л., Казанцев А.А., Тихоненко А.В., Травлеев А.А., Уваров А.А. Саморегулируемый реактор сверхмалой мощности для теплоснабжения «МАСТЕР-ИАТЭ // Известия Вузов.

Ядерная энергетика, № 3, 2003. – С. 63-71.

3. Тихоненко А.В. Комплекс программ для расчета модели реактора теплоснабжения малой мощности в прикладном математическом пакете MATHCAD // Известия Вузов. Ядерная энергетика, № 2, 2007. – С. 44-52.

4. Тихоненко А.В. Модельный анализ температурных полей в многослойном сферически-симметричном устройстве. Известия Вузов. Ядерная энергетика, № 2, 2007. – С. 110-118.

5. Конобеев А.Ю., Коровин Ю.А., Наталенко А.А., Осыкин С.А., Пильнов Г.Б., Станковский А.Ю., Тихоненко А.В. Нейтронная библиотека активационных файлов «IEAF-2005» в энергетическом диапазоне от 150 МэВ до 1 ГэВ // Известия Вузов. Ядерная энергетика, № 2, 2007. – С. 8-15.

6. Пильнов Г.Б., Тихоненко А.В. Статистический анализ данных для нейтронной библиотеки активационных файлов «IEAF-2005» // Известия Вузов. Ядерная энергетика, № 3, в.1, 2007. – С. 109-119.

7. Тихоненко А.В. Задача рассеяния для модельных потенциалов деления тяжелых ядер и комплексы программ в прикладных математических пакетах // Известия Вузов. Ядерная энергетика, № 3, в.2, 2007. – С. 110-119.

8. Морозов М.Ю., Тихоненко А.В. Квантование свободных полей в пространстве де Ситтера // Изв. АН ТССР. Сер.физ.-тех,хим., и геол.наук.1987, № 6, с.11–17.

9. Гальцов Д.В., Нуньес С.Д., Морозов М.Ю., Тихоненко А.В. Безмассовые поля в статическом пространстве де Ситтера и проблема динамического нарушения симметрии. Препринт физ. ф-та МГУ, 1986, № 4, 1987, 4 с.

10. Гальцов Д.В. Морозов М.Ю. Нуньес С.Д., Тихоненко А.В. Квантовые поля в инфляционной Вселенной // Гравитация и фунд. взаим. М. 1988. С. 68-49.

в периодическом научном электронном журнале, выпускаемом в РФ:

11. Тихоненко А.В. Решение дифференциальных уравнений в частных производных методом функционального программирования в MAPLE. 1. Уравнения параболического и гиперболического типов // Эл. журнал «Исследовано в России», 046, с. 502-511, 2007. http://zhurnal.ape.relarn.ru/articles/2007/046.pdf.

12. Тихоненко А.В. Решение дифференциальных уравнений в частных производных методом функционального программирования в MAPLE. 2. Двумерные уравнения эллиптического типа // Эл. журнал «Исследовано в России», 061, с.

652-661, 2007. http://zhurnal.ape.relarn.ru/articles/2007/061.pdf.

13. Тихоненко А.В. Тензорное исчисление и его приложения в MAPLE и MATHCAD // Электронный журнал «Исследовано в России», 069, с. 720-729, 2007.

http://zhurnal.ape.relarn.ru/articles/2007/069.pdf.

14. Тихоненко А.В. Интегрирование уравнений движения заряженных частиц методом функционального программирования в MAPLE // Эл. журнал «Исследовано в России», 089, с. 923-933, 2007.

http://zhurnal.ape.relarn.ru/articles/2007/089.pdf.

15. Тихоненко А.В. Реализация математических моделей теплопроводности и комплекс программ в Maple для расчета тепловых полей в многослойном устройстве с тепловыделяющим элементом // Эл. журнал «Исследовано в России», 25, стр. 290-300, 2008 г. http://zhurnal.ape.relarn.ru/articles/2008/024.pdf в трудах Всероссийских и международных научных конференций:

16. Морозов М.Ю. Тихоненко А.В. Безмассовые поля в конформно-конических пространствах. Материалы 7-ой Всес. конф. «Совр. теор. и эксп. проблемы теории относительности и грав.». ЕГУ. Ереван. 1988, – С. 310–311.

17. Тихоненко А.В. Информационные технологии в физическом образовании.

«Физика в системе инженерного образования стран ЕврАзЭС». Тезисы докладов. Москва: Авиаиздат, 2006. с. 272-274.

18. Тихоненко А.В. Модульная система обучения в курсе общей физики. «Физика в системе инженерного образования стран ЕврАзЭС». Тезисы докладов.

Москва: Авиаиздат, 2006. с. 272-276.

19. Тихоненко А.В. Компьютерный практикум по общей физике. «Физика в системе инженерного образования стран ЕврАзЭС». Тезисы докладов.

Москва: Авиаиздат, 2006. с. 276-278.

20. Тихоненко А.В. Прикладные математические пакеты в курсах общей и теоретической физики. «Физика в системе инженерного образования стран ЕврАзЭС». Тезисы докладов. Москва: Авиаиздат, 2006. с. 278-280.

21. Конобеев А.Ю., Коровин Ю.А., Наталенко А.А., Осыкин С.А., Пильнов Г.Б., Станковский А.Ю., Тихоненко А.В. О возможности использования каскадноэкситонных моделей CEM2K, CEM03.01, CEM03.02 и CASCAD/INPE для расчета протонных активационных данных в энергетическом диапазоне от 1МэВ до 1 ГэВ. Безопасность АЭС и подготовка кадров. X Межд. конф. Тезисы докладов. Часть 1. – Обнинск: ИАТЭ, 2007. С. 48-49.

22. Pilnov G. B., Tikhonenko A.V. Statistical analysis data for neutron library of activation files «IEAF-2005». NPP and Safety and personal training. X Int. Conf. Abstracts. Part 2. – Obninsk: INPE, 2007. P. 32.

23. Saenko A.V., Tikhonenko A.V., Artisyuk V.V. The quantitive approach to the problem of the denaturating plutonium: temperature profile in implosion device. NPP and Safety and personal training. X Int. Conf. Abs.. P. 1. – Obninsk: INPE, 2007. P. 107–108.

24. Тихоненко А.В. Реализация математических моделей теплопроводности и анализ тепловых полей в многослойном устройстве. Тр. III межд. конф. «Математические идеи П.Л. Чебышёва и их приложение к современным проблемам естествознания». Обнинск, 14-18 мая 2006 г. – Обнинск: ИАТЭ, 2008. с. 151-160.

25. Korovin Yu., Fischer U., Konobeyev A., Natalenko A., Pilnov G., Stankovskiy A., Tikhonenko A., Evaluation of activation nuclear data in the energy region 150 MeV to 1 GeV, Proc. Int. Conf. On Nuclear Data for Science and Technology ND20(2007), April 23 - 27, Nice, France, p.1175-1178.

26. Yu. Korovin, I. Kuptsov, A.Natalenko, A.Stankovskiy, A.Tikhonenko, A. Konobeyev.

Development of a New Code to Simulate Radiation Damage and Gas Accumulation in the Structural Materials of ADS. Proc. First Workshop on Accelerator Radiation Induced Activation, October 13-17, 2008. Paul Scherrer Institut, Switzerland. P. 48-54.

в научных монографиях:

27. Тихоненко А.В. Интегрирование уравнений движения заряженных частиц в MAPLE. – Обнинск: ИАТЭ, 2007. – 204 с.

28. Тихоненко А.В. Тензорное исчисление и его приложения в прикладных математических пакетах. – Обнинск: ИАТЭ, 2007. – 204 с.

в учебных пособиях:

29. Тихоненко А.В. Компьютерные математические пакеты в курсе общей физики. Обнинск: ИАТЭ, 2003. – 84 с.

30. Тихоненко А.В. Компьютерный практикум по общей физике. Часть 1. Классическая Механика. Обнинск: ИАТЭ, 2003. – 84 с.

31. Тихоненко А.В. Компьютерный практикум по общей физике. Часть 2. Механические колебания и волны. Термодинамика и молекулярная физика.

Обнинск: ИАТЭ, 2004. – 84 с.

32. Тихоненко А.В. Компьютерный практикум по общей физике. Часть 3. Электричество и магнетизм. Обнинск: ИАТЭ, 2004. – 84 с.

33. Тихоненко А.В. Компьютерный практикум по общей физике. Часть 4. Оптика.

Обнинск: ИАТЭ, 2004. – 80 с.

34. Тихоненко А.В. Компьютерный практикум по общей физике. Часть 5. Квантовая физика. Обнинск: ИАТЭ, 2004. – 80 с.

35. Тихоненко А.В. Компьютерные математические пакеты в курсе «Линейные и нелинейные уравнения физики». Часть 1. Гиперболические уравнения в MAPLE. Обнинск: ИАТЭ, 2005. – 80 с.

36. Тихоненко А.В. Компьютерные математические пакеты в курсе «Линейные и нелинейные уравнения физики». Часть 2. Параболические уравнения в MAPLE. Обнинск: ИАТЭ, 2005. – 80 с.

37. Тихоненко А.В. Решение краевых задач для двумерного уравнения Лапласа методом разделения переменных в MAPLE. Обнинск: ИАТЭ, 2005. – 80 с.

38. Тихоненко А.В. Компьютерные математические пакеты в курсе «Электродинамика». Часть 1. Электромагнитные поля в проводящих средах и поля движущегося заряда. Обнинск: ИАТЭ, 2005. – 80 с.

39. Тихоненко А.В. Решение уравнения Шредингера для одномерного рассеяния в MAPLE и MATHEMATICA. Обнинск: ИАТЭ, 2005. – 80 с.

40. Тихоненко А.В. Векторный анализ в прикладных математических пакетах.

Обнинск: ИАТЭ, 2006. – 80 с.

41. Тихоненко А.В. Решение уравнений колебаний в прикладных математических пакетах. Обнинск: ИАТЭ, 2008. – 64 с.

42. Тихоненко А.В. Компьютерные аналитические методы решения задач электростатики и магнитостатики. – Обнинск: ИАТЭ, 2008. – 48 с.

на официальных сайтах производителей математических пакетов:

MATHSOFT ENGINEERING & EDUCATION INC 43. Tikhonenko A.V. Visual Quantum Mechanics of the Isotropic Oscillator and Hydrogen-Like Atoms with MATHCAD. Electronic Book in MATHCAD. Mathsoft Engineering & Education Inc. 2002.

http://www.ptc.com/appserver/wcms/standards/fileothumbredirect.jsp?&im_dbkey= 62175&icg_dbkey=888.

44. Tikhonenko A.V. Solar System Mechanics. Mathsoft Engineering & Education Inc. 2003.

http://www.MATHCAD.com/resources/MATHCAD_files/DisplayCategory.asp?c=48.

45. Tikhonenko A.V. Interference Effects of Monochromatic and Non-monochromatic light. Mathsoft Engineering & Education Inc. 2003.

http://www.MATHCAD.com/resources/MATHCAD_files/DisplayCategory.asp?c=48.

46. Tikhonenko A.V. One-Dimensional Motion of Quantum Particles: Tunneling Through a Square Well. Mathsoft Engineering & Education Inc. 2004.

http://www.MATHCAD.com/resources/MATHCAD_files/DisplayCategory.asp?c=48.

47. Tikhonenko A.V. Wave Functions of the Hydrogen-like Atoms. Mathsoft Engineering & Education Inc. 2004.

http://www.MATHCAD.com/resources/MATHCAD_files/DisplayCategory.asp?c=48.

WOLFRAM RESEARCH INC 48. Tikhonenko A.V. Quantum Mechanics of the Bound Isotropic Systems with MATHEMATICA: Analysis and Visualization of Quantum States. Electronic Book.

Wolfram Research Inc. 2003.

http://library.wolfram.com/infocenter/MathSource/4285/ WATERLOO MAPLE INC 49. Tikhonenko A.V. Vibrations: Complete Set of Lessons. Electronic Book in MAPLE. Waterloo MAPLE Inc. 2003.

http://www.maplesoft.com/applications/view.aspx?SID=4750. Tikhonenko A.V. Equations of mathematical Physics in MAPLE. Waterloo MAPLE Inc. 2005.

http://www.MAPLEsoft.com/applications/app_center_view.aspx?AID=1867.

51. Tikhonenko A.V. Two-dimensional partial elliptic differential equations in MAPLE.

Waterloo MAPLE Inc. 2007.

http://www.maplesoft.com/applications/app_center_view.aspx?AID=2052. Tikhonenko A.V. Tensor algebra in the four-dimensional pseudo-euclidean spacetime. Waterloo MAPLE Inc. 2007.

http://www.maplesoft.com/applications/app_center_view.aspx?AID=2053. Tikhonenko A.V. Tensors of relativistic electrodynamics with MAPLE. Waterloo MAPLE Inc. 2007.

http://www.maplesoft.com/applications/app_center_view.aspx?AID=2054. Tikhonenko A.V. Quantum scattering by the one-dimensional potential barrier in MAPLE. Waterloo MAPLE Inc. 20http://www.maplesoft.com/applications/app_center_view.aspx?AID=2055. Tikhonenko A.V. Quantum scattering by the one-dimensional potential step in MAPLE. Waterloo MAPLE Inc. 2007.

http://www.maplesoft.com/applications/app_center_view.aspx?AID=20







© 2011 www.dissers.ru - «Бесплатная электронная библиотека»

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