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

книги / Нестационарные структуры и диффузионный хаос

..pdf
Скачиваний:
6
Добавлен:
12.11.2023
Размер:
14.86 Mб
Скачать

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

Типичный пример возникновения локализованной неста­

ционарной структуры в поле сходящейся цилиндрической СВЧ-

волны показан

на рис.

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

Соседние файлы в папке книги