Biological pattern formation in cell dynamics under cross-diffusion: An Isogeometric analysis perspective

• Published in: Mathematical biosciences Vol. 384; p. 109444
• Main Author: Asmouh, Ilham
• Format: Journal Article
• Language: English
• Published: United States Elsevier Inc 01.06.2025
• Subjects: Algorithms; Animals; BDF2; Cell dynamics; Diffusion; Diffusion-reaction systems; Energy dissipation; Isogeometric analysis; Models, Biological; Morphogenesis - physiology; Nonlinear Dynamics; Pattern formation; Isogeometric analysis; BDF2; Pattern formation; Energy dissipation; Cell dynamics; Diffusion-reaction systems
• ISSN: 0025-5564; 1879-3134; 1879-3134
• DOI: 10.1016/j.mbs.2025.109444

This note presents an efficient numerical method based on isogeometric analysis (IgA) and an operator splitting approach for solving nonlinear reaction–diffusion systems with cross-diffusion. Such problems are often used in mathematical modeling of developmental biology and are subject to highly rigid reactive and diffusive terms. Similarly, the interactions between substances produce complex morphologies (Roth, 2011) [1]. In this note we present two different types of solutions. Mainly, the Turing patterns and the traveling waves, which are a direct result of the presence of linear diffusion and/or cross-diffusion in the dynamical system. To deal with the multiphysical nature of the nonlinear system, we propose a time-splitting method. The spatial discretization is performed using IgA-based Non-Uniform Ratinal B-spline (NURBS) functions, where the semidiscrete problem is integrated using an implicit scheme. The nonlinear terms are treated by an adaptive fourth-order Runge–Kutta method. The well-known FitzHugh–Nagumo and Gray–Scott models are used to study the performance of the new method. The results obtained demonstrate the ability of our algorithm to accurately maintain the shape of the solution in the presence of complex patterns arising from biological cells on complex geometries. Furthermore, the energy dissipation in the Allen-Cahn equation is analyzed and the new method clarifies the effect of the geometry on the formed patterns and on the energy decay for the considered benchmarks. •Strang splitting separates reaction and diffusion for efficient computation.•Isogeometric analysis ensures accurate diffusion discretization on complex geometries.•BDF2 integrates diffusion, while Runge–Kutta–Fehlberg handles reaction adaptively.•Energy decay in Allen–Cahn is validated with analytical comparisons.•The approach works well on irregular geometries.