Анализ линейной модели движения малоконцентрированной суспензии монодисперсных стоксовских частиц в плоском канале Текст научной статьи по специальности «Физика»

УДК 517.9

DOI: 10.14529/mmpl40405

АНАЛИЗ ЛИНЕЙНОЙ МОДЕЛИ ДВИЖЕНИЯ МАЛОКОНЦЕНТРИРОВАННОЙ СУСПЕНЗИИ МОНОДИСПЕРСНЫХ СТОКСОВСКИХ ЧАСТИЦ В ПЛОСКОМ КАНАЛЕ

В.И. Ряжских, A.A. Богер, A.B. Ряжских

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

Ключевые слова: уравнение конвективной диффузии; краевая задача; преобразование Лапласа; аналитическое решение; толщина осадки; седиментация; степень перемешивания; поглощательная способность стенки.

Введение

Исследование гидромеханики суспензий началось с работ Эйнштейна, Каннингема и Смолуховского [1-3], в которых были синтезированы феноменологические модели эффективной вязкости, единичной ячейки и отражения соответственно. Примененная ими схема реализовывала гипотезу взаимодействующих континуумов дисперсионной и дисперсной фаз при ламинарном режиме течения, что давало возможность рассматривать поведение фаз отдельно друг от друга, а их взаимное влияние учитывать через так называемые представительные объемы. Это позволяло использовать совместно уравнения Навье - Стокса и переноса частиц твердой фазы в Эйлеровом представлении модифицированными тензорными соотношениями реологических компонентов. Попытка [4] распространить такой подход на турбулентный режим течения суспензии в рамках энергетического подхода Кармана не привела к ожидаемым результатам: оставалось не ясным, каким образом получить корректную информацию о локальной концентрации твердой фазы в задачах с непроницаемыми границами и о распределении осадка на смоченных поверхностях. Параллельно с феноменологическим подходом стал развиваться и стохастический [5], однако основная проблема решений различного спектра задач динамики дисперсных систем методами статистической физики состоит в согласовании локальных характеристик фаз для однозначного поведения суспензии в целом. Все перечисленные трудности остаются и в случае малоконцентрированных суспензий, для которых влиянием твердой фазы на несущую среду можно пренебречь.

Необходимость совершенствования технологий получения композиционных материалов, используемых в энергетике, ракетно-космической технике, химической и пищевой промышленности и др., в последнее время вновь пробудило интерес к этой тематике, подробный анализ работ приведен, например, в [6-9]. Однако, судя по полученным пока результатам, проблема все еще далека от своего разрешения.

Альтернативным подходом для описания гидромеханики дисперсных систем является использование конвективно-диффузионных модельных представлений о сепарации дисперсной фазы в суспензиях, который подтвердил свою состоятельность при решении ряда важных практических задач [10, 11]. Уравнение транспорта моттодисперспых частиц имеет следующий вид

ди

— + (и ■V) и = V (БУи), (1)

где т - текущее время; и - локальная счетная концентрация; и = V + щ V - скорость движения дисперсионной среды; щ - седиментационная скорость частиц; О - коэффициент диффузии (перемешивания). В общем случае и и О являются функциями времени, координат, концентрации и размера частиц. На поверхности осадка Г, как правило, формулируется граничное условие [12]

(щи + Оgradп)|г = [к ■ и (1 — и/и3)]\г , (2)

где п3 - счетная концентрация частиц взвеси в осадке, к - кинетический коэффициент, характеризующий скорость встраивания частиц в структуру осадка. Соотношения (1) и (2), дополненные начальным условием

и\т=0 = и0, (3)

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

1. Постановка задачи и ее физическая линеаризация

Рассматривается однонаправленное движение малокопцептрироваппой суспензии моттодисперспых стоксовских частитт, в плоском слое, тта входе которого осуществляется подача суспензии с постоянной концентрацией (рис. 1).

//////////////////

о о о —- V > о а о ° о 0 ,

Рис. 1. Расчетная схема

В силу малокопцептрироваппости суспензии влиянием частитт, тта гидродинамическую структуру течения можно пренебречь, так что определение и возможно из несопряжен-ттых с полем концентрации частитт, уравнений Навье - Стокса или уравнений Рейттольдса. Известно [13], что довольно сложттое решение гидродинамической задачи может быть заменено гидродинамическим режимом идеального вытеснения, который в первом приближении

вполне адекватно описывает, например, турбулентное течение маловязких сред, откуда следует, что V ~сопз^ Далее, так как частицы имеют стоксовскую природу, то и ад ~соп8^ Кроме того, будем считать, что О ~соп8^ а также не учитывать эффект скольжения фаз. Таким образом, в выбранной системе координат (рис. 1) уравнение (1) запишется

дп (х,х,т) дп (х,х,т) дп (х,х,т) ^

-V = -V- ' + ш- ' + О

дт дх дх

д2п (х,х,т) д2п (х,х,т) дх2 дх2

(4)

В силу квазистационарности рассматриваемой задачи д/дт = 0, а диффузионный перенос частиц в направлении ог существенно меньше конвективного, получим из (4)

дп (х,х) дп (х,х) о д2п (х,х) дх дх дх2 '

Граничные условия (2) и (3) с учетом того, что п ^ п3, модифицируются в следующий вид:

п (х, 0) = по; (6)

шп (0, х) + Вдп(°х'Х) = кьп (0, х); (7)

шп (Н, х) + = -кип (Н,х), (8)

дх

кь ки

в осадок на нижней и верхней стенках канала соответственно.

Для сокращения числа независимых переменных в системе (6) - (8) приведем ее к безразмерной форме записи:

дм (X, z) = дм (Х,г) -1 д2м (X, z) д2 дх + дх2 '

м (X, 0) = 1; (10)

М (0, Z) + В о-1 дМ^) = КьМ (0, Z); (11)

М (1^)+В о-1 дМ ^) = -Ки М (1^) , (12)

где Z = шх/^Н)] X = х/Н; Во = шН/О - число Боденштейна; Кь = кь/щ Ки = ки/ш. Отметим, что краевая задача (9) - (12) является линейной, и, следовательно, ее решение может быть получено классическими методами, например, с помощью одностороннего преобразования Лапласа [14].

2. Аналитическое решение

Изображение системы (9) - (12) по Лапласу относительно переменной Z таково: й2Мь (X, в) „ <1МЬ (X, в) „ лт _ .

-X2 + Во--Во^ 8 • Мь (X, в) = -Во; (13)

МЬ (0, в) + В о-1 ^^^ = КьМь (0, в) ; (14)

МЬ (1,в)+Во-1 ^^^ = КиМЬ (1,в) , (15)

ГД6 в ; N. - изображения Z и N соответственно. Общее решение (13) - (15) есть

N. (Х, в) = 1+ [Ох эИ (вХ) + С2 еИ (вХ)] ехр X

2

где в = ^/в°2/4 + в - Во; Сх, С2 - константы интегрирования, которые определены из условий (14) и (15). Поэтому

1 Ф (Х в)

(16)

N. (X, в) = - + ^^, в ф (в)

Ф (Х, в) = < <| (1 + Ки )(* - К,) ехр - (1 - Кь)

(2 + К^ - еИв+

+Во"хвэИ в]} + (1 - Кь)

эИ (вХ) _в

эИ в

-В о"1 (1 + Ки )ехр — +

(?)

-+ К^ ^ + В<Г1еЪ в

в

еЪ(вХ)^ехр х) ;

ф (в) = в^В о"1 (Ки + Кь) еИ в +

Во"2в2 -( - - К,) (2 + Ки

эИ в\

ТУ

Нетрудно показать, что ф (Х, в) и ф (в) являются обобщенными полиномами, причем степень ф (в) больше степени ф (Х, в), поэтому воспользуемся второй теоремой обращения Лапласа [15]. Для этого найдем корни уравнения ф (в) = 0, которыми будут: 1) в = 0; 2) действительный корень Л уравнения

Ш Л = рЛ/ (д - Л2) ; (17)

и 3) бесконечное счетное множество корней ¡п уравнения

tg ¡п = ри-п! (д - ¡П), п = 1, ^

(18)

р = Во (Ки + Кь); д = В о2 - - Кь

(- - Кь) (2+ Ки).

При р > 0 д > 0и р < д уравнение (17) имеет один корень в интервале (0, ^/д), а при р > 0, д > 0 и р > д или когда р > 0 д < 0 у уравнения (17) действительных корней нет. При р > 0, д > 0 и р < д уравнение (18) имеет по одному корню в интервалах периодичности функции тангенс, причем в интервале (0, п/2) действительных корней нет. При р > 0 д > 0и р > д уравнение (18) имеет по одному корню в интервалах периодичности функции тангенс, а также корень в интервале (0, п/2). И, наконец, при р > 0 д < 0 уравнение (18) на одном из интервалов периодичности функции тангенс имеет два корня. Номер этого интервала определяется величиной ^[д.

Из вышесказанного следует, что оригинал (16) есть

N (Х, Z) =

М (Х, Z) =

Мл (Х, Z) + N (Х, Z), при р > 0, д> 0 и р< д; (Х, Z), при р > 0, д > 0, р > д или при р > 0, д < 0,

(19)

Ф эИ (ЛХ) + Б л еИ (ЛХ) Л

ехр

Л2

-т +и - т

)z /ал,

N. (Х,2) = £

п=1

А

й1п (рпХ) + В. ссй (рпХ)

п

ехр

^X-( Ф- + 2 \Во 4 ) .

Во

т

с

.1

Ал = (1 + Ки) - К^ ехр ^ - (1 - КЬ)

Вх = -Во-1 Л (1 + Ки) ехр ( ^ ) + (1 - Кь)

2 + Ки) сИ Л + Во-1Л йИ Л 2+ Ки) ^ +ВО-1 сИ Л

2

Л2 Во\ ( 1

Сл = ( ^ - |2л (Ки + Кь)8ЬЛ + Во-1Л-1 8ИЛ+

Во

2л2

Во2

-(2 - *)(

2 + Ки

(сИ Л - Л-1 йИ Л)} ;

А. = (1 + Ки) ^ - К^ ехр - (1 - Кь) ^ + Ки) ссй¡п - Во-1^ й1п¡п

В. = -Во-1 (1 + Ки) ехр ) + (1 - Кь)

1 \ й1п ¡п -1

2 + Ки )--+ В о ссй ¡¡п

2 / ¡п

(ё+?){

С. = - ^ ^ (Ки + Кь)

й1п ¡п Во

1

+

¡п

й1п ¡п+

+

Во

2Лп

¡к + (!_ Во2 + \ 2

(1 - *) (2+*,)](■

сСЙ ¡п

Й1П ¡¡п ¡п

Осредненное по поперечному сечению значение локальной концентрации (19) таково:

N (^) = / Nл (2) + N. (2), при р> 0, д> 0 и р < д; ^

\ N. (2), при р > 0, д > 0 и р > д или при р > 0, д < 0,

^(2 )= С-Ч{ АЛ + Вл

+ (Л - ВО [йЫ"4- (л + Т)

Л

ехр

1

2Л

ехР ( Л - ^ ) - 1

+

1

ехр

( ^ - —) 2

/Л2_Во VBo 4

N. (2) = Е

п=1

А.

¡п

/Во

¡п V йт ¡п + ¡п сей ¡п

_Во\

+В.

хехр

)6ХР( 2 ) ) ^ (-Т) (Ё + 2\1\.( *)])■

+

Во /Во

— - I — сой ¡п - ¡п Й1П ¡п

Наиболее важной характеристикой движения суспензии в канале является распределение толщины осадка частиц на смоченных поверхностях. Расчетные соотношения для толщины осадка на нижней стенке канала получены из следующих соображений. Из (7) следует, что локальный штучный поток частиц на поверхность осаждения есть

3 (0, г) = кьи (0, г), откуда их число на элементарной площади dzdy будет

¿ш (0, г) = кьи (0, г) dтdzdy,

где - элемент длины в направлении плоскости оху. С другой стороны, объем одной ча-

стицы равен Ур = ку I3, где ку - коэффициент формы частицы с характерным размером I, тогда общий объем частиц на элементарной площади составит

(У (0, х) = Ур(ш (0, х) = ку¡'3кьи (0, х) (твхву.

Учитывая порозность осадка через е, получим величину толщины осадка на элементарной

(б (0,х) = кьп (0, х) (21)

Соотношение (21) в безразмерном дифференциальном виде таково

КЬМ {0,2) , (22)

где А(0,х) = б (0,х)/Ыс); к которому добавляется очевидное условие

А (0, го) = 0. (23)

Из (22) и (23) окончательно получим

А (0, 2) = КьIN (0, 2)(2. (24)

те

Аналогичное соотношение может быть записано для безразмерной толщины осадка на верхней стенке канала

А (1, 2) = Ки IN (1, 2)(2. (25)

те

3. Результаты расчетов

Целью расчетов являлось установление закономерностей кинетики осаждения твердой фазы в суспензии в зависимости от интенсивности перемешивания дисперсионной среды и влияние поглощательной способности стенок канала при образовании

осадка.

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

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

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

локальная концентрация частиц болытте у верхней стенки. Этим объясняется и то, что толщина осадка на верхней смоченной поверхности болытте, чем тта нижней (рис. 4 б).

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

Выявленные закономерности позволяют управлять процессом седиментации твердой фазы суспензии.

а б

Рис. 2. Локальная концентрация (а) при различных Z: 1-0; 2-1,0; 3-3,0; 4-4,0; 5-10,0; средняя концентрация по длине канала и толщины осадков на верхней (Ли) и нижней (Дь) стенках (б) для Во=0,1; Ки = 0,1; Кь = 0,1

Рис. 3. Локальная концентрация (а) при различных Z: 1-0; 2-0,1; 3-1,0; 4-2,0; средняя концентрация по длине канала и толщины осадков на верхней (Ли) и нижней (Ль) стенках (б) для Во=10,0; Ки = 0,1; Кь = 0,1

Рис. 4. Локальная концентрация (а) при различных Z: 1-0; 2-0,1; 3-3,0; 4-4,0; 5-10,0; средняя концентрация по длине канала и толщины осадков на верхней (Ди) и нижней (Дь) стенках (б) для Во=0,1; Ки = 10 Кь = 10

Рис. 5. Локальная концентрация (а) при различных Z: 1-0; 2-0,1; 3-0,3; 4-0,6; 5-0,8; средняя концентрация по длине канала и толщины осадков на верхней ( Ди) и нижней (Дь) стенках (б) для Во=10; Ки = 10 Кь = 10

Работа проводилась при финансовой поддержке РФФИ по гранту Ж'1 18-08-00261.

Литература

1. Einstein, А. Eine neue Bestimmung der Molekuldimensionen ./' A. Einstein /7 Annalen der Physik. - 1906. - V. 19. - P. 289-306.

2. Cunningham, E. On the Velocity of Steady Fall of Spherical Particles Through Fluid Medium / E. Cunningham // Proc. R. Soc. Loud. A. - 1910. - V. 83. - P. 357-365.

3. Smoluchowski, M. Uber die Wechselwirkung von Kugeln, die sich in einer zähen Flüssigkeit bewegen / M. Smoluchowski /'/ Bull. Int.. Acad. Sei. Cracovie. Ser. 1A. - 1911. - P. 28-39.

4. Einstein, A. Effects of Heavy Concentration Near the Bed on the Velocity and Sediment Distribution / A. Einstein, N. Chien. - MRD Sediment Series 8, University of California, Berkeley, 1955.

5. Покровский, В.П. Статистическая механика разбавленных суспензий / В.П. Покровский. - М.: Наука, 1978. - 135 с.

6. Junke, G. Turbulent Velocity Profiles in Sediment-Laden Flows / G. Junke, Y.J. Pierre // Journal of Hydraulic Research. - 2001. - V. 39, № 1. - P. 11-23.

7. Ventura, A.F. Modeling the Turbulent Flow of Pulp Suspensions / A.F. Ventura, A.P. Garcia, P.J. Ferreira, M.G. Rasteiro // Ind. Eng. Chem. Res. - 2011. - V. 50, issue 16. - P. 9735-9742.

8. Shanliang, Z. Numerical Research on the Fiber Suspensions in a Turbulent T-Shaped Branching Channel Flow / Z. Shanliang, L. Jianzhong, Z. Weifeng // Chinese Journal of Chemical Engineering. - 2007. - V. 15, № 1. - P. 30-38.

9. Гуськов, С.Б. Метод самосогласованного поля применительно к динамике вязких суспензий / С.Б. Гуськов // Прикладная математика и механика - 2013. - Т. 77, вып. 4. -С. 557-572.

10. Нигматулин, Р.Н. Основы механики гетерогенных сред / Р.Н. Нигматулин. - М.: Наука, 1978. - 336 с.

11. Процессы и аппараты химической технологии. Явления переноса, макрокинетика, подобие, моделирование, проектирование: В 5 т. Т. 1. Основы теории процессов химической технологии / Под ред. A.M. Кутепова. - М.: Логос, 2000. - 480 с.

12. Харин, В.М. К теории осаждения / В.М. Харин, В.И. Ряжских // Теоретические основы химической технологии. - 1989. - Т. 23, № 5. - С. 651-658.

13. Цветков, ф.ф. Тепломассообмен / ф.ф. Цветков, Б.А. Григорьев. - М.: Изд-во МЭИ, 2005. - 550 с.

14. Дёч, Г. Руководство к практическому применению преобразования Лапласа и z-преобразования / Г. Дёч. - М.: Физматгиз, 1971. - 288 с.

15. Беляев, Н.М. Методы теории теплопроводности. В 2-х ч. Ч. 1 / Н.М. Беляев, А.А. Рядно. - М.: Высш. шк., 1982. - 327 с.

Виктор Иванович Ряжских, доктор технических наук, профессор, кафедра «Прикладная математика:», Воронежский государственный технический университет (г. Воронеж, Российская Федерация), ryazhskih@scientist.com.

Андрей Александрович Богер, кандидат технических наук, доцент, кафедра «Математика», Военный учебно-научный центр Военно-воздушных сил «Военно-воздушная академия имени профессора Н.Е. Жуковского и Ю.А. Гагарина» (г. Воронеж, Российская Федерация), a-boger@yandex.ru.

Александр Викторович Ряжских, кандидат физико-математических наук, доцент,

«»

ский государственный технический университет (г. Воронеж, Российская Федерация), ryazhskihav@bk.ru.

Поступила в редакцию 18 сентября 2014 г.

Bulletin of the South Ural State University. Series "Mathematical Modelling, Programming & Computer Software",

2014, vol. 7, no. 4, pp. 65-75.

MSC 35K60 DOI: 10.14529/mmpl40405

A Linear Model of the Motion of a Low-Concentration Suspension of Monodisperse Stokes Particles in a Flat Channel

V.l. Ryazhskih, Voronezh State Technical University, Voronezh, Russian Federation, ryazhskih@scientist.com,

A.A. Boger, Military training and scientific center of the air force "Air force Academy named after Professor N.E. Zhukovsky and Y.A. Gagarin", Voronezh, Russian Federation, a-boger@yandex.ru,

A.V. Ryazhskih Voronezh State Technical University, Voronezh, Russian Federation, ryazhskihav@bk.ru

In the framework of the convection-diffusion approach to the monodisperse low-concentration solid phase sedimentation in a suspension moving in a flat horizontal channel, we obtain a linear boundary value problem for a parabolic equation on the local concentration of particles. We impose boundary conditions of the third kind using the condition that the flux of particles on the wetted surface is proportional to their concentration at the wall. Integral transformations yield an analytical solution of the stated boundary value problem, which we use to find the relations to determine the thickness of sediment on the bottom and top walls of the channel. Simulations show that the kinetics of the solid phase settling from a flowing suspension, as well as the sediment formation rate and its distribution on the botton and top walls of the flat channel, substantially depend on the degree of mixing of the dispersion medium and the absorption capacity of wet surfaces. We establish that for walls with low absorption capacity a decrease in the mixing intensity reduces the rate of particle sedimentation on the walls, but increases it in the case of high absorption capacity.

Keywords: convection-diffusion equation; boundary value problem; Laplace transform; analytical solution; sediment thickness; sedimentation; degree of mixing; wall absorption capacity.

References

1. Einstein A. Eine neue Bestimmung der Molekuldimensionen. Annalen der Physik, 1906, vol. 19, pp. 289-306. DOI: 10.1002/andp.l9063240204

2. Cunningham E. On the Velocity of Steady Fall of Spherical Particles Through Fluid Medium. Proc. R. Soc. Lond. A, 1910, vol. 83, pp. 357-365. DOI: 10.1098/rspa.l910.0024

3. Smoluchowski M. Uber die Wechselwirkung von Kugeln, die sich in einer zähen Flüssigkeit bewegen. Bull. Int. Acad. Sei. Cracovie. Ser. 1A, 1911, pp. 28-39.

4. Einstein H.A., Chien N. Effects of Heavy Concentration Near the Bed on the Velocity and Sediment Distribution. MRD Sediment Series 8, University of California, Berkeley, 1955.

5. Pokrovskiy V.P. Statisticheskaya mehanika razbavlennyh suspenziy [Statistical Mechanics of Dilute Suspensions]. Moscow, Nauka, 1978. 135 p.

6. Junke G., Pierre Y.J. Turbulent Velocity Profiles in Sediment-Laden Flows. Journal of Hydraulic Research, 2001, vol. 39, no. 1, pp. 11-23. DOI: 10.1080/00221680109499798

7. Ventura A.F., Garcia A.P., Ferreira P.J., Rasteiro M.G. Modelling the Turbulent Flow of Pulp Suspensions. Ind. Eng. Chem. Res., 2011, vol. 50, issue 16, pp. 9735-9742. DOI: 10.1021/iel01500h

8. Shanliang Z., Jianzhong L., Weifeng Z. Numerical Research on the Fiber Suspensions in a Turbulent T-Shaped Branching Channel Flow. Chinese Journal of Chemical Engineering, 2007, vol. 15, no. 1, pp. 30-38. DOI: 10.1016/S1004-9541(07)60030-5

9. Gus'kov S.B. [Self-consistent Field Method is Applied to the Dynamics of Viscous Suspensions]. Prikladnaja Matematika i Mehanika [Applied Mathematics and Mechanics], 2013, vol. 77, no. 4, pp. 557-572. (in Russian).

10. Nigmatulin R.N. Osnovy mehaniki geterogennyh sred [Fundamentals of Mechanics of Heterogeneous Media]. Moscow, Nauka, 1978. 336 p.

11. D.A. Baranov, A.V. Vjaz'min, A.A. Tuhman Processy i apparaty himicheskoy tehnologii. Javleniya perenosa, makrokinetika, podobie, modelirovanie, proektirovanie. T. 1. Osnovy teorii processov himicheskoy tehnologii [Processes and Devices of Chemical Technology. Transport Phenomena, Macrokinetics, Likeness, Modelling, Design. Fundamentals of the Theory of Chemical Technology Processes]. Moscow, Logos, 2000. 480 p.

12. Harin V.M., Ryazhskih V.l. [On the Theory of Deposition]. Teoreticheskie osnovy himicheskoy tehnologii [Theoretical Foundations of Chemical Engineering], 1989, vol. 23, no. 5, pp. 651-658. (in Russian)

13. Cvetkov F.F., Grigor'ev B.A. Teplomassoobmen [Heat and Mass Transfer]. Moscow, Izd-vo MJel, 2005. 550 p.

14. Dötsch G. Anleitung zum praktischen gebrauch der Laplace-transformation und der z-transformation, dritte Auflage. Wien, 1967.

15. Belyaev N.M., Ryadno A.A. Metody teorii teploprovodnosti [Methods of the Teory of Heat Conduction. Ch. 1]. Moscow, Vysshaya shkola, 1982. 327 p.

Received September 18, 2014