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

Математическое моделирование

УДК 517.958:536.2

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

А. Ю. Шупляк, С. П. Шкаруппа, А. М. Штеренберг

Самарский государственный технический университет,

Россия, 443100, Самара, ул. Молодогвардейская, 244.

E-mails: [email protected], [email protected], [email protected]

Разработана математическая модель процесса испарения многокомпонентных жидкостей для открытых систем с учётом массопереноса веществ в жидкой фазе. Особенностью данной модели является математическое описание процесса одномерного испарения многокомпонентных жидких смесей, содержащих плохо испаряющийся растворитель и n растворённых веществ. Точное аналитическое решение краевой задачи с использованием метода функций Грина получено путём приведения системы дифференциальных уравнений с движущейся границей для n растворённых веществ к n дифференциальным уравнениям с неподвижными границами.

Ключевые слова: многокомпонентная жидкость, концентрация, диффузия, испарение, математическая модель, метод функций Грина, задача Стефана.

Процесс испарения многокомпонентных жидкостей зависит от многих факторов: состава жидкой фазы, изменения ее объёма, наличия градиентов температур и концентраций, взаимного влияния процессов теплообмена и массо-обмена и т.д. Многообразие факторов, влияющих на процесс, обуславливает многообразие подходов к данной проблеме [1-6].

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

Анатолий Юрьевич Шупляк, инженер, каф. химической технологии и промышленной экологии. Светлана Петровна Шкаруппа (к.х.н., доц.), доцент, каф. химической технологии и промышленной экологии. Александр Моисеевич Штеренберг (д.ф.-м.н., проф.), заведующий кафедрой, каф. общей физики и физики нефтегазового производства.

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

В связи с этим в задачу настоящей работы входило изучение закономерностей процесса испарения многокомпонентных жидкостей в открытых системах при постоянной температуре и скорости движения газовой фазы. Особенность данного исследования заключается в разработке математической модели процесса одномерного испарения многокомпонентных жидких смесей с высокой концентрацией одного из компонентов (далее по тексту — растворитель) и с низкой концентрацией всех остальных N компонентов (далее по тексту — растворённые вещества) при постоянном движении воздушного потока над границей раздела фаз. Будем рассматривать два случая: когда растворитель испаряется и не испаряется.

При разработке математической модели будем исходить из следующих допущений:

- в начальный момент времени концентрации веществ распределены рав-

- испарение веществ с поверхности жидкости подчиняется закону Герца— Кнудсена;

- парциальное давление насыщенного пара каждого компонента над раствором подчиняется закону Рауля;

- массоперенос растворённых веществ в жидкой фазе осуществляется за счет диффузионных процессов, подчиняющихся законам Фика при постоянном коэффициенте диффузии;

- массоперенос растворителя в жидкости отсутствует;

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

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

Математическое описание испарения многокомпонентной жидкости при постоянной температуре растворённых веществ и отсутствии испарения растворителя имеет следующий вид:

номерно;

(1)

Сі(х, і) = С і = со^ при і = 0;

(2)

$ д Сі(х, і)$і^(і) + А—Сі(х,і) = .1і(і) при х = Л,(і);

(3)

Сі(х, і) = Сі = со^ при х = те;

(4)

(5)

Н(і) = 0 при і = 0;

(6)

СО = аУз (*)рУ & 1т ; (7)

/ м V1

у (£) = [Ск (ж, С\(х,г) при х = Л(^), (8)

\ и—(Л '

х = О

а; = /г(^)

Рис. 1. Схема движения границы: 1 — слой испаряющейся жидкости; 2 — слой жидкости, который испарится за время Ь

где £ — время, с; х — координата в жидкой фазе, перпендикулярная границе раздела фаз (начало координат соответствует нахождению границы раздела фаз в начальный момент времени, см. рис. 1), м; г — порядковый номер компонента (г = 0 соответствует растворителю), г = 1, 2,..., N; Сг, С — начальная концентрация и концентрация в произвольный момент времени (г/м3) г-того компонента; А, Зг, ^, Рг, рг —коэффициент диффузии (м2/с), интенсивность испарения (г/(см2)), молярная масса (г/моль), давление насыщенных паров (Па) и плотность (г/м3) г-того компонента соответственно; У — мольная доля компонента в растворе на границе раздела фаз; Л(£) — функция зависимости координаты границы раздела фаз от времени, м; N — число растворённых веществ в растворителе; а — эмпирический коэффициент, зависящий от движения воздушных масс над поверхностью жидкости.

В задаче (1)-(8) соотношение (2) представляет начальное условие, которым задаётся равномерное распределение концентрации всех компонентов в растворе в момент времени £ = 0. Соотношение (3) представляет граничное условие, согласно которому на границе раздела фаз сумма диффузионного потока, определяемого первым законом Фика [7], и потока, возникающего за счёт движения границы, г-того испаряющегося растворённого вещества равна его интенсивности испарения. Соотношение (4) характеризует равенство концентрации г-того компонента начальному значению при х ^ те. Величина йНС^/йЬ — скорость конвективного потока (стефановский поток [7]), возникающего от поступления веществ к поверхности раздела фаз за счёт изменения объёма жидкости. Выражение (6) задаёт координату границы раздела фаз «газ-жидкость» в начальный момент времени. Выражение (7) представляет собой совмещение законов Герца—Кнудсена и закона Рауля, где Ур — парциальное давление насыщенного пара испаряющегося компонента над раствором.

Краевая задача (1)-(8) представляет систему N дифференциальных уравнений с нелинейным граничным условием третьего рода (совокупность выражений (3), (7) и (8)) и соотношениями (5)-(6), описывающими движение верхней границы.

Для упрощения процесса получения аналитического решения задачи (1)-(8) перейдём к подвижной системе координат и приведём граничное условие (3), (7), (8) к линейному виду. Для перехода к подвижной системе координат с фиксированной верхней границей введём новую независимую переменную г = х — Л(£) [9]. Считая, что концентрации растворённых веществ малы по

сравнению с концентрацией растворителя, соотношение (7) в новой системе координат приведём к виду

Зг(-) = кСг(г, -) при г = 0, (9)

где к — константа скорости (м/с) испарения компонента, определяемая по формуле

к = а^ргЦС ; (10)

^о, С о — молярная масса (г/моль) и начальная концентрация (г/м3) растворителя.

Введём величину Ь как среднее значение стефановского потока за некоторый промежуток времени -тах:

1 Г ^тах У 1 ^ / 1 Г ^тах \

Ь = --- — й(-)й- = ----- — М^м). (11)

-тах ,)0 -тах ^ ^ \ рг .)0

С учётом введения подвижной системы координат и выражений (9) и (11), систему уравнений (1)—(8) можно представить в виде

д д2 д

—Ci(z, t) = Di—2Ci(z, t) + b—Ci(z, t) при 0 < z < to, (12)

dt dz2 dz

Ci(z, t) = Ci = const при t = 0, (13)

д

Di~dz,Ci(z, t) = (ki - b)Ci(z, t) при z = 0, (14)

Ci(z,t) = Ci = const при z = to, (15)

где i = 1, 2,... ,N, т.е. в виде N независимых дифференциальных уравнений

(12) с граничными условиями (13)—(15) для каждого растворённого вещества.

Решение дифференциального уравнения (12) при заданных граничных условиях в интегральной форме приведено в [9]:

р ГО

Ci(z,t) = C i Gi(z,^,t)d^,

Jo

где Gi(z,^,t) — функция Грина, определяемая выражением

G„(;.i,.t) = * exp(Ь(ф - z) - ) Lp(-(z - Ф)2) +

n ’ 2^ПЩ v 2Di 4Di) у Pv 4Dit )

( (z + ф)2) 2ki - b f™ ( (z + ф + n)2 2ki - b ) '

+ exp( - -UDDT) - — Jo exn------------------------m--------25Г nr>

После преобразований решение дифференциального уравнения (12) примет вид

Ci(z,t) = I К-Ш^ exp(-D) erfc(^) +

+ 2ki - b / (ki - b)(kit + z) \ erf/ z + (2ki - b)t \

ki

Di

2jD~t

(16)

Подстановкой выражения (16) при z = 0 в выражение (9) находим функцию интенсивности испарения от времени при условии, что стефановский поток постоянен (b = const):

Ji(t) = C i

- erfcf--AP- V L2 V 2\l DiJ + 2ki - be(ki-b)kit/Di erfc

, b ki - 2 J

(17)

Следуя [7], преобразуем второе слагаемое выражения (17) для вычисления при больших значениях времени: применим разложение функции вй^ж) в асимптотический ряд. После упрощений получим выражение

Ji(t) = C i

b erfcf - + ,/ —e

2 V 2 V Di) V nt

Di - ^ ~

i 4 Di

E

n=0

2nDn

(2ki - b)2ntn

П(! - 2j)

j=0

При больших значениях времени t имеем [7]

Е

n=0

2nDn

(2ki - b)2ntn

П(1 - 2j)

j=0

erfc(-0,5Ьд/t/Dj) & 2, exp(-0,25b2t/Di) & 1.

(18)

(19)

С учётом выражений (18), (19) получаем приближённую формулу для расчёта интенсивности испарения при больших значениях времени Ь:

Ji (t) & bCi + CiVDi/(nt).

(20)

Схематический график зависимости логарифма функции испарения (17) от времени представлен на рис. 2.

Анализ результатов, представленных на рис. 2, позволяет заключить, что в начальный момент времени (Ь = 0) интенсивность испарения растворённых веществ определяется только константой испарения (кг) (по аналогии с химическими процессами эту область можно назвать кинетической); при больших значениях времени (Ь ^ то) процесс испарения растворённых веществ можно рассматривать как стационарный:

Ji(t) = kiCi при t = 0, lim (Ji(t)) = bCi = J*.

t

(21)

(22)

где кг — константа испарения, определяемая по формуле (10); 3* — установившаяся с течением времени интенсивность испарения г-того компонента.

Равновесные концентрации растворённых веществ на границе раздела фаз (г/м3) при установившемся стационарном режиме определяются подстановкой в формулу (9) выражения (22):

C* = J*/ki = bCi/ki.

(23)

1

Выражения (17)-(22) выведены при условии постоянства стефановского потока (формула (11)), но при больших значениях времени осреднения средняя скорость стефановского потока может стремиться к нулю. Рассмотрим принципиальную возможность установления стационарного состояния, для которого характерно при больших значениях времени равенство значений функции её среднему значению за этот же промежуток времени.

При подстановке в формулу (11) зависимости (20) получаем следующее соотношение:

Ь = Ііт

гшах

1 Л

^шах V рі ^0

ЬС і + с п/Б V пі

ді

N

=

і_1

ЬС1

Рі

(24)

Из выражения (24) получаем условие перехода процесса в стационарное состояние:

N

1=1

Здесь уI = СI/р{ — объёмная доля г-того растворённого вещества на границе раздела фаз в начальный момент времени (£ = 0).

Сумма объёмных долей испаряющихся веществ при наличии неиспаряю-щегося испарителя не может быть равна единице, следовательно, предположение о возможности установления стационарного состояния без участия в испарении растворителя было ложным. В этом случае при больших значениях времени осреднения величина Ь ^ 0, следовательно, для бесконечного слоя при отсутствии испарения растворителя интенсивность испарения каждого растворённого вещества будет асимптотически стремиться к выражению, характеризующему диффузионную область протекания процесса:

Зі (і) к С і^ Бі/(пі).

(25)

При участии в испарении плохо испаряющегося растворителя его концентрация на границе раздела фаз будет увеличиваться, следовательно, задачу

Рис. 2. Схематический график зависимости логарифма интенсивности испарения от времени для слоя бесконечной глубины

г

1

шах

(1)—(8) необходимо дополнить выражением для определения объёмной концентрации растворителя через концентрации растворённых веществ:

N

yo(x,t) = Co(x,t)/po = і - ^2 yi(x,t)

(26)

i_1

где уі(х,і) —объёмные доли растворённых веществ (при і = 1, 2,...,Ы) и растворителя (при і = 0) соответственно в момент времени, отличный от начального при координате х.

С учётом испарения растворителя выражения (5) и (11) примут вид

N

d h(t) = Е[“■ Ji(t)

i_o pi

b=

dt

і

tm

N

Е

i_0

і f

_pi J0

Ji(t)dt

(27)

(28)

где Jo(t) —интенсивность пспаренпя растворителя (г/(с м2)); po — плотность растворителя (г/м3); отношение J0(t)/p0 — вклад в скорость стефановского потока процесса испарения растворителя (м/с).

Рассуждая аналогично предыдущему случаю, с учётом выражений (27) и (28) получаем модифицированное выражение (24) в виде формулы

N

b=b

+ koC*0/po,

(29)

i_1

где С0 — равновесная концентрация растворителя на границе раздела фаз (г/м3); член коС0/ро характеризует испарение растворителя в кинетической области.

Определяя отношение С0/ро из выражения (26) и подставляя его в выражение (29), находим

N г N

i_1

i_1

(ЗО)

Для нахождения объёмных долей растворённых веществ разделим концентрацию из выражения (23) для каждого компонента на его плотность:

Уі = byi/ki.

(зі)

Подставляя (3l) в (30), после преобразований получаем уравнение относительно величины стефановского потока:

N N

і - Е y і + ko Е {yi/k:

i_1 i_1

= ko.

Из выражения (32) получаем значение стефановского потока:

b=

N

Е

Ч_0

У:

max

b

При

I»£(к

г=0

скорость стефановского потока определяется константой испарения растворителя: Ь « к0.

Для слоя конечной глубины Н имеем задачу, аналогичную рассмотренной задаче (1)—(8), где вместо граничного условия (4) необходимо использовать соотношение

д

—С(х, £) = 0 при х = Н.

дх

Проведя аналогичные операции, получим задачу N независимых дифференциальных уравнений (12) с граничными условиями (13)—(15) для каждого растворённого вещества, где вместо граничного условия (15) принимается соотношение

д

—С (г, £) = 0 при г = Н — Н(Ь). (33)

Из граничного условия (33) следует, что для ограниченной области жидкой фазы при подвижной системе координат необходимо учитывать движение нижней границы. Между тем для начальной стадии процесса испарения слой можно считать бесконечным. С течением времени концентрации растворённых веществ на поверхности раздела фаз снижаются: уменьшается скорость стефановского потока. При достижении изменения концентрации г-того компонента влияние движения границы может оказаться незначительным. Введём параметр эффективной толщины слоя

N 1

г=0

7г(£)^£

рг Jо

(34)

С учётом выражения (34) граничное условие (33) приводится к соотношению

д

—С(г, £) = 0 при г = Heff.

Таким образом, с помощью упрощающего предположения (34) приводим N дифференциальных уравнений (12) с подвижной нижней границей к классическому виду.

Решение дифференциального уравнения (12) при заданных граничных условиях в интегральной форме приведено в [9]:

__ г нед

Сг(г,г) = С г Сг(г,ф,г)йф,

о

где Сг(г,ф,Ь) — функция Грина, определяемая так:

Сг(г,ф,г) =ехр(Ь(^2— г) — 4^)£ /га(г)/га(^)ехр(—АСО;

г г п=0 п

/п(г) = ео8(6>г) + (2кг — Ь) 8т(0г)/(2А0п);

Ъ

тах

в =_Ь ЩвП + (2кг — Ь)2 + 2кг — Ь + Ня Н^и — Ь)2

п 40#* 40291 + Ь2 4Д02 2 ЗО2^ ;

9п — множество положительных корней трансцендентного уравнения

tg(Heffвn) = 40г(кг — Ь)

9п = 402^п + Ь(2кг — Ь).

После преобразований решение дифференциального уравнения (12) примет вид

_ ~ en LHffj 1 '‘s™ LHeffJ ~ ■ "■

Ci(z,t) = 2Cie Heff 72 I (72 I \2 )(72 I \2) /(A , (35

7n - AbAs I (7n I As)(7n I \V(As - Ab)

7n

7n cos

zgr

Heff

I As sin

Z?r

e-Dit(^r+^2)/H2ff

_Heff

_ _ 's)(7n

где As, Ab, Ak — критериальные числа, определяемые выражениями

Л£ = Лк — Ль; Ль = Pe/2;

Afc = kiHeff/Di; (36)

Pe = bHeff/Di —диффузионный критерий Пекле; £n — множество положительных корней трансцендентного уравнения

tg(Cn) = £n(ЛЕ - ЛьV(СП + ЛЬЛТ,) (37)

на отрезках [пп; пп + п/2], п = 0,1, 2,____

Подставляя значение концентрации вещества в жидкости на границе раздела фаз, определяемой из выражения (35) при z = 0, в выражение (9), после преобразований получаем выражение для расчёта интенсивности испарения:

Т(t) = 2 DiCi V ЛкСП eXP ( - Dit[£n + Л1]/Н2еш) (38)

Heff n= - ЛЬЛЕ + (cn + Л|)(СП + Л2)/(Л£ - ЛЬ)

С течением времени ряд (38) сходится к первому члену (п = 0) и функция интенсивности испарения приводится к экспоненте (см. рис. 3):

Ji(t) = Ai exp(-A2t),

1 2С^Лк Со(ЛЕ — ЛЬ)

где

At =

Heff (As - Ab)(70 - AbAs) I (70 I As)(70 I a2)’ 222

А2 = А(Со + Аь)/НеА?.

Для того чтобы математическая модель, выраженная формулами (37) и (38), имела смысл, необходимо, чтобы у трансцендентного уравнения (37) на отрезке [0,п/2] был корень. Запишем условие нахождения корня на данном отрезке:

(Ац — Аь)/АьАт, > 1.

Рис. 3. Зависимость логарифма интенсивности испарения от времени для слоя глубины Н

Подставляем в это условие вместо А^ выражение (36) и решаем неравенство относительно Аь с учётом, что Ак и Аь больше нуля и Аь < Ак:

Аь < 1 + 0,5Ак — 0,5^1 Ак + 4. Последнее выражение можно переписать в виде

Pe < 2 I A k - A 2k I 4 .

(39)

В предельном случае Иш [\к — д Лк + 4] = 0 и Ре < 2.

Ак V

При невыполнении условия (39) процесс не выходит на экспоненту, изменение интенсивности испарения вещества во времени будет описываться выражением (17), как для слоя жидкости бесконечной глубины. Если растворитель не испаряется или значение числа Пекле настолько мало, что им можно пренебречь (Ль ^ 0), то выражение для расчёта интенсивности испарения примет вид

, 2 A 2 DiCi О exp(DiC2t/Hffr)

Ji(i)-2Akiff n=o"^IAkTa- ’

где Z определяется из решения трансцендентного уравнения

tg(C) = Ak/Z.

(40)

(41)

Из формул (17), (22), (25), (37), (38), (39), (40) и (41) для расчёта испарения следует, что если растворитель не испаряется, то конвективный (сте-фановский) поток, возникающий за счёт поступления веществ к границе раздела фаз, с течением времени прекращается, и испарение будет стремиться

к диффузионному режиму (выражения (40) и (41), при этом в начальные моменты времени, пока концевой эффект слоя не проявляется, для расчёта можно использовать выражение (25)). При испарении растворителя для веществ, диффузионный критерий которых соответствует условию (39), может установиться квазистационарное состояние, выраженное формулой (22). Длительность стационарного состояния обуславливается соблюдением условия (39). Когда условие (39) нарушится, для некоторого растворённого вещества его испарение будет описываться формулами (37) и (38), а процесс будет стремиться к стационарному испарению без участия этого компонента.

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

Исходя из модели при рассмотрении процесса испарения многокомпонентных жидкостей можно выделить три группы веществ: растворитель, растворённые вещества, для которых выполняется условие (39), и растворённые вещества, для которых условие (39) не выполняется.

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

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

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

БИБЛИОГРАФИЧЕСКИЙ СПИСОК

1. J. F. Richardson, “The evaporation of two-component liquid mixtures” // Chemical Engineering Science, 1959. Vol. 10, no. 4. Pp. 234-242.

2. А. В. Кульчицкий, Математическое моделирование процессов испарения в многокомпонентных средах: Дис. ... канд. физ.-мат. наук. М., 1999. 168 с. [A. V. Kul’chitskiy, Mathematical simulation of evaporation processes in multicomponent media: Ph.D. Thesis (Phys. & Math.). Moscow, 1999. 168 pp.]

3. R. Chebbi, M. S. Selim, “The Stefan problem of evaporation of a volatile component from a binary liquid mixture” // Heat and Mass Transfer, 2006. Vol. 42, no. 3. Pp. 238-247.

4. А. Д. Галлеев, С. И. Поникаров, “Прогнозирование зон токсической опасности и пожароопасности при авариях на объектах хранения нефтепродуктов” // Безопасность жизнедеятельности, 2009. №5. С. 29-34. [A. D. Galleyev, S. I. Ponikarov, “Prediction of toxic

and flammable hazard zones at accidents on oil product storage facilities” // Bezopasnost’ Zhiznedeyatel’nosti, 2009. no. 5. Pp. 29-34].

5. М. А. Долгова, А. Д. Галлеев, С. И. Поникаров, “Моделирование процесса испарения жидкости с поверхности аварийного пролива” // Безопасность жизнедеятельности, 2009. №6. С. 45-48. [M. A. Dolgova, A. D. Galleyev, S. I. Ponikarov, “Modeling of process the evaporation liquids from the surface accident spills” // Bezopasnost' Zhiznedeyatel'nosti, 2009. no. 6. Pp. 45-48].

6. А. Д. Галеев, С. И. Поникаров, “Оценка влияния профиля местности на формирование зон взрывоопасных концентраций при испарении нефти с поверхности аварийного пролива” // Электронный научный журнал “Нефтегазовое дело”, 2010. №2, 39. 17 с. [A. D. Galeev, S. I. Ponikarov, “Estimation of influence of terrain profile on generation of explosive cloud during evaporation from oil surface spill” // Electronic scientific journal “Oil and Gas Business”, 2010. no. 2, 39. 17 pp.]

7. D. A. Frank-Kamenetskii, Diffusion and heat transfer in chemical kinetics. New York: Plenum Press, 1969. 574 pp.; русск. пер.: Д. А. Франк-Каменский, Диффузия и теплопередача в химической кинетике. М.: Наука, 1987. 502 с.

8. Э. М. Карташов, Аналитические методы в теории теплопроводности твердых тел. М.: Высш. шк., 2001. 550 с. [E. M. Kartashov, Analytical methods in the theory of thermal conductivity of solids. Moscow: Vyssh. Shk., 2001. 550 pp.]

9. А. Д. Полянин, Справочник по линейным уравнениям математической физики. М.: Физ-матлит, 2001. 576 с.; англ. пер.: A. D. Polyanin, Handbook of Linear Partial Differential Equations for Engineers and Scientists. Boca Raton: Chapman & Hall/CRC Press, 2002. xviii+781 pp.

Поступила в редакцию 20/III/2013; в окончательном варианте — 27/VI/2013.

MSC: 80A22; 76R10, 65M80

MATHEMATICAL MODELLING OF MASS-TRANSFER UNDER THE EVAPORATION OF MULTICOMPONENT LIQUIDS

A. Yu. Shuplyak, S. P. Shkaruppa, A. M. Shterenberg

Samara State Technical University,

244, Molodogvardeyskaya st., Samara, 443100, Russia.

E-mails: [email protected], [email protected], [email protected]

The mathematical model of evaporation of multicomponent liquids is developed, for open systems with consideration of mass-transfer substances in a liquid phase. Feature of the given model is the mathematical description of one-dimensional evaporation of the multicomponent liquid mixes containing poorly volatile solvent and n of diluted substances. Exact analytical solution of the boundary value problem based on Green’s functions method is received by reduction of the differential equations system with moving boundary of the n diluted substances to the n differential equations with fixed boundary.

Key words: multicomponent liquid, concentration, diffusion, evaporation, mathematical model, Green’s functions method, Stefan problem.

Original article submitted 20/III/2013; revision submitted 27/VI/2013.

Anatoliy Yu. Shuplyak, Engineer, Dept. of Chemical Engineering & Industrial Ecology. Svetlana P. Shkaruppa (Ph.D. Chem.), Associate Professor, Dept. of Chemical Engineering & Industrial Ecology. Alexander M. Sterenberg (Dr. Phys. & Math. Sci.), Professor, Head of Dept., Dept. of General Physics and Physics of Oil and Gas Production.