CC BY
22Вестник Сыктывкарского университета. С ер Л. Вып. 13.2011
УДК 519.6
ДИССИПАТИВНАЯ СТРУКТУРА И ОБЛАСТЬ СВЕРХАНОМАЛИИ КУЭТТОВСКОГО ТЕЧЕНИЯ СТРУКТУРИРОВАННОЙ ЖИДКОСТИ В ПЛОСКОМ
ЗАЗОРЕ 1
Н. А. Беляева, К. П. Кузнецов
Проведено бифуркационное исследование куэттовского течения структурированной жидкости в плоском зазоре в области сверханомалии. Получены бифуркационные диаграммы и определены значения параметров, соответствующих области сверханомалии. Бифуркационный метод позволил получить аналитическое приближение стационарного неоднородного решения в окрестности точки бифуркации. Проведено численное моделирование течения.
Ключевые слова: бифуркационный анализ, структурированная жидкость, параметрический анализ, область сверханомалии, ку-эттовское течение
1. Постановка задачи. Пусть структурированная вязкая несжимаемая жидкость заполняет полосу между плоскостями ( = 0 и ( = й, а её течение происходит в направлении оси г/ и вызывается движением плоскости £ = К. Обозначим через й - вектор скорости жидкости в произвольной точке полосы в момент времени тогда й — (0, г^, 0), причем иц — пБудем предполагать, что жидкость является смесью двух компонент А1 и А2, которые под действием приложенного механического поля, характеризуемого напряжением а и скоростью деформации
ди
1 Работа выполнена при финансовой поддержке Министерства образования и науки РФ ГК № 02.740.11.0618
© Беляева Н. А., Кузнецов К. П., 2011.
могут взаимно превращаться друг в друга. Суммарная скорость превращения определяется соотношением:
ф(а, 7) = к2( 1 - а) - kia, (0.1)
ki = к0 ехр (р0сг + §о72) •
Здесь а — степень структурных превращений (доля А1 в общей смеси), ко, qo~ параметры жидкости.
Будем предполагать, что жидкость подчиняется реологическому уравнению состояния вида
а = Д7, (0.2)
где вязкость ¡л = ц(а) задается соотношением
¡i~l(a) — ¡i[la + /^(1 — a)i
/ii,/i2 ~~ вязкости компонент А1 и А2, соответственно.
Модель течения жидкости описывается с помощью диффузионно-кинетического уравнения относительно степени структурных превращений а и уравнения движения [1,2]:
д2а
да =
dt эе
Ф(а,7)
du д ( . ,ди\
"lit = äe
(0.3)
р — плотность, И — коэффициент диффузии.
Граничные и начальные условия рассматриваемой задачи имеют вид:
да t = 0 : а\
да
Ж
= о, щ=о = 0;
Ü=h
= 0.
0 s^i<h
(0.4)
(0.5)
10^^/1 ~ а
В случае псевдопластической жидкости считаем, что подвижная плоскость движется с постоянной скоростью ^о:
и\^=н = > о. (0.6)
В случае дилатантной жидкости — на подвижной скорости задана постоянная скорость деформирования 70:
дп
= 7о-
(0.7)
i=h
2. Кривая течения. Рассмотрим однородные стационарные состояния системы (2)-(6), то есть решения вида:
a = const, 7 = const.
Тогда система (0.3) сводится к уравнению
Ф(а,7) =0. (0.8)
Из соотношений (0.2) и (0.8) заменой переменных [2]
к
у = сф0, х = рот, 5 = 9о(роМ2)~2, X = "Г (°'9)
получим уравнение кривой течения или реологической кривой, связывающей кинематическую и динамическую характеристики:
у = х( 1 + fc)(A + 1 + к), к = хехр(у + &Г2) (0.10)
На рис.1 изображены реологические кривые для различных значений параметра Л.
Рис. 1. Реологические кривые 0: Л = 0 — ньютоновская жидкость, 1-3: — 1 < Л < 0 — псевдопластическая жидкость, 4-6: Л > 0— дилатантная жидкость
Явление падения напряжения при увеличении скорости деформирования для псевдопластической жидкости (рис.1 кривая 3) и падение скорости деформации при увеличении напряжения для дилатант-ной жидкости (рис.1 кривая 6) названо сверханомалией вязкости. В
работе [2] проведено исследование поведения кривой (0.10) в зависимости от значений параметра х — Получены граничные значения X*, которым соответствуют кривые 2 и 5, разделяющие область монотонных кривых у(х) и х(у) и область кривых с двумя экстремумами. В настоящей работе анализируется поведение степени структурных превращений в области сверханомалии вязкости в зависимости от изменений некоторых характерных параметров модели.
3. Параметрический анализ. Исследование зависимости решений уравнений от параметров позволяет находить условия возникновения критических явлений: множественность стационарных состояний, потеря устойчивости этих состояний. Введем безразмерные переменные и параметры:
рО
'
рк2И2 №
(0.11)
V
и
Щ
для псевдопластической жидкости;
V =
и
То
—для дилатантнои.
Тогда из системы (0.3)-(0.7) получим безразмерную модель:
(0.12)
с начальными
Г = 0 : = а1о^1 = а0
и граничными условиями
Будем рассматривать стационарные решения модели. Тогда, выразив из первого уравнения градиент скорости
и подставив во второе, систему (0.12) приведем к следующему уравнению второго порядка:
Стационарное уравнение (0.14) запишем в виде нормальной системы:
Исследуем однородные решения системы (0.15) для псевдопластической жидкости. Произведем замену:
у =рС,5 = —.
р1
Считаем у и 8 варьируемыми в определенных пределах параметрами, остальные параметры фиксируем. Однородные решения системы (0.15) удовлетворяют уравнению
Ф(а, у, 8) = ~ [1 - а - аХехр (у + 8у2 (1 + Аа)2)] = 0. (0.16)
Нелинейное уравнение (0.16) может иметь несколько решений, что приводит к множественности однородных состояний: при значениях параметров у и 8 из области сверханомалии существует [1] три решения степени структурирования а.
При варьировании одного из параметров существуют его особые значения — бифуркационные значения, при которых меняется число и устойчивость однородных состояний. Варьирование второго параметра приводит к тому, что на плоскости (у,8) бифуркационные значения описывают некоторые кривые, которые являются бифуркационными кривыми. На плоскости существует две основных бифуркационных кривых: кривая кратности однородных состояний 1/д и кривая нейтральности
6У_
Лг
с( 1 + Аа), с =
Си &
¡3— = -к [1 - о - ах ехр (рс + дс2 (1 + Аа)2)] . (0.14)
(Iг2
(0.15)
Сформируем матрицу Якоби:
з = (У «12 ) = (0.17)
ехр [у + 5у2 (1 + Аа)2] (Х + 2дХХ6у2 (1 + Ла)))
/5 0
Для построения линии кратности Ьд уравнение (0.16) необходимо дополнить условием А(7) = «11022 — О12О21 = 0. Тогда границы области множественности стационарных однородных состояний определяются решением системы:
(0.18)
[ Д(«7) = аца2 2 аХ2а2Х = 0.
Решив данную систему, получим выражение для кривой 1/д в параметрическом виде:
2Ха (а — 1)
о А — -2>
(1 + Ха) (2Ха2 1п - Ла - 2Ла 1п -2а2 -Ла-2Ла1п(^) -
Уа
2Ла (а - 1)
Кривая кратности делит параметрическую плоскость (у, 5) на области, в каждой из которых число стационарных однородных решений остается постоянным. Изменение значений параметров у и 6 с переходом через линию кратности приводит к изменению числа однородных решений на два. Точка на линии кратности соответствует точке поворота на диаграмме решений.
В явном виде получить решение а = а(у) уравнения (0.16) не удается. Однако можно получить функцию у (а) обратную искомой параметрической зависимости. Тогда график функции а (у) получается простой перестановкой координат у и а. На рис. 2 приведены несколько различных диаграмм решений с указанием точек поворота.
а
Рис. 2. Диаграмма стационарных однородных решений а(у). Сплошные линии — устойчивые решения, штриховые — неустойчивые решения, Ьр - точки поворота;
Л = -0.7, к = 0.1, х = 1е - 5, = 0.001
Исследуем устойчивость стационарных однородных решений уравнения (0.14). Для этого рассмотрим характеристический многочлен системы (0.15):
ср2 - (ап + а22)^ + (аца22 - а^г) (0.20)
Учитывая, что а = ац + а22 = 0, получим характеристическое уравнение:
+ Д = 0, (0.21)
корни которого определяют устойчивость решений. Область с тремя стационарными однородными решениями, одно из которых неустойчиво, характеризуется условием А < 0. При А > 0 уравнение (0.21) имеет два чисто мнимых сопряженных корня. В работе [2] показано, что существующее решение в данном случае устойчиво. Случай А = 0 определяет границу области устойчивых решений системы. Таким образом, система (0.18) задает не только множество точек поворота — кривая кратности £д), но и кривую точек комплексной бифуркации или бифуркации Хопфа — кривая нейтральности Ьа. Изменение параметров у и 6 с переходом через кривую нейтральности приводит к изменению устойчивости одного из имеющихся стационарных однородных решений задачи.
Анализ расположения линий 1/д, Ьа на плоскости (у, 5) позволяет определить параметрический портрет системы — рис. 3. В области,
отмеченной цифрами 1, существует одно устойчивое стационарное однородное решение. В области, отмеченной цифрой 3 и соответствующей области сверханомалии, существует три решения уравнения (0.14): а\ < а2 < аз. При этом среднее решение а2 является неустойчивым, а крайние а\ и аз устойчивы [1]. Несколько расчетных значений параметра у, задающих область сверханомалии при 6 = 0.2, приведены в таблице 1.
При 6 —>> 0 ветви кривой сближаются и исчезают — реологическая кривая течения в этом случае монотонна.
б
0.0-
10 У
Рис. 3. Бифуркационная диаграмма псевдопластической жидкости. Л = -0.7, к = 0.1, Х=1е-5,Р = 0.001
Схема исследования однородных решений системы (0.15) дилатант-ной жидкости аналогична исследованию псевдопластической жидкости, приведенному выше. Используя замены (0.9) и (0.11), получим следующую систему уравнений:
±А Лг
к
1
а — ах ехр
х
(1 + Л а)
+ 5х'
Аа Лг
Якобиан которой имеет вид:
1 ~
ах А
(1 + А а)
.1
0
(0.22)
(0.23)
Решая соответствующие системы уравнений получаем соотношения, параметрически задающие совпадающие кривые кратности и нейтральности. Кривая кратности (хд,5д) определяется соотношениями:
хА = --
1 + 2Ла + А 2 а2 (а — 1) Аа
(1п (а ~ Х)2Ла + (Ла + Х)(а ~ !)) Ла (1 + Аа)4
(0.24)
б
0.10'
0.05-
0.00-
—г-
12
16
20 у
Рис. 4. Бифуркационная диаграмма дилатантной жидкости. Цифрами 1,3 отмечено количество стационарных однородных решений. (А = 0.8, к = 0.1, х = 1е - 6, /? = 0.001)
В отличие от псевдопластической жидкости область множественности однородных состояний, изображенная на рис.4, замкнута в случае дилатантной жидкости: сохраняется при условии 5 = 0.
В работе [1] показано, что в области монотонности реологической кривой однородное решение устойчиво, а в области сверханомалии одно из решений (среднее) неустойчиво.
4. Бифуркационный анализ. Запишем уравнение (0.14) в следующем виде:
лч/ \
Рдг2=ф(а>У)>
(0.25)
Выберем в качестве бифуркационного параметра
7 дФ
где а0 — неустойчивое однородное стационарное решение; остальные параметры считаем фиксированными. Существуют критические значения ¿о, при которых возможны бифуркации неоднородных структур из однородного распределения а^. Находится из условия обращения в ноль собственного значения оператора задачи (0.25), линеаризованного на постоянном решении ао:
При значениях с1 близких к происходит рождение неоднородного стационарного решения.
Бифуркационный анализ состоит в представлении решения вблизи критического значения параметра в виде асимптотического ряда по степеням некоторого параметра е, характеризующего отклонение бифуркационного параметра от его критического значения, и нахождения коэффициентов этого разложения. Известный метод Ляпунова-Шмидта [3] определяет рекуррентную процедуру для получения коэффициентов разложения. На каждом шаге этой процедуры необходимо решать линейные уравнения с введенным выше оператором Ь. Правые части уравнений определяются разложением нелинейных членов уравнения (0.25) и зависят от решений, полученных на предыдущих шагах алгоритма. В соответствии со сказанным представим искомое стационарное решение уравнения (0.25) в виде:
Подставляем выражение (0.27) в стационарное уравнение, раскладывая нелинейные члены в ряд Тейлора в окрестности точки ао- В полученном разложении слагаемые при е,е2,ег приравниваем к нулю — получаем уравнения для определения коэффициентов а^, Таким образом, слагаемые первого порядка дают уравнение относительно а\\
дг2
(0.26)
Соответствующие значения
¿о = тг2п2/3, п — 1,2,...
а — а0 + ах{г)е + а2(г)е2 + а3(г)б:3 + 0{еА) в, = ¿о + ¿хе + (12е2 + 3 + 0(б4)
(0.27)
1/(а0, = 0.
Отсюда следует, что а\ является собственной функцией оператора L, соответствующей нулевому собственному значению. Эта функция имеет вид:
а\ — A cos (тгг), где А = const. Для остальных членов ряда получаем рекуррентные соотношения:
L(a0, d0)di = Фг(аг_ь = 2, 3, (0.28)
где функции Ф^ определяются из разложения нелинейных членов уравнения (0.25) в ряд Тейлора в окрестности а0. Так как оператор L(a0, dQ) имеет нулевое собственное значение, необходимым условием разрешимости этих уравнений относительно является ортогональность правых частей уравнения (0.28) к собственным функциям, соответствующим нулевому собственному значению оператора L. Из этого условия находится di-i:
) cos (тгг) Фг (аг_ь di-г) dr = 0 (0.29)
о
Уравнения (0.28) и (0.29) позволяют рекуррентно вычислить В
целях удобства описанная выше процедура была проведена с использованием математического пакета Maple. С использованием обозначения dl Ф
получены следующие результаты:
а=ао
dl = 0
. _ A2(5F2 - 3F37T2/3)
d2 ~ 24^5
ах = A cos (тгг)
, . A2F2(cos (2irr) — 3) a2 = В cos (тгг) + — v v ; ;
F =
da1
+
12тг 2/3
+ F37T'-192tt
32tt2/5F2A В (cos (2vrr) - 3)
, , A3 (F22 + F3Tr2f3) cos (З-лт) , a3 = С cos (тгг) +-^^-+
192тг4/32
В разложении (0.27) для разности (с? — do) первый отличный от нуля член содержит е2 и, следовательно, е = \Jid- ¿о) /с?2 + О (в). Поэтому бифуркация односторонняя: если с/2 > 0, то неоднородное решение вблизи точки бифуркации существует для значений d > do — надкритическая бифуркация; при < 0 неоднородное решение существует
для значений А < ¿о — бифуркация докритическая. В общем виде выражение для ¿2 громоздко и трудно поддается анализу.
Рис. 5. Диаграмма стационарных однородных решений псевдопластической жидкости; Вр —точки бифуркации; штриховая линия: ¿2 = 0;
Л = -0.7, к = 0.1, * = 1е - 4, = 0.001
Проведенный численный анализ показал, что для всех рассмотренных значений параметров бифуркации надкритические (рис. 5).
В работе [1] численно показано: в области сверханомалии в результате потери устойчивости однородного стационарного решения формируется устойчивое неоднородное решение — диссипативная структура, являющаяся монотонно возрастающей или монотонно убывающей функцией. Одной из задач данной работы являлось: проведение аналитических исследований на предмет определения вида формирующегося устойчивого неоднородного состояния. Сравнение численного решения задачи (0.12) с аналитическим решением на основе бифуркационного анализа, представленным в данной работе, подтверждает результаты численного анализа: качественно оба решения совпадают. На рисунке 6 представлен случай монотонно убывающей диссипативной структуры.
а
1.0-
0.5-
0.0
0.0
0.5
1.0
г
Рис. 6. График стационарных неоднородных решений задачи (0.12)-(0.13).
1 - численное решение, 2 - аналитическое решение (0.27) стационарной задачи
(Л = -0.7, ß = 0.001, к = 0.1, х = 1е - 5, у = 7, S = 0.2, р = 5.35, q = 5.73, а0 = 0.43, А — \,В — \,С — 500, г = 0.078)
Выполненный параметрический и бифуркационный анализ области немонотонности реологической кривой куэттовского течения структурированной жидкости, позволяет дать более широкую характеристику области сверханомалии: определить значения параметров ей соответствующих, охарактеризовать тип бифуркации в окрестности неустойчивого однородного состояния, дает аналитический метод исследования формирующегося устойчивого неоднородного состояния.
Литература
1. Беляева H.A. Неоднородное течение структурированной жидкости // Математическое моделирование. 2006. Т. 18. № 6. — С. 3-14-
2. Худяев С.И. Пространственная неоднородность и автоколебания при течении структурированной жидкости // Математическое моделирование. 2002. Т. Ц. № 7.- С. 53-73.
3. Иванова А.Н., Тарнопольский B.J1. О бифуркации неоднородных стационарных состояний в задаче о фазовом разделении бинарной полимерной смеси // Журнал вычислительной математики и математической физики Т.49. № 1999. С. 653-659.
4. Холодниок М., Клич А., Кубичек М., Марек М. Методы анализа нелинейных динамических моделей. М.: Мир, 1991. 368с.
Summary
Belyaeva N. A., Kuznetsov К. P. The dissipative structure and domain of anomaly structural liquid Couette flow in a flat clearance
The bifurcation study of structural liquid Couette flow in a flat clearance in the superanomaly area was conducted. Bifurcation diagrams and the values of parameters corresponding to the superanomaly area were obtained. Bifurcation method allowed to obtain an analytical approximation of the stationary inhomogeneous solution in the neigborhood of the bifurcation point. A numerical simulation of the flow was conducted.
Keywords: bifurcation analysis, couette flow, structural liquid, superanomaly area, parametrical analysis
Сыктывкарский государственный университет Поступила 01.04-2011
