Global existence, blow-up behavior and numerical simulations for a class of chemotaxis-driven fish-mussel systems

Resumen Técnico: Existencia Global, Comportamiento de Explosión y Simulaciones Numéricas para una Clase de Sistemas de Peces y Mejillones Impulsados por Quimiotaxis

Planteamiento del Problema
Este trabajo investiga un modelo de ecosistema impulsado por quimiotaxis que describe la dinámica espaciotemporal de nutrientes ($u$), poblaciones de peces ($v$) y poblaciones de mejillones ($w$) dentro de un dominio acotado $\Omega \subset \mathbb{R}^n$. El modelo se formula como un sistema acoplado de ecuaciones en derivadas parciales (PDE) sujeto a condiciones de contorno de Neumann homogéneas. El sistema incorpora la difusión de nutrientes, la quimiotaxis de los peces hacia regiones ricas en nutrientes y el papel de filtración de los mejillones en el reciclaje de nutrientes. Las ecuaciones gobernantes (1.1) incluyen un parámetro $\tau \in \{0, 1\}$, que distingue entre un sistema parabólico-parabólico-parabólico ($\tau=1$) y un sistema parabólico-elíptico-parabólico ($\tau=0$). El estudio tiene como objetivo establecer la existencia global y la acotación de soluciones clásicas, analizar las condiciones que conducen a la explosión en tiempo finito (particularmente en tres dimensiones) y desarrollar un esquema numérico para validar estos hallazgos teóricos.

Metodología
Los autores emplean una combinación de técnicas analíticas y simulaciones numéricas:

1. Marco Analítico:

   • Existencia Local: La existencia de una solución clásica única en el tiempo local se establece mediante métodos de semigrupos y el teorema del punto fijo de Banach, apoyándose en la teoría de regularidad elíptica y parabólica.
   • Estimaciones A Priori: Para probar la existencia global, los autores derivan cotas $L^\infty$ y $L^p$. Para el caso $\tau=1$, se utiliza una función de prueba especializada $\phi(x)$ (construida mediante el Lema 2.5) para manejar el término quimiotáctico $-\chi\nabla \cdot (v\nabla u)$. Esto permite la derivación de estimaciones $L^p$ para la población de peces $v$ bajo restricciones específicas de la sensibilidad quimiotáctica $\chi$ y los límites de nutrientes.
   • Criterios de Extensibilidad: La prueba de la existencia global se basa en demostrar que la norma $L^\infty$ de la solución permanece acotada, evitando la explosión. Esto se logra combinando las estimaciones $L^p$ con las propiedades de suavizado del semigrupo de calor de Neumann.
   • Análisis de Explosión (Blow-Up): Para el caso tridimensional ($n=3$), los autores derivan una cota inferior para el tiempo de explosión $T_{max}$. Esto implica la construcción de un funcional de energía auxiliar $\Psi_\tau$ (definido en 1.8) y la derivación de una desigualdad diferencial para su evolución. Al integrar esta desigualdad, se obtiene una cota inferior explícita para el tiempo de explosión.

2. Enfoque Numérico:

   • Discretización: Se desarrolla un método de elementos finitos (FEM) utilizando elementos tetraédricos de dimensión $P1$ para la discretización espacial y el esquema de Euler hacia atrás para la integración temporal.
   • Implementación: El sistema no lineal totalmente acoplado se resuelve utilizando el resolvedor Newton line-search de PETSc SNES, con sistemas linealizados manejados por GMRES precondicionado con factorización LU. Las simulaciones se realizan utilizando FEniCSx y se visualizan en MATLAB.
   • Validación: El esquema numérico se valida mediante estudios de convergencia de malla y comparando la dinámica de largo plazo contra un modelo de ODE espacialmente homogéneo resuelto con ode45 de MATLAB.

Contribuciones Clave y Resultados

• Existencia Global y Acotación (Teorema 1.1): El artículo establece la existencia global de soluciones clásicas únicas para el sistema en dimensiones $n \geq 3$. Específicamente, si los datos iniciales satisfacen ciertas condiciones de regularidad y los parámetros cumplen la desigualdad:
  $$\max \left\{ \tau\|u_0\|_{L^\infty(\Omega)}, \frac{M}{\alpha} \right\} < \min \left\{ \frac{\pi}{\chi}\sqrt{\frac{2}{n}}, \frac{\delta_2}{\sigma_2} \right\}$$
  entonces la solución $(u, v, w)$ es uniformemente acotada en el tiempo. Este resultado es válido tanto para $\tau=0$ como para $\tau=1$.

• Estimaciones del Tiempo de Explosión (Teorema 1.2): En el caso donde las soluciones explotan en tiempo finito (específicamente en un dominio convexo 3D), los autores derivan una cota inferior computable para el tiempo de explosión $T_{max}$. Esta cota se expresa como una integral que involucra el valor inicial del funcional de energía $\Psi_\tau(0)$ y los parámetros del sistema. También se proporciona una cota inferior simplificada explícita:
  $$T_{max} \geq \frac{1}{2E_\tau \Psi_\tau(0)^2}$$

• Convergencia Numérica: Se demuestra que el esquema de elementos finitos propuesto logra una convergencia de segundo orden tanto en las normas $L^2(\Omega)$ como $L^\infty(\Omega)$, consistente con el uso de elementos lineales.

• Simulaciones Numéricas:

  • Estabilidad: Las simulaciones con parámetros moderados demuestran que el sistema converge a un estado estacionario de coexistencia, validando la estabilidad de largo plazo predicha por el análisis de ODE.
  • Comportamiento de Explosión: Experimentos numéricos en un dominio cúbico 3D con altas concentraciones iniciales y elecciones de parámetros específicas ilustran el crecimiento rápido y descontrolado de la densidad de la población de peces $v$, confirmando la posibilidad teórica de una explosión en tiempo finito. Las simulaciones muestran que mientras $v$ explota, las densidades de nutrientes $u$ y mejillones $w$ permanecen acotadas, lo cual es consistente con las cotas teóricas derivadas.

Significancia y Reivindicaciones
El artículo afirma proporcionar una base matemática rigurosa para comprender la estabilidad y la potencial inestabilidad de los ecosistemas de acuicultura que involucran peces quimiotácticos y mejillones filtradores. Al extender los modelos previos basados en ODE a un marco de PDE, el trabajo destaca el papel crítico de la heterogeneidad espacial y el movimiento dirigido (quimiotaxis) en la dinámica de los ecosistemas.

Los autores posicionan su trabajo como un puente entre el análisis teórico y el modelado computacional. Demuestran que, si bien la existencia global puede garantizarse bajo regímenes de parámetros específicos (evitando el colapso inducido por la eutrofización), el sistema es susceptible de una explosión en tiempo finito en dimensiones superiores si estas condiciones se violan. El desarrollo de un esquema de elementos finitos robusto permite la visualización de estas dinámicas complejas, ofreciendo una herramienta para verificar los criterios teóricos de explosión y simular la transición de la coexistencia estable a la explosión poblacional. El estudio no propone nuevas estrategias de gestión ecológica, sino que proporciona las herramientas analíticas y numéricas necesarias para evaluar la estabilidad de tales sistemas bajo diversas condiciones ambientales.