Введение
Монография написана по новому подходу к расчету методом конечных элементов (МКЭ) массообменных процессов таких как ректификация, абсорбция, теплообменных процессов, процессов в химических реакторах с учетом протекания химических процессов.
1. Касаткин известной работе [1] отметил, что критериальные зависимости вводятся из-за невозможности решения дифференциальных уравнений гидродинамики. Проблема решения систем уравнений заключалась в отсутствии мощных вычислительных средств, имеющихся в настоящее время.
2. Нобелевский лауреат академик Ландау в работе [2,с.12] указывал о недостаточном физическом обосновании критериальных методик по сравнению с решением дифференциальных уравнений методами вычислительной гидродинамики.
Существующие работы по технологическому расчету аппаратов и процессам переработки нефти и газа [3], [4], [5], [6], [7], [8], [9] химической технологии содержат критериальные методики расчета таких процессов как теплообмен, ректификация, абсорбция. Отдельно рассматриваются реакторы с химическими процессами [10]. На определенном историческом этапе эти работы были востребованы в связи с отсутствием вычислительной техники. Этим и обосновывается широкое применение критериальных методик. Но в настоящее время методики, приведенные в этих работах являются устаревшими. Так как с применением современной вычислительной техники и программных пакетов получают решение дифференциальных уравнений численным методом без необходимости применения теории подобия и критериальных методик. Решение численным методом дифференциальных уравнений является более обоснованным физически по сравнению с решением критериальных уравнения, содержащих эмпирические данные. Расчет процесса в общем виде заключается в решении системы уравнений гидродинамики, химической кинетики, теплопередачи, массообмена и др. в узлах сетки с получением точного решения. . Кроме того, результат численного расчета представляется в виде цветовой диаграммы со шкалой, то есть предоставляет более наглядные данные, которые в том числе можно сравнить с фотографиями и данными результатов эксперимента. Расчет численными методами позволяет разработать цифровой двойник изделия.
На современном уровне балансы объектов, соответствующих аппаратам, на технологической схеме рассчитываются в специализированных программах, таких как пакет Aspen. В этих программах существуют модули технологического расчета аппаратов по критериальным методикам с получением технологических характеристик оборудования и геометрических размеров. Или используются отдельные программы, в результате которых получают баланс аппарата и его эскиз с размерами. Эскиз передается в чертежную программу или пакет 3D-моделирования для проектирования аппарата. Так обеспечивается сквозное проектирование. Отметим, что в процесс разработки чертежей аппарата, как и других машиностроительных изделий может быть включена передача геометрии для расчета методом конечных элементов и затем возврата полученных результатом в разработку проекта конструкции аппарата. Расчет численным методом можно выполнить для результатов расчета по критериальным методикам и затем передать данные на проектирование конструкции или без выполнения расчета по критераильным методикам с использованием данных балансов технологической схемы.
В технологии нефти и газа, в химической технологии, процессы в аппаратах являются физико-химическими, которые можно разделить на составляющие гидродинамические, тепловые, массообменные и химические. Задача расчета физико-химического процесса является междисциплинарной так как требует совместного решения всех процессов, из которых состоит физико-химический процесс, протекающий в аппарате. Итак, в общем случае составляется система из дифференциальных уравнений, описывающих каждый процесс, и находится численное решение в узлах сетки.
В расчетах таких процессов как сгорание топлива в камере анализируется по термодинамическим параметрам кинетика процесса (механизмы реакций, то есть последовательность промежуточных частиц от сырья до продуктов). Отметим, что для установления механизма реакции по наименьшей энергии промежуточных частиц можно использовать отдельную программу для квантово-химических расчетов. Однако, для процесса необходимо описывать кинетику на макроскопическом уровне с учетом гидродинамики, то есть использовать химическую гидродинамику. Подход в современных пакетах с анализом пути реакции и моделями эквивалентных реакторов по-видимому уступает расчету по уравнениям химической гидродинамики.
Система уравнений из приведенных выше уравнений гидродинамики расширяется дифференциальными уравнениями химической кинетики или процессов тепломассообмена. Затем эта система решается методом конечных объемов (численным методом) в узлах расчетной сетки аналогично тому, как в этих узлах может решаться исходная система уравнений гидродинамики.
1 Теория гидродинамики
Из теоретической гидродинамики используются уравнения, из которых составляется система. Затем эта система уравнений может быть дополнена дифференциальными уравнениями теплообмена, массообмена, химической гидродинамики. Система решается численным методом в узлах сетки.
В гидродинамике описание движения потока жидкости производится описанием вектора скорости и двух термодинамических величин (давление и плотность) в зависимости от координат и времени. Система уравнений должна содержать 5 уравнений, в том числе уравнение неразрывности, уравнение Навье-Стокса (уравнение Эйлера для идеальной жидкости), уравнение переноса тепла (уравнение сохранения энтропии для идеальной жидкости).
– уравнение неразрывности потока:
– уравнение Навье-Стокса
для сжимаемой жидкости (уравнение движения вязкой среды):
для несжимаемой жидкости при :
– уравнение переноса тепла
Закон сохранения для идеальной жидкости (при отсутствии вязкости и теплопроводности правая часть уравнения становится равной нулю и получится уравнение сохранения энтропии):
Уравнения Навье-Стокса в настоящее время решены для нескольких простых случаев, например, течения Пуазейля, течения Куэта.
Анри Навье вывел уравнения используя представления, то есть расчетную модель, о молекулярных взаимодействиях. Эта модель не отражает картины турбулентного потока, предложенной Колмогоровым и указанной справедливой Ландау.
Академик А.Н. Колмогоров записал уравнения для турбулентного движения [11]:
По модели Колмогорова А.Н. Колмогоров в работе [2,с.294] происходит наложение на осредненный поток различных по масштабу турбулентных пульсаций. Наибольшим масштабом является масштаб L пути перемешивания, наименьшим масштабом является масштаб λ, на котором вязкость уже оказывает влияние. Пульсации передаются от крупного масштаба L к меньшим и на самом низком масштабе λ происходит рассеяние энергии за счет вязкости.
2 О решении уравнений Навье-Стокса для пространства R3 в постановке математического института Клея
В работе Ефанова [12] показана попытка доказательства невозможности решения уравнений Навье-Стокса.
Для этого Ефановым рассмотрены расчетные модели и применена теорема Курта Геделя о неполноте непротиворечивых систем.
Применяя эту теорему, модель Колмогорова включает в себя представления Анри Навье и является расширенной системой по отношению к расчетной модели Анри Навье.
И следовательно, средствами базовой системы нельзя получить решение для расширенной системы.
А. Эйнштейн, А.Н. Колмогоров (в кн. «О профессии математика») писали, что теория строится на гипотезах, наглядных представления. Конкретно сопоставление основания уравнений Навье-Стокса с возможностью их решения выполнено К.В. Ефановым и указано о большем обосновании модели Колмогорова и отсутствие решения уравнений Навье-Стокса, применительно к пространству R3.
Точное решение численным методом DNS уравнений Навье-Стокса приводит в интегральному масштабу L, который похож на масштаб L в модели Колмогорова. И именно по этой причине возможно решение методом DNS уравнений Навье-Стокса. Однако, это решение теоретически строго не корректно, так как требуется решать уравнения, составленные на модели Колмогорова, но не на представлениях Анри Навье о молекулярных взаимодейсвтиях.
Таким образом, по-видимому проблема уравнений Навье-Стокса решена Ефановым К.В. путем обоснования некорректности теоретического основания уравнений для пространства R3.
3 Химическая гидродинамика
Химическая гидродинамика в дифференциальных уравнениях для макроскопических масштабов позволяет уточнить влияние на гидродинамические процессы протекания физико-химических процессов таких как химические реакции, диффузия. Химическая гидродинамика может использоваться также для анализа химической кинетики в условиях процесса.
Для химической гидродинамики существует в части переноса вещества в потоке аналогия с процессами переноса тепла конвективным теплообменом (для химических превращений присутствует большее количество параметров и граничных условий).
Диффузионная кинетика используется для анализа гетерогенных химических реакций в направлении изучения истинной кинетики на границе раздела и молекулярной и конвективной (турбулентной) диффузии вещества к границе раздела.
Используются 3 подхода:
– решение дифференциальных уравнений диффузии в граничных условиях с заданием кинетики реакции на границе раздела, граничные условия [13]:
Строгое решение получают для ламинарного потока, для случая турбулентного потока применяются приближенные методы турбулентности. Аналитические задачи диффузионной кинетики решают применением преобразования Лапласа.
– применение теории подобия к процессам переноса и к химическому процессу. Возникает проблема совместимости подобия процессов переноса и химических процессов.