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

L'effet simultané de la chimiotaxie et de la taxie d'alarme sur l'existence globale et la stabilité d'un système prédateur-proie

Informations de base

  • ID de l'article: 2510.14728
  • Titre: L'effet simultané de la chimiotaxie et de la taxie d'alarme sur l'existence globale et la stabilité d'un système prédateur-proie
  • Auteurs: Gnanasekaran Shanmugasundaram (NIT Tiruchirappalli), Jitraj Saha (NIT Tiruchirappalli), Rafael Díaz Fuentes (Université de Cagliari)
  • Classification: math.AP (Analyse des EDP)
  • Date de soumission: 16 octobre 2025
  • Lien de l'article: https://arxiv.org/abs/2510.14728

Résumé

Cette étude examine un système complètement parabolique de chimiotaxie-taxie d'alarme prédateur-proie avec conditions aux limites de Neumann homogènes sur un domaine borné ΩRn\Omega \subset \mathbb{R}^n avec frontière lisse Ω\partial\Omega. Sous certaines conditions de paramètres, nous démontrons l'existence d'une unique solution classique globale bornée du système. La convergence des solutions est établie par la construction d'une fonctionnelle de Lyapunov appropriée. De plus, le comportement asymptotique des solutions est vérifié par simulation numérique. Les résultats soulignent l'importance des coefficients de chimiotaxie et de taxie d'alarme dans la détermination de l'existence et de la stabilité des modèles prédateur-proie.

Contexte et motivation de la recherche

Contexte biologique

  1. Mécanisme d'appels d'alarme: L'étude s'appuie sur l'hypothèse du « vol d'alarme » dans les écosystèmes marins, selon laquelle les proies émettent des signaux pour attirer des prédateurs secondaires, dissuadant ainsi les prédateurs primaires. Un exemple typique est la bioluminescence émise par les dinoflagellés en réponse aux perturbations causées par les copépodes, attirant les poissons comme prédateurs secondaires.
  2. Phénomènes de chimiotaxie: Le comportement de mouvement des microorganismes en réponse aux gradients chimiques, largement appliqué en écologie, immunologie, biologie et médecine.

Besoins de modélisation mathématique

Les problèmes fondamentaux que cette recherche vise à résoudre sont:

  • Comment modéliser mathématiquement un système d'interaction à trois espèces incluant simultanément la chimiotaxie et la taxie d'alarme
  • Sous quelles conditions la solution de ce système existe-t-elle globalement et reste-t-elle bornée
  • Quel est le comportement asymptotique à long terme du système

Limitations de la recherche existante

La littérature existante se concentre principalement sur:

  • L'étude de mécanismes de taxie uniques (chimiotaxie pure ou taxie d'alarme pure)
  • Les résultats locaux en dimension deux
  • L'absence d'analyse mathématique rigoureuse considérant simultanément les deux mécanismes de taxie

Contributions principales

  1. Première établissement d'une théorie mathématique rigoureuse pour l'existence globale, la bornitude et la convergence d'un système prédateur-proie incluant simultanément la chimiotaxie et la taxie d'alarme
  2. Construction d'une fonctionnelle de Lyapunov innovante pour analyser la stabilité asymptotique globale du système
  3. Preuve de l'existence d'une solution classique globale dans l'espace de dimension n1n \geq 1 et sa bornitude uniforme
  4. Établissement d'un cadre complet d'analyse de stabilité, incluant le comportement asymptotique des états de coexistence, d'extinction et de semi-coexistence
  5. Fourniture de vérifications numériques validant les résultats théoriques en dimensions deux et trois

Explication détaillée des méthodes

Modèle mathématique

L'étude considère le système couplé suivant d'équations paraboliques non linéaires:

{ut=d1Δuχ1(u(vw))+μ1u(1u+a1v+a2w)vt=d2Δvχ2(vz)+μ2v(1va3u+a4w)wt=d3Δw+ξ(wz)+μ3w(1wa5ua6v)zt=d4Δz+αv+βwγz\begin{cases} 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}

où:

  • u(x,t)u(x,t): densité du prédateur secondaire
  • v(x,t)v(x,t): densité du prédateur primaire
  • w(x,t)w(x,t): densité de la proie
  • z(x,t)z(x,t): concentration de l'attractif chimique

Signification des paramètres clés

  • χ1\chi_1: coefficient de taxie d'alarme (mouvement du prédateur secondaire selon le gradient de densité du prédateur primaire-proie)
  • χ2\chi_2: coefficient de chimiotaxie (mouvement du prédateur primaire selon le signal chimique)
  • ξ\xi: coefficient de répulsion chimique (la proie s'éloigne du signal chimique)
  • did_i: coefficients de diffusion
  • μi\mu_i: coefficients de croissance logistique
  • aia_i: coefficients d'interaction interspécifique

Points d'innovation technique

1. Stratégie de preuve de bornitude globale

Approche multi-étapes:

  • Première étape: Établissement d'estimations de normes Lp(Ω)L^p(\Omega) (Lemme 3.1)
  • Deuxième étape: Utilisation de la technique du semi-groupe de chaleur de Neumann pour obtenir la bornitude L(Ω)L^{\infty}(\Omega) (Lemme 3.2)
  • Troisième étape: Établissement d'estimations de gradients W1,(Ω)W^{1,\infty}(\Omega) (Lemmes 3.4-3.6)

2. Construction de la fonctionnelle de Lyapunov

Construction de fonctionnelles d'énergie appropriées pour différents états d'équilibre:

État de coexistence: E1(t)=Ω(uuulnuu)+Γ1Ω(vvvlnvv)+Γ2Ω(wwwlnww)+12Ω(zz)2E_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

Γ1=μ1a1μ2a3\Gamma_1 = \frac{\mu_1a_1}{\mu_2a_3}, Γ2=μ1a2μ3a5\Gamma_2 = \frac{\mu_1a_2}{\mu_3a_5}

Résultats théoriques principaux

Théorème 1.1 (Existence globale)

Il existe μ>0\mu > 0 tel que lorsque min{μ2,μ3}>μ\min\{\mu_2,\mu_3\} > \mu, pour des données initiales non négatives satisfaisant les conditions, le système (1.1) admet une unique solution classique (u,v,w,z)(u,v,w,z) uniformément bornée:

u(,t)L(Ω)+v(,t)W1,(Ω)+w(,t)W1,(Ω)+z(,t)W1,q(Ω)C\|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

Résultats de stabilité

L'étude établit quatre classes de comportements asymptotiques:

  1. État de coexistence (Théorème 1.2): Lorsque les paramètres satisfont certaines conditions, la solution converge exponentiellement vers l'état d'équilibre de coexistence (u,v,w,z)(u^*,v^*,w^*,z^*)
  2. État d'existence exclusive du prédateur secondaire (Théorème 1.3): La solution converge vers (1,0,0,0)(1,0,0,0)
  3. États de semi-coexistence (Théorème 1.4):
    • État d'extinction de la proie: (uˉ,vˉ,0,zˉ)(ū,v̄,0,z̄)
    • État d'extinction du prédateur primaire: (u^,0,w^,z^)(û,0,ŵ,ẑ)

Expériences numériques

Configuration expérimentale

  • Domaine spatial: [0.5,0.5]n[-0.5,0.5]^n (n=2,3n=2,3)
  • Maillage: 101×101101×101 en 2D, 50×50×5050×50×50 en 3D
  • Conditions initiales: Fonctions en cloche à symétrie radiale
  • Paramètres: di=χ1=χ2=ξ=μi=α=β=1d_i=\chi_1=\chi_2=\xi=\mu_i=\alpha=\beta=1, γ=2\gamma=2

Résultats numériques de vérification

Quatre exemples typiques vérifient les prédictions théoriques:

  1. Exemple 5.1: Avec ai=0.5a_i=0.5, la solution converge vers l'état de coexistence (1.286877,0.428973,0.144368,0.286793)(1.286877, 0.428973, 0.144368, 0.286793)
  2. Exemple 5.2: Sous certains paramètres, la solution converge vers (1,0,0,0)(1,0,0,0)
  3. Exemples 5.3-5.4: Vérification de la convergence des deux états de semi-coexistence

Défis techniques et solutions

Principales difficultés techniques

  1. Traitement des termes non linéaires: Complexité du terme de taxie d'alarme χ1(u(vw))\chi_1\nabla \cdot (u\nabla(vw))
  2. Analyse du système couplé: Nature fortement couplée des quatre équations
  3. Estimations de gradients: Complexité technique en dimension supérieure

Outils techniques clés

  1. Inégalité de Gagliardo-Nirenberg
  2. Théorie du semi-groupe de chaleur de Neumann
  3. Régularité maximale de Sobolev
  4. Inégalités de Young et Hölder

Comparaison avec les travaux connexes

Développement historique

  • Modèle de Keller-Segel (1970): Travail fondateur sur la chimiotaxie
  • Modèle de Haskell-Bell (2021): Première modélisation mathématique de la taxie d'alarme
  • Extensions existantes: Principalement concentrées sur des mécanismes uniques ou des cas de faible dimension

Innovations du présent article

Comparé aux travaux existants, cet article est le premier à:

  • Traiter rigoureusement l'action simultanée de deux mécanismes de taxie
  • Établir une théorie globale en dimension supérieure
  • Fournir un cadre complet d'analyse de stabilité

Conclusions et discussion

Conclusions principales

  1. Sous des conditions de paramètres appropriées, le système chemo-alarm-taxis admet une solution classique globale bornée
  2. Le comportement à long terme du système est déterminé par les relations entre paramètres, pouvant converger vers quatre états d'équilibre distincts
  3. Les coefficients de chimiotaxie et de taxie d'alarme jouent un rôle clé dans la détermination de la stabilité du système

Signification biologique

  • Valide la faisabilité mathématique de l'hypothèse du « vol d'alarme »
  • Explique la dynamique complexe des interactions à trois espèces dans la nature
  • Fournit une base théorique pour les stratégies de conservation écologique

Limitations

  1. Conditions de paramètres: Les résultats de stabilité nécessitent des contraintes de paramètres relativement fortes
  2. Analyse linéarisée: L'analyse de stabilité locale est basée sur la linéarisation
  3. Vérification numérique: Effectuée uniquement dans une plage de paramètres spécifique

Directions futures

  1. Relâcher les restrictions sur les conditions de paramètres
  2. Considérer des termes de réaction plus complexes
  3. Étudier la formation de motifs spatiaux
  4. Étendre à des environnements stochastiques

Évaluation approfondie

Points forts

  1. Rigueur théorique: Fournit un cadre complet de preuves mathématiques
  2. Innovation méthodologique: Combine habilement plusieurs techniques d'analyse
  3. Complétude des résultats: Couvre l'existence, la bornitude et la stabilité
  4. Vérification numérique: Les expériences numériques en 2D et 3D soutiennent les résultats théoriques

Insuffisances

  1. Complexité des conditions de paramètres: Vérification difficile des paramètres dans les applications pratiques
  2. Réalisme biologique: Certaines hypothèses mathématiques peuvent être trop idéalisées
  3. Complexité computationnelle: Défis du calcul numérique en dimension supérieure

Impact

  • Contribution théorique: Fournit un nouveau cadre d'analyse pour la recherche sur les systèmes de taxie
  • Valeur applicative: Fournit des outils importants pour la modélisation mathématique écologique
  • Valeur méthodologique: La méthode de construction de la fonctionnelle de Lyapunov a une valeur généralisable

Scénarios applicables

  • Modélisation des écosystèmes marins
  • Transduction de signaux en biologie cellulaire
  • Analyse de la dynamique du système immunitaire
  • Élaboration de stratégies de conservation écologique

Références

L'article cite 46 références pertinentes, couvrant plusieurs domaines incluant la théorie de la chimiotaxie, la modélisation de la taxie d'alarme, la théorie des équations aux dérivées partielles et d'autres travaux importants, fournissant une base théorique solide pour la recherche.


Évaluation globale: Cet article est un travail de haute qualité ayant une importance théorique significative dans le domaine de la biologie mathématique. Il établit pour la première fois une théorie mathématique complète du système chemo-alarm-taxis, avec des méthodes rigoureuses et des résultats profonds, jouant un rôle important dans la promotion du développement des domaines connexes.