Каковы ограничения уравнений состояния SRK и PR и как их можно устранить?

Стандартные уравнения SRK и PR систематически завышают плотность жидкости вблизи критической точки и не могут воспроизвести правильный второй вириальный коэффициент при очень высоких температурах — особенно для легких газов, таких как водород и гелий. Для термодинамических экспериментов в химической технологии это означает, что вы не можете полагаться на «сырые» данные SRK или PR при проектировании сосудов опытной установки, подборе насосов или моделировании потоков, богатых водородом при высоких температурах. Наиболее надежные исправления делятся на три категории: применение объемного сдвига, ограничение уравнения состояния известными PVT-данными или полная замена кубического уравнения специальной эмпирической корреляцией.

Хотя уравнения Соаве-Редлиха-Квонга и Пенга-Робинсона остаются превосходными универсальными инструментами для фазового равновесия, их универсальная форма содержит два жестких ограничения: фиксированный критический коэффициент сжимаемости, искажающий объемы жидкости, и член притяжения, который расходится неправильно при экстремальных температурах. Глубокая потребность заключается в закрытии этого разрыва, чтобы данные опытной установки и моделирование процессов оставались достоверными там, где это наиболее важно — при определении размеров сосудов, выборе насосов и работе с водородом при высоких температурах.

Фундаментальные недостатки стандартных кубических уравнений состояния

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

Каждое кубическое уравнение принудительно задает универсальный критический коэффициент сжимаемости Zc.
SRK фиксирует Zc на уровне 1/3, PR — на уровне 0,3074. Реальные жидкости обычно концентрируются вокруг Zc ≈ 0,27. Поскольку критический объем масштабируется прямо пропорционально Zc, стандартные кубические уравнения завышают критический объем на 10–20%, и эта ошибка распространяется на всю жидкостную ветвь.

Несоответствие наиболее сильно выражено для насыщенных жидкостей, близких к критической температуре.
В ректификационной колонне опытной установки или в эксперименте со сжиженным газом это приводит к заниженной плотности, которая может ввести в заблуждение при расчетах объема удержания, определении размеров флегмовой емкости и оценках NPSH насоса. Стандартные уравнения SRK и PR могут легко показывать ошибки, превышающие 5–10% для насыщенных жидкостей, в то время как паровая фаза остается приемлемой.

Плохой прогноз — это не неудача кубической формы, это неудача предположения о постоянном Zc.
Параметр притяжения a и совокупный объем b калибруются по давлению пара и условиям критической точки, а не по плотности жидкости. Без дополнительной степени свободы уравнение просто не может изогнуть жидкостную изохору, чтобы соответствовать реальности.

Сбой вириального коэффициента при высоких температурах

Второй вириальный коэффициент B(T) всех простых жидкостей при очень высокой температуре стремится к малому положительному предельному значению.
Этот предел происходит от вклада отталкивания (твердых сфер). Однако члены притяжения в стандартных уравнениях SRK и PR содержат температурно-зависимую альфа-функцию, которая никогда не исчезает полностью. Их подразумеваемый B(T) становится отрицательным и продолжает расти по модулю с повышением температуры — что физически невозможно.

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

Корневая причина заключается в альфа-функциях типа Соаве и Пенга–Робинсона.
Они были разработаны для улучшения давления пара в докритической области; при приведенных температурах выше 10 их алгебраическая форма экстраполируется в физически неверном направлении. Именно поэтому дополнительные справочные материалы подчеркивают, что стандартные кубические уравнения дают надежные объемные прогнозы при высоких температурах и низких давлениях — где второй вириальный член мал и доминирует закон идеального газа.

Устранение ограничений плотности жидкости в экспериментальной работе

Коррекция объемным сдвигом (Пенелу)

Простейшее и наиболее широко применяемое исправление — это постоянный объемный сдвиг.
Вычтите небольшой поправочный член c из мольного объема, предсказанного SRK или PR. Поскольку сдвиг не зависит от давления и температуры, он смещает всю кривую плотности жидкости, не изменяя расчеты парожидкостного равновесия. Многие коммерческие симуляторы реализуют это как сдвиг Пенелу или основанный на Ракетте.

В эксперименте на опытной установке это означает, что вы сохраняете те же параметры бинарного взаимодействия и качество VLE, приводя плотность жидкости в соответствие с экспериментальными данными в пределах 1–3%.
Для криогенного сепаратора или установки с СПГ эта точность часто достаточна для определения размеров сосудов и контроля запасов, при условии, что рабочая точка находится вдали от непосредственной близости к критической точке.

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

Эмпирические корреляции для проектирования процессов

Когда даже кубические уравнения с объемным сдвигом недостаточны, инженеры обращаются к специальным корреляциям плотности жидкости.
Дополнительные ссылки указывают на семейства Yu-Lu, Spencer-Danner и Chiu-Hsi-Lu, которые могут прогнозировать плотность точки кипения со средним отклонением менее 0,3% для сжиженных газов и углеводородных смесей.

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

Эти корреляции являются эмпирическими; они должны быть подобраны к химической системе и диапазону давлений.
Использование корреляции Chiu-Hsi-Lu #2 для тяжелой нефтяной фракции, например, даст плохие результаты, если параметры корреляции были регрессированы для легких потоков NGL. Всегда проверяйте корреляцию по нескольким измерениям плотности на лабораторном стенде, прежде чем доверять ей в опытном эксперименте.

Ограничение констант экспериментальными PVT-данными

Когда у вас уже есть высококачественные PVT-данные для вашей смеси, вы можете регрессировать пользовательский набор констант уравнения состояния.
Это обычная практика в исследовательских лабораториях: измерьте плотность жидкости при нескольких точках температура-давление, а затем скорректируйте параметр объемного сдвига c или даже совокупный объем b, пока уравнение состояния не воспроизведет эти точки точно.

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

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

Укрощение поведения при высоких температурах для легких газов

Модифицированные функции члена притяжения

Для исправления вириального коэффициента при высоких температурах альфа-функция должна стремиться к конечной положительной константе при T → ∞.
Уравнения, такие как Пател–Тейя, или оптимизированные формы, такие как альфа-функция Доберта-Грабовски, включают третий параметр, который позволяет вкладам отталкивания и притяжения правильно балансироваться в пределе.

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

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

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

Когда экспериментальная матрица включает доли водорода выше 20 мол.% при температурах превышающих 400 К, даже модифицированные кубические уравнения могут испытывать трудности.
Дополнительные ссылки косвенно указывают на вириальное уравнение, усеченное на втором или третьем коэффициенте, или на многопараметрические уравнения состояния, такие как BWR-11, которые routinely используются в промышленном проектировании СПГ.

Для опытной установки это не означает, что вы должны полностью отказаться от кубических уравнений.
Хорошо работает гибридная стратегия: используйте кубическое уравнение состояния для VLE и энтальпии, но переопределите расчеты фугитивности паровой фазы и объема отдельной точной корреляцией для компонента легкого газа. Многие процессные симуляторы позволяют вам определять компоненты, подчиняющиеся «закону Генри», или присоединять второе уравнение состояния к конкретному потоку.

Этот подход добавляет сложность, но гарантирует физическое поведение.
Это стандартное решение при моделировании реакторов высокотемпературного гидрирования или сосудов с гелиевой подушкой, где стандартные уравнения SRK/PR незаметно внесли бы ошибку 2–4% в коэффициент сжимаемости, которая перейдет в показания массового расхода.

Понимание компромиссов

Баланс точности и простоты

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

Не существует универсального «наилучшего» решения — есть только наилучшее соответствие вашей экспериментальной цели.
Студенческая лаборатория по VLE может позволить себе ошибку сырого SRK, так как фокус на равновесии, а не на определении размеров. Колонна промывки LN2 опытного масштаба не может, так как ошибка плотности в 5% изменяет прогнозируемый удерживаемый объем и гидравлику колонны.

Когда критическую область невозможно избежать

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

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

Правильный выбор для вашего эксперимента

Ваш термодинамический эксперимент попадет в одну из нескольких типичных категорий. Сопоставьте свой подход с доминирующей потребностью:

• Если ваш основной фокус — определение размеров оборудования для работы с жидкостью (резервуары, насосы, флегмовые емкости): Примените версию SRK или PR с объемным сдвигом с корреляцией типа Yu-Lu для окончательной проверки. Это дает быстрые, удобные для симулятора ответы с приемлемыми запасами безопасности.
• Если ваш основной фокус — работа с легкими газами при высоких температурах (водород или гелий выше 400 К): Замените альфа-функцию выражением типа Доберта-Грабовски или Патела-Тейя, или переключитесь на вириальную коррекцию второго этапа. Полагаться на стандартные SRK/PR означает ввести ошибку отрицательной сжимаемости.
• Если ваш основной фокус — академические исследования, где физическая корректность является paramount (обучение ограничениям классических уравнений состояния): Начните со стандартных SRK/PR, покажите сбой в прогнозируемой плотности жидкости и втором вириальном коэффициенте, а затем продемонстрируйте пошаговое улучшение, достигаемое с помощью объемного сдвига и модифицированных альфа-функций. Это создает глубокое понимание, необходимое вашим студентам.
• Если ваш основной фокус — опытная установка для многокомпонентного разделения, где важны и плотность жидкости, и поведение легкого газа: Ограничьте модифицированное кубическое уравнение вашими собственными экспериментальными PVT-данными с помощью регрессии по нескольким свойствам, сохраняя точность давления пара нетронутой, настраивая сдвиг объема и асимптотику высоких температур.

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

Итоговая таблица:

Оптимизируйте ваши термодинамические эксперименты с LABPARK

Вы проектируете лаборатории химической технологии или масштабируете операции опытной установки? Точное термодинамическое моделирование критично для вашего успеха.

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

Свяжитесь с LABPARK сегодня, чтобы обсудить, как наши решения для опытных установок могут повысить эффективность ваших исследований и обучения!