UK 
EnglishТипова
Polski🌊 Хвиля Популяції — Просторове Поширення Виду
Симулюйте просторове поширення виду за Fisher-KPP на 2D сітці. Фронт хвилі рухається зі швидкістю v=2√(Dr). Ефект Аллі.
Схожі симуляції
Про симуляцію
Ця симуляція моделює просторове поширення інвазивного виду за допомогою рівняння реакції-дифузії Фішера-КПП, ∂u/∂t = D∇²u + ru(1−u), яке розв'язується чисельно на сітці 150×150. Невелика центральна колонія розширюється назовні у вигляді біжучої хвилі, швидкість якої збігається до теоретичного мінімуму v = 2√(Dr), а необов'язковий доданок ефекту Аллі може сповільнити або зупинити вторгнення нижче порогової щільності.
Щільність популяції u(x,y) еволюціонує на сітці 150×150 через явні скінченні різниці з п'ятиточковим шаблоном лапласіана та нульовим потоком (умова Неймана) на межах. Синій колір позначає порожнє середовище (u=0), зелений — половину ємності (u=0,5), а червоний — насичену колонію (u=1). Смужка під сіткою показує загальну популяцію в часі, а живі показники порівнюють виміряну швидкість хвилі з теоретичним значенням 2√(Dr). Регулюйте коефіцієнт дифузії D і швидкість росту r, щоб змінити темп розсіювання й розмноження особин — теоретична швидкість хвилі оновлюється миттєво. Увімкніть ефект Аллі та підвищуйте поріг a, щоб побачити, як вторгнення сповільнюється або повністю зупиняється. Перемикач «Неоднорідний ландшафт» рандомізує локальну дифузійність на ±50%, створюючи нерівний, пальцеподібний фронт замість гладкого кола.
Та сама формула v = 2√(Dr), що керує цією змодельованою колонією, підтверджена на реальних вторгненнях — ондатри, що поширювалися Європою після інтродукції 1905 року, просувалися приблизно на 1 км/рік, що відповідає хвилі Фішера-КПП із D≈51 км²/рік. Та сама математика описує й поширення каштанового раку, ВІЛ у лімфовузлах і навіть історичне поширення землеробства неолітичною Європою.
Поширені запитання
Які повзунки керують симуляцією і що вони змінюють?
Коефіцієнт дифузії D (0,1–3,0) задає, наскільки швидко особини розсіюються за рахунок випадкового руху; швидкість росту r (0,1–2,0) задає логістичну швидкість розмноження; а поріг Аллі a (0–0,5) задає мінімальний штраф щільності, що застосовується, коли увімкнено ефект Аллі. І D, і r напряму впливають на теоретичну швидкість хвилі, показану як 2√(Dr), тож зміна будь-якого повзунка миттєво оновлює цей показник, а «Скинути» пересіває свіжу колонію радіусом 8 клітин у центрі сітки 150×150.
Як симуляція насправді обчислює хвилю щокадру?
Кожен кадр виконує чотири явні кроки Ейлера рівняння Фішера-КПП. Лапласіан наближається п'ятиточковим шаблоном (сума чотирьох сусідів мінус чотири рази центральна клітина, поділена на Δx²=1), крок часу обмежується dt = min(0,1; 1/(4D)) для дотримання умови стійкості CFL двовимірної дифузії, а ріст використовує або ru(1−u), або, з увімкненим ефектом Аллі, ru(u−a)(1−u). Межі мають нульовий потік, тож щільність не може покинути сітку.
Чим виміряна швидкість хвилі відрізняється від показника «Теорія 2√(Dr)»?
Симуляція відстежує радіус, на якому щільність вперше падає до 0,5, скануючи назовні від центру, а потім ділить зміну цього радіуса на витрачений змодельований час, щоб отримати виміряну швидкість, яка оновлюється приблизно кожні пів одиниці часу. Показник «Теорія 2√(Dr)» — це замкнена формула передбачення Фішера, обчислена безпосередньо з поточних значень повзунків D і r. На початку прогону виміряна швидкість нижча за теоретичну, поки фронт стає крутішим зі свого початкового гладкого профілю; після короткого перехідного періоду обидва значення сходяться.
Що насправді змінює увімкнення ефекту Аллі в моделі?
З вимкненим перемикачем Аллі ріст описується стандартним логістичним доданком ru(1−u), і будь-яка додатна щільність зрештою заповнює всю сітку. Увімкнення змінює ріст на ru(u−a)(1−u): щойно локальна щільність падає нижче порогу Аллі a, цей доданок стає від'ємним, тож ділянки нижче a скорочуються замість росту. Підвищення a повзунком робить ефективний поріг вторгнення суворішим, а достатньо високе значення може повністю зупинити фронт або призвести до колапсу колонії замість поширення.
Що відбувається, коли вмикається «Неоднорідний ландшафт»?
Увімкнення будує карту дифузійності для кожної клітини, де кожна з клітин сітки 150×150 отримує D, помножений на випадковий коефіцієнт від 0,5 до 1,5, згенерований один раз при скиданні. Під час обчислення кроків симуляція використовує це локальне значення замість однорідного D усюди, де активний перемикач, тож фронт хвилі просувається швидше через ділянки з високою дифузійністю та відстає в ділянках з низькою дифузійністю — створюючи той самий тип нерівного, пальцеподібного поширення, що спостерігається, коли реальні вторгнення слідують дорогами, річками чи лісовими коридорами, а не рівною однорідною місцевістю.