Реакція Grey-Scott: рівняння та параметри
Дві хімічні речовини, дві швидкості дифузії й два числа, які можна крутити вручну — з цього оманливо простого рецепту модель Грея-Скотта (Пітер Грей та Стівен Скотт, 1983) породжує цілий зоопарк самоорганізованих патернів: пульсуючі плями, гілки коралу, лабіринтоподібні мазі та повзучі "черв'яки". Вона стала стандартною точкою входу для кожного, хто досліджує реакційно-дифузійні системи в дусі Тюрінга — саме тому, що її простір параметрів достатньо малий, щоб досліджувати вручну, і водночас достатньо багатий, щоб ніколи не вичерпати сюрпризи.
1. Рівняння Грея-Скотта
Модель описує реактор із безперервною подачею й гарним перемішуванням, що містить дві хімічні речовини U та V, пов'язані єдиною незворотною автокаталітичною реакцією:
Записана як зв'язані реакційно-дифузійні диференціальні рівняння в частинних похідних на 2D-області:
Кожен доданок має пряме фізичне тлумачення. У рівнянні для U: Dᵤ·∇²u розповсюджує U за рахунок дифузії; −u·v² споживає U в автокаталітичній реакції; F·(1 − u) — стала поповнення (подача), яка повертає u до 1 там, де воно було виснажене. У рівнянні для V: D_v·∇²v розповсюджує V; +u·v² — той самий реакційний доданок, тепер утворюючи V; −(F+k)·v видаляє V як через відтік подачі, так і через додаткову швидкість знищення k.
2. Чотири параметри
| Параметр | Типовий діапазон | Роль |
|---|---|---|
| F (швидкість подачі) | 0,01–0,09 | Як швидко U поповнюється ззовні; контролює загальну густину патерну |
| k (швидкість знищення) | 0,03–0,07 | Додаткова швидкість видалення V; разом із F обирає режим патерну |
| Dᵤ | ≈ 0,16 | Коефіцієнт дифузії U — зазвичай приблизно вдвічі більший за D_v |
| D_v | ≈ 0,08 | Коефіцієнт дифузії V — має бути меншим за Dᵤ, щоб патерни утворювалися |
Ключова асиметрія — D_v < Dᵤ: V ("активатор", що каталізує власне утворення) дифундує повільніше, ніж U (постачання "субстрату/інгібітора"). Ця комбінація — локальна самоактивація з довшим за радіусом виснаженням субстрату — це саме той рецепт, який Алан Тюрінг визначив у 1952 році як здатний порушувати просторову симетрію.
3. Карта параметрів F-k
Змінюючи F та k при фіксованих Dᵤ та D_v, ми виявляємо виразні іменовані області, кожна з яких відповідає іншому довгостроковому атрактору динаміки:
| Патерн | Прибл. (F, k) | Поведінка |
|---|---|---|
| Плями (перлини) | 0,035, 0,065 | Ізольовані, приблизно круглі V-плями на фоні U |
| Мітоз | 0,028, 0,062 | Плями видовжуються та розділяються надвоє під час росту — нагадує клітинний поділ |
| Черв'яки | 0,078, 0,061 | Видовжені рухомі сегменти-смуги, що блукають і час від часу розділяються |
| Смуги / лабіринт | 0,029, 0,057 | Стаціонарні лабіринтоподібні, просторово заповнюючі смугасті патерни |
| Корал / хаос | 0,082, 0,059 | Розгалужений, деревоподібний ріст із хаотичним розгалуженням кінчиків |
| Хвилі | 0,014, 0,054 | Спіральні та цільові хвилі, що поширюються, не осідаючи |
Замкненої формули для меж між цими областями немає — площину F-k досліджують емпірично та через лінійний аналіз стійкості поблизу меж. Невеликі зміни F чи k (усього 0,002–0,005) можуть перевести систему від однієї якісної поведінки до іншої, що частково пояснює, чому карту параметрів так цікаво досліджувати інтерактивно.
Симуляція реакції-дифузії
Змінюйте F та k наживо й дивіться, як плями перетворюються на черв'яків, лабіринти й корал у реальному часі на GPU-прискореній сітці.
4. Дискретизація та стійкість
На сітці з кроком h Лапласіан наближають 5-точковим (або ізотропнішим 9-точковим) кінцево-різницевим шаблоном, а ДЧП просувають явним методом Ейлера в часі:
u[i,j] ← u[i,j] + Δt·( Dᵤ·∇²u − u·v² + F·(1−u) ) v[i,j] ← v[i,j] + Δt·( D_v·∇²v + u·v² − (F+k)·v )
Явний Ейлер для дифузії стійкий лише коли крок часу задовольняє умову Куранта-Фрідріхса-Леві (CFL) для доданку дифузії:
Перевищення цієї межі породжує чисельні коливання, що зростають без обмеження — сітка помітно "вибухає" в шум за кілька десятків кроків. Через це обмеження більшість реалізацій виконують кілька підкроків на кожен відрендерений кадр замість одного великого кроку, а GPU-реалізації (через WebGL ping-pong-фреймбуфери) паралелізують обчислення шаблону для кожного пікселя одночасно, щоб підтримувати реальночасові частоти кадрів навіть за високої роздільної здатності.
5. Чому виникають патерни: нестійкість Тюрінга
Лінеаризація рівнянь навколо однорідного стану та застосування збурення з хвильовим числом k (просторова частота, не плутати з параметром швидкості знищення) дає дисперсійне співвідношення σ(k), що описує, наскільки швидко збурення цієї довжини хвилі зростає чи згасає. Дифузійно-обумовлена (тюрінгова) нестійкість вимагає, щоб швидкість зростання σ(k) була додатною для смуги ненульових хвильових чисел, навіть якщо мода k = 0 (однорідна) стійка. Необхідна умова — щоб інгібувальна речовина (тут, ефективно виснаження U) дифундувала достатньо швидко відносно самоактивувальної речовини V — що відповідає вимозі D_v < Dᵤ з розділу 2.
Хвильове число з піковою швидкістю зростання задає характерну відстань між елементами результуючого патерну — більші співвідношення Dᵤ/D_v та нижчі F зазвичай дають грубіші, ширше розставлені елементи, тоді як менші співвідношення дають тонші текстури.
6. Реальні аналоги
Візерунки шкіри тварин
Реакційно-дифузійні системи цього загального типу активатор-інгібітор — провідне пояснення періодичних візерунків на шкірі тварин: плями леопарда, смуги зебри, візерунки риби-їжака. Сама модель Грея-Скотта є спрощеною карикатурою; для кількісного порівняння з біологією використовують детальніші моделі морфогенів (наприклад, системи активатор-інгібітор Мейнгардта), але Грей-Скотт відтворює той самий якісний зоопарк патернів з набагато меншою кількістю параметрів.
Хімічні осцилятори
Хімія Грея-Скотта (хоч і гіпотетична в буквальному записі) належить до тієї самої широкої родини, що й реакція Білоусова-Жаботинського: автокаталітична стадія, пов'язана з виснажуваним субстратом. Обидві системи ілюструють, як прості локальні хімічні правила породжують глобальний просторовий порядок далеко від рівноваги.
Еталон для чисельних методів
Оскільки Грей-Скотт є жорсткою системою (швидка реакція, повільніша дифузія), вона є поширеним еталоном для перевірки неявних і напівнеявних розв'язувачів ДЧП, схем адаптивного кроку часу та стратегій GPU-паралелізації, перш ніж застосовувати їх до складніших реакційно-дифузійних систем чи систем рідин.
Часті запитання
Чому D_v має бути меншим за Dᵤ, щоб патерни утворювалися?
Якби V дифундував так само швидко чи швидше за U, будь-який локальний надлишок V розповсюдився б раніше, ніж автокаталітична реакція встигла б його підсилити, і система просто повернулася б до однорідного стану. Коли D_v менший, локальна пляма V залишається концентрованою достатньо довго, щоб спожити навколишній U і рости, тоді як виснажена область U розповсюджується назовні швидше за V, створюючи самопідсилювану пляму, оточену "ровом" низького U. Це механізм короткодіючої активації та далекодіючого пригнічення в основі формування патернів Тюрінга.
Чому в площині F-k існує так багато іменованих областей (плями, черв'яки, лабіринти)?
Нелінійний реакційний доданок u·v² означає, що локальна динаміка системи бістабільна або збудлива залежно від того, як саме F та k врівноважують подачу й видалення. Поблизу одних меж плями стають нестійкими до видовження (утворюючи черв'яків); поблизу інших плями, що ростуть, розділяються замість простого розширення (мітоз); поблизу ще інших сам однорідний стан стає повністю нестійким, породжуючи просторово заповнюючі смугасті лабіринти замість ізольованих елементів. Єдиного механізму немає — кожна іменована область відповідає окремій біфуркації локальної реакційної кінетики в поєднанні з патерноутворювальною нестійкістю.
Чи є модель Грея-Скотта реалістичною хімічною реакцією?
Не напряму — це мінімальна, ідеалізована модель, розроблена як найпростіша можлива система з багатим зоопарком патернів, а не буквальний опис конкретної хімічної суміші. Її цінність педагогічна й обчислювальна: маючи лише два види й чотири параметри, її можна вичерпно дослідити, водночас відтворюючи ті самі якісні явища (плями, смуги, хвилі), що спостерігаються в справжній реакційно-дифузійній хімії та біології.
Як Грей-Скотт пов'язаний із реакцією Білоусова-Жаботинського?
Обидві є реакційно-дифузійними системами, побудованими навколо автокаталітичної стадії, але БЖ — реальна, експериментально реалізована хімічна реакція з добре охарактеризованим механізмом (модель Орегонатора), що породжує часові осциляції та спіральні хвилі, що поширюються в збудливому середовищі. Грей-Скотт — суто математична модель без прив'язки до конкретної реальної хімії, налаштована на дослідження максимально широкого діапазону стаціонарних і квазістаціонарних патернів тюрінгового типу, а не осцилюючої поведінки.
Що станеться, якщо збільшити роздільну здатність сітки?
Збільшення роздільної здатності (більше точок сітки на одиницю довжини) за фіксованого фізичного розміру області змінює крок сітки h. Оскільки довжина хвилі патерну задається фізичними коефіцієнтами дифузії та швидкостями реакції, а не сіткою, вища роздільна здатність розкриває тонші деталі в межах того самого патерну, а не змінює сам патерн — за умови, що крок часу перемасштабовано відповідно до умови CFL (Δt ∝ h²), щоб зберегти чисельну стійкість.
Чи можна обчислити патерни Грея-Скотта в 3D?
Так — те саме ДЧП безпосередньо поширюється на три виміри з 3D-шаблоном Лапласіана (7-точковим чи вищого порядку). 3D Грей-Скотт породжує аналогічні структури: плями стають кулями, смуги — листами чи трубками, а мітозоподібне розділення відбувається в трьох вимірах. 3D-симуляція значно дорожча обчислювально (сітка N³ проти N² у 2D) і зазвичай візуалізується через вилучення ізоповерхонь чи об'ємний рендеринг, а не пряму колірну карту.
Хто винайшов модель Грея-Скотта і навіщо?
Пітер Грей та Стівен Скотт, хімічні інженери з Лідського університету, представили модель у серії статей 1983-84 років, аналізуючи гіпотетичний реактор безперервного перемішування (CSTR) із кубічною автокаталізою. Їхній оригінальний інтерес стосувався стійкості й осциляцій хімічного реактора, а не формування патернів як такого; вражаючу здатність моделі формувати просторові патерни під дією дифузії детальніше дослідили пізніші науковці, зокрема стаття Джона Пірсона 1993 року в журналі Science, яка картографувала простір параметрів F-k і назвала багато класів патернів, що використовуються досі.
Чим Грей-Скотт відрізняється від моделі Фітцг'ю-Нагумо?
Фітцг'ю-Нагумо — спрощена модель збудливого середовища (виведена з рівнянь нейрона Ходжкіна-Хакслі), налаштована на утворення хвиль на кшталт потенціалу дії та спіральних хвиль, з швидким активатором і повільною змінною відновлення. Грей-Скотт натомість налаштований на стаціонарні й квазістаціонарні патерни Тюрінга з хімічної точки зору. Обидві є двозмінними реакційно-дифузійними ДЧП і мають спільну математику нестійкості Тюрінга в основі, але наголошують на різних областях якісного простору патернів (збудливі хвилі проти статичного патерноутворення).