Як це працює
Модель Грея-Скотта описує дві хімічні речовини U та V на дискретній сітці. U дифундує швидше та поповнюється із зовнішнього резервуару; V дифундує повільніше та виробляється реакцією U + 2V → 3V. Конкуренція між локальною активацією (V каталізує власне виробництво) та глобальним інгібуванням (U вичерпується швидше, ніж дифундує назад) створює самоорганізовані просторові патерни.
Кожен крок: обчислюємо лапласіан кожної клітинки за 9-точковим зваженим шаблоном, потім одночасно оновлюємо U та V методом явного Ейлера. За один кадр анімації виконується кілька кроків для прискорення формування патернів.
∂v/∂t = Dv∙∇²v + u∙v² - (f+k)∙v
∇²u[i,j] = 0.2∙(N+S+E+W) + 0.05∙діаг - u[i,j]
Колір кожного пікселя кодує концентрацію хімічної речовини V: низька V — темно-синя, висока — яскраво-бурштинова або біла залежно від вибраної кольорової карти. Клік по канвасу засіває невелику пляму з високою концентрацією V, запускаючи ріст нового патерну з цієї точки.
Часті запитання
Що таке модель реакції-дифузії Грея-Скотта?
Модель Грея-Скотта описує дві хімічні речовини U та V, що реагують і дифундують. U перетворюється на V при зустрічі: U + 2V → 3V, і V розпадається зі швидкістю k. U поповнюється зі швидкістю подачі f. Залежно від f і k система утворює плями, смуги, спіралі або лабіринти.
Що таке патерни Тюрінга?
Патерни Тюрінга — просторові патерни, що виникають у системі реакції-дифузії, де активатор і інгібітор дифундують з різними швидкостями. Алан Тюрінг запропонував у 1952 р., що цей механізм пояснює утворення біологічних патернів: плями на шкірі тварин, смуги риб, розвиток органів.
Що контролюють швидкість подачі f та швидкість знищення k?
Швидкість подачі f визначає, як швидко U поповнюється із зовнішнього резервуару. Швидкість знищення k визначає, як швидко V видаляється. Різні пари (f,k) дають різні патерни: малі f/k — плями, середні — смуги, великі — лабіринти.
Які ДРЧ описують модель Грея-Скотта?
∂u/∂t = Du∙∇²u - u∙v² + f∙(1-u); ∂v/∂t = Dv∙∇²v + u∙v² - (f+k)∙v. Тут Du і Dv — коефіцієнти дифузії (Du > Dv), ∇² — оператор Лапласа.
Як обчислюється лапласіан на дискретній сітці?
На 2D сітці 9-точковий шаблон з ваговими коефіцієнтами 0.05 для діагоналей та 0.2 для кардинальних напрямів: ∇²u ≈ 0.2∙(N+S+E+W) + 0.05∙(NE+NW+SE+SW) - u[i,j]. Це дає плавніші результати.
Чому реакція-дифузія утворює патерни шерсті тварин?
Під час ембріонального розвитку морфогени дифундують крізь тканину та реагують за принципом активатор-інгібітор. Принцип локальної активації та далекодіючого інгібування природно створює правильні відстані між плямами чи смугами, що відповідають реальним тваринним забарвленням.
Що таке морфогенез?
Морфогенез — біологічний процес, завдяки якому організм набуває форми. Системи реакції-дифузії — одна математична модель, яка пояснює, як прості хімічні градієнти можуть створювати складні просторові патерни без жодного плану зверху.
Як коефіцієнти дифузії впливають на патерн?
Відношення Du/Dv є критичним. Якщо обидві речовини дифундують однаково, патерни не утворюються. Активатор (V) повинен дифундувати повільніше за інгібітор (U). Типово Du ≈ 2∙Dv. Більше відношення дає дрібніші патерни.
Чи можуть системи реакції-дифузії демонструвати хаос?
Так. Поблизу точок біфуркації системи можуть проявляти просторово-часовий хаос — патерни, що безперервно руйнуються та відновлюються непередбачуваним чином. Це пов'язано з хімічним хаосом у реакції Бєлоусова-Жаботинського.
Що таке реакція Бєлоусова-Жаботинського?
Реакція Бєлоусова-Жаботинського (БЖ) — реальний хімічний осцилятор, де концентрації коливаються між окисленим та відновленим станами. У тонкому шарі рідини вона утворює спіральні хвилі — яскраву демонстрацію динаміки реакції-дифузії.