МОДЕЛЬ МЕЖФАЗНОЙ ПОВЕРХНОСТИ

Уравнения (24) и (25) являются общими для любой плоской межфазной области между двумя несжимаемыми ньютонов­скими жидкостями, если только к описанию локальных свойств применимо условие аксиальной симметрии. В частности, они тогда должны быть справедливыми для любимой модели, исполь­зуемой специалистами по макроскопической динамике жидко­сти, а именно: бесконечно тонкая межфазная поверхность. В такой поверхности локальные свойства меняются скачком при пересечении плоскости 2 = 0. Запишем в явном виде:

Ц = и" = и’ = 4" + и11) + 4- * W (ип — и[2])

О 1 (ду1х дК1 1 (dVx Wl Г

Х N о

О

О


Здесь h (z) — функция Хэвисайда: h (z) = —1 для z < 0 и h (z) = 1 для z > 0. Из уравнений (26), а также рис. 3 сле­дует, что в этой модели межфазной поверхности плотность, коэффициенты вязкости и статический тензор давлений раз­рывны. Составляющие вектора скорости непрерывны, но произ­водные от тангенциальных составляющих dVJdz и ‘ dVJdz Разрывны, в то время как для нормальной составляющей раз­рывна вторая производная d2VJdz2. Наконец, величина 0° есть скорость дилатации поверхности разрыва

< Дх

К

Ду

00 =

= vs-v°

И условие несжимаемости (9) выполнено.

Подстановкой уравнений (26) в (24) и (25) все интегриро­вания для такой модели межфазной поверхности могут быть формально закончены (интересующийся читатель может про­следить за одним из этих интегрирований в Приложении II, Но нас не интересуют явные результаты). Для наших целей достаточно допущения, что все эти интегралы имеют поря-

Док 0 (Az) или меньше, т. е. они являются произведениями таких величин объемных фаз, как р1, р11 и т. д., на интервал Az. Наша цель при построении модели межфазной поверхности состояла в сравнении положения этой модели с реальной меж­фазной поверхностью. Этого мы достигаем, прибавляя и вычи­тая результаты для модельных интегралов к уравнениям (24) и (25).

В качестве примера одной из этих операций рассмотрим интеграл j Ptdz, который появляется в уравнениях (25). За­пишем

Dz-f

+I [~Т(pI + рП) + ~Тh (z) (рП ~ pl)] ^ = — V + о (Дz)

У = J [4 + рП) + 4" Л (Z) {PU ~ Pl)Pt Dz

IЛ л = | [Л — 4" ^ + PU)—Т h <г) {рП — р1)

Где

(27)

Является, очевидно, избыточной капиллярной величиной, тесно связанной с поверхностным натяжением, определенным урав­нением Баккера (5) для статической межфазной поверхности.

В действительности я использовал для обеих величин одно и то же обозначение у, но вопрос о согласованности уравне­ний (5) и (27) рассмотрю в следующем разделе.

Аналогичным образом можно привести интегралы j рVxdz И j р‘Vxdz к следующим:

J 9Vx Dz = J V® [ — j — (p1 + p11) + ~R h & (p" ~ рг) ] } dz + 0

(28)

+ 0(Az) (29)

Дальнейшие упрощения я обосновываю доводами, что поле скорости в экспериментальной системе очень слабо меняется на толщине Az рассматриваемой нами межфазной поверхности. Конечно, производные dVJdz и DVy!Dz значительно меняются на этом расстоянии, но я тем не менее предполагаю, что тол­щина Az может — быть выбрана столь малой, что в хорошем при-

Ближении V V1 ^ Vй везде внутри межфазной области. Если бы это было не так, экспериментатор наблюдал бы столь большие изменения в поле скорости по обе стороны межфазной поверхности, что был бы вынужден допустить существование

МОДЕЛЬ МЕЖФАЗНОЙ ПОВЕРХНОСТИ

А б

Рис. 4. Сравнение диффузной экспериментальной межфаз ной поверхности с моделями, приведенными на рис. с

А — плотность р монотонно меняется в анизотропной облает» избыточная капиллярная плотность массы Г, определенна уравнением (30), является суммой заштрихованных площаде} подходящим выбором разделяющей поверхности г = О, можн обратить Г в нуль;

Б — сдвиговая вязкость изображена с резко выраженным мак симумом: поверхностная сдвиговая вязкость ц’, определенна уравнением (31), равна заштрихованной площади; для приведер ной разделяющей поверхности 2=0, вязкость г) положительш но т| может быть сделана равной нулю достаточно большим см( щением плоскости г — 0 вправо.

В его системе плоскости скольжения, совпадающей с оптиче ски наблюдаемой межфазной поверхностью. Эксперименталь ных свидетельств о существовании такой плоскости сколь жения нет; поэтому-то я и выбираю величину в уравнения: (28) и (29), являющуюся ^-составляющей вектора скорост] межфазной поверхности, и полагаю Vx = Vx во всей облает] интегрирования. Вынося V°x из-под знака интеграла, получае1 из уравнений (28) и (29)

JPVxdz= rvj-f 0(Az);

J n‘Vxdz = r)V°x-t-0(Az)

Где

Г = J [р———- g — (р1 + р11)——— LTh(z) (р11 р1)] dz (30)

Dz (31)

11=1 К -4- W+4л (2) ^ ^

Избыточные капиллярные плотность массы и сдвиговая вяз­кость.

Уравнения (30) и (31) — основные в том классе явлений капиллярной реологии, рассмотрение которых является целью этой статьи. Возможные формы подынтегральных выражений изображены на рис. 4. Без помощи строгого статистико-механи — ческого рассмотрения установить аналитический вид подын­тегральных выражений невозможно, но выбранная форма, по крайней мере, согласуется с известными из опыта свойствами реальных жидкостей. Мы вернемся к этому вопросу в следую­щем разделе. Сейчас же достаточно отождествить интегралы (30) и (31) с заштрихованной площадью на рис. 4. Очевидно также, что значения этих интегралов зависят от положения раз­деляющей поверхности Z 0, которая в нашей модели является поверхностью разрыва функции Хэвисайда H (Z).

Последовательное приложение этих методов к уравнениям баланса импульса (24) и (25) позволяет привести их к виду:

Нормальная составляющая —

П

О = | Ргг + 0 (AZ) (32)

I

Тангёнциальные составляющие —

= | + "Ж+ * IF + IVK + 0 (А.)

(33)

II

— TV0 = Dt VY

I

Где у, Г и г] определены уравнениями (27), (30) и (31), а

* = ] И- —I — Л (г) (fx11 — f^1)] ^ (34)

Избыточная капиллярная дилатационная вязкость.

Наряду с у, Г и ц величина K не определена точно, пока не сделан выбор положения разделяющей поверхности Z = 0.

Комментирование и размещение ссылок запрещено.

Комментарии закрыты.