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

In this work, we investigate a chemotaxis-driven fish–mussel ecosystem model described by a coupled system of partial differential equations subject to homogeneous Neumann boundary conditions. Under suitable assumptions on the system parameters and initial data, we establish the global existence of classical solutions by employing semigroup methods together with a priori estimates. We also examine the possible blow-up behavior of solutions in a three-dimensional domain. To support the theoretical analysis, a finite element method is developed for the numerical approximation of the system, and convergence studies based on mesh refinement are carried out to verify the accuracy and stability of the proposed numerical scheme. Furthermore, numerical simulations illustrating the blow-up behavior of solutions in the computational domain are presented.