Belgorod, Belgorod, Russian Federation
GRNTI 61.01 Общие вопросы химической технологии и химической промышленности
BBK 35 Химическая технология. Химические производства
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
nonstationary diffusion, Fick's law, layered structures, nonstationary boundary conditions of the third kind, boundary conditions of conjugation of the fourth kind.
Введение. Известно, что массообменные (диффузионные) процессы характеризуются переносом одного или нескольких веществ исходной смеси из одной фазы в другую через поверхность раздела фаз. В группу процессов, которые рассматриваются, как массообменные входит молекулярная диффузия распределяемого вещества. Молекулярная диффузия определяет процессы абсорбции, перегонки (ректификации), экстракции из растворов, растворение и экстракцию из пористых тел, кристаллизацию, адсорбцию и сушку. Химические (реакционные) процессы протекают со скоростью, определяемой законами химической кинетики. В тоже время химическим реакциям сопутствует перенос массы и энергии. Поэтому скорость реакций подчиняется законам макрокинетики и определяется наиболее медленным из последовательно протекающих химического взаимодействия и диффузии [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).
Предложенная форма решения имеет явный вид и благодаря рекуррентной форме записи основных соотношений может быть полезной при численных расчетах и анализе кинетики нестационарнойдиффузии в многослойных средах.
1. Rudobashta S.P., Kartashov E.M. Diffusion in chemical-technological processes. M.: KolossS, 2010, 478 p.
2. Kasatkin A.G. The main processes and apparatuses of chemical technology. M.: OOO TID "Alliance", 2004, 753 p.
3. Romankov P.G., Frolov V.F., Flisyuk OM Methods for calculating the processes of chemical technology apparatuses. SPb.: Himizdat, 2009, 544 p.
4. Romankov P.G., Frolov V.F. Mass transfer processes of chemical technology. Systems with a dispersed solid phase. L.: Chemistry, 1990. 388 p.
5. Vendin S.V. To the calculation of non-stationary heat conduction in multilayer objects with boundary conditions of the third kind. IFZH, 1993, vol. 65, no. 2, pp. 249-251.
6. Vendin S.V. There are also a number of problems in the field of thermal engineering and thermophysics, 1993, vol. 65, no. 2, pp. 823.
7. Vendin S.V., Scherbinin I.A. To solving problems of non-stationary heat conduction in layered media. Bulletin of Belgorod State Technological University V.G. Shukhov, 2016, № 3, pp. 96-99.
8. Vendin S. Conduction in Layered Media. International Journal of Environmental and Science Education, 2016, vol. 11, no. 18, pp. 12253-12258.
9. Vendin S.V., MamontovA.Yu. Calculation of the power of additional heat sources for heating biomass in a biogas reactor. Bulletin of BSTU named after V.G. Shukhov, 2017, no. 7, pp. 97-99.
10. Vendin S.V. Towards solving some boundary problems of unsteady heat conduction in layered media by separating variables. In collection: Actual problems of drying and thermo-moisture treatment of materials in various industries and the agro-industrial complex; Lykov. Moscow, 2015, pp. 78-80.
11. Vendin S.V. Solution of the problem of non-stationary heating of layered materials. In the collection: Problems and prospects of innovative development of agricultural technologies. Proceedings of the XX International Scientific and Production Conference. Maysky: Belgorod State Agrarian University, 2016, pp. 15-16.
12. Vendin, S.V. Theory and mathematical methods for the analysis of thermal processes in the microwave treatment of seeds. M..: OJSC Central Library Collector BIBKOM, TRANSLOG LLC, 2016, 143 p.
13. Vendin S.V. International Journal of Applied Engineering Research On Solving the International. Journal of Applied Engineering Research, 2017, vol. 12, no. 22, pp. 12272-12274
14. Kartashov E.M. Analytical methods in the theory of thermal conductivity of solids: Proc. Benefit. Ed. 3rd, Pererab. and add. M.: Higher School, 2001, 550 p.
15. Korn G., Korn T. Mathematics Handbook for Scientists and Engineers. M.: Nauka, 1984, 835 p.
