Bifurcations and Dynamics Emergent From Lattice and Continuum Models of Bioactive Porous Media

We study the dynamics emergent from a two-dimensional reaction–diffusion process modeled via a finite lattice dynamical system, as well as an analogous PDE system, involving spatially nonlocal interactions. These models govern the evolution of cells in a bioactive porous medium, with the evolution of the local cell density depending on a coupled quasi-static fluid flow problem. We demonstrate differences emergent from the choice of a discrete lattice or a continuum for the spatial domain of such a process. We find long-time oscillations and steady states in cell density in both lattice and continuum models, but that the continuum model only exhibits solutions with vertical symmetry, independent of initial data, whereas the finite lattice admits asymmetric oscillations and steady states arising from symmetry-breaking bifurcations. We conjecture that it is the structure of the finite lattice which allows for more complicated asymmetric dynamics. Our analysis suggests that the origin of both types of oscillations is a nonlocal reaction–diffusion mechanism mediated by quasi-static fluid flow.

[1]  R. Seydel.  Practical Bifurcation and Stability Analysis , 1994 .

[2]  David H. Sattinger,et al.  Bifurcation and symmetry breaking in applied mathematics , 1980 .

[3]  Philippe Baveye,et al.  Computational pore network modeling of the influence of biofilm permeability on bioclogging in porous media , 2008, Biotechnology and bioengineering.

[4]  D. Hutmacher,et al.  Scaffolds in tissue engineering bone and cartilage. , 2000, Biomaterials.

[5]  P. Grassberger,et al.  Symmetry breaking bifurcation for coupled chaotic attractors , 1991 .

[6]  Ian Stewart,et al.  Some Curious Phenomena in Coupled Cell Networks , 2004, J. Nonlinear Sci..

[7]  Dirk Drasdo,et al.  Individual-based and continuum models of growing cell populations: a comparison , 2009, Journal of mathematical biology.

[8]  Lenya Ryzhik,et al.  On the nonlocal Fisher–KPP equation: steady states, spreading speed and global bounds , 2013, 1307.3001.

[9]  Hiroshi Matano,et al.  Convergence of solutions of one-dimensional semilinear parabolic equations , 1978 .

[10]  Changsheng Liu,et al.  Biomimetic porous scaffolds for bone tissue engineering , 2014 .

[11]  J. Mathiot,et al.  Modelling by percolation theory of the behaviour of natural coral used as bone substitute , 2003, Physics in medicine and biology.

[12]  Hermann Riecke,et al.  Symmetry-breaking Hopf bifurcation in anisotropic systems , 1992 .

[13]  Jesús María Sanz-Serna,et al.  Convergence of method of lines approximations to partial differential equations , 1984, Computing.

[14]  Hal L. Smith,et al.  Monotone Dynamical Systems: An Introduction To The Theory Of Competitive And Cooperative Systems (Mathematical Surveys And Monographs) By Hal L. Smith , 1995 .

[15]  Fernando Fraternali,et al.  A discrete-to-continuum approach to the curvatures of membrane networks and parametric surfaces , 2014 .

[16]  Mark Newman,et al.  Networks: An Introduction , 2010 .

[17]  Steven Robert McDougall,et al.  The application of network modelling techniques to multiphase flow in porous media , 1995, Petroleum Geoscience.

[18]  R. Baker,et al.  Discrete and Continuum Approximations for Collective Cell Migration in a Scratch Assay with Cell Size Dynamics , 2018, Bulletin of mathematical biology.

[19]  T. Randolph Beard,et al.  The Global Organ Shortage: Economic Causes, Human Consequences, Policy Responses , 2013 .

[20]  S. McDougall,et al.  Mathematical modelling of flow through vascular networks: Implications for tumour-induced angiogenesis and chemotherapy strategies , 2002, Bulletin of mathematical biology.

[21]  Richard S. Graham,et al.  A continuum model of cell proliferation and nutrient transport in a perfusion bioreactor. , 2013, Mathematical medicine and biology : a journal of the IMA.

[22]  Martin Golubitsky,et al.  Patterns in Square Arrays of Coupled Cells , 1997 .

[23]  Guanrong Chen,et al.  Spatio-temporal patterns in a square-lattice Hodgkin-Huxley neural network , 2006 .

[24]  P. Mañas,et al.  Morphological and Physiological Changes Induced by High Hydrostatic Pressure in Exponential- and Stationary-Phase Cells of Escherichia coli: Relationship with Cell Death , 2004, Applied and Environmental Microbiology.

[25]  L Preziosi,et al.  A review of mathematical models for the formation of vascular networks. , 2013, Journal of theoretical biology.

[26]  Helen M. Byrne,et al.  Continuum Modelling of In Vitro Tissue Engineering: A Review , 2012 .

[27]  H. Schättler,et al.  Local bifurcations and feasibility regions in differential-algebraic systems , 1995, IEEE Trans. Autom. Control..

[28]  Boris A. Malomed,et al.  Spontaneous symmetry breaking in Schrödinger lattices with two nonlinear sites , 2011, 1105.0266.

[29]  Liesbet Geris,et al.  Computational modeling in tissue engineering , 2013 .

[30]  Ragnar Fleischmann,et al.  Exponentially fragile PT symmetry in lattices with localized eigenmodes. , 2009, Physical review letters.

[31]  Matthias Wolfrum,et al.  The Turing bifurcation in network systems: Collective patterns and single differentiated nodes , 2012 .

[32]  B. A. Malomed,et al.  Spontaneous symmetry breaking in photonic lattices: Theory and experiment , 2004, cond-mat/0412381.

[33]  Kambiz Vafai,et al.  Porous media : applications in biological systems and biotechnology , 2010 .

[34]  Fordyce A. Davidson,et al.  Spatio-temporal pattern formation in a nonlocal reaction-diffusion equation , 2001 .

[35]  R. G. Casten,et al.  Instability results for reaction diffusion equations with Neumann boundary conditions , 1978 .

[36]  Ioannis G Kevrekidis,et al.  From discrete to continuum models of three-dimensional deformations in epithelial sheets. , 2015, Biophysical journal.

[37]  A. Duncan,et al.  Noise-induced multistability in chemical systems: Discrete versus continuum modeling. , 2014, Physical review. E, Statistical, nonlinear, and soft matter physics.

[38]  A. V. Shapovalov,et al.  Asymptotics semiclassically concentrated on curves for the nonlocal Fisher–Kolmogorov–Petrovskii–Piskunov equation , 2016 .

[39]  E. C. Childs.  Dynamics of fluids in Porous Media , 1973 .

[40]  Julian R Jones,et al.  Quantifying the 3D macrostructure of tissue scaffolds , 2009, Journal of materials science. Materials in medicine.

[41]  Sophie C Cox,et al.  3D printing of porous hydroxyapatite scaffolds intended for use in bone tissue engineering applications. , 2015, Materials science & engineering. C, Materials for biological applications.

[42]  M. Golubitsky,et al.  Bifurcation on the hexagonal lattice and the planar Bénard problem , 1983, Philosophical Transactions of the Royal Society of London. Series A, Mathematical and Physical Sciences.

[43]  T. Janssen,et al.  Bifurcations of lattice structures , 1983 .

[44]  L. Preziosi,et al.  Modelling Solid Tumor Growth Using the Theory of Mixtures , 2001, Mathematical medicine and biology : a journal of the IMA.

[45]  Shui-Nee Chow,et al.  DYNAMICS OF LATTICE DIFFERENTIAL EQUATIONS , 1996 .

[46]  M. Chaplain,et al.  Continuous and discrete mathematical models of tumor-induced angiogenesis , 1998, Bulletin of mathematical biology.

[47]  Hiroko Kamei,et al.  The Existence and Classification of Synchrony-Breaking bifurcations in Regular Homogeneous Networks Using Lattice Structures , 2009, Int. J. Bifurc. Chaos.

[48]  David M. Winterbottom.  sqrt(5):2 Mode Interactions on a Square Lattice , 2004, Int. J. Bifurc. Chaos.

[49]  Cleo Choong,et al.  Three-dimensional scaffolds for tissue engineering applications: role of porosity and pore size. , 2013, Tissue engineering. Part B, Reviews.

[50]  Nicola Elvassore,et al.  Mathematical Modeling of Three-Dimensional Cell Cultures in Perfusion Bioreactors , 2006 .

[51]  E. Knobloch,et al.  Symmetry and Symmetry-Breaking Bifurcations in Fluid Dynamics , 1991 .

[52]  Kangsun Lee,et al.  Design of pressure-driven microfluidic networks using electric circuit analogy. , 2012, Lab on a chip.

[53]  John Billingham,et al.  Dynamics of a strongly nonlocal reaction–diffusion population model , 2004 .

[54]  Spencer Brown,et al.  Tissue Engineering (Academic Press Series in Biomedical Engineering) , 2010 .

[55]  J P Whiteley,et al.  The influence of hydrostatic pressure on tissue engineered bone development. , 2016, Journal of theoretical biology.

[56]  Nikolaos S. Christodoulou,et al.  Discrete Hopf bifurcation for Runge-Kutta methods , 2008, Appl. Math. Comput..

[57]  Sundararajan V. Madihally,et al.  Applications of Computational Modelling and Simulation of Porous Medium in Tissue Engineering , 2016, Comput..

[58]  Chris Antonopoulos,et al.  Chaos in a Near-integrable Hamiltonian Lattice , 2002, Int. J. Bifurc. Chaos.

[59]  Edgar Knobloch,et al.  Hopf bifurcation on a square lattice , 1991 .

[60]  M. Blunt.  Flow in porous media — pore-network models and multiphase flow , 2001 .

[61]  Andrew L. Krause,et al.  Lattice and continuum modelling of a bioactive porous tissue scaffold. , 2017, Mathematical medicine and biology : a journal of the IMA.