2025-11-18T00:37:13.255978

The simultaneous effect of chemotaxis and alarm-taxis on the global existence and stability of a predator-prey system

Shanmugasundaram, Saha, Fuentes
This study examines a fully parabolic predator-prey chemo-alarm-taxis system under homogeneous Neumann boundary conditions in a bounded domain $Ω\subset \mathbb{R}^n$ with a smooth boundary $\partialΩ$. Under specific parameter conditions, it is shown that the system admits a unique, globally bounded classical solution. The convergence of the solution is established through the construction of an appropriate Lyapunov functional. In addition, numerical simulations are presented to validate the asymptotic behaviour of the solution. The results highlight the significant role of chemotaxis and alarm-taxis coefficients in determining the existence and stability of predator-prey models, as discussed in the literature.
academic

Одновременное влияние химиотаксиса и тревожного таксиса на глобальное существование и устойчивость системы хищник-жертва

Основная информация

  • ID статьи: 2510.14728
  • Название: The simultaneous effect of chemotaxis and alarm-taxis on the global existence and stability of a predator-prey system
  • Авторы: Gnanasekaran Shanmugasundaram (NIT Tiruchirappalli), Jitraj Saha (NIT Tiruchirappalli), Rafael Díaz Fuentes (University of Cagliari)
  • Классификация: math.AP (Анализ дифференциальных уравнений в частных производных)
  • Дата подачи: 16 октября 2025 г.
  • Ссылка на статью: https://arxiv.org/abs/2510.14728

Аннотация

В данном исследовании рассматривается полностью параболическая система химиотаксиса-тревожного таксиса хищник-жертва с однородными граничными условиями Неймана на ограниченной области Ω⊂Rn\Omega \subset \mathbb{R}^n с гладкой границей ∂Ω\partial\Omega. При определённых параметрических условиях доказано существование единственного глобально ограниченного классического решения системы. Сходимость решений установлена путём построения подходящего функционала Ляпунова. Кроме того, асимптотическое поведение решений проверено численным моделированием. Результаты подчёркивают важную роль коэффициентов химиотаксиса и тревожного таксиса в определении существования и устойчивости модели хищник-жертва.

Исследовательский контекст и мотивация

Биологический фон

  1. Механизм тревожных сигналов: Исследование основано на гипотезе «воровского тревожного сигнала» в морских экосистемах, согласно которой жертва издаёт сигналы, привлекающие вторичных хищников, тем самым отпугивая основного хищника. Типичный пример — биолюминесцентные сигналы динофлагеллят в ответ на возмущение копеподами, привлекающие рыб в качестве вторичных хищников.
  2. Явление химиотаксиса: Движение микроорганизмов в ответ на химические градиенты, имеющее широкое применение в экологии, иммунологии, биологии и медицине.

Потребность в математическом моделировании

Основные проблемы, которые должно решить исследование:

  • Как математически моделировать трёхвидовую систему взаимодействия, одновременно включающую химиотаксис и тревожный таксис
  • При каких условиях решения такой системы существуют глобально и остаются ограниченными
  • Каково долгосрочное асимптотическое поведение системы

Ограничения существующих исследований

Существующая литература в основном сосредоточена на:

  • Исследовании отдельных механизмов таксиса (чистый химиотаксис или чистый тревожный таксис)
  • Локальных результатах в двумерном пространстве
  • Отсутствии строгого математического анализа, одновременно учитывающего оба механизма таксиса

Основные вклады

  1. Впервые установлена строгая математическая теория глобального существования, ограниченности и сходимости системы хищник-жертва, одновременно содержащей химиотаксис и тревожный таксис
  2. Построен инновационный функционал Ляпунова для анализа глобальной асимптотической устойчивости системы
  3. Доказано существование глобального классического решения в пространстве размерности n≥1n \geq 1 с равномерной ограниченностью
  4. Установлена полная система анализа устойчивости, включающая асимптотическое поведение состояний сосуществования, вымирания и полусосуществования
  5. Предоставлена численная верификация результатов теории в двумерном и трёхмерном пространстве

Подробное описание методов

Математическая модель

Исследование рассматривает следующую систему связанных нелинейных параболических дифференциальных уравнений в частных производных:

u_t = d_1\Delta u - \chi_1\nabla \cdot (u\nabla(vw)) + \mu_1u(1-u+a_1v+a_2w) \\ v_t = d_2\Delta v - \chi_2\nabla \cdot (v\nabla z) + \mu_2v(1-v-a_3u+a_4w) \\ w_t = d_3\Delta w + \xi\nabla \cdot (w\nabla z) + \mu_3w(1-w-a_5u-a_6v) \\ z_t = d_4\Delta z + \alpha v + \beta w - \gamma z \end{cases}$$ где: - $u(x,t)$: плотность вторичного хищника - $v(x,t)$: плотность основного хищника - $w(x,t)$: плотность жертвы - $z(x,t)$: концентрация химического аттрактанта ### Значение ключевых параметров - $\chi_1$: коэффициент тревожного таксиса (движение вторичного хищника по градиентам плотности основного хищника и жертвы) - $\chi_2$: коэффициент химиотаксиса (движение основного хищника по химическому сигналу) - $\xi$: коэффициент химического отталкивания (удаление жертвы от химического сигнала) - $d_i$: коэффициенты диффузии - $\mu_i$: коэффициенты логистического роста - $a_i$: коэффициенты межвидового взаимодействия ### Технические инновации #### 1. Стратегия доказательства глобальной ограниченности Применяется многошаговый подход: - **Первый шаг**: Установление оценок норм $L^p(\Omega)$ (лемма 3.1) - **Второй шаг**: Использование техники полугруппы Неймана для получения ограниченности $L^{\infty}(\Omega)$ (лемма 3.2) - **Третий шаг**: Установление оценок градиентов $W^{1,\infty}(\Omega)$ (леммы 3.4-3.6) #### 2. Конструкция функционала Ляпунова Для различных равновесных состояний строятся соответствующие энергетические функционалы: **Состояние сосуществования**: $$E_1(t) = \int_\Omega \left(u-u^*-u^*\ln\frac{u}{u^*}\right) + \Gamma_1\int_\Omega \left(v-v^*-v^*\ln\frac{v}{v^*}\right) + \Gamma_2\int_\Omega \left(w-w^*-w^*\ln\frac{w}{w^*}\right) + \frac{1}{2}\int_\Omega (z-z^*)^2$$ где $\Gamma_1 = \frac{\mu_1a_1}{\mu_2a_3}$, $\Gamma_2 = \frac{\mu_1a_2}{\mu_3a_5}$ ## Основные теоретические результаты ### Теорема 1.1 (Глобальное существование) Существует $\mu > 0$ такое, что при $\min\{\mu_2,\mu_3\} > \mu$ для неотрицательных начальных данных, удовлетворяющих условиям, система (1.1) имеет единственное классическое решение $(u,v,w,z)$, которое равномерно ограничено: $$\|u(\cdot,t)\|_{L^{\infty}(\Omega)} + \|v(\cdot,t)\|_{W^{1,\infty}(\Omega)} + \|w(\cdot,t)\|_{W^{1,\infty}(\Omega)} + \|z(\cdot,t)\|_{W^{1,q}(\Omega)} \leq C$$ ### Результаты устойчивости Исследование устанавливает четыре класса асимптотического поведения: 1. **Состояние сосуществования** (теорема 1.2): При выполнении определённых параметрических условий решение экспоненциально сходится к равновесному состоянию сосуществования $(u^*,v^*,w^*,z^*)$ 2. **Состояние единственного вторичного хищника** (теорема 1.3): Решение сходится к $(1,0,0,0)$ 3. **Состояния полусосуществования** (теорема 1.4): - Состояние исчезновения жертвы: $(ū,v̄,0,z̄)$ - Состояние исчезновения основного хищника: $(û,0,ŵ,ẑ)$ ## Численные эксперименты ### Установка экспериментов - **Пространственная область**: $[-0.5,0.5]^n$ ($n=2,3$) - **Сетка**: $101×101$ для 2D, $50×50×50$ для 3D - **Начальные условия**: Радиально симметричные функции в форме колокола - **Параметры**: $d_i=\chi_1=\chi_2=\xi=\mu_i=\alpha=\beta=1$, $\gamma=2$ ### Результаты численной верификации Четыре типичных примера подтверждают теоретические предсказания: 1. **Пример 5.1**: При $a_i=0.5$ решение сходится к состоянию сосуществования $(1.286877, 0.428973, 0.144368, 0.286793)$ 2. **Пример 5.2**: При определённых параметрах решение сходится к $(1,0,0,0)$ 3. **Примеры 5.3-5.4**: Подтверждена сходимость двух типов состояний полусосуществования ## Технические вызовы и решения ### Основные технические трудности 1. **Обработка нелинейных членов**: Сложность члена тревожного таксиса $\chi_1\nabla \cdot (u\nabla(vw))$ 2. **Анализ связанной системы**: Сильная связанность четырёх уравнений 3. **Оценки градиентов**: Техническая сложность в пространстве высокой размерности ### Ключевые технические инструменты 1. **Неравенство Гальярдо-Ниренберга** 2. **Теория полугруппы Неймана** 3. **Максимальная регулярность Соболева** 4. **Неравенства Янга и Гёльдера** ## Сравнение с существующими работами ### Историческое развитие - **Модель Келлера-Сегела** (1970): Основополагающая работа по химиотаксису - **Модель Хаскелла-Белла** (2021): Первое математическое моделирование тревожного таксиса - **Существующие расширения**: В основном сосредоточены на отдельных механизмах или низкомерных случаях ### Инновации данной работы По сравнению с существующими работами, данная статья впервые: - Строго рассматривает одновременное действие обоих механизмов таксиса - Устанавливает глобальную теорию в пространстве высокой размерности - Предоставляет полную систему анализа устойчивости ## Заключение и обсуждение ### Основные выводы 1. При надлежащих параметрических условиях система хемо-тревожного таксиса имеет глобально ограниченное классическое решение 2. Долгосрочное поведение системы определяется параметрическими соотношениями и может сходиться к четырём различным равновесным состояниям 3. Коэффициенты химиотаксиса и тревожного таксиса играют ключевую роль в определении устойчивости системы ### Биологическое значение - Подтверждает математическую осуществимость гипотезы «воровского тревожного сигнала» - Объясняет сложную динамику трёхвидового взаимодействия в природе - Предоставляет теоретическую основу для стратегий охраны окружающей среды ### Ограничения 1. **Параметрические условия**: Результаты устойчивости требуют достаточно строгих параметрических ограничений 2. **Линеаризованный анализ**: Анализ локальной устойчивости основан на линеаризации 3. **Численная верификация**: Проведена только в определённом диапазоне параметров ### Направления будущих исследований 1. Ослабление ограничений на параметры 2. Рассмотрение более сложных реакционных членов 3. Исследование формирования пространственных структур 4. Расширение на стохастическую среду ## Глубокая оценка ### Преимущества 1. **Теоретическая строгость**: Предоставляет полную систему математических доказательств 2. **Методологическая инновативность**: Искусно объединяет множество аналитических техник 3. **Полнота результатов**: Охватывает существование, ограниченность и устойчивость 4. **Численная верификация**: Численные эксперименты в 2D и 3D поддерживают теоретические результаты ### Недостатки 1. **Сложность параметрических условий**: Трудность верификации параметров при практическом применении 2. **Биологическая реалистичность**: Некоторые математические предположения могут быть чрезмерно идеализированы 3. **Вычислительная сложность**: Вызовы численных вычислений в случае высокой размерности ### Влияние - **Теоретический вклад**: Предоставляет новую аналитическую систему для исследования систем таксиса - **Прикладная ценность**: Предоставляет важный инструмент для математического моделирования в экологии - **Методологическая ценность**: Метод конструкции функционала Ляпунова имеет потенциал для обобщения ### Применимые сценарии - Моделирование морских экосистем - Анализ передачи сигналов в клеточной биологии - Анализ динамики иммунной системы - Разработка стратегий охраны окружающей среды ## Библиография Статья цитирует 46 связанных источников, охватывающих важные работы в нескольких областях, включая теорию химиотаксиса, моделирование тревожного таксиса и теорию дифференциальных уравнений в частных производных, обеспечивая прочную теоретическую основу для исследования. --- **Общая оценка**: Это высококачественная статья с важным теоретическим значением в области математической биологии, впервые устанавливающая полную математическую теорию системы хемо-тревожного таксиса, с строгими методами и глубокими результатами, имеющими важное значение для развития соответствующих областей.