Добавил:
tg: @Yr66gi4 Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

экзамен / Билеты по синергетике

.pdf
Скачиваний:
0
Добавлен:
08.01.2026
Размер:
13.78 Mб
Скачать

Турбулентность в её обычном понимании возникает в пристеночных слоях слабовязких жидкостей или газов либо на некотором удалённом расстоянии за плохо обтекаемыми телами.

Для теоретического описания турбулентности применяются различные подходы. При статистическом подходе считается, что турбулентность порождает случайно изменяющаяся совокупность вихревых элементов различных размеров. Другим подходом является метод спектрального анализа, который дополняет статистический подход.

Рис BEE_4. Рождение турбулентности при обтекании шара потоком жидкости:

а) ламинарное течение (Re=10-2);

б) появление вихрей в кормовой части (Re=20); в) развитие вихрей (Re=102);

г) развитая турбулентность (Re>104)

Рассмотрим вначале цилиндр, ось которого перпендикулярна скорости V движущейся жидкости. На рис. 6,а) схематически показаны линии тока жидкости при малой скорости ее движения. Характер этих линий зависит не

только от скорости, но и от кинематической вязкости

 

( - вязкость

 

 

 

87

 

 

Tmax
Tmin

жидкости - ее средняя плотность) и от диаметра d цилиндра. Эти числа объединяют в безразмерный комплекс Рейнольдса, Re Vd , который более

полно, чем одна скорость описывает картину обтекания цилиндра жидкостью. Итак, при малых числах Re 20 линии тока стационарны, т. е. не меняются со временем. Но после того, как скорость превысит некий порог, появляются рециркуляционные вихри в следе за цилиндром (рис. BEE_4,б). Стационарный режим исчезает, уступая место цепочке вихрей, вращающихся попеременно то в одну, то в другую сторону. Это явление носит название вихревой дорожки Бенера - Кармана. На рис. BEE_4 изображена эволюция

вихрей для различных значений 20 Re 106 . При Re > 20 появляется пара вихрей, при Re > 102 вихри осциллируют. При еще более высокой скорости (Re > 106) появляется нерегулярная картина - турбулентный поток. В последней можно также усмотреть появление новой картины самоорганизации - порядка в хаосе. Подчеркнем универсальный характер описываемых здесь периодических явлений: след турбулентных явлений можно обнаружить в структуре облаков, создаваемых ветром над городом и за Большим красным пятном на Юпитере.

Возвратимся в нашем рассмотрении у явлению Бенара (рис. BEE_1). Управляющим параметром, определяющим смену режимов поведения жидкости при нагреве, является значение числа Рэлея R

 

 

R

Tmax Tmin

.

 

(BEE.1)

 

 

 

 

 

 

 

 

 

 

 

 

В этом

выражении

 

 

- кинематическая вязкость,

-вязкость

 

 

 

 

 

 

 

 

жидкости,

- средняя плотности жидкости, - коэффициент

теплопроводности, определяющий

соотношение

между потоком тепла W

и градиентом

температуры

T . Соотношение,

связывающее

поток тепла

W , градиент температуры T и коэффициент теплопроводности, имеет вид W T .

Вывод. В эксперименте Бенара, при котором жидкость в поле тяготения подогревается снизу при малых разностях температур

преобладает вязкость. Жидкость покоится и тепло переносится постоянной теплопроводностью. Такое состояние становится неустойчивым при критическом значении числа Рэлея Rcr . Такое критическое значение для

данной жидкости достигается увеличением разности температур Tmax Tmin . Если R Rcr в жидкости появляются стационарные конвективные валы,

88

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

направлением

поля скоростей

 

движения жидкости. В таком хаотическом

режиме присутствуют

характерные возмущения

с

длинами

волн

,

которые могут пробегать целый диапазон значений

min max .

 

 

Количественное описание явлений Бенара

 

 

 

Движущаяся жидкость

описывается

полем

скоростей

V r ,t

и

скалярным

полем

температур

T r ,t .

Основными

уравнениями,

описывающими эту систему являются уравнения

Навье-Стокса (BEE.2),

уравнение теплопроводности (BEE.3) и уравнения непрерывности (BEE.4)

 

 

dV

F P V ;

 

 

 

 

(BEE.2)

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

T ;

 

 

 

 

(BEE.3)

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

cp

 

 

 

 

 

 

 

 

 

 

div

V

0 .

 

 

 

 

(BEE.4)

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В этих уравнения cp -

удельная

теплоемкость,

r ,t -

плотность,

P

давление, -вязкость, F - внешняя сила, действующая

на элемент объема.

Если такой силой является сила тяжести, то F r ,t g , где

g -ускорение

свободного падения. Заметим, что нелинейность в гидродинамики связана с полной производной по времени

 

dVi r ,t

 

Vi

r ,t

 

 

V dx

k

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

dt

 

 

t

xk

 

dt

(BEE.5)

В векторной форме это соотношение можно представить в виде

 

 

 

 

 

dV

 

V

V

V .

 

 

 

 

 

 

 

t

 

(BEE.6)

 

 

 

 

dt

 

 

Чтобы упростить все вычисления предполагается, что система обладает трансляционной симметрией по оси Y. Учтем, что зависимостью всех постоянных нашей задачи, входящих в уравнения (BEE.2)-(BEE.4), от температуры T можно пренебречь. Единственной величиной, которая

89

Tmax , а

зависит от температуры T , является плотность жидкости T , которая

линейным образом зависит от температуры. Такое приближение называется приближением Буссинеска. Требования симметрии определяют зависимость искомых величин от своих аргументов (рис. BEE_5). Ширина плоской пластины по оси z , нижняя часть которой нагрета до температуры

верхняя часть которой имеет температуру Tmin ,

равна

h .

Ширина

конвективных валов по оси

x равна

h

. Симметрия рисунка подсказывает,

a

 

 

 

 

 

 

что искомыми являются

следующие

компоненты

скорости

Vx x, z,t ,

Vz x, z,t и температуры

T x, z,t .

Соответствующими

собственными

ax

функциями по координате x и по координате z являются sin иh

sin z .

h

Рис. BEE_5. Конвективные валы жидкости.

 

 

Оставляя

 

 

 

только

главные пространственные

Фурье

компоненты

Vx x, z,t X t ;

Vz x, z,t Y t ;

T x, z,t Z t

получаем три

обыкновенных

 

дифференциальных уравнения для

коэффициентов Фурье

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

t

,Y

t

, Z

t

 

. Вывод этих соотношений приведен в книге Г. Шустера

 

 

 

 

 

 

 

«Детерминированный хаос».

90

X X Y

Y rX Y XZ

(BEE.7)

Z XY bZ

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

коэффициентов Фурье

 

X

t

,Y

t

, Z

t

 

 

 

 

с

постоянными

значениями

 

 

 

 

 

 

 

 

 

 

,b, r

в точности соответствует

 

системе

 

дифференциальных уравнений

Э.Лоренца, которая приводит к явлению,

называемому странный аттрактор.

Это явление мы обсуждали

в

 

главе 1, параграф 1. Заметим, что параметр

r R Tmax T min будет

 

 

пропорционален

управляющему

параметру

Рэлея

R (BEE.1).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Численный

анализ

 

этой простой

 

системы

нелинейных

уравнений

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

показывает, что

величины

 

 

 

 

X

t

,Y

t

, Z

t

 

могут

проявлять свои

 

 

 

 

 

 

 

 

 

 

хаотические поведения, если параметр r

 

увеличивается и превышает свое

критическое значение rc . Компоненте

 

X t

соответствует

 

 

Фурье

 

 

 

 

 

 

 

 

 

 

компоненты X

f , зависящая от частоты

 

f , измеряемой в Герцах. На рис.

 

 

 

 

 

 

 

f

 

2 ,

BEE_6 приведен спектр мощности

 

 

компоненты

P f

X

 

соответствующей компоненте

 

скорости

 

Vx x, z,t X t . Такие

спектры

мощности измеряются при рассеянии света с помощью изучения

эффекта

Доплера.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

91

Рис. BEE_6. Спектр мощности конвективного потока в эксперименте

Бенара. При увеличении относительного числа Рэлея R* наблюдаются следующие состояния: a) верхний рисунок. периодическое движение с одной

частотой и ее

небольшими

гармониками. Относительное число Рэлея

R* 31.0; б)

квазипериодическое движение с двумя несоизмеримыми

частотами и их комбинациями

R* 35.0 ; в) непериодическое хаотическое

движение с несколькими узкими линиями в спектре R* 46.8; г) хаос,

R* 65.4

92

Вывод. Критерием хаотического поведения системы в задаче Бенара является сложный спектр Фурье исследуемой величины. При хаотическом поведении величины Vx x, z,t спектр мощности характеризуется широким

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

93

2.6. Реакция Белоусова-Жаботинского

Реакция Белоусова Жаботинского — класс химических реакций,

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

Реакция Белоусова-Жаботинского

https://www.youtube.com/watch?v=muYgGjNxB8c ( 0-30 сек; 8 min-13 min; 17 мин)

https://www.youtube.com/watch?v=GZnwjyVjH_I

В настоящее время под этим названием объединяется целый класс родственных химических систем, близких по механизму, но различающихся

используемыми катализаторами (Ce3+,

Mn2+ и комплексы

Fe2+,

Ru2+),

органическими восстановителями (малоновая кислота CH2 (COOH )2 ,

броммалоновая кислота, лимонная

кислота, яблочная кислота и

др.)

и

окислителями (броматы, иодаты и др.).

Реагентами являются Ce2

SO4

3 ;

NaBrO3 ; H2SO4 . В реакции

Белоусова-Жаботинского

органические

молекулы (например, молекулы малоновой кислоты) окисляются бромат

ионами при

катализации окислительно-

восстановительной системы

(Ce4 / Ce3 ) .

Напомним, что элемент Ce

находится в Периодической

системе элементов Д.И.Менделеева в группе лантаноидов, его порядковый номер равен 58.

Все

реагенты в реакции Белоусова-Жаботинского участвуют в N =18

элементарных реакциях.

Рассмотрим вектор концентраций

компонентов

x c1;c2

;...cd , где величины ck

имеют размерность концентрации,

 

 

dx

F x, .

 

(BZH.1)

 

 

dt

 

 

 

 

 

 

 

В этой системе из

N уравнений векторная функция F

представляет

собой нелинейную функцию,

зависящую от всех компонент концентраций

x c1;c2

;...cd . Величина

представляет собой внешний

управляющий

параметр.

От его значения зависит

характер развития реакций. В качестве

такого управляющего параметра может выступать характерное время

пребывания компонент реакции Белоусова-Жаботинского

в проточном

химическом реакторе.

 

 

 

 

При определенных условиях эти системы могут демонстрировать очень

сложные

формы

поведения

от

регулярных

периодических

 

 

 

94

 

 

до хаотических колебаний и являются важным объектом исследования универсальных закономерностей нелинейных систем. В частности, именно в реакции Белоусова — Жаботинского наблюдался первый экспериментальный странный аттрактор в химических системах и была осуществлена экспериментальная проверка его теоретически предсказанных свойств.

История открытия колебательной реакции Б. П. Белоусовым, экспериментальное исследование её и многочисленных аналогов, изучение механизма, математическое моделирование, историческое значение приведены в коллективной монографии. Борис Павлович Белоусов проводил исследования цикла Кребса, пытаясь найти его неорганический аналог. В

результате одного

из экспериментов

в 1951 году, а

именно

окисления лимонной

кислоты броматом

калия в кислотной

среде в

присутствии катализатора — ионов церия Ce+3, он обнаружил автоколебания.

Течение реакции менялось со временем, что проявлялось периодическим изменением цвета раствора от красного (избыток Ce+3) к синему (избыток

Ce+4) и обратно. Эффект ещё более заметен в присутствии индикатора ферроина. Сообщение Белоусова об открытии было встречено в советских научных кругах скептически, поскольку считалось, что автоколебания в химических системах невозможны. Статью Белоусова дважды отклоняли в редакциях советских журналов, поэтому опубликовать результаты исследований колебательной реакции он смог только в сокращённом виде спустя 8 лет в ведомственном сборнике, выходившем небольшим тиражом. Впоследствии эта статья стала одной из самых цитируемых в данной области, а реакция получила название реакции Белоусова.

Дальнейшее развитие исследований этой реакции произошло, когда

профессор С.Э.Шноль предложил

своему

аспиранту,

в

будущем

лауреату Ленинской премии А.М.

Жаботинскому, исследовать

механизм

реакции. От приглашения проводить совместные исследования Белоусов отказался, хотя выражал удовлетворение тем, что его работа продолжена. Жаботинский провёл подробные исследования реакции, включая её различные варианты, а также составил её первую математическую модель (1964). Основные результаты были изложены в книге Жаботинского

«Концентрационные колебания». В 1969 году Жаботинский с коллегами обнаружили, что если реагирующую смесь разместить тонким плоским слоем, в нём возникают волны изменения концентрации, которые видны невооружённым глазом в присутствии индикаторов. Сейчас известно

довольно

много

реакций

типа

Белоусова —

Жаботинского,

например, реакция Бриггса — Раушера.

 

 

Таким образом, реакция Белоусова-Жаботинского представляет собой окисление малоновой кислот броматом в кислой среде. Изначально Белоусов

95

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

Реакция Белоусова Жаботинского катализируется такими ионами металлов

или металлокомплексами, как ионы церия Ce3 , Mn2 , Fe2 . Для наблюдения волн или структур обычно применяются ионы Ce , за которыми удобно следить оптически, так как их оптический спектр очень хорошо известен. Концентрация ионов церия может быть найдена по селективному поглощению света этими ионами.

Рис. BZH_1. Концентрические и спиральные волны в реакции Белоусова-Жаботинского.

При определенных составах наблюдается периодическое изменение окраски раствора от красного (избыток ионов церия Ce3 ) для синего

(избыток ионов Ce4 ). В других случаях после периодических колебаний цвета: красный, синий, красный, ... с периодом около 4 минут в пробирке возникают неоднородности концентрации и образуются устойчивые красные и синие слои, которые можно наблюдать примерно в течение получаса В замкнутой системе колебания цвета постепенно затухают. Но, если непрерывно подводить реагенты и отводить конечные продукты, то колебания могут продолжаться неограниченно долго. Если реакционную смесь не перемешивать, то возникают пространственные структуры, либо подобные тем, что описаны выше, либо в виде концентрических волн или спиралей видных в широкой кювете.

96