Аналитическое исследование и численное моделирование контрастных диссипативных структур в поле флуктуаций динамических переменных Текст научной статьи по специальности «Физика»

Изв. вузов «ПНД», т. 17, № 6, 2009

УДК 536.75, 519.6

АНАЛИТИЧЕСКОЕ ИССЛЕДОВАНИЕ И ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ КОНТРАСТНЫХ ДИССИПАТИВНЫХ СТРУКТУР В ПОЛЕ ФЛУКТУАЦИЙ ДИНАМИЧЕСКИХ ПЕРЕМЕННЫХ

С.Е. Курушина

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

Ключевые слова: Диссипативные структуры, незатухающие моды, поле флуктуаций динамических переменных, численное моделирование.

Введение

В настоящее время ведутся интенсивные исследования явлений самоорганизации в различных неравновесных системах, заключающихся в возникновении эволюции упорядоченных пространственно-временных структур. Примером последних могут служить автоволны [1,2], которые формируются в так называемых возбудимых средах в ответ на внешнее возмущение. Существует множество примеров возбудимых сред: нервные и мышечные ткани [3], колонии микроорганизмов [4], ряд химических растворов и гелей [5,6], магнитные сверхпроводники с током [7], некоторые твердотельные системы [8].

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

= fk (Xk ,t) + DXk Ax-, k = 1, 2, 3 ... (1)

Одним из типов автоволновых процессов являются диссипативные структуры (ДС), теория которых в настоящее время активно строится в приложениях к моделям различной природы (биологическим, химическим, физическим и т. д.) [1-3,11-13].

125

Особенности формирования ДС в автономных системах типа реакция-диффузия к настоящему времени изучены достаточно подробно [11-13]. Автономная система типа (1) аналитически исследовалась в [12,14,15] и других работах. Численное моделирование формирования ДС при различных граничных условиях проведено в работах [12,16,17].

В реальной ситуации та или иная система часто находится под действием случайных внешних полей. Известно, что при этом картина поведения систем типа реакция-диффузия может качественно меняться [18,19]. Однако теоретическое изучение динамики систем вида (1) при этом чрезвычайно усложняется, и в данной области существует весьма ограниченное число работ, обычно касающихся линейной стадии образования ДС. В работе [13] в аддитивном случайном поле вблизи точки бифуркации Тьюринга методом Ланжевена проведено аналитическое исследование линеаризованной системы морфогенеза Гирера-Майнхардта [14]. Численное исследование модели, рассматриваемой в [13], проведено в [20]. При этом случайное внешнее воздействие полагалось S-коррелированным по времени и пространству. Но в [13] и [20] не учитывался нелинейный характер автокаталитических реакций и не проводилось аналитическое исследование влияния флуктуаций на условия возникновения ДС.

В работе [19] проанализировано поведение системы, состоящей из двух потребителей, ориентированных на один ресурс. В частности, показано, что при сверхкритической интенсивности шума становится возможным стационарное сосуществование обоих потребителей, в отсутствие внешних воздействий запрещаемое теоремой Гаузе [21]. Однако в работе [19] анализ проводился при условии малости отклонений концентраций видов от стационарных значений и при определенных ограничениях на интенсивность шума. Кроме того, при анализе использовался целый ряд дополнительных ограничений, серьезно сужающих область применимости полученных результатов.

В этой работе исследуется влияние аддитивного однородного изотропного поля гауссовых флуктуаций динамических переменных на образование ДС при мягком режиме их возбуждения в отсутствие указанных выше ограничений.

В качестве конкретной модели используется модель морфогенеза Гирера-Майнхардта [14]. Пониманию процессов биологического формообразования мы обязаны, прежде всего, работам Б.Н. Белинцева [13,22-24] и Г. Майнхардта [14,16,17,25].

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

да

dt

dh

dt

а2

р + k~h - Ца + DaAa + Fi(r,t), ca2 — vh + DhAh + F2(r, t).

(2)

Здесь р - скорость образования автокаталитической переменной, ц - ее постоянная распада, Da - коэффициент диффузии активатора; v - постоянная скорости распада демпфирующей переменной; член ca2 задает скорость образования ингибитора; Dh - коэффициент диффузии ингибитора. Слагаемое ka2/h учитывает, что увеличение концентрации ингибитора снижает скорость автокаталитического процесса.

126

Случайные функции Fi(r,t) определяют пространственные и временные гауссовы флуктуации скорости изменения концентраций активатора и ингибитора, соответственно, с корреляционными функциями (Fi(r,t)Fi(r', т)) = 0j/j(|r — r'|) S(t — т) и нулевыми средними значениями. Взаимная корреляция случайных функций Fi(r,t) отсутствует. Величины 0i характеризуют интенсивность флуктуаций; /i(|r — r'|) -некоторая функция, определяющая пространственную зависимость корреляций однородного и изотропного случайного поля [26]; 8-коррелированность во времени фактически означает, что время корреляции случайного поля гораздо меньше всех характерных времен задачи.

Флуктуации вызываются различными причинами. В частности, концентрация молекул может изменяться из-за тепловой и турбулентной диффузии через виртуальные или вполне реальные стенки элементарного объема (например, оболочки клеток, внутриклеточные мембраны). Необходимо также учитывать флуктуации за счет стохастического характера взаимодействия между молекулами [10]. В этом случае динамические переменные модели являются случайными функциями.

1. Аналитическое решение

В данной работе с использованием концепции параметра порядка [12] аналитически исследовано поведение системы (2) и осуществлено численное моделирование эволюции системы при двух различных типах граничных условий.

Вводя новые переменные, r' = гд/v/Da, t' = vt, a' = ca/k, h' = vch/k2, получим

„/2

= p' + h — Ц a' + Aa' + F (r, t), dh'

dF = a'2 — h' + D'Ah' + F'2(r', t'),

(3)

где p' = pc/vk, ц' = ц/v, D = Dh/Da. Далее штрихи опускаем.

Исследование устойчивости системы (3) при Fi = F2 = 0 проведено в работе [12]. Показано, что неустойчивость типа мягкой моды возникает при двух критических значениях волновых чисел kc = ± 4ц/D.

Чтобы учесть влияние случайного воздействия, перепишем (3) в виде

д_

dt

где

K (A) =

/

V

2

Р + 1

Ц (p+1»

— 1 + A

K(A)j q = g + F,

\

(4)

(p + 1)2 —1 + DA

; qi = a—ao, q2 = h—ho;

qi, q2 описывают малые отклонения концентраций относительно стационарных состояний h0 = (p + 1)2 /ц2; a0 = (p + 1)/ц. Вектор g содержит квадратичные и кубические нелинейности, полученные разложением в ряд правой части (3), и определяется так:

g = ((q2/ho — (q2q2 + 2aoqiq^ /h2+ (2aoqiq| + a§q|) /h0 — a^q^/^), q]2

2

ц

Ц

q

127

Представим вектор q в виде суперпозиции

q(r,t) = E O(j)lk\t)eikr,

k,j

(5)

где O(j) - собственные векторы линейного оператора K(А) задачи (4) (j = 1, 2).

Считая, что неустойчивость возникает при k = kc, после некоторых преобразований и проведения процедуры адиабатического исключения устойчивых мод система уравнений для амплитуд неустойчивых мод ^^(t) с точностью до слагаемых третьего порядка малости включительно принимает вид

4)

dt

Mfcc)4 = Е c(kc)444kC + kC' - kc) +

+ Е d(kc,k)i^/)lw S(kc+kc'+k'c - kc)+

f J I'-c IN'C IN'C

+

kC,k

X(kc,k)[ j O^ke ikrFi(r,t)d r+ O^ke ikrF2(r,t)dr ) +

+ e(kc,k)[ j O(2ke-ikrFi(r,t)dr + I O^ke“*krF2(r,t)dr

)(2) -ikr ;

+ f Ol}Ce ikcrF1(r,t)dr + J O2kce ikcrF2(r,t)dr.

48(k + k - kc)+

(6)

Здесь Oj - компоненты собственных векторов оператора, сопряженного K(А). Выражения для ~kj(k) и компонент векторов O(j) и O(j) [27] приведены ниже:

/

O(j)

O(j)

2(Р + 1) j >

.(1 + Dk2 + 7j) )

(-1)j

1 + Dk2 + 7,-

(-1)j

7,2 — 7i

ц (1 + Dk2 + 7i) (1 + Dk2 + X2) 2(p + 1)(7i — 72)

. a(k) a2(k)

7i,2 = -^ ±y — P(k);

a(k) = —(D + 1)k2 + 2ц/(р + 1) — ц — 1;

|3(k) = (k2 + p)(Dk2 + 1) — 2pDk2/(p + 1).

k/c,k/c

k/c,k//,k/c

1

128

Функции 6(k' + k" — k) и 6(k' + k" + kw — k) дают «правила отбора», которые важны для определения возникающих пространственных структур. Моды ^^(i) служат параметрами порядка. Их кооперация или конкуренция определяет вид возникающих структур.

В уравнениях (6) введены следующие обозначения:

е(кс)

0И

(1+Р)2

0(1)

01,kc

+

^2 f0(1) ]2___2^~0(1) 0(1) 1+0(1) (о(1) ]2 •

(1 + р)2 \O2,kc j (1 + р) O1>fcc02,kcJ + 02,kc \O1,kc) •

%(кс,к) = —

1

h(k)

2г>(1)

2 ^2O

1,kc

(1+ P)2

O(i o^

+ _____O(1) O1)

(i + p)2 2’kc 2k

— O0 O“) + 241OH Ok

e(kc ,k) = —

(1 + p) 11 2,k

1

Ык)

2р2о(1)

d(kc,k) = —

(1 + P)

^о(1 (1+P)4

^2

_______kc (O( 1 ) O(2) +__^ o( 1 ) O(2) __

(1 + p)2 v 1kc 1k (1 + p)2 2kc 2k

(Ol o22k+O0 Ok >) +2<52‘kc oi‘kc Ol

(о(ч V2O(d | p2 (O(1) v_ (O(1 n2O(d

\0 1 O2kc +(1+p)2 y^kc) (1+ p) yO2,kc) 01 k

1,kc ) 2,kc ' (1 + p)2 V 2,kc,

i(1) , „(u k)^(2)'

(1+p)

+

+ (Trp^ + е(к:,кЩ) X

x ((o1i) + (o2i)— (irpy <<) +

+x(kc,k)021k

O(1)

O1,kc

2

+ e(kc, k)0

(2) 2 ,kk

O(1)

O1,kc

2

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

Предположим, что реакция протекает в двумерном слое, размеры которого l1 и I2, и примем периодические граничные условия. Волновой вектор определяется как

( n/h \

k = 2п , n,m = 0, ±1, ±2...

I m/l2 )

Заменим вектор kc его модулем кс и углом ф, который этот вектор образует с фиксированной осью. Этот угол изменяется от 0 до п. В первой сумме уравнений (6) может остаться только слагаемое, для которого kC + k^ = kc. Это возможно, если рассматриваемая тройка векторов образует правильный треугольник. Аналогично во второй сумме уравнений (6) kc + kc' + k^ = kc. Следовательно,

129

О = — kC; k'C = kc. Модуль и угол в волнового вектора устойчивых мод k легко выражается через kc, ф и ф' (ф' - угол, который вектор kC образует с выбранной осью): k = \/2kC\J 1 — cosO — ф); tg в = — ctg(ф/ + ф)/2. В результате систему уравнений (6) для двумерного слоя можно переписать так:

d£(1)

кс’ф х (к )t(1)

dt kl(kc)'^*

c(kc)l

(1)

kc ,ф+п/3

fc(1)

^с,ф — п/3

+

где Ю1(|ф — ф'|) ®2(|ф — ф'О

+ Е d(|Ф

ф'

+ Y Н(1)

1 ^ ^с,ф

ф'

: °{i,lx(l ф : °{2,1 Х(|ф

ф'^Йф

fc(1)

0,ф'

' z2k (t) ^^2 (| ф — ф'

— ф'|) + °дИф

— ф'|) + е(|ф

+ Е ?kC)»' z1,k (tW^

+ <O1,kc z1,kc (t) + °2л

(1)

2

— ф'|);

— ф'|);

— ф'|)+

z2,kc(t),

(7)

zj,k(t) = / e ikrFj(r,t)dr

- компоненты случайного векторного поля z(t), j и k - индексные аргументы этого поля. Корреляционные функции для компонент поля z(t) имеют вид: (zj,k(t)zj,k' (т))= = 20jФj(к) S(k + k')6(t — т). Очевидно, являются функционалами компонент поля z(t).

Усредним уравнения системы (7), используя развитый в [28] подход, и, принимая во внимание ее формальные решения и раскрывая коррелятор Zj,k(t)^ с

помощью многомерного обобщения формулы Фурутцу-Новикова [26], окончательно получим

dt

40 + c(kc) (?

(1)

kc, ф+п/3

I

(1)

kc,ty — n/3

+

+ Е Оф ф'

ф'|Н ?Й)ф

(1)

kc ,ф'

+

2

(8)

+ Е (®2(|ф — ф'|)elфl(k)+ ®2(|ф — ф'|)02ф2О) (Й1^)>

где

X = h(kc) + 01 ю1(0)Ф1(0) + 02®2(О)Ф2(О), (9)

Ф1,2(0) > 0; ®1 2(0) - в общем случае, комплексные величины. Знак о над суммой означает, что в нее не входит слагаемое с ф' = ф. Из выражения (9) для собственных чисел X неустойчивых мод видно, что их действительная часть пропорциональна интенсивности флуктуаций.

В одномерном случае система (6) после процедуры усреднения примет еще более простой вид:

d( ?k1)

dt

= X

О)

о

+ 2d(kc, 0)

О)

О

fc(1)

о

2

(10)

130

Пусть для простоты 01 = 02 = 0, fi = exp(-kf |r — r'l). Здесь kf - величина, обратная радиусу корреляции. Для одномерной среды Ф = 2kf/(k2 + kf), для двумерной Ф = 2nkf/(k2 + kf)3/2, для трехмерной Ф = 4nkf /(k2 + kf )2.

Численно исследовалась зависимость действительной части X от волнового числа k при различных интенсивностях шума 0. Результаты численных расчетов для двумерной среды в случае, когда становится неустойчивой только одна мода kc (маргинальная), при 0 = 0 и для 01 < 02 < 03 < 04 приведены на рис. 1. Из рисунка видно, что при увеличении интенсивности флуктуаций концентраций активатора и ингибитора область неустойчивых мод, для которых Re(X) > 0, увеличивается.

Рис. 1. Зависимость действительной части X от волнового числа

2. Численное моделирование

Для двумерной среды было проведено численное моделирование эволюции системы (3) с граничными условиями вида [10]:

da da da da dh dh dh dh

dx x=0,y dx x=h,7 dy x,y=0 dy , dx X,y = l2 x=0,y dx x=h,7 dy x,y=0 dy

0

- I тип

и

a(0,y,t) = a(li,y,t); a(x, 0,t) = a(x,l2,t); h(0,y,t) = h(l1,y,t); h(x, 0,t) = h(x,l2,t);

da

dx

da da da ; dh dh ; dh dh

x=0,y dx x=h,f dy x,y=0 dy l ’ dx x,y = l2 x=0,y dx x=h,f dy x,y=0 dy

- II тип.

Вдали от точки бифуркации Тьюринга установившиеся ДС пичкового типа в отсутствие флуктуаций на слое большой площади при разных граничных условиях показаны на рис. 2. Процесс установления стационарных ДС пичкового типа в

Рис. 2. Установившиеся ДС на слое большой площади (400x400 единиц длины) при отсутствии флуктуаций динамических переменных. Параметры модели: р = 0.1; ц =1; kf = 1; D = 400. а - граничные условия II типа; б - граничные условия I типа. Слева концентрация активатора, справа - ингибитора

131

отсутствие флуктуаций при граничных условиях II типа и при тех же параметрах модели представлен на рис. 3.

Формирование пичков при наличии шума различной интенсивности и при разных граничных условиях представлено на рис. 4, 5, а, 6, а. Видно, что при увеличении интенсивности флуктуаций процесс формирования стационарной ДС происходит быстрее. Кроме того, из сравнения рис. 3 и 4, 5, а можно заметить, что при наличии флуктуаций формируются ДС других периодов, то есть флуктуации играют роль фактора отбора, выбирая на стадии формирования ДС наиболее устойчивые. Последний результат согласуется с выводами, сделанными по этому поводу в работе [10].

Рис. 3. Процесс образования ДС пичкового типа из однородного состояния в отсутствие флуктуаций вдали от критической точки неустойчивости Тьюринга (граничные условия II типа). Параметры модели: р = 0.1; ц =1; kf = 1; D = 400. Слой 200 х 200 единиц длины. В середине указаны соответствующие моменты времени; слева концентрация автокаталитической переменной, справа - демпфирующей

Рис. 4. Формирование ДС концентрации активатора и ингибитора при наличии флуктуаций вдали от точки бифуркации (граничные условия II типа). Параметры модели: р = 0.1; ц = 1; D = 400; kf = 1; 0 = 0.01. Слой 200 х 200 единиц длины. В середине указаны соответствующие моменты времени; слева концентрация автокаталитической переменной, справа - демпфирующей

132

Также численно исследовались процессы изменения концентраций реагентов при р = 0.8 в отсутствие и при наличии флуктуаций. Результаты моделирования вблизи неустойчивости Тьюринга показаны на рис. 5, б и 6, б. При выбранных параметрах модели D и ц (см. подписи к рисункам) критическая скорость образования автокаталитической переменной рс = 0.814.

Все численные эксперименты производились на слое такого размера, чтобы в результате образовывалась статистически стационарная ДС. При этом среднее амплитудное значение концентрации активатора оказалось равным приблизительно 51.9, а ингибитора - 27.1.

Рис. 5. Формирование ДС концентраций реагентов при наличии флуктуаций: а - вдали от точки бифуркации (параметры модели: р = 0.1; ц = 1; D = 400; kf = 1; 0 = 0.1); б - вблизи точки бифуркации (р = 0.8; ц = 1; D = 400; kf = 5; 0 = 0.2). Слой 200 х 200 единиц длины (граничные условия II типа). В середине указаны соответствующие моменты времени; слева концентрация автокаталитической переменной, справа - демпфирующей

133

Интересен тот факт, что средние амплитудные значения концентраций автокаталитической и демпфирующей переменных практически не зависят от интенсивности аддитивного шума, однако, как показывает эксперимент, сильно зависят от размеров слоя 1\ и При уменьшении 1\ и I2 (до значений 15 х 15 единиц длины

Рис. 6. Формирование ДС концентрации активатора при наличии флуктуаций (граничные условия I типа). Слой 50 х 50 единиц длины. а - Вдали от точки бифуркации. Параметры модели: р = 0.1, р =1, D = 400, kf =5, 0 = 0.01 (левый вертикальный ряд); р = 0.1, р =1, D = 400, kf = 5, 0 = 0.1 (средний вертикальный ряд). б - Вблизи точки бифуркации. Параметры модели: р = 0.8, р =1, D = 400, kf = 5, 0 = 0.2

134

Рис. 7. Изменение поля флуктуаций концентрации активатора во времени при переходе системы от однородного состояния к ДС. Параметры модели: р = 0.1; ^ = 1; D = 400, kf = 1, 0 = 0.1. Слой 200 х 200 единиц длины (граничные условия II типа)

при данных параметрах модели; см. подпись к рис. 3) амплитуда ДС уменьшается, кроме того, сами ДС становятся неустойчивыми и при большом интервале модельного времени можно наблюдать их неоднократные образования и разрушения в поле флуктуаций (образуется «пульсирующая» ДС).

Численно исследовался вопрос об изменении интенсивности флуктуаций в системе при ее переходе от однородного состояния к ДС. На рис. 7 показаны поля флуктуаций концентрации активатора в разные моменты времени. По оси аппликат отложено отношение величины флуктуации концентрации активатора о = \J(a2) — (а)2 в данной точке слоя к А - среднему значению амплитуды ДС. Усреднение производилось по ансамблю реализаций.

На рис. 8 приведены усредненные по поверхности слоя значения о/A в зависимости от времени. Эти графики иллюстрируют тот факт, что в области перехо-

Рис. 8. Усредненные по поверхности слоя относительные флуктуации концентраций активатора (а) и ингибитора (б) при различных интенсивностях аддитивного шума вдали от точки бифуркации Тьюринга (верхний ряд, р = 0.1) и вблизи этой точки (нижний ряд, р = 0.8). Параметры модели: ^ = 1; D = 400, kf = 1. Слой 200 х 200 единиц длины (граничные условия (II типа)

135

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

Заключение

Таким образом, в модели Гирера-Майнхардта вблизи маргинального состояния аддитивные флуктуации динамических переменных приводят к увеличению области неустойчивых мод (см. рис. 1), а в области контрастных структур вдали от бифуркации Тьюринга способствуют более быстрому образованию последних (см. рис. 4, 5, а, 6, а). Из рис. 5, б, 6, б видно, что вблизи точки бифуркации аддитивный шум замедляет процесс установления ДС. Аналитически показано, что действительная часть собственных чисел неустойчивых мод пропорциональна интенсивности флуктуаций и некоторой степени радиуса корреляции, зависящей от размерности пространства системы реакция-диффузия.

Приведенный анализ будет качественно справедлив и для других распределенных моделей типа реакция-диффузия в присутствии аддитивных шумов. Однако весьма широкий класс задач требует анализа поведения динамических систем в присутствии мультипликативных флуктуаций. Такие ситуации имеют место при изучении влияния «естественных» (внутренних) шумов, роль которых существенно возрастает в системах с малым числом взаимодействующих объектов или в системах, близких к своим точкам бифуркации и т.д. [1]. Кроме того, когда между неравновесной системой и окружающей ее средой нет обратной связи и среда не претерпевает систематических изменений во времени, флуктуации среды (внешние шумы) приводят к случайному поведению параметров системы [18], которые в этом случае можно разложить на постоянную и случайную компоненты. Это, в свою очередь, приводит к появлению в уравнениях, описывающих эволюцию системы, слагаемых, содержащих шумы, величины которых пропорциональны динамическим переменным. Задачи такого рода требуют отдельного рассмотрения.

Библиографический список

1. Васильев В.А., Романовский Ю.М., Яхно В.Г. Автоволновые процессы. М.: Наука, 1987.

2. Кринский В.И., Михайлов А.С. Автоволны. М.: Знание, 1984.

3. Иваницкий Г.Р., Кринский В.И., Сельков Е.Е. Математическая биофизика клетки. М.: Наука, 1978.

4. Gerisch G. // Wilhelm Roux Archiv Entwickelungsmech Organizmen, 1965, 156, 127.

5. Белоусов Б.П. // Сборник рефератов по радиационной медицине за 1958 г. М.: Медгиз, 1959, 145; то же в кн. Автоволновые процессы в системах с диффузией. ИПФ АН СССР, Горький, 1981, 176.

6. Жаботинский А.М. Концентрационные автоколебания. М.: Наука, 1974.

7. Буздин А.И., Михайлов А.С. //ЖЭТФ. 1986. Т. 90. С. 294.

136

8. Скотт Э. Волны в активных нелинейных средах в приложении к электронике. М.: Сов. радио, 1977.

9. Cohen D., White A. // SIAMJ. Appl. Math. 1991. Vol. 51. P. 472.

10. Романовский Ю.М., Степанова Н.В., Чернавский Д.С. Математическое моделирование в биофизике (Введение в теоретическую биофизику). Москва; Ижевск: ИКИ, 2004.

11. Эбелинг В. Образование структур при необратимых процессах. Введение в теорию диссипативных структур. Москва; Ижевск: ИКИ НИЦ РХД, 2004.

12. ХакенГ. Синергетика. М.: Мир, 1980.

13. Белинцев Б.Н. Динамические коллективные свойства развивающихся систем. Дисс. .. .канд. физ.-мат. наук. М.: МФТИ, 1979.

14. Meinhardt H., Gierer A. Generation and regeneration of sequences of structures during morphogenesis // J. Theor. Biol. 1980. Vol. 85. P. 429.

15. Keener I.P. Activaters and ingibitors in pattern formation // Stadies and Applied Mathematics. 1978. Vol. 59. P. 1.

16. Gierer A., Meinhardt H. Biological pattern formation involving lateral inhibition // Lectures on Mathematics in the Life Sciences. 1974. Vol. 7. P. 163.

17. Meinhardt H., Gierer A. Applications of a theory of biological pattern formation based on lateral inhibition // Journ. Cell. Sci. 1974. Vol. 15. P. 321.

18. Хорстхемке В., Лефевр Р. Индуцированные шумом переходы: теория и применение. М.: Мир, 1987.

19. Михайлов А.С., Упоров И.В. Критические явления в средах с размножением, распадом и диффузией // УФН. 1984. Т. 144. Вып. 1. С. 79.

20. СоляникГ.И., Чернавский Д.С. Математические модели морфогенеза. Препринт ФИАН, 1980, № 8.

21. ГаузеГ.Ф. Борьба за существование. Москва; Ижевск: Изд. РХД. 2000. 234 С.

22. Белинцев Б.Н. Физические основы биологического формообразования. М.: Наука. Гл. ред. физ.-мат. лит., 1991.

23. Белинцев Б.Н. Диссипативные структуры и проблема биологического формообразования//УФН. 1983. Т. 141. Вып.1. С. 55-101.

24. Белинцев Б.Н. Элементарные процессы формирования надклеточной организации при морфогенезе. Автореферат дисс... докт. физ.-мат.наук. М., 1986.

25. Meinhardt H. The Algorithmic Beauty of Sea Shells. Berlin, Heidelberg, New York: Springer-Verlag, 1999.

26. Кляцкин В.И. Стохастические уравнения глазами физика. М.: Физматлит, 2001.

27. Абрамов Е.И., Курушина С.Е. Влияние флуктуаций динамических переменных на образование диссипативных структур в модели морфогенеза Гирера - Майн-хардта // Материалы международной междисциплинарной научной конференции «III Курдюмовские чтения. Синергетика в естественных науках». Тверь, 2007. С. 48.

28. Свирижев Ю.М., Логофет Д.О. Устойчивость биологических сообществ. М.: Наука, 1978.

Поступила в редакцию 27.10.2008 После доработки 17.06.2009

137

ANALYTICAL RESEARCH AND NUMERICAL SIMULATION OF CONTRAST DISSIPATIVE STRUCTURES IN THE FIELD OF FLUCTUATIONS OF DYNAMICAL VARIABLES

S.E. Kurushina

The influence of additive homogeneous isotropic field of Gauss fluctuations of dynamical variables of Gierer-Meinhardt model to formation of dissipative structures in soft mode regime was investigated. The system of equations for description of undamped modes interaction was received. It was shown that fluctuations of dynamical variables are widening the instability region. The numerical simulation of considered model with different boundary condition was performed. It was shown that far from Turing bifurcation fluctuations further to acceleration of dissipative structures forming and change the order of location of separate peaks.

Keywords: Dissipative structures, undamped modes, field of fluctuations of dynamical variables, numerical simulation.

Курушина Светлана Евгеньевна - родилась в Самаре (1971), окончила Самарский государственный университет (1994). Работает в Самарском государственном аэрокосмическом университете. Защитила диссертацию на соискание ученой степени кандидата физико-математических наук в СГАУ (2002) по специальности 01.04.05 - Оптика. Опубликовала 30 научных работ. Область научных интересов: оптика, теория колебаний, волн, структур; динамика нелинейных стохастических систем; теория самоорганизации; прикладная математика; стохастические уравнения.

443086 Самара, Московское шоссе, 34 Самарский государственный аэрокосмический университет им. академика С.П. Королева E-mail: [email protected]

138