CC BY
37
Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки. 2012. № 3(28). С.81—87
Математическое моделирование
УДК 517.957:519.21:519.62
ПЕРЕХОД «БЕСПОРЯДОК - ПОРЯДОК - БЕСПОРЯДОК»
В БИОФИЗИЧЕСКОЙ СИСТЕМЕ РЕАКЦИОННОДИФФУЗИОННОГО ТИПА
В. В. Максимов
Самарский государственный университет путей сообщения,
443066, Россия, Самара, 1-й Безымянный пер., 18.
E-mail: [email protected]
Аналитически исследована эволюция пространственных диссипативных структур, возникающих в биофизической системе реакционно-диффузионного типа во внешних шумах. Изучены поведение плотности вероятности параметра порядка, его среднего и наиболее вероятного значений, а также восприимчивости и кумулянта второго порядка в зависимости от интенсивности внешнего шума в статистически стационарном состоянии. Определена граница перехода «порядок - беспорядок». Показано, что в 'рассматриваемой системе возникает последовательность шумоиндуцированных упорядочивающего и разупорядочивающего переходов.
Ключевые слова: система «реакция - диффузия», неравновесные фазовые переходы, внешние шумы, параметры порядка.
Введение. В последние годы наблюдается все возрастающий интерес к изучению влияния внешних шумов на открытые нелинейные пространственно распределенные системы, что подтверждается появлением в ведущих научных изданиях и журналах сравнительно большого числа работ, посвященных этим исследованиям, например [1—7]. Актуальность такого рода исследований объясняется широчайшим спектром приложений таких систем к различным областям знаний, а также тем очевидным фактом, что шумы являются неотъемлемой частью окружающего нас мира и часто играют решающую роль в изменении режимов поведения динамических систем.
В работе [8] авторами развита теория, позволяющая с единой точки зрения концепции параметров порядка провести последовательное и детальное изучение процессов образования пространственных диссипативных структур (паттернов), спонтанно возникающих в открытых нелинейных пространственно распределенных системах с внешними шумами как в окрестности, так и вдали от точки перехода.
Настоящая работа является продолжением исследований, начатых авторами в [8—11], и посвящена изучению шумоиндуцированных упорядочивающего и разупорядочивающего фазовых переходов, возникающих в конкретной биофизической системе в области параметров детерминированной системы, соответствующих бифуркации Тьюринга.
Валерий Владимирович Максимов (к.т.н., доц.), доцент, каф. высшей математики.
81
Максимов В. В.
1. Математическая модель. Математическая модель изучаемой биофизической системы описывается уравнениями [12-14]
dxi
~Ж
дх2
~Ж
rx i(1 — Xi) — axix2
1 + bxi
mx2
axix2 . n „2
——----h AV X\,
1 + bxi
tSsS + D^'
(1)
где xi, x2 — функции состояния; параметры r, a, b, m, g, h, f, Di и D2 описаны подробно в [12, 14]. Исследование локальной динамики и бифуркационный анализ системы (1) представлены в [13, 14].
Введём безразмерные время т = rt, координатв1 х' = хл/r/D\ и представим параметры m/r и a/r в виде
m/r = (mo/ro )(1 + fi(x',r)), a/r = (ao/ro )(1 + f2(x',r)).
Здесь m0, r0, a0 — пространственно-временные средние соответствующих параметров, случайные однородные изотропные гауссовы поля fi(x/, т) определяют пространственно-временные флуктуации этих параметров и имеют нулевые средние и корреляционные функции вида
K [fj(r,t),fj/(r/,t/)] = (|r — r/|)exp(—ktj\t — ^|) j,
где Фj(|r — r/1) = 9j exp(—kfj|r — r/|), 9j — интенсивности флуктуаций, kfj и ktj определяют их характерные пространственный и временной масштабы. Принимая во внимание внешний шум, получим
dxi
Ifr
дх2
Ifr
ao / xix2
+ DiV/2 xi,
«оn f f , \\ xix2
— (1 + /2(x ,r))———
ro 1 + bxi
- —(1 + /1 (х',т))ж2-
ro
22
g2x2 n „/2
l + h2x2J + 2V Ж2‘
(2)
2. Результаты аналитического исследования. В работе [8] было получено уравнение Фоккера—Планка для критического параметра порядка стохастических систем реакционно-диффузионного типа. В этом разделе представлены результаты аналитического исследования системы (2), полученные на основании этого уравнения при следующих значениях параметров: ro = 1, ao = 8, g = 1,434, f = 0,093, h = 0,857, b = 11,905, mo = 0,490, r/i =
= r/2 = 1, $i = 2,4 ■ 10-5, критическое значение контрольного параметра D = 135. Известно, что в докритической области состояние любой системы является однородным статистически стационарным (беспорядок) и этому состоянию соответствует одномодальная плотность распределения вероятности. Возникновение неоднородного статистически стационарного состояния (порядок) проявляется в расщеплении максимума плотности распределения вероятности на два симметричных. На рис. 1 и 2 представлено изменение плотности стационарного распределения вероятности амплитуды критической моды £kc с увеличением интенсивности шума $2 при переходе через точку бифуркации детерминированной системы (1).
82
Переход «беспорядок — порядок — беспорядок» ...
Wat
60
40
20
0
-0,4
-0,2
0,2
0,4 £кс
Рис. 1. Стационарная плотность распределения вероятности значений амплитуды критической моды системы (2) в докритической области для двух значений интенсивности шума в2: а) в2 = 6,4 • 10-5, б) в2 = 3,2 • 10-5
Рис. 1 иллюстрирует плотность распределения вероятности амплитуды критической моды в докритической области. При малых шумах стационарная плотность вероятности близка к ^-функции, и среднее и наиболее вероятное значения модуля амплитуды критической моды (параметра порядка) совпадают и равны нулю, т.е. однородное статистически стационарное состояние системы является наиболее вероятным (см. рис. 1, а). При увеличении интенсивности шума (см. рис. 1, б) происходит деформация кривой стационарной плотности вероятности: максимум, не смещаясь, значительно уменьшается, при этом основание кривой расширяется. Так как полученное распределение вероятности не является Гауссовым, среднее значение модуля амплитуды становится отличным от наиболее вероятного. Таким образом, несмотря на отсутствие расщепления одномодальной плотности вероятности на бимодальную, среднее значение параметра порядка становится отличным от нуля и следует ожидать возникновения неоднородного статистически стационарного состояния (порядка), вероятность которого мала.
Очевидно, что в рассматриваемом случае возникновение нового состояния носит случайный характер. Объяснить это можно так. Формирование структур в докритической области может происходить на сильных неоднородностях среды. Такую неоднородность может создать сильная (крупномасштабная) флуктуация, вероятность которой невелика и зависит от параметров шума. Таким образом, при длительном наблюдении за эволюцией системы и подходящих параметрах внешнего шума можно ожидать формирования таких случайных неоднородностей, которые вызовут возникновение структур в области их нахождения.
Рис. 2 демонстрирует стационарную плотность вероятности амплитуды критической моды системы (2) в закритической области при различных интенсивностях шума. Бимодальные распределения соответствуют существованию диссипативных структур. При этом наиболее вероятное значение и математическое ожидание параметра порядка становится отличным от нуля. Рис. 2 ясно показывает, что при увеличении интенсивности шума происходит
83
Максимов В. В.
Рис. 2. Стационарная плотность распределения вероятности значений амплитуды критической моды системы (2) в закритической области для пяти значений интенсивности шума £2: а) в2 = 0, б) в2 = 6,4 • 10-5, в) в2 = 6,4 • 10-3, г) в2 = 1,6 • 10-1, д) в2 = 6,4 • 10-1
постепенное слияние максимумов, и при некотором критическом значении интенсивности шума вновь возникает одномодальная плотность. При этом |£kcmp| = 0, а (|£kc|) = 0. Система (2) переходит в состояние сильно нерегулярного поведения (беспорядок). Таким образом, полученное изменение плотности стационарного распределения вероятности амплитуды критической моды свидетельствует о существовании в системе (2) фазового перехода «беспорядок - порядок - беспорядок». Обсуждавшееся выше изменение статистически стационарных среднего и наиболее вероятного значений модуля амплитуды
84
Переход «беспорядок — порядок — беспорядок» ...
критической моды, соответствующих изменению плотностей, изображенных на рис. 2, демонстрирует рис. 3.
На рис. 4 представлена граница шумоиндуцированного фазового перехода «порядок - беспорядок» для системы (2), предсказанная с применением подхода, развитого в [8]. Следует отметить, что при приближении к детерминированной точке перехода даже очень малые флуктуации будут способствовать потере устойчивости неоднородного состояния и вызывать неупорядоченное состояние.
Определим относительные флуктуации параметра порядка (восприимчивость) в соответствии с работой [15] как х = [(£kc) — (£кс)2]/$2 и его кумулянт второго порядка к2 = (£kc)/(£kc)2. На рис. 5 изображены графики х и к2 как функции интенсивности шума $2. Наличие максимумов восприимчивости ясно показывает увеличение флуктуаций в окрестности двух критических точек. Качественно вид кривых К2 и х, полученных для системы (2), совпадает в соответствующей области с аналогичными кривыми, полученными численно в работе [15] для системы, в которой наблюдается чисто индуцированный шумом переход. Это свидетельствует о том, что в системе (2) шумоиндуцированный переход «порядок - беспорядок» имеет ту же природу.
(I£kc|)> |fk cmp|
Рис. 3. Статистически стационарные среднее (|£kc|) (сплошная линия) и наиболее вероятное |£kcmP | (штриховая линия) значения модуля амплитуды критической моды в зависимости от интенсивности шума в2 (закритическая область)
1(Г2, D
Рис. 4. Граница перехода «порядок - беспорядок» на плоскости параметров D и в2
85
Максимов В. В.
«2
1,е
1,4
1,2
1,0
0,8
Рис. 5. Восприимчивость х (сплошная линия) и кумулянт второго порядка к2 (штриховая линия) как функции интенсивности шума 02
Заключение. В работе исследована эволюция двухкомпонентной биофизической системы реакционно-диффузионного типа в поле внешних флуктуаций, имеющая важное теоретическое и практическое значение. Изучено поведение статистически стационарных характеристик первого и второго порядков амплитуды критической моды этой системы. Определена граница перехода «порядок - беспорядок». Исследование восприимчивости и кумулянта второго порядка абсолютного значения амплитуды критической моды позволяет сделать вывод о том, что шумоиндуцированные упорядочивающий и разупорядочивающий переходы, связанные с образованием диссипативных структур в этой системе, аналогичны чисто индуцированным шумом переходам и являются неравновесными фазовыми переходами второго рода. Результаты численного моделирования эволюции системы (2) [11] имеют удовлетворительное количественное соответствие результатам аналитического исследования в окрестности точки перехода «беспорядок - порядок».
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Horsthemke W., Lefever M. Noise-induced Transition. Berlin: Springer, 1984. 318 pp.
2. Garcia-Ojalvo J., Sancho J. M. Noise in Spatially Extended Systems. New York: Springer Verlag, 1999. 307 pp.
3. Lindner B., Garcia-Ojalvo J., Neimand A., Schimansky-Geiere L. Effects of noise in excitable systems// Phys. Rept., 2004. Vol. 392, no. 6. Pp. 321-424.
4. Haken H. Synergetics. Berlin: Springer, 2004. 758 pp.
5. Gardiner C. Handbook of Stochastic Methods: for Physics, Chemistry and the Natural Sciences. Berlin: Springer, 2004. 424 pp.
6. Valenti D., Schimansky-Geier L., Sailer X., Spagnolo B. Moment equations for a spatially extended system of to competing species // Eur. Phys. J. B, 2006. Vol. 50, no. 1-2. Pp. 199203.
7. Hutt A., Longtin A., Schimansky-Geier L. Additive noise-induced Turing transitions in spatial systems with application to neural fields and the Swift-Hohenberg equation // Physica D, 2008. Vol. 237, no. 6. Pp. 755-773.
8. Курушина С.Е., Громова Л. И., Максимов В. В. Стохастические уравнения и уравнение Фоккера—Планка для параметров порядка в исследовании динамики шумоиндуцированных пространственных диссипативных структур // Изв. вузов. Прикладная нелинейная динамика, 2011. Т. 19, №5. С. 45-63. [Kurushina S.E., Gromova L.I., Maksimov V. V. Stochastic equations and the Fokker-Planck equation for order parameters
86
Переход «беспорядок — порядок — беспорядок» ...
for studying of dynamics noise-induced spatial dissipative structures // Izv. Vuzov. Prikladnaya Nelineynaya Dinamika, 2011. Vol. 19, no. 5. Pp. 45-63].
9. Курушина С.Е., Желнов Ю.В., Завершинский И.П., Максимов В. В. Образование диссипативных структур в двухкомпонентных системах типа реакция-диффузия во флуктуирующей среде// Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2010. №1(20). С. 143-153. [Kurushina S.E., Zhelnov Yu. V., Zavershinskii I.P., Maximov V. V. Pattern formation in two-component reaction-diffusion systems in fluctuate environment // Vestn. Samar. Gos. Tekhn. Univ. Ser. Fiz.-Mat. Nauki, 2010. no. 1(20). Pp. 143-153].
10. Курушина С.Е. Иванов А. А., Желнов Ю. В. Завершинский И. П., Максимов В. В. Автоволновые структуры во внешней флуктуирующей среде // Изв. СНЦ РАН, 2010. Т. 12, №4. С. 41-50. [Kurushina S.E., Ivanov A. A., Zhelnov Yu. V., Zavershinskii I.P., Maximov V. V. Autowave structures in fluctuating external environment // Izv. SNC RAN, 2010. Vol. 12, no. 4. Pp. 41-50].
11. Курушина С.Е., Иванов А. А., Желнов Ю.В., Завершинский И.П., Максимов В. В. Моделирование пространственно-временных структур в системе хищник-жертва во внешней флуктуирующей среде// Матем. моделирование, 2010. Т. 22, №10. С. 317. [Kurushina S.E., Ivanov A. A., Zhelnov Yu. V., Zavershinskii I.P., Maksimov V. V. Modeling spatial-temporal structures in a predator-prey system in an external fluctuating environment// Matem. Modelirovanie, 2010. Vol. 22, no. 10. Pp. 3-17].
12. Scheffer M. Fish and nutrients interplay determines algal biomass: a minimal model // OIKOS, 1991. Vol. 62, no. 3. Pp. 271-282.
13. Malchow H. Motional instabilities in prey-predator systems// J. Theor. Biol., 2000. Vol. 204, no. 4. Pp. 639-647.
14. Malchow H. Spatiotemporal pattern formation in nonlinear non-equilibrium plankton dynamics // Procc. R. Soc. Lond. B, 1993. Vol. 251, no. 1331. Pp. 103-109.
15. Van den Broeck C., Parrondo J. M. R. Toral R., Kawai R. Nonequilibrium phase transitions induced by multiplicative noise// Phys. Rev. E, 1997. Т. 55, №4. С. 4084-4094.
Поступила в редакцию 20/III/2012; в окончательном варианте — 21 /VIII/2012.
MSC: 60H15,35R60
TRANSITION DISORDER-ORDER-DISORDER IN REACTION-DIFFUSION BIOPHYSICAL SYSTEM
V. V. Maksimov
Samara State Transport University,
18, First Bezimyanniy per., Samara, 443066, Russia.
E-mail: [email protected]
The evolution of spatial pattern formation, which arises in the biophysical system of reaction-diffusion type in external noise, is researched analytically. The behavior of the probability density of the order parameter, its mean and the most probable values, susceptibility and second-order cumulant as a function of external noise intensity are studied in the statistically steady state. The boundary of transition “order-disorder” is defined. It is shown that there is a sequence of noise-induced ordering and disordering transitions in this system.
Key words: reaction-diffusion system, nonequilibrium phase transitions, external noise, order parameters.
Original article submitted 20/III/2012; revision submitted 21/VIII/2012.
Valerii V. Maksimov (Ph. D. (Techn.)), Dept. of Higher Mathematics.
87
