A space time conservation element and solution element method for solving two-species chemotaxis model

Abstract This work is related to the numerical investigation of two species chemotaxis models in both one and two dimensions. This system in its simpler form is a set of non-linear partial differential equations. First equation represents the dynamics of cell densities and the others are responsible for chemoattractant concentration. The system is known to produce delta-type singularities in finite time. Therefore, a less order accurate numerical schemes are not capable of handling such a complicated spiky solution. A Space Time Conservation Element and Solution Element (CE/SE) method is proposed for solving such systems. This method has distinct attributes which includes treatment of space and time variables in unified fashion. Furthermore, numerical diffusion is reduced on the basis that both conserved quantities and derivatives are anonymous. Due to this feature, diffusion in numerical schemes is inherently reduced. The Nessyahu–Tadmor (NT) central scheme is also implemented for validation and comparison. Several one and two dimensional case studies are carried out. Numerical results obtained through proposed numerical method are analyzed in comparison to Nessyahu–Tadmor (NT) central scheme. Moreover, effects of the different densities and concentration functions are explored to see generic applicability of the proposed method for current system of equations. It is observed that the CE/SE and NT central method has the capability to capture the delta type singularities in the solution profile. However, CE/SE resolves the solution peaks better as compare to NT central scheme.

[1]  Waqas Ashraf,et al.  A SPACE–TIME CE/SE METHOD FOR SOLVING HYPERBOLIC HEAT CONDUCTION MODEL , 2014 .

[2]  Sin-Chung Chang.  The Method of Space-Time Conservation Element and Solution Element-A New Approach for Solving the Navier-Stokes and Euler Equations , 1995 .

[3]  Francis Filbet,et al.  Approximation of Hyperbolic Models for Chemosensitive Movement , 2005, SIAM J. Sci. Comput..

[4]  Sin-Chung Chang,et al.  Regular Article: The Space-Time Conservation Element and Solution Element Method: A New High-Resolution and Genuinely Multidimensional Paradigm for Solving Conservation Laws , 1999 .

[5]  J. Christiansen.  Numerical Simulation of Hydrodynamics by the Method of Point Vortices , 1997 .

[6]  Raluca Eftimie,et al.  Hyperbolic and kinetic models for self-organized biological aggregations and movement: a brief review , 2012, Journal of mathematical biology.

[7]  C. Patlak.  Random walk with persistence and external bias , 1953 .

[8]  Alexander Kurganov,et al.  ON A CHEMOTAXIS MODEL WITH SATURATED CHEMOTACTIC FLUX , 2012 .

[9]  Muhammad Yousaf,et al.  The space-time CESE method for solving special relativistic hydrodynamic equations , 2012, J. Comput. Phys..

[10]  Shamsul Qamar,et al.  Application of central schemes for solving radiation hydrodynamical models , 2013, Comput. Phys. Commun..

[11]  R. LeVeque.  Numerical methods for conservation laws , 1990 .

[12]  Sin-Chung Chang,et al.  A space-time conservation element and solution element method for solving the two- and three-dimensional unsteady euler equations using quadrilateral and hexahedral meshes , 2002 .

[13]  H. Huynh.  Accurate upwind methods for the Euler equations , 1995 .

[14]  Jianxian Qiu,et al.  Hermite WENO schemes and their application as limiters for Runge-Kutta discontinuous Galerkin method: one-dimensional case , 2004 .

[15]  Alexander Kurganov,et al.  Numerical study of two-species chemotaxis models , 2013 .

[16]  E. Tadmor,et al.  Non-oscillatory central differencing for hyperbolic conservation laws , 1990 .

[17]  T. Hillen.  HYPERBOLIC MODELS FOR CHEMOSENSITIVE MOVEMENT , 2002 .

[18]  L. Segel,et al.  Model for chemotaxis. , 1971, Journal of theoretical biology.

[19]  Wai Ming To,et al.  Application of the space-time conservation element and solution element method to one-dimensional convection-diffusion problems , 2000 .

[20]  Waqas Ashraf,et al.  A space–time CE/SE method for solving single and two-phase shallow flow models , 2014 .

[21]  Nicolas Vauchelet,et al.  Traveling Pulses for a Two-Species Chemotaxis Model , 2016, PLoS Comput. Biol..

[22]  M. Liu,et al.  The Direct Aero-Acoustics Simulation of Flow around a Square Cylinder Using the CE/SE Scheme , 2007 .

[23]  K. Painter,et al.  A User's Guide to Pde Models for Chemotaxis , 2022 .

[24]  Xijun Yu,et al.  The high order control volume discontinuous Petrov-Galerkin finite element method for the hyperbolic conservation laws based on Lax-Wendroff time discretization , 2015, Appl. Math. Comput..

[25]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[26]  Shamsul Qamar,et al.  On the application of a variant CE/SE method for solving two-dimensional ideal MHD equations , 2010 .

[27]  L. Segel,et al.  Initiation of slime mold aggregation viewed as an instability. , 1970, Journal of theoretical biology.

[28]  Alexander Kurganov,et al.  High-order positivity-preserving hybrid finite-volume-finite-difference methods for chemotaxis systems , 2018, Adv. Comput. Math..