Дополнительные темы:Устойчивость систем с n переменными. Критерии устойчивости. Функция Ляпунова.

↑

▷ Содержание

⇐ ПредыдущаяСтр 25 из 26Следующая ⇒

10.2. Катастрофа «складка»

Пример 1. Рассмотрим систему с одним управляющим параметром. Пусть потенциал системы можно выразить в форме

                                           ,                                      (10.4)

где a – единственный управляющий параметр.

Такого типа потенциал мы имели в случае нелинейных химических систем. Форма потенциала при изменении параметра a меняется, как показано на рис.10.1.

Рис.10.1. Задача с одним управляющим параметром в теории катастроф, форма потенциала

Определим, присутствует ли в этой системе катастрофа, т.е. скачкообразное изменение формы потенциала, и найдем вид бифуркационного множества.

Для этого определим геометрическое место точек стационарного состояния на плоскости параметров q – a . Если dV/dq = 0,то q² + a = 0.

Полученное уравнение описывает геометрическое место особых бифуркационных (стационарных) точек (точек, соответствующих решению полученного уравнения). Уравнение имеет решение в отрицательной полуплоскости (см. рис.10.2). В правой положительной полуплоскости особых точек нет. Переход от отсутствия решений к двум решениям происходит, как видно, в точке a = 0, которую называют двукратно-вырожденной. Положение двукратно-вырожденной точки находят, решая уравнение d²V/dq² = 0, из которого следует, что q = 0 и a = 0. Эта особая двукратно-вырожденная бифуркационная точка представляет собой искомое бифуркационное множество в пространстве управляющего параметра a (рис.10.2б). Таким образом, в системе имеет место катастрофа, которая называется катастрофой складки – по форме кривой стационарных состояний в плоскости параметров.

Рис.10.2. Катастрофа типа «складка»: а – геометрическое место особых (стационарных) точек;
б – диаграмма решений и точка бифуркации.

10.3. Катастрофа «сборка»

Пример 2. Задача с двумя управляющими параметрами.

Пусть потенциал системы имеет вид

                              ,                          (10.4)

где a и b – управляющие параметры.

Построим стационарную поверхность в пространстве параметров  q – a – b , для чего решим стационарное уравнение

                                           .                             (10.5)

Внешний вид поверхности равновесия (стационарности), определяемый этим уравнением показан на рис. 10.3.

Рис.10.3. Катастрофа типа сборки в двухпараметрической задаче: а – поверхность равновесия в пространстве параметров; б – бифуркационное множество (сепаратриса) на плоскости
a – b.

Предположим, что в системе присутствует катастрофа, Проверим это. Определим вид геометрического места двукратновырожденных особых точек в прастранстве параметров a – b. Для этого решим уравнение

                                  .                  (10.6)

Подставив полученное значение а в уравнение стационарной (равновесной) поверхности, получим b = 2q² . Система уравнений

                                                                  a = - 3q²

                                                 b = 2q³                                                              (10.7)

определяет положение бифуркационного множества на плоскости управляющих параметров a – b. Это множество представляет собой линию, которая называется сепаратрисой, состоящую из двух симметрично расположенных ветвей, каждая из которых разделяет на плоскости параметров области с разным типом потенциала. Точка схода ветвей сепаратрисы носит название трехкратно-вырожденной особой точки. Ее положение (а она, как видно, находится в начале координат) определяется уравнением

                                                    ,                                      (10.8)

откуда q = 0 и a = b = 0.

Катастрофа канонической формы потенциала в данном примере также, видимо, из-за внешнего сходства поверхности равновесия с портновской сборкой, носит название катастрофы сборки.

Итак, если система градиентная, то с помощью теории катастроф можно исследовать структурную устойчивость системы, найти точки, линии или поверхности бифуркационного множества, а значит определить, где (при каких значениях управляющих параметров) расположена граница, разделяющая различные неравновесные фазы (состояния) системы. Для этого надо привести потенциал системы к каноническому виду и определить все n – кратновырожденные особые точки (n = 1,2…k, где k – общее количество параметров системы).

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

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

Решая стационарные уравнения, нужно определить особые точки, которые соответствуют стационарным или равновесным состояниям.

Если модель описывается градиентной системой уравнений, исследовать устойчивость стационарных состояний по форме потенциала или путем линейного анализа устойчивости (можно использовать локальный критерий Гурвица, см. приложение 1); если полное число переменных не более пяти можно воспользоваться методами теории катастроф; для не градиентных систем следует использовать глобальный критерий Ляпунова (см. приложение).

11. АКТИВНЫЕ СРЕДЫ

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

Это относится не только к химическим системам. Аналогичные явления имеют место в непрерывных средах, в которых под влиянием внешних воздействий происходят диссипативные процессы. Среды, в которых процессы диссипации зависят от процессов переноса, называются активными.

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

                                                       .                                           (11.1)

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

                                                                                     (11.2)

Для среды, в которой присутствуют потоки диффузионного типа, полученное уравнение приобретет вид

                                                    ,                          (11.3)

где D – коэффициенты переноса (коэффициенты диффузии), которые не зависят от пространственных координат.

Отклик среды на внешнее воздействие будет зависеть от характера процессов, происходящих в каждом элементе, т.е. от вида функции F_j и от интенсивности внутренних взаимодействий, мерой которых являются коэффициенты переноса.

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

11.1. Бистабильные среды

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

Ячейка горения

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

                                ,                                    (11.4)

где Q(T) – теплота горения (Дж/с); С – теплоемкость; γ – коэффициент теплообмена; Т₁ – температура окружающей среды.

Рис.11.1. Вид динамической функции F(T) для ячейки горения.

Схема работы ячейки горения

Ячейка горения имеет три стационарных состояния (рис.11.1): Т₁ – горение отсутствует, Т₃ – устойчивое горение, когда тепловыделение равно теплоотводу, и Т₂ – появляется как математическое следствие двух устойчивых стационарных состояний при Т₁ и Т₃ (вспомним, что устойчивость стационарных состояний можно определить по знаку тангенса угла наклона функции F(T) в особых точках).

Элемент среды может быть подвержен двум типам жесткого внешнего воздействия: его можно перевести из состояния Т₁ в состояние Т₃ (т.е. зажечь) и из состояния Т₃ в состояние Т₁ (т.е. погасить). Оба воздействия носят характер жесткого возбуждения.

Если учесть взаимодействие таких ячеек, то уравнение теплопроводности следует записать следующим образом:

                                               ,                               (11.5)

где - коэффициент теплопроводности.

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

Модель Шлегля

Рассмотрим химический процесс, состоящий из двух звеньев:

                                                             А + 2Х 3Х

                                                                   Х В

Кинетическая модель этого процесса имеет вид

                                   (11.6)

и носит название модели Шлегля.

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

                                            ,                                  (11.7)

где Х – координата, вдоль которой распространяется взаимодействие элементов среды.

Функция F(x) имеет 3 стационарных состояния и геометрически выглядит также, как F(Т) для ячейки горения. Реакцию можно жестким возбуждением запустить, подав в реакционную зону количество реагента х, большее некоторого порогового значения, или остановить, выведя весь реагент х из реакционной зоны.

Запишем общий для рассмотренных трех случаев вид модели бистабильной среды

                                                  .                                 (11.8)

Обозначим стационарные состояния модели следующим образом: q₁ и q₃ – устойчивые, а q₂ – неустойчивое стационарное состояние. Рассмотрим поведение бистабильной среды при наложении на нее внешнего возбуждения. Предположим, что стационарное состояние q₁ является более устойчивым, чем состояние q₃ . Тогда элементы среды, подверженные внешнему воздействию будут стремиться к состоянию q₁ и потянут за собой своих соседей. По среде пройдет волна перехода q₃ → q₁ , которая называется волной переключения.

Рассмотрим распространение волны вдоль оси х. В уравнении заменим переменную t на :

                                                   ,                                             (11.9)

где с – скорость распространения волны переключения, t – время. Тогда . Введем граничные условия:

                        при ;  при ,             (11.10)

что означает, что сначала вся система была в состоянии q₃ , а в конце перешла в состояние q₁ , которое устойчивее. Уравнение среды в новых переменных запишется так

                                           .                                  (11.11)

Введем потенциал , тогда уравнение среды примет вид

                                               .                                 (11.12)

Полученное уравнение по форме совпадает с уравнением движения частицы в потенциальном поле V с вязким трением, пропорциональным скорости (роль которой играет ). Переведем обозначения на язык механики: D (коэффициент диффузии) – масса «частицы»; с (скорость волны переключения) – динамическая вязкость; q – пространственная координата;  - время. Адекватность динамических уравнений свидетельствует, как известно, о подобии явлений. Поэтому, проследив, как будет двигаться частица, попробуем сделать выводы о поведении бистабильной среды.

Если вязкое трение отсутствует, т.е. с = 0, то величина  - энергия «частицы». «Частица», получив жесткое возбуждение, достаточное для преодоления потенциального барьера при q₂ , падая вниз, проскочит яму при q₁ поднимется вверх по склону горы на уровень, соответствующий запасу ее энергии, а затем обратно перескочит потенциальный барьер и вернется в положение q₃ (см. рис.11.2).

Рис.11.2. К вопросу определения скорости волны переключения.

Если значение вязкости с велико, то частица под действием возбуждения не сможет проскочить потенциальный барьер и перейти в состояние q₁ – она застрянет на склоне. Существует единственное значение вязкости с = с₀ , при котором частица, получив жесткое возбуждение, попадает из q₃ в состояние q₁ . При с₀ потеря энергии на трение должна быть в точности равна разности потенциала в точках q₁ и q₃ : . Это значение с₀ и определяет скорость волны переключения.

Скорость волны переключения с₀ определяется только характеристиками среды (а не величиной внешнего воздействия). С уменьшением  с₀ убывает, при  волна пойдет в другую сторону. Для произвольного вида функции F(q) нет аналитических способов расчета с₀ . В частном случае, если F(q) представляет собой полином третьей степени:

               ,            (11.13)

можно получить точное решение:

                                         .                          (11.14)

Рассмотрим поведение активной среды в зависимости от начальных условий. Допустим, что в начальный момент времени среда неоднородна, т.е. часть ее находится в состоянии q₁ , а часть – в q₃ . Проследим за эволюцией начального состояния, при этом рассмотрим несколько основных случаев.

1. Начальное состояние среды незначительно отклонено от одного из стационарных состояний (см. рис. 11.3а). Тогда со временем вся среда переходит в ближайшее стационарное состояние.

2. Часть среды находится в состоянии q₁ , часть – в состоянии q₃ . Граница будет двигаться так, чтобы во всей системе установилось более устойчивое стационарное состояние (рис.11.3б).

3. Любое другое начальное состояние можно получить как комбинацию случаев 1 и 2 (рис.11.3в). Всегда в среде происходит релаксация к одному из стационарных состояний. Поэтому две столкнувшиеся волны гасят друг друга – аннигилируют.

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

Рис.11.3. Эволюция начальных распределений в бистабильной среде:   а – случай 1; б – случай 2;
в – случай 3; г – случай 4 (образование критического зародыша)

При анализе поведения среды обнаруживается аналогия с фазовыми превращениями первого рода. Пусть фазы соответствуют состояниям q₁ и q₃. Фаза q₁ отвечает абсолютному минимуму потенциала. При фазовых превращениях 1 рода это означает, что фаза стабильна. Состояние q₃ отвечает локальному минимуму потенциала и соответствует метастабильной фазе в равновесии. Образовавшийся зародыш стабильной фазы достаточно большого размера начинает расти и дает начало двум разбегающимся волнам переключения (фронт кристаллизации, например), после чего вся среда переходит в стабильное состояние.

Критический размер зародыша определяется, как и в кинетике фазовых превращений, из баланса двух факторов: уменьшения (термодинамического) потенциала системы при переходе в более устойчивое состояние и наличия возмущения, которое повышает потенциал.

 Если оба состояния равнозначны ( ) , то возможно стационарное сосуществование двух фаз с плоским переходным слоем (аналог границы зерна в поликристалле), несмотря на то. Что это повышает энергию системы. Когда размер системы много больше толщины переходного слоя, система может разбиться на произвольное количество областей – доменов, соответствующих разным фазам.

11.2. Возбудимые среды

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

Рассмотрим вновь ячейку горения. Предположим, что в ней присутствует ингибитор, т.е. вещество, ухудшающее условия горения. Тогда тепловой эффект горения будет зависеть не только от температуры Т, но и от концентрации ингибитора - c_in (чем больше c_in , тем меньше Q). В этом случае имеем следующую модель процесса горения:

                                 ,                     (11.15)

а модель среды

                                     .                                (11.16)

Если концентрация ингибитора не меняется в процессе горения, т.е. c_in играет роль внешнего параметра, среда остается бистабильной. При малых значениях c_in и Т > T₁ (см. рис.11.4) в среде распространяется волна зажигания; при тех же условиях, если c_in велико, в среде распространяется волна гашения.

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

                                        ,                                   (11.17)

где – характеристическое время релаксации, которое велико по сравнению с временем перехода от холодного состояния к горячему;  - равновесная концентрация ингибитора для данной температуры, монотонно растет с ростом Т.

Уравнение для ингибитора вместе с моделью среды описывает одиночную волну (уединенный бегущий импульс см. рис.11.4).

Рис.11.4. Распространение волны в возбудимой среде: а – изменение температуры в волне;
б – изменение концентрации ингибитора в волне

Действительно, при загорании происходит резкий рост температуры – это фронт импульса, волна загорания. Затем происходит постепенное накопление ингибитора, и температура снижается, пока не опустится до предельного значения, при котором еще возможно горение: Т = Т₂ . Тогда наступает резкий спад температуры – волна гашения – и состояние среды возвращается к исходному. Концентрация ингибитора меняется плавно, без скачков, что является следствием большого времени релаксации – .

БИБЛИОГРАФИЧЕСКИЙ СПИСОК

1. Пригожин И., Кондепуди Д. Современная термодинамика. От тепловых двигателей до диссипативных структур. – М.: Мир, 2009. – 461 с.

2. Хакен Г. Синергетика. – М.: Мир, 1980. – 404 с.

3. Эбелинг В. Образование структур при необратимых процессах. – М.: Мир, 1979. – 279 с.

4. Жуховицкий А.А., Шварцман Л.А. Физическая химия. – М.: Металлургия, 2000. – 688 с.

5. Бокштейн Б.С. Диффузия в металлах. – М.: Металлургия, 1978. – 248 с.

6. Гилмор Р. Прикладная теория катастроф. Кн.1. – М.: Мир, 1984. – 350 с.

7.  Лоскутов А.Ю., Михайлов А.С. Введение в синергетику. – М.: Наука, 1990. – 272 с.

8. Петелин А.Л. Основы синергетики для металлургов: Курс лекций. – М.: МИСиС, 1993. – 118 с.

ПРИЛОЖЕНИЕ

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

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

                                                      .                        (П.1)

Так же, как и раньше, определим стационарные значения q_j^(s) из уравнений

                                                        F_j(q₁,q₂,...q_n) = 0.                                     (П.2)

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

                                                  ,                                           (П.3)

где

Решение этой системы имеет вид:

         .          (П.4)

Для нахождения p₁ – p_n запишем характеристическое уравнение этой системы

                                                   ,                       (П.5)

или

                                                                             (П. 6)

Полученное алгебраическое уравнение n – ой степени в некоторых частных случаях (приведенное кубическое, биквадратное и т.д.) можно решить, найдя значения всех n корней. Однако общие выводы относительно устойчивости системы можно сделать не решая данного уравнения. Они заключаются в следующем:

- если действительные части всех (!) корней p_i меньше нуля – Re p_i < 0, то отклонения со временем затухают, т.е. соответствующая особая точка (стационарное состояние) устойчива;

- если хотя бы один из корней имеет Re p_i > 0 , то все возмущения будут со временем неограниченно возрастать, особая точка неустойчива;

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

Познавательные статьи:



Психологические особенности спортивного соревнования

Приготовление дезинфицирующих растворов различной концентрации

Занятость населения и рынок труда

Социальный статус семьи и её типология