Scientific journal
International Journal of Applied and fundamental research
ISSN 1996-3955
ИФ РИНЦ = 0,593

THE RELAXING FLOW OF A VISCOUS LIQUID BETWEEN SOLID WALLS

Sennitskiy V.L. 1, 2
1 Lavrentiev Institute of Hydrodynamics SB RAS
2 Novosibursk State Universuty
The problem on a plane flow of a viscous incompressible liquid between absolutely solid walls is stated and solved. Force interactions supporting the liquid motion are absent; the liquid motion is relaxing. The exact solution of the problem is found. Integral parameters of the liquid flow are determined. The asymptotic formulas are obtained which characterize the hydro-mechanical system under consideration at large time. Problems on plane flows of a viscous liquid in the presence of solid walls are invariably in the number of actual problems of hydro-mechanics. The results of investigations of plane flows are of an independent interest and also can serve as important orienting points under the study of spatial flows. In particular in this connection it is actual an effective search of new adequate hydro-mechanical problems. One of perspective directions of viscous fluid mechanics involving such search consists in a study of the evolution of hydro-mechanical systems under a cessation of force influences supporting their motion. An especial role in such investigations belongs to exact solutions finding of the problems of viscous fluid mechanics. The actuality of works in this scientific field in particular is caused by their obvious applied significance and the presence of their connection with a wide spectrum of phenomena and processes.
asymptotic formulas
remaining mass of liquid
solid walls
stationary periodical relaxing flows,

Выявлению закономерностей движения жидкости в присутствии твердых стенок посвящено значительное число исследований (см., например, [1, 2], а также [3–5], и представленную там литературу).

В данной работе рассматривается гидромеханическая система, состоящая из вязкой несжимаемой жидкости и абсолютно твердых стенок Ξ1, Ξ2. Стенки покоятся относительно инерциальной прямоугольной системы координат XYZ. Стенка Ξ1 ограничена плоскостью Γ1: Y = 0; стенка Ξ2 – плоскостью Γ2: Y = H (H > 0 – постоянная). Жидкость заполняет область Ω: 0 < Y < H. В начальный момент времени t, при t = 0, течение жидкости является симметричным относительно плоскости Y = H / 2; жидкость движется со скоростью V0 = {V0(Y), 0, 0} (V0(Y) = V0(H – Y); V0 = 0 на Γ1, Γ2).

Цель данной работы состоит в определении скорости жидкости V = {V(Y, t), 0, 0}, а также силового взаимодействия жидкости и стенок, средней скорости жидкости, массопереноса жидкости – при t > 0.

Общая задача

Пусть y = Y / H; ssen01.wmf,
если ssen02.wmf, ssen03.wmf sup│V0│, если ssen04.wmf; τ = ssen05.wmft / H; v0 = V0/ssen06.wmf; v = V/ssen07.wmf;
ν – кинематический коэффициент вязкости жидкости; Re = H ssen08.wmf – число Рейнольдса.

Уравнение Навье – Стокса и условия, которые должны выполняться на границах стенок и в начальный момент времени, имеют следующий вид:

ssen09.wmf; (1)

v = 0 при y = 0; (2)

v = 0 при y = 1; (3)

v = v0 при τ = 0. (4)

Отметим, что при t > 0 какие-либо воздействия на жидкость, поддерживающие ее движение, отсутствуют.

Наряду с задачей (1)–(4) будем рассматривать также вспомогательную задачу

ssen10.wmf; (5)

w = 0 при y = 0; (6)

w = 0 при y = 2; (7)

w = w0 при τ = 0. (8)

Здесь w0 – периодическая с периодом 2 функция y, которая при 0 ≤ y ≤ 2 определяется формулой

ssen11.wmf (9)

Обратимся к задаче (5)–(7). Применяя метод разделения переменных, найдем, что уравнение (5) имеет последовательность решений

ssen12.wmf (10)

каждое из которых удовлетворяет условиям (6), (7). Используя (10), построим следующее решение задачи (5)–(7):

ssen13.wmf (11)

где cm (m = 1, 2, …) – постоянные.

Рассмотрим условие (8). Представим w0(y) в виде ряда Фурье. С учетом (9) и соотношения

v0(y) = v0(1 – y) (12)

получим

ssen14.wmf (13)

Здесь

ssen15.wmf

Используя (11), (13) и формулу

ssen16.wmf

найдем

cm = ssen17.wmf(m = 1, 2, … ; n = 1, 2, …). (14)

Из (11), (14) следует, что w (решение задачи (5) – (8)) удовлетворяет условию

w = 0 при y = 1. (15)

Таким образом, согласно (1)–(9), (11), (14), (15) задача (1)–(4) имеет решение

ssen18.wmf (16)

Формулой (16), в частности, демонстрируется, что изучаемое течение жидкости является затухающим.

Пусть ρ – плотность жидкости; γ1 – площадка площадью Sγ, принадлежащая плоскости Γ1; γ2 – площадка площадью Sγ, принадлежащая плоскости Γ2; Flw1 = {Flw1, 0, 0} – тангенциальная сила, действующая со стороны жидкости на часть ξ1 стенки Ξ1, граничащую с жидкостью на площадке γ1 ssen19.wmf Fwl1 = {Fwl1, 0, 0} – тангенциальная сила, действующая со стороны части ξ1 стенки Ξ1 на жидкость; Flw2 = {Flw2, 0, 0} – тангенциальная сила, действующая со стороны жидкости на часть ξ2 стенки Ξ2, граничащую с жидкостью на площадке γ2 ssen20.wmf; Fwl2 = {Fwl2, 0, 0} – тангенциальная сила, действующая со стороны части ξ2 стенки Ξ2 на жидкость; flwk = HFlwk / (ρssen21.wmf ,
fwlk = HFwlk / (ρssen22.wmf (k = 1, 2); σ – площадка: X = X*, 0 ≤ Y ≤ H , – Z*/2 ≤ Z ≤ Z*/2 (X*, Z* > 0 – постоянные); Sσ = H Z*;

ssen23.wmf (17)

– среднее значение скорости V по координате Y; ssen24.wmf

Используя (16), (17), получим

ssen25.wmf (18)

(fwl1 = – flw1, fwl2 = – flw2);

ssen26.wmf (19)

Формулой (18) определяется силовое взаимодействие жидкости и стенок (вопрос о нормальном силовом взаимодействии жидкости и стенок является тривиальным).

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

ssen27.wmf (20)

– масса жидкости, которая протекает через площадку σ из области X < X* в область X > X* за промежуток времени 0÷t (eX = {1, 0, 0}) . Используя (17), (19), (20), найдем

ssen28.wmf (21)

Из (16), (18), (19), (21) следуют асимптотические формулы, характеризующие рассматриваемую гидромеханическую систему на больших временах

ssen29.wmf при τ → ∞; (22)

ssen30.wmf при τ → ∞; (23)

ssen31.wmf при τ → ∞; (24)

ssen32.wmf при τ → ∞. (25)

Здесь

ssen33.wmf (26)

– полная остаточная масса жидкости (предел Μ при t → ∞). Отметим, что ряд в (26) является абсолютно сходящимся (данный ряд мажорируется сходящимся рядом ssen34.wmf.

Первая частная задача

Пусть жидкость совершает движение со скоростью U = {U, 0, 0}, не изменяющейся со временем; задача о течении жидкости имеет вид

ssen35.wmf (27)

U = 0 при Y = 0; (28)

U = 0 при Y = H. (29)

Здесь η = – (1/ρ) ∂P /∂X – постоянная (P – давление в жидкости; без умаления общности может быть принято, что η > 0). Использование (27)–(29) приводит к формуле

ssen36.wmf (30)

Положим

V0 = U. (31)

Отметим, что

ssen37.wmf

удовлетворяет соотношению (12).

Выполнение (31) соответствует тому, что задачей (1)–(4), определяемым ею течением жидкости моделируется происходящее при t > 0 остаточное, не поддерживаемое силовыми воздействиями затухающее течение вязкой жидкости, совершающей при t ≤ 0 движение с не изменяющейся со временем скоростью

ssen38.wmf

С учетом (16), (17), (19)–(21), (26), (30), (31) для данного (моделирующего) течения жидкости, в частности, имеем

ssen39.wmf

ssen40.wmf

Здесь

ssen41.wmf, ssen42.wmf

Вторая частная задача

Пусть жидкость совершает движение со скоростью U2 = {U2, 0, 0}, периодически изменяющейся со временем; задача о течении жидкости имеет вид

ssen43.wmf (32)

U = 0 при Y = 0; (33)

U = 0 при Y = H. (34)

Здесь η = ξ [1 + sin(2πt / T + φ)] (ξ, T > 0. 0 ≤ φ < 2π – постоянные; без умаления общности может быть принято, что ξ > 0). Использование (32)–(34) приводит к формуле

ssen44.wmf (35)

где χ = {sh[(1 + i)ϰY] + sh[(1 + i)ϰ(H – Y)] – sh[(1 + i)ϰH]} / sh[(1 + i)ϰH]; ssen45.wmf

Положим

ssen46.wmf (36)

Отметим, что

ssen47.wmf

удовлетворяет соотношению (12) .

Выполнение (36) соответствует тому, что задачей (1) – (4), определяемым ею течением жидкости моделируется происходящее при t > 0 остаточное, не поддерживаемое силовыми воздействиями затухающее течение вязкой жидкости, совершающей при t ≤ 0 движение с периодически изменяющейся со временем скоростью

ssen48.wmf

С учетом (16), (17), (19)–(21), (26), (35), (36) для данного (моделирующего) течения жидкости, в частности, имеем

ssen49.wmf

ssen50.wmf

Здесь

ssen51.wmf

ssen52.wmf (λ2n – 1 = 2H2 / [(2n – 1)2πνT]; 0 < φ2n – 1 < π /2 (n = 1, 2, …) – углы, удовлетворяющие соотношениям ssen53.wmf

Заключение

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