Неравновесное испарение жидкостей Леннарда-Джонса: теория Энског-Власова и модель Герца-Кнудсена

Основная информация

• ID статьи: 2508.15785
• Название: Non-equilibrium evaporation of Lennard-Jones fluids: Enskog-Vlasov theory and Hertz-Knudsen model
• Авторы: Shaokang Li, Livio Gibelli, Yonghao Zhang
• Классификация: cond-mat.soft physics.flu-dyn
• Дата публикации: 14 октября 2025 г. (arXiv v2)
• Ссылка на статью: https://arxiv.org/abs/2508.15785

Аннотация

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

Исследовательский контекст и мотивация

1. Исследуемая проблема

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

2. Важность проблемы

Испарение в нанометровом масштабе играет ключевую роль в нескольких передовых технологических областях:

• Повышение эффективности разделения веществ в наноразмерных пористых мембранах
• Эффективное рассеивание тепла в высокопроизводительных электронных устройствах
• Практические приложения в проточном охлаждении

3. Ограничения существующих методов

• Соотношение Герца-Кнудсена: предполагает локальное равновесие в газовой фазе, игнорируя кинетическую структуру в слое Кнудсена
• Уравнение Энског-Власова: использует потенциал Сазерленда, который значительно отличается от реального потенциала Леннарда-Джонса
• Метод молекулярной динамики: имеет чрезмерно высокую вычислительную стоимость, что затрудняет применение в инженерной практике

4. Исследовательская мотивация

Разработка единой схемы, способной:

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

Основные вклады

1. Предложена улучшенная кинетическая модель Энског-Власова: специально разработана для жидкостей Леннарда-Джонса путём коррекции корреляционной функции для точного воспроизведения уравнения состояния реальной жидкости
2. Установлен точный метод калибровки термодинамических параметров: путём одновременного согласования параметров критической точки и тройной точки обеспечивается согласованность модели со свойствами реальной жидкости
3. Подтверждена способность модели предсказывать равновесные свойства: достигнуто высокое соответствие с экспериментальными данными по кривой жидко-газового сосуществования, коэффициентам переноса, давлению насыщенного пара и коэффициенту поверхностного натяжения
4. Выявлены характеристики распределения скоростей при неравновесном испарении: подтверждено отклонение распределения скоростей в газовой фазе от распределения Максвелла, что обеспечивает теоретическое обоснование ограничений соотношения Герца-Кнудсена

Подробное описание методов

Определение задачи

Исследование процесса испарения жидкостей Леннарда-Джонса при неравновесных условиях с акцентом на анализ характеристик распределения молекулярных скоростей вблизи жидко-газового раздела и оценку применимости классического соотношения Герца-Кнудсена.

Архитектура модели

1. Базовое кинетическое уравнение

Модель основана на модифицированном уравнении Энског-Власова:

$\frac{\partial f}{\partial t} + \boldsymbol{\xi} \cdot \frac{\partial f}{\partial \boldsymbol{x}} + \frac{\boldsymbol{F}}{m} \cdot \frac{\partial f}{\partial \boldsymbol{\xi}} = \Omega$

где:

• $f(\boldsymbol{x}, \boldsymbol{\xi}, t)$: функция распределения молекулярных скоростей
• $\boldsymbol{F}$: самосогласованное силовое поле, создаваемое потенциалом Сазерленда
• $\Omega$: упрощённый столкновительный член

2. Упрощённый столкновительный член

Столкновительный член использует разложение Тейлора: $\Omega = J_S + J_e$

где $J_S$ — член модели Шахова, $J_e$ — избыточный член, обеспечивающий правильное число Прандтля и объёмную вязкость.

3. Коррекция уравнения состояния

Используется обобщённая форма Ван-дер-Ваальса: $p = nk_BT(1 + \rho b\chi) - an^2$

Корреляционная функция $\chi$ корректируется для согласования с термодинамическими свойствами жидкости Леннарда-Джонса.

Технические инновации

1. Эмпирическая коррекция корреляционной функции

Введена эмпирическая функция $\Phi$: $\Phi = \rho b + a_1(\rho b)^2 + a_2(\rho b)^3 + a_3(\rho b)^4 + a_4(\rho b)^5$

где $\chi = \Phi' = d\Phi/d(\rho b)$

2. Калибровка с множественными ограничениями

Одновременное удовлетворение условиям:

• Критической точки: $\left(\frac{\partial p}{\partial n}\right)_{T=T_c,n=n_c} = 0$, $\left(\frac{\partial^2 p}{\partial n^2}\right)_{T=T_c,n=n_c} = 0$
• Тройной точки: равенство давления и химического потенциала

3. Коррекция коэффициентов переноса

Учёт эффектов плотности на вязкость и теплопроводность: $\mu = \frac{\mu^*}{\chi}\left(1 + \frac{2}{5}\rho b\chi\right)^2 + \frac{3}{5}\varpi$$\kappa = \frac{\kappa^*}{\chi}\left(1 + \frac{3}{5}\rho b\chi\right)^2 + c_v\varpi$

Экспериментальная установка

Набор данных

• Эталонная жидкость: аргон как типичный представитель жидкостей Леннарда-Джонса
• Экспериментальные данные: измеренные значения кривой жидко-газового сосуществования, коэффициентов переноса, давления насыщенного пара и коэффициента поверхностного натяжения
• Данные моделирования МД: для сравнительной проверки

Показатели оценки

1. Равновесные свойства:
   • Плотность жидко-газового сосуществования
   • Сдвиговая вязкость и теплопроводность
   • Давление насыщенного пара
   • Коэффициент поверхностного натяжения
2. Неравновесные характеристики:
   • Степень отклонения функции распределения скоростей от распределения Максвелла
   • Структура слоя Кнудсена

Методы сравнения

• Исходное уравнение Энског-Власова
• Экспериментальные данные
• Результаты моделирования молекулярной динамики

Детали реализации

• Длина расчётной области: $L$, с периодическими граничными условиями (равновесие) или поглощающими граничными условиями (неравновесие)
• Диапазон температур: от тройной точки до критической точки
• Параметры аргона: $a = 162.5 \text{ Дж·л·моль}^{-2}$, $n_c = 13.41 \text{ моль·л}^{-1}$, $T_c = 150.69 \text{ К}$

Результаты экспериментов

Основные результаты

1. Проверка равновесных свойств

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

2. Свойства ненасыщенной жидкости

При условиях вдали от кривой насыщения:

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

Анализ неравновесного испарения

1. Распределение макроскопических параметров

• Температура уменьшается от центра к границе
• Числовая плотность слегка увеличивается из-за эффекта охлаждения
• Объёмная скорость увеличивается в области раздела и газовой фазе

2. Характеристики функции распределения скоростей

Анализ в различных положениях ($x = 11, 13, 14.5, 15, 17$) показывает:

• Жидкая фаза и область раздела: распределение скоростей близко к распределению Максвелла
• Газовая фаза: явное отклонение в диапазоне отрицательных скоростей из-за недостаточного числа столкновений молекул
• Область вблизи границы: значительный эффект объёмного движения с явным смещением распределения

Экспериментальные находки

1. Ограничения соотношения Герца-Кнудсена: немаксвеллевский характер распределения скоростей в газовой фазе подтверждает неприменимость соотношения ГК при сильно неравновесных условиях
2. Два ключевых фактора:
   • Отклонение распределения скоростей в газовой фазе вблизи жидко-газового раздела от распределения Максвелла
   • Возрастающее влияние объёмного движения на газовой стороне

Связанные работы

Основные направления исследований

1. Классическая теория: соотношение Герца-Кнудсена, уравнение Шраге
2. Кинетические методы: уравнение Больцмана, метод DSMC, теория Энског-Власова
3. Молекулярное моделирование: молекулярная динамика, методы Монте-Карло

Преимущества данной работы

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

Заключение и обсуждение

Основные выводы

1. Успешно разработана улучшенная кинетическая модель для жидкостей Леннарда-Джонса
2. Модель показывает высокое соответствие с экспериментальными данными при прогнозировании равновесных свойств
3. Выявлены немаксвеллевские характеристики распределения скоростей при неравновесном испарении
4. Обоснованы теоретически ограничения соотношения Герца-Кнудсена

Ограничения

1. Область применения: модель оптимизирована в основном для насыщенных свойств, точность ограничена при ненасыщенных условиях
2. Типы жидкостей: в настоящее время применима только к инертным газам, расширение на многоатомные молекулы требует учёта внутренних энергетических мод
3. Геометрическая сложность: основная проверка проведена для плоского раздела

Направления будущих исследований

1. Расширение на другие инертные газы
2. Разработка моделей для сложных многоатомных жидкостей
3. Исследование приложений при сложных геометрических условиях

Глубокая оценка

Достоинства

1. Методологическая инновативность: искусное сочетание теории Энског и приближения Власова с эмпирической коррекцией для точного описания реальных жидкостей
2. Полнота проверки: модель проверена с нескольких углов, включая равновесные свойства и неравновесные характеристики
3. Теоретический вклад: обеспечивает важные понимания для изучения неравновесного испарения, особенно в отношении выявления ограничений соотношения ГК
4. Практическая ценность: значительно повышает точность прогнозирования при сохранении вычислительной эффективности

Недостатки

1. Зависимость от параметров: модель содержит несколько эмпирических параметров, что может влиять на её универсальность
2. Теоретическое обоснование: некоторые коррекции не имеют строгого теоретического вывода
3. Диапазон проверки: проверка проведена в основном на аргоне, применимость к другим жидкостям требует дальнейшего подтверждения

Влияние

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

Применимые сценарии

• Проектирование систем испарительного охлаждения в нанометровом масштабе
• Теплопередача при фазовых переходах в пористых средах
• Управление тепловыми потоками в микро- и наноустройствах
• Оптимизация процессов дистилляции и разделения

Библиография

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