Совместное движение бинарной смеси и вязкой теплопроводной жидкости под действием термоконцентрационных сил

Вычислительные технологии

Том 18, № 3, 2013

Совместное движение бинарной смеси и вязкой

теплопроводной жидкости под действием

>к

термоконцентрационных сил*

М. В. ЕФИМОВА

Институт вычислительного моделирования СО РАН, Красноярск, Россия

e-mail: efmavi@icm.krasn.ru

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

Ключевые слова: начально-краевая задача, преобразование Лапласа, бинарная смесь, термоконцентрационный эффект.

Введение

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

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

* Работа выполнена при финансовой поддержке РФФИ (грант № 11-01-00283) и проектов № 38, № 116 фундаментальных исследований СО РАН.

1. Постановка задачи

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

а(9, с) = ао - ^(9 - 9о) - «2(с - со),

где а0 = а(90,с0), > 0 — температурный коэффициент, ж2 — концентрационный коэффициент (обычно ж2 < 0, поскольку поверхностное натяжение увеличивается с ростом концентрации).

Будем считать, что в начальный момент времени смесь заполняет слой -11 < у < 0, а вязкая жидкость — слой П2: 0 < у < 12. Система находится в покое, и при Ь = 0 во всем пространстве мгновенно создаются поля температур 9у = Лу х и концентрации с1 = Вх, где Лу ,В — постоянные, ] = 1, 2. Термоконцентрационный эффект вызывает однонаправленное движение системы, в которой поверхностью раздела останется плоскость у = 0. В такой интерпретации решение задачи будем искать в виде

Щ = Щ (У,Ь), 9 у = Л у х + Ту (у,Ь), с1 = Вх + К1(у,Ь).

Функции щ , Ту, К1 можно назвать возмущениями состояния покоя смеси и жидкости.

В данном предположении система уравнений термодиффузионного движения имеет следующий вид:

щг = ц Щуу,

ТИ = ХзТууу - ЛЩ,

Кц =

при -11 < у < 0 (] = 1), 0 < у < 12 (] = 2). Здесь V — коэффициент кинематической вязкости, Ху — коэффициент температуропроводности, ( 1 — коэффициент диффузии, а1 — параметр термодиффузии, индекс у означает дифференцирование по поперечной координате у.

На поверхности раздела выполнены условие равенство скоростей и динамическое условие

щ

где р1, р2 — плотности в слое П1, П2 соответственно. Условие непрерывности температур и равенства потоков тепла имеет вид

№у

здесь к1,к2 — коэффициенты теплопроводности в слое , П2 соответственно. Ещё одно условие — отсутствие потока вещества через границу раздела

К1у

В равенстве (2) Л = Л1 = Л2 (следствие равенства температур при у = 0).

К условиям (4)-(6) необходимо добавить условия на твёрдых стенках у = -1\, у = /2:

щ{-1ъг) = 0, И2

условия прилипания,

г1

задание нулевых возмущений температуры,

=0

у=-к

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

щ

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

Уравнения (1)-(3) с граничными условиями (4)-(10) образуют три последовательно решаемые задачи для функций (щ, «2), (Т1,Т2), К1.

2. Стационарное течение

Предположим, что все параметры, характеризующие течение в слоях, не зависят от времени, тогда начальные условия (10) не учитываются. В этом случае краевая задача (1)—(9) значительно упрощается. Её решение показывает, что в системе реализуется течение типа течения Куэтта и распределение температуры в слоях описывается полиномами третьего порядка

Щ = + у), «0 = - £),

Т0(у) = - (£ + у2) + а1у + «2,

где

Х1 \611 2

Т0

Х2 \ 6/2 2 /

Н/

а =--т-тг", Н = -(^А + Ж2 В1),

+ 1)

аА^/2 - х) аА/^к/ + х)

а1 = ~т;—71->\—, а2 —

3Х1(к + /) ' 2 3х1(к + 1) постоянные, найденные из граничных условий. Здесь введены обозначения / = /1//2, х

Х1/Х2, к = ^1/^2, ^ = ^1/^2.

Интегрирование уравнения для определения концентрации даёт зависимость

(В - ОА) а (£ + £) + * + Ь2.

а б в

Рис. 1. Стационарные профили скоростей

Из граничного условия на поверхности раздела определим, что Ъ1 = — а\а\. Граничное условие на твёрдой стенке выполняется только при В = 0. Значит, стационарное распределение концентрации в слое имеет место только при отсутствии её градиента в начальный момент времени в направлении движения. При В = 0 распределение концентрации может быть только нестационарным. Постоянную Ъ2 определим, если зададим среднюю концентрацию в сечении х = 0

о

У ко Шу = 0.

-1г

Таким образом, стационарное распределение концентрации (В = 0) в слое описывается многочленом третьего порядка

о

^о «1А / У3 , У2^ , агаА12 [ 1 I2 — х ,

К = — НАЧ У1 + -2) — а1а1У + "ХТ" 11— «№+) 1. (13

Для построения стационарных профилей скоростей, температур и концентрации перейдём к безразмерным переменным, определяемым формулами

о

« = У • й' = — Т • Т = К = — ОААк'0

Стационарное распределение скоростей, температур и концентрации представлено на рис. 1. Для расчётов была выбрана система этанол — вода и бензол с физическими свойствами х = 1.104, I = 1, V = 2.965, р = 1.104, а1 = 0.26, к = 3.082.

3. Решение нестационарной задачи

Для исследования нестационарных процессов при решении линейных начально-краевых задач в математике, механике и технике хорошо зарекомендовали себя операционные

методы на основе интегрального преобразования Лапласа, которое определяется формулой

сю

f

0

где f (p) — изображение функции f (t), f (t) — оригинал для f (p). Как правило, перевод в область изображений производится по переменной времени.

Основное свойство преобразования Лапласа, используемое далее, заключается в следующем предельном соотношении.

Теорема. Если f (t) является оригиналом и f (p) — изображение f (t), то

lim pf (p) = f (0),

p^-c

г^е p —^ то внутри угла | argp| < n/2 — # и f (0) = ^lim^ f (t); если, кроме того, существует lim f(t) = f(то), то

lim pf

г^е p — 0 внутри того же угла.

Доказательство утверждения и условие применимости формулы (15) описано в [5].

3.1. Определение нестационарного поля скоростей

В результате применения преобразования Лапласа к равенствам (1), (4), (7) приходим

к краевой задаче для изображений скорости Uj

p

u/--ut = 0, —Zt < y < 0, Ui(—Zi,p) = 0,

vi

p

Uo--U2 = 0, 0 < y < I2; U2(Z2,p) = 0,

H

U"

p

Здесь и далее штрих означает дифференцирование по у. Решение задачи (16) находится в виде

л/^ H rp>~

Ui(y,p) =--ATw/ 2 , sh\/ V7 (li + y), —li <y< 0,

ß2 Vp3 Wi (p) ch vV- i11 V vi 4 Vv2 H t^V/pVrT li 1 IT n Ч П/

u2(y,p) =--^ ^v v (l2 — У), 0 < у < l2,

ß2 VT3 Wi (p) s^pv2 T l2 V v2

где

Wi(p) = -ß= + th A /^li cth A /pz2. V v V Vi V V2

С учётом асимптотических представлений shx ~ x + x3/6, ch ж ~ 1 + x2/2 при ж — 0 и равенства pv = ß получим

lim pui(y,p) = — - T

p^0 ' ß2

1+y

= ui

Ит рЙ2(у,р) =--,Я11, л

Р^0 + 1)

1 - У

¿2

= «2 (у).

:18)

Таким образом, согласно теореме 1, с ростом времени решение задачи (16) выходит на стационарный режим и0(у), и°(у), определяемый равенством (11).

3.2. Определение температурных возмущений

Применение преобразования Лапласа к формулам (2), (5), (8) приводит к краевой задаче для изображений возмущений температур

Т" - — Т Х1

р ~ А«1

х1

, —¿1 <у< 0,

гТ// Р Т _

Т 2--Т 2 =

2 х2

А«2(у,Р) Х2

, 0 < у < ¿2,

^(0,р) = ТТ2(0,р), кТ (0,р) = Т2 (0,р), Т1(-/1,р) = 0, Т2 (¿2, р) = 0. Решение задачи (19) определяется в виде

Т(у,р) = ¿1 бЬ,/ —у + ¿2 сЬА/ — у+ х1 х1

49)

А

Vрх-1

х

и1(г,р) йЬ

-¿1

— (у - г) х1

Тг(у,р) = ¿3 вЬА/ —у + ¿4 сЬА/ — у+ х2 х2

А

Х2\[рХ

-1 2

и2 (г,р) йЬ

'— (у - г) х2

Коэффициенты = ¿¿(р), г = 1, 4, находятся из граничных условий и задаются формулами

= с'(р)- в2<р>, = и й, Мь, ь3 = -к=ь1 - сь

ИЪ(р) V Х1 -X

¿4 — ¿2 -

-, [ и1(г,р) бЬ , /-^г^г.

х^рх- I Vх1

-¿1

В формулах (22) введены обозначения

кА / х2 (* I р

в1 (р) =--4 — и1(г,р) сЬ , — г^г,

х1 V р J V х1

-¿1

в2(р) = - А^^¿2 /«1(г,р) бЬ

х^л/рхТ1 I V х1

У

У

о

A

X2V/pX-T sWpx-1 l

u2(z,p) sh

p

— (l2 — z) X2

dz,

W2(p) = + th, /^li cth, /^¿2. л/х V xi V X2

Покажем, что lim Tj

t

p—0

В качестве примера приведём расчёты для ^ = 1. С этой целью, используя (22), несколько преобразуем выражение (20):

Т1(У,р = , ^ 1—/—г, — (У +11)+ ^(р) с^^РХТ11г УХ1

A

Xiv7^"1p

ut(z,p) sh

-¿i

'—

Xi

dz.

Далее, учитывая равенства (17), (18), (23) при p — 0 (shx ~ x, chx ~ 1, x — 0) получим

W2(p) ~ k +l, pGi(p) ~ — kA

pG2

л/Х A

XiVX- tp

u0(z) dz,

-ii

XiWpx-1

0 ¿2 — J^ u0(z)zdz + xJ u0(z)(l2 — z) dz -i1 0

Подстановка эквивалентных выражений в (24) даёт

f ( )= AHlil2(3fcl + lil + 2X)(y + li) AHli Г + y2 + yli + ¿t p— i(y,p) 6ß2Xi(ß + l)(k +1) ß2 Xi(ß + ZK6li 2 2 6

что совпадает с TT0(y) из формулы (12).

Аналогично показывается, что limpT2(y,p) = T2 (y). Таким образом, температурные

p—^0

возмущения с ростом времени выходят на стационарный режим, определяемый формулами (12).

2

y

0

3.3. Определение возмущения концентрации

При анализе задачи для определения возмущения концентрации (3), (6), (9) удобно ввести замену Р = К1 + а1Т1. Тогда начально-краевая задача преобразуется к виду

Р = ^ + а Туу - Вм1,

Ру (0,*) = 0, Ру (-/ь*) = 0,

р (у, 0) = 0.

Решение задачи в изображениях по Лапласу будет следующим:

y

~ /Р /Р d i /

F(y,p) = D 1 sh w — y + D2 ch w — y +W— h(z,p) sh

V d 1 V d 1 V p j

' -ii

'dt

dz,

где

h = -apTi + 1, Dj = Dj(p), j = 1, 2. d d

После подстановки (25) в граничные условия находим

D1(p) = — /— f h(z,p) ch,/-^zdz, pd1 ' -ii

I— _ о _

d1 p p D2(p) = — W — cthw — l1 I h(z,p) chw — zdz. p d1 d1 ' -ii

Далее для исследования функции F(y,p) рассмотрим limp_P1(y,p). Для этого перепишем (25) в виде

~ [Z СЧ di (y +11) о rf

F(y, p) = — W----=- h(z,p) ch л — z dz+

Vp sh, /dil i Vd1

p 1 -ii

dp i h(z,p) sh

'dt

dz.

Учитывая limpu 1

lim pF1

p^0

a 1aAl2 a 1a1l 1

"8X7

,

+ a1a2, если В = 0, если B = 0.

Постоянные a,a1,a2 определены выше.

Так как K1 = F — a1T1, то получим, что с ростом времени возмущение концентрации выходит на стационарный режим при условии отсутствия её градиента в начальный момент времени, т.е. limpKK1 (y,p) = K°(y), где K°(y) определяется формулой (13).

4. Результаты численных расчётов

Полученные формулы в изображениях по Лапласу были использованы для численного нахождения полей скорости, температур и концентрации. Численные расчёты проводились в безразмерных переменных (14). При этом безразмерное время вводилось

соотношением т = -ртс.

¿2

-ii

Рис. 2. Нестационарные профили скорости (а), температуры (б) и концентрации (в) в разные моменты времени

Результаты численного обращения преобразования Лапласа, приведённые на рис. 2, подтверждают выход решения рассматриваемой задачи на стационарный режим. При этом время выхода для скорости, температуры и концентрации различно. Для скорости оно составляет т3 = 8, для температуры — т3 = 68, для концентрации — т3 = 1000.

Автор выражает благодарность д-ру физ.-мат. наук, профессору В.К. Андрееву за помощь в постановке задачи и постоянное внимание к работе.

Список литературы

[1] Андреев В.К., Захватаев В.Е., Рявицкий Е.А. Термокапиллярная неустойчивость. Новосибирск: Наука, 2000. 280 с.

[2] Андреев В.К. Эволюция совместного движения двух вязких теплопроводных жидкостей в плоском слое под действием нестационарного градиента давления // ПМТФ. 2008. T. 49, № 4. C. 94-107.

[3] Совачкина Н.Л. О совместном движении бинарной смеси и вязкой жидкости в теплоизолированной цилиндрической трубе // Вычисл. технологии. 2011. T. 16, № 4. C. 120-133.

[4] Лемешкова Е.Н. Прямая и обратная задача о совместном движении трёх вязких жидкостей в плоских слоях // Журн. СФУ. Математика и физика. 2011. T. 4, № 3. C. 363-370.

[5] Лаврентьев М.А., Шават Б.В. Методы теории функций комплексного переменного. М.: Наука, 1973. 736 с.

Поступила в 'редакцию 25 февраля 2013 г., с доработки — 22 апреля 2013 г.