Генеруйте патерни хутра тварин — плями, смуги, лабіринти — за моделлю реакційної дифузії Тюрінга (варіант Грея-Скотта). Два хімічних види — активатор та інгібітор — взаємодіють для самоорганізованих просторових структур. Натисніть на полотно для посіву нових патернів.
Клікайте / тягніть по полотну для малювання активатора. Патерни виникають за 5–30 секунд.
Кроки: 0
F: 0.035
k: 0.065
Du/Dv: 0.21 / 0.10
Фізика
Модель Грея-Скотта: ∂u/∂t = Du·∇²u − u·v² + F(1−u); ∂v/∂t = Dv·∇²v + u·v² − (F+k)v. Речовина u (синя) надходить зі швидкістю F і витрачається у реакції u·v². Речовина v (червона) утворюється в реакції і знищується зі швидкістю k.
Нестійкість Тюрінга: малі збурення ростуть, коли інгібітор дифундує значно швидше за активатор. Різні F/k дають плями, смуги, лабіринти — ті самі патерни, що є на леопардах, рибах-зебра та мушлях.
Про патерни реакції-дифузії Тюрінга
У своїй знаменитій статті 1952 року «Хімічна основа морфогенезу» Алан Тюрінг запропонував, що два хімічних види, які взаємодіють — активатор і швидший за дифузією інгібітор — можуть мимовільно порушувати просторову симетрію та утворювати періодичні патерни. Цей механізм, тепер відомий як нестабільність Тюрінга, вважається причиною пігментації хутра тварин (плями леопарда, смуги зебри, малюнки риб) і формування скелету хребетних. Модель Грея-Скотта, використана тут, описує ту саму фізику: хімічна речовина u подається зі швидкістю F і споживається реакцією u·v²; речовина v продукується тією самою реакцією і знищується зі швидкістю k.
Змінюйте швидкість подачі F і швидкість знищення k, щоб досліджувати простір параметрів і спостерігати, як полотно самоорганізується в плями, смуги, лабіринти або черв'яки протягом кількох десятків секунд. Клацайте або перетягуйте по полотну, щоб засіяти нові плями активатора й спостерігати, як патерни поширюються від кожного збурення.
Часті запитання
Що таке нестабільність Тюрінга і чому вона є несподіваною?
Нестабільність Тюрінга виникає, коли просторово рівноважний стан стійкий без дифузії, але стає нестійким з її додаванням — протилежно до звичайного очікування, що дифузія все згладжує. Ключова умова: інгібітор дифундує значно швидше за активатор (Dᵥ > Dᵤ). Дрібні випадкові флуктуації потім зростають у просторовий патерн із характерною довжиною хвилі, яка визначається відношенням коефіцієнтів дифузії.
Що контролюють параметри F і k?
У моделі Грея-Скотта F визначає, як швидко свіжий активатор (u) надходить у систему, а k — як швидко продукт (v) видаляється. Різні пари (F, k) дають якісно різні патерни: мале F і помірне k — ізольовані плями; більше F — смуги або лабіринти; певні комбінації — хвилі або «черв'яки». Межа між типами патернів називається простором Тюрінга.
Які рівняння описують цю симуляцію?
Модель Грея-Скотта: ∂u/∂t = Dᵤ·∇²u − u·v² + F(1 − u) і ∂v/∂t = Dᵥ·∇²v + u·v² − (F + k)v. Доданок u·v² — автокаталітична реакція: v каталізує власне утворення за рахунок u. Тут Dᵤ = 0,2097 і Dᵥ = 0,1050, що дає відношення дифузії 2:1, достатнє для запуску нестабільності в правильній ділянці (F, k).
Чи справді патерни Тюрінга зустрічаються у тварин?
Є переконливі докази механізмів типу Тюрінга у біології. Рівномірний розподіл волосяних фолікулів у мишей потребує взаємодій, що відповідають кінетиці активатор-інгібітор. Риби-ангели перемальовують свої смуги в міру зростання — відповідно до механізму Тюрінга. Пігментація зебрафіш включає два типи клітин — ксантофори (активатор) і іридофори (інгібітор) — поведінка яких після видалення та трансплантації відповідає передбаченням Тюрінга.
Як симуляція дискретизує рівняння?
Полотно — це сітка 200×200 з вузловими (замкненими) граничними умовами. На кожному кроці лапласіан ∇²u апроксимується 5-точковою схемою скінченних різниць: ∇²u ≈ u(x−1,y) + u(x+1,y) + u(x,y−1) + u(x,y+1) − 4u(x,y). Потім явно застосовується оновлення реакції-дифузії з кроком Δt = 1. Кілька кроків на кадр анімації (повзунок швидкості) прискорює формування патерну без зміни кінцевого стану.
Яка характерна довжина хвилі патерну Тюрінга?
Довжина хвилі патерну визначається коефіцієнтами дифузії і швидкостями реакцій. З аналізу лінійної стійкості найбільш підсилене хвильове число q* задовольняє q*² ≈ √(fᵤhᵥ/DᵤDᵥ), де fᵤ і hᵥ — часткові похідні реакційних доданків в рівноважній точці. Збільшення відношення Dᵥ/Dᵤ зменшує довжину хвилі (дрібніші плями); зменшення може повністю пригнітити патерн.
Які ще моделі продукують схожі патерни?
Модель Гірера-Майнхардта (1972) — це класична формуляція активатор-інгібітор у біології розвитку. Модель Фіцхью-Нагумо продукує подібні до Тюрінга патерни і є спрощеною моделлю збудливості нейронів. Моделі Шнакенберга і Брюсселятора широко використовуються в теоретичній біології. Усі вони поділяють ключову умову: локальна позитивна зворотна реакція активатора у поєднанні з далекодіючою негативною зворотною реакцією інгібітора, який дифундує швидше.
Чому патерни з'являються через 5–30 секунд?
Нестабільність Тюрінга зростає з мікроскопічного шуму, тому потрібно багато кроків для підсилення просторових мод до видимої амплітуди. Швидкість зростання σ найбільш підсиленої моди зазвичай є малою часткою часового масштабу реакції, що означає тисячі ітерацій до насичення патерну. Повзунок швидкості множить кроки на кадр, дозволяючи обрати між повільним природним виникненням або прискореним досягненням кінцевого стану.
Що відбувається на межі простору Тюрінга?
На краях параметричної ділянки (F, k), що підтримує патерни, система зазнає біфуркацій. Поблизу нижньої межі плями нуклеюють повільно і залишаються ізольованими. Поблизу верхньої межі патерни стають нерегулярними або колапсують у однорідний стан. Певні значення параметрів дають «самореплікаційні плями», що розщеплюються при надмірному зростанні — явище, пов'язане з солітонними розв'язками моделі Грея-Скотта.