Table of Contents

Введение в моделирование мороза в COMSOL Multiphysics

Формирование и таяние мороза — сложные мультифизические явления, играющие важнейшую роль во многих промышленных и природных системах. В аэрокосмической технике мороз на поверхностях самолётов может нарушать воздушный поток и снижать подъём. В холодильной технике наращивание мороза на катушках испарителя ухудшает эффективность теплопередачи и увеличивает потребление энергии. В климатологии накопление мороза на растительности и почве влияет на местные энергетические и влагобалансы. Точное моделирование этих процессов требует сопряженной обработки теплопередачи, массового транспорта и фазового изменения. COMSOL Multiphysics с её гибким модулем CFD и встроенной мультифизической связью обеспечивает среду, в которой инженеры и исследователи могут моделировать динамику мороза с высокой точностью. Эта статья расширяется на фундаментальную настройку, численные методы, подходы валидации и практические применения моделирования образования мороза и плавления в COMSOL.

Физические принципы образования и таяния мороза

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

Уравнения для управляющей фазы жидкости

В области воздушного потока уравнения Рейнольдса-среднего Навье-Стокса (RANS) обычно решаются для моделирования турбулентного потока вокруг подверженных морозу поверхностей. Уравнение непрерывности обеспечивает сохранение массы, в то время как уравнения импульса включают закрытие турбулентности, часто используя модели k-ε или k-ω SST. Уравнение энергии включает в себя конвективный теплообмен и, в сочетании с переносом влаги, эффект скрытого теплоотдачи или поглощения. Для поля влаги решается уравнение конвекции-диффузии для фракции массы водяного пара, с термином источника или поглотителя, представляющим сублимацию или осаждение. Граничное условие на поверхности обморожения соединяет поток пара с тепловым балансом через отношение Клаузиуса-Клапейрона, которое определяет равновесное давление пара над льдом. Эти нелинейные связи требуют надежного численного подхода; мультифронтальный растворитель COMSOL или итеративный растворитель GMRES с предварительной кондиционацией обычно используется

Модель изменения фазы для мороза

Изменение фазы от пара к льду (осаждение) и от льда к жидкости (плавление) моделируется с использованием метода энтальпийной пористости или эквивалентной теплоемкости. В COMSOL изменение фазы может быть реализовано с помощью функции сглаженного шага, которая определяет фракцию массы льда как функцию температуры и концентрации пара. Скрытое тепло включено в качестве объемного исходного термина в энергетическом уравнении, пропорциональное временной производной фракции массы льда. Для роста мороза подмодель отслеживает пористость и теплопроводность слоя мороза, оба из которых развиваются с плотностью и толщиной. Эмпирические корреляции для теплопроводности мороза (например, из работы Хаяси или О'Нила и Древа) часто вводятся в качестве функций плотности, которая сама зависит от скорости осаждения и времени. Эти подмодели имеют решающее значение, потому что плотный слой мороза с низкой пористостью проводит тепло иначе, чем пушистый, высокопористый.

Настройка модели COMSOL шаг за шагом

Создание геометрии и Mesh-расчеты

Начните с определения геометрии во встроенной среде САПР COMSOL или импорта из внешнего инструмента. Типичная модель включает плоскую пластину, представляющую поверхность плавника или крыла теплообменника, помещенную в канал потока, который определяет воздушную область. Для аксисимметричной геометрии (например, охлаждающей трубки) 2D аксисимметричная установка может снизить вычислительную стоимость. Сетка должна быть очищена вблизи поверхности обморожения для разрешения резких градиентов температуры, концентрации пара и скорости. Сетка пограничного слоя с 5-10 слоями призмы и первой толщиной ячейки, такая, что y + ≈ 1 рекомендуется для моделей турбулентности, которые разрешают вязкий подслой. Исследования независимости сетки должны проводиться путем сравнения результатов толщины мороза и коэффициента теплопередачи на более грубых, средних и тонких сетках.

Определение физических интерфейсов

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

  • Передача тепла в жидкостях (ht) — для распределения температуры в воздухе и твердом веществе.
  • Перенос тепла в твердых телах (ht) — для стенки или подложки.
  • Транспорт разбавленных видов (tds) — для фракции массы водяного пара в воздухе.
  • Турбулентный поток, k-ε или k-ω (spf) — для движения движения.
  • Мультифизические связи — неизотермическая связь потока (псевдосвязь) и пользовательская связь «Отложение мороза», которая связывает поток пара с тепловым балансом через пограничное условие.

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

Свойства материала

Точные данные материала необходимы. Для воздуха плотность следует закону идеального газа, теплопроводность и вязкость зависят от температуры (используются встроенные выражения). Диффузивность водяного пара в воздухе может быть взята из эмпирических формул (например, из Cussler или Perry’s Handbook). Для слоя мороза обычно составляет от 50 до 300 кг/м3 в зависимости от условий осаждения. Для выражения, такого как ρ frost = 650 − 2.0×106/(T frost + 273.15) или корреляция от Hayashi (1977) часто задается k frost = 0,024 + 7,52×10−4 ρ frost + 4,1×10−9 ρ frost3. Скрытое тепло сублимации составляет примерно 2,83×106 J/кг, а для плавления составляет 3,34×105 J/кг.

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

Для стандартного моделирования обморожения:

  • Впуск: постоянная скорость (1-5 м/с, типичная для охлаждения), температура (например, от -5 до 10 °C) и относительная влажность (50-90%).
  • Пустошь: давление нулевого калибра или состояние оттока.
  • Поверхность для обморожения: стенка с предписанной температурой (например, -15°C для холодной пластины) или соединенная с энергетическим балансом, который включает конвекцию и скрытое тепло. Концентрация пара на стене устанавливается на значение насыщения льда при этой температуре стенки (через функцию или таблицу поиска).
  • Начальные условия: Весь домен потока установлен на температуру и влажность входа, нулевую фракцию массы льда. Моделирование зависит от времени; типичные сроки выполнения охватывают от минут до часов в физическом времени, с временными шагами порядка 1-10 секунд.

Моделирование образования мороза: численный рабочий процесс

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

  • Толщина мороза как функция времени (пост-обработка движущегося интерфейса или поля фракции массы льда).
  • Коэффициент теплопередачи поверхности (вычисляется из общего теплового потока через стену, деленного на разницу температур).
  • Общая масса мороза на единицу площади.
  • Профили температуры и скорости вблизи стены.

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

Моделирование циклов таяния мороза и размораживания

Моделирование плавления требует изменения обратной фазы: слой льда переходит в жидкую воду, которая затем течет или испаряется. В COMSOL один подход заключается в том, чтобы включить модель плавления тонкого слоя, которая активируется, когда температура поверхности превышает 0°C. Для многих инженерных применений может быть отслежена с помощью метода заданного уровня или интерфейса, но для многих инженерных применений достаточно упрощенного метода теплобаланса. Введите тепловой поток или повышенное температурное граничное условие на стене (например, 30 °C для электрического обогревателя мороза) и решите для времени, необходимого для полного расплавления заданной толщины мороза. Моделирование должно учитывать энергию, потребляемую скрытым теплом плавления, а также разумный нагрев пористого мороза и подложки. Дополнительная физика может включать поток водяной пленки, приводимый в действие гравитацией и сдвиг от воздушного потока, который может быть смоделирован с интерфейсом Thin-Film Flow , соединенный с передачей тепла

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

Валидация и лучшие практики

Валидация моделей мороза по экспериментальным данным имеет решающее значение. Опубликованные ориентиры включают работу Лю и Якоби (2013) по росту мороза на плоских пластинах в принудительной конвекции и исследования Падхманабана и Шерифа (2000) по морозным свойствам. Как правило, проводятся сравнения для:

  • Толщина мороза vs время
  • Температура поверхности мороза
  • Общее снижение коэффициента теплопередачи из-за мороза

Пользователи должны откалибровать эмпирические корреляции для теплопроводности и плотности мороза до их конкретных диапазонов температуры и влажности. Анализ чувствительности по этим параметрам помогает количественно оценить неопределенность. Исследования сетчатой уточнения должны быть нацелены на неопределенность толщины мороза <5%. Additionally, for turbulent flows, ensure that the y+ value at the first cell is checked—the built‑in wall functions in the k-ε model may require the first cell to lie in the log‑law region (y+ > 30 ), тогда как модель с низкой турбулентностью (например, k-ω) требует y + ~ 1. Несоответствие может привести к ошибочным прогнозам тепло- и переноса массы.

Расширенные соображения: турбулентность, сцепление и вычислительный спрос

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

Полностью сопряженное решение переноса потока, тепла и влаги с движущимся фронтом мороза может быть вычислительно тяжелым. Для длительных физических раз (например, часов роста мороза) рекомендуется адаптивное времяшагово с допуском 0,001. Пользователи могут воспользоваться эффективными параметрическими разметками COMSOL для оценки нескольких условий входа. Для больших 3D-моделей грубая сетка на незамерзающих областях и тонкая сетка только вблизи поверхности может уменьшить количество ячеек. Общее время моделирования для 2D-модели (несколько тысяч элементов) с 30 минутами физического времени может занять 1-4 часа на рабочей станции, в то время как 3D-модель плавникового и трубчатого теплообменника может занять 24 часа и более.

Приложения и промышленная значимость

Аэрокосмические противообледенительные системы

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

Производительность холодильных и тепловых насосов

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

Энергоэффективные строительные системы

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

Защита от льда ветровых турбин

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

Проблемы и новые направления исследований

Несмотря на мощь COMSOL, моделирование мороза остается сложным из-за сложности микромасштабного зарождения и стохастической природы роста кристаллов мороза. Многие современные модели полагаются на эмпирические корреляции, которые проверяются только в узких диапазонах параметров. Будущие исследования движутся к прямому численному моделированию (DNS) взаимодействий замораживания капель и образования кристаллической решетки, хотя они являются вычислительно дорогостоящими для макромасштабных проблем. Еще одной развивающейся тенденцией является использование машинного обучения для суррогатной части модели мороза — например, замена корреляции пористой теплопроводности нейронной сетью, обученной на экспериментальных данных. Способность COMSOL интегрироваться с продуктами MATLAB и LiveLink облегчает такие гибридные подходы.

Заключение

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