Derivation of the Amplitude Equation for Reaction-Diffusion Systems via Computer-Aided Multiple-Scale Expansion

The amplitude equation describes a reduced form of a reaction–diffusion system, yet still retains its essential dynamical features. By approximating the analytic solution, the amplitude equation allows the examination of mode instability when the system is near a bifurcation point. Multiple-scale expansion (MSE) offers a straightforward way to systematically derive the amplitude equations. The method expresses the single independent variable as an asymptotic power series consisting of newly introduced independent variables with differing time and space scales. The amplitude equations are then formulated under the solvability conditions which remove secular terms. To our knowledge, there is little information in the research literature that explains how the exhaustive workflow of MSE is applied to a reaction–diffusion system. In this paper, detailed mathematical operations underpinning the MSE are elucidated, and the practical ways of encoding these operations using MAPLE are discussed. A semi-automated MSE computer algorithm Amp_solving is presented for deriving the amplitude equations in this research. Amp_solving has been applied to the classical Brusselator model for the derivation of amplitude equations when the system is in the vicinity of a Turing codimension-1 and a Turing–Hopf codimension-2 bifurcation points. Full open-source Amp_solving codes for the derivation are comprehensively demonstrated and available to the public domain.

[1]  Katharina Krischer,et al.  Spatio-Temporal Pattern Formation , 1999 .

[2]  M. Segev,et al.  Spontaneous pattern formation in a cavity with incoherent light. , 2004, Optics express.

[3]  De Wit A,et al.  Spatiotemporal dynamics near a codimension-two point. , 1996, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[4]  M. Shahabadi,et al.  MULTIPLE-SCALE ANALYSIS OF PLANE WAVE REFRACTION AT A DIELECTRIC SLAB WITH KERR-TYPE NONLINEARITY , 2006 .

[5]  J. Langer,et al.  Pattern formation in nonequilibrium physics , 1999 .

[6]  Derivation of amplitude equations for nonlinear oscillators subject to arbitrary forcing. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[7]  Amplitude equation for the generalized Swift–Hohenberg equation with noise , 2014 .

[8]  Pei Yu,et al.  COMPUTATION OF NORMAL FORMS VIA A PERTURBATION TECHNIQUE , 1998 .

[9]  Moira L Steyn-Ross,et al.  Gap junctions mediate large-scale Turing structures in a mean-field cortex driven by subcortical noise. , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.

[10]  Ryan A Chisholm,et al.  Critical slowing down as an indicator of transitions in two-species models. , 2009, Journal of theoretical biology.

[11]  Milos Dolnik,et al.  Spatial resonances and superposition patterns in a reaction-diffusion model with interacting Turing modes. , 2002, Physical review letters.

[12]  Pei Yu,et al.  Analysis on Double Hopf Bifurcation Using Computer Algebra with the Aid of Multiple Scales , 2002 .

[13]  Li Li,et al.  Pattern formation induced by cross-diffusion in a predator–prey system , 2008 .

[14]  Competing spatial and temporal instabilities in a globally coupled bistable semiconductor system near a codimension-two bifurcation , 2000, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[15]  István Szalai,et al.  An Experimental Design Method Leading to Chemical Turing Patterns , 2009, Science.

[16]  I. Prigogine,et al.  Formative Processes. (Book Reviews: Self-Organization in Nonequilibrium Systems. From Dissipative Structures to Order through Fluctuations) , 1977 .

[17]  Swinney,et al.  Pattern formation in the presence of symmetries. , 1994, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[18]  Grégoire Nicolis,et al.  Self-Organization in nonequilibrium systems , 1977 .

[19]  On Mode Interactions in Reaction Diffusion Equation with Nearly Degenerate Bifurcations , 1980 .

[20]  L. Pismen,et al.  ON AUTOMATED DERIVATION OF AMPLITUDE EQUATIONS IN NONLINEAR BIFURCATION PROBLEMS , 1996 .

[21]  A..  SPATIAL PATTERNS AND SPATIOTEMPORAL DYNAMICS IN CHEMICAL SYSTEMS , 2007 .

[22]  M. Cartmell,et al.  A computerised implementation of the multiple scales perturbation method using Mathematica , 2000 .

[23]  L. Pismen,et al.  New amplitude equation for Boussinesq convection and nonequilateral hexagonal patterns , 1995 .

[24]  M. Haragus,et al.  Local Bifurcations, Center Manifolds, and Normal Forms in Infinite-Dimensional Dynamical Systems , 2010 .

[25]  A. M. Turing,et al.  The chemical basis of morphogenesis , 1952, Philosophical Transactions of the Royal Society of London. Series B, Biological Sciences.

[26]  D. Walgraef.  Spatio-Temporal Pattern Formation: With Examples from Physics, Chemistry, and Materials Science , 2012 .

[27]  D. Alistair Steyn-Ross,et al.  Gap junctions modulate seizures in a mean-field model of general anesthesia for the cortex , 2012, Cognitive Neurodynamics.

[28]  L. M. Pismen,et al.  COMPUTER TOOLS FOR BIFURCATION ANALYSIS: GENERAL APPROACH WITH APPLICATION TO DYNAMICAL AND DISTRIBUTED SYSTEMS , 1999 .

[29]  A. Zhabotinsky,et al.  Pattern formation arising from interactions between Turing and wave instabilities , 2002 .

[30]  B. Peña,et al.  Selection and competition of Turing patterns , 2000 .

[31]  Rui Peng,et al.  Pattern formation in the Brusselator system , 2005 .

[32]  B. Peña,et al.  Stability of Turing patterns in the Brusselator model. , 2001, Physical review. E, Statistical, nonlinear, and soft matter physics.

[33]  Bifurcation Equations Through Multiple-Scales Analysis for a Continuous Model of a Planar Beam , 2005 .

[34]  P. Yu,et al.  SYMBOLIC COMPUTATION OF NORMAL FORMS FOR RESONANT DOUBLE HOPF BIFURCATIONS USING A PERTURBATION TECHNIQUE , 2001 .