Моделирование теплопередачи в ядерных топливных стручках с помощью Cfd в Ansys Fluent

Введение в CFD в ядерной инженерии

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

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

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

Уравнения для передачи тепла и потока жидкости

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

В ANSYS Fluent эти уравнения дискретизируются с использованием метода конечного объема. Инженеры должны выбирать соответствующие модели турбулентности и схемы решения для баланса точности и вычислительной стоимости. Для ядерных применений реализуемая модель k-ε или модель SST (SST:2) для точного захвата переноса тепла вблизи стенки. Более продвинутые методы, такие как моделирование большого Эдди (LES), иногда используются для исследований, но остаются слишком дорогими для рутинных проектных исследований.

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

Настройка моделирования в ANSYS Fluent

Геометрия и сетчатая генерация

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

Сетка имеет решающее значение для точных предсказаний теплопередачи. Структурированная гексахедральная сетка предпочтительна для цилиндрической геометрии, поскольку она выравнивается с направлением потока и уменьшает численную диффузию. Сетка ANSYS или CFD ICEM может генерировать такие сетки с уточнениями пограничного слоя вблизи облицовочной стены для захвата крутых температурных градиентов в вязком подслое. Хорошая практика заключается в том, чтобы обеспечить безразмерное расстояние стенки (y +) около 1 для моделей турбулентности, которые разрешают вязкий подслой. Для сопряженной теплопередачи сетка должна быть непрерывной на жидкостно-твердом интерфейсе, или для обеспечения энергосбережения необходимо конформное соответствие.

Свойства материалов и граничные условия

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

Типичные граничные условия включают:

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

Физика моделей и настройки Solver

ANSYS Fluent предлагает широкий спектр физических моделей. Для однофазного охлаждающего агента (большинство PWR и BWR нормальной работы) характерны следующие варианты:

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

Конъюгатная теплопередача в топливных стержнях

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

Существует два подхода к CHT в Fluent:

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

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

Моделирование CHT показывает температурные градиенты по всей грануле, облицовке и в охлаждающую жидкость. Максимальная температура на центральной линии топлива является критическим параметром безопасности, который должен оставаться ниже точки плавления (например, ~ 2850 ° C для UO2). Высокоточная CHT помогает определить, как изменения скорости потока охлаждающей жидкости, уровня мощности или толщины облицовки влияют на этот запас.

Соображения для моделирования турбулентности

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

Для реакторов с водяным охлаждением число Рейнольдса в подканалах обычно находится в диапазоне от 104 до 105, что указывает на полностью турбулентный поток. Стандартная модель k-ε с усиленной обработкой стенки является общим выбором, поскольку она надежна и вычислительно эффективна. Однако она может недооценивать теплообмен в потоках с сильной кривизной или разделением. Модель SST k-ω сочетает в себе формулу k-ω в области ближнего стенка с k-ε в дальнем поле, часто давая лучшие результаты для теплообмена во внутренних потоках. Более продвинутые модели, такие как Reynolds Stress Model (RSM) , редко используются для моделирования промышленных топливных стержней из-за их высокой стоимости и численной жесткости.

Важно проверить выбор модели турбулентности на основе экспериментальных данных (например, эксперименты по смешиванию субканальных каналов PWR или корреляции теплопередачи стержневого пучка). Многие руководства по наилучшей практике из Комиссии по ядерному регулированию США рекомендуют использовать модель SST k-ω для CFD стержневого пучка.

Анализ результатов и оптимизация дизайна

Распределение температуры и определение точки доступа

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

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

Паттерны охлаждающего потока и эффективность теплового удаления

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

Коэффициенты теплопередачи (HTC) рассчитываются на основе моделирования путем деления потока тепла стенки на разницу между температурой стенки и местной температурой объемной жидкости. Эти HTC сравниваются с корреляциями (например, Dittus-Boelter, Gnielinski) для проверки модели. Расхождения могут указывать на необходимость усовершенствованного моделирования сетки или лучшего моделирования турбулентности.

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

Преимущества использования CFD в конструкции ядерного реактора

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

Практические проблемы и решения

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

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

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

Будущие тенденции в области CFD для ядерной инженерии

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

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

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

Заключение

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

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