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

Понимание PWR и необходимость моделирования

Реактор с водой под давлением является наиболее распространенным типом реактора с ядерной энергией во всем мире, работающим в более чем 300 реакторах по всему миру. В PWR вода поддерживается под высоким давлением (около 15,5 МПа или 2250 пси) внутри первичной петли охлаждающей жидкости, предотвращая ее кипение даже при температурах, превышающих 300 ° C. Эта вода под давлением поглощает тепло от топливных сборок в ядре реактора и транспортирует его в парогенератор, где она передает тепло во вторичную петлю, производя пар для привода турбины. Первичная петля содержит реакторный сосуд, насосы охлаждающей жидкости, герметизатор и трубчатый пучок парогенератора. Вторичная петля физически отделена, гарантируя, что радиоактивные материалы в первичной петле не достигают турбины или окружающей среды.

Сложность систем PWR вводит многочисленные потенциальные сценарии аварий — как аварий на проектной основе (DBA), так и аварий вне проектной основы (BDBA), которые должны быть поняты для обеспечения безопасной работы. Такие аварии, как авария с потерей охлаждающей жидкости (LOCA), авария с включением реактивности (RIA) или отключение станции, включают в себя каскад взаимозависимых физических процессов: поток жидкости через сложные геометрии, теплообмен от топлива до охлаждающей жидкости, нейтронная кинетика в ядре и поведение материала при экстремальных температурах и давлении. Экспериментальное воссоздание этих условий в полном масштабе является непомерно дорогим и часто невозможным из-за ограничений безопасности. Вычислительные модели заполняют этот пробел, позволяя аналитикам прогнозировать прогрессирование несчастных случаев, оценивать производительность инженерных функций безопасности (например, системы аварийного охлаждения ядра, системы сдерживания распыления) и информировать регулирующие решения.

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

Типы используемых вычислительных моделей

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

Тепло-гидравлические системные коды

Коды теплогидравлической системы являются рабочими лошадками анализа аварий PWR. Они имитируют общую систему охлаждающей жидкости реактора (RCS) и ее реакцию на переходы и аварии с использованием одномерного или квазидвумерного представления путей потока (например, трубы, насосы, сосуды). Коды, такие как RELAP5 (разработанные Комиссией по ядерному регулированию США), TRACE и ATHLET , решают управляющие уравнения для массы, импульса и энергосбережения для двухфазной жидкости (], включают модели критического потока (захлеб), режимы теплопередачи (от кипения с субхолодным ядром до кипения пленки), нагрева ядра и поведение систем пассивной безопасности, таких как аккумуляторы и аварийные охлаждающие баки ядра. Системные коды относительно быстры в вычислительном отношении, что делает их пригодными для

Вычислительная динамика жидкости (CFD)

В то время как системные коды обеспечивают представление на системном уровне, модели CFD захватывают трехмерные, переходные потоки и локальные эффекты, которые имеют решающее значение для определенных аварий. Инструменты CDF, такие как ANSYS Fluent , STAR-CCM + и OpenFOAM, решают уравнения Навье-Стокса с моделями турбулентности (например, k-ε, k-ω SST) и многофазными моделями (например, Eulerian-Eulerian, Volume of Fluid).

  • Термическое расслоение в системе охлаждения реактора во время естественной циркуляции
  • Смешивание аварийной охлаждающей воды с первичным охлаждающим веществом
  • Кипение бассейна и критический тепловой поток (CHF) на поверхностях топливных стержней
  • Распределение водорода в контейнере после тяжелой аварии
  • Эффективность распылительного сопла в охлаждении в условиях содержания

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

Нейтронные транспортные симуляции

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

  • Детерминированные коды (например, PARCS, SIMULATE-3): Решите уравнение диффузии нейтронов или уравнение транспортировки в стационарном и переходном режимах.
  • Коды Монте-Карло (например, MCNP, Змеиный ): Используйте случайную выборку для моделирования отдельных историй нейтронов, обеспечивая очень точное распределение потоков и скорости реакции.В то время как вычислительно интенсивные методы Монте-Карло необходимы для бенчмаркинга детерминированных кодов и для анализа сложных геометрий, таких как топливные сборки с горючими ядовитыми стержнями.

Связанные теплогидравлические/нейтронные коды позволяют аналитикам моделировать обратную связь с реакционной способностью во время аварий (например, расширение Доплера в топливе, охлаждение, обратная связь с модератором температуры). Такая связь имеет решающее значение для анализа РИА и ожидаемого переходного без скрема (ATWS).

Структурная целостность и химические модели

Прогрессирование аварии может налагать экстремальное механическое и химическое напряжение на компоненты реактора. Коды анализа конечных элементов (FEA), такие как ANSYS Mechanical и ABAQUS , оценивают предельные значения напряжения, деформации и отказа для сосуда давления реактора, трубопроводов и сдерживания под тепловыми и давлением нагрузками. Модели химических реакций, часто встроенные в коды тяжелых аварий, такие как MAAP или MELCOR , предсказывают такие явления, как:

  • Окисление Циркалоя и образование водорода
  • Взаимодействие с топливной облицовкой
  • Выпуск и транспортирование продуктов деления (йод, цезий и т. Д.)
  • Расплавленное ядро-конкретное взаимодействие во время сценариев расплава ядра

Применение в анализе сценариев несчастных случаев

Вычислительные модели используются в различных сценариях аварий для информирования о проектировании, лицензировании и готовности к чрезвычайным ситуациям. Стандарты безопасности МАГАТЭ требуют, чтобы детерминированный анализ безопасности проводился с использованием проверенных кодов для всех аварий на проектной основе. Ниже приведены основные категории сценариев и используемые конкретные подходы к моделированию.

Аварии с потерей охлаждающей жидкости (LOCA)

LOCA предполагает разрыв трубы в системе первичной охлаждающей жидкости, приводящий к разгерметизации и потере запасов охлаждающей жидкости. Сценарий далее подразделяется на крупноразрывную LOCA (LBLOCA, область разлома > 0,093 м2) и мелкоразрывную LOCA (SBLOCA). Для LBLOCA коды теплогидравлической системы имитируют фазу выдувания, пополнения и повторного затопления ядра системами аварийного охлаждения ядра. Ключевые явления включают критический двухфазный поток через разрыв, вызванное разгерметизацией вакуумирование и распространение затухания передней части, когда вода поднимается со дна, чтобы повторно намокнуть топливные стержни. Надежные модели CCFL (ограничение потока тока) необходимы для прогнозирования того, может ли вода достичь ядра. CDF часто используется для анализа эффективности впрыска прижимного устройства. Для SBLOCA доминирует естественная циркуляция и петлевое поведение, требующее системных кодов с подробной нодаляцией первичных и вторичных петель.

Реактивность инициированных аварий (РИА)

События РИА, такие как выброс управляющего стержня или упавший стержень, вызывают внезапное включение реактивности в ядро, что приводит к быстрому отводу энергии. Температура топлива может резко повышаться, вызывая механическое взаимодействие в гранулах (PCMI) и, в тяжелых случаях, фрагментацию и рассеивание топлива. Эти аварии обычно моделируются путем соединения кода нейтронной кинетики с тепловой гидравлической моделью, которая учитывает обратную связь температуры топлива и охлаждающей жидкости. Высокоточность, тепловая гидравлическая система подканального уровня может использоваться для отслеживания местных горячих каналов и обеспечения того, чтобы пиковая энтальпия топлива оставалась ниже нормативных пределов. OECD / NEA RIA бенчмарк сыграла важную роль в проверке этих связанных кодов.

Отключение станции и тяжелые аварии

Отключение станции (SBO) связано с потерей всей мощности переменного тока, что приводит к отказу активных систем охлаждения. Завод полагается на пассивные системы (или аварийные дизель-генераторы) и, в конечном итоге, система охлаждающей жидкости реактора нагревается и оказывает давление. Если охлаждение не восстанавливается, ядро может раскрывать, перегреваться и плавиться. Тяжелые коды аварий, такие как MELCOR или ]MAAP используются для моделирования всего перехода от разгерметизации системы к повреждению ядра, отказу судна и реакции сдерживания. Эти коды включают модели поведения кория, высвобождения продукта деления и термогидравлических систем сдерживания. Результаты определяют разработку руководящих принципов управления серьезными авариями (SAMGs) и стратегий смягчения последствий аварии после Фукусимы.

Другие важные сценарии

  • Предусмотренный переход без скрема (ATWS) : Неисправность управляющих стержней для вставки во время ожидаемого события (например, поездка турбины, потеря подачей воды). Связанные нейтронные / термально-гидравлические коды имитируют вставку реактивности из разбавления бора или охлаждения жидкости и оценивают, может ли растение выжить без скрема через внутреннюю обратную связь.
  • Основной разрыв паровой линии (MSLB) : разрыв во вторичной стороне, который увеличивает нагрузку на парогенератор, потенциально переохлаждая первичную систему и вставляя положительную реактивность. Для обеспечения того, чтобы ядро оставалось подкритическим, необходимы подробные системные коды и переходный нейтронный анализ.
  • Реакция на содержание и безопасность водорода: После тяжелой аварии водород от окисления циркония может накапливаться и представлять риск дефлагации или детонации. Модели CFD смешивания атмосферы сдерживания в сочетании с моделированием химической реакции используются для проектирования систем смягчения воздействия водорода (например, воспламенители, пассивные автокаталитические рекомбинеры).

Преимущества и ограничения вычислительного моделирования

Преимущества

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

Ограничения и вызовы

Несмотря на свою мощь, вычислительные модели имеют неотъемлемые ограничения. Основная проблема заключается в неопределенной количественной оценке . Параметры ввода — такие как начальные условия, граничные условия и свойства материала — неопределенны, а структуры моделей упрощают реальность. Пользователи кода должны применять методы наилучшей оценки плюс неопределенность (BEPU), как это предусмотрено правилами, такими как USNRC 10 CFR 50.46 для LOCA. Кроме того, мультифизическая связь страдает от численных ошибок и проблем синхронизации по времени. Вычислительная сетка должна захватывать масштабы от миллиметров (поверхности топливных стержней) до метров (сдерживание), что делает моделирование высокого разрешения чрезвычайно дорогим. Валидация остается критической: модель только так хороша, как экспериментальные данные, против которых она оценивается. Пробелы существуют для таких явлений, как распространение расплава или химия продуктов деления в тяжелых авариях, где эксперименты скудны. Следовательно, одобрение регулирующих органов часто требует «глубинной защиты», которая сочетает в себе

Будущие направления

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

  • Высокоточная мультифизическая связь : Усилия, такие как CASL (Консорциум по усовершенствованному моделированию легких водных реакторов) Министерства энергетики США, направлены на то, чтобы связать нейтронику, тепловую гидравлику и характеристики топлива с беспрецедентным разрешением, используя 3-D CFD и нейтронику Монте-Карло на экзафлопсных компьютерах. Эти подробные модели уменьшат консерватизм в запасах безопасности.
  • Дирижируемое данными и AI-ассистируемое моделирование: Машинное обучение применяется для ускорения распространения неопределенности, разработки суррогатных моделей для моделирования аварий и улучшения мониторинга состояния растений. Например, нейронные сети могут быстро предсказать эволюцию ключевых параметров безопасности во время сценария на основе предварительно вычисленной базы данных моделирования.
  • Цифровые двойники: цифровые двойники PWR в реальном времени, сочетающие данные датчиков с живым моделированием, могут обеспечить раннее предупреждение об аномальных условиях и поддержку принятия оператором решений во время аварий.
  • Улучшенное моделирование тяжелых аварий : Улучшенные модели поведения кория, расплавленного топлива-конкретные взаимодействия и аэрозольной физики разрабатываются для поддержки следующего поколения анализов безопасности, особенно для растений поколения III+.

По мере того, как вычислительная мощность продолжает расти, а данные экспериментальных объектов (например, тестов ОЭСР/НЕА ]ATLAS , LOBI и PKL ) становятся все более доступными, вычислительное моделирование будет играть все более центральную роль в обеспечении работы PWR с самыми высокими стандартами безопасности и надежности.

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