🎨 3D Патерни Реакції-Дифузії

Модель Грея-Скотта описує дві хімічні речовини U та V, що реагують і дифундують. U перетворюється на V при зустрічі: U + 2V → 3V, а V розпадається зі швидкістю k. U поповнюється зі швидкістю подачі f. Залежно від f і k, система утворює плями, смуги, спіралі або лабіринтові патерни.

Патерни Тюрінга — це просторові патерни, що виникають у системі реакції-дифузії, де активатор та інгібітор дифундують з різними швидкостями. Алан Тюрінг у 1952 році припустив, що цей механізм може пояснити формування біологічних патернів, таких як забарвлення шкіри тварин, смугастість риб та розвиток органів.

Швидкість подачі 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), ∇² — лапласіан, f — швидкість подачі, а k — швидкість знищення.

На 2D-сітці дискретний лапласіан u в комірці (i,j) наближається як: ∇²u ≈ u[i-1,j] + u[i+1,j] + u[i,j-1] + u[i,j+1] - 4·u[i,j] (4-сусідній шаблон). 9-точковий шаблон з діагональними вагами 0.05 та кардинальними вагами 0.2 дає плавніші результати з меншими артефактами сітки.

Під час ембріонального розвитку морфогени дифундують крізь тканину і реагують за кінетикою активатор-інгібітор. Принцип «локальна активація — далекодіюче пригнічення» означає: пляма активує сусідні клітини приєднатися до патерну, водночас пригнічуючи віддалені клітини, природно створюючи регулярний інтервал у плямах чи смугах.

Морфогенез — це біологічний процес, за яким організм розвиває свою форму. Системи реакції-дифузії — одна з математичних моделей того, як прості хімічні градієнти можуть створювати складні, відтворювані просторові патерни без будь-якого зверху заданого плану — виключно з локальної хімії.

Співвідношення Du/Dv є критичним. Якщо обидві речовини дифундують однаково, патерни не можуть утворитися. Активатор (V) повинен дифундувати повільніше за інгібітор (U). Зазвичай Du ≈ 2·Dv. Більше співвідношення дає тонші патерни; менше робить патерни грубішими і більш плямоподібними.

Так. У відповідному діапазоні параметрів (особливо поблизу точок біфуркації) системи реакції-дифузії можуть демонструвати просторово-часовий хаос — патерни, що постійно розпадаються і перебудовуються непередбачуваним чином. Це відрізняється від стаціонарних патернів Тюрінга і пов'язане з хімічним хаосом у реакції Білоусова-Жаботинського.

Реакція Білоусова-Жаботинського (БЖ) — це реальний хімічний осцилятор, у якому концентрації реагентів коливаються між окисненим і відновленим станами. У тонкому шарі рідини вона утворює хвилі спіралей, що розширюються — яскрава фізична демонстрація динаміки реакції-дифузії.

Пресети лише змінюють значення швидкості подачі та швидкості знищення, що подаються в ті самі рівняння Грея-Скотта, а цей параметричний простір відомий своєю чутливістю — невеликі зсуви штовхають ту саму реакцію в абсолютно різні сімейства патернів.

Клік засіває невелику кругову ділянку клітин навколо курсора свіжою хімічною концентрацією, яка потім росте і взаємодіє з наявним патерном.

Виконання оновлення реакції-дифузії п'ять разів перед кожним перемальовуванням дозволяє патерну помітно еволюціонувати та стабілізуватися значно швидше, ніж дозволив би один крок за кадр.

U повинен дифундувати швидше за V, щоб локальна активація через V могла випереджати далекодіюче пригнічення через U; якби обидві швидкості дифузії були рівними, хімічні речовини змішалися б у плоску, безпатернову рівновагу.

Він надає більшу вагу чотирьом кардинальним сусідам, ніж чотирьом діагональним, даючи плавніше, більш ізотропне наближення неперервного лапласіана.