Про симуляцію фронту реакції
Рівняння Фішера-КПП — названо на честь Рональда Фішера (1937) та незалежно Колмогорова, Петровського й Піскунова — є одним із найфундаментальніших диференціальних рівнянь у математичній біології та хімії. Воно поєднує два фізичні процеси: дифузію, яка поширює речовину в просторі, та логістичний ріст, що моделює автокаталітичне підсилення, яке насичується при вичерпанні реагенту. Результатом є мандрівна хвиля, що зберігає свою форму нескінченно довго, поширюючись з постійною швидкістю, визначеною виключно двома параметрами D та r.
Ця симуляція розв'язує 1D рівняння Фішера-КПП числовим методом кінцевих різниць з явним методом Ейлера. Просторова область має N=600 комірок з граничними умовами Неймана (нульовий потік). Крок за часом адаптивно вибирається для задоволення критерію стійкості dt ≤ dx²/(2D), що гарантує обмеженість числового розв'язку. На кожному кадрі анімації виконуються три підкроки для покращення часової роздільної здатності при збереженні 60 кадрів за секунду.
На головному полотні відображається поточний профіль концентрації u(x) у вигляді заповненої площі з накладеною лінією. Вставка-каймограф у верхньому правому куті показує просторово-часову діаграму: кожен горизонтальний рядок є знімком u(x), захопленим у послідовні моменти часу, що прокручується вниз. Діагональна смуга фронту хвилі на каймографі безпосередньо кодує швидкість хвилі — крутіша діагональ означає швидше поширення. Порівняйте виміряну швидкість хвилі з теоретичним передбаченням v = 2√(D·r) для верифікації формули Фішера-КПП.
Часті запитання
Що таке рівняння Фішера-КПП?
Рівняння Фішера-КПП (Колмогорова–Петровського–Піскунова) ∂u/∂t = D∂²u/∂x² + ru(1-u) описує поширення вигідного гена (Фішер, 1937) або автокаталітичної хімічної реакції. Член D∂²u/∂x² моделює дифузію, а ru(1-u) — логістичний ріст, який насичується при u=1.
Що таке мандрівний фронт хвилі?
Мандрівний фронт хвилі — це розв'язок u(x,t) = f(x - vt), який зберігає форму, рухаючись з постійною швидкістю v. У рівнянні Фішера-КПП будь-який початковий стан з компактним носієм еволюціонує до хвилі, що рухається з мінімальною швидкістю v = 2√(Dr).
Що означає формула швидкості хвилі v=2√(Dr)?
Швидкість хвилі v = 2√(Dr) показує, що більша дифузія (D) і швидший ріст (r) обидва збільшують швидкість поширення. Ця формула випливає з лінеаризації рівняння біля u=0: фронт 'тягнеться' своїм переднім краєм, де u мале і ріст приблизно експоненційний.
Що таке автокаталіз?
Автокаталіз — хімічна реакція, де продукт каталізує власне утворення. Якщо існує невелика кількість продукту, він прискорює подальше виробництво, спричиняючи експоненційне зростання до вичерпання реагентів — саме така поведінка ru(1-u) моделюється тут.
Чим бістабільна хвиля відрізняється від хвилі Фішера-КПП?
Бістабільна хвиля використовує кубічний член реакції r·u·(u-a)·(1-u) з порогом a. На відміну від Фішера-КПП, яка поширюється з будь-якого позитивного зародку, бістабільні хвилі поширюються лише якщо початкове збурення перевищує поріг a.
В чому аналогія з нервовим імпульсом?
Нервовий імпульс Ходжкіна-Хакслі є мандрівною хвилею в системі реакції-дифузії, де напруга мембрани u дифундує вздовж аксона, а іонні канали забезпечують бістабільну нелінійність. Модель Фішера-КПП — спрощена версія, що захоплює основну фізику поширення хвилі.
Що таке принцип найбільшого спуску для вибору швидкості хвилі?
Якщо початковий стан спадає швидше ніж e^(-λx) для λ = √(r/D), хвиля Фішера-КПП встановлює мінімальну швидкість v_min = 2√(Dr). Крутіші початкові умови породжують швидші хвилі; швидкість 2√(Dr) є найповільнішою можливою стабільною хвилею.
Що таке просторово-часова діаграма (каймограф)?
Вставка показує каймограф: 2D зображення, де горизонтальна вісь — простір x, а вертикальна — час t (прокручується вниз). Кожна горизонтальна смуга — знімок u(x) в один момент часу. Діагональна смуга — це фронт мандрівної хвилі, нахил якого безпосередньо дає швидкість хвилі.
Чому фронт реакції зберігає свою форму?
Форма зберігається тому, що дифузія розмиває фронт (робить ширшим), а член реакції протидіє цьому (прискорений ріст на фронті загострює його). Ці два ефекти балансуються на характерній ширині ~√(D/r), надаючи хвилі стабільний профіль.
Де в природі зустрічаються фронти реакції-дифузії?
Мандрівні хімічні хвилі з'являються в реакції Белоусова-Жаботинського (спіральні хвилі в чашці), кальцієвій сигналізації в клітинах (хвилі запліднення), поширенні інвазивних видів (оригінальна мотивація Фішера), потенціалах дії в нейронах та поширенні епідемій серед населення.