Scientific journal
Modern high technologies
ISSN 1812-7320
"Перечень" ВАК
ИФ РИНЦ = 0,940

ON A MODEL OF THE MOTION OF A VISCOELASTIC FLUID WITH WALL SLIP

Baranovskiy E.S. 1 Artemov M.A. 1
1 Voronezh State University
This paper is concerned with a mathematical model describing unsteady flows of an incompressible viscoelastic fluid of the Oldroyd kind in a three-dimensional bounded domain with solid impermeable boundaries. In the constitutive law, we use the total derivative. An important feature of the problem under consideration is that Navier’s slip law (the free-slip condition) at the solid surface is used instead of the classical no-slip condition. This law gives a relation between the slip velocity and the force acting on the fluid from the interphase boundary of liquid-solid. For the corresponding evolution problem we introduce the concept of generalized solutions in a class of suitable function spaces. The paper discusses an algorithm for finding solutions by the Faedo–Galerkin method. Our main result is an analogue of the Sather–Serrin theorem on the uniqueness, regularity and energy estimates of global generalized solutions.
mathematical models
non-Newtonian fluids
the Oldroyd model
the effect of slip
the Faedo–Galerkin method
global solutions
regular solutions
energy estimates

Как известно, многие материалы, используемые на практике, при определенных условиях проявляют текучесть, для которой характерно отклонение от закона трения Ньютона [1]. Для решения задач о течении неньютоновских жидкостей успешно применяется метод математического моделирования, который предполагает построение адекватной математической модели течения и изучение её, например, с помощью реализуемых на компьютерах вычислительных алгоритмов. При этом требуется тщательное обоснование используемых алгоритмов. Такое обоснование тесно связано с изучением качественных свойств модели, в том числе c рассмотрением вопросов о существовании, единственности и регулярности решений соответствующей гидродинамической задачи.

В настоящей статье исследуется одна математическая модель, описывающая движение вязкоупругой жидкости типа Олдройда [11] в трёхмерной ограниченной области с непроницаемыми стенками. В уравнении состояния Олдройда используется полная производная. Важной особенностью рассматриваемой модели является то, что на границе области течения вместо классического условия прилипания задано условие проскальзывания Навье (обзор основных краевых условий, которые используются в механике жидкостей, приводится в [5]). Это условие определяет связь между скоростью скольжения жидкости и силой, действующей на жидкость со стороны межфазной границы «жидкость – твердое тело».

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

Основной результат работы представлен в разделе 3, где сформулирована и доказана теорема о единственности, регулярности и энергетической оценке глобального, т.е. заданного на заранее выбранном промежутке времени [0, T], обобщенного решения. Эта теорема обобщает известные результаты Дж. Сэзера и Дж. Серрина [12] для трёхмерных эволюционных уравнений Навье–Стокса. Как и в случае системы Навье–Стокса, единственность решения начально-краевой задачи для модели Олдройда удается установить в классе, более узком, чем пространство, для которого известно существование решений. Тем не менее, полученные результаты демонстрируют связь между регулярностью и единственностью, а также дают возможность определить энергетические оценки решений.

Модель движения вязкоупругой жидкости в ограниченной области

Рассматривается математическая модель, описывающая движение несжимаемой вязкоупругой среды типа Олдройда в ограниченной трёхмерной области Ω с неподвижной непроницаемой границей Г класса С2 при краевом условии проскальзывания Навье:

baran01.wmf (A)

где QT = Ω× (0, T) – пространственно-временной цилиндр c T > 0, ГT = Гт(0, T), v = v(x, t) – скорость, S = S(x, t) – упругая часть избыточных напряжений, p = p(x, t) – давление, E = E(v) – тензор скоростей деформации, E(v) = (∇v + (∇v)T)/2, f = f(x, t) – плотность внешних сил, D/Dt – полная (субстанциональная) производная,

baran02.wmf

η, λ и α – положительные постоянные (η – вязкость, λ – время релаксации, α = 1 – λ-1κ – – параметр ретардации, где κ – время запаздывания, 0 < κ < λ), n = n(x) – единичный вектор внешней нормали к Г, k = k(x) – коэффициент проскальзывания, k(x) > 0, [·]tan – касательная составляющая вектора. Символы v0 и S0 обозначают соответственно поле скоростей и упругую часть избыточных напряжений при t = 0. Более подробное описание указанных физических величин приводится в [1].

Рассмотренную в этом разделе начально-краевую задачу для краткости будем называть задачей (A). Неизвестными в системе (A) являются вектор-функции v, S и функция p, а все остальные величины считаются заданными.

Обобщённая постановка задачи

Ниже мы будем использовать пространства Лебега и Соболева. Определения и свойства этих пространств систематически изложены в [8]. Скалярное произведение в L2(Ω) и L2(Г)) будем обозначать через (·, ·) и (·, ·)Г соответственно. Положим,

baran03.wmf

baran04.wmf

baran05.wmf

baran06.wmf

Предполагая, что функция k принадлежит пространству L∞(Г) и k(x) > k0 > 0, введем скалярное произведение в X(Ω) по следующей формуле:

baran07.wmf

Из неравенства Корна следует, что норма, соответствующая этому скалярному произведению, эквивалентна норме, наследованной из пространства H1(Ω)

Следуя [6], будем отождествлять пространство H(Ω) c его сопряженным: H*(Ω) ≡ H(Ω). Поэтому имеем цепочку вложений: X(Ω)⊂H(Ω) ≡ H*(Ω)⊂X*(Ω).

Наконец, обозначим через baran08.wmf пространство симметрических 3/3-матриц.

Теперь мы можем сформулировать определение обобщенного решения задачи (A). Пусть

baran09.wmf

Определение. Обобщенным решением задачи (A) назовем пару вектор-функций

baran10.wmf baran11.wmf,

такую, что

baran12.wmf

baran13.wmf

выполнены начальные условия: v(·,0) = v0 и S(·,0) = S0, а также следующие равенства

baran14.wmf

baran15.wmf

baran16.wmf

baran17.wmf

в смысле распределений на (0, T).

Замечание 1. Определение обобщенного решения задачи (А) в целом соответствует восходящей к работам Ж. Лерэ концепции слабых, или турбулентных, решений уравнений Навье–Стокса, которые с предельной тщательностью обсуждаются в фундаментальных трактатах О.А. Ладыженской [4] и Р. Темама [6].

Замечание 2. Как обычно, при рассмотрении обобщенных решений гидродинамических уравнений из определения решения «выпадает» давление p, которое после нахождения остальных неизвестных определяется с помощью результатов де Рама (см., например, [6, гл. I]).

Замечание 3. Для построения обобщенного решения может быть использован метод Фаэдо–Галеркина. Приведем основную идею данного метода. В подходящем функциональном пространстве выбирается базис. Затем вводится приближенное решение, которое записывается в виде суммы первых m базисных функций с коэффициентами, зависящими от t. Использование приближенных решений позволяет свести исходную начально-краевую задачу к задаче Коши для системы обыкновенных дифференциальных уравнений, решение которой не представляет существенных сложностей. Используя энергетические оценки и теоремы о компактности, можно выделить из последовательности приближенных решений подпоследовательность, которая сходится при m> к решению исходной задачи. На основе приведенной выше схемы в [7, 9] строятся решения различных типов начально-краевых задач для модели Олдройда; ряд стационарных задач о течении жидкости Олдройда рассмотрен в работах [2, 3, 10].

Основной результат работы

Введем обозначения:

M1 = L8(0, T; L4(Ω)),

baran18.wmf

Теорема. Пусть пара вектор-функций (v, S) – обобщенное решение задачи (A), принадлежащее классу M = M1/M2. Тогда это решение (если оно существует) единственно в классе M, удовлетворяет включениям

baran19.wmf baran20.wmf

и энергетической оценке

baran21.wmf

baran22.wmf

baran23.wmf

baran24.wmf

Доказательство. Пусть (v1, S1) и (v2, S2) – обобщенные решения из класса M. Покажем, что v1 = v2 и S1 = S2.

Проверим сначала, что

baran25.wmf (1)

Введем операторы:

baran26.wmf

baran27.wmf

baran28.wmf

baran29.wmf

baran30.wmf

baran31.wmf

Из определения обобщенного решения следует, что

baran32.wmf

Очевидно, что имеют место включения:

baran33.wmf (2)

Покажем, что

baran34.wmf (3)

Используя неравенство Гёльдера, находим, что

baran35.wmf

baran36.wmf

Здесь и ниже через Сi, i = 1, 2,…, обозначаются константы. Из полученной оценки и условий теоремы следует включение (3), которое вместе с (2) влечет (1).

Докажем теперь, что

baran37.wmf (4)

Введем операторы:

baran38.wmf

baran39.wmf

baran40.wmf

baran41.wmf

baran42.wmf

baran43.wmf

Из определения обобщенного решения следует, что

baran44.wmf

Очевидно, что

baran45.wmf. (5)

Используя неравенство Гёльдера, получаем, что

baran46.wmf

baran47.wmf

откуда, с учетом условий теоремы, вытекает, что

baran48.wmf

Это включение вместе с (5) влечет выполнение (4).

Принимая во внимание (1), (4) и применяя лемму 1.2 из [6, гл. III], получаем, что

baran49.wmf

baran50.wmf

и п.в. на [0,T] выполнены равенства:

baran51.wmf (6)

baran52.wmf (7)

где u = v1 – v2 и F = S1 – S2. Так как (v1, S1) и (v2, S2) – обобщенные решения задачи (A), то нетрудно вывести равенства:

baran53.wmf

baran54.wmf

baran55.wmf

baran56.wmf

baran57.wmf

Полагая φ = u(t), Ф = (αη)-1F(t) и складывая полученные равенства, с учетом (6), (7) находим

baran58.wmf

baran59.wmf

baran60.wmf

baran61.wmf (8)

Применяя неравенства Гёльдера, Юнга и теоремы о вложении пространств Соболева, находим, что

baran62.wmf

baran63.wmf

baran64.wmf

baran65.wmf

baran66.wmf (9)

Нам потребуется также следующая оценка:

baran67.wmf

baran68.wmf

baran69.wmf

baran70.wmf. (10)

С учетом (9) и (10), мы выводим из равенства (8) следующую оценку:

baran71.wmf

baran72.wmf

c функцией

baran73.wmf

Поскольку u(0) = 0 и S(0) = 0, то, применяя лемму Гронуолла–Беллмана, получаем, что u = 0 и S = 0. Таким образом, мы установили, что в классе M не может быть более одного обобщенного решения задачи (A). Из приведенного выше доказательства видно, что если такое обобщенное решение существует, то оно принадлежит пространству C([0, T]; H(Ω))×C([0, T]; L2(Ω)). Энергетическая оценка обобщенного решения выводится аналогично оценке (8) с последующим интегрированием в пределах от 0 до t.

Работа выполнена при финансовой поддержке РФФИ, грант № 16-31-00182 мол_а.