Анализ процессов горения с помощью Ansys Fluent: советы и хитрости для получения точных результатов

Создание надежного фундамента: предварительный процесс

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

Первым шагом в любом исследовании CFD является создание точной вычислительной области. Для приложений сгорания даже небольшие геометрические детали – такие как держатели пламени, насадки инжектора или охлаждающие проходы – могут значительно влиять на структуру потока и стабильность пламени. Начните с импорта чистой модели CAD, удаления ненужных функций и упрощения некритических секций без ущерба для точности. Затем сосредоточьтесь на сетке. Всегексаэдрическая сетка или хорошо структурированная гибридная сетка часто дает превосходные результаты для потоков горения, потому что гексаэдрические элементы лучше выравниваются с направлением потока и уменьшают численную диффузию. В регионах, где находится пламя, таких как зона реакции и смежные слои сдвига, вы должны совершенствовать сетку, чтобы захватывать крутые градиенты температуры, концентрации видов и скорости. Используйте местные средства управления размерами, чтобы обеспечить общее количество ячеек по толщине пламени. Общее руководство заключается в обеспечении по меньшей мере 10 – 15 ячеек по толщине пламени; если пламя тонкое,

Избегать сильно перекошенных ячеек (перекос ниже 0,9) и высоких пропорций сторон в интересующих областях. метрика ортогонального качества в Fluent является надежным показателем; значения выше 0,15 в целом приемлемы, но для сгорания следует стремиться к 0,2 или выше вблизи зон пламени. Используйте буферизованный сетчатый подход: грубо огрубляйте выше и ниже по течению от зоны реакции постепенно, чтобы снизить вычислительные затраты без ущерба для точности. Для переходных симуляций с движущимися пламенями (например, в поршневых двигателях или пульсирующих горелках) рассмотрите динамическую адаптацию сетки на основе температуры или градиентов скорости реакции.

Граничные условия и инициализация

Реалистичные граничные условия являются источником жизненной силы надежного моделирования. Для впускных отверстий указывают общую температуру, скорость или массовый расход и интенсивность турбулентности (обычно 5–10% для воздуха сгорания). Обратите пристальное внимание на фракции массы видов: если используется упрощенный топливный отвод (например, метан), установите топливный отвод с соответствующим составом. Для воздушных отверстий используют 23,3% O2 по массе (или 21% по объему) и 76,7% N2. Если система включает зоны рециркуляции (например, в печи), выход давления с заданной температурой обратного потока и видовой состав помогает стабилизировать конвергенцию. Настенные граничные условия требуют тщательной термической обработки – изотермический, тепловой поток или в сочетании с конъюгированной теплообменной передачей – в зависимости от физической проблемы. Для изолированных стен используют условие нулевого теплового потока; для сжигателей с водяным охлаждением более целесообразным является постоянный температурный или коэффициент теплопередачи.

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

Выбор правильной модели турбулентности

Турбулентность и горение тесно связаны; смешивание топлива и окислителя, морщины пламени и выделение тепла зависят от турбулентных вихрей. Ansys Fluent предоставляет ряд моделей турбулентности на основе RANS. Ниже приведено практическое руководство по выбору среди наиболее распространенных вариантов применения горения.

RNG k-ε и реализуемый k-ε

Обе модели являются промышленными рабочими лошадками. Модель RNG k-ε включает дополнительный термин в уравнении рассеивания, который повышает точность для закручивания потоков и умеренной рециркуляции. Модель Realizable k-ε удовлетворяет математическим ограничениям на нормальные напряжения и лучше работает для плоских и круглых струй, которые являются общими в конструкциях горелок. Для несмешанных пламен с умеренным закручиванием (число вихрей <0,6) любая модель может дать удовлетворительные средние температуры и профили видов. Если поток включает сильную обтекаемую кривизну (например, в сливных горючих материалах), часто предпочтительнее RNG k-ε.

k-ω SST (транспорт сдвига)

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

Стрессовая модель Рейнольдса (RSM)

RSM решает транспортные уравнения для каждого компонента напряжения Рейнольдса, что делает его наиболее физически полной моделью RANS. Это вычислительно дорого (5–10 раз медленнее, чем модели с двумя уравнениями), но может быть необходимо для сильно закрученных потоков (число вихрей > 0,7), сильно анизотропной турбулентности или потоков со сложной кривизной. Для большинства проблем промышленного сгорания RSM редко оправдан, если данные проверки не показывают, что более простые модели терпят неудачу.

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

Химия горения: балансировка деталей и вычислительные затраты

Химическая кинетика регулирует скорость пламени, задержку зажигания, образование загрязняющих веществ (NOx, CO, сажа) и пределы исчезновения. Ansys Fluent предлагает несколько подходов к модели химии горения, каждый из которых имеет компромиссы.

Концепция диссипатии Эдди (EDC)

EDC — это конечная модель химии, которая предполагает, что реакции происходят в тонких турбулентных структурах, где молекулярное смешивание интенсивно. Он подходит для турбулентного премиксированного и непремиксированного пламени, где реакция быстрая по отношению к смешиванию. EDC может обрабатывать подробные химические механизмы (например, GRI 3.0 для метана) без необходимости ламинарных предположений о фламлетах. Однако для практических симуляций многие инженеры используют сокращенные механизмы (например, 15–30 видов), которые точно воспроизводят скорость пламени, воспламенение и поведение вымирания. Хорошей отправной точкой является механизм Петерсена-Макгуайра для метана или модель CRECK для более тяжелых углеводородов. Валидируйте выбранный вами механизм против экспериментальных данных о скорости пламени ламинар, доступных из источников, таких как NIST Chemical Kinetics Database .

Для непремиксированного пламени, модель устойчивый ламинарный фламлет (SLFM) является высокоэффективной. предполагает, что пламя тонкое и химические временные шкалы намного короче, чем турбулентные шкалы смешивания. Предварительно составленная библиотека скаляров (температура, фракции массы вида) в качестве функции фракции смеси и скорости рассеивания скаляров генерируется с использованием 1D-решителя пламени. Затем флюент просматривает эти значения во время моделирования. Этот подход на порядки быстрее, чем EDC и хорошо работает для диффузионного пламени без частичного вымирания. Однако он борется с поднятым пламенем, задержкой воспламенения и случаями, когда потери тепла или излучение являются доминирующими.

Частично смешанное горение

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

Запуск симуляции: настройки Solver и конвергенция

Сцепление скоростей давления

Для моделирования сгорания в устойчивом состоянии используйте растворитель Fluent, соединенный с FLT:0. Он обеспечивает более быструю конвергенцию для сильно связанных потоков (например, реагирующих потоков с большим выделением тепла) по сравнению с сегрегированными растворителями. Алгоритм FLT:2 SIMPLE работает для простых геометрий, но может бороться с изменениями рециркуляции и плотности в пламени. Сначала устанавливайте факторы пониженного расслабления для плотности и температуры до низких значений (0,3 –0,5), чтобы предотвратить расхождение, а затем постепенно увеличивайтесь по мере стабилизации раствора.

Схемы дискретизации

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

Мониторинг конвергенции

Не полагайтесь исключительно на масштабированные остатки (например, 1e-5 для непрерывности). Моделирование горения часто демонстрирует осцилляционные остатки из-за динамики пламени. Вместо этого, контролируйте интегральные величины: температуру выхода, дисбаланс массового потока и общее выделение тепла. Когда они остаются стабильными в течение сотен итераций, решение, вероятно, сходится. Кроме того, проверьте, что средняя температура в ключевых плоскостях не дрейфует. Для переходного моделирования используйте этап времени, который захватывает самую быструю химическую шкалу времени (обычно от 1e-6 до 1e-5s для углеводородного пламени) и работает в течение 3–5 времени прохождения потока.