Isogeometric Analysis of the Gray-Scott Reaction-Diffusion Equations for Pattern Formation on Evolving Surfaces and Applications to Human Gyrification

We propose a numerical scheme based on the principles of Isogeometric Analysis (IgA) for a geometrical pattern formation induced evolution of manifolds. The development is modelled by the use of the Gray-Scott equations for pattern formation in combination with an equation for the displacement of the manifold. The method forms an alternative to the classical finite-element method. Our method is based on partitioning the initially spherical geometry into six patches, which are mapped onto the six faces of a cube. Major advantages of the new formalism are the reconstruction of the manifold based on bicubic spline-functions used for the computation of the concentrations as well as the evolution of the mapping operator. These features give a smooth representation of the manifold which, in turn, may lead to more realistic results. The method successfully reproduces the smooth but complicated geometrical patterns found on the surface of human brains.

[1]  D. V. van Essen,et al.  A tension-based theory of morphogenesis and compact wiring in the central nervous system. , 1997, Nature.

[2]  J. Murray,et al.  On pattern formation mechanisms for lepidopteran wing patterns and mammalian coat markings. , 1981, Philosophical transactions of the Royal Society of London. Series B, Biological sciences.

[3]  Alvaro L. G. A. Coutinho,et al.  Control strategies for timestep selection in simulation of coupled viscous flow and heat transfer , 2002 .

[4]  Q Ouyang,et al.  Pattern Formation by Interacting Chemical Fronts , 1993, Science.

[5]  C Marsault,et al.  Are the brains of monozygotic twins similar? A three-dimensional MR study. , 1998, AJNR. American journal of neuroradiology.

[6]  R. A. Barrio,et al.  The Effect of Growth and Curvature on Pattern Formation , 2004 .

[7]  Giancarlo Sangalli,et al.  IsoGeometric Analysis: A New Paradigm in the Numerical Approximation of PDEs , 2016 .

[8]  Arnold Neumaier,et al.  Introduction to Numerical Analysis , 2001 .

[9]  S. Rosenberg The Laplacian on a Riemannian Manifold: An Introduction to Analysis on Manifolds , 1997 .

[10]  Jean-Francois Mangin,et al.  A Reaction-Diffusion Model of Human Brain Development , 2010, PLoS computational biology.

[11]  B. Nagorcka,et al.  The role of a reaction-diffusion system in the initiation of primary hair follicles. , 1985, Journal of theoretical biology.

[12]  R. Raghavan,et al.  A CONTINUUM MECHANICS-BASED MODEL FOR CORTICAL GROWTH , 1997 .

[13]  V. Caviness,et al.  Mechanical model of brain convolutional development. , 1975, Science.

[14]  J. E. Pearson Complex Patterns in a Simple System , 1993, Science.

[15]  A J Barkovich,et al.  Syndromes of bilateral symmetrical polymicrogyria. , 1999, AJNR. American journal of neuroradiology.

[16]  T. Sachs Patterned Differentiation in Plants , 1978 .

[17]  Y. Burnod,et al.  A morphogenetic model for the development of cortical convolutions. , 2005, Cerebral cortex.

[18]  Steven J. Ruuth,et al.  Implicit-explicit Runge-Kutta methods for time-dependent partial differential equations , 1997 .

[19]  T. Ivancevic,et al.  Applied differential geometry : a modern introduction , 2007 .

[20]  D. V. Essen,et al.  A tension-based theory of morphogenesis and compact wiring in the central nervous system , 1997, Nature.

[21]  J. Lefévre,et al.  A reaction-diffusion model of the human brain development , 2010, 2010 IEEE International Symposium on Biomedical Imaging: From Nano to Macro.

[22]  Steven J. Ruuth,et al.  Implicit-explicit methods for time-dependent partial differential equations , 1995 .

[23]  T. Hughes,et al.  Isogeometric analysis : CAD, finite elements, NURBS, exact geometry and mesh refinement , 2005 .

[24]  Pasko Rakic,et al.  Genetic Control of Cortical Convolutions , 2004, Science.