Вычислительные методы Опции

[AtomInfo.Ru]

А речь, видимо, о том, что у М.Н. есть в наличии только некоторые данные пусковых измерений Калинина-3, результаты обработки которых он несколько раз докладывал.

Пропустил.

А что, уже есть другие тесты? Калину-3 и Армянскую станцию мы считали и сами, когда ещё работали по специальности. Появилось уже что-то другое?

[AtomInfo.Ru]

- почему нейтронно-физические коды у нас создаются почти ремесленными методами разрозненными, почти герметичными группами?

Ну потому что мы же все умные

Простейший пример. Оболочки и интерфейсы. Не поверите, но у меня был собственный аналог ShIPR. Назывался Barsic, написан на TurboPascal, потом на Delphi. С языком управления заданиями, на минуточку Там было всё, что положено, в этом языке - циклы, переменные, подпрограммы и т.д. И даже с закосом на объектоориентированность.

Кто-нибудь, кроме меня самого, им пользовался? По факту - нет. Или под моим непосредственным присмотром. Даже мой лучший друг, а теперь соучредитель AtomInfo - впихнуть его код была война: "А чё это мой код будет под твоей оболочкой работать? У него своя оболочка есть". А уж согласование форматов данных - это холодная война на долгие годы.

Ввести в практику чужой расчётный код? Ага, его автор(ы) на тестовом примере вылизали. А чуть изменишь исходные данные, так программа делит на нуль. А в тех же случаях, когда программу хоть как-то тестировали, выяснялось, что она слишком простая для решаемых задач. Например, мифисты замахнулись в 90-ые на российский WIMS (если помните такую - GETERA). Так в ней даже Bj-тиков не было.

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

ИБРАЭ сейчас взялся за единый национальный код для быстрых реакторов. К С-нёву у меня вопросов нет. Но посмотрите, как недовольно забурчало расчётчицкое сообщество - мол, а почему это код будет собирать ИБРАЭ?

- почему не появляется хорошей отечественной площадки для встречи наших расчетчиков и создателей кодов? (у Нейтроники, по-моему, не получается)

Обижаете. Нейтроника - вполне приличный семинар. Я на нём отдыхаю. И не только потому, что приезжают одногруппники, одно/сокурсники и т.д. Он и вправду приятно отличается от многих других подобных мероприятий.

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

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

Украинцы вообще никуда не ездят. На их фоне мы выигрываем с разгромным счётом. Тот же С-нёв ездит на Физор

Я не в восторге от забугорных митингов по расчётчицким делам - много откровенной фигни. Что до кооперации - ну, давайте, для начала попробуем хотя бы просто купить MCNP для коммерческого использования. Рассказать, как DoE одной рукой требует для совместного проекта использовать "признанный в мире" код, а потом другой отвергает заявку на приобретение. т.к. россияне не имеют права задействовать MCNP в коммерческих проектах? Какая кооперация может быть в этом случае? Только с Минском по константам.

- где российские коды - законченные рыночные продукты (типа Helios, и вообще Studsvikовские, или MCNP, т.п.)?

MCU не так плох.

Касательно рыночных продуктов, ждал на этой Нейтронике, что расскажет С-нёв и ИБРАЭ. А они чего-то не приехали. Разрабатывают, видимо

[alex_bykov]

Ну потому что мы же все умные

Простейший пример. Оболочки и интерфейсы. Не поверите, но у меня был собственный аналог ShIPR. Назывался Barsic, написан на TurboPascal, потом на Delphi. С языком управления заданиями, на минуточку Там было всё, что положено, в этом языке - циклы, переменные, подпрограммы и т.д. И даже с закосом на объектоориентированность.

А у нас смешнее, расчётов идёт слишком много - до написания user-friendly оболочки руки не доходят. Пока код только спортирован по форматам файлов с нашей системой ВРК, чтобы возни меньше было...

[AtomInfo.Ru]

А у нас смешнее, расчётов идёт слишком много - до написания user-friendly оболочки руки не доходят. Пока код только спортирован по форматам файлов с нашей системой ВРК, чтобы возни меньше было...

И это тоже всегда было.

С одной стороны, была в моё время куча народу, которая дурью маялась - ну, оболочки и пр. А с другой, был великий Унитазус которым реально считался ВВЭР, но у которого ввод был форматный - типа I3, E12.0.

[MVS]

Было бы полезно, я думаю.

Программа была сериальная, стояла задача распараллелить. По модели распараллелить с MPI (OpenMP было менее выгодно) и сразу вылезли глюки - выходные циферки сильно отличаются от референса. Ошибок в MPI нет, но не работает даже с -O0.

Дальше обезьянья работа - но итог: в коде есть суммирования и результат зависит от порядка слагаемых. MPI ориентируется на топологию кластера и складывает (AllReduce) не в том порядке, как в сериальном коде...

Насчет OpenMP - это не "перспективная" техноголия, ей уже давным-давно пользуются. В гибридных запусках ей цены нет. Глюки - да, бывают, но на то и программисты чтобы с ними бороться.

Но сегодня вы выключаете OpenMP, а завтра MIC или другие карточки будете выкидывать? Или просто хасвелл какой-нибудь придет с новым компилятором и MPI - тоже отмахиваться и сидеть на старых бинарях и 286-х процессорах? Грубо говоря, чтобы считать год на тысячи ядрах, когда за бугром затратят неделю на сотне...

Ремарка: имелись ввиду процессорные ядра, а не ядра атомов

Сообщение отредактировал MVS - 30.11.2012, 21:11

[asv363]

NRC разрешила французам использовать комплекс кодов ARCADIA в США

Группа AREVA сообщила о том, что регуляторы США одобрили использование комплекса кодов ARCADIA для расчётного анализа активных зон и процедуры лицензирования легководных реакторов под давлением в Соединённых Штатах.

Сообщение об этом было опубликовано на сайте американского филиала AREVA 13 марта 2013 года.

Решение регуляторов "существенно расширит" арсенал средств, применяемых группой AREVA при получений лицензий в США, отмечается в пресс-релизе.

Основными компонентами ARCADIA являются спектральный/ячеечный код APOLLO-2A и код для расчёта активной зоны ARTEMIS.

Код APOLLO-2A основывается на известном коде APOLLO-2, созданном в комиссариате по атомной энергии Франции. Группа AREVA адаптировала его для нужд расчётов легководных реакторов под давлением.

Библиотека констант для APOLLO-2A насчитывает 281 группу. В ней, среди прочего, имеются данные по 26 делящимся изотопам и 131 осколку деления. Групповая структура хорошо описывает область энергий ниже 23 эВ, что позволяет точно учитывать резонансную самоблокировку сечений.

Код ARTEMIS предназначен для расчётов активных зон. ARTEMIS - трёхмерный код (нодальный метод, учёт выгорания, возможности восстановления потвэльных полей энерговыделения). Это комплексный код, в котором, кроме нейтроники, учитываются также теплогидравлика и термомеханика.

В состав комплекса входит также ряд вспомогательных кодов, в том числе, CLARUS (подготовка файлов исходных данных для генерации библиотек констант для ARTEMIS) и HERMES (формирование библиотек констант для ARTEMIS по результатам работы APOLLO-2A). ... (блок-схему скопировать не могу)
http://www.atominfo.ru/newsd/k0729.htm

Маленькая ремарка: Для произведения расчетов по реактору, активной зоне и топливу проекта АР-1000 использовались двумерные коды, 26-групповой метод Монте-Карло, и какие-то собственные топливные коды. Цитирую по памяти, в переводе: "С 1969 года этод метод (код Вестингауза) стал стандартом в индустрии." Оригинал изложу в другой теме, тут мусорить желания нет.

[AtomInfo.Ru]

26-групповой метод...

Вопрос о числе групп не такой простой как кажется.

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

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

В СССР, например, для борьбы с этим явлением использовалась процедура "бэжитиков" - каждый расчёт выполнялся дважды, после первого расчёта делалась попытка восстановить спектр нейтронов внутри группы, переусреднялись сечения упругого рассеяния и делался окончательный расчёт.

Переход на систему, в которой порядка 300 групп, позволяет проблему упругого рассеяния снять (это я видел своими глазами).

Но! Вышесказанное важно для реакторов быстрых и с промежуточным спектром. Для тепловых реакторов проблема корректного усреднения групповых сечений упругого рассеяния практически не имеет значения. И в первом общемировом стандарте расчёта ячеек тепловых реакторов (британская программа WIMS) диапазон от 10 МэВ до 4 эВ спокойно разбивался всего на 27 групп.

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

Есть ещё одна область, где может потребоваться много групп - тепловая область энергий. Её можно определять по-всякому - например, в WIMS она начиналась от 4 эВ и шла до нуля, были варианты, когда тепловая область начиналась от 2 эВ. Общее у всех способов - тепловая область должна начинаться ниже первого резонанса U-238 (6,67 эВ, если помню верно) и доходить до нуля энергии.

Как показала практика, на эту область (её называют термализацией, и её характерная особенность - здесь нейтрон может при столкновении не только потерять, но и приобрести энергию, а матрица групповых уравнений перестаёт быть треугольной) достаточно порядка 40-50 групп. У того же WIMS - 42 группы. У советского ТЕРМАК'а было 45 групп (по памяти). Можно взять, допустим, 75 групп и это будет с огромным запасом по детализации тепловой области энергии.

Отсюда вывод. Для хороших расчётов тепловых реакторов достаточно иметь порядка 100 или менее групп по энергии (20-30 в области выше термализации и до 75 в области термализации). Расчёты же с 300 группами выглядят... м-м-м, некоторым излишним архитектурным украшательством. Вот если бы было 3000 групп, это был бы качественный скачок вверх.

И о Вестингаузе. Скорее всего, там были не просто 26-групповые расчёты, а расчёты в 26 группах плюс некоторый механизм в области термализации. Какой механизм - не знаю. Но, скажем, разработчики советской БНАБ заменяли 26-ую группу на тепловую точку и очень грамотно присваивали ей сечения. В результате, в 26 группах можно было с хорошей точностью считать ВВЭР-ы.

[asv363]

Возник вопрос о средствах и методах расчёта равновесной концентрации изотопов и их активности в различных зонах, помещениях ЭБ и вспомогательных строений. Аналогично, интересует вопрос о расчёте будущих индивидуальных и коллективных доз. Связано это с пристальным исследованием одного конкурирующего американо-японского проекта, не отработавшего ни одного дня. Формулы, в принципе, приведены в 12.2.2.4 для каждого изотопа, и в общем знакомы, однако, нет никаких упоминаний о кодах и методах расчёта транспорта. Для данного проекта, это является редкостью.

Глава 12.2 страница 12.2-7 - две формулы, далее таблицы.
AP12.2

Глава 12.4, начиная с 12.4-6 - таблицы.
AP12.4

Если внимательно посмотреть, то они постулируют некую течь, объем помещения, время для достижения равновесной концентрации в 100 дней и прочее, но чем считают, не ясно.

Сообщение отредактировал asv363 - 9.4.2013, 0:11