🔥 Фронт Реакції — Хімічна Мандрівна Хвиля
Симулюйте мандрівний фронт хвилі за рівнянням Фішера-КПП. Спостерігайте поширення автокаталітичних реакцій зі швидкістю v = 2√(Dr).
Схожі симуляції
Про симуляцію «Фронт реакції»
Рівняння Фішера-КПП — назване на честь Рональда Фішера (1937) та незалежно Колмогорова, Петровського і Піскунова — одне з найфундаментальніших диференціальних рівнянь у математичній біології та хімії. Воно поєднує два фізичні процеси: дифузію, яка поширює речовину в просторі, і логістичне зростання, що моделює автокаталітичне підсилення, яке насичується, коли реагент вичерпується. Результатом є мандрівна хвиля, що безкінечно зберігає свою форму, поширюючись з постійною швидкістю, яка повністю визначається двома параметрами D і r.
Ця симуляція чисельно розв'язує одновимірне рівняння Фішера-КПП за допомогою методу скінченних різниць з явним методом Ейлера за часом. Просторова область складається з 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) — найповільніша можлива стабільна хвиля, яка обирається при t→∞.
Що таке просторово-часова діаграма (кімограф)?
Вставка показує кімограф: 2D-зображення, де горизонтальна вісь — це простір x, а вертикальна вісь — час t (прокручується вниз). Кожна горизонтальна смуга — це знімок u(x) в один момент часу. Діагональна смуга — це мандрівний фронт хвилі, нахил якої безпосередньо дає швидкість хвилі.
Чому фронт реакції зберігає свою форму?
Форма зберігається тому, що дифузія прагне розширити фронт (робить його ширшим), тоді як реакційний член протидіє цьому розширенню (швидше зростання на фронті робить його крутішим). Ці два ефекти врівноважуються на характерній ширині ~√(D/r), надаючи хвилі стабільний профіль.
Де в природі зустрічаються хвилі реакції-дифузії?
Мандрівні хімічні хвилі зустрічаються в реакції Білоусова–Жаботинського (спіральні хвилі в чашках), кальцієвій сигналізації в клітинах (хвилі запліднення), поширенні інвазивних видів (первинна мотивація Фішера), потенціалах дії в нейронах та поширенні епідемій серед популяцій.