автореферат диссертации по информатике, вычислительной технике и управлению, 05.13.16, диссертация на тему:Анализ точности моделирования параметров электроядерных установок

^ гф

о

/

^ЪЕДИНЕННЫЙ ИНСТИТУТ ЯДЕРНЫХ ИССЛЕДОВАНИЙ

па правах рукописи УДК 621.039

БАТЯЕВ Вячеслав Феликсович

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

Специальность: 05.13.16 - применение вычислительной техники, математического моделирования и математических методов в научных

исследованиях

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата физико-математических наук

ДУБНА - 1999

Работа выполнена в ГНЦ РФ Институте теоретической и экспериментальной Физ ки, г. Москва

Научный руководитель:

кандидат физико-математических наук Титарепко Ю.Е.

Официальные оппоненты: доктор физико-математических наук,

профессор Барашенков B.C.

доктор физико-математических наук, Севастьянов JI.A.

Ведущая организация:

ГНЦ РФ Физико-энергетический институт, г. Обнинск.

Защита диссертации состоится " " 2000г. в ^^.часов на зас

дании диссертационного совета Д047.01.04 при Лаборатории вычислительной технш и автоматизации ОИЯИ по адресу: г.Дубна Московской области, ЛВТА ОИЯИ.

С диссертацией можно ознакомится в библиотеке ОИЯИ.

Автореферат разослан "'"iiZl" ______г* 1999г.

Ученый секретарь Совета, кандидат физико-математических наук З.М.Иванченко

\Ш-01Ис!и. О

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

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

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

В рамках поставленной цели были решены следующие задачи:

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

2) определен набор параметров, наиболее пригодных для верификации расчетных зограмм;

3) разработана методика сравнения экспериментальных результатов по выходу 1ер-продуктов с результатами, полученными с помощью различных используемых в ютоящее время программ;

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

5) создан взаимосвязанный самосогласованный комплекс программ для расчета [.ернофизических параметров толстых мишеней;

6) выполнены расчеты параметров мишенных устройств ЭлЯУ, проектируемых ИТЭФ с целью их оптимизации.

Научная повизна

1. Впервые изучена точность используемых программ для моделирования пара\ тров мишенных устройств из тяжелых материалов.

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

3. Проведена верификация основных используемых в насточщее время nporpaj путем сравнения расчетных результатов с экспериментальными по выходам яде продуктов из материалов, предполагаемых для использования в качестве мишенш и конструкционных материалов ЭлЯУ.

4. Рассчитаны основные ядерно-физические параметры проектируемых в ИТЕ мишенных устройств ЭлЯУ.

Практическая ценность работы

Практическая ценность работы состоит в исследованной точности расчетов ост пых ядерно-физических параметров мишенных устройств ЭлЯУ. Особенность дан» работы, отличающей ее от других подобных работ (напр. состоит как в использое нии широкого набора программ, применяемых для расчета ЭлЯУ, так и в использован! набора экспериментальных данных, требуемых сугубо для проектирования ЭлЯУ.

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

Точность продемонстрирована на примере широкого пакета программ, включающ] широкий спектр моделей внутриядерного взаимодействия и активно использующих для расчетов параметров ЭлЯУ, как в РФ, так и зарубежом.

Точнось программ продемонстрирована по: (1) независимым и кумулятивным в: ходам ядер-продуктов в тонких мишенях из мишенных и конструкционных материал! ЭлЯУ; и (2) скоростям реакций на поверхности толстой W-мишени, облученной прот нами с энергией 1 ГэВ.

1R. Michel, Р. Nagel, Intereiiational Codes and Model Intercomparison for Intermediate Energy Activati Yields, NSC/DOC(97)-l, NEA/P&T No 14.

1. Тонкие мишени. Для корректного сравнения расчетных результатов с экспе-шментальными создана методика сличения экспериментальных и расчетных резуль-атов по выходам независимых и кумулятивных ядер-продуктов в тонких мишенях, 1азирующаяся на методике экспериментального получения непосредственно в ИТЭФ :умулятпвных и независимых выходов ядер продуктов.

Используя более 720 экспериментальных значений выходов ядер-продуктов на [драх-мишенях конструкционных материалов ЭлЯУ, измеренных в ИТЭФ (облучались 09Bi,™'Pb, 208РЬ, 207РЬ, 206РЬ, 63Си, 65Си, и 59Со протонами низких (0.07 - 0.2 ГэВ) [ высоких (1.2 и 1.5ГэВ) энергиях); было выполнено сравнение с ними расчетных ре-ультатов, полученных по программам: СЕМ95, CASCADE, INUCL, НЕТС, LAHET, 3NASH, ALICE, QMD, NUCLEUS, YIELDX, и полуэмпиричсской формуле Foshina et .1.

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

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

2. Толстая W-мишень. Используя программы LAHET, СЕМ95, INUCL и CASCADE для моделирования скоростей реакций на поверхности толстой W-мишени 1Ыло проведено сравнение расчетных значений скоростей пороговых реакций актива-(ионных детекторов (In, Au, А1) с экспериментальными значениями, полученными в 1ТЭФ.

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

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

Апробация работы

Основные положения работы докладывались на: International Workshop on Nuclea Methods of Transmutation of Nuclear Wastes: Problems, Perspectives, Cooperativ Research, Dubna, Russia, May 29-31,1996; Second Intternational Conference on Acceleratoi Driven Transmutation Technologies and Applications, Kalmar, Sweden, June 3-7, 199G; X Международном семинаре по точным измерениям в ядерной спектрометрии, 2-6 сени бря 1996, Саров; International Conference on Nuclear Data for Science and Technology, Ma 19-24, 1997, Trieste, Italy; 3rd Specialists' Meeting on High Energy Nuclear Data, Marc 30-31, 1998 at JAERI, Tokai, Ibaraki-ken, Japan; Second International Topical Meetin on Nuclear Applications of Accelerator Technology (AccApp;98), Gatlinburg, TN, US/ September 20-23, 1998; 3rd International Conference on Accelerator-Driven Transmutatio Technologies and Applications, Praha, June 7-11,1999; Third International Topical Meetin on Nuclear Applications of Accelerator Technology (AccApp'99), Long Beach, CA, US/ November 14-18, 1999.

Основные результаты работы опубликованы в сборниках указанных конференци£ в виде препринтов ИТЭФ и LANL, а также в журналах ВАНТ Сер. "Физика Ядерны Реакторов" и "Nucl. Instr. and Meth.".

Структура и объем работы.

Диссертация состоит из введения, трех глав, заключения, списка литературы. Он изложена на 147 страницах, включая 43 рисунка, 44 таблицы и список литературы и 169 наименований.

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

Во Введении описано назначение электроядерных установок и дан краткий обзор [роектируемых мишенных устройств для реализации ЭлЯУ, который показывает наиболее часто предлагаемые мишенные материалы, энергии и мощности протонных пуч-;ов. Представлены необходимые для проектирования параметры мишенных устройств, I также основные методы и программы, используемые для их моделирования. Поста-¡лена задача по необходимости исследования точности расчетов активно используемх I настоящее время программ.

В Главе 1 проведено исследование экспериментальной и расчетной точности основ-[ых параметров мишеней. Показана удовлетворительная (~15%), по большинству нро-рамм, точность расчета полного нейтронного выхода и полного энерговыделения. Ис-юльзуя экспериментальные данные показано, что данная точность соответсвует по юрядку величины экспериментальной точности. Показана менее удовлетворительная ~30%) точность расчетов локального энерговыделения, в первую очередь, за счет проектов ядерных взаимодействий. Показана неудовлетворительная точность (фактор 2 и ыше) расчета выходов ядер-продуктов. Это позволило сделать выбор в пользу выхода дер-продуктов как параметра для верификации программ. Показана необходимость спользования изотопных мишеней для накопления ядерных данных для верификации рограмм. Проведен статистический анализ адронного каскада в толстой мишени, в ре-ультате показаны диапазоны энергий нуклонов, принимающих наибольшее участие в [ежъядерном каскаде, и, следовательно, наиболее важные для верификации программ.

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

В связи с тем, что большая часть экспериментальных даиых содержит кумулятив-ые выходы, была разработана методика сравнения результатов различных программ экспериментальными значениями по независимым и кумулятивным выходам ядер-родуктов. Методика согласована с разработанной в ИТЭФ экспериментальной мето-икой определения выходов ядер-продуктов.

Используя эту методику было проведено сравнение экспериментальных и расчетных анных по выходам ядер-продуктов.

В сравнении использовались:

1. Экспериментальные результаты по выходам ядер-продуктов в: 209Bi, 208Pb, 207Pb, 206РЬ, облученных протонами с энергией 0.13 и 1.5 ГэВ; ""'РЬ, облученном с энергией 1.5 ГэВ; 59Со, облученном протонами с энергией 0.07, 0.10, 0.13 и 1.2 ГэВ; 63Cu, 65Си, облученных протонами с энергией 0.13, 1.2 и 1.5 ГэВ.

2. Результаты моделирования экспериментально измеренных ядер-продуктов используя программы:

• каскадно-экситонную СЕМ95,

• каскадно-испарительно-делительную CASCADE,

• каскадно-предравновсспо-испаритсльпо-делительную INUCL,

• каскадно-испарительную НЕТС,

• каскадно-испарительно-делительную LAHET,

• предравновесно испарительную GNASH,

• предравновесно-испарительную ALICE,

• квантово-молекулярно-динамическую QMD,

. • каскадно-испарительно-делительную NUCLEUS,

• полуфеномеиологическую YIELDX,

• полуфеноменологическую формулу Foshina.

Результаты сравнения представлены в виде таблиц и рисунков понуклидного срав нения, таблиц и рисунков статистических отклонений "расчет/эксперимент" и в вид| рисунков массовых выходов. Сделаны выводы о сходимости расчетных и эксперимент тальных результатов по каждой программе, группе ядер и энергии. Рассчитаны сред неквадратичные отклонения расчетных результатов от экспериментальных по четыре», основным наборам ядра-мишени/энергия и сделаны выводы только о качественном со ответствии большинства значений выходов моделируемых по всем программам ядер продуктов экспериментальным результатам. Показано, например, что в случае тяже лых ядер (свинец, внесут) облучаемых при высоких энергиях, наилучшей предсказа тельной силой обладает программа СЕМ95, предсказывающая выходы ядер продукте! со среднеквадратичным отклонением 1.92, т.е. с точностью на уровне всего лишь

В Главе 3 представлен взаимосвязанный самосогласованный комплекс програмх для расчета ядерно-физических параметров толстых мишеней, используя программь

1ЕТС, МСИР, БРЕСТ, ЕСтЯ. Т\УСЮА1ЧТ, баз данных ENDF/B и БНАБ. Произведе-[а модификация программы БРЕСТ для расчета требуемых параметров мишенных стройств, включая нейтронный выход, потоки каскадных частиц, полное и локальное нерговыделение. Разработана и реализована методика расчета парциального энерго-ъщеления в компонентах смеси, заданной гомогенно.

Используя эти программы проведены оптимизационные расчеты параметров мишен-[ых устройств, нредставленых для проектирования ЭлЯУ. Были проанализированы две кпдкосвинцовьге, свшгоцово-висмутовая, вольфрамовая и свинцово-тяжеловодная ми-пени, предложенные для проектной реализации в ходе ряда проектов по исследованию ЭлЯУ. В ходе оптимизаций получены:

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

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

- зависимость выхода нейтронов свинцово-тяжеловодной композиции в зависимости >т пористости;

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

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

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

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

1. Создан взаимосвязанный комплекс программ для расчета нейтронофизически: параметров толстых мишеней, используя программы НЕТС, МСКР, БРЕСТ, НС£ TWODANT, баз данных ЕКОР/В и БНАБ (рис. 1). Проведена модификация программ* БРЕСТ для расчета требуемых параметров мишенных устройств, включая нейтропны! выход, потоки каскадных частиц и полное и локальное энерговыделение (рис. 2). Разра ботана и реализована методика расчета парциального эперговыделеиия в компонента; смеси, заданной гомогенно.

Выполнена модеринизация программы НЕТС путем включения в нее альтернатив ных модулей из программ СЕМ95 и ШиСЬ, моделирующих адрон-ядерные взаимодей ствия.

2. Показана удовлетворительная (~15%), по большинству программ, точность рас чета полного нейтронного выхода и полного энерговыделения. Используя эксперимен тальные данные показано, что данная точность соответствует по порядку величин! экспериментальной точности. Показана менее удовлетворительная (~30%) точност! расчетов локального энерговыделения, в первую очередь, за счет продутков ядерньс взаимодействий, а также неудовлетворительная точность (фактор порядка 2 и выше расчета выходов ядер-продуктов.

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

4. Проведен статистический анализ адронного каскада в толстой мишени, в ре зультате показаны диапазоны энергий нуклонов, принимающих наибольшее участш в межъядерном каскаде, и, следовательно, наиболее приоритетные для верификацщ программ. Получено, что наиболее существенный вклад в динамику адронного каска да играют диапазоны энергий: Е0-0.25 ГэВ - ЕО (ЕО - исходная энергия налетающей протона) и 0.1 - 0.25 ГэВ (рис. 3).

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

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

Зис. 1: Схема применения программ для вычисления параметров мишенных устройств ЭлЯУ.

7. Проведено сравнепие экспериментальных я расчетных данных по 11 программам [о выходам ядер-продуктов.

Результаты сравнения представлены в виде таблиц и рисунков понуклидного срав-[ения (рис. 4), таблиц и рисунков статистических отклонений "расчет/эксперимент" табл. 1 и рис. 5) и в виде рисунков массовых выходов (рис. 6). Сделаны выводы о ходимости расчетных и экспериментальных результатов по каждой программе, груп-ie ядер и энергии. Рассчитаны среднеквадратичные отклонения расчетных результатов it экспериментальных по четырем основным наборам ядра-мишени/энергия (табл. 2 и >ис. 7, 8) и сделаны выводы только о качественном соответствии большинства моде-шруемых по программам выходов ядер-продуктов экспериментальным результатам. В [астностп, показано, что в случае тяжелых ядер (свинец, висмут), облучаемых протона-га при высоких энергиях, наилучшая сходимость, представленная программой СЕМ95, гаходится чуть лучше фактора 2.0 Только в отдельно взятых случаях (про-

рамма CASCADE при 208РЬ-0.13ГэВ и СЕМ95 при 63Си-0.13ГэВ) программы показали федсказательную силу на уровне требуемой точности для проектных работ (отклоне-шя менее 30%).

Рис. 2:11-2 распределение удельного энерговыделения (в кэВ/(р*см2)) мишени на осно] исевдосжиженной РЬ-020 композиции при пористости 0.7.

Взаимодействия по налетающей Взаимодействия по налетающей энергии энергии вторичных частиц

250 500 750 1000 Энергия взаимодействия (МэВ)

Нейтронный выход по налетающей энергии

200 400

Энергия взаимодействия (МэВ)

Нейтронный выход по налетающей энергии вторичных частиц

250 500 750 1000 Энергия взаимодействия (МэВ)

200 400

Энергия взаимодействия (МэВ)

Рис. 3: Статистика адрон-ядерных взаимодействий и генераций низкоэнергичных (£„<20МэВ) нейтронов в свишюво-висмутовой мишени размером 02О*6Осм, в зависимости от энергии взаимодействия внутри мищени. Исходная энергия протонов - 1 ГэВ. Расчет с помощью программы НЕТС с встроенным модулем ШиСЬ. Число исходных протонов - 10000.

выходы продуктов (мбарн)

Ю' 10 1

-1

10 10

10' 10 1

СЕМ95

СЕМ95

вац^вавввйдв^Р!

Л 0) О Ф <Л ГО «¿¿ГГ 3 о>зг о>3 3 г л г £ ОД^Р-ОСт ОГ7 V т т т т О О

' ^ л! 9

ю ю 1

•1

ю

-2

10

ю2 10 1

-1

10

•2

10

иНЕТ

ЬАНЕТ

<0 Ф Ф <1>

к а: сс сп

<р о> -г- см

Г4« г- со со

1Л СО .Ь ОТ ^

оыо«

см со аз ц") о и со да »■

^ т- т- ^ г- су су ГЧ «00

Г\1

? т т со ю ,9 Р

3 п .

Рис. 4: Сравнение экспериментальных (закрашенные символы) и расчетных значений (1 мбарн) по программам СЕМ95 и ЬАНЕТ (открытые сиволы) выходов ядер-продукто! в 209В1, облученного протонами с энергией 1.5ГэВ.

НЕТС

1 1 1 . | 1 | 1 | 1-

■А 1

У1Е1-ОХ

г1—I—1—г

_Ъва]_I_|_

0.01 0.1 1 10 100 РобЖпэ е( а1.

0.01 0.1

аю

СТса1с,/°(

ехр,1

10 100 ®са1с,/®ехр,1

Рис. 5: Статистика отношений расчетных и экспериментальных значений по выходам продуктов в 2МВ1, облученным протонами с энергией 1.5 ГэВ

Выход (мбари)

Т—I—I—I—|—I—I—I—I—]—I—I—I—I—р-

—I—|—1—I—1—I—I—I—I—I—I—I—I—I—I—1—I—I—I—I—I—I—I—г

10'

10

10

10

А

Л

1 А ^¿Г

ш

о о о %

^ 1 _д*

о»

о

/

°о ♦

0« о о

о о

о

<0 о

л о

Ехрептеп! СЕМ95 1.АНЕТ

о мио.

□ НЕТС А У1ЕШХ О САЭСАОЕ

о *

° ®й о

<Ьо

И) № дд

о

ао д о д

ф® о д

J_1_1_I_I_I_I_I_I_I_1_1_1_1_

_1_1_1_1_I_I_I_1_1_I_I_I_I_и

25 50 75 100 125 150 175 200

Масса ядср-продукто!

А

д

+

+

■

о

О

0 <00

Рис. 6: Массовое распедедеиие продуктов реакции 209В1(р,х) при облучении протонам! с энергией 1.5 ГэВ.

Таблица 1: Статистика сравнений чксперимента пьных и расчетных результатов выходов ядер-нродукгов в 209В1, облученном протонами с энергией 0.13 и 1.5 ГэВ; А/т -число экспериментальных выходов; N0 - число выходов, отобранных для сравнения; дг5 - число выходов, промоделированных программой; Л^.,, - число промоделированных выходов, отличающихся от экспериментальных значений не более чем в 1.3 и в 2.0 раза, соответственно: (Я) - среднеквадратичный фактор отклонения расчетных значений выходов ядер-продуктов от экспериментальных; 8({Я)) - дисперсия фактора отклонения.

Программа

А'т=26.

.13 ГэВ

N0=24

Ер = 1.5 ГэВ ЛГГ=99, Дз=88

М^^/Г^/^^^т) N^/N^/N5 (Н) $({Н))

СЕМ95 14/22/24 1.46 1.30 25/48/65 2.38 2.00

шись 3/6/21 13. 7.99 29/48/85 3.46 2.72

3/6/19* 3.22 1.95

САЭСАИЕ 12/22/24 1.63 1.48 27/48/66 2.68 2.31

НЕТС 8/19/24 2.38 2.05 23/42/61 2.92 2.39

ЬАНЕТ(18АВЕЬ) 16/23/24 1.39 1.26

' ЬАНЕТ(ВеП1ш) 16/22/24 1.41 1.27 23/63/87 2.22 1.78

АИСЕ 13/21/24 1.76 1.59

С КАБЫ 2/6/9 2.35 1.78

QMD 8/16/24 2.11 1.67 20/63/87 2.36 1.83

ПЕЫЭХ 13/21/24 1.52 1.34 18/41/83 3.52 2.47

с исключением 200В1, 201В1, 202В1, г02Ро, и 203Ро

Таблица 2: Общая статистика сравпений экспериментальных и расчетных результатов по выходам ядер-продуктов. Обозначения Л^, Л^, Л^о, (Я), 5((Я)) даны в

оглавлении табл. 1.

63'65Си, Бр=0.13ГэВ 59Со,63,65Си, Ер=1.2, 1.5ГэВ

59Со, Ер=0.07, 0.10, 0.13, 0.20ГэВ Программа ЛГТ=9 2, Nc=71 ЛГТ=193, Л'с=161

XC,JNC2JNs (H) 5«Я» NC,JNC2JNs <H> да»

СЕМ95 27/46/64 2.24 1.88 46/96/149 2.22 1.69

LAHET 27/45/69 2.33 1.86 45/112/151 1.90 1.47

INUCL 15/39/71 2.51 1.74 38/89/153 2.20 1.65

НЕТС 8/21/49 5.47 3.47 19/43/96 4.47 2.93

CASCADE 5/10/33" 3.97 2.23 20/58/152 3.40 2.09

ALICE(Fermi) 7/15/34" 3.17 2.02 - - -

ALICE(Kataria) 12/21/34" 2.15 1.64 - - -

YIELDX 24/45/71 2.28 1.78 64/142/160 1.59 1.33

NUCLEUS 8/14/30" 2.85 2.02 16/54/89* 2.18 1.51

Foshina et al. 14/36/71 2.72 1.90 43/108/156 2.25 1.81

QMD 11/24/32" 1.98 1.60 17/50/78* 2.45 1.81

206,207,208 pj^ 209Bii =0.13ГэВ I 200,207^08,natpjj 2"9gj JJ =1.5ГэВ

Программа ЛГг=90, =76 JVr=471, 7VC=418

NCJNc,JNs (Я> 5({Я» Nc,JNc>JNs <Я> 5({Я»

CEM95 37/59/76 2.16 1.94 154/257/327 1.92 1.66

LAHET 42/70/76 1.47 1.28 103/312/417 2.02 1.60

INUCL 23/48/71 2.14 1.68 123/247/414 3.02 2.26

HETC 23/55/76 2.21 1.78 92/215/308 3.29 2.80

CASCADE 34/68/76 1.61 1.40 123/246/330 2.20 1.85

ALICE(Fermi) 39/69/76 1.62 1.44 - - -

YIELDX 39/68/76 1.49 1.29 77/199/413 3.23 2.20

QMD 13/40/76 2.32 1.58 88/284/415 2.33 1.80

GNASH 2/6/9»" 2.35 1.78 - - -

* здесь ЛГг = 110 и Иа = 93; " здесь Л^ = 41 и Л^ = 34. *" было промоделировано только 9 ядер-продуктов

10 9 8 7 6 5 4

3

1.3(30%) 1

I 1 I 1 1 1 1 1 1 1

: е,= 0.07, 0 10, 0.13, 0.20 ГэВ, Число сравнений = 71

- 5.47

■ . '.:'.. ' 3.97

Г98 2.24 • „«-V 2.72 2-85 • - . -•-■ --'=•*-'■',•' 2.28 2.33 .

> i i .......

Q

о

LU

O

_J

<

X Q

LU

Ш X

3

o

Э

со

LU —I

O 3

Ш

o

o <

o

со <

o

Ш I

20

3 2

1.3(30%) 1

Ep = 1.2, 1.5 ГэВ, Число сравнений ■

161

4Í47

3.40

1.90

2.18 2.20 2.22 2.25

2,45

í.59

X I- со ш га

Q ш э о О с

-J X ш 3

Ш >- 5 _1 о э Z Z ш о ел о LL

О О

ш

Q <

О

со <

о

ш I

4

L

Рис. 7: Общий фактор отклонения по всем отношениям расчетных и экспериментальных шачений для "легких" (конструкционных) ядер-мишеней (5ЭСо, 63,60 Си)

1 Г

1-г

Ер = 0.13 ГэВ, Число сравнений. = 76 \ : ......

2.14 2.16 2-21

2.32 2-35

1.47 1-49

1.61 1.62

1.3(30%)-

10 9 8 7 6 5 4

1.3(30%) 1

ш х

3

х о

ш

J_1_

ш

о <

о

сл <

о

ш

о

о э 2

5 ш О

J_|_

ш X

о О

х

ю <

2 а

Ер = 1.5 ГэВ, Число сравнений — 418

3.02

3.23 ': 3.29

1.92

2.02'

2 20

2 33

ш 1- ш а

о ш а о

5 X < п э

ш < о у г

О _1 С/5

<

О

х о

о 1-

3

Рис. 8: Общий фактор отклонения по всем отношениям расчетных и экспериментальны: значений для "тяжелых" (мишенных) ядер-мишеней (209В1, 208,207,20б,тш1рЬ)

8. Проведено сравнение расчетных результатов по четырем программам (LAHET, СЕМ95, INUCL, CASCADE) с экспериментальными данными по скоростям реакций на поверхности толстой W-мишени. Полученные значения отклонений расчетных результатов от экспериментальных сходны с отклонениями в случае тонких мишеней (факторы 1.6-г2.4). Сделаны предварительные выводы о корректности расчетных спектров и используемых сечений.

9. Проведены расчеты параметров мишенных устройств, предлагаемых для проектирования ЭлЯУ с целью их оптимизации. Были проанализированы две жидкосвинцо-вые, свипоцово-висмутовая, вольфрамовая и свинцово-тяжеловодная мишени. В ходе оптимизаций получено:

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

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

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

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

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

Таблица 3: Программы, результаты которых отличаются от экспериментальных со средним отклонением менее 2.0.

Мишени Энергия протонов (ГэВ)

<0.20 1.2, 1.5

"Мишенные" ЬАНЕТ (1.47)

материалы У1ЕЬВХ (1.49) СЕМ95 (1.92)

206,207,208,паф^ САЭСАБЕ (1.61)

209 Ш АЫСЕ (1.62)

"Конструкционные" У1ЕЬОХ (1.59)

материалы С^МБ (1.98) ЬАНЕТ (1.90)

мСо, 63'65Си

35-1 30 25 20-| ^ 151050

—■— Нейтронорождаемость —•— Полный вылет нейтронов —— Вылет через боковую поверхность —▼— Поглощение нейтронов

▼— 1

т--_

1-1-1-'-1-'-1---1-1-1—

0.0 0.2 0.4 0.6 0.8 1.0

пористость

▼

Рис. 9: Нейтронные параметры мишени на основе пссвдосжилесиной РЬ-БгО композиции в зависимости от пористости.

ОСНОВНОЕ СОДЕРЖАНИЕ ДИССЕРТАЦИИ ИЗЛОЖЕНО В РАБОТАХ

L. V.F. Batyacv, "Nuclear Data requirements in Intermediate Energy Range for Improvement of Calculations of ADTT Target Process", "Proc. of Int. Workshop on Nucl. Methods of Transmutation of Nuclear Wastes: Problems, Perspectives, Cooperative Research," Dubna, Russia, May 29-31, 1996, World Scientific, M.Khankhasayev, H.Plendl, Zh.B.Kurmanov (1997), pp. 166-172; В.Ф.Батяев, "Потребности в ядерных данных промежуточных энергий для верификации расчетов мишенных процессов электроядерных установок", Материалы XI Международного семинара по точным измерениям в ядерной спектрометрии, 2-6 сентября

1996, Саров, Вопросы атомной науки и техники, Сер. Физика Ядерных реакторов,

1997, ТИЯС-XI, Спец. Выпуск, 179-183.

!. Yu. Е. Titarenko, О. V. Shvedov, V. F. Batyaev, Е. I. Karpikhin, V. М. Zhivun, А. В. Koldobsky, М. М. Igumnov, I. S. Sklokin, R. D. Mulambetov, A.N. Sosnin, H. Yasuda, II. Takada, S. Chiba, Y. Kasugai, S. G. Mashnik, R. E. Prael, M. B. Chadwick, T. A. Gabriel, M. Blann, "Experimental and Computer Simulation Study of Radionuclide Yields in the ADT Materials Irradiated with Intermediate Energy Protons", "Proceedings of the third Specialists' Meeting on High Energy Nuclear Data", March 30-31, 1998 at JAERI, Tokai, Ibaraki-ken, Japan, JAERI-Conf 98-016 pp.125135; Yu. E. Titarenko, О. V. Shvedov, V. F. Batyaev, E. I. Karpikhin, V. M. Zhivun,

A. B. Koldobsky, M. M. Igumnov, I. S. Sklokin, R. D. Mulambetov, A.N. Sosnin, H. Yasuda, H. Та Takada, S. Chiba, Y. Kasugai, S. G. Mashnik, R. E. Prael, M.

B. Chadwick, T. A. Gabriel, M. Blann, "Experimental and Computer Simulation Study of Radionuclide Formation in the ADT Materials Irradiated with Intermediate Energy Protons", "Second Int. Topical Meeting on Nuclear Applications of Accelerator Technology (AccApp'98)", Gatlinburg, TN, Sept. 20-23, 1998, pp. 164-171; Yu. E. Titarenko, О. V. Shvedov, V. F. Batyaev, E. I. Karpikhin, V. M. Zhivun, R. D. Mulambetov, A.N. Sosnin, S. G. Mashnik, R. E. Prael, T. A. Gabriel, M. Blann, "Experimental and Computer Simulation Study of Radionuclide Production in Heavy Materials Irradiated by Intermediate Energy Protons", "Topical Meeting - Nuclear Applications of Accelerator Technology (AccApp'99)", Long Beach, CA, November 1418, 1999, pp. 212-221.

I. Yu. E. Titarenko, О. V. Shvedov, V. F. Batyaev, E. I. Karpikhin, V. M. Zhivun, A.B. Koldobsky, M. M. Igumnov, A.N. Sosnin, S. G. Mashnik, R. E. Prael, M. B.

Chadwick, T. A. Gabriel, M. Blann, "Exerimental and Theoretical Study of the Yield: of Radionuclides Produced in 209Bi thin target Irradiated by 1500 MeV and 130MeV Protons" Nucllnstr. and Meth. A 414 (1998) 73-99; Los Alamos Preprint LA-UR-97-3787; nucl-th/9709056.

4. P.P.Blagovolin, V.D.Kazaritsky, V.F.Batyaev, E.A.Zolotareva, E.S.Nikolaevsky V.T.Gorshkov, V.A.Shulyndin, "Accelerator-Driven Molten-Fluoride Reactor witl Modular Heat Exchangers on Pb-Bi Eutectic", "Proc. of Int. Workshop on Nucl Methods of Transmutation of Nuclear Wastes: Problems, Perspectives, Cooperativt Research," Dubna, Russia, May 29-31, 1996, World Scientific, M.Khankhasayev H.Plendl, Zh.B.Kurmanov (1997), pp. 235-240.

5. V.F. Batyaev, "Neutronics and Power Deposition Parameters of the Targets Proposec in the ISTC Project 17", "Proc. of Int. Workshop on Nucl. Methods of Transmutatioi of Nuclear Wastes: Problems, Perspectives, Cooperative Research" Dubna, Russia May 29-31,1996, World Scientific, M.Khankhasayev, H.Plendl, Zh.B.Kurmanov (1997) pp. 216-221.

6. V.D.Kazaritsky, M.L.Okhlopkov, V.R.Mladov, A.N.Legotsky, S.S.Rozhkov B.V.Sharov, M.M.Kats, V.V.Seliverstov, V.F.Batyaev, "A concept for a higt flux spallation neutron source as a radiation damage test facility", Nucl.Instr. ant Meth. A 414 (1998) 43-48.

7. V.D.Kazaritsky, P.P. Blagovolin, M.L.Okhlopkov, V.R.Mladov, V.F. Batyaev, V.V Seliverstov, E.B. Strakhov, V.l. Volk, A.Y. Vakhrushin, M.I. Zavadsky, Z.G. Ilina S.V. Petrov, S.V. Petrin, E.F. Fomushkin, "A conceptual design study and efficiency assessment of an accelerator driven reactor system for the transmutation of nuclear waste", Nucl.Instr. and Meth. A 414 (1998) 21-27.

ВВЕДЕНИЕ

1 Назначение ЭлектроЯдерных Установок (ЭлЯУ).

2 Основные разработки мишенных частей ЭлЯУ.

3 Обзор параметров мишенных устройств, необходимых для проектирования ЭлЯУ.

4 Обзор методов моделирования ядерных процессов мишенных частях ЭлЯУ.

4.1 Основные модели нуклон-ядерных взаимодействий в промежуточной области энергий (выше 20МэВ).

4.2 Программные комплексы для моделирования адрон-ядерного каскада в протяженной геометрии и оценки основных параметров мишенной части ЭлЯУ.

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

6 Постановка задачи

7 Актуальность и новизна

1 ПОТРЕБНОСТИ В ВЕРИФИКАЦИИ РАСЧЕТНЫХ ПРОГРАММ

1.1 Оценка точности расчета параметров мишенных устройств ЭлЯУ.

1.1.1 Нейтронный выход

1.1.2 Параметры энерговыделения.

1.1.3 Выход вторичных ядер-продуктов.

1.2 Выбор параметра для тестирования расчетных программ.

1.3 Статистика адронных взаимодействий в толстых мишенях.

1.4 Выводы.

2 ВЕРИФИКАЦИЯ РАСЧЕТНЫХ ПРОГРАММ.

2.1 Особенности верификации прогррамм по выходам ядер-продуктов.

2.2 Обзор имеющихся экспериментальных данных.

2.3 Используемые в верификации расчетные программы.

2.4 Систематизация сравнений расчетных и экспериментальных данных по выходам ядер-продуктов.

2.4.1 Сравнение экспериментальных и расчетных результатов.

2.4.1.1 209Bi, облученный протонами с энергией 0.13 и 1.5 ГэВ.

2.4.1.2 20б,207,208р^ облученные протонами с энергией 0.13 ГэВ.

2.4.1.3 206,207,208,natpj^ облученные протонами с энергией 1.5 ГэВ.

2.4.1.4 63'65Си, облученные протонами с энергией 0.13, 1.2 и 1.5 ГэВ и 59Со, облученный протонами с энергией 0.07, 0.10, 0.13,

0.20 и 1.2 ГэВ

2.4.1.5 Объединенное сравнение.

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

2.5 Верификация программ по интегральным экспериментам.

2.6 Анализ сравнений по интегральным экспериментам.

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

3 НЕЙТРОНО-ФИЗИЧЕСКИЕ ПАРАМЕТРЫ И ПАРАМЕТРЫ ЭНЕРГОВЫДЕЛЕНИЯ ИССЛЕДОВАННЫХ МИШЕННЫХ ЧАСТЕЙ ЭлЯУ.

3.1 Методика расчета.

3.1.1 Используемые программы.

3.1.2 Расчет нейтронного выхода.

3.1.3 Расчет энерговыделения

3.1.3.1 Расчет парциальных энерговыделений.

3.1.3.2 Трудности расчетов локального энерговыделения

3.2 Жидкосвинцовая и свинцово-висмутовая мишень.

3.3 Псевдо-сжиженная свинцово-тяжеловодная композиция.

3.4 Газоохлаждаемая вольфрамовая мишень.

3.5 Выбор окна.

3.6 Выводы.

1 Назначение ЭлектроЯдерных Установок (ЭлЯУ).

Несмотря на сложившееся после Чернобыльской аварии негативное, в целом, отношение к себе мировой общественности, ядерная энергетика продолжает оставаться одним из наиболее перспективных источников энергии и имеет серьезные предпосылки увеличить свой вклад в энергообеспечение планеты. Основными причинами сдерживания развития ядерной энергетики до сих пор остаются проблемы (1) безопасного функционирования АЭС, (2) переработки радиоактивных отходов, и (3) нераспространения делящихся материалов [1]. Для государств, обладающих ядерным оружием, в последние годы появилась новая проблема - уничтожение или невоенное использование (конверсия) избыточного количества оружейного плутония и высокообогащенного урана. Существуют два основных направления решения проблемы радиоактивных отходов и конверсии плутония. Первое, традиционное направление, заключается в долговременном захоронении долгоживущих отходов и плутония и, соответственно, откладывания проблемы их уничтожения на поколения вперед. Второе направление заключается в стремлении переработать (трансмутировать) отходы и плутоний в настоящее время или в ближайшем будущем, тем самым избавляя последующие поколения как от самой проблемы отходов, так и от проблемы нераспространения, из-за риска экологических катастроф от накопленных отходов наработки ядерной энергетики в объеме, по крайней мере, в десятки тысяч реакторо-лет

Как указано в [2], к концу 1995 года общемировой опыт ядерной энергетики составил 7696 реакторо-лет. Существующие 437 ядерноэнергетические установки в более чем 30 странах ежегодно выдают более 10,000 тонн отработанного топлива, из которых только 3000 тонн перерабатываются [3] При этом, согласно оценке [4], из 100000 тонн отработанного топлива, находящегося во временных хранилищах, примерно 1000 тонн составляет плутоний. Благодаря накоплению в отработанном топливе значительной части продуктов деления, данный плутоний не влияет, на первый взгляд, на проблему нераспространения. Однако, так как отсутствуют принципиальные проблемы по химическому выделению плутония из отходов имеетпозволяет трансмутировать изотопы, считающиеся в обычных реакторах отходами, т.е. работать с топливом, непригодным в обычных критических реакторах.

Как отмечается в [8], кроме основных, ЭлЯУ имеют и другие потенциальные достоинства:

1. Полное выжигание топлива;

2. Отказ в топливном цикле от процесса обогащения, экстракции и рефабрикации;

3. Отсутствие накапления плутония;

4. Выжигание текущих отходов, а также отходов, полученных в прошлом;

5. Повышенный уровень безопасности;

6. Защищенность от диверсий и несанкционированного распространения;

7. Сравнимость финансовых затрат с альтернативными концепциями;

8. Использование готовых (разработанных) технологий.

В [9], наряду с достоинствами ЭлЯУ в сравнении с критическим системами, обозначены также и некоторые недостатки, которые обусловлены введением нетрадиционных для ядерныой энергетики элементов:

1. Проблемы согласования ускорителя с реактором;

2. Высокое удельное тепловыделение в мишени.

3. Тепловые нагрузки в аварийных ситуациях в модулях ускорителя и в транспортных каналах;

4. Случайное повышение тока приводит к случайной повышенной энергонапряженности в б л анкете;

5. Нестабильность работы ускорителя является причиной частого срабатывания СУЗ;

6. Повышенная выдача продуктов деления,

7. Повышенный тепловой выброс в атмосферу.

Кроме трансмутации отходов, рассматривается использование ЭлЯУ в исследованиях, направленных на:

1. Создание ядерных энергетических установок, основанных на использовании ториево-го цикла, отходы которого обладают значительно меньшим уровнем экологической токсичности [10]. Имеется ряд проектов производства энергии с помощью ускорителя (Accelerator-Driven Energy Production (ADEP)) [11, 12, 13], которые предполагают не только извлечение энергии деления из тория, но и уничтожение долгоживущих отходов, причем без производства материалов оружейного характера.

2. Создание мощных нейтронных генераторов для специализированных задач, таких как производство трития (например, проект Accelerator Production of Tritium (APT) [14]), использование нейтронов в материаловедческих исследованиях (например, в проекте Spallation Neutron Source (SNS) [15]), легирования материалов[16], производство изотопов [17] и т.д.

В ЭлЯУ происходят принципиально новые с точки зрения традиционной реакторной физики процессы взаимодействия высокоэнергичных ( с энергией порядка 1ГэВ) протонов с материалами специального устройства - мишенью [18]. К мишенному устройству ЭлЯУ, которое служит нейтронным источником в подкритической системе, предъявляются следующие требования [19]:

1. Высокая нейтронная производительность;

2. Способность работать при высоких энерговыделениях порядка нескольких десятков МВт;

3. Компактные размеры для введения его внутрь активной зоны;

4. Надежное функционирование, не требующее интенсивного обслуживания;

5. Безопасность;

6. Малый уровень образования радиоактивных отходов.

Стремление к увеличению выхода нейтронов из мишени предполагает использование в качестве мишенных материалов тяжелые элементы, как неделящиеся - свинец, свинецвисмут [20], ртуть [21]2, так и материалы, включающие соли делящихся актинидов [23, 24]. Возможны и комбинированные мишени с добавками продуктов деления, подлежащих трансмутации - 90Sr, 137Cs [25].

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

В последнее десятилетие во многих странах проведены десятки крупных международных конференций, посвященных проблемам ЭлЯУ, управляемым ускорителем [27]-[43].Работы по нескольким трансмутационным проектам уже идут во многих ядерных центрах США, Европы, России, Японии, Китая и Кореи (что отражено в обзорах [44, 45, 19] и в специальных докладах на последних конференциях по трансмутации, проведенных в 1996г. в Швеции [36] и в 1999г. в Чехии [37]).

Среди этих проектов особое место занимают следующие:

1. Проект OMEGA, разрабатываемый в Японском институте Атомной Энергии (JAERI).

2. Американские проекты, разрабатываемые в: a) Лос-Аломосской лаборатории (LANL) - проекты ATW, APT, ADEP, ABC; b) Брук-Хевинской лаборатории (BNL) - проекты ADAPT, PHOENIX; c) Оак-Риджской лаборатории (ORNL) - проект SNS.

3. Европейские проекты SINQ, ESS и проект под руководством К.Руббиа (C.Rubbia) [12].

4. Российские исследования проблем трансмутации наиболее полно продемонстрированы в завершившимся в 1996г. проекте МНТЦ-017 [46], результатом которого явились,

2Неделящихся в реакторной области энергий. При высоких энергиях (десятки МэВ и выше) эти материалы, так же как и актиниды, испытывают деление (см., например, [22]). в частности, конструкторские разработки нескольких вариантов мишенных частей ЭлЯУ (МЧЭлЯУ). Кроме того, важнейшие ядерно-физические и технологические аспекты МУПС прорабатывались (и/или прорабатываются) в проектах МНТЦ-157, 176, 183, 477, 540, 554, 559, 839, 1145, 1314 и др.

Различные проблемы трансмутации радиоактивных отходов обсуждались практически на всех последних международных [47, 48, 49, 50] и национальных (см. например, [51, 52, 53, 54, 55, 56, 57, 58, 59] конференциях по ядерным данным.

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

Так же как и в случае ядер-продуктов на тонких мишенях, необходимо сделать вывод, что на данный момент программы не в состоянии достаточно точно предсказать все необходимые параметры. Полученные расхождения (в среднем, в 2 раза) с экспериментальными значениями свидетельствуют о том, что использовать данные программы уместно только для качественных(оценочных) расчетов параметров МЧЭлЯУ.

Глава 3

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

3.1 Методика расчета. 3.1.1 Используемые программы

Для расчета нейтронофизических параметров и параметров энерговыделения МЧЭлЯУ, исследованных в ИТЭФ, были задействованы следующие расчетные программы:

1. НЕТС - программа, генерирующая адронный каскад в мишени;

2. MCNP - программа, ведущая транспорт медленных нейтронов и фотонов;

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

4. SPECT - вспомогательная программа к НЕТС, рассчитывающая параметры непосредственно адронного каскада;

5. Т^ЮБАМТ - двумерная многогрупповая программа, рассчитывающая перенос нейтронов Эдг-методом.

Общая схема стыковки программ представлена на Рис. 3.1.

Рис. 3.1: Схема применения программ для вычисления параметров МЧЭлЯУ.

Как видно их схемы, основной программой является программа НЕТС, которая моделирует адронный-ядерный каскад с энергией частиц выше 20МЭВ1, и выдает три файла событий, описывающих моделируемый каскад:

1. файл рождений низкоэнергичных нейтронов (Еп <20МэВ) и фотонов;

2. файл рождений 7г°-мезонов;

3. общий файл адрон-ядерных взаимодействий.

Первый файл считывается программой МСКР с помощью специально внедренного в нее модуля. Программа МСКР моделирует транспорт нейтронов и фотонов с целью опре

1 Величина энергии, при которой транспорт частиц заканчивается (энергия обрезания), задается опционально для каждого типа частиц. Как правило, для протонов и нейтронов задавалось 20МэВ, для 7г-мезонов - 1МэВ. деления выхода(вылета) низкоэнергичных нейтронов, их спектров, а также энерговыделения за их счет. Файл рождений 7г°-мезонов использовался программой ЕОБ для определения пространственного распределения энерговыделения от жестких фотонов, образующихся при распадах 7г°-мезонов. Общий файл взаимодействий использовался программой БРЕСТ. Авторская версия программы включает расчет пространственного распределения энерговыделения за счет взаимодействий с энергией выше 20МэВ, потоки и спектры адро-нов, выходы вторичных ядер-продуктов (включая спектры легких ядер). Детально расчет нейтронного выхода и энерговыделения описан в разд. 3.1.2 и 3.1.3.

В качестве альтернативы программе МСКР применялась программа Т\УОБАМТ. Используя пространственно-энергетическое распределение рожденных низкоэнергичных нейтронов, программа вычисляет нейтронные потоки, поглощения и выходы. Ввиду относительно высокой скорости счета, применение программы TWODANT предпочтительнее МСЫР, в случаях проведения значительного числа оптимизационных расчетов.

Имеющаяся версия программы НЕТС была модернизирована включением в нее программ Ш11СЬ[131] и СЕМ95[129], которые, в качестве коренному блоку ВЕИИШ программы НЕТС, моделируют адрон-ядерные взаимодействия.

Программа БРЕСТ отличается от стандартной из-за включения в нее дополнительного модуля, который разделяет парциальные энерговыделения в гомогенно заданных компонентах (смесях) (см. п. 3.1.3). Это позволяет рассчитывать гетерогенную среду (например, свинцово-водную композицию) в виде гомогенной среды, получая парциальные энерговыделения (для воды и свинца отдельно). Это исключает трудности с заданием сложных геометрий и сокращает время расчета.

Остальные программы используются в стандартном виде без изменений.

При проведении расчетов поперечная форма пучка была задана в соответствии с рас-счетным моделированием пучка после поворотного магнита [46] (Рис. 3.2).

Геометрия мишеней была задана в соответствии с результатами их предварительной конструкторской проработки (см. Рис. 3.3, 3.5, 3.6 и 3.10).

3.1.2 Расчет нейтронного выхода

Обычно, при расчете нейтронного выхода из "голой" (ничем неокруженной) мишени, нейтронный выход рассчитывается подсчетом числа смоделированных нейтронов (как быI о X I о о.

§1

50 Ч о ь-о с 0 5

10 15 20

Радиус(см)

25

Рис. 3.2: Плотность протонного пучка в зависимости от радиуса. Нормировка на ток ЮОмА. стрых, с £,„>20МэВ, так и медленных), пересекающих внешнюю поверхность мишени.

Однако, такой подход не отражает нейтронофизические свойства мишени в случае, если она окружена каким-либо, особенно, размножающим материалом По причине того, что нейтронные параметры мишени сильно зависят от внешних условий, а именно, от того какая среда окружает мишень. Например, если говорить о свинцовых мишенях, то, в силу наличия у свинца существенного теплового поглощения, они, находясь в тяжеловодном бланкете с мягким нейтронным спектром, будут поглощать значительно больше нейтронов, чем в жидкосолевом бланкете с промежуточным или быстром спектром. Поэтому в конкретных случаях целесообразно говорить не о параметрах мишени как таковой, а о параметрах всей системы " мишень-бланкет". Для того, чтобы сравнивать различные мишени, следует определить параметр, который являлся бы принципиальным для всей установки. Такой параметр должен быть тесно связан с эффективностью ЭлЯУ, под которой, например, для трансмутационной установки следует понимать, число делений загруженных в б л анкет актинидов. Поэтому полезно было бы проводить расчеты параметров всей системы "мишень-бланкет" в целом, вычисляя число делений на один исходный протон. Проводить такие расчеты возможно с использованием программ Монте-Карло типа МСи или МС№, но при расчете в размножающей среде это требует практически недостижимой статистики. Расчет можно значительно упростить, если исключить выход нейтронов при делении актинидов, т.е. моделировать каждое деление как поглощение (без нейтронного выхода). Так как при фиксированном кец соотношение числа делений и поглощений фиксировано для конкретного бланкета, то число захватов, равное сумме числа делений без выхода нейтронов и поглощений, будет пропоционально числу делений, т.е. пропорционально эффективности установки и, в том числе, мишени. Поэтому расчет нейтронного выхода мишени, непосредственно окружаемой размножающим бланкетом, целесообразно проводить, вычисляя число указанных захватов в бланкете, что и было сделано с применением программы МСКР с параметром 1ЧС)Ми=0 (см. описание программы в [83] стр. 3-90), при котором программа трактует деление как захват нейтрона с выходом только фотонов.

3.1.3 Расчет энерговыделения

При расчете энерговыделения учитываются следующие процессы:

• ионизационные потери заряженных частиц (кроме легких ядер) - протонов, заряженных 7г-мезонов, и продуктов их распада (мюонов);

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

• энерговыделение за счет распада нейтральных пионов.

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

1. отрезок делился на равные подотрезки, в количестве, при котором в каждой области получалось несколько подотрезков2;

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

3. значения ионизационных потерь на под отрезках нормировались на полное значение ионизационных потерь на отрезке.

2Минималыюе количество подотрезков в 11^-области задается опционально в задании для программы БРЕСТ.

Энерговыделение от ядерных взаимодействий относилось к месту нахождения соответствующего ядра (считается, что пробеги образующихся легких ядер много меньше по сравнению с протяженность адронного каскада). Пространственно распределение энерговыделение от гамма-квантов, образующихся в результате распада 7г°-мезонов моделировалось программой ЕОБ.

3.1.3.1 Расчет парциальных энерговыделений.

Разделение энерговыделения на парциальные составляющие основывается на том, что:

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

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

3. при разделение энерговыделения от медленных нейтронов (£^<20МэВ) использовались пространственные распределения групповых потоков и групповые терма-факторы материалов компонент.

Качественного распределение энерговыделения за счет нейтральных 7г°-мезонов по компонентам не проводилось. Упрощенно, энерговыделение приписывалось той компоненте, в ядерном взаимодействии с ядром которого рождался конкретный нейтральный пион. Считалось, что такое упрощение не значительно снижает точность расчетов, т.к. максимальная доля локального энерговыделения за счет 7г°-мезонов никогда не поднимаШсг выше 10-=-12% (см. Табл. 1.3).

3.1.3.2 Трудности расчетов локального энерговыделения.

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

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

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

3.2 Жидкосвинцовая и свинцово-висмутовая мишень.

Этот вариант мишени разработан для использования в жидкосолевом бланкете [60]. Материалом мишени служит расплавленный свинец, отделенный от ускорительного вакуума стальным окном конической формы.

С целью оптимизации диаметров мишени и пучка протонов был сделан ряд расчетов жидкосолевого варианта при различных значениях диаметров мишени и протонного пучка (Рис. 3.3).

Результаты представлены в Таблице 3.1 и на Рис. 3.4, из которых следует, что наибольшим числом поглощений в бланкете обладает вариант с диаметром мишени 40см и диаметром пучка 30см. При больших размерах мишени увеличивается число паразитных захватов нейтронов мишенью, а при меньших диаметрах пучка максимальное удельное энерговыделение в мишени становится крайне высоким (>1кВт/г). Поэтому данное соотношение диаметров следует считать наилучшим. Полезный выход нейтронов при нем составляет ЗО.Зп/р при энергии протонов 1 ГэВ и является практически максимально возможным при данной энергии протонного пучка. Столь высокий выход нейтронов стал возможен благодаря отсутствию в системе сильнопоглощающих конструкционных материалов и солевом бланкету, поддерживающему достаточную жесткость нейтронного спектра, благодаря чему обеспечивается низкий уровень нежелательного поглощения свинцом мишени.

Другой вариант мишени, на основе свинцово-висмутовой эвтектики, был предложен женных выше, представлены в Таблицах 3.1 и 3.2.

ЗАКЛЮЧЕНИЕ

1. Создан взаимосвязанный комплекс программ для расчета нейтронофизических параметров толстых мишеней, используя программы НЕТС, МСКР, БРЕСТ, ЕвБ, ТЭДШАОТ, баз данных ЕШР/В и БНАБ.

2. Показана удовлетворительная (~15%), по большинству программ, точность расчета полного нейтронного выхода и полного тепловыделения. Используя экспериментальные данные показано, что данная точность соответствует по порядку величины экспериментальной точности. Показана менее удовлетворительная (порядка 30%) точность расчетов локального тепловыделения, в первую очередь, за счет продутков ядерных взаимодействий. Показана неудовлетворительная точность (фактор 2 и выше) расчета выходов ядер-продуктов.

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

4. Проведен статистический анализ адронного каскада в толстой мишени, в результате показаны диапазоны энергий, принимающих наибольшее участие в межьядерном каскаде, и, следовательно, наиболее приемлемые для верификации программ. Получено, что наиболее существенный вклад в динамику адронного каскада играют диапазоны энергий: Е0-0.25 ГэВ - Е0 (Е0 - исходная энергия налетающего протона) и 0.1 - 0.25 ГэВ (Рис. 2).

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

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

7. Проведено сличение экспериментальных и расчетных данных по 11 программам по выходам ядер-продуктов. Полученные среднестатистические отклонения расчетных значений выходов ядер-продуктов от экспериментальных данных на уровне фактора 2.0 позволяют сделать вывод о недостаточно высоком уровне предсказательной способности используемых программ и, как следствие, о необходимости их совершенствования, прежде чем они станут надежным инструментом предсказаний параметров ЭлЯУ. Пример программы CASCADE, в случае 208РЬ - 0.13ГэВ, и программы СЕМ95, в случае 63Си -0.13ГэВ, когда расчетные результаты "совпали" с экспериментальными с требуемой точностью 30% - подтверждает возможность достижения требуемого уровня предсказательной способности.

8. Произведены расчеты параметров мишенных устройств, предлагаемых для проектирования подкритических систем с целью их оптимизации. Были проанализированы две жидкосвинцовые, свиноцово-висмутовая, вольфрамовая и свинцово-тяжеловодная мишени, предложенных для проектной реализации в ходе ряда проектов МНТЦ. В ходе оптимизаций получено:

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

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

- показана нейтронопроизводительность свинцовой мишени охлаждаемой тяжелой водой, как функция ее пористости (Рис. 8);

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

Также, исследована проблема первой стенки (окна) по тепловыделению. Показаны зависимости удельных тепловыделений при использовании конусного окна.

1. J.Kupitz, V.Arkhipov, "Emerging Nuclear Systems for Energy Production and Transmutation", Ссылка 36], 132-143.

2. S. Bjurstrôm, "Radioactive Waste Management in Sweden" , Ссылка 36], 116-123.

3. C.Mileikowsky, "Eliminating of risk of nuclear bombs from civil plutonium", Ссылка 36], 249-252.

4. V.V. Orlov, "Lead-cooled fast reactor Concept", Ссылка 36], 124-131.

5. T.Mukaiyama, T.Ogawa, M.Mizumoto, T.Takizuka, R.Hino, and Y.Oyama, "Omega Program & Neutron Science Project for Development of Accelerator Hubryd System at JAERI", Ссылка 40]; N.Watanabe," JAERI 5 MW Spallation Source Project", Ссылка [42] pp. 3-10.

6. F.Vennery, "The Physics Design of Accelerator-Driven Transmutation Systems ", Ссылка 33], 117-137.

7. P.Wydler, E.Curti, S.Pelloni, G.Yonirou, "Role of Hybrid Systems in Waste Transmuation", Ссылка 40]

8. С. Bowman, "Basis and Objectives of Los Alamos Accelerator-Driven Transmutation Technology Project", Ссылка 33], 22-43

9. A. N. Sisakian, "Elaboration of Physical Conception of Electronuclear Power Generation," p. T/B 27 in Ref. 36].

10. Modular Plant Design and Supporting Engineering Development and Demonstration Work", Ссылка 42. pp. 186-193; см. также http://strider.lansce.lanl.gov.

11. J.R. Alonso, "The Spallation Neutron Source (SNS) Project: Introduction", Ссылка 42] pp. 197-201.

12. V.G.Alpatov, A.V.Davydov et al. "Excitation of long living isomers 107m, 109Ag in the fast neutron inelastic scattering reactions", 1-st Int. Workshop on induced 7-emission, Predeal, Romania, 17-20 августа 1997.

13. G.J.Russel, E.J.Pitcher and L.L.Daemen, "Introduction to Spallation Physics and Spallation Target Design", Ссылка 33], 93-104.

14. W. Gudowski, "Accelerator-Driven Systems Survey of the research programs in the world", Ссылка 38], 3-25.

15. B.F. Gromov, E.I.Yefimov, A.A.Veremeev, M.P.Leonchuk, V.T.Gorshkov, "Precompound Design of Liquid Metal Targets with Diafragm for Accelerator-Based System", Ссылка 36], 822-827.

16. G.S.Bauer, "Research and Development for Molten Salt Metal Targets", Ссылка 36], 803-814

17. V.P.Eismont, A.I.Obukhov, A.V.Prokofiev, A.N.Smirnov, "An Experimental Database on Proton-Induced Fission Cross sections of Tantalum, Tungsten, Lead, Bismuth, Thorium and Uranium", Ссылка 36], 592-598.

18. T.Takizuka, T.Nishida, M.Mizumoto and H.Yoshida, "JAERI R&D on Accelerator-Based Transmutation Under OMEGA Programm", Ссылка 33], 64-73.

19. T.Takizuka, T.Nishida, T.Sasa, "Conceptual Design Study by Accelerator-Driven Systems for Nuclear Waste Transmutation", Ссылка 36], 179-185.

20. V.G.Semenov, N.M.Sobolevsky, "Computer Study of 90-Sr and 137-Cs Transmutation by Proton beam", Ссылка 36], 353-359.

21. J. P. Schapira, "Transmutation of Nuclear Wastes, A working report to Nupecc," I.P.N. Report IPNO DRE 94-04, Orsay (1994).

22. P.P.Blagovolin, V.D.Kazaritsky, V.F.Batyaev, E.A.Zolotareva, E.S.Nikolaevsky, V.T.Gor-shkov, V.A.Shulyndin, "Accelerator-Driven Molten-Fluoride Reactor with Modular Heat Exchangers on Pb-Bi Eutectic", Ссылка 38], 235-240.

23. A.V.Lopatkin, V.M. Matyushechkin, I.T.Tretyakov, P.P.Blagovolin, V.D.Kazaritsky, V.I.Kostenkov, D.Yu.Chuvilin, "Lead-Cooled Heterogeneous Accelerator Driven Molten-Fluoride Blanket for Inceneration of Long-Lived Radiation Waste", Ссылка 36], 878-884.

24. В.Т.Горшков, В.А.Шулындин, отчет по Проекту МНТЦ- 017.

25. G.S.Bauer, F.Atchison, T.A.Broome, H.M.Conrad, "A Target Development Program for Beamhole Spallation Neutron Sources in the Megawatt Range", Ссылка 33], 105-116.

26. G.J. Van Tuyle et al. "Topical Report on a Preconceptual Design for the Spallation-Induced Lithium Conversion (SILC) Target for the Accelerator Production of Tritium", BNL Report BNL-52401 (1993).

27. G.S.Bauer, "The European Spallation Source Study, ESS", Ссылка 36], 159-168; G.S.Bauer, "The ESS as the 21st Century User Facility for Neutron Scattering in Europe", Ссылка [42].

28. В.Д.Казарицкий, П.П.Благоволин, Е.А.Золотарева, В.Ф. Батяев, "Бланкет электроядерной установки для дожигания долгоживущих радиоактивных отходов ядерной энергетики", Ссылка 43].

29. С.В.Керсновский, В.И.Костенков, В.М.Новиков, А.В.Пузырев, Д.Ю.Чувилин, П.П.Благоволин, В.Д.Казарицкий, Ю.В.Кочевалин, " Электроядерне установки для пережигания актинидов на основе расплавов солей", ИАЭ-5989/3 (1996).

30. Д.Ю.Чувилин, В.И.Костенков. "Интегральный нейтроно- физический эксперимент в обоснование проекта бланкета электроядерной установки на расплавах фтористых солей". ИАЭ-5979/2 (1996).

31. V.R.Mladov, M.L.Okhlopkov, V.F.Batyaev,"Prospects of designing a proton target with ftuidized lead particle bed in heavy water", Ссылка 36], 885-890.

32. G.J. Van Tuyle, "Technologies Using Accelerator-Driven Targets under Development at BNL", Ссылка 33], 74-82.

33. Westinghouse Hanford Company, "Final Safety Analysis Report for the Fast Flux Test Facility", HEDL-TI-75001, December 1975.

34. G.J.Van Tuyle, M.Todosow, M.J.Geiger, A.L.Aronson, H.Takahashi, Accelerator-Driven Subcritical Target Concept for transmutation of Nuclear Wastes, Ссылка 28], 98-119.

35. M.A.Lone, "Data Needs for Accelerator-Based Neutron Radiography Sources", Intermediate Energy Nuclear Data for Applications, (Proc. of the Advisory Group Meeting organized by the International Atomic Energy Agency) INDC(NDS)-245 (1991), 57-60.

36. D.Filges, P.Cloth, R.-D.Neef, H.Schaal, "Monte-Carlo Simulation of Accelerator Reactor System", Ссылка 38], 117-131.

37. Б.С.Сычев А.Я.Серов, Б.В.Манько, "Аналитическая аппроксимация дифференциальных сечений образования вторичных частиц в неупругих нуклон-ядерных взаимодействиях при энергиях выше 20 МэВ", Препринт МРТИ-799 (1979).

38. В.П. Крючков, О.В. Суманеев, MOSKIT1 программа для рещения задачи переноса частиц (n,p,7r,j) в веществе методом монте-карло. Препринт ИФВЭ 92-132.

39. Д.В. Горбатков, В.П. Крючков, В кн.: Тринадцатое совещание по ускорителям заряженных частиц. Аннотации докладов. Дубна, 1992, с. 195.

40. T.W. Armstrong and К.С. Chandler, Nucl. Sci. Eng. 49 (1972) 110.

41. P.A.Aarnio et al., Report CERN TIS-RP 168 (1986).

42. M.B. Emmett, Report ORNL 4972 (1975).

43. J.F. Briesemeister, "MCNP: A General Monte-Carlo N-Particle Transport Code, Version 4A", LA-12625-M (1993).

44. Jl. Абагян и др. "Программа MCU для расчетов ядерных реакторов. Верификация", Препринт ИАЭ 5751/5 (1994).

45. R.D.O'Dell and F.W.Brinkley, Jr. and D.R.Marr, "User's Manual for ONEDANT: A Code Package for One-Dimensional, Diffusion-Accelerated, Neutral-Particle Transport", LANL manual LA-9184-M (1982).

46. R.E.Alcouffe, F.W.Brinkley, Jr., D.R. Marr, and R.D. O'Dell, "User's Guide for TWODANT: A Code Package for Two-Dimensional Diffusion-Accelerated, Neutral-Particle Transport", LANL manual LA-1049-M, Rev. (1984).

47. R.E.Alcouffe et al., Report LA-10049-11 (1990).

48. W.R.Nelson et al., Report SLAC-265 (1985).

49. T.A.Gabriel et al., Report ORNL/TM-IIO6O (1989).

50. P.Cloth et al., "HERMES: A Monte-Carlo Program Sytem for Beam-Materials Interaction Studies", KFA Report Jiiel-2203 (1989).

51. R.E.Prael and H.Lichtenstein, "User Guide to LCS: The LAHET Code System", Report LA-UR-89-3014 (1989).

52. F.Atchison, Report NEA NSC/DOC/95.2 (1995).

53. O.Bersillon et al., "TIERCE: A code system for Particles and Radiation Transport in Thick Targets", Ссылка 36], 520-526.

54. A.V.Dementyev, N.M. Sobolevsky, "SHIELD a Monte-Carlo Hadron Transport Code" Proc. of Specialists Meeting "Intermediate Energy Nuclear Data: Models and Codes", Paris May 30 - June 1, 1994, NEA OECD, Paris, p.237; Препринт ИЯИ 0874, M. 1994.

55. A.V.Daniel, "Transport Code SITHA", Khlopin Radium Institute Preprint 181 (1984).

56. W.A. Coleman and T.W. Armstrong, Nucl. Sei. Eng. 43 (1971) 353

57. T.Nishida, T.Sasa, H.Takada, T.Takizuka, "Development of the Code System ACCEL for Accelerator Based Transmutation Research", Ссылка 36], 668-674.

58. A.J.Koning, "Requirements for an evaluated nuclear data file for accelerator-based transmutation", ECN-C-93-005, ECN-C-93-041 (1993).

59. P.G. Young, W.B.Wilson and M.B.Chadwick, "Nuclear Data Requirements for Accelerator-Driven Transmutation Systems", Ссылка 33], 703-709.

60. R.Michel, "Nuclide Production at Intermediate Energies", 8eme Journees Saturne, 5-6 May 1994, Saclay

61. R. Michel and P. Nagel, International Codes and Model Intercomparison for Intermediate Energy Activation Yields, NEA/OECD, Paris, 1997, NSC/DOC(97)-l; см. также htt p : / / www. nea. fr/html/science/pt / ieay.

62. Vasilkov R.G. "Neutron Emission from an Extended Lead Target under Action of Light Ions in the GeV region". Report KEK-90-25 v.I p.340.

63. N.M.Sobolevsky, "International Code Comparison for Intermidiate Energy Nuclear Data. Thick-target benchmark for Lead and Tungsten". Report NSC/DOC(95)2.

64. P.R.Tunnicliffe, B.H.Chidley, J.S.Eraser, "High Current Proton Linear Accelerators and Nuclear Power", Int. Conf. on Accelerators in Chalk River. Ontario, 1976.

65. В.С.Барашенков, "Ядерно-Физические Аспекты Электроядерного Метода", Физика элементарных частиц и атомного ядра, 9 вып. 5 (1978) с. 871-921.

66. M.Gloris, R.Michel, U.Herpers, F.Sudbrock, D.Filges, "Production of residual nuclei from irradiation of Thin Pb targets with protons upto 1.6 GeV", Nucl. Instr. and Meth., ВИЗ (1996) 429-433.

67. Belyakov-Bodin V.I. et al. "Heat deposition in targets bombarded by medium-energy protons", Nucl. Instr. and Meth. in Phys. Res. A335 30-36 (1993).

68. Yu.E. Titarenko, S.G. Mashnik, V.D.Kazaritsky, V.F. Batyaev et al. "Exerimental and Theoretical Study of the Yields of Radionuclides Produced in 209Bi thin target Irradiated by 1500 MeV and 130MeV Protons", Nucl. Instr. and Meth. A414(l) (1998) 73-99.

69. C.L.Dunford and T.W.Burrows, "National Nuclear Data Center online nuclear data service", BNL report NNDC/BNL-92/08.

70. NEA report DBG-006.2, "EXFOR network service" (1989); H.D.Lemmer (editor) NDS EXFOR Manual, IAEA-NDS3 Vienna, IAEA (1984).

71. O. Bersillon, T. Gabriel and S.G. Mashnik, "Cascade-exiton model detailed analysis of proton spallation at energies from lOMeV to 5 GeV", Bruyere-le-Chatel (1996)

72. M. Blann, H. Gruppelar, P. Nagel, and J. Rodens, International Code Comparison for Intermediate Energy Nuclear Data, NEA OECD, Paris (1994).

73. P. Nagel, J. Rodens, M. Blann, and H. Gruppelar, "Intermediate Energy Nuclear Reaction Code Intercomparison: Application to Transmutation of Long-Lived Reactor Wastes," Nucl. Sei. Eng., 119, 97 (1995).

74. R. Michel, R. Bodemann, H. Busemann, R. Daunke, M. Gloris, B. Klug, A. Krins, H. -J. Lange, I. Leya, M. Lüpke, S. Neumann, H. Reinhard, U. Herpers, Th. Schiekel, F.

75. Bowman et al., "Nuclear Energy Generation and Waste Transmutation Using an Accelerator-Driven Inense Thermal Neutron Source". Nucl. Instr. and Meth. in Phys. Res. A320, p. 336 (1992).

76. H.Takahashi, Survey of codes Relevant to Design, Engineering and Simulation of Transmutation of actinides by Spallation, OECD/NEA report 1993.

77. K.K. Gudima, S.G. Mashnik and V.D. Toneev, Nucl. Phys. A 401 (1983) 329; JINR Communications P2-80-774 and P2-80-777, Dubna (1980).

78. G.A. Lobov, N.V. Stepanov, A.A. Sibirtsev and Yu.V. Trebukhovskii, ITEP Preprint ITEP-91, Moscow (1983); A.A. Sibirtsev, N.V. Stepanov and Yu.V. Trebukhovskii, ITEP

79. Preprint ITEP-129, Moscow (1985); N.V. Stepanov, ITEP Preprint ITEP-81, Moscow (1987); N.V. Stepanov, ITEP Preprint ITEP-55-88, Moscow (1988) (in Russian).

80. M. Blann, Phys. Rev., C54, No.3 (1996) 1341-1349.

81. P.G. Young, E.D. Arthur and M.B.Chadwick, Los Alamos National Laboratory Report LA-12343-MS (1992); M.B. Chadwick and P.G. Young, Phys. Rev. С 47 (1993) 2255.

82. Т. Nishida, Y. Nakahara, and T. Tsutsui, "Development of a Nuclear Spallation Simulation Code and Calculations of Primary Spallation Products", JAERI-M 86-116 (1986).

83. K. Niita et al., Phys. Rev., C52, 2620 (1995).

84. M.Foshina, J.B.Martins, O.A.P.Tavares and V.di Napoli, "Systematics of Spallation Yields with Four-Parameter Formula", Radiochimica Acta, 35 (1984) 121-131.

85. M. Blann, Phys. Rev. С 54 (1996) 1341.

86. В.С.Барашенков, В.Д.Тонеев, Взаимодействие высокоэнергичных частиц и ядер с атомными ядрами, Атомиздат, Москва, 1972.

87. В.С.Барашенков, А.С.Ильинов, Н.М.Соболевский, В.Д.Тонеев, Успехи Физ. Наук 109 (1973) 91.

88. A.S. Iljinov, M.V. Kazarnovsky, E.Ya. Paryev, Intermediate-Energy Nuclear Physics (CRC Press, Boca Raton, Florida, 1994).

89. H.W. Bertini and M.P. Guthrie, Nucl. Phys. A 169 (1971) 670; H.W. Bertini, Phys. Rev. 188 (1969) 1711.

90. K. Chen, Z. Frankel, G. Fridlandder, J.R. Grover, J.M. Miller and Y. Shimamoto, Phys. Rev. 166 (1968) 948; K. Chen, G. Ffiedlander, G.D. Harp and J.M. Miller, Phys. Rev. С 4 (1971) 2234.144145146147148149150151152153154155

91. V.S. Barashenkov, H.W. Bertini, K. Chen, G. Friedlander, G.D. Harp, A.S. Iljinov, J.M. Miller, V.D. Toneev, Nucl. Phys. A 187 (1972) 531.

92. К.К.Гудима, Г.А.Осоков, В.Д.Тонеев, Ядерная Физика 21 (1975) 260 Sov. J. Nucl. Phys. 21 (1975) 138.

93. S.G. Mashnik and V.D. Toneev, JINR Communication P4-8417, Dubna (1974).

94. S.G. Mashnik, A.J. Sierk, 0. Bersillon and T. Gabriel, Los Alamos National Laboratory Report LA-UR-97-3176 (1997).

95. H.J. Krappe, J.R. Nix and A.J. Sierk, Phys. Rev. С 20 (1979) 992. A.G.W. Cameron, Can. J. Phys. 35 (1957) 1021.

96. V.S. Barashenkov, A.S. Iljinov, V.D. Toneev, F.G. Gereghi, Nucl. Phys. A 206 (1973) 131.

97. A.S. Iljinov, M.V. Mebel, N. Bianchi, E.De Sanctis, C. Guardalo, V. Lucherini, V. Muccifora, E, Polli, A.R. Reolon and P. Rossi, Nucl. Phys. A 543 (1992) 517.

98. J.W. Truran, A.G.W. Cameron and E. Hilf, in: Proc. Int. Conf. on the Properties of Nuclei Far From the Region of Beta-Stability, Leysin, Switzerland, 1970, v. 1, p. 275.

99. V.S. Barashenkov, A.S. Iljinov, V.D. Toneev and F.G. Gereghi, Nucl. Phys. A 222 (1974) 204.

100. P. Fong, Statistical theory of nuclear fission (Gorgon and Breach Science Publishers, Ney York, 1969).

101. M.P. Guthrie, ORNL-TM-3119, Oak Ridge National Laboratory (1970).

102. Radiation Scielding Information Center, "HETC Monte Carlo High-Energy Nucleon-Meson Transport Code", Report CCC-178, Oak Ridge National Laboratory (August 1977).

103. Y. Yariv and Z. Frankel, Phys. Rev. C 20 (1979) 2227; Phys. Rev. C 24 (1981) 488.

104. M.B. Chadwick, P.G. Young, P. Oblozinski and A. Marcinkowski, Phys. Rev. C 49 (1994) R2885.

105. M.V. Mebel, A.S. Iljinov, C. Grandi, G. Reffo and M. Blann, submitted Nucl. Inst. Meth. (1997); Proc. Int. Conf. Nucl. Data for Science and Tech., Trieste, Italy May, 1997.

106. Н.Г.Гусев, П.П.Дмитриев "Радиоактивные цепочки. Справочник", 2-е издание. Москва, Энергоатомиздат, 1988

107. J.R.Letaw et al., Astrophys.J.Suppl. 51, 271 (1983).

108. S.Pearlstein, Astrophys.J. 346, 1049 (1989).

109. A.V. Prokofiev, S.G.Mashnik, A.J.Sierk, " Cascade-Exciton Model Analisys of Nucleon-Induced Fission Cross Sections of Lead and Bismuth at Energies from 45 to 500 MeV", LA-UR-98-0418 PrePrint (1998), submitted to Nucl. Sci and Eng., nucl-th/9802027.

110. R.Kinsey (compiler), "ENDF/B Sumary Documentation", BNL report BNL-NCS-17541 (ENDF-201), 3rd edition (ENDF/B5) (1979).

111. Л.П.Абагян, Н.О.Базазянц, М.Н.Николаев, A.M.Цибуля, "Групповые нейтронные константы для расчетов реакторов и защит", Атомиздат, Москва, 1981.

112. V.F. Batyaev, "Neutronics and Power Deposition Parameters of the Targets Proposed in the ISTC Project 17", Ссылка 38], 216-221.