книги / Нестационарные структуры и диффузионный хаос
..pdfтаком подходе учитываются в изобарическом приближении. Для функции Е выполнены условия излучения.
Типичный пример возникновения локализованной неста
ционарной структуры в поле сходящейся цилиндрической СВЧ-
волны показан |
на рис. |
2.11. |
Видно, |
что |
рост |
амплитуды |
распределения |
сопровождается |
сокращением |
его |
полуширины |
||
(рис. 2.11,а). На рис. 2.11,6 |
представлен |
рост |
концентра |
|||
ции ng в центре волны. |
Видно, |
что на |
определенном времен |
ном интервале она выглядит так, как если бы процесс шел в режиме с обострением. И действительно, представления, развитые при анализе тепловых структур, эффективно исполь зуются и в теории СВЧ-пробоя.
Режимы с обострением характерны и для других модель ных уравнений. Например, недавно было показано, что упро
щенной |
моделью многих |
явлений в |
гидродинамике, |
астро |
|
физике, |
физике |
плазмы и |
нелинейной |
оптике является |
уравне |
ние |
|
|
|
|
|
|
Р, + |
(pvK = ° ’ |
vt + vvx = Clm(PVm)x • |
|
Это уравнение движения газа с отрицательной сжима
емостью описывает не бегущие волны, типичные для обычных
газов, |
а стоячие |
возмущения, нарастающие со временем в |
|
режиме с обострением. Одна из первых моделей такого |
типа |
||
была |
предложена |
русским ученым С.А. Чаплыгиным в |
1896- |
1902гг. Однако общность этой модели была осознана только в последние годы [383]. Интересно, что с помощью преобразо вания годографа эта модель сводится к линейному уравнению.
Локализованы могут быть не только стационарные дисси пативные структуры или структуры, развивающиеся в режиме с обострением, но и более сложные колебательные или хаоти ческие процессы. Вероятно, в ближайшие годы анализ таких локализованных процессов будет активно развиваться.
61
§2.2. Диссипативные структуры в средах
стриггерными свойствами
Модель тепловых структур и другие модели, связанные с одним параболическим уравнением, не обладают двумя важными особенностями, характерными для многих нелинейных сред. В
силу принципа максимума в таких системах не |
могут |
возни |
||||
кать |
новые |
экстремумы, |
и, |
следовательно, |
не |
могут |
появляться новые структуры. Кроме того, в рассмотренных средах устойчивы только простые структуры, для создания сложной упорядоченности нужно задать специальным образом начальные данные.
Описание сложных устойчивых структур, процессов их
возникновения требует перехода к системе уравнений. В ряде
случаев и здесь удается выяснить, каким образом из простых структур могут быть построены сложные, как возникает
упорядоченность в таких средах, как меняется ход процессов
при переходе от одномерных к многомерным моделям. Эти воп росы были исследованы для систем с триггерными свойствами [128]. Обратим внимание на некоторые из полученных
результатов. |
|
|
|
|
|
|
|
|
Пусть |
при />1 |
= |
0, |
£>2 = |
0 |
(когда |
нет диффузии), |
систе |
ма может |
находиться |
в |
двух |
устойчивых состояниях |
и |
|||
и |
в 0ДН0М |
неустойчивом |
(0,0). |
Такие системы |
полу |
чили название бистабильных или триггерных. Триггерные сис темы с диффузией возникают в химической кинетике, в нели нейной оптике, при моделировании морфогенеза [33].
Для простоты возьмем следующую модель:
и(х,0) = uQ(x), v(x,0) = vQ(x), |
0 |
£ х < /, |
(2.24) |
||
u j О,/) |
= ux(l,t) = vjfi.t) = vx(l,t) = |
0; |
|
||
и и v |
здесь |
можно рассматривать как отклонения от равнове |
|||
сия концентраций каких-либо веществ. |
Понятно, что и и v |
||||
могут |
быть |
как положительны, так |
и |
отрицательны. |
Будем |
62
считать, что собственные значения матрицы А - комплексно
сопряженные числа с положительными действительными частя ми. (Это является достаточным условием неустойчивости ну левого решения задачи (2.24).) Матрица В определена отри цательно, что является достаточным условием ограниченности решения при 0 < t < оо. Задача (2.24) инвариантна относи тельно преобразования {u,v} —> {-и, -о}, поскольку правые
части уравнений содержат только нечетные степени функций,
поэтому |
ы2 |
= |
- |
ы, = - |
и* v2 |
= — vx = |
— v*, (Qt(u.,Vj) = 0; |
|
t = 1,2; |
j |
= |
1,2). |
Будем |
рассматривать такие |
значения |
||
и |
/, |
при |
которых неустойчивость |
Тьюринга |
не возни |
|||
кает. |
|
|
|
|
|
|
|
|
Рис. 2.12. Пример |
элементарной |
структуры в |
три герной среде. Параметры |
||||
расчета: Л . = 0,18; |
= 0,576; |
Q. = и |
- |
v |
+(4v |
- и)(и 2 + о2); Q = и + |
|
|
+ v |
- |
(4о |
+ |
ы)(ы2 + |
у2) |
Вкачестве начальных данных вначале зададим
«ступеньку» (uQ |
= и* vQ = |
v* при |
0 £ |
х s а, и = -и* |
v = |
= -v* при а £ |
х ^ /). В расчетах наблюдается быстрый выход |
||||
на стационарное |
решение, |
типичный |
вид |
которого показан |
на |
рис. 2.12. Возникает переходная область, слева от которой
находится один устойчивый |
фон (u*v*), справа |
другой (- и * |
-V *). Такое решение мы |
назовем элементарной |
структурой |
(смысл названия мы поясним далее). Можно в большом интер вале менять длину, области и параметр а, возникающая эле ментарная структура при этом практически не меняется. Она эффективно локализована.
63
|
Обратим |
внимание |
на |
то, |
что |
в |
модели |
тепловых |
|||||
структур |
была |
возможна |
строгая |
локализация: T(x,t) |
= |
0 вне |
|||||||
(Уд |
при |
0 |
< |
( < ({. Здесь есть эффективная локализация: |
|||||||||
можно задать |
малое |
число е |
и указать |
область |
G, |
вне |
кото |
||||||
рой |
шах |
[(ы - |
u 'f |
+ (v |
- |
i>*)2] |
< е |
или |
шах |
[(и |
+ |
и*)2 + |
|
|
*€G |
|
|
|
|
|
|
|
|
x€ G |
|
|
|
+ (v + i/*)2] |
< |
e. |
|
|
|
|
|
|
|
|
|
0,5
О
-0 ,5
-1,0
Рис. 2.13. Пример сложной стационарной структуры в три герной среде. Па раметры средыте же, что на рис. 2.12
|
В отличие от модели тепловых структур, здесь явление |
|||||||||||
локализации связано |
не |
с |
источниками, |
а |
со |
стоками, |
так |
|||||
как |
в окрестности |
каждого |
из |
устойчивых |
фонов |
Q^(u,v) и |
||||||
Q2(« ,P) ведут себя |
как |
стоки. |
Кроме того, |
если |
в |
триггер |
||||||
ной |
среде |
возникли |
структуры, |
то далее |
решение |
от |
времени |
|||||
не зависит. |
|
|
|
|
|
|
|
|
|
|
|
|
|
Поскольку элементарные структуры локализованы, из иих |
|||||||||||
можно построить более сложные конфигурации, задав в |
ка |
|||||||||||
честве начальных данных несколько «ступенек». |
Расчеты |
по |
||||||||||
казывают, |
что их действительно |
можно |
построить, |
если |
рас- |
|||||||
|
|
|
|
|
|
|
|
|
|
Л / |
/V |
|
стояние между «ступеньками» не слишком мало (L > 1,1 опре деляется свойствами нелинейной среды). Пример решения, возникающего в этом случае, показан на рис. 2.13. Все ос тальные структуры в этой среде состоят из множества эле ментарных, находящихся на разном расстоянии друг от друга. Метод сшивания, оказавшийся очень полезным при исследова нии модели тепловых структур, и в этом случае дает хорошие
64
зультады, т. е. приближенные решения, соответствующие эле
ментарным структурам [161].
± =1,55
10 |
2 0 |
30 се* 4 |
Рис. 2.14. Последовательное появление структур в соседних пространствен ных областях: = 9,18; Z>2 = 0,576; А = I I ; В - 1“]_ 4 1
65
Одним из |
самых интересных свойств этой нелинейной |
среды является |
возможность рождения диссипативных структур |
и их самодостройки. Для того чтобы создать сложную упоря
доченность, совсем не обязательно задавать ее извне,
«навязать» нелинейной среде, задав нужное число ступенек. Посмотрим, что произойдет с малым возмущением, поставлен ным на неустойчивом фоне. Пример такого расчета показан на рис. 2.14. Видно, что на соседних участках последовательно
возникают структуры. При этом в расчетах наблюдается чет
кая граница области, где происходит воспроизведение дисси
пативных структур, что позволяет говорить о волновом ха
рактере процесса и отождествлять эту границу с фронтом распространяющейся волны.
Такое поведение характерно и для многих других двух компонентных систем, когда в начале имеются малые возмуще
ния вблизи потерявшей устойчивость термодинамической вет
ви. Например, для одной из наиболее известных двухкомпо
нентных |
систем - |
модели |
брюсселятора, |
описывающей некото |
||||
рый |
класс |
химических реакций |
в открытых системах, |
Q. = А - |
||||
- (В |
+ |
1)ы |
+ и v, |
Q2 = |
Ви - |
и v [151] |
или модели |
Гирера - |
Мейнхардта, возникающей при математическом описании морфо
генетических |
процессов |
Q1 |
= р |
+ |
ku2v~* - ци, |
Q2 = |
си2 - |
||||
- vv. |
(Здесь |
А В, |
р, |
к, |
ц, |
с, |
v |
- постоянные |
величины, |
||
являющиеся параметрами |
модели.) |
|
|
|
|
|
|||||
Триггерные среды рассмотренного типа обладают |
инте |
||||||||||
ресной |
особенностью - |
диссипативные |
структуры |
в |
одномерном |
и многомерном случае в таких системах качественно различа
ются. В многомерном случае возникает большой набор локали
зованных нестационарных структур, время существования ко
торых |
конечно. |
|
|
|
|
|
|
|
Будем |
рассматривать |
двумерный |
аналог |
уравнения |
||
(2.24), в котором члены D^u |
и &<ухх |
заменены |
соответст |
||||
венно |
на div(Z>1gradu) |
и div(Z>2gradt/), |
в |
квадратной области |
|||
со стороной, |
равной /, |
при |
условии |
отсутствия |
потоков на |
границе. Поставим в качестве начальных данных в правом нижнем углу области квадрат со стороной /, в пределах ко-
66
торого |
и = |
и*, |
v = |
v*, |
в |
остальной |
области |
пусть и |
= |
- |
и*, |
||||
v = |
- |
* |
Величину |
/ |
выберем |
|
таким |
образом, |
чтобы |
в |
одно |
||||
v . |
|
||||||||||||||
мерной |
задаче |
при |
|
« 0(*) |
= |
и*, |
v0(x) |
= |
v*(0 - |
х |
- |
0.* |
|||
uQ(x) |
= |
- |
и*, |
VQ(x) |
= |
- |
v* |
(/ |
^ х |
£ /) |
возникали |
стацио |
нарные структуры. За время т возникает переходная область, структура которой остается постоянной для различных точек границы, несмотря на то, что кривизна контура Г0, ограни
чивающего |
область G0, где и(х,у) > 0, меняется от точки к |
точке (см. |
рис. 2.15). |
Область GQ медленно снмметризуется, стремясь к кругу, радиус которого уменьшается со временем. Увеличение пара
метра / практически не меняет процесса, выход на стацио
нарное решение не наблюдается ни при каких |
значениях /. |
||
Независимость |
переходной области от начальных данных, эф |
||
фективная локализация, большие |
времена существования, |
||
крторые могут |
на несколько порядков |
превышать |
т, позволяют |
рассматривать этот процесс как эволюцию нестационарной диссипативной структуры.
Эти качественные представления позволяют построить приближенную модель, описывающую эволюцию контура Г(*).
Можно убедиться, что если переписать двумерный аналог сис
темы |
(2.24) в |
цилиндрической^ системе координат, |
то |
в нем |
||||
появятся члены |
~ 1 |
и ~ — |
Расчеты показывают, |
что |
||||
вдоль |
контура |
Г(^) |
профиль практически |
такой |
же, |
как |
в |
|
случае |
элементарной |
структуры |
(которая |
описывается |
уравне- |
67
ниями (2.24)). Поэтому именно слагаемые ~ |
Л |
л |
д у и ~ |
Vr |
приводят к деформации контура. Следовательно, для точек контура Г можно записать уравнение
dr (х,у, t) |
сп |
дг |
г |
где п - внешняя нормаль к контуру, г - радиус кривизны в данной точке, с определяется свойствами среды. Такая упро
щенная модель хорошо согласуется с |
результатами |
ряда |
|||||||
численных расчетов. |
|
|
|
|
|
|
|||
|
В двумерной задаче малое возмущение на неустойчивом |
||||||||
фоне |
вначале, |
как |
и в |
одномерном |
случае, |
приводит |
к |
после |
|
довательному |
возникновению структур (см. |
на |
рис. |
2.16 |
|||||
*2’ |
*3^' °Д »ако |
адесь |
структуры |
после |
их |
возникновения |
медленно перестраиваются. Радиусы структур, имеющих форму
круга, медленно уменьшаются (см. на рис. 2.11 *4, |
<5> t&). |
|||
На первый взгляд кажется, что в двумерном случае в |
||||
триггерных средах |
вообще |
нет |
стационарных, не |
зависящих от |
времени структур. |
Однако |
это |
не так. Построить |
такие реше |
ния позволяют соображения симметрии. В самом деле, уравне
ние |
не меняется при замене х на - |
х, |
у |
на |
- |
у, u,v |
на - и, |
- v. |
Будем обозначать область, |
где |
и |
< |
0, |
черным |
цветом, |
где и > 0, белым. Эволюция одной структуры, имеющей форму
круга, |
будет выглядеть как сокращение черного |
круга на |
|||
белом фоне или белого круга на черном. |
|
|
|
||
Теперь поставим «черное» и «белое» в одинаковые усло |
|||||
вия. |
Например так, как показано |
на рис. |
2.17, а |
справа. |
|
Если бы начали уменьшаться нерные квадраты, |
то |
это |
означа |
||
ло бы |
нарушение симметрии между |
черным и |
белым, |
которое |
присутствует в уравнении (2.24). Расчеты показывают, что в
этом случае действительно возникает стационарная структура
«крест» (см. |
рис. |
2.17, а слева). |
Резонансное возбуждение в |
|||||
такой |
среде |
определяется |
симметрией, |
конкретный |
вид |
|||
начальных данных, в отличие от |
модели тепловых структур, |
|||||||
здесь |
оказывается |
не |
важен. |
|
|
|
|
|
|
Интересно, что метод сшивания позволяет с высокой |
|||||||
точностью построить |
приближенное |
решение, |
соответствующее |
68
Рис. 2.16. Возникновение |
структур и |
их последующий |
распад |
в двумерной |
задаче. |
Параметры |
расчета те же, что |
на рис. |
2.14 |
69
Рис. 2.17. Примеры двумерных стационарных структур в триггерной среде
70