Table of Contents

Введение в криогенную динамику жидкости в космических системах

Криогенные жидкости — сжиженные газы, хранящиеся при температурах ниже -150 ° C — являются источником жизненной силы современных систем космического движения и жизнеобеспечения. Жидкий водород (LH2) и жидкий кислород (LOX) питают основные двигатели ракет, такие как Raptor SpaceX, BE-4 Blue Origin и система космического запуска НАСА. Управление этими жидкостями в условиях микрогравитации и экстремальных тепловых средах представляет собой уникальные инженерные проблемы. Изменение фазы, стратификация, нестабильность гейзера и откачка должны быть точно предсказаны для обеспечения безопасности миссии и топливной эффективности.

Инструменты вычислительной динамики жидкости (CFD), такие как COMSOL Multiphysics с модулем CFD , предлагают надежную основу для моделирования криогенного поведения жидкости. Благодаря соединению потока жидкости, теплопередачи, термодинамики и фазового изменения в одной среде инженеры могут тестировать конструкции резервуаров, линии подачи и стратегии изоляции до того, как будет построено оборудование. Эта статья предоставляет авторитетное, готовое к производству руководство по моделированию криогенной динамики жидкости в COMSOL CFD для космических применений.

Почему криогенные жидкости требуют специального моделирования

В отличие от обычных жидкостей, криогены демонстрируют сильные изменения свойств с температурой и давлением. Жидкий водород, например, имеет плотность всего около 70 кг/м3 в точке кипения и теплопроводность, которая резко меняется вблизи критической точки. Режимы кипения — ядро, переход и кипение пленки — могут сосуществовать в одном резервуаре из-за локализованного ввода тепла от двигателя или солнечного излучения. В микрогравитации доминирует поверхностное натяжение, создавая потоки, управляемые капиллярами, которые влияют на позиционирование жидкости и теплообмен.

Эти сложности делают эмпирические корреляции ненадежными для проектирования. Мультифизическое моделирование должно одновременно разрешать:

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

Мультифизический подход Comsol позволяет инженерам объединять эти явления в одну модель, а не повторяться между отдельными инструментами.

Основные возможности COMSOL CFD для криогенного моделирования

Управляющие уравнения и модели турбулентности

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

Изменение фазы и многофазный поток

COMSOL предлагает два основных подхода для изменения фазы:
Набор уровней и Фазовое поле Методы явно отслеживают интерфейс жидкостного пара. Они идеально подходят для задач с свободной поверхностью, таких как выплескивание резервуара или струйное сопротивление. Модель смеси лучше подходит для дисперсных двухфазных потоков (пузырь или капля). Для кипения интерфейс Изменение тепловой фазы может использовать пользовательские функции для скорости передачи массы на основе температурных градиентов.

Неидеальная термодинамика

Функция COMSOL Real Gas позволяет импортировать фундаментальные уравнения состояния, такие как Пэн-Робинсон или Спан-Вагнер, для высокоточных предсказаний свойств. Для водорода и кислорода эти уравнения необходимы вблизи критической точки, где предположения идеального газа дают ошибки, превышающие 30 процентов.

Пользовательские библиотеки материалов

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

Пошаговый процесс моделирования в COMSOL

1. Определение геометрии

Начните с 2D аксисимметричной или полной 3D-геометрии криогенного резервуара, кормовой линии или головки двигателя. Для типичного сферического топливного бака включите толщину стенки, впускные/выходные порты, противовихрные перегородки и слои изоляции. Упростите некритические филе и болтовые отверстия, чтобы поддерживать управляемость сетки.

2. Присвоение свойств материала

Создать пользовательский материал для криогена с использованием таблиц NIST. Например, для LH2 при 20 К - 30 К, плотность импорта в зависимости от температуры и давления, и установить динамическую вязкость с использованием модели Леннарда-Джонса. Для стенки резервуара (например, алюминий 2219 или нержавеющая сталь 316L), определить теплопроводность, плотность и коэффициент теплового расширения. Добавить изоляционный материал, такой как многослойная изоляция (MLI) с эффективной теплопроводностью, которая зависит от вакуумного давления.

3.Физические интерфейсы и связи

Добавьте следующие интерфейсы из Мастера Модели:

  • Ламинарный поток или Турбулентный поток (k-ω SST) для движения жидкости.
  • Передача тепла в жидкостях и Передача тепла в твердых телах для сопряженного теплопередачи.
  • Поле фазы (или набор уровней) для отслеживания интерфейса, если исследовать двухфазную динамику.
  • Неизотермический поток мультифизическая связь для учета температуры-зависимой плотности и плавучести.

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

4.Границы условий

Типичные условия для бака с топливом во время береговой фазы:

  • Наружная поверхность стенки резервуара: тепловой поток от солнечного излучения (например, 1,3 кВт/м2 при 1 а.е. плюс инфракрасный луч Земли).
  • Внутренняя поверхность стенки резервуара: нет-скользящая, конъюгированный тепловой поток.
  • Интерфейс жидкостного пара: непрерывность температуры, баланс массы с частотой изменения фазы.
  • Выходная (линия подачи): граница давления, равная давлению в корпусе танка минус головка.
  • Вход (для заполнения): граничное условие скорости потока массы.

Для моделирования микрогравитации используйте Сила объема для поверхностного натяжения (на основе модели континуумной поверхностной силы) и установите ускорение до нуля или небольшого значения микро-g (10-4 г).

5.Стратегия сетки

Высококачественная сетка имеет решающее значение для криогенного моделирования. Используйте сетку пограничного слоя (по крайней мере, 5-10 призмовых слоев) на стенке резервуара для захвата крутых температурных градиентов. В области двухфазного интерфейса уточните элементы размером около 1 мм (для резервуара диаметром 1 м) для разрешения кривизны интерфейса. Используйте сетку свободного треугольного угла для 2D аксисимметричных моделей или свободную тетраэдрическую сетку для 3D.

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

6. Настройки растворителя и шаг времени

Криогенные проблемы часто переходные. Используйте Формулу обратной дифференциации (BDF) ] решатель порядка 1 или 2 с начальным шагом времени 0,01 с. Для сопряженной физики, включите Полностью связанный решатель с коэффициентом демпфирования 0,8 для улучшения конвергенции. Если модель работает медленно, переключитесь на Сегрегированный подход, сначала решая поток, затем теплообмен, затем фазовое поле на каждом шаге времени.

Конвергенция контролируется путем проверки остатков (ниже 10-4) и консервированных количеств (баланс массы в пределах 0,1 %).

7.Постпроцессинг и анализ

Ключевые выходы для извлечения:

  • Температурные контуры на стенках резервуара и флюидной области.
  • Скорость отключения (кг/с) интегрированы по интерфейсу.
  • Эволюция давления в объеме ульяжа (важно для конструкции резервуара).
  • Стремовые линии или линии пути , показывающие зоны рециркуляции, которые могут вызывать термическое расслоение.
  • Положение интерфейса как функция времени — критическое для оценки устройств приобретения жидкости.

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

Пример: моделирование выхлопа в криогенном резервуаре для топлива

Рассмотрим сферический резервуар LH2 (радиус 1,5 м) с алюминиевыми стенками толщиной 1 мм и 5 см MLI. Танк заполнен на 50%, при давлении 3 бар, с жидкостью при температуре насыщения (около 23 К). Внешний тепловой поток составляет 150 Вт/м2 от тепловой среды космического корабля. Цель состоит в том, чтобы предсказать скорость откипела в течение 24-часового периода побережья.

После настройки модели, как описано выше, моделирование показывает, что MLI уменьшает утечку тепла примерно до 8 Вт/м2. Скорость кипения стабилизируется на 0,25 г/с после начального переходного периода в 200 с. Параметрические прострелы показывают, что удвоение толщины MLI уменьшает откипь на 45 %, но весовой штраф может быть неприемлемым. Модель также прогнозирует температурный градиент 2 К между днищем и верхом резервуара, что помогает определить места, где может произойти конденсация пара, что потенциально вызывает распад давления в резервуаре и трудности с зажиганием.

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

Приложения в жизненном цикле космического корабля

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

Криогенные линии подачи должны передавать жидкость без кавитации или двухфазного потока, который может привести к голоду двигателя. Модели COMSOL линии подачи, включая изгибы, соединения и клапаны, могут прогнозировать падения давления и обнаруживать области, где кипение может произойти из-за локальных горячих точек. Модель смеси с скоростью скольжения захватывает разделение фаз в микрогравитационных каналах, помогая инженерам проектировать устройства получения жидкости (LAD), такие как экранные каналы или лопасти.

Проверка системы термозащиты

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

Запуск Pad Chill-Down и заполнение операций

Заполнение резервуара криогеном из условий окружающей среды приводит к сильному кипению и тепловому удару. Переходные модели в COMSOL могут оценить время, необходимое для охлаждения стенок резервуара до безопасных температур (часто ниже 100 К) перед основным заполнением. Эти модели также помогают оптимизировать скорость заполнения, чтобы избежать чрезмерного образования пара, который может повредить клапаны.

Управление в полете

Слошинг в условиях микрогравитации вызывает движение жидкости, которое может дестабилизировать космический корабль. Совместные модели гидроструктуры в COMSOL позволяют инженерам проектировать перегородки, диафрагмы или резервуары поверхностного натяжения, которые удерживают жидкость. Метод FLT:0 Фазовое поле фиксирует движущийся интерфейс во время ожогов двигателя или натяжных двигателей, обеспечивая силы на стенках резервуара для структурного анализа.

Корреляции наземных испытаний

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

Внешние ресурсы для продвинутых тем

Для инженеров, стремящихся углубить свои криогенные знания о CFD, рекомендуется использовать следующие ресурсы:

Будущие направления: цифровые близнецы и машинное обучение

Следующим рубежом в криогенном моделировании жидкости является создание цифровых двойников в реальном времени. Модели с пониженным порядком (ROM) могут быть построены из параметрических проверочных приборов COMSOL, а затем встроены в программное обеспечение для полета транспортного средства для оценки состояния. Например, ROM, который прогнозирует скорость отката в зависимости от утечки тепла, уровня заполнения и ускорения, может работать в миллисекундах, что позволяет измерять бортовое топливо.

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

LiveLink от COMSOL для MATLAB позволяет легко взаимодействовать с данными моделирования с помощью таких инструментов ИИ, позволяя инженерам итеративно улучшать свои ПЗУ с помощью свежих данных высокой точности по мере продвижения миссии.

Заключение

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

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