доктор технических наук;
doctor of technical sciences;
Рассмотрены вопросы нестационарной диффузии в слоистых структурах. При разработке конструкций аппаратов для реализации массообменных процессов необходимо учитывать, свойства вещества и характер протекаемых процессов. Сроки проектирования значительно сокращаются, а КПД аппаратов получается выше, если удается построить хорошую физическую модель и применить математический анализ с учетом кинетики процессов. Трудности теоретического анализа и расчета массопереноса определяются сложностью механизма переноса к границе раздела фаз и от нее. Поэтому применяют упрощенные модели процессов массопереноса в которых механизм массоотдачи характеризуется сочетанием молекулярного и конвективного массопереноса. Многие важные практические задачи предполагают расчет нестационарной диффузии (второго закона Фика) для определенного объема вещества (веществ). Для качественной оценки процессов, в случае симметрии, объемные задачи можно рассматривать как одномерные задачи, т.е. зависящие от одной координаты. Предложено общее решение уравнения нестационарной диффузии для слоистых сред. При этом рассматривался случай нестационарных граничных условий третьего рода на внешней поверхности и граничных условий сопряжения четвертого рода для соприкасающихся слоев. Решение получено методом разделения переменных Фурье по собственным функциям задачи с применением интеграла Дюамеля. Предложенная форма решения имеет явный вид и благодаря рекуррентной форме записи основных соотношений может быть полезной при численных расчетах.
The issues of nonstationary diffusion in layered structures are considered. When designing the devices for implementing mass transfer processes, it is necessary to take into account the properties of the substance and the nature of the processes. Design time reduces significantly and the efficiency of the devices is higher if a good physical model is built and a mathematical analysis with kinetics of the processes is applied. The difficulties of theoretical analysis and calculation of mass transfer are determined by the complexity of the transfer mechanism to and from the phase boundary. Therefore, simplified models of mass transfer processes are used in which the mass transfer mechanism is characterized by a combination of molecular and convective mass transfer. Many important practical problems involve the calculation of nonstationary diffusion (Fick's second law of diffusion) for a certain volume of substance (substances). For qualitative evaluation of processes, in the case of symmetry, volumetric issues can be considered as one-dimensional tasks, i.e. dependent on one coordinate. The general solution of the non-stationary diffusion equation for layered environments is proposed. The case of non-stationary boundary conditions of the third kind on the external surface and boundary conditions of the fourth kind conjugation for contiguous layers has been considered. The solution is obtained by separating the Fourier variables by the eigenfunctions of the problem using the Duhamel integral. The proposed solution is explicit and due to the recurrent form of the basic relations can be useful in numerical calculations
Введение. Известно, что массообменные (диффузионные) процессы характеризуются переносом одного или нескольких веществ исходной смеси из одной фазы в другую через поверхность раздела фаз. В группу процессов, которые рассматриваются, как массообменные входит молекулярная диффузия распределяемого вещества. Молекулярная диффузия определяет процессы абсорбции, перегонки (ректификации), экстракции из растворов, растворение и экстракцию из пористых тел, кристаллизацию, адсорбцию и сушку. Химические (реакционные) процессы протекают со скоростью, определяемой законами химической кинетики. В тоже время химическим реакциям сопутствует перенос массы и энергии. Поэтому скорость реакций подчиняется законам макрокинетики и определяется наиболее медленным из последовательно протекающих химического взаимодействия и диффузии [1-4].При разработке конструкций аппаратов для реализации указанных процессов необходимо учитывать, свойства вещества и характер протекаемых процессов. Сроки проектирования значительно сокращаются, а КПД аппаратов получается значительно выше, если удается построить хорошую физическую модель и применить математический анализ с учетом кинетики процессов.Трудности теоретического анализа и расчета массопереноса определяются сложностью механизма переноса к границе раздела фаз и от нее. Поэтому применяют упрощенные модели процессов массопереноса в которых механизм массоотдачи характеризуется сочетанием молекулярного и конвективного массопереноса. Среди моделей массопереноса следует выделить пленочную модель, модель диффузионного пограничного слоя, модель обновления поверхности фазового контакта и модифицированные модели обновленной поверхности. В любом случае основу всех моделей составляет основное уравнение массообмена, связывающее изменение концентрации вещества во времени с координатой точки в объеме вещества, а также условиями массобмена (диффузии) на свободных поверхностях и на границе соприкасающихся поверхностей с различной концентрацией вещества. Приходится идеализировать и упрощать рассматриваемую задачу, делая определенные допущения относительно начального распределения концентрации вещества и коэффициентов диффузии.Аппараты реализующие массообменные процессы представляют собой объемные тела различной конфигурации, для которых требуется правильно подобрать геометрические размеры и технические параметры с учетом скорости протекания процессов и производительности оборудования.Методология. Многие важные практические задачи предполагают расчет нестационарной диффузии (второго закона Фика) [1] для определенного объема вещества (веществ). Заметим, что для качественной оценки процессов, в случае симметрии, многие объемные задачи можно рассматривать как одномерные задачи, т.е. зависящие от одной координаты. В более ранних работах автора были рассмотрены аналитические решения однородной задачи нестационарной теплопроводности в слоистых структурах [5-13]. С учетом методического подхода, изложенного в этих работах, методом разделения переменных Фурье можно получить решение уравнения нестационарной диффузии для слоистых сред.Основная часть. Ниже приведено общее решение уравнения нестационарной диффузии для слоистых сред. При этом рассматривался случай нестационарных граничных условий третьего рода на внешней поверхности и граничных условий сопряжения четвертого рода для соприкасающихся слоев.Одномерная задача нестационарной диффузии для слоистых сред в математической постановке должна определяться следующей системой дифференциальных уравнений в частных производных: ,xi-1 ≤ r ≤ xi , i = 1, 2,…n, (1)где Сi(r,t), – концентрация вещества в i-м слое; Di, – коэффициенты диффузии в i-м слое; x0, – координаты нижней геометрической (свободной) поверхности объекта; xn – верхней геометрической (свободной) поверхности объекта.На внешних (свободных) поверхностях r = xo, r = xn граничные условия определим как нестационарные граничные условия третьего рода: , (2) Граничные условия сопряжения четвертого рода концентрационных полей и концентрационных потоков для соприкасающихся слоевв общем виде определяются следующим образом: , i = 12,…n-1, (3) где Ki – константа растворимости i- го слоя, Распределение концентрационных полей в начальный момент времени в каждом слое имеет вид: , i = 1, 2,…n, (4)где – функция начального распределения концентрационных полей.Искомое решение задачи (1) представим в виде суммы , (5)где – функции, которые являются решением задачи с нулевыми начальными условиями и удовлетворяют уравнениям (1)-(3).Функции можно определить интегралом Дюамеля [14, 15]: , при t>0 (6)где – решение задачи при условии, что τявляется параметром.После преобразования функции окончательно приобретают вид: , (7)где – собственные функции задачи; , i = 1,2,…n, (8) , (9) , (10)i = 2,3,…n. ,i = 2,3,…n. (11) Собственные функции задачи и определяются как корни трансцендентного уравнения: m = 0,1,2,… (12) Коэффициенты при собственных функциях задачи определяются условиями ортогональности функций, которые выполняются следующим образом: , (13) , (14) , , i = 2,3,…n (15) ,i = 1, 2,…n. (16) (17) Конкретный вид функций и , а также весовая функция , в различных системах координат определяются выражениями:а). Декартова (прямоугольная) система координат: (18)б). Сферическая система координат: (19)в). Цилиндрическая система координат: (20) Следует заметить, что при решении задач для шара или цилиндра полученное решение требует ограниченности в центре шара или на оси цилиндра. Тогда для граничных условий на свободных поверхностях необходимо вместо (2) использовать другую форму записи: , (21) Кроме того, в общем решении вместо (9) и (14) также следует полагать: , i=1,2...n (22)Остальные расчеты проводятся в соответствии с основным решением.Выводы. В заключение отметим, что представленные выше выражения определяют общее решение уравнения нестационарной диффузии для слоистых сред при нестационарных граничных условиях третьего рода на внешней поверхности и граничных условиях сопряжения четвертого рода для соприкасающихся слоев.На практике могут встречаться различные случаи физической и математической постановки задачи нестационарной диффузии для слоистых сред. Однако для одномерного случая различные частные решения могут быть сразу же определены с учетом условий (2), а также выражений (4), (5), (7) и (21) – (22).Предложенная форма решения имеет явный вид и благодаря рекуррентной форме записи основных соотношений может быть полезной при численных расчетах и анализе кинетики нестационарнойдиффузии в многослойных средах.